matlab编程矢量衍射积分的程序

以下是matlab编写的矢量衍射积分程序：

function E = diffraction_integral(lambda, a, b, L, X, Y, Z, theta, phi)
% lambda: 波长
% a, b: 矩形孔的宽度和高度
% L: 矩形孔到屏幕的距离
% X, Y, Z: 屏幕上的点的坐标
% theta, phi: 入射光线的极角和方位角

k = 2 * pi / lambda; % 波数
N = 100; % 数值积分的分段数
delta = a / N; % 积分步长
E = 0;
for m = -N/2:N/2-1
    for n = -N/2:N/2-1
        x = m * delta;
        y = n * delta;
        r = sqrt((X-x)^2+(Y-y)^2+L^2); % 点(X,Y)到矩形孔的距离
        cos_theta = sin(theta) * cos(phi) * (X-x)/r + sin(theta) * sin(phi) * (Y-y)/r + cos(theta) * L/r;
        E = E + exp(1i * k * (X * sin(theta) * cos(phi) + Y * sin(theta) * sin(phi) + Z * cos(theta))) * exp(1i * k * (x * sin(theta) * cos(phi) + y * sin(theta) * sin(phi))) * (cos_theta * exp(-1i * k * r) / r^2) * delta^2;
    end
end
E = E * (k / (2 * pi));

该程序采用数值积分法计算矩形孔的衍射积分，通过对矩形孔进行分段求和，最终得到屏幕上每个点的光强分布。