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

有会相位恢复算法的吗？能帮我看看下面的代码错哪儿了吗？不胜感激 FFX-kS  
67/JsL  
2d centered Fourier transform，居中的傅里叶变换 l{<+V)  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 15dhr]8E  
% Citation for this code and algorithm: Ro3C(aRx  
% Tatiana Latychevskaia and Hans-Werner Fink 9oBK(Sf@^  
% "Practical algorithms for simulation and reconstruction of digital in-line holograms", ~A^E_  
% Appl. Optics 54, 2424 - 2434 (2015) ~SUrbRaY>  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% '0q.zzv|_  
% The code is written by Tatiana Latychevskaia, 2002 "g27|e?y  
% The version of Matlab for this code is R2010b 36"-cGNr{  
4'u +%6+__  
function [out] = FT2Dc(u0); IlN: NS  
0%J0.USkM7  
[Nx, Ny] = size(1i,1j); VF[$hs
  
ZD!?mR+-  
f1 = zeros(Nx,Ny); *g7BR`Bt]z  
l@:Tw.+
/9  
for ii = 1:Nx vDI$ QUMD6  
    for jj = 1:Ny 'kU5  
         f1(ii,jj) = exp(1i*pi*(ii + jj)); ov;1=M~RF  
    end ayF+2(vch)  
end I<E~=  
NzS(,F  
FT = fft2(f1.*in); ]M3V]
m  
D!7-(3R  
out = f1.*FT; >s dT=6v  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

求大神指教