蒲丰投针matlab代码

蒲丰投针实验（Poisson's needle problem）是一个经典的统计学问题，用来估计圆周率π的值，通过随机投掷细针并计算落在圆内的次数。在MATLAB中，你可以使用循环和随机数生成来模拟这个实验。以下是一个简单的MATLAB代码示例：

% 定义常量
num_trials = 100000; % 投掷次数
radius = 1; % 圆的半径
circle_area = pi * radius^2; % 圆的面积

% 初始化变量
inside_circle = 0; % 落在圆内的针的数量

% 进行投针实验
for i = 1:num_trials
    % 生成随机点 (x, y) 在单位正方形内
    x = rand();
    y = rand();
    
    % 计算点到圆心的距离
    distance = sqrt(x^2 + y^2);
    
    % 判断针是否落在圆内
    if distance <= radius
        inside_circle = inside_circle + 1;
    end
end

% 估算π
approx_pi = 4 * inside_circle / num_trials;

% 打印结果
fprintf('Estimated value of pi: %f\n', approx_pi);