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经验模态分解残差量分析的详细内容，更多关于PyEMD经验模态分解残差量的资料请关注脚本之家其它相关文章！
相关文章
这篇文章主要介绍了Python返回真假值（True or False）小技巧,本文探讨的是最简洁的条件判断语句写法,本文给出了两种简洁写法,需要的朋友可以参考下2015-04-04
这篇文章主要介绍了Python-OpenCV的基本使用和相关函数介绍，主要包括图像的读取保存图像展示问题，结合实例代码给大家介绍的非常详细,需要的朋友可以参考下2022-05-05
浅析Python 3 字符串中的 STR 和 Bytes 有什么区别
Python2的str相当于Python3的Bytes，而Unicode相当于Python3的Bytes。这篇文章主要介绍了Python 3 字符串中的 STR 和 Bytes 究竟有什么区别？需要的朋友可以参考下2018-10-10
这篇文章主要介绍了Python实现Keras搭建神经网络训练分类模型教程，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧2020-06-06
读写文件是最常见的IO操作,Python内置了读写文件的函数,下面这篇文章主要给大家介绍了关于Python文件读写操作的相关资料,文中通过实例代码介绍的非常详细,需要的朋友可以参考下2023-05-05
Python之lambda匿名函数及map和filter的用法
今天小编就为大家分享一篇关于Python之lambda匿名函数及map和filter的用法，小编觉得内容挺不错的，现在分享给大家，具有很好的参考价值，需要的朋友一起跟随小编来看看吧2019-03-03
本篇文章主要介绍了python+selenium识别验证码并登录的示例代码，小编觉得挺不错的，现在分享给大家，也给大家做个参考。一起跟随小编过来看看吧2017-12-12
这篇文章主要介绍了深入了解Python 中线程和进程区别，一个进程由一个或者多个线程组成，线程是一个进程中代码的不同执行路线。切换进程需要的资源比切换线程的要多的多，下面来了解更多的详细内容，需要的小伙伴可以参考一下2022-03-03
Flask中最方便用的数据库框架是flask_sqlalchamy,是对 SQLAlchamy 在 Flask 中的扩展,它主要在于简化Flask 中 sqlalchamy的使用,本篇文章给大家介绍Python Flask的数据库操作使用方法,感兴趣的朋友一起看看吧2024-02-02
这篇文章主要为大家详细介绍了python实现多图像叠置输出，文中示例代码介绍的非常详细，具有一定的参考价值，感兴趣的小伙伴们可以参考一下2021-04-04
最新评论