epicure 发表于 2007-6-14 10:11

求助!EMD中的IMF

小弟初学HHT ,有小白的话不要笑话我。
我合成了一个信号 从10.8HZ开始,以7HZ递增
k=[1,1,1,1,0,1,...
0,1,0,0,0,0,...
]
s=0;
for i=1:1:12
s=s+k(i)*sin(pi*(10.8+(i-1)*7)*t+pha(i)*pi);
end
我使用EMD分解,但是得出IMF只有3层,而我在一些论文中看到对频率不规则的信号可分为很多层
imf=emd(s);
=hhspectrum(imf);
=toimage(A,f);
disp_hhs(im);
colormap(flipud(gray))

是我合成的信号的频率差太低,还是该信号的各个分量有可能出现倍数关系,使内够模态信号只有3层??
哪位大哥解答一下
谢谢

[ 本帖最后由 epicure 于 2007-6-14 10:13 编辑 ]

zhlong 发表于 2007-6-14 10:22

回复 #1 epicure 的帖子

你换个别的合成信号吧,你这信号分量也不多,k里面那么多都是0。
本版有HHT的例子

zhangnan3509 发表于 2007-6-14 10:26

回复 #2 zhlong 的帖子

每次增加的7HZ与10.8相比大于0.5小于2倍,没有别的的措施,看这个程序主要是还是根据论坛上zhlong求时频图的程序来的,所以这是肯定分不开的。

epicure 发表于 2007-6-14 10:34

那如果改变递增的频率,想分出结果是不是要大于10.8的2倍呢

zhangnan3509 发表于 2007-6-14 10:36

回复 #4 epicure 的帖子

等于也是可以的,咱们版面以前有人做过例子

zhlong 发表于 2007-6-14 10:49

回复 #4 epicure 的帖子

其实不仅和频率之比相关,还与幅值比、相位差等有关系。

zhangnan3509 发表于 2007-6-14 10:53

回复 #6 zhlong 的帖子

总之不是紧靠这几行程序,就能把信号分解得到好的结果

epicure 发表于 2007-6-15 11:08

谢谢   
:lol :lol :lol :lol

继续学习中

brightfay 发表于 2007-6-15 15:00

刚才看见大家讨论,说分解结果受到频率比,幅值比,相位比的影响,但是实际信号并不是我们能控制的啊。我现在刚接触EMD,说实话做的比较郁闷:'(

zhlong 发表于 2007-6-15 16:19

回复 #9 brightfay 的帖子

所以EMD并不是适合处理任何信号的!
页: [1]
查看完整版本: 求助!EMD中的IMF