高尔顿钉板试验模拟程序

1、高尔顿钉板试验模拟（程序）这是我2005年12的课程设计中程序的核心部分，写完后自己非常得意，等着老师表扬。等啊等, 等待现在也没等到：em16:现将它献给大家（若有版权那遵守BSD吧）注1:程序以前是用Matlab写的现用Java重写注2：原程序中galton返回值为int口 grid、没有“输出结果”部分 public void galton(int sumOfGrid, int sumOfBall)int口 grid = new intsumOfGrid;int number = 0;/一个小球从顶端落下过程中向右偏移的总次数int rand ;/随机数，取值范围为0,1,为0、为1的概

2、率相等for( int counter_ball = 1; counter_ball <= sumOfBall; counter_ball+ )/<核心/ (sumOfGrid - 1)为钉板的层数for( int times = 1; times <= ( sumOfGrid - 1 ); times+ ) rand = (int)( Math.random()*2 );number += rand;gridnumber+;number = 0;/</核心>/输出结果System.out.println("小球的总数为"+sumOfBall+

3、"t 格子的个数为"+sumOfGrid );for( int index = 0; index < grid.length; index+ )System.out.println( (index+1)+”号格子中的小球数为：t"+gridindex);/end of metod galton补充：（谢谢2楼提醒：-D ）高尔顿钉板试验：自板上端放入一小球，任其自由落下.在下落过程中，当小球碰到钉子时，从左边落 下与从右边落下的机会相等.碰到下一排钉子也是如此.自板上端放入n（n自行输入）个小球，观察小球 落下后呈现曲线并统计小球落入各个格子的频率 .钉子一

4、小球格子高尔顿钉板试验可见概率论（复旦大学李贤平）当小球数量少时分布无明显特征，当小球数量多时（>100）分布近似正态分布。（即两边对称：-D）当时为了证明服从正态分布投1千万个小球（计算机模拟：-D）Galton BoardThe Galton board, also known as a quincunx or bean machine, is a device for statistical experiments named after English scientist Sir Francis Galton. It consists of an upright board wi

5、th evenly spaced nails (or pegs) driven into its upper half, where the nails are arranged in staggered order, and a lower half divided into a number of evenly-spaced rectangular slots.The front of the device is covered with a glass cover to allow viewing of both nails and slots. In the middle of the

6、 upper edge, there is a funnel into which balls can be poured, where the diameter of the balls must be much smaller than the distance between the nails. The funnel is located precisely above the central nail of the second row so that each ball, if perfectly centered, would fall vertically and direct

7、ly onto the uppermost point of this nail's surface (Kozlov and Mitrofanova 2002). The figure above shows a variant of the board in which only the nails that can potentially be hit by a ball dropped from the funnel are included, leading to a triangular array instead of a rectangular one.Each time

8、 a ball hits one of the nails, it can bounce right (or left) with some probability p (and q = 1 p ).For symmetrically placed nails, balls will bounce left or right with equal probability, so p = q = 1 / 2 . If therows are numbered from 0 to N -1, he path of each falling ball is a Bernoulli trial con

9、sisting of N steps. Each ball crosses the bottom row hitting the nth peg from the left (where0±n三 N -1) if it has taken exactly right turns, which occurs with probability= 尸This process therefore gives rise to a binomial distribution of in the heights of heaps of balls in the lower slots.If the

10、 number of balls is sufficiently large and一 "一 1 : , then according to the weak law of large numbers, the distribution of the heights of the ball heaps will approximate a normal distribution.Some care is needed to obtain these idealized results, however, as the actual distribution of balls depe

11、nds on physical properties of the setup, including the elasticity of the balls (as characterized by their coefficient of restitution), the radius of the nails, and the offsets of the balls over the funnel's opening when they are dropped (Kozlov and Mitrofanova 2002).AXL RRLLLR LLRL LftLL RLLLLLRR RLRL LRLR RRLLLRQL RLLRRRRRPfnbdbilici«sC.