浦丰投针试验
一，试验内容
试验步骤是：
1） 取一张白纸，在上面画上许多条间距为 d 的平行线。
2） 取一根长度为 l（l<d） 的针，随机地向画有平行直线的纸上掷 n
次，观察针与直线相交的次数，记为 m
3）计算针与直线相交的概率．
其概率的值 freq=2*l/( π*d)
根据这一公式和大量测试结果，即可估算出圆周率 π =2*l/(d*freq)。
下面我们将用 MATLAB 来模拟这一试验，估算出圆周率的值，并探讨
影响其准确率的因素。
二，试验过程
试验程序如下：
>> clear
>> d=3;
>> r=1;
>> NoI=0;
>> n=100;
>> x=d/2*rand(1,n);
>> p=pi*rand(1,n);
>> for i=1:n
if x(i)<r*sin(p(i))
NoI=NoI+1
end
end;
>> freq=NoI/n;
>>
estpi=4*r/(d*freq);
>> delta=abs(estpi-
pi)/pi;
>> freq
freq =
0.4700
>> estpi
estpi =
2.8369
>> delta=
0.0970
注：d 为两平行线之间的距离，r 为针的长度的一半，NoI（numbers of intersction）
表示相交次数,n 表示投针次数，相交概率为 freq，估算出的 π 的值为 estpi，最后是计算
出来的相对误差。
资源评论
