假设数据是定期抽样的，因此所有数据的抽样间隔都是相同的，我通常会这样做:
y =.．.；%信号(矢量或矩阵)
t =.．.；%时间矢量
Fs = 1/mean(diff(t));采样频率
Fn = Fs/2;奈奎斯特频率
Ym = y -均值(y);%从所有列中减去列的含义(消除D-C偏移效果)
L = numel(t);信号长度%
f = fft(ym)/L;%傅里叶变换(缩放长度)
Fv = linspace(0,1, fix(L/2)+1)*Fn;%频率矢量
Iv = 1:numel(Fv);%索引向量
数字
情节(阵线、abs (FTy(四世:))* 2)
网格
它计算并绘制了一个单侧傅里叶变换。你将需要供应
“t”
而且
“y”
．代码应该相当健壮。
．
