2021-05-27

clear

%模拟光栅，光栅常数d=10mm，缝宽为3mm，缝数为10。
f=zeros(100,100);          
for i=0:2;                 
    f(8+i:10:98+i,:)=1;    
end                                          
d=0.01;

y=fft2(f);      %对光栅进行二维傅里叶变换，得到光栅场分布的空间频谱。
lamda=5.00e-7;  %波长           
k=2*pi/lamda;  %波数              
z=2*d^2/lamda; %Talbot距离              
H=exp(j*k*z); %菲涅尔衍射传递函数              
r=H*y;       %观察平面场分布的频谱               
f0=ifft2(r);    %二维傅里叶逆变换，得到观察平面场分布的复振幅分布。           
subplot(1,2,1),imshow(f)    
subplot(1,2,2),imshow(f0)   
 

D-15