matlab使用經驗模式分解emd 對信號進行去噪

原文鏈接 : http://tecdat.cn/?p=2567

 

對於這個例子，考慮由具有明顯頻率變化的正弦波組成的非平穩連續信號。手提鑽的振動或煙花聲是非平穩連續信號的例子。

以采樣頻率加載非平穩信號數據fs，並可視化混合正弦信號。

1.  
   load（'sinusoidalSignalExampleData.mat'，'X'，'fs'）;
2.  
    
3.  
   xlabel（'Time（s）'）;

觀察到混合信號包含具有不同幅度和頻率值的正弦波。

為了創建希爾伯特譜圖，您需要信號的IMF。執行經驗模式分解以計算信號的固有模式函數和殘差。由於信號不平滑，請指定' pchip'作為Interpolation方法。

[imf，residual，info] = emd（X，'Interpolation'，'pchip'）;

1.  
   目前的IMF | #Sift Iter | 相對Tol | 停止標准命中
2.  
   1 | 2 | 0.026352 | SiftMaxRelativeTolerance
3.  
   2 | 2 | 0.0039573 | SiftMaxRelativeTolerance
4.  
   3 | 1 | 0.024838 | SiftMaxRelativeTolerance
5.  
   4 | 2 | 0.05929 | SiftMaxRelativeTolerance
6.  
   5 | 2 | 0.11317 | SiftMaxRelativeTolerance
7.  
   6 | 2 | 0.12599 | SiftMaxRelativeTolerance
8.  
   7 | 2 | 0.13802 | SiftMaxRelativeTolerance
9.  
   8 | 3 | 0.15937 | SiftMaxRelativeTolerance
10.  
    9 | 2 | 0.15923 | SiftMaxRelativeTolerance
11.  
    分解停止是因為殘留信號的極值數小於'MaxNumExtrema'。

在命令窗口中生成的表指示每個生成的IMF的篩選迭代次數，相對容差和篩選停止標准。此信息也包含在info。您可以通過指定Display為隱藏表0。

HHT（IMF，FS）;

頻率對時間圖是一個稀疏圖，其中垂直顏色條表示IMF中每個點的瞬時能量。該圖表示從原始混合信號分解的每個分量的瞬時頻譜。從該圖中可以觀察到三個IMF，其頻率在1s處有明顯變化。

可視化信號的殘余和內在模式功能

 

對於這個例子，考慮由具有明顯頻率變化的正弦波組成的非平穩連續信號。手提鑽的振動或煙花聲是非平穩連續信號的例子。

加載非平穩信號數據，並可視化混合正弦信號。

load（'sinusoidalSignalExampleData.mat'，'X'，'fs'）; 
 

觀察到混合信號包含具有不同幅度和頻率值的正弦波。

執行經驗模式分解以繪制固有模式函數和信號殘差。由於信號不平滑，請指定' pchip'作為Interpolation方法。

emd（X，'Interpolation'，'pchip'）;

1.  
   目前的IMF | #Sift Iter | 相對Tol | 停止標准命中
2.  
   1 | 2 | 0.026352 | SiftMaxRelativeTolerance
3.  
   2 | 2 | 0.0039573 | SiftMaxRelativeTolerance
4.  
   3 | 1 | 0.024838 | SiftMaxRelativeTolerance
5.  
   4 | 2 | 0.05929 | SiftMaxRelativeTolerance
6.  
   5 | 2 | 0.11317 | SiftMaxRelativeTolerance
7.  
   6 | 2 | 0.12599 | SiftMaxRelativeTolerance
8.  
   7 | 2 | 0.13802 | SiftMaxRelativeTolerance
9.  
   8 | 3 | 0.15937 | SiftMaxRelativeTolerance
10.  
    9 | 2 | 0.15923 | SiftMaxRelativeTolerance
11.  
    分解停止是因為殘留信號的極值數小於'MaxNumExtrema'。

生成具有原始信號，前3個IMF和殘差的交互式圖。在命令窗口中生成的表指示每個生成的IMF的篩選迭代次數，相對容差和篩選停止標准。您可以通過指定Display為隱藏表0。

右鍵單擊圖中的空白區域以打開IMF選擇器窗口。使用IMF選擇器有選擇地查看生成的IMF，原始信號和殘差。

 

 

從列表中選擇要顯示的IMF。選擇是否在圖上顯示原始信號和殘差。

 

 

選定的IMF現在顯示在圖上。

 

 

使用該圖可視化從原始信號中分解的各個組件以及殘差。請注意，殘差是根據IMF總數計算的，並且不會根據IMF選擇器窗口中選擇的IMF進行更改。

非常感謝您閱讀本文，有任何問題請在下方留言！