作者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