高尔顿钉板R语言实验.docxVIP

高尔顿钉板试验【实验目的】1、加强对正态分布的理解2、了解独立同分布的中心极限定理3、掌握R在计算机模拟中的应用【实验要求】1、了解R程序文件的建立和运行，理解循环等控制语句的应用。2、了解R的程序设计，掌握用R处理实际问题的能力。【实验内容】高尔顿钉板试验，这个试验是英国科学家高尔顿设计的，具体如下：自板上端放一个小球，任其自由下落。在其下落过程中，当小球碰到钉子时从左边落下的概率为p，从右边落下的概率为1-p，碰到下一排钉子又是如此，最后落到底板中的某一格子，因此任意放入一球，则此球落入哪个格子事先难以确定（设横排共有m=20排钉子，每一排钉子等距排列，下一排每个钉子恰好在上一排两相邻钉子中间）。（1）分别取p=0.15,0.5,0.85，自板上端放入n个小球，取n=5000，观察n个小球落下后呈现的曲线（直方图）。（2）固定p=0.3，分别取n=1000,10000,100000，观察小球落下后呈现的曲线的变化。【实验思路】令k表示某一个小球在第k次碰到钉子后向左或向右落下这一随机现象相联系的随机变量（=1表示向右落下，=-1表示向左落下），令n=,其中k（k=1,2,…,n）相互独立，则n表示这个小球第n次碰钉后的位置。【实验步骤】实验需要在1和-1之间随机取20次数，并计算其和为横坐标，随机5000次，统计落在各坐标的点的个数做为纵坐标画出图像。在1和-1之间随机取20次数，这相当于随机生成数问题，可以参考二项分布的随机生成，由于二项分布随机生成的为0和1，而本实验想要的是1和-1，则在利用二项生成一组01随机数后，统计其中0的个数作为我所需要的-1的个数。将求和后的数据做为一个向量，将其做成直方图即可。【实验程序及结果】⑴d-NAfor(i in 1:5000){ a-rbinom(20,1,0.15) b-sum(a) d-c(d,(b-10))}hist(d)⑵d-NAfor(i in 1:5000){ a-rbinom(20,1,0.15) b-sum(a) d-c(d,(b-10))}hist(d)⑶d-NAfor(i in 1:5000){ a-rbinom(20,1,0.85) b-sum(a) d-c(d,(b-10))}hist(d)⒉d-NAfor(i in 1:1000){ a-rbinom(20,1,0.3) b-sum(a) d-c(d,(b-10))}hist(d)d-NAfor(i in 1:10000){ a-rbinom(20,1,0.3) b-sum(a) d-c(d,(b-10))}hist(d)d-NAfor(i in 1:100000){ a-rbinom(20,1,0.3) b-sum(a) d-c(d,(b-10))}hist(d)【实验结论】当取定小球数时，概率为0.5时整体图像大致为正态分布图，当概率小于0.5时图像最高点向左偏移，大于0.5时向右便宜。当概率去定时，随着小球数目的增多，图像和正态分布图的拟合程度越来越高，但当小球数超过10000时，变化不明显。