python来模拟高尔顿钉板实验（带下落过程的）

1.逻辑的思考

 首先，要用turtle把模型画出来

 然后，把turtle变成蓝色的小球，来模拟小球下落的轨迹。

 定义一个draw函数来控制小球的左右移动

 再定义一个box函数来模拟下落小球的堆积

2.代码的实现

模型板块，用循环即可画出来

draw函数里，x作为入参数，来控制turtle的方向和画轨迹的方向。在0,1里面取随机数，0就左走，1就右走。

而box函数里，定义一个9个元素的数组array，初始值都为0，用count来标记小球进了哪个坑（一共九个），小球下落后清零。而数组里面的数则控制画的长方形的高度（即模拟小球的堆积）下落一个小球对应的数就加一。

import numpy as np
import turtle as t
import random
from turtle import Terminator
win = t.Screen()
win.setup(1000,1000)
t.speed(1000000)
for i in range(1,9):
    t.circle(10)
    t.penup()
    t.goto(-15*i,-15*(3**(1/2)*i))
    t.pendown()
for i in range(1,8):
    t.penup()
    t.goto(15 * i, -15 * (3 ** (1 / 2) * i))
    t.pendown()
    t.circle(10)
    t.penup()
    t.pendown()
t.left(90)
t.penup()
t.forward(20)
t.pendown()
t.left(90)
for i in range(1,7):
    t.penup()
    t.forward(30)
    t.pendown()
    t.circle(10)
    t.pendown()
t.left(180)
for i in range(1,6):
    t.penup()
    t.goto(-15*(i-1),-15*(i+1)*(3**(1/2)))
    t.pendown()
    t.circle(10)
    t.penup()
    t.pendown()
for i in range(1,6):
    t.penup()
    t.goto(15*(i-1),-15*(i+1)*(3**(1/2)))
    t.pendown()
    t.circle(10)
    t.penup()
    t.pendown()
for i in range(1,4):
    t.penup()
    t.goto(15*(i-1),-15*(i+3)*(3**(1/2)))
    t.pendown()
    t.circle(10)
    t.penup()
    t.pendown()
for i in range(1,4):
    t.penup()
    t.goto(-15*(i-1),-15*(i+3)*(3**(1/2)))
    t.pendown()
    t.circle(10)
    t.penup()
    t.pendown()
t.penup()
t.goto(0,(-90*(3**(1/2))))
t.pendown()
t.circle(10)
t.penup()
t.right(90)
for i in range(10):
    t.penup()
    t.goto(-135+30*i,-105*(3**(1/2)))
    t.pendown()
    t.forward(200)
t.penup()
t.goto(-135,-105*(3**(1/2)-200))
t.left(90)
t.pendown()
t.forward(9*30)

#变成小球
t.penup()
t.speed(1000)
t.shape("circle")
t.color("blue")
t.shapesize(0.2)
t.goto(0,20.8)
t.pendown()

#随机
count = 0
array=np.zeros(9)
def draw(x):
    t.penup()
    t.seth(0 if x<0 else 180)
    t.pendown()
    r = 18
    t.circle(x*r,116)
    t.penup()
    t.pendown()
#底下的格子
def box():
    t.penup()
    t.goto(135-30*count,-200-105*(3**(1/2)))
    t.pendown()
    t.begin_fill()
    t.fillcolor("black")
    t.seth(90)
    t.fd(0.3*array[8-count])
    t.seth(180)
    t.fd(30)
    t.seth(270)
    t.fd(0.3*(array[8-count]))
    t.seth(360)
    t.fd(30)
    t.end_fill()
try:
    for z in range(100):
        for i in range(1,9):
            x = random.randint(0,1)
            if x == 0:
                draw(1)
                count = count+1

            else:
                draw(-1)
        array[8-count]+=1
        t.seth(270)
        t.forward(200)
        box()
        count = 0
        t.penup()
        t.goto(0,20.8)
        t.pendown()
except Terminator:pass
t.mainloop()