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

有会相位恢复算法的吗？能帮我看看下面的代码错哪儿了吗？不胜感激 ~i>'3j0@k  
A$;*O)  
2d centered Fourier transform，居中的傅里叶变换 :z|$K^)7Z  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% uD ;T   
% Citation for this code and algorithm: U-{3HHA  
% Tatiana Latychevskaia and Hans-Werner Fink Kf$6D 79#  
% "Practical algorithms for simulation and reconstruction of digital in-line holograms", jDj=a->e^  
% Appl. Optics 54, 2424 - 2434 (2015) yC!>7@m  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 2{g&9  
% The code is written by Tatiana Latychevskaia, 2002 8. %g&%S  
% The version of Matlab for this code is R2010b yDpv+6(a  
EvA8<o  
function [out] = FT2Dc(u0); (;HO3Z".q$  
VHL[Y  
[Nx, Ny] = size(1i,1j); 
a^VI)  
;B:\e8  
f1 = zeros(Nx,Ny); 7CX5pRNL  
-}8r1jQH;  
for ii = 1:Nx u=(.}  
    for jj = 1:Ny M?['HoRo  
         f1(ii,jj) = exp(1i*pi*(ii + jj)); x3jjtjf  
    end 0zA:
?}  
end %fh-x(4v  
&G3$q,`H  
FT = fft2(f1.*in); %@C$xM"  
D{4]c)>  
out = f1.*FT; z34+1d  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

求大神指教