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

有会相位恢复算法的吗？能帮我看看下面的代码错哪儿了吗？不胜感激 C4~;yhz  
u%$Zqee  
2d centered Fourier transform，居中的傅里叶变换 3b+d"`Y^S  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +eFFSt  
% Citation for this code and algorithm: ev#;t@^
 
% Tatiana Latychevskaia and Hans-Werner Fink ,!7 H]4Qx  
% "Practical algorithms for simulation and reconstruction of digital in-line holograms", 2\7`/,U6  
% Appl. Optics 54, 2424 - 2434 (2015) d*~ICir7  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% ]cGA~d  
% The code is written by Tatiana Latychevskaia, 2002 z#]Jv!~EPE  
% The version of Matlab for this code is R2010b (fCXxyZrr  
W!6qqi{  
function [out] = FT2Dc(u0); :yeq(oK,  
vfq%H(  
[Nx, Ny] = size(1i,1j); WC.t_"@  
BbgnqzU  
f1 = zeros(Nx,Ny); )0Me?BRp  
<-,gAk)u  
for ii = 1:Nx IFDZfx  
    for jj = 1:Ny Yx eOI#L  
         f1(ii,jj) = exp(1i*pi*(ii + jj));
uoXAQ6k  
    end $D}{]MN.  
end "QWq_R  
UmQ?rS8d  
FT = fft2(f1.*in); )e a:Q?  
{3.r6ZwCn  
out = f1.*FT; +eX@U;J,g  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

求大神指教