有会相位恢复算法的吗?能帮我看看下面的代码错哪儿了吗?不胜感激 }NwmZw>_
w9GY/]
2d centered Fourier transform,居中的傅里叶变换 O4f9n
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 0'r}]Mws
% Citation for this code and algorithm: A\rY~$Vr
% Tatiana Latychevskaia and Hans-Werner Fink yC5>k;/6#K
% "Practical algorithms for simulation and reconstruction of digital in-line holograms", 5OGwOZAj52
% Appl. Optics 54, 2424 - 2434 (2015) V.u^;gr3
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% [ _xOz4`%
% The code is written by Tatiana Latychevskaia, 2002 3{co.+
% The version of Matlab for this code is R2010b sq#C|v/
85io%>&0
function [out] = FT2Dc(u0); uGXvP(Pg'
hl**G4z9q
[Nx, Ny] = size(1i,1j); J+ :3==,
XS@iu,uO
f1 = zeros(Nx,Ny); ?^!J:D?
HI D6h!
for ii = 1:Nx >BO!jv!a
for jj = 1:Ny _?{KTgJ G
f1(ii,jj) = exp(1i*pi*(ii + jj)); {)r[?%FMgV
end :%mlsNw
end \)\n5F:Zu
C'A
D[`p
FT = fft2(f1.*in); 8b,Z)"(U3
8cY5:plK
out = f1.*FT; #oEtLb@O
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%