用python写一个双矩孔衍射的代码，要求其波长、衍射距离、光强幅值、矩孔大小可调

以下是一个用Python编写的双矩孔衍射的代码示例： ```python import numpy as np import matplotlib.pyplot as plt def double_slit_diffraction(wavelength, distance, intensity, slit_width, slit_distance, screen_width, screen_distance): # 计算衍射角 theta = np.arctan(screen_width / (2 * screen_distance)) # 计算衍射公式中的常数 k = 2 * np.pi / wavelength # 生成屏幕上的坐标 x = np.linspace(-screen_width/2, screen_width/2, 1000) # 计算衍射光强 diffraction_intensity = intensity * np.sinc(k * slit_width * np.sin(theta) * x / distance)**2 \ * np.cos(k * slit_distance * np.sin(theta) * x / distance)**2 # 绘制衍射图像 plt.plot(x, diffraction_intensity) plt.xlabel('Position on Screen') plt.ylabel('Intensity') plt.title('Double Slit Diffraction Pattern') plt.show() # 设置参数 wavelength = 500e-9 # 波长（单位：米） distance = 1 # 衍射距离（单位：米） intensity = 1 # 光强幅值 slit_width = 50e-6 # 矩孔宽度（单位：米） slit_distance = 200e-6 # 矩孔间距（单位：米） screen_width = 0.1 # 屏幕宽度（单位：米） screen_distance = 1 # 屏幕距离（单位：米） # 调用函数绘制衍射图像 double_slit_diffraction(wavelength, distance, intensity, slit_width, slit_distance, screen_width, screen_distance) ``` 在这个代码中，我们使用了numpy库进行数值计算和matplotlib库进行图像绘制。通过调整参数，你可以自由地改变波长、衍射距离、光强幅值、矩孔大小等参数，从而观察到不同的双矩孔衍射图案。