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