有会相位恢复算法的吗?能帮我看看下面的代码错哪儿了吗?不胜感激 Fr2kbQTg;
gY%OhYtF2
2d centered Fourier transform,居中的傅里叶变换 3? };
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% @ "0uM?_)-
% Citation for this code and algorithm: /Bh>
% Tatiana Latychevskaia and Hans-Werner Fink Ej{+U
% "Practical algorithms for simulation and reconstruction of digital in-line holograms", 4?v$<=#21*
% Appl. Optics 54, 2424 - 2434 (2015) e) ]RA?bF
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% D/cg7
% The code is written by Tatiana Latychevskaia, 2002 blUY.{NN3
% The version of Matlab for this code is R2010b m[W/j/$A+x
ovB=Zm
function [out] = FT2Dc(u0); L,WkJe3
%uj[ `
[Nx, Ny] = size(1i,1j); 9@Q&B+!
B#QL M^
f1 = zeros(Nx,Ny); :Jm!=U%'Z
x\6];SXX
for ii = 1:Nx <o~t$TH
for jj = 1:Ny 4`$5
_}
j!
f1(ii,jj) = exp(1i*pi*(ii + jj)); b'&LBT7
end `|&\e_"DE
end WFRsSp2
}M^_Z#|,
FT = fft2(f1.*in); xC'mPcU8
zf`5>h|
out = f1.*FT; aXX,Zu^
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%