有会相位恢复算法的吗?能帮我看看下面的代码错哪儿了吗?不胜感激 /7De.O~H
DZ?>9W{
2d centered Fourier transform,居中的傅里叶变换 Y.sf^}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +V4BJ/H
% Citation for this code and algorithm: 41>Bm*if
% Tatiana Latychevskaia and Hans-Werner Fink 1!/cd;{B
% "Practical algorithms for simulation and reconstruction of digital in-line holograms", 'D#}ce)s#
% Appl. Optics 54, 2424 - 2434 (2015) 0I* ^VGZ
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% #.?DsK_:@
% The code is written by Tatiana Latychevskaia, 2002 H6 ( ~6Bp5
% The version of Matlab for this code is R2010b %?J\P@
KS}Ci-
function [out] = FT2Dc(u0); j )Zi4<./
P1]ucu_y,
[Nx, Ny] = size(1i,1j); pOS:/~I3
"tBdz V
f1 = zeros(Nx,Ny); 7qt<CLJ
%1 <No/
for ii = 1:Nx ?q1&(g]qO
for jj = 1:Ny %gSqc
}v*
f1(ii,jj) = exp(1i*pi*(ii + jj)); us\@n"
end 6*cm
end 3}= .7qm
w~afQA>
FT = fft2(f1.*in); =H.l/'/Z
V7%G?
out = f1.*FT; YiB]}/
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%