MATLAB夫琅禾费衍射矩形孔代码，还附带有3D的图

下面是一个用MATLAB实现夫琅禾费衍射的矩形孔代码，并生成对应的3D图像：

% 定义常数
lambda = 632.8e-9; % 波长（单位：米）
a = 10e-3; % 孔径尺寸（单位：米）
z = 1; % 距离（单位：米）
k = 2*pi/lambda; % 波数

% 定义坐标系
N = 1000; % 离散点数
x = linspace(-a/2, a/2, N);
y = linspace(-a/2, a/2, N);
[X,Y] = meshgrid(x, y);

% 计算衍射场
r = sqrt(X.^2 + Y.^2 + z^2);
theta = k*r;
field = (exp(1i*theta) / (1i*lambda*z)) * (a^2/2) * (abs(besselj(1,k*a*X./r))./r).^2;

% 绘制夫琅禾费衍射图像
figure;
subplot(1,2,1);
imagesc(x, y, abs(field).^2);
colormap('gray');
title('Fraunhofer Diffraction');
xlabel('X (m)');
ylabel('Y (m)');

% 绘制夫琅禾费衍射3D图像
subplot(1,2,2);
mesh(x, y, abs(field).^2);
colormap('gray');
title('Fraunhofer Diffraction (3D)');
xlabel('X (m)');
ylabel('Y (m)');
zlabel('Intensity');

这段代码与之前的代码类似，只是在最后增加了绘制3D图像的部分。使用`mesh`函数可以生成夫琅禾费衍射的3D图像，其中`x`、`y`是坐标，`abs(field).^2`是衍射场的强度。通过调整视角，可以更清楚地观察夫琅禾费衍射的形态。