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

有会相位恢复算法的吗？能帮我看看下面的代码错哪儿了吗？不胜感激 bj 8pqw|;  
w^ofH-R/  
2d centered Fourier transform，居中的傅里叶变换 f<nK;  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% =_g#I  
% Citation for this code and algorithm: `,/5skeJ  
% Tatiana Latychevskaia and Hans-Werner Fink Y}s6
__  
% "Practical algorithms for simulation and reconstruction of digital in-line holograms", /L@o.[H  
% Appl. Optics 54, 2424 - 2434 (2015) 3XQe? 2:<  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -e_TJA  
% The code is written by Tatiana Latychevskaia, 2002 fkX86  
% The version of Matlab for this code is R2010b $]}K;  
}Y!s:w#  
function [out] = FT2Dc(u0); .m>Qlh  
O'#;Ge/,  
[Nx, Ny] = size(1i,1j); 8\:>;XG6f  
?[>Y@w
e  
f1 = zeros(Nx,Ny); matm>3n  
T"C.>G'[B  
for ii = 1:Nx 5WxN
H}{  
    for jj = 1:Ny w2/
3[VZ}l  
         f1(ii,jj) = exp(1i*pi*(ii + jj)); fO^s4gWTg  
    end /38I(0  
end }U5$~,*p  
qTwl\dcncC  
FT = fft2(f1.*in); @(E6P;+{  
F`(;@LO  
out = f1.*FT; \T<F#a  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

求大神指教