有会相位恢复算法的吗?能帮我看看下面的代码错哪儿了吗?不胜感激 h!GixN?
qNL~m'
2d centered Fourier transform,居中的傅里叶变换 zxdO3I
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% _,~zy9{,
% Citation for this code and algorithm: DL_\luh
% Tatiana Latychevskaia and Hans-Werner Fink KIJ[ cIw
% "Practical algorithms for simulation and reconstruction of digital in-line holograms", 2jC\yY |PN
% Appl. Optics 54, 2424 - 2434 (2015) ~T9[\nU\
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% |;yb *
% The code is written by Tatiana Latychevskaia, 2002 l si8?91
% The version of Matlab for this code is R2010b SbivW5|61
=1capix 1r
function [out] = FT2Dc(u0); pC8i&_A
hBu=40K
[Nx, Ny] = size(1i,1j); !6C d.fpWL
i-b7
f1 = zeros(Nx,Ny); &pz8vWCk
` 454=3H
for ii = 1:Nx Wz]S+IpY
for jj = 1:Ny {{,%p#/b
f1(ii,jj) = exp(1i*pi*(ii + jj)); ]"6<"1)
end L+L9)8FJ
end a06DeRCej
.YRSd
FT = fft2(f1.*in); C-Z,L#
i$bBN$<b<
out = f1.*FT; LAeX e!y
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%