用matlab编写一段代码：用波长0.1mm的单色平面波照明一个焦距为20mm的薄透镜，设透镜平面为边长为10mm的方形，在透镜前放置一个边长为1mm的正方形小孔，计算此时透镜焦平面上的光强分布

以下是一段Matlab代码，用于计算上述情况下透镜焦平面上的光强分布：

% 光学参数定义
wavelength = 0.1; % 波长（mm）
f = 20; % 焦距（mm）
L = 10; % 透镜平面边长（mm）
D = 1; % 入射光小孔边长（mm）

% 网格参数定义
N = 1000; % 网格点数
dx = L/N; % 网格尺寸
x = linspace(-L/2, L/2, N); % 网格

% 定义透镜传输函数
k = 2*pi/wavelength;
[X,Y] = meshgrid(x);
R = sqrt(X.^2 + Y.^2);
T = exp(1i*k/2/f*(R.^2));
T(abs(R) > D/2) = 0; % 加上小孔

% 定义入射光波前
E_in = ones(N,N);

% 计算透镜焦平面上的光强分布
E_out = fftshift(fft2(ifftshift(E_in.*T)));
I_out = abs(E_out).^2;

% 绘制结果
figure;
imagesc(x, x, I_out);
axis equal tight;
colormap('gray');
colorbar;
title('透镜焦平面上的光强分布');
xlabel('x (mm)');
ylabel('y (mm)');

这段代码首先定义了光学参数（波长、焦距、透镜平面边长、入射光小孔边长），然后定义了网格参数和网格。接下来，根据透镜传输函数的定义，计算了透镜传输函数，并在传输函数上加上了小孔。然后定义了入射光波前，并使用传输函数计算了透镜焦平面上的光强分布。最后，使用Matlab的图形函数绘制了结果。

希望这段代码能对您有所帮助。如果您有任何其他问题，请随时问我。