[其他] 希爾伯特-黃的訊號拆解後的還原

看板MATLAB作者 (嗶波)時間12年前 (2013/11/20 11:39), 編輯推噓0(009)
留言9則, 3人參與, 最新討論串1/1
各位先進日安,小弟是matlab的新手有個問題想請教各位 在eemd進行HHT拆解後,得到分離出來的9個IMFseries 在確認出所謂noise後要如何針對剩下的幾個IMF進行還原(reconstruction)? 看著matlab的書弄了好久就是沒有成功 是否能請板上的高手指點我一下..萬分感謝 之前弄出來的是 function [x]=reconstruction(c) % xt is orignal signal % c is IMF [m,n]=size(c) % c is 9*307, so xt will be 9*307; m is number of IMF x=zeros(1,n) for i=1:n for j=1:m x(i)=x(i)+c(j,i) end end 但是就是動不了,┌(_#_)┐懇求指導 -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 123.192.235.106 ※ 編輯: rosso0922 來自: 123.192.235.106 (11/20 11:47)

11/20 13:14, , 1F
你denoise前能還原嗎? 如果可以,那就是你後來送給recons
11/20 13:14, 1F

11/20 13:14, , 2F
的input結構有問題.
11/20 13:14, 2F

11/20 13:16, , 3F
如果你的emd是從黃鍔院士的網頁上抓的,程式碼是ok的.
11/20 13:16, 3F

11/20 14:43, , 4F
這個是我從中央大學的研究中心網頁抓來的
11/20 14:43, 4F

11/20 14:44, , 5F
eemd應該是沒有問題,但事實上denoise前我也還原不回去
11/20 14:44, 5F

11/20 14:45, , 6F
我第一次寫還原的式子..,這個跟fastica的還原法似乎有
11/20 14:45, 6F

11/20 14:45, , 7F
所不同,另外謝謝Y大也請您繼續指點
11/20 14:45, 7F

11/20 15:03, , 8F
不是直接把IMF的向量相加嗎?雙重迴圈的作用也是如此吧?
11/20 15:03, 8F

11/20 17:31, , 9F
阿阿! 我馬上試試
11/20 17:31, 9F
文章代碼(AID): #1IZ2xaxU (MATLAB)
文章代碼(AID): #1IZ2xaxU (MATLAB)