clfftf.m

时间域有限差分法（TD-FDM）源代码

function [FT,omega] = clfftf(ft,N,Ts)

FT1 = Ts*fft(ft,N);
FT = fftshift(FT1);
fs = 1/N/Ts;
omega = 2*pi*linspace(-N/2,N/2-1,N)*fs;

end