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

有会相位恢复算法的吗？能帮我看看下面的代码错哪儿了吗？不胜感激 {;.q?mj  
E#3KWp#
M  
2d centered Fourier transform，居中的傅里叶变换 qFEGV+  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% zO#{qF+~;  
% Citation for this code and algorithm: q;co53.+P)  
% Tatiana Latychevskaia and Hans-Werner Fink -<rQOPH%  
% "Practical algorithms for simulation and reconstruction of digital in-line holograms", |GPR3%9  
% Appl. Optics 54, 2424 - 2434 (2015) 7w}D2|+  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% {ctEjgiE  
% The code is written by Tatiana Latychevskaia, 2002 ~x<nz/^  
% The version of Matlab for this code is R2010b 52+;j[ ]/O  
gq9D#B  
function [out] = FT2Dc(u0); 6u-@_/O5R3  
/*g9drwaa  
[Nx, Ny] = size(1i,1j); Ki#({~  
jG($:>3a@  
f1 = zeros(Nx,Ny); jDI)iW`P  
=gI;%M\'  
for ii = 1:Nx :eaqUW!Y  
    for jj = 1:Ny >7@kwj-f)  
         f1(ii,jj) = exp(1i*pi*(ii + jj)); $@m)8T  
    end Hya  ";'  
end QF^_4Yn  
~ou1{NS  
FT = fft2(f1.*in); ^/;W;C{4
 
BR8W8nRb  
out = f1.*FT; +<p?i]3CHe  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

求大神指教