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

有会相位恢复算法的吗？能帮我看看下面的代码错哪儿了吗？不胜感激 #s-iy+/1oN  
gk0.zz([  
2d centered Fourier transform，居中的傅里叶变换 bzpFbfb  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% knp>m,w  
% Citation for this code and algorithm: A;XOT6jv?  
% Tatiana Latychevskaia and Hans-Werner Fink p zw8T  
% "Practical algorithms for simulation and reconstruction of digital in-line holograms", Nh?|RE0t  
% Appl. Optics 54, 2424 - 2434 (2015) u
YI@9U  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% DbI!l`Vn4  
% The code is written by Tatiana Latychevskaia, 2002 F0Jx(  
% The version of Matlab for this code is R2010b ;F:fM!l=  
hS [SRa'.  
function [out] = FT2Dc(u0); XKOUQ
c4!R  
Njc%_&r  
[Nx, Ny] = size(1i,1j); lJ;J~>  
p&p.Q^"ok  
f1 = zeros(Nx,Ny); (46 {r}_O  
b,H[I!. %  
for ii = 1:Nx %V!iQzL1  
    for jj = 1:Ny "d/s5sP|S  
         f1(ii,jj) = exp(1i*pi*(ii + jj)); '=E9En#@  
    end 5/><$06rq  
end F t&+vS  
gwkb!#A  
FT = fft2(f1.*in); \`.F\Z  
E`\8TqO  
out = f1.*FT; M+>`sj  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

求大神指教