伽尔顿板程序设计.PDF

伽尔顿板程序设计 刘慧琳 北京师范大学系统科学学院 学号：201821250021 email: liuhl@mail.bnu.edu.cn 一、 简要介绍 伽尔顿板是在一块竖直木板的上部规则地钉上铁钉，木板的下部用竖直隔板 隔成等宽的狭槽，从顶部中央的漏斗形入口处可以投入小球，板前覆盖玻璃使小 球不致落到槽外。小球从入口处投入，在下落过程中将与铁钉发生多次碰撞，最 后落入某一槽中。具体实验模型见图 1。 图 1 伽尔顿板实验 伽尔顿板实验通常是分别多次投入单个小球或者同时投入许多小球，观察比 较小球在各个槽中的分布。实验结果发现：投入单个小球，小球与铁钉碰撞后落 入哪个槽中完全是偶然的或者随机的。大量小球同时投入或单个小球分别多次投 入，最终落入中间部位槽中的小球总是较多，而落入两侧槽中的小球总是较少。 多次重复实验发现各槽中小球数目分布基本不变，但又不是绝对相同。 二、 程序设计 我们可以将实验看作多次投入单个小球，每次在与铁钉碰撞后，小球都完全 随机地选择向左偏还是向右偏。我们在程序中可以设置左偏为取值−1，右偏为 取值1，那么将一个小球每次偏转对应的取值相加的和就代表该小球最终的落点 位置。具体程序编写思路见图 2： 开始 设置小球数N， 设置钉子排数M 确定小球每 一步的偏转 计算小球最终落点位置 绘制所有小球 落点分布图 结束 图 2 程序思路 但是在确定小球每一步的偏转时，有一个需要关注的问题，即小球向左偏和 向右偏的概率是否相等？如果不相等，又取什么值呢？ 对此问题，我们将本程序分为以下几种版本。其中主要差别在于子程序—— 确定小球每一步的偏转。 1、原始版本 当小球左右偏的概率相同，即皆为1/2 时， 程序如下： % step1 设置 N=100000; %小球数 M=10; %钉子排数 % step2 确定小球每一步的选择 step=randsrc(N,M); %1/2 的几率向左 % step3 得到小球最终落点 point=sum(step,2); %行求和，得落点 % step4 画出所有小球落点分布图 len(1,:)=-M:M; for k=-M:M len(2,k+M+1)=length(find(point==k)); end relen=len(:,1:2:end); figure; plot(relen(1,:)./2,relen(2,:)) 得到的一次分布图见图 3： 图 3 原始版本小球分布图 2、改进版本1：可设置偏转概率 当小球左右偏的概率不同，设置小球向左偏的概率为p，则向右偏的概率为 1-p，所得程序如下： % step1 设置 N=100000; %小球数 M=10; %钉子排数 % step2 确定小球每一步的选择 p=0.3; %设定向左偏转的概率 alphabet=[-1,1]; prob=[p,1-p]; for i=1:N for j=1:M step(i,j)=randsrc(1,1,[alphabet;prob]); end end % step3 得到小球最终落点 point=sum(step,2); %行求和，得落点