fs = 2000;
t = 0:1/fs:2-1/fs;
q = chirp(t-2,4,1/2,6,'quadratic',100,'convex').*exp(-4*(t-1).^2);
plot(t,q)
emd を使用して、固有モード関数 (IMF) と残差を可視化します。
emd(q)
信号の IMF を計算します。'Display' の名前と値のペアを使用して、各 IMF のふるい分け反復の数、相対許容誤差、およびふるい分け停止基準を示すテーブルを出力します。
imf = emd(q,'Display',1);
Current IMF | #Sift Iter | Relative Tol | Stop Criterion Hit
1 | 2 | 0.0063952 | SiftMaxRelativeTolerance
2 | 2 | 0.1007 | SiftMaxRelativeTolerance
3 | 2 | 0.01189 | SiftMaxRelativeTolerance
4 | 2 | 0.0075124 | SiftMaxRelativeTolerance
Decomposition stopped because the number of extrema in the residual signal is less than the 'MaxNumExtrema' value.
計算された IMF を使用して、二次チャープのヒルベルト スペクトルをプロットします。周波数範囲を 0 Hz ~ 20 Hz に制限します。
信号の最初の 3 つの IMF を計算します。'Display' の名前と値のペアを使用して、各 IMF のふるい分け反復の数、相対許容誤差、およびふるい分け停止基準を示すテーブルを出力します。
imf = emd(whale,'MaxNumIMF',3,'Display',1);
Current IMF | #Sift Iter | Relative Tol | Stop Criterion Hit
1 | 1 | 0.13523 | SiftMaxRelativeTolerance
2 | 2 | 0.030198 | SiftMaxRelativeTolerance
3 | 2 | 0.01908 | SiftMaxRelativeTolerance
Decomposition stopped because maximum number of intrinsic mode functions was extracted.
計算された IMF を使用して、信号のヒルベルト スペクトルをプロットします。周波数範囲を 0 Hz ~ 1400 Hz に制限します。
[1] Huang, Norden E, and Samuel S P Shen. Hilbert–Huang Transform and Its Applications. 2nd ed. Vol. 16. Interdisciplinary Mathematical Sciences. WORLD SCIENTIFIC, 2014. https://doi.org/10.1142/8804.
[2] Huang, Norden E., Zhaohua Wu, Steven R. Long, Kenneth C. Arnold, Xianyao Chen, and Karin Blank. “ON INSTANTANEOUS FREQUENCY.” Advances in Adaptive Data Analysis 01, no. 02 (April 2009): 177–229. https://doi.org/10.1142/S1793536909000096.