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)
D15