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

有会相位恢复算法的吗？能帮我看看下面的代码错哪儿了吗？不胜感激 {j!jm5  
,2R7AHk  
2d centered Fourier transform，居中的傅里叶变换 8F;f&&L"y  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% |RA|nu   
% Citation for this code and algorithm: keMfK]9  
% Tatiana Latychevskaia and Hans-Werner Fink B'kV.3t  
% "Practical algorithms for simulation and reconstruction of digital in-line holograms", ylo/]pVs  
% Appl. Optics 54, 2424 - 2434 (2015) XP|qY1  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% [l7 G9T}/[  
% The code is written by Tatiana Latychevskaia, 2002 &{5v[:$  
% The version of Matlab for this code is R2010b l)m]<EX  
RaBq@r*(
 
function [out] = FT2Dc(u0); MB:VACCr  
VOY#Y*)g  
[Nx, Ny] = size(1i,1j); `-J$7)d@  
^G*zFqa+`  
f1 = zeros(Nx,Ny); v1m'p:7uGB  
itpljh  
for ii = 1:Nx -<jd/ 5  
    for jj = 1:Ny @8;0p  
         f1(ii,jj) = exp(1i*pi*(ii + jj)); y. A]un1  
    end ,o{|W9  
end D#pZN,'  
r@.3.Q  
FT = fft2(f1.*in); | 
WN9&  
yE6EoC^  
out = f1.*FT; B/n/bi8T  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

求大神指教