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

有会相位恢复算法的吗？能帮我看看下面的代码错哪儿了吗？不胜感激 [l-zU}u&v  
d=oOMXYa   
2d centered Fourier transform，居中的傅里叶变换 ]1)@.b;QR  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% >, 234ab=d  
% Citation for this code and algorithm: }h
+a8@  
% Tatiana Latychevskaia and Hans-Werner Fink nKu`Ta*fX  
% "Practical algorithms for simulation and reconstruction of digital in-line holograms", @!0j)5%  
% Appl. Optics 54, 2424 - 2434 (2015) 3a!/EP  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% t<ftEJU"'w  
% The code is written by Tatiana Latychevskaia, 2002 R:E6E@T  
% The version of Matlab for this code is R2010b -fL|e/  
6|T{BOW!d  
function [out] = FT2Dc(u0); **! lV]/  
)P|%=laE8  
[Nx, Ny] = size(1i,1j); :{_Or'L  
AA>5h<NM  
f1 = zeros(Nx,Ny); .Z%7+[  
^b-H  
for ii = 1:Nx _q6+]  
    for jj = 1:Ny 2*K0~ b`  
         f1(ii,jj) = exp(1i*pi*(ii + jj)); _\@i&3hkx  
    end e!tgWYN  
end AU87cqq  
*XniF~M  
FT = fft2(f1.*in); m9#u.Q*  
Qi=rhN`  
out = f1.*FT; D.*JG7;=Z  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

求大神指教