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

有会相位恢复算法的吗？能帮我看看下面的代码错哪儿了吗？不胜感激 %;#^l+UB  
Nn>'^KZNG  
2d centered Fourier transform，居中的傅里叶变换 F~ Lx|)0M  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% ~>9_(L  
% Citation for this code and algorithm: t6v/sZ{F  
% Tatiana Latychevskaia and Hans-Werner Fink tabT0  
% "Practical algorithms for simulation and reconstruction of digital in-line holograms", S/;bU:  
% Appl. Optics 54, 2424 - 2434 (2015) -OSa>-bzNx  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% o1nURJ!  
% The code is written by Tatiana Latychevskaia, 2002 O-)-YVU  
% The version of Matlab for this code is R2010b IK*07h/!  
r$)w7Gk<  
function [out] = FT2Dc(u0); n~VD uKn9  
1+?N#Fh  
[Nx, Ny] = size(1i,1j); ~  T>U  
OtFh,}E  
f1 = zeros(Nx,Ny); pW4 cX  
`est|C '+  
for ii = 1:Nx !!Z?[rj  
    for jj = 1:Ny  Q1@A2+ c  
         f1(ii,jj) = exp(1i*pi*(ii + jj)); 1}~(Yj@f%  
    end =B.F;40  
end ftH:r_"O#  
 TX  
FT = fft2(f1.*in); ]qhPd_$?D'  
+SJd@y@fR  
out = f1.*FT; jt+iv*2N>  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

求大神指教