菲涅耳衍射的S-FFT计算.docVIP

菲涅耳衍射的S-FFT计算

a=512;

b=512;

forn=1:a

form=1:b

if(n106|n406|m106|m406)

A(n,m)=0;

else

A(n,m)=1;

end

end

end

X=A;

z0=1500%----衍射距离(mm),可按需要修改

z0=input(衍射距离z0=?(mm));

h=0.532e-3;%----波长(mm),可按需要修改

k=2*pi/h;

[M,N]=size(X);

N=min(M,N);

U0=double(X);%初始场复振幅,修改为实际光场复振幅可进行实际衍射场的计算

L0=sqrt(h*z0*N)%FFT计算时同时满足振幅及相位取样条件的物光场宽度

L0=input(初始场宽度L0=?(mm));

Uh=[0:N-1]-N/2;Vh=[0:N-1]-N/2;

[mh,nh]=meshgrid(Uh,Vh);

figstr=strcat(初始图像宽度=,num2str(L0),mm);

figure(1),imshow(X,[]),colormap(gray);ylabel(衍射计算);xlabel(figstr);title(S-FFT方法计算衍射);

%---------------菲涅耳衍射的S-FFT计算起始

n=1:N;

x=-L0/2+L0/N*(n-1);

y=x;

[yy,xx]=meshgrid(y,x);

Fresnel=exp(i*k/2/z0*(xx.^2+yy.^2));

f2=U0.*Fresnel;%S-FFT计算菲涅耳衍射时的傅里叶变换函数

Uf=fft2(f2,N,N);%对N*N点的离散函数f2作FFT计算

Uf=fftshift(Uf);%将FFT计算结果进行整序

L=h*z0*N/L0;%FFT计算后观测屏的物理宽度

x=-L/2+L/N*(n-1);

y=x;

[yy,xx]=meshgrid(y,x);

phase=exp(i*k*z0)/(i*h*z0)*exp(i*k/2/z0*(xx.^2+yy.^2));%菲涅耳衍射积分号前方的相位因子

Uf=Uf.*phase;

T=L0/N;%空域取样间隔

Uf=Uf*T*T;%二维离散变换量值补偿

%---------------菲涅耳衍射的S-FFT计算结束

If=Uf.*conj(Uf);%形成衍射场强度分布

figstr=strcat(衍射场图像宽度=,num2str(L),mm);

%-figure,imshow(A);xlabel(重建Y);ylabel(figstr);

figure(2),imshow(abs(Uf),[]),colormap(gray);ylabel(衍射计算);xlabel(figstr);title(S-FFT计算衍射);

1、衍射距离z0=0且初始场宽度L0=4mm时

2、衍射距离z0=100mm且初始场宽度L0=4mm时

3、衍射距离z0=1000mm且初始场宽度L0=4mm时

4、衍射距离z0=20000mm且初始场宽度L0=