高尔顿板
从对称高尔顿板到逐层不对称概率:观察弹珠自由下落、碰钉反弹与落槽,比较实际频率和理论分布。
游戏规则
- 设置每层可独立设置左移概率,右移概率自动互补。
- 观察单颗查看碰撞,批量比较分布,支持暂停和继续。
- 比较对称、整体偏置、逐层异质三种情形均可探索。
- 统计只统计已落槽弹珠,改变参数后重新实验。
落点分布 · 共同概率纵轴
浅蓝柱:实际频率 钢蓝虚线:理论概率;落槽区深蓝实线:实际频率平滑曲线
实验结果怎么看
这张表与上方实验实时联动,用同一行直接比较理论预测与已经落槽的实际结果。
| 槽位 | 理论概率 | 期望颗数 | 实际颗数 | 实际频率 | 绝对误差 |
|---|
期望颗数已完成弹珠数
理论概率;绝对误差
实际频率
理论概率|,单位为百分点。尚无完成样本时显示“—”。
1.从碰钉到槽位:先把物理过程翻译成随机变量
设钉板有 层。第
层向左的概率记为
,向右的概率就是
。把“向右”记为
,把“向左”记为
,那么
是一次伯努利试验。
每向右一次,最终位置就比“全向左”的路径多向右一格。因此,从左到右把槽位编号为 0、1、…、 后,最终槽位
本模型假设:同一颗弹珠在不同层的选择相互独立;不同弹珠之间也相互独立;一轮实验中各层概率保持不变。这些是假设,不是所有真实钉板都会自动满足的物理定律。
2.为什么对称板中间高、两边低:多条路径会汇入同一槽
层一共有
条左右路径,却只有 n+1 个槽。最左槽要求“每层都向左”,最右槽要求“每层都向右”,各自只有一条路径;中间槽往往可由许多条路径到达。在对称板上,每条路径等可能,所以中间槽累积的概率更大。
注意:当各层左右概率不同时,各条路径未必等可能,路径多的槽也未必概率最大。我们可以把到达同一槽位的路径概率逐条相加;程序采用下一节的逐层递推,结果相同但计算更高效。
3.对称、不对称,以及它们各自对应的分布
| 各层设置 | 最终槽位 | 分布形状 |
|---|---|---|
| 所有 | 二项分布 | 关于 |
| 所有 | 二项分布 | 可整体偏向一侧 |
| 各层 | 泊松二项分布 | 形状由每一层共同决定 |
这里 表示左移概率,而
统计右移次数,所以二项分布中的“成功概率”是
。当所有层的
都等于
时:
如果各层概率不同,不能把它们先求平均再代入二项分布。这样做虽然能保留均值,却通常会改变方差和各槽概率。名称相近也要注意:泊松二项分布不等于泊松分布。
4.不用枚举所有路径:逐层更新槽位概率
令 表示“经过前
层后,恰好向右
次”的概率。开始时还没经过任何层,所以
。
经过第 层后要处在
:
- 可以先有
次右移,这一层再向左,概率贡献为
;
- 也可以先有
次右移,这一层再向右,概率贡献为
。
两种情况互斥,因此相加:
两层例子:左移概率依次为 70% 和 40%。于是 0 号槽为 ;1 号槽包含“左→右”和“右→左”两条路径,概率为
;2 号槽为
。三者相加正好为 1。
5.均值定位置,方差定分散程度
因为 是各层
选择之和,所以
均值告诉我们落点中心大约在哪里:向右概率整体增大,平均槽位就右移。方差告诉我们分布有多分散:某层的 接近
时最不确定,对方差贡献较大;接近 0 或 1 时几乎确定,对方差贡献较小。
上面的数值会随层数和各层概率实时更新。改变这些参数会清空本轮样本,避免把来自不同理论分布的数据混在一起。
6.实验结果怎么看:频率会靠近概率,但不会整齐地靠近
只统计已经完成下落的 颗弹珠。第
槽的实际颗数记为
,则实际频率
;理论上预计约有
颗落入该槽。
| 槽位 | 理论概率 | 期望颗数 | 实际颗数 | 实际频率 | 绝对误差 |
|---|
绝对误差是 ,表中以百分点表示:例如 30% 与 28% 相差 2 个百分点。槽位计数满足
;图中的柱和理论虚线共用同一概率纵轴。平滑曲线只是在视觉上连接离散频率,不是连续概率密度。
因此,“投得多”只会让频率更稳定,并不会把任意偏置分布变成对称钟形。若每层选择都完全确定( 只能取 0 或 1),所有弹珠会落入同一槽,方差为 0,更谈不上正态近似。