www.pudn.com > work.rar > fresnell.m

%菲涅耳积分傅立叶变换方法 
function[f1,dx1,dy1,x1,y1]=fresnell(f0,M,N,dx0,dy0,z,lambda) 
k=2*pi/lambda; 
x0=ones(N,1)*[-M/2:M/2-1]*dx0;y0=[-N/2:N/2-1]'*ones(1,M)*dy0; 
g=f0.*exp(j*0.5*k*(x0.^2+y0.^2)/z); 
G=fftshift(fft2(fftshift(g))); 
du=1./(M*dx0);dv=1./(N*dy0); 
dx1=lambda*z*du;dy1=lambda*z*dv; 
x1=ones(N,1)*[-M/2:M/2-1]*dx1;y1=[-N/2:N/2-1]'*ones(1,M)*dy1; 
f1=G.*exp(i*0.5*k*(x1.^2+y1.^2)/z); 
f1=f1./(i*lambda*z);