zhlong
发表于 2007-5-8 21:54
eight版主,还有一个问题请教,图(1)是我从一篇学位论文上看到的信号及其频谱,说明信号有基频和二倍频的存在;图(2)是该文中EMD分解结果,说明信号不能分解,已满足IMF的定义。图(3)是信号的HHT时频图。
按道理,这个信号确实是满足了IMF的定义不能分解。但Fourier谱又说明信号有基频和其二倍频的存在,这本是EMD的一个缺陷吧(二倍频的幅值比基频幅值小很多的时候,两个信号相加的波形和基频波形区别不大,EMD把它当作是一个单频率分量了)。但是该文作者却根据不正确的分解得到的不正确的HHT时频图说该信号存在调相(调频)现象。不知道我的理解是否正确,请您指教!
[ 本帖最后由 zhlong 于 2007-5-8 21:55 编辑 ]
eight
发表于 2007-5-8 22:05
前面的理解都正确,但是对于以下句子,由于我没有看过作者的文章,所以不知道是他的原话,还是你的评价。
但是该文作者却根据不正确的分解得到的不正确的HHT时频图说该信号存在调相(调频)现象
何谓“不正确的分解”?你的意思是他运用EMD能够分解出几个而非IMF?调频(调相)现象个人认为应该与分解得到多少个IMF无关吧,因为就算只得到一个IMF(正如你说的例子),也不能排除原信号存在调频(调相)现象的可能。我对信号处理也不是十分了解,仅供参考,需要静候两位版主(yangzj,zhangnan3509)及其他高手(如songzy41, pengzk 等等)的高见了,呵呵
zhlong
发表于 2007-5-8 22:10
我指的错误的分解是指 根据频谱该信信号由基频和其二倍频两个分量组成, 但EMD却没能将该信号分解为基频和二倍频两个IMF分量
[ 本帖最后由 zhlong 于 2007-5-8 22:12 编辑 ]
eight
发表于 2007-5-8 22:12
原帖由 zhlong 于 2007-5-8 22:10 发表 http://forum.vibunion.com/forum/images/common/back.gif
我指的错误的分解是指 根据频谱该信信号由基频和其二倍频两个分量组成, 但EMD却没能将该信号分解维基频和二倍频两个分量
我认为EMD对此应该无能为力,除非结合其他方法(例如小波,或者mask等)
zhlong
发表于 2007-5-8 22:14
EMD 一定不能分解基频和其倍频分量的叠加吗?
zhlong
发表于 2007-5-8 22:15
原帖由 eight 于 2007-5-8 22:12 发表 http://forum.vibunion.com/forum/images/common/back.gif
我认为EMD对此应该无能为力,除非结合其他方法(例如小波,或者mask等)
那这种分解得到的HHT时频图可信吗?我要是有原始信号,我想用STFT分析得到的时频图应该不是这样的。
但作者也没有将这个结果和别的方法相比
[ 本帖最后由 zhlong 于 2007-5-8 22:17 编辑 ]
zhlong
发表于 2007-5-8 22:20
下面图中的信号也是由 基频和二倍频组成,EMD结果显示能够分解
eight
发表于 2007-5-8 22:25
原帖由 zhlong 于 2007-5-8 22:14 发表 http://forum.vibunion.com/forum/images/common/back.gif
EMD 一定不能分解基频和其倍频分量的叠加吗?
我的意思是指频率很近的两个信号。你第一个例子原信号明显就是一个IMF,第二个例子原信号明显就有骑波
zhlong
发表于 2007-5-8 22:29
两个例子都是基频和二倍频合成,即 高频与低频的频率比为2
那两个信号中频率接近程度应该是一样的啊,根据信号的频谱图,造成它们之间区别的原因是不是基频和二倍频幅值比例的问题呢
[ 本帖最后由 zhlong 于 2007-5-8 22:40 编辑 ]
rc-hw-0002
发表于 2007-5-8 22:54
原帖由 zhlong 于 2007-5-8 22:29 发表 http://forum.vibunion.com/forum/images/common/back.gif
两个例子都是基频和二倍频合成,即 高频与低频的频率比为2
那两个信号中频率接近程度应该是一样的啊,根据信号的频谱图,造成它们之间区别的原因是不是基频和二倍频幅值比例的问题呢
可能是幅值比原因吧,在一篇论文上看到A1*w1>A2*w2(w1>w2)是两分量能正确分解的充分条件
[ 本帖最后由 rc-hw-0002 于 2007-5-8 23:08 编辑 ]
zhlong
发表于 2007-5-8 22:56
原帖由 eight 于 2007-5-8 22:25 发表 http://forum.vibunion.com/forum/images/common/back.gif
我的意思是指频率很近的两个信号。你第一个例子原信号明显就是一个IMF,第二个例子原信号明显就有骑波
还请教eight版主一个问题,IMF的两条定义是不是就是要使得基于Hilbert变换的“瞬时频率”有意义,即将信号分解为多个单频率分量之和。
显然我举的第一个例子中信号包含两个频率分量,所以直接对这样的信号进行hilbert变换得到的HHT时频图反应的信号频率与时间之间的关
系是不正确的。
zhlong
发表于 2007-5-8 22:58
原帖由 rc-hw-0002 于 2007-5-8 22:54 发表 http://forum.vibunion.com/forum/images/common/back.gif
可能是幅值比原因吧,在一篇论文上看到A1*w1>A2*w2是两分量能正确分解的充分条件
非常感谢提供线索,能告知哪篇文章吗?呵呵,我这方面文献看的太少了
eight
发表于 2007-5-8 23:00
原帖由 zhlong 于 2007-5-8 22:29 发表 http://forum.vibunion.com/forum/images/common/back.gif
两个例子都是基频和二倍频合成,即 高频与低频的频率比为2
那两个信号中频率接近程度应该是一样的啊,根据信号的频谱图,造成它们之间区别的原因是不是基频和二倍频幅值比例的问题呢
的确有可能是幅值的问题,这也是不容忽视的
rc-hw-0002
发表于 2007-5-8 23:02
原帖由 zhlong 于 2007-5-8 22:58 发表 http://forum.vibunion.com/forum/images/common/back.gif
非常感谢提供线索,能告知哪篇文章吗?呵呵,我这方面文献看的太少了
多频信号经验模态分解的理论研究及应用
沈国际,陶利民,陈仲生
zhlong
发表于 2007-5-8 23:03
原帖由 rc-hw-0002 于 2007-5-8 22:54 发表 http://forum.vibunion.com/forum/images/common/back.gif
可能是幅值比原因吧,在一篇论文上看到A1*w1>A2*w2是两分量能正确分解的充分条件
t=1:500;
t=t*1/2000;
x=sin(2*pi*20*t);
y=0.2*sin(2*pi*40*t);
plot(t,x,t,x+y)
A1*w1=20>A2*w2=8
但是从图中看,显然这两个信号也不能分解!