以下是一个用Python实现的高中物理波动演示程序，使用`matplotlib`库来绘制波的传播和振动方向。该程序允许用户控制波的传播方向、传播速度以及观察固定点的振动方向。 ```python import numpy as np import matplotlib.pyplot as plt from matplotlib.animation import FuncAnimation # 参数设置 fig, ax = plt.subplots() x = np.linspace(0, 10, 500) # 波的空间范围 t = 0 # 初始时间 wave_speed = 1 # 波速 (可调整) direction = 1 # 波传播方向：1为正向，-1为反向 # 固定点位置 fixed_point = 3.0 # 观察的固定点位置 # 波函数 y = A * sin(kx - ωt) A = 1 # 振幅 k = 2 * np.pi # 波数 omega = wave_speed * k # 角频率 line, = ax.plot(x, np.sin(k * x - omega * t)) # 初始波形 point, = ax.plot([], [], 'ro') # 固定点的振动显示 text = ax.text(0.02, 0.95, '', transform=ax.transAxes) # 显示时间信息 def init(): ax.set_xlim(0, 10) ax.set_ylim(-1.5, 1.5) ax.set_xlabel('Position') ax.set_ylabel('Amplitude') ax.set_title('Wave Propagation and Vibration Direction') return line, point, text def update(frame): global t t += 0.1 # 更新时间 y = A * np.sin(k * x - direction * omega * t) # 更新波形 line.set_data(x, y) # 更新固定点的振动方向 fixed_y = A * np.sin(k * fixed_point - direction * omega * t) point.set_data(fixed_point, fixed_y) # 显示当前时间 text.set_text(f'Time = {t:.2f} s') return line, point, text ani = FuncAnimation(fig, update, frames=range(100), init_func=init, blit=True, interval=50) plt.show() ``` ### 上述代码解释： 1. **波函数**：我们使用了标准的波动方程 `y = A * sin(kx - ωt)`，其中： - `A` 是振幅。 - `k` 是波数（与波长有关）。 - `ω` 是角频率（与周期有关）。 - `t` 是时间。 - 方向由参数 `direction` 控制（正向或反向传播）。 2. **动画更新**：通过 `FuncAnimation` 动态更新波形和固定点的位置。 - 每次调用 `update` 函数时，时间 `t` 增加 0.1 秒，并重新计算波形和固定点的振动值。 3. **固定点振动方向**：在特定位置（如 `fixed_point = 3.0`），显示该点随时间的振动情况。 4. **交互性**：可以通过修改变量 `wave_speed` 和 `direction` 来控制波速和传播方向。 --- ###