菲涅尔衍射仿真matlab程序

% 本程序用于计算宽度为0.2mm的狭缝的菲涅尔衍射
clear all;close all;
lamda = 0.650e-3;   % 波长
k = 2*pi/lamda;     % 波数
d = 500;            % 传播距离
N = 4096;           % 采样数
HP=30;
no=2.0;
ho=8e-5;
x = linspace（-11N）;   % 衍射屏和观察屏
slit = zeros（1N）;      % 狭缝
slit（abs（x）<=1） = 1;  % 宽度为0.2mm的狭缝
for point=（N/2-HP/2）:（N/2+HP/2-1）
    slit（point）=exp（1i*k*（no-1）*ho）;
end 
for pitch=1:N/4/HP-1
    for HighPhase=1:HP
        slit（N/2+HP/2*3+（pitch-1）*2*HP+（