有会相位恢复算法的吗?能帮我看看下面的代码错哪儿了吗?不胜感激 jZk dTiI
P N(<=v&E
2d centered Fourier transform,居中的傅里叶变换 T"e"?JSRJ
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 79uAsI2-Y
% Citation for this code and algorithm: Jq:Wt+a
% Tatiana Latychevskaia and Hans-Werner Fink QsPZ dC
% "Practical algorithms for simulation and reconstruction of digital in-line holograms", "td ,YVK
% Appl. Optics 54, 2424 - 2434 (2015) B$G9#G6pZ
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% r.' cjUs
% The code is written by Tatiana Latychevskaia, 2002 {>d\
% The version of Matlab for this code is R2010b ;?{OX
Yq:+.UU
function [out] = FT2Dc(u0);
8#|PJc
t;
@T~%
[Nx, Ny] = size(1i,1j); v807)JwS
pYj}
f1 = zeros(Nx,Ny); wc ;^C?PX
Me-H'Mp~
for ii = 1:Nx -e.ygiK.`S
for jj = 1:Ny 'kekJ.wJ;
f1(ii,jj) = exp(1i*pi*(ii + jj)); }q)dXFL=I#
end qSQsY:]j0
end 7Pr5`#x#
(X|lK.W y
FT = fft2(f1.*in); {WuUzq`
(o8?j^ -v
out = f1.*FT; 51puR8AG>
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%