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