傅里叶变换
傅里叶变换函数
T2F是一个傅里叶变换函数,输入时间样本t与幅度样本st,输出频域频率样本f以及幅频特性sf。
function[f,sf]=T2F(t,st)
dt=t(2)-t(1);
T=t(end);
df=1/T;
N=length(st);
f=-N/2*df:df:N/2*df-df;
sf=fft(st);
sf=T/N*fftshift(sf);
注意:
- sf是一个复数,所以输出的sf要画图需要使用abs();
- fftshift()的原理及用法我觉得这篇文章解释地很好。
傅里叶反变换
傅里叶反变换函数
function[t,st]=F2T(f,sf)
df=f(2)-f(1);
Fmx=(f(end)-f(1)+df); %s