用matlab模拟拉盖尔高斯光束圆孔衍射

以下是使用MATLAB模拟拉盖尔高斯光束通过圆孔进行衍射的示例代码：

% 定义参数
wavelength = 0.5; % 波长
k = 2*pi/wavelength; % 波数
a = 1; % 圆孔半径
f = 10; % 焦距
z = 100; % 衍射距离
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);
theta = atan2(Y, X);
U = sqrt(2/pi) * (wavelength/(pi*a))^0.25 * exp(-R.^2/a^2) .* exp(1i*k*z) .* exp(1i*theta);

% 绘制光强分布
figure;
imagesc(x, y, abs(U).^2);
colormap('hot');
colorbar;
axis square;
title('拉盖尔高斯光束圆孔衍射');

% 相关问题: