六路高斯光 经过透镜聚焦 matlab代码

以下是一个简单的 Matlab 代码，可以用于计算六路高斯光经过透镜聚焦后的光强分布：

% 定义六路高斯光的参数
w0 = 1; % 光斑半径
zR = pi * w0^2 / 2; % 焦深
k = 2 * pi / 1064e-9; % 波矢量

% 定义透镜的参数
f = 10e-3; % 焦距
D = 10e-3; % 孔径直径

% 定义计算区域
L = 10e-3; % 计算区域大小
N = 512; % 离散点数
x = linspace(-L/2, L/2, N);
y = x;
[X, Y] = meshgrid(x, y);

% 计算光波复振幅
u = zeros(N);
for n = 1:6
    u = u + exp(-(X.^2 + Y.^2) / w0^2) .* exp(-1i * k * (X * cos((n-1)*pi/3) + Y * sin((n-1)*pi/3)));
end

% 计算透镜传输函数
kx = 2*pi/L * (-N/2:N/2-1);
ky = kx;
[KX, KY] = meshgrid(kx, ky);
H = exp(-1i * pi * lambda * f * (KX.^2 + KY.^2));

% 计算透镜聚焦后的光强分布
u_f = ifft2(fft2(u) .* H);
I = abs(u_f).^2;

% 显示结果
figure;
subplot(1,2,1);
imagesc(x, y, abs(u).^2);
xlabel('x (m)');
ylabel('y (m)');
title('初始光强分布');
axis square;
subplot(1,2,2);
imagesc(x, y, I);
xlabel('x (m)');
ylabel('y (m)');
title('透镜聚焦后的光强分布');
axis square;

在这个代码中，我们首先定义了六路高斯光的参数，包括光斑半径、焦深和波矢量。然后定义了透镜的参数，包括焦距和孔径直径。接着定义了计算区域的大小和离散点数，并生成了网格点坐标矩阵。然后计算了光波复振幅，并通过傅里叶变换计算了透镜传输函数。最后通过傅里叶变换和逆傅里叶变换计算了透镜聚焦后的光强分布，并将结果显示出来。

需要注意的是，这个代码只是一个简单的示例，实际的光学系统可能会更加复杂，需要根据具体情况进行修改。