求助!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 编辑 ]
回复 #1 epicure 的帖子
你换个别的合成信号吧,你这信号分量也不多,k里面那么多都是0。本版有HHT的例子
回复 #2 zhlong 的帖子
每次增加的7HZ与10.8相比大于0.5小于2倍,没有别的的措施,看这个程序主要是还是根据论坛上zhlong求时频图的程序来的,所以这是肯定分不开的。 那如果改变递增的频率,想分出结果是不是要大于10.8的2倍呢回复 #4 epicure 的帖子
等于也是可以的,咱们版面以前有人做过例子回复 #4 epicure 的帖子
其实不仅和频率之比相关,还与幅值比、相位差等有关系。回复 #6 zhlong 的帖子
总之不是紧靠这几行程序,就能把信号分解得到好的结果 谢谢:lol :lol :lol :lol
继续学习中 刚才看见大家讨论,说分解结果受到频率比,幅值比,相位比的影响,但是实际信号并不是我们能控制的啊。我现在刚接触EMD,说实话做的比较郁闷:'(
回复 #9 brightfay 的帖子
所以EMD并不是适合处理任何信号的!
页:
[1]