方孔菲涅尔衍射代码

tic;
lambda=589.3e-9;
N=2000; %设定输入面上的采样点阵N×N
L=1e-3;    %设定孔径大小
z=200*L;    %设定观察屏至衍射屏距离
xi=linspace（-LLN）;
xom=1.5*L;   %（重要参数）输出面采样区域边界;
         %-xom≤x≤xom-xom≤y≤xom
M=401;      %设定输出面上的采样点阵M×M
xo=linspace（-xomxomM）;
Xo=meshgrid（xoxo）;
Ix=zeros（si