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

有会相位恢复算法的吗？能帮我看看下面的代码错哪儿了吗？不胜感激  f4Xk,1Is  
PcA2/!a  
2d centered Fourier transform，居中的傅里叶变换 i.3=!6z  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Y8m1M-#w  
% Citation for this code and algorithm: d/BM&r  
% Tatiana Latychevskaia and Hans-Werner Fink BQ=PW|[  
% "Practical algorithms for simulation and reconstruction of digital in-line holograms", E^vJ@O  
% Appl. Optics 54, 2424 - 2434 (2015) CDP U\ZG  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 3QG7C{  
% The code is written by Tatiana Latychevskaia, 2002 r[4tPk  
% The version of Matlab for this code is R2010b \6R,Nq
 
`%+Wz0(K  
function [out] = FT2Dc(u0); Xc?&_\. +  
RloK,bg  
[Nx, Ny] = size(1i,1j); $wo?!gt  
GT0Of~?f  
f1 = zeros(Nx,Ny); y+"X~7EX  
fA+,TEB~d  
for ii = 1:Nx \QYFAa  
    for jj = 1:Ny ~]nSSD)\  
         f1(ii,jj) = exp(1i*pi*(ii + jj)); ti I.W  
    end !HnXXVW  
end B"!l2  
R)QC
)U  
FT = fft2(f1.*in); .P[ _<8  
S/9DtXQ  
out = f1.*FT; -'t)=YJ  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

求大神指教