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

lz1736

　　对于给定的信号,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 ?

　　请高人指教。

　　另附上程序，请高手帮忙看看为什么运行出错。谢谢!!