python布朗运动模拟程序

以下是一个简单的Python布朗运动模拟程序的例子：

import numpy as np
import matplotlib.pyplot as plt

# 设置参数
num_particles = 1000  # 粒子数量
num_steps = 1000  # 步数
step_size = 0.1  # 步长

# 初始化粒子位置
particles = np.zeros((num_particles, 2))

# 模拟运动
for i in range(num_steps):
    # 生成随机步长
    steps = np.random.normal(scale=step_size, size=(num_particles, 2))
    # 更新粒子位置
    particles += steps
    # 绘制粒子位置
    plt.clf()
    plt.scatter(particles[:, 0], particles[:, 1], s=1)
    plt.xlim(-50, 50)
    plt.ylim(-50, 50)
    plt.pause(0.01)

plt.show()

该程序使用numpy和matplotlib库来生成随机步长并绘制粒子位置。程序中的参数可以根据需要进行调整，以获得所需的模拟效果。