1. 程式人生 > >matlab使用經驗模式分解emd 對訊號進行去噪

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

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

以取樣頻率載入非平穩訊號資料fs,並可視化混合正弦訊號。

load('sinusoidalSignalExampleData.mat''X''fs'); 
 
xlabel('Time(s)');

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

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

[imf,residual,info] = emd(X,'Interpolation'
'pchip');
目前的IMF | #Sift Iter | 相對Tol | 停止標準命中   
      1 | 2 | 0.026352 | SiftMaxRelativeTolerance 
      2 | 2 | 0.0039573 | SiftMaxRelativeTolerance 
      3 | 1 | 0.024838 | SiftMaxRelativeTolerance 
      4 | 2 | 0.05929 | SiftMaxRelativeTolerance 
      5 | 2 | 0.11317 | SiftMaxRelativeTolerance 
      6 | 2 | 0.12599 | SiftMaxRelativeTolerance 
      7 | 2 | 0.13802 | SiftMaxRelativeTolerance 
      8 | 3 | 0.15937 | SiftMaxRelativeTolerance 
      9 | 2 | 0.15923 | SiftMaxRelativeTolerance
分解停止是因為殘留訊號的極值數小於'MaxNumExtrema'。

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

HHT(IMF,FS);

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

視覺化訊號的殘餘和內在模式功能

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

載入非平穩訊號資料,並可視化混合正弦訊號。

load('sinusoidalSignalExampleData.mat'
'X''fs');

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

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

emd(X,'Interpolation''pchip');
目前的IMF | #Sift Iter | 相對Tol | 停止標準命中   
      1 | 2 | 0.026352 | SiftMaxRelativeTolerance 
      2 | 2 | 0.0039573 | SiftMaxRelativeTolerance 
      3 | 1 | 0.024838 | SiftMaxRelativeTolerance 
      4 | 2 | 0.05929 | SiftMaxRelativeTolerance 
      5 | 2 | 0.11317 | SiftMaxRelativeTolerance 
      6 | 2 | 0.12599 | SiftMaxRelativeTolerance 
      7 | 2 | 0.13802 | SiftMaxRelativeTolerance 
      8 | 3 | 0.15937 | SiftMaxRelativeTolerance 
      9 | 2 | 0.15923 | SiftMaxRelativeTolerance
分解停止是因為殘留訊號的極值數小於'MaxNumExtrema'。

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

右鍵單擊圖中的空白區域以開啟IMF選擇器視窗。使用IMF選擇器選擇地檢視生成的IMF,原始訊號和殘差。

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

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

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

大資料部落——中國專業的第三方資料服務提供商,提供定製化的一站式資料探勘和統計分析諮詢服務

統計分析和資料探勘諮詢服務 :y0.cn/teradat(諮詢服務請聯絡官網客服

點選這裡給我發訊息QQ:3025393450

【服務場景】  

科研專案;  公司專案外包 ;線上線下一對一培訓 ;資料採集;學術研究;報告撰寫;市場調查。

【大資料部落】提供定製化的一站式資料探勘和統計分析諮詢服務

【大資料部落】大資料部落提供定製化的一站式資料探勘和統計分析諮詢服務

分享最新的大資料資訊,每天學習一點資料分析,讓我們一起做有態度的資料人【大資料部落】大資料部落提供定製化的一站式資料探勘和統計分析諮詢服務

微信客服號:lico_9e

QQ交流群:186388004 大資料部落

【大資料部落】r語言電商網站爬蟲


歡迎關注微信公眾號,瞭解更多資料乾貨資訊!

【大資料部落】大資料部落提供定製化的一站式資料探勘和統計分析諮詢服務