Spectrogram使用小记

短时傅立叶变换spectrogram最详细的说明当时是看matlab官方的help文件


格式:
spectrogram(x,window,noverlap,nfft)

说明:

x:为输入向量
window:为使用的窗函数,默认为hamming窗,长度为nfft长度
noverlap:为每一段重叠的样本数,一般使用窗函数的50%
nfft:为fft的长度

例子:

产生一个线性调频信号x,采样率1kHz,采样时间为2s。起始频率为100Hz,终止频率为200Hz,时间1s。

t = 0:0.0001:2;
x = chirp(t,100,1,200,'quadratic');

使用spectrogram函数之间绘制频谱图

spectrogram(x,128,120,128,1e3)
  • 使用的是128阶的hamming窗
  • 重叠长度为120
  • 频率分量有[128/2+1]=65
    时间分量有[length(x)-120/(128-120)] = 235
    screen-shot-2016-11-24-at-3-48-35-pm
    在最后添加‘yaxis’可以将X轴变为时间轴,Y轴变为频率
spectrogram(x,128,120,128,1e3,'yaxis')

screen-shot-2016-11-24-at-3-49-14-pm

发表评论

电子邮件地址不会被公开。