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

有会相位恢复算法的吗？能帮我看看下面的代码错哪儿了吗？不胜感激 Fr2kbQTg;  
gY%OhYtF2  
2d centered Fourier transform，居中的傅里叶变换 3?  };
 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% @"0uM?_)-  
% Citation for this code and algorithm: /Bh>  
% Tatiana Latychevskaia and Hans-Werner Fink Ej{+U  
% "Practical algorithms for simulation and reconstruction of digital in-line holograms", 4?v$<=#21*  
% Appl. Optics 54, 2424 - 2434 (2015) e) ]RA?bF  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% D/cg7  
% The code is written by Tatiana Latychevskaia, 2002 blUY.{NN3  
% The version of Matlab for this code is R2010b m[W/j/$A+x  
ovB
=Zm  
function [out] = FT2Dc(u0); L,WkJe3  
%uj[`  
[Nx, Ny] = size(1i,1j); 9@Q&B+!  

B#QL
M^  
f1 = zeros(Nx,Ny); :Jm!=U%'Z  
x\6];SXX  
for ii = 1:Nx <o~t$TH  
    for jj = 1:Ny 4`$5 _} j!  
         f1(ii,jj) = exp(1i*pi*(ii + jj)); b'&LBT7  
    end `|&\e_"DE  
end WFRsSp2  
}M^_Z#|,  
FT = fft2(f1.*in); xC'mPcU8  
zf`5>h|  
out = f1.*FT; aXX,Zu^  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

求大神指教