有会相位恢复算法的吗？能帮我看看下面的代码错哪儿了吗？不胜感激 &xC5Mecb*
O+DYh=m*p
2d centered Fourier transform，居中的傅里叶变换 EUPc+D3
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% a}k5[)et
% Citation for this code and algorithm: oBPm^ob4
% Tatiana Latychevskaia and Hans-Werner Fink 0w2<2grQ
% "Practical algorithms for simulation and reconstruction of digital in-line holograms", HErG%v]nw
% Appl. Optics 54, 2424 - 2434 (2015) *'@T+$3s
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% *MCkezW7{
% The code is written by Tatiana Latychevskaia, 2002 K [-G2
% The version of Matlab for this code is R2010b `Z@qWB<
2tp95E`(O
function [out] = FT2Dc(u0); eN TKX
>/-Bg:
[Nx, Ny] = size(1i,1j); c5eimA%`
2 ) Q/cH\g
f1 = zeros(Nx,Ny); x)kp*^/
~MK% ^5y?
for ii = 1:Nx &x[V<Gq
for jj = 1:Ny -0*z"a9<p8
f1(ii,jj) = exp(1i*pi*(ii + jj)); S]c& T`jx
end jpiBHi]5+
end .Jc<Gg
@ w,O1Xwj
FT = fft2(f1.*in); kbx4I?
Lnltt86
out = f1.*FT; QIcc@PGT9a
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
