[求助]相位恢復(fù)算法 [復(fù)制鏈接]

有會(huì)相位恢復(fù)算法的嗎？能幫我看看下面的代碼錯(cuò)哪兒了嗎？不勝感激 8 F2|  
T/c<23i  
2d centered Fourier transform，居中的傅里葉變換 iJv48#'ii  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Xt{*N-v\  
% Citation for this code and algorithm: FVB;\'/  
% Tatiana Latychevskaia and Hans-Werner Fink ;uqx@sx ;  
% "Practical algorithms for simulation and reconstruction of digital in-line holograms", Uz608u  
% Appl. Optics 54, 2424 - 2434 (2015) zf.-I  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% ~"lJ'&J}  
% The code is written by Tatiana Latychevskaia, 2002 h6%[q x<  
% The version of Matlab for this code is R2010b 'q>2t}KG  
ExSO|g]%
 
function [out] = FT2Dc(u0); ? b[n|^wS  
2oZ9laJO  
[Nx, Ny] = size(1i,1j); e8h,,:l3j  
T*I?9d{k  
f1 = zeros(Nx,Ny); E