有会相位恢复算法的吗？能帮我看看下面的代码错哪儿了吗？不胜感激 xoA\^AA
QA+qFP
2d centered Fourier transform，居中的傅里叶变换 Vdjf
F&q
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 66,?f<b
% Citation for this code and algorithm: 9@nX 6\ ,
% Tatiana Latychevskaia and Hans-Werner Fink [Z` q7ddd^
% "Practical algorithms for simulation and reconstruction of digital in-line holograms", APsd^ J
% Appl. Optics 54, 2424 - 2434 (2015) w(]Q `
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 9\0
% The code is written by Tatiana Latychevskaia, 2002 % TyR8
%
% The version of Matlab for this code is R2010b cP`o?:
1;i[H[hNY
function [out] = FT2Dc(u0); LBk1Qw}-
P`v%<
9~
[Nx, Ny] = size(1i,1j); ^!3Sz1
5lu620o
f1 = zeros(Nx,Ny); ^ D/:[
l
Ztq_* Fl
for ii = 1:Nx B-
N
for jj = 1:Ny q*SX.A>YR
f1(ii,jj) = exp(1i*pi*(ii + jj)); ]6v6&YV
end ~Yv "=
end }Gqx2 )H
'
^^]Or
FT = fft2(f1.*in); l}&egq
DC
M~t S
*
out = f1.*FT; Q7gBxp
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
