菲涅尔衍射matlab仿真

在MATLAB中，可以使用Fresnel衍射公式进行菲涅尔衍射的仿真。下面是一个简单的MATLAB代码示例，用于进行菲涅尔衍射的仿真：

lambda = 632.8e-9;   % 波长（单位：米）
D = 1e-3;   % 光阑尺寸（单位：米）
L = 1;   % 传播距离（单位：米）
N = 1024;   % 采样点数
dx = D/N;   % 采样间隔

x = (-N/2:N/2-1) * dx;   % 采样点位置

% 产生光阑函数
A = double(abs(x)<D/2);   % 方形光阑

% 计算衍射场
U0 = A.*exp(1i*pi/lambda/L*(x.^2));   % 初始场
U1 = fftshift(fft(U0));   % 进行傅里叶变换
U2 = exp(1i*pi*lambda*L*(x.^2));   % 相因子
U3 = ifft(fftshift(U1.*U2));   % 进行逆傅里叶变换

% 绘制结果
figure;
subplot(2,2,1);
plot(x, abs(U0).^2);
title('初始场');
xlabel('位置（米）');
ylabel('强度');

subplot(2,2,2);
plot(x, abs(U1).^2);
title('傅里叶变换后');
xlabel('位置（米）');
ylabel('强度');

subplot(2,2,3);
plot(x, abs(U2).^2);
title('相因子');
xlabel('位置（米）');
ylabel('强度');

subplot(2,2,4);
plot(x, abs(U3).^2);
title('逆傅里叶变换后');
xlabel('位置（米）');
ylabel('强度');

上述代码中，使用方形光阑作为初始场，在傅里叶域中进行变换并添加相因子，然后进行逆傅里叶变换得到衍射场。最后，绘制了初始场、傅里叶变换后、相因子和逆傅里叶变换后的结果。

你可以根据自己的需求修改代码中的参数和光阑函数，以得到不同条件下的菲涅尔衍射仿真结果。