高尔顿钉板试验模拟.docx

高尔顿钉板试验模拟(程序)

...这是我2005年12的课程设计中程序的核心部分，写完后自己非常得意，等着老师表扬。等啊等，等待现在也没等到：em16:

现将它献给大家...(若有版权那遵守BSD吧)

注1：程序以前是用Matlab写的现用Java重写

注2：原程序中galton返回值为int[]grid、没有输出结果”部分

publicvoidgalton(intsumOfGrid,intsumOfBall){

int[]grid=newint[sumOfGrid];

intnumber=0; 〃一个小球从顶端落下过程中向右偏移的总次数

intrand; //随机数，取值范围为{0,1}，为0、为1的概率相等

for(intcounter_ball=1;counter_ball=sumOfBall;counter_ball++){

//核心

//(sumOfGrid-1)为钉板的层数

for(inttimes=1;times=(sumOfGrid-1);times++){

rand=(int)(Math.random()*2);

number+=rand;

}

grid[number]++;

number=0;

///核心

}

//输出结果

System.out.println(小球的总数为+sumOfBall+\t格子的个数为+sumOfGrid);

for(intindex=0;indexgrid.length;index++)

System.out.println((index+1)+号格子中的小球数为：\t+grid[index]);

}

}//endofmetodgalton

补充：（谢谢2楼提醒:-D）

高尔顿钉板试验：自板上端放入一小球，任其自由落下.在下落过程中，当小球碰到钉子时，从左边落下与从右边落下的机会相等.碰到下一排钉子也是如此.自板上端放入n（n自行输入）个小球，观察小球落下后呈现曲线并统计小球落入各个格子的频率.

1?■I■.■4

1

?■

I■

.

■4

.

■

.

1■.

.

?

?

.

I■,

.

.

?

.

,?

.

i?

.

?

?

r—

i

.

.

.

?

.

■■

?

*4

4

?.

?

?

.

.

?a

.

■

11

.

*9?

0

(P

o%Io

o

o

o

O

o

%o

o

o

钉子

小球

格子

In

高尔顿钉板试验可见《概率论》（复旦大学李贤平）

当小球数量少时分布无明显特征，当小球数量多时（＞100）分布近似正态分布。（即两边对称:-D）当时为了证明服从正态分布投1千万个小球（计算机模拟:-D）

GaltonBoard

TheGaltonboard,alsoknownasaquincunxorbeanmachine,isadeviceforstatisticalexperimentsnamedafterEnglishscientistSirFrancisGalton.Itconsistsofanuprightboardwithevenlyspacednails(orpegs)drivenintoitsupperhalf,wherethenailsarearrangedinstaggeredorder,andalowerhalfdividedintoanumberofevenly-spacedrectangularslots.Thefrontofthedeviceiscoveredwithaglasscovertoallowviewingofbothnailsandslots.Inthemiddleoftheupperedge,thereisafunnelintowhichballscanbepoured,wherethediameteroftheballsmustbemuchsmallerthanthedistancebetweenthenails.Thefunnelislocatedpreciselyabovethecentralnailofthesecondrowsothateachball,ifperfectlycentered,wouldfallverticallyanddirectlyontotheuppe