声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1733|回复: 4

[HHT] 求助:EMD的滤波停止条件

[复制链接]
发表于 2007-3-15 17:06 | 显示全部楼层 |阅读模式

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

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

x
我自己编的EMD程序。其中的零均值条件是自己设定的:上下包络线的均值小于某一个特定的数值。不知道这样对分解效果会有什么影响?
对分解后的IMF求瞬时频率时得到的结果,我觉得不对,频率并没有减小。

[ 本帖最后由 zhlong 于 2007-6-4 17:24 编辑 ]
回复
分享到:

使用道具 举报

发表于 2007-3-15 17:28 | 显示全部楼层
我觉得中止条件对HHT至关重要,而楼主在自己编写的时候,有些随意了,建议好好看看Rilling G的程序以及他们所提出的中止条件在程序中是怎么应用的。如果还是有疑问 可以看看 On empirical mode decomposition and its algorithms 这篇文章
 楼主| 发表于 2007-3-16 15:47 | 显示全部楼层
[imf,nbits]=emd2(x);
for k=1:nbits
c(k,:)=imf(k,:)+i*hilbert(imf(k,:));
[fnormhat(:,k)]=instfreq(c(k,:)');
end
这段求瞬时频率的程序对吗?可为什么求的结果不好呢?
发表于 2007-3-21 11:12 | 显示全部楼层
原帖由 xxf 于 2007-3-16 15:47 发表
=emd2(x);
for k=1:nbits
c(k,:)=imf(k,:)+i*hilbert(imf(k,:));
=instfreq(c(k,:)');
end
这段求瞬时频率的程序对吗?可为什么求的结果不好呢?



请看清楚 hilbert 函数的输出,c=hilbert(a) 的结果是: c = a + b i ,并非 c = b
 楼主| 发表于 2007-3-21 11:45 | 显示全部楼层
我明白了,是我理解错了
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-28 16:37 , Processed in 0.054126 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表