function?Ud=D_FFT（U0zL01L02lambda）
%參數預處理
k=2*pi/lambda;
[MN]=size（U0）;
%計算頻域網格坐標
m=1:M;n=1:N;
t1=-M/L01/2+1/L01*（m-1）;
t2=-N/L02/2+1/L02*（n-1）;
[xy]=meshgrid（t2t1）;
%計算菲涅爾衍射
trans=exp（1i*k*z*（1-lambda^2.*（x.^2+y.^2）/2））;
Uf=fftshift（fft2（U0MN））;
Ud=ifft2（Uf.*transMN）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????312??2016-04-23?19:14??D_FFT.m

?????文件???????1522??2016-04-25?12:48??SLM_IG_simulation.m

-----------?---------??----------?-----??----

?????????????????1834????????????????????2