详解Python中PyEMD残差量的两种实现形式是怎样的-群英

上一篇：python selenium模拟点击的过程及实现是怎样下一篇：python中logging模块用来做什么，简单的用法怎样

在实际应用中，我们有时候会遇到“详解Python中PyEMD残差量的两种实现形式是怎样的”这样的问题，我们该怎样来处理呢？下文给大家介绍了解决方法，希望这篇“详解Python中PyEMD残差量的两种实现形式是怎样的”文章能帮助大家解决问题。

PyEMD是经验模态分解 (EMD)及其变体的Python实现，EMD最流行的扩展之一是集成经验模态分解 (EEMD)，它利用了噪声辅助执行的集成。

顾名思义，这个包中的方法获取数据（信号）并将其分解为一组组件。所有这些方法理论上都应该将信号分解为同一组分量，但实际上有很多细微差别和不同的方法来处理噪声。无论采用何种方法，获得的分量通常称为本征模态函数(IMF)，以强调它们包含固有（自身）属性，即特定振荡（模态）。（以上来自官方文档）

两种实现形式

最近尝试实现CEEMDAN，CEEMADN也是EMD的一种变体。按照官方API，有以下两种形式的写法：

大部分博客采用的是第一种写法：

ceemdan = CEEMDAN()
ceemdan.ceemdan(load)
imfs, res = ceemdan.get_imfs_and_residue()
vis = Visualisation()
vis.plot_imfs(imfs, res)

这种写法得到的图为：

第二种写法，也是官方的写法：

ceemdan = CEEMDAN()(load)
imfs, res = ceemdan[:-1], ceemdan[-1]
vis = Visualisation()
vis.plot_imfs(imfs, res)

得到的图示为：

可以发现，两张图最大的区别在于Res。一般论文中给出的图示是第二种。

我们尝试输出第一种方法中的Res：

[ 0.00000000e+00 -2.84217094e-14 0.00000000e+00 ... -2.84217094e-14
0.00000000e+00 0.00000000e+00]

可以发现其量级特别小，还原时我们可以不加上这一部分。而方法二中的Res显然量级是最大的，还原时必须加上。

因此，方法一中的Res是我们一般人所理解的残余量，在真正进行建模时可以不考虑。

方法一中画图时如果我们不包含残余量，即：

vis.plot_imfs(imfs=imfs, residue=res, include_residue=False)

我们将得到：

在github上经过交流后，得到如下结论：

方法一中的Res是真正意义上的残余量，或许叫残差更合适一点，也就是分解之后不能再分解的部分。在PyEMD的源码中被定义为：

S * scale_s - np.sum(self.C_IMF, axis=0)

因此，ceemdan.get_imfs_and_residue()实际上得到的是最终的IMF和重建误差，而不是残差。

残差的正确获取方式是ceemdan[-1]。

“详解Python中PyEMD残差量的两种实现形式是怎样的”的内容就介绍到这里了，感谢大家的阅读。如果想了解更多行业技术相关的知识可以关注群英网络网站，小编每天都会为大家更新不同的知识。

标签： python

免责声明：本站发布的内容（图片、视频和文字）以原创、转载和分享为主，文章观点不代表本网站立场，如果涉及侵权请联系站长邮箱：mmqy2019@163.com进行举报，并提供相关证据，查实之后，将立刻删除涉嫌侵权内容。

上一篇：python selenium模拟点击的过程及实现是怎样下一篇：python中logging模块用来做什么，简单的用法怎样