python实现PyEMD经验模态分解残差量分析

 更新时间:2022年05月10日 10:55:32   作者:Cyril_KI  
这篇文章主要为大家介绍了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经验模态分解残差量的资料请关注脚本之家其它相关文章!

您可能感兴趣的文章:

相关文章

  • 从零开始的TensorFlow+VScode开发环境搭建的步骤(图文)

    从零开始的TensorFlow+VScode开发环境搭建的步骤(图文)

    这篇文章主要介绍了从零开始的TensorFlow+VScode开发环境搭建的步骤(图文),文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧
    2020-08-08
  • Python超有用的多版本管理工具pyenv

    Python超有用的多版本管理工具pyenv

    使用不同的Python版本,如果我们把需要的不同版本的Python都下载到服务器上,管理起来会非常困难,多版本并存又容易互相干扰,接下来就来介绍一个Python环境管理工具:pyenv,就可轻松的在多个版本的Python之间自由切换,需要的朋友可以参考下
    2021-09-09
  • Python使用Opencv实现边缘检测以及轮廓检测的实现

    Python使用Opencv实现边缘检测以及轮廓检测的实现

    这篇文章主要介绍了Python使用Opencv实现边缘检测以及轮廓检测的实现,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧
    2020-12-12
  • python if三元表达式实例用法详解

    python if三元表达式实例用法详解

    在本篇文章里小编给大家整理的是一篇关于python if三元表达式实例用法详解内容,有兴趣的朋友们可以学习下。
    2021-04-04
  • python实现单线程多任务非阻塞TCP服务端

    python实现单线程多任务非阻塞TCP服务端

    这篇文章主要为大家详细介绍了python实现单线程多任务非阻塞TCP服务端的相关资料,具有一定的参考价值,感兴趣的小伙伴们可以参考一下
    2017-06-06
  • Python之os操作方法(详解)

    Python之os操作方法(详解)

    下面小编就为大家带来一篇Python之os操作方法(详解)。小编觉得挺不错的,现在就分享给大家,也给大家做个参考。一起跟随小编过来看看吧
    2017-06-06
  • Python中的函数参数(位置参数、默认参数、可变参数)

    Python中的函数参数(位置参数、默认参数、可变参数)

    这篇文章主要介绍了Python中的函数参数:位置参数、默认参数、可变参数、关键字参数和命名关键字参数,需要的小伙伴可以参考下面文章内容
    2021-09-09
  • 亲手教你用Python打造一款摸鱼倒计时界面

    亲手教你用Python打造一款摸鱼倒计时界面

    前段时间在微博看到一段摸鱼人的倒计时模板,感觉很有意思,于是我用了一个小时的时间写一个页面出来,下面小编把实现过程分享给大家,对Python摸鱼倒计时界面感兴趣的朋友一起看看吧
    2021-12-12
  • OpenCV哈里斯(Harris)角点检测的实现

    OpenCV哈里斯(Harris)角点检测的实现

    这篇文章主要介绍了OpenCV哈里斯 (Harris)角点检测,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧
    2020-01-01
  • python spotlight库简化交互式方法探索数据分析

    python spotlight库简化交互式方法探索数据分析

    这篇文章主要为大家介绍了python spotlight库简化的交互式方法探索数据,有需要的朋友可以借鉴参考下,希望能够有所帮助,祝大家多多进步,早日升职加薪
    2024-01-01

最新评论