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

有会相位恢复算法的吗？能帮我看看下面的代码错哪儿了吗？不胜感激 =z=$S]qN  
_A8x{[$  
2d centered Fourier transform，居中的傅里叶变换 qTG
Ei  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 3jH8pO^  
% Citation for this code and algorithm: z x7fRd$  
% Tatiana Latychevskaia and Hans-Werner Fink 
'h?;i2[  
% "Practical algorithms for simulation and reconstruction of digital in-line holograms", P^1+;dL,D  
% Appl. Optics 54, 2424 - 2434 (2015) s'4S,  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 6$d3Ap@Gl  
% The code is written by Tatiana Latychevskaia, 2002 m' LRP:9v  
% The version of Matlab for this code is R2010b FIB 9W@oao  
uk8vecj  
function [out] = FT2Dc(u0); yNwYP%"y  
(A6-9g>  
[Nx, Ny] = size(1i,1j); yAge2m]<B  
Tug}P K   
f1 = zeros(Nx,Ny); 5<?O S &B  
d[b(+sHp a  
for ii = 1:Nx >&1
MD}  
    for jj = 1:Ny 7y4!K$c$  
         f1(ii,jj) = exp(1i*pi*(ii + jj)); Anpx%NVo  
    end v=(L>
gg  
end 3c#CEuu  
INm2
1MS$  
FT = fft2(f1.*in); i[gq8%  
' 9K4A'2[  
out = f1.*FT; ERql^Yr  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

求大神指教