matlab模拟蒲丰掷针实验计算pi的值

matlab模拟蒲丰掷针实验计算pi的值

%投掷次数为n,针和直线相交的次数为m

%分析知真和直线相交的概率为p=2l/pid

%则有pi=2nl/md

clear all

clc

d=1;% 设置两条平行线之间的距离,相当于在一张纸上画很多间距为d的平行线

l=0.6;% 投针的长度（l

n=1000000;% n为投掷次数

x=unifrnd(0,d/2,1,n);%产生n个(0,d/2)之间均匀分布的随机数，这里d/2是投针的中点到最近的平行线的距离

alpha=unifrnd(0,pi,1,n);% 产生n个(0,pi)之间均匀分布的随机数，这里pi是投针到最近的平行线的角度

m=0;%记针与平行线相交的次数初始值为0

for i=1:n

if x(i)

m=m+1;

end

end

p=m/n; % 计算相交的频率，即相交次数比总次数

Vpi=2*l/(d*p) % 从相交的频率总求的Value pi,即pi的值

%运行结果：

%Vpi =

% 3.140710093614099

%因每次产生的随机数都不一样，每次运算出来的pi的值都不同。n越大，越精确。