有会相位恢复算法的吗?能帮我看看下面的代码错哪儿了吗?不胜感激 |X~T</{8i
PrDvRWM
2d centered Fourier transform,居中的傅里叶变换 q<,?:g$k
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% TnPd pynP
% Citation for this code and algorithm: Ds&)0Iwf
% Tatiana Latychevskaia and Hans-Werner Fink .6(Bf$E
% "Practical algorithms for simulation and reconstruction of digital in-line holograms", s@^GjA[6+
% Appl. Optics 54, 2424 - 2434 (2015) W{=>c/
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% }P-9\*hlm
% The code is written by Tatiana Latychevskaia, 2002 ;G;vpl
% The version of Matlab for this code is R2010b 1a/@eqF''
Ndx.SOj
function [out] = FT2Dc(u0); g#9KG
HTf7r-
[Nx, Ny] = size(1i,1j); sJt&`k Z
~"
}t8`vP1
f1 = zeros(Nx,Ny); <ynmA
JAmv 7GL'6
for ii = 1:Nx k{y@&QNj
for jj = 1:Ny ToDNBt.u{+
f1(ii,jj) = exp(1i*pi*(ii + jj)); P[#V{%f*5
end Zhz.8W
end Zo-s_6uC
qJ5Y}/r
FT = fft2(f1.*in); vRRi"bo
x8b w#
out = f1.*FT; .~ZNlI {K
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%