### 夫琅禾费衍射动画的Python实现 以下是一个使用Python实现夫琅禾费衍射动画的代码示例。该代码基于`numpy`和`matplotlib`库，能够生成动态的衍射图案并以动画形式展示。 ```python import numpy as np import matplotlib.pyplot as plt from matplotlib.animation import FuncAnimation # 参数设置 wavelength = 500e-9 # 光波长 (米) distance = 1.0 # 衍射屏到观察屏的距离 (米) slit_width = 0.0001 # 缝隙宽度 (米) screen_size = 0.01 # 观察屏尺寸 (米) num_points = 500 # 屏幕上采样点的数量 steps = 100 # 动画帧数 # 计算屏幕上的位置 x = np.linspace(-screen_size / 2, screen_size / 2, num_points) # 夫琅禾费衍射强度公式 def fraunhofer_intensity(x, slit_width, wavelength, distance): k = 2 * np.pi / wavelength beta = k * slit_width * x / (2 * distance) intensity = (np.sinc(beta / np.pi) ** 2) * (slit_width ** 2) return intensity # 创建图形对象 fig, ax = plt.subplots() line, = ax.plot([], [], lw=2) ax.set_xlim(-screen_size / 2, screen_size / 2) ax.set_ylim(0, 1.2 * fraunhofer_intensity(0, slit_width, wavelength, distance)) ax.set_xlabel("Position on Screen (m)") ax.set_ylabel("Intensity (a.u.)") ax.set_title("Fraunhofer Diffraction Animation") # 初始化函数 def init(): line.set_data([], []) return line, # 更新函数 def update(frame): global slit_width slit_width += frame * 1e-6 # 动态调整缝隙宽度 intensity = fraunhofer_intensity(x, slit_width, wavelength, distance) line.set_data(x, intensity) return line, # 创建动画对象 ani = FuncAnimation(fig, update, frames=np.linspace(0, 1e-4, steps), init_func=init, blit=True, interval=50) # 显示动画 plt.show() ``` #### 代码说明 1. **参数设置**：定义了光波长、衍射屏到观察屏的距离、缝隙宽度等物理参数[^1]。 2. **夫琅禾费衍射公式**：通过数学公式计算不同位置上的衍射强度，其中使用了`np.sinc`函数来表示正弦平方项[^2]。 3. **动画生成**：利用`FuncAnimation`类将衍射强度随时间的变化以动画形式展示，动态调整缝隙宽度以观察衍射图案的变化[^3]。 --- ###