EMD 经验模态分解

1. EMD 的起因

在信号处理方面我们了解了时域处理方法(如有效值、峭度)、频域处理方法(如频谱、功率谱)以及一些时频域处理方法(如小波分析).

EMD(Empirical Mode Decomposition)作为时频域的处理方法,相对于同样是时频域方法的小波分析有什么好处呢?

EMD最显著的特点,就是其克服了基函数无自适应性的问题。啥意思呢?

回忆小波分析部分的内容,我们会知道小波分析是需要选定某一个小波基的,小波基的选择对整个小波分析的结果影响很大,一旦确定了小波基,在整个分析过程中将无法更换,即使该小波基在全局可能是最佳的,但在某些局部可能并不是,所以小波分析的基函数缺乏适应性。

通俗的说,用EMD有什么好处呢?对于一段未知信号,不需要做预先分析与研究,就可以直接开始分解。这个方法会自动按照一些固模式按层次分好,而不需要人为设置和干预。

再通俗一点,EMD就像一台机器,把一堆混在一起的硬币扔进去,他会自动按照1元、5毛、1毛、5分、1分地分成几份。

2. 各个本征函数的分量

涵模态分量(Intrinsic Mode Functions, IMF)就是原始信号被EMD分解之后得到的各层信号分量。EMD的提出人黄锷认为,任何信号都可以拆分成若干个内涵模态分量之和。而内涵模态分量有两个约束条件:

1)在整个数据段内,极值点的个数和过零点的个数必须相等或相差最多不能超过一个。

2)在任意时刻,由局部极大值点形成的上包络线和由局部极小值点形成的下包络线的平均值为零,即上、下包络线相对于时间轴局部对称。

啥意思?

用不严谨的语言和灵魂画师来解释一下:

1.图线要反复跨越x轴,像这样:

在整个数据段内,极值点的个数和过零点的个数必须相等或相差最多不能超过一个

在这里插入图片描述
而不能像这样某次穿过零点后出现多个极点:

在这里插入图片描述

  1. 包络要对称,像这样:
    包络线对称

在这里插入图片描述
而不能像这样:

在这里插入图片描述

看个正常点的例子吧:

在这里插入图片描述

上图由7张图片组成,其中第1张为原始信号,后边依次为EMD分解之后得到的6个分量,分别叫做IMF1~IMF5,最后一张图为残差,每一个IMF分量代表了原始信号中存在的一种内涵模态分量。可以看出,每个IMF分量都是满足这两个约束条件的。

3. EMD 分解步骤

EMD的分解过程是简单直观的:

1)根据原始信号上下极值点,分别画出上、下包络线。
上、下包络线
在这里插入图片描述

2)求上、下包络线的均值,画出均值包络线。均值包络线。
在这里插入图片描述
3)原始信号减均值包络线,得到中间信号。

原始信号减均值包络线
在这里插入图片描述
4)判断该中间信号是否满足IMF的两个条件,如果满足,该信号就是一个IMF分量;如果不是,以该信号为基础,重新做1)~4)的分析。IMF分量的获取通常需要若干次的迭代。

不满足约束2,需要继续迭代
在这里插入图片描述

使用上述方法得到第一个IMF后,用原始信号减IMF1,作为新的原始信号,再通过1)~4)的分析,可以得到IMF2,以此类推,完成EMD分解。

在这里插入图片描述

上述例子中的图来自http://perso.ens-lyon.fr/patrick.flandrin/emd.ppt

这个ppt中用了几十页的图片完整地展示了EMD分解的过程,看一下该文档相信你会对该过程有一个更清楚的了解。

原文出自这里

https://zhuanlan.zhihu.com/p/40005057
https://zhuanlan.zhihu.com/p/121803211

  • 1
    点赞
  • 5
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
EMD经验模态分解算法是一种局部平稳性分解技术,用于将时间序列分解成若干个组成部分,每个部分中包含不同的时间尺度和不同频率的信息,以更好地描述时间序列的特征。该算法能够充分利用时间序列数据的内在特征结构,适用性广泛。在EMD分解实验中,可以通过导入PyEMD工具库,构建信号并提取IMFs(内禀模态函数)和剩余信号。通过绘制IMFs和剩余信号的瞬时频率,可以更好地理解时间序列的特征。\[1\]\[3\]该模型还可以结合机器学习相关算法如神经网络或支持向量回归(SVR)等进行时间序列预测,以提高预测的准确性。\[1\]\[2\] #### 引用[.reference_title] - *1* *2* [基于EMD经验模态分解算法Python程序](https://blog.csdn.net/qq_41728700/article/details/129288324)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^control_2,239^v3^insert_chatgpt"}} ] [.reference_item] - *3* [Python-经验模态分解EMD分解及测试用例](https://blog.csdn.net/weixin_41990278/article/details/123127868)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^control_2,239^v3^insert_chatgpt"}} ] [.reference_item] [ .reference_list ]

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值