有谁给我看看我这个信号分解成IMF第一个分量对不对啊？

flq103 · 发表于 2011-3-7 12:18

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

上面一个是IMF第一个分量，下面是差残，总觉得不对啊，论坛中的大牛帮忙指正一下啊

nameshenzaiyang · 发表于 2011-3-7 14:56

不对，是不是你程序中残差把IMF给覆盖了造成的？

Crocodilesky · 发表于 2011-3-7 21:35

本帖最后由 Crocodilesky 于 2011-3-7 21:38 编辑

EMD分解是从高到低。从你的残差看，似乎分解不对。

hitdely · 发表于 2011-3-7 21:38

你的输入信号是什么？只有知道了输入信号的成分才能知道你的IMF对不对

flq103 · 发表于 2011-3-8 08:48

输入是这样的，大家给看一下吧，是怎么回事啊，我对HHT刚研究不久！

mni987654321 · 发表于 2011-3-10 14:44

回复 5 # flq103 的帖子

你这是什么信号啊？

Cena马 · 发表于 2011-3-10 15:23

你这个是不是分段余弦信号啊，要是的话分解的没错，直接满足IMF条件，不用EMD了~，直接做Hilbert变化就行~

dsp2008 · 发表于 2011-3-11 08:01

什么信号？送过来瞧一瞧。

gczhang · 发表于 2011-3-11 15:27

我觉得你第一个图，是将5个IMF分量画在一起了，这个残差可能是信号和那5个IMF想减之后得到的。

flq103 · 发表于 2011-3-12 15:36

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)上面是该程序，高手们看看吧！！

gczhang · 发表于 2011-3-14 10:55

本帖最后由 gczhang 于 2011-3-14 10:55 编辑

我试了一下，结果供参考.

gczhang · 发表于 2011-3-14 10:59

本帖最后由 gczhang 于 2011-3-14 10:59 编辑

总残差图，非常小。

flq103 · 发表于 2011-3-17 10:36

本帖最后由 flq103 于 2011-3-17 10:38 编辑

回复 12 # gczhang 的帖子

那我分析出来残差怎么是那样的啊，高手帮忙分析一下啊，还有你分析出来的瞬时幅值和瞬时频率怎么样啊

gczhang · 发表于 2011-3-17 11:35

本帖最后由 gczhang 于 2011-3-17 13:20 编辑

我感觉你的那个不是残差，更像是最后一阶 IMF8。看我得到的IMF8，和你的所谓残差基本一致。
不知怎么的，图片传不上来。

flq103 · 发表于 2011-3-17 21:36

回复 13 # flq103 的帖子

是这样啊，我用的是法国那个人编程序，你的HHT程序边界效应是不是处理过？

账号		自动登录	找回密码
密码			我要加入

[FFT] 有谁给我看看我这个信号分解成IMF第一个分量对不对啊？

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

点评