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

[rosso0922 (嗶波)]

各位先進日安，小弟是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)
→ YoursEver:你denoise前能還原嗎? 如果可以,那就是你後來送給recons
→ YoursEver:的input結構有問題.
→ YoursEver:如果你的emd是從黃鍔院士的網頁上抓的,程式碼是ok的.
→ rosso0922:這個是我從中央大學的研究中心網頁抓來的
→ rosso0922:eemd應該是沒有問題，但事實上denoise前我也還原不回去
→ rosso0922:我第一次寫還原的式子..，這個跟fastica的還原法似乎有
→ rosso0922:所不同，另外謝謝Y大也請您繼續指點
→ aia:不是直接把IMF的向量相加嗎？雙重迴圈的作用也是如此吧？
→ rosso0922:阿阿! 我馬上試試