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

有会相位恢复算法的吗？能帮我看看下面的代码错哪儿了吗？不胜感激 ?GGBDql  
95*=&d  
2d centered Fourier transform，居中的傅里叶变换 @!B%ynrG  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \N.Bx
  
% Citation for this code and algorithm: Tz{-L%*#  
% Tatiana Latychevskaia and Hans-Werner Fink xd!GRJ<I  
% "Practical algorithms for simulation and reconstruction of digital in-line holograms", A[H"(E#k  
% Appl. Optics 54, 2424 - 2434 (2015) \iAs  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% s}?QA cC  
% The code is written by Tatiana Latychevskaia, 2002 0>yuBgh  
% The version of Matlab for this code is R2010b i.''\  
ms#|Yl1/|  
function [out] = FT2Dc(u0); n
8o(>?Kw  
,p6o "-  
[Nx, Ny] = size(1i,1j); u(g9-O  
\'+P5,  
f1 = zeros(Nx,Ny); Ex5LhRe>=  
5$c*r$t_RK  
for ii = 1:Nx -~]]%VJP|  
    for jj = 1:Ny <h*$bx]9 +  
         f1(ii,jj) = exp(1i*pi*(ii + jj)); lz(}N7SLa  
    end A5,(P$@k  
end gCaxZ~o  
aA-
s{af  
FT = fft2(f1.*in); R!2E`^{Wl  
S{UEV7d:n0  
out = f1.*FT; <,Fj}T-  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

求大神指教