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. 說明默認的窗口長度不是最佳窗口長度.