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

有会相位恢复算法的吗？能帮我看看下面的代码错哪儿了吗？不胜感激 4Smno%jq  
sOa`Tk  
2d centered Fourier transform，居中的傅里叶变换 Wy.2*+5FX0  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% >k<.bEx(A  
% Citation for this code and algorithm: -_[n2\|we)  
% Tatiana Latychevskaia and Hans-Werner Fink tbJB0T|G  
% "Practical algorithms for simulation and reconstruction of digital in-line holograms", 2P`hdg  
% Appl. Optics 54, 2424 - 2434 (2015) J_ ?;On5  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% =SA@3)kHH  
% The code is written by Tatiana Latychevskaia, 2002 x/{  
% The version of Matlab for this code is R2010b yBXdj`bV  
A7>0Pn%D3  
function [out] = FT2Dc(u0); H^_[nL
 
7G.IGXK$  
[Nx, Ny] = size(1i,1j); g05:A0X#  
kLU$8L  
f1 = zeros(Nx,Ny); )?^0<l#s  
7FW!3~3A_  
for ii = 1:Nx
Ytmt+9  
    for jj = 1:Ny 5rRYv~+  
         f1(ii,jj) = exp(1i*pi*(ii + jj)); |k$[+53A  
    end ..UmbJJ.u  
end R!0O[i  
%k_R;/fjW  
FT = fft2(f1.*in); }_u1'  
/zMiy?  
out = f1.*FT; tH,}_Bp  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

求大神指教