有会相位恢复算法的吗?能帮我看看下面的代码错哪儿了吗?不胜感激 :pqUUZ6x&
k||t<&`Ze
2d centered Fourier transform,居中的傅里叶变换 tg/UtE`V
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% `x9Eo4(/
% Citation for this code and algorithm: M Irx,d
% Tatiana Latychevskaia and Hans-Werner Fink ,RV
qYh(-|
% "Practical algorithms for simulation and reconstruction of digital in-line holograms", }aVZ\PDg
% Appl. Optics 54, 2424 - 2434 (2015) ,_Z(!|
rW
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 5QMra5N k
% The code is written by Tatiana Latychevskaia, 2002 s{Z)<n03
% The version of Matlab for this code is R2010b 5 8bW
}'$PYAf6
function [out] = FT2Dc(u0); ZD]1C~)
EO&Q
[Nx, Ny] = size(1i,1j); <W"W13*j!
^a4z*#IOr
f1 = zeros(Nx,Ny); ^}_Ka //k
dI{DiPho
for ii = 1:Nx lM6pYYEq=
for jj = 1:Ny ;.A}c)b
f1(ii,jj) = exp(1i*pi*(ii + jj)); M3pE$KT0x
end 2*|T)OA`m,
end +WH\,E
%;,fI'M
FT = fft2(f1.*in); FE~D:)Xj'?
fK J-/{|
out = f1.*FT; ;)(g$r^_i
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%