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

有会相位恢复算法的吗？能帮我看看下面的代码错哪儿了吗？不胜感激 Q=dw 6  
;B,nzx(L  
2d centered Fourier transform，居中的傅里叶变换 r{bg
TG  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% &9S8al 8"  
% Citation for this code and algorithm:
xq8}6Q  
% Tatiana Latychevskaia and Hans-Werner Fink p|xs|O6{  
% "Practical algorithms for simulation and reconstruction of digital in-line holograms", {(8U8f<'=y  
% Appl. Optics 54, 2424 - 2434 (2015) A&xab  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 't||F1X~J  
% The code is written by Tatiana Latychevskaia, 2002 v/QEu^C  
% The version of Matlab for this code is R2010b )zo#1$C-  
^U R-#WaQ  
function [out] = FT2Dc(u0); &Xj{:s#  
oUnq"]  
[Nx, Ny] = size(1i,1j); kq-mr  
#i'C  
f1 = zeros(Nx,Ny); 7[(Lrx.pM  
L{4),65  
for ii = 1:Nx 3U :YA&K(  
    for jj = 1:Ny ^\xCqVk_R  
         f1(ii,jj) = exp(1i*pi*(ii + jj)); "it`X B.  
    end ZJiuj!  
end GLbc/qs  
{o Q(<&Aw  
FT = fft2(f1.*in); 4f_ZY5=  
T)
$6H}[c  
out = f1.*FT; 9;;]q?*  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

求大神指教