有会相位恢复算法的吗？能帮我看看下面的代码错哪儿了吗？不胜感激 g0;6}n
4`RZ&w;1H2
2d centered Fourier transform，居中的傅里叶变换 X .K*</(g
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 8Vt' X2
% Citation for this code and algorithm: RgM=g8}M
% Tatiana Latychevskaia and Hans-Werner Fink u'Hh||La"
% "Practical algorithms for simulation and reconstruction of digital in-line holograms", mOLz(0
% Appl. Optics 54, 2424 - 2434 (2015) +#X+QG
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% ZR{YpLFQ
% The code is written by Tatiana Latychevskaia, 2002 Y2g% {keo
% The version of Matlab for this code is R2010b vn@sPT
* $1F|G
function [out] = FT2Dc(u0); 8e`HXU(A
6W;`}'ap
[Nx, Ny] = size(1i,1j); Z
6][9o
v`6vc)>8
f1 = zeros(Nx,Ny); YXdo&'Q<qX
2IkyC`
for ii = 1:Nx "`Xbi/i
for jj = 1:Ny C!^;%VQ}d
f1(ii,jj) = exp(1i*pi*(ii + jj)); 9u3P>a~b
end C}9|e?R[Rz
end Upg8t'%{op
Az>r}*F Gr
FT = fft2(f1.*in); X @)'E9g5:
$8 [JL \
out = f1.*FT; r,L`@A=v
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
