%?本程序用于計(jì)算寬度為0.2mm的狹縫的菲涅爾衍射
clear?all;close?all;
lamda?=?0.650e-3;???%?波長(zhǎng)
k?=?2*pi/lamda;?????%?波數(shù)
d?=?500;????????????%?傳播距離
N?=?4096;???????????%?采樣數(shù)
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+（