请教高手（emd法求信号的逼近分量）

lz1736 · 发表于 2006-5-18 11:16

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

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

x

对于给定的信号,Huang 介绍的EMD 方法是:首先找到信号的极大值和极小值,通过3 次样条拟合,从而获得信号的上包络曲线和下包络曲线,计算上下包络曲线在每一点上的平均值,从而获得一平均值曲线m1 ,设分析信号为x ( t) ,则
x ( t) - m1 = C1
　　从理论上讲, C1 即为第一阶IMF 分量,然后,从原始信号中减去C1 即可获得信号的逼近分量R1 。
x ( t) - C1 = R1
　　对R1 重复上面的过程,就可以获得第二阶IMF 分量。通过EMD 方法对信号的一次次的筛分,就可以获得信号的多个IMF 分量和一个逼近分量Rn ,从而信号可由下式表示:
x ( t) = Σni =1Ci + Rn
请问以上过程matlab有现成的工具包吗？如果没有，找到信号的极大值和极小值,通过3 次样条拟合,从而获得信号的上包络曲线和下包络曲线,计算上下包络曲线在每一点上的平均值,从而获得最后的逼近分量Rn ？
请高人指教。
另附上程序，希望高手帮忙看看为什么运行出错，谢谢！！

suffer · 发表于 2006-5-18 11:35

程序没附上，请重新上传

suffer · 发表于 2006-5-18 11:36

matlab有现成的emd工具包，你可以去信号处理版找

qiuqia17 · 发表于 2006-5-18 12:15

不要在 LOCAL 直接写入信号

你在LOCAL定义了 T=[

]；

在建立新M文件，将T 还成 X 再用 [imf,ort,nbits] = emd_local(x）就可以得到 IMF

是不是这个意思

lz1736 · 发表于 2006-5-18 15:43

不好意思，我太笨了，没看懂楼上这为仁兄是什么意思，麻烦好人做到底，能不能写成m文件传上来？
小弟先谢过了！！

happy · 发表于 2006-5-18 15:50

emd工具包：http://forum.vibunion.com/thread-10456-1-1.html

账号		自动登录	找回密码
密码			我要加入

[HHT] 请教高手（emd法求信号的逼近分量）

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

回复：（lz1736）请教高手（emd法求信号的逼近分量）