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

有会相位恢复算法的吗？能帮我看看下面的代码错哪儿了吗？不胜感激 X>I)~z}9#  
KZbR3mi,  
2d centered Fourier transform，居中的傅里叶变换 x-'~Bu  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 68iV/7  
% Citation for this code and algorithm: d$Y7u  
% Tatiana Latychevskaia and Hans-Werner Fink j.i#*tN//  
% "Practical algorithms for simulation and reconstruction of digital in-line holograms", m,R Dr  
% Appl. Optics 54, 2424 - 2434 (2015) IhiGP {  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% @[6,6:h|  
% The code is written by Tatiana Latychevskaia, 2002 &5%dhc4&!&  
% The version of Matlab for this code is R2010b &6j<ca  
H2r8,|XL  
function [out] = FT2Dc(u0); #n  
phYDs9-K  
[Nx, Ny] = size(1i,1j); OI0B:()  
j+8TlVur  
f1 = zeros(Nx,Ny); gFPi7 o1  
}K(o9$V ^!  
for ii = 1:Nx 7dX/bzUVz8  
    for jj = 1:Ny oA tsUF+a  
         f1(ii,jj) = exp(1i*pi*(ii + jj));  +<.\5+  
    end q+a.G2S  
end kLS(w??T  
e#/kNHl  
FT = fft2(f1.*in); _}mK!_`  
4.K'
\S  
out = f1.*FT; v8I{XU@%  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

求大神指教