蒲丰投针实验模拟.doc

概率论与数理统计实验蒲丰投针与蒙特卡罗法班级 应数12级01班学号 2012444086姓名 张旭东蒲丰投针与蒙特卡罗法张旭东 2012444086（重庆科技学院 数学与应用数学 ，重庆 沙坪坝） 【摘 要】 通过设计一个投针实验使这个事件的概率和未知量有关，然后通过重复实验，以频率估计概率，即可求得未知参数 的近似解。这种方法称为随机模拟法，也称为蒙特卡罗法。一般来说，实验次数越多所得的近似值就越接近真值。可以利用MATLAB来大量重复地模拟所设计的随机实验。 【关键词】 随机模拟；投针实验；重复实验1 引言 蒲丰投针问题是由法国科学家蒲丰（Buffon）在1777年提出的，它是概率中非常有代表性的问题，它是第一个用几何形式表达概率问题的例子，其结论具有很强的理论与实际意义。蒲丰针问题的解决不仅较典型的反应了集合概率的特征及处理方法，而且还可以由此领略到从“概率土壤”上开出的一朵瑰丽的鲜花蒙特卡洛（Monte-Carlo）方法。蒙特卡罗(Monte Carlo)方法，也称计算机模拟方法，是一种基于“随机数”的计算方法，大数定律为近年来发展迅速的随机计算机和随机模拟方法提供了理论基础。MATLAB是一个适合多学科，具有多种工作平台的功能强大的大型软件。MATLAB已经成为线性代数、自动控制理论、数理统计、数字信号处理、时间序列分析、动态系统仿真等高级课程的进本教学工具，Matlab随机数发生器的种类丰富且用法简便。本文介绍了利用随机模拟方法和大数定律的相关理论解决蒲丰投针问题计算的近似值。2 有关数学实验的有关基础 定理（贝努力大数定律） 设是n重贝努力实验中事件A出现的次数，P是事件A每次实验中出现的概率，即P(A)=p,则对任意的0，有 3 实验 蒲丰投针问题 在平面上画有等距离的一些平行线，平行线间的距离为a(a0)，向平面上随机投一长为l(la)的针。求此针与平行线相交的概率？解 以x表示针的中点与最近一条平行线的距离，又以表示针与此直线间的交角，见图1.2.易知样本空间满足 由这两式确定平面上的夜歌矩形，这就是样本空间，其面积为。这时针与平行线相交（记为事件A）的充要条件是由这个不等式表示的区域是图1.2中的阴影部分。图1.1a/2图1.2由于针是向平面任意投掷的，所以由等可能性知这是一个概率的问题。由此得如果l,a为已知，则以的值代入上式即可计算得P（A）之值。反之如果已知P(A)的值，则也可以利用上式去求，而关于P(A)的值，可以从实验中获得频率去近似它：即投掷掷其中针与平行线相交n次，则频率n/N可作为P(A)的估计值，于是由可得历史上有一些学者曾亲自做过这个实验，下表记录了他们的实验结果实验者 年份 针长 投掷次数 相交次数 的近似值Wolf 1850 0.8 5000 2532 3.1596Smith 1855 0.6 3204 1218.5 3.1554DeMorgan.c 1860 1.0 600 382.5 3.137Fox 1884 0.75 1030 489 3.1595Lazzerini 1901 0.83 3408 1808 3.1415929Reina 1925 0.5419 2520 859 3.1795历史上有一些学者曾亲自做过这个实验，下表记录了他们的实验结果可以采用MATLAB软件进行模拟实验，即用MATLAB编写程序来进行“蒲丰投针实验”。MATLAB编程 clear (n)clear(a)clear(x)clear(f)clear (y)clear (m)disp(本程序用来进行投针实验的演示，a代表两线间的宽度，针的长度l=a/2，n代表实验次数);a=input(请输入a：);n=input(请输入n：);x=unifrnd(0,a/2,n,1);f=unifrnd(0,pi,n,1);y=x0.25*a*sin(f);m=sum(y);PI=vpa(a*n/(a*m)实验数据(部分程序截屏见后)anPI第一次3100003.123048第二次3100003.177629第三次31000003.140703第四次31000003.168266第五次310000003.139954第六次310000003.138160第七次3100000003.140658第八次3100000003.141319第九次3200000003.141348第十次3300000003.1421684、结束语从上述数据分析可知，随着模拟次数的越来越多，PI的值逐渐稳定在值附近，即越来越趋近于，故蒲丰投针实验确