### 圆孔衍射的Python模拟 为了实现圆孔衍射的模拟，可以利用菲涅尔-基尔霍夫衍射积分理论来描述光场分布。以下是基于 Python 的一种方法，该方法使用 `numpy` 和 `matplotlib` 库完成数值计算和可视化。 #### 参数说明 - **波长 (`wavelength`)：** 衍射现象中的光源波长。 - **孔径直径 (`diameter`)：** 圆形孔洞的实际物理尺寸。 - **观察距离 (`distance`)：** 孔到屏幕的距离。 这些参数直接影响远场（Fraunhofer）或近场（Fresnel）区域内的衍射图案特性[^1]。 下面是一个完整的代码示例： ```python import numpy as np import matplotlib.pyplot as plt def simulate_circular_aperture_diffraction(wavelength, diameter, distance, screen_size=0.01, resolution=500): """ 模拟圆形孔口的衍射效应 :param wavelength: 波长 (m) :param diameter: 圆孔直径 (m) :param distance: 屏幕与孔之间的距离 (m) :param screen_size: 屏幕大小的一半长度 (m)，默认为 0.01 m :param resolution: 计算分辨率，默认为 500 像素 :return: None """ k = 2 * np.pi / wavelength # 波数 x = np.linspace(-screen_size, screen_size, resolution) # X轴坐标网格 y = np.linspace(-screen_size, screen_size, resolution) # Y轴坐标网格 xx, yy = np.meshgrid(x, y) r = np.sqrt(xx**2 + yy**2) # 距离中心点的径向距离矩阵 q = np.arctan2(r, distance) # 极角 # Fraunhofer衍射强度公式 intensity = ((2 * jv(1, k * diameter / 2 * np.sin(q))) / (k * diameter / 2 * np.sin(q)))**2 # 处理零除错误的情况 intensity[np.isnan(intensity)] = 1 # 可视化结果 plt.figure(figsize=(8, 6)) plt.imshow(intensity, extent=[-screen_size*1e3, screen_size*1e3, -screen_size*1e3, screen_size*1e3], cmap='viridis') plt.colorbar(label="Intensity") plt.title("Circular Aperture Diffraction Pattern", fontsize=14) plt.xlabel("X Position (mm)", fontsize=12) plt.ylabel("Y Position (mm)", fontsize=12) plt.show() # 示例调用 simulate_circular_aperture_diffraction( wavelength=500e-9, # 波长设为500 nm diameter=1e-3, # 圆孔直径设为1 mm distance=1, # 观察屏距离设为1米 screen_size=0.01 # 屏幕范围设置为±1 cm ) ``` 此函数实现了对给定波长、孔径直径以及观测距离下的圆孔衍射图样的计算并绘制其强度分布图像[^2]。 --- ### 关键概念解释 1. **Bessel 函数的应用** 上述程序中使用的 Bessel 函数来源于光学中的 Airy 公式，用于描述理想情况下圆孔产生的衍射斑点形状[^3]。 2. **远场条件假设** 当满足 \( d \gg \lambda L/d \) （其中 \( d \) 是孔径直径，\( L \) 是观察距离），则可认为处于远场区，适用 Fraunhofer 衍射模型[^4]。 ---