有会相位恢复算法的吗?能帮我看看下面的代码错哪儿了吗?不胜感激 6i\b&
[w~1e)D
2d centered Fourier transform,居中的傅里叶变换 femAVx}go
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% g;bkVq
% Citation for this code and algorithm: X;p,Wq#D'
% Tatiana Latychevskaia and Hans-Werner Fink ?Zoq|Q+
% "Practical algorithms for simulation and reconstruction of digital in-line holograms", A\.M/)Qo
% Appl. Optics 54, 2424 - 2434 (2015) YKUs>tQ!
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% ]}G(@9
% The code is written by Tatiana Latychevskaia, 2002 rYq8OZLi
% The version of Matlab for this code is R2010b D1lHq/
37!}8
function [out] = FT2Dc(u0); YrgwR
<
v_ ?}
[Nx, Ny] = size(1i,1j); diqG8KaK
M>_S%V4a
f1 = zeros(Nx,Ny); Z!|nc.
OY`G _=6!N
for ii = 1:Nx e v?Hz8Q;(
for jj = 1:Ny JUU0Tx:`9)
f1(ii,jj) = exp(1i*pi*(ii + jj)); e[n T'e
end z/rN+ ,
end S3j/(BG
XUUS N
FT = fft2(f1.*in); Q.6pmaXrb
f+$/gz
out = f1.*FT; g[@]OsX
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%