MATLAB R2018a中信号处理相关的新函数 : pkurtosis与kurtogram

pkurtosis: Spectral kurtosis from signal or spectrogram

用来判断是否为平稳的(stationary)白噪信号, 如果是的话, 取值在各个频率上都是较低的(在置信区间内).

如果某些频率上在置信区间外面的话, 说明在对应的频率下, 存在非平稳信号或者是非白噪信号.

基于STFT的结果, 进一步计算频域的峰度.

 

 

一个300Hz到400Hz的非平稳信号 + 白噪声做实验.

 

很清晰看到, 检测到300Hz到400Hz是非平稳或非白噪

 

kurtogram: Visualize spectral kurtosis

相关函数为pkurtosis.

 

pkurtosis使用了一个窗口长度来计算STFT(进一步计算频域峰度) , 而kurtogram使用了很多个窗口长度来计算STFT(进一步计算频域峰度)

 

kurtogram可以用来观察不同窗口长度下, 频域峰度的变化, 还可以用来确定最佳窗口长度.

最佳指的是频域峰度最大化.

 

 

使用最佳的窗口长度来计算pkurtosis

 

可以看出, 最佳窗口长度下, 频域峰度最大到了10, 而使用默认的窗口长度下, 频域峰度最大才到6. 说明默认的窗口长度不是最佳窗口长度.