[求助]相位恢复算法 [复制链接]

有会相位恢复算法的吗？能帮我看看下面的代码错哪儿了吗？不胜感激 <I?f=[  
iP9Dr<P  
2d centered Fourier transform，居中的傅里叶变换 hsK(09:J
 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% pJo4&Ff  
% Citation for this code and algorithm: @P=n{-pIW  
% Tatiana Latychevskaia and Hans-Werner Fink h9nh9a(2  
% "Practical algorithms for simulation and reconstruction of digital in-line holograms", 0G'v4Vj0'  
% Appl. Optics 54, 2424 - 2434 (2015) ]N1,"W}  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% h,/3}  
% The code is written by Tatiana Latychevskaia, 2002 'I[xZu/8yg  
% The version of Matlab for this code is R2010b mkgL/h*  
S4E@wLi  
function [out] = FT2Dc(u0); pUgas?e&  
~jWpD7px  
[Nx, Ny] = size(1i,1j); mndEB!b  
6Tq2WZ}<'  
f1 = zeros(Nx,Ny); Ghe@m6|D  
ILHn~d IC  
for ii = 1:Nx :19s=0  
    for jj = 1:Ny kWbY&]ZO  
         f1(ii,jj) = exp(1i*pi*(ii + jj)); E*v+@rv  
    end \ov]Rn  
end IJ{VCzi  
1"1ElH  
FT = fft2(f1.*in); On+0@hh  
zHW&i~  
out = f1.*FT; 4LJOT_  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

求大神指教