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

有会相位恢复算法的吗？能帮我看看下面的代码错哪儿了吗？不胜感激 lF!Iu.MM 9  
k'+Mc%pg4E  
2d centered Fourier transform，居中的傅里叶变换 X3rvM8  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 6w^Fee`>]  
% Citation for this code and algorithm: T13Jno  
% Tatiana Latychevskaia and Hans-Werner Fink x)o`w"]al  
% "Practical algorithms for simulation and reconstruction of digital in-line holograms", xGymQ|y84  
% Appl. Optics 54, 2424 - 2434 (2015) JV9Ft,xk  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% A+F@JpV  
% The code is written by Tatiana Latychevskaia, 2002 8VZLwhj  
% The version of Matlab for this code is R2010b lDhuL;9e  
*|k/lI  
function [out] = FT2Dc(u0); p*(]8pDC  
f}F   
[Nx, Ny] = size(1i,1j); a,*~wmg  
2u'h,on?  
f1 = zeros(Nx,Ny); $qj||zA  
?BnjtefIe  
for ii = 1:Nx 4 g^oy^~  
    for jj = 1:Ny ?]u=5gqUU  
         f1(ii,jj) = exp(1i*pi*(ii + jj)); %1VfTr5  
    end -dsE9)&8DX  
end ZtqN8$[6n  
0^rDf L  
FT = fft2(f1.*in); B>W!RyH8o  
E`>u*D$un~  
out = f1.*FT; Qape DU;  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

求大神指教