好的，我可以为您提供一个用Python编写双孔菲涅尔衍射的代码示例。下面是代码： ```python import numpy as np import matplotlib.pyplot as plt # 定义计算函数 def fresnel_diffraction(A, wavelength, distance, size, num_points): x = np.linspace(-size/2, size/2, num_points) y = np.linspace(-size/2, size/2, num_points) X, Y = np.meshgrid(x, y) r1 = np.sqrt((X - distance/2)**2 + Y**2) r2 = np.sqrt((X + distance/2)**2 + Y**2) phase1 = 2 * np.pi * r1 / wavelength phase2 = 2 * np.pi * r2 / wavelength amplitude = A * np.sinc(phase1) + A * np.sinc(phase2) intensity = np.abs(amplitude)**2 return intensity # 参数设置 A = 1.0 # 幅度 wavelength = 0.5 # 波长 distance = 1.0 # 双孔间距 size = 5.0 # 平面尺寸 num_points = 500 # 离散点数 # 计算双孔菲涅尔衍射 intensity = fresnel_diffraction(A, wavelength, distance, size, num_points) # 绘制结果 plt.imshow(intensity, cmap='gray') plt.colorbar() plt.title('Double Slit Fresnel Diffraction') plt.xlabel('x') plt.ylabel('y') plt.show() ``` 在这个代码示例中，我们定义了一个名为`fresnel_diffraction`的函数来计算双孔菲涅尔衍射的强度分布。该函数接受幅度（A）、波长、双孔间距、平面尺寸和离散点数作为输入，并返回计算得到的强度分布。 您可以根据需要调整参数的值来进行实验和观察不同情况下的双孔菲涅尔衍射现象。希望对您有帮助！如果您有任何其他问题，请随时提出。