使用MATLAB生成菲尼尔透镜

以下是使用MATLAB生成菲涅耳透镜的示例代码：

% 定义透镜参数
r1 = 25; % 内半径
r2 = 50; % 外半径
f = 100; % 焦距
m = 20; % 阶数

% 定义网格参数
N = 500; % 网格大小
x = linspace(-r2, r2, N);
y = linspace(-r2, r2, N);
[X, Y] = meshgrid(x, y);
R = sqrt(X.^2 + Y.^2);

% 计算透镜相位
phase = zeros(size(R));
for n = 1:m
    phase = phase + (-1)^(n+1) * (2*pi*f*R).^n / (factorial(n) * (r2-r1)^n) .* (R >= r1 & R <= r2);
end

% 绘制透镜相位
figure;
imagesc(phase);
axis equal off;
colormap gray;

% 计算透镜衍射图像
lambda = 0.5; % 波长
k = 2*pi/lambda; % 波数
field = exp(1i*k*phase);

% 绘制透镜衍射图像
figure;
imagesc(abs(field).^2);
axis equal off;
colormap gray;

这段代码首先定义了透镜的参数，包括内外半径、焦距和阶数。然后定义了网格参数，并利用网格计算了透镜的相位。接着计算了透镜的衍射图像，并将结果绘制出来。可以根据需要修改透镜参数和网格参数，得到不同形状和性质的菲涅耳透镜。