Fraunhofer三角孔衍射matlab仿真，代码显示给我

以下是一个简单的 MATLAB 代码，可以用于模拟 Fraunhofer 三角孔衍射：

% 定义三角孔的大小和形状
a = 0.05; % 孔的半径
b = 0.5; % 孔到屏的距离
n = 512; % 屏幕像素数目

% 生成坐标网格
xmin = -1;
xmax = 1;
ymin = -1;
ymax = 1;
x = linspace(xmin, xmax, n);
y = linspace(ymin, ymax, n);
[X, Y] = meshgrid(x, y);

% 计算衍射场
k = 2 * pi / 0.5; % 光的波数
r = sqrt(X.^2 + Y.^2); % 到中心点的距离
theta = atan2(Y, X); % 到中心点的角度
H = zeros(n); % 衍射场
for i = 1:n
    for j = 1:n
        if r(i, j) < a
            H(i, j) = 1;
        else
            H(i, j) = 0;
        end
    end
end
H = H .* exp(-1i * k * b * sqrt(1 - (a ./ r).^2)); % 衍射场

% 计算衍射图样
I = abs(fftshift(fft2(H))).^2;

% 显示衍射图样
figure;
imagesc(x, y, I);
colormap(gray);
axis('equal');
axis('off');

运行以上代码后，可以得到一个 Fraunhofer 三角孔衍射的图样。