作者rosso0922 (嗶波)
看板MATLAB
標題[其他] 希爾伯特-黃的訊號拆解後的還原
時間Wed Nov 20 11:39:14 2013
各位先進日安，小弟是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)
1F：→ YoursEver:你denoise前能還原嗎? 如果可以,那就是你後來送給recons 11/20 13:14
2F：→ YoursEver:的input結構有問題. 11/20 13:14
3F：→ YoursEver:如果你的emd是從黃鍔院士的網頁上抓的,程式碼是ok的. 11/20 13:16
4F：→ rosso0922:這個是我從中央大學的研究中心網頁抓來的 11/20 14:43
5F：→ rosso0922:eemd應該是沒有問題，但事實上denoise前我也還原不回去 11/20 14:44
6F：→ rosso0922:我第一次寫還原的式子..，這個跟fastica的還原法似乎有 11/20 14:45
7F：→ rosso0922:所不同，另外謝謝Y大也請您繼續指點 11/20 14:45
8F：→ aia:不是直接把IMF的向量相加嗎？雙重迴圈的作用也是如此吧？ 11/20 15:03
9F：→ rosso0922:阿阿! 我馬上試試 11/20 17:31