有谁给我看看我这个信号分解成IMF第一个分量对不对啊？
上面一个是IMF第一个分量，下面是差残，总觉得不对啊，论坛中的大牛帮忙指正一下啊
不对，是不是你程序中残差把IMF给覆盖了造成的？ 本帖最后由 Crocodilesky 于 2011-3-7 21:38 编辑
EMD分解是从高到低。从你的残差看，似乎分解不对。 你的输入信号是什么 ？只有知道了输入信号的成分才能知道你的IMF对不对 输入是这样的，大家给看一下吧，是怎么回事啊，我对HHT刚研究不久！
回复 5 # flq103 的帖子
你这是什么信号啊？ 你这个是不是分段余弦信号啊，要是的话分解的没错，直接满足IMF条件，不用EMD了~，直接做Hilbert变化就行~ 什么信号？送过来瞧一瞧。 我觉得你第一个图，是将5个IMF分量画在一起了，这个残差可能是信号和那5个IMF想减之后得到的。 p=zeros(1,20001);
t=0:0.0001:2;
aa=find(t>0.325&t<=0.556);
p(aa)=0.9;
aa=find (t>1.358&t<=1.647);
p(aa)=1.1;
aa=find ((t>0.556&t<=1.358)|(t>0&t<=0.325)|(t>1.647&t<=2));
p(aa)=1;
s=p.*cos(2*pi*50*t).*(t>0&t<=2) +0.1*sin(2*pi*1.5*50*t).*(t>0&t<=0.5)+0.25*sin(2*pi*3*50*t).*(t>=0.5&t<=1)+0.2*sin(2*pi*4*50*t).*(t>=1&t<=2);
plot(t,s)上面是该程序，高手们看看吧！！{:{13}:}
本帖最后由 gczhang 于 2011-3-14 10:55 编辑
我试了一下，结果供参考.
本帖最后由 gczhang 于 2011-3-14 10:59 编辑
总残差图，非常小。 本帖最后由 flq103 于 2011-3-17 10:38 编辑
回复 12 # gczhang 的帖子
那我分析出来残差怎么是那样的啊，高手帮忙分析一下啊，还有你分析出来的瞬时幅值和瞬时频率怎么样啊 本帖最后由 gczhang 于 2011-3-17 13:20 编辑
我感觉你的那个不是残差，更像是最后一阶 IMF8。看我得到的IMF8，和你的所谓残差基本一致。
不知怎么的，图片传不上来。
回复 13 # flq103 的帖子
是这样啊， 我用的是法国那个人编程序，你的HHT程序边界效应是不是处理过？
页:
[1]
2