### Python 实现伽尔顿板模拟 以下是基于 Python 的伽尔顿板模拟程序，该程序允许用户自定义左右移动的概率以及小球的数量： ```python import random import matplotlib.pyplot as plt def galton_board(num_balls, num_slots, prob_left): """ 模拟伽尔顿板实验。 参数: num_balls (int): 小球总数。 num_slots (int): 插槽数量（即钉子层数加一）。 prob_left (float): 向左偏移的概率，范围为0到1。 返回: list: 记录每个插槽中小球数量的结果列表。 """ if not (0 <= prob_left <= 1): raise ValueError("prob_left 必须在 0 和 1 之间")[^2] results = [0] * num_slots for _ in range(num_balls): position = 0 for level in range(num_slots - 1): # 随机决定向左还是向右 direction = random.random() if direction < prob_left: position += 1 # 向右移动 else: position += 0 # 继续留在原位置 results[position] += 1 return results def plot_results(results): """ 使用柱状图展示伽尔顿板实验结果。 参数: results (list): 结果数据列表。 """ slots = list(range(len(results))) plt.bar(slots, results, color='blue', edgecolor='black') plt.xlabel('插槽编号') plt.ylabel('落入的小球数目') plt.title('伽尔顿板模拟结果') plt.show() if __name__ == "__main__": try: number_of_balls = int(input("请输入小球的总数量: ")) number_of_slots = int(input("请输入插槽数量(通常等于钉子层数加一): ")) probability_left = float(input("请输入向左移动的概率(介于0和1之间): ")) result_distribution = galton_board(number_of_balls, number_of_slots, probability_left) print(f"最终分布情况: {result_distribution}") plot_results(result_distribution) except Exception as e: print(f"发生错误: {e}")[^3] ``` #### 解析 上述代码实现了伽尔顿板的核心逻辑并提供了可视化功能。 - `galton_board` 函数负责执行核心计算过程，接受三个参数：小球数量 (`num_balls`)、插槽数量 (`num_slots`) 和向左移动的概率 (`prob_left`)。如果输入的概率不在合法范围内，则会抛出异常提示用户修正输入值[^2]。 - `plot_results` 是用于绘制直方图的一个辅助函数，它通过调用 Matplotlib 库来直观显示每次试验后的统计结果[^3]。 此脚本还包含了简单的命令行界面让用户可以灵活设置不同的条件来进行多次测试观察不同设定下的效果变化。