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

有会相位恢复算法的吗？能帮我看看下面的代码错哪儿了吗？不胜感激 B#EF/\5  
o~es>
;  
2d centered Fourier transform，居中的傅里叶变换 9Cf^Q3)5o  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% #Kn7 xn[  
% Citation for this code and algorithm: hh:)"<[  
% Tatiana Latychevskaia and Hans-Werner Fink ax7 M  
% "Practical algorithms for simulation and reconstruction of digital in-line holograms", IFpmf0;^  
% Appl. Optics 54, 2424 - 2434 (2015) ]*zF#Voc  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% QfwGf,0p  
% The code is written by Tatiana Latychevskaia, 2002 3
]}W  
% The version of Matlab for this code is R2010b  pLyX9C  
pR\etXeLd  
function [out] = FT2Dc(u0); dF1Bo  
fCEd :Kr  
[Nx, Ny] = size(1i,1j); 8eOQRC33  
ZTt%7K"L  
f1 = zeros(Nx,Ny); gfr``z=>O  
@X\nY</E#M  
for ii = 1:Nx |C-B=XE;3  
    for jj = 1:Ny eTrGFe!8w  
         f1(ii,jj) = exp(1i*pi*(ii + jj)); @}fnR(fS  
    end $V5Ol6@2  
end {MYlW0)~  
BAt2m-  
FT = fft2(f1.*in); uc7Y8iO  
WQ"ZQ  
out = f1.*FT; )]qFI"B7  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

求大神指教