以下是一个用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库进行图像绘制。通过调整参数，你可以自由地改变波长、衍射距离、光强幅值、矩孔大小等参数，从而观察到不同的双矩孔衍射图案。