码迷,mamicode.com
首页 > 其他好文 > 详细

信号分解 - PyEMD

时间:2020-04-17 09:51:07      阅读:304      评论:0      收藏:0      [点我收藏+]

标签:5*   param   信号   using   evo   fine   sim   python实现   content   

安装:pip install EMD-signal

文档:https://pyemd.readthedocs.io/en/latest/

 

我大概看了原理,看的不是特别明白,只记录一下我觉得有用的地方

 

EMD,经验模态分解,是一种信号分解的技术;

它提出了一个概念叫 基本模态分量 IMF,

EMD 用于处理非平稳信号,可用于任意数据,基于数据本身进行分解;

EMD 把一个信号分解成 多个 IMF,每个 IMF 具有线性和非线性的特点,还有一个 信号残余分量,常常代表信号的直流分量或者信号的趋势;

EMD 分解得到的 IMF,频率逐渐降低,尺度各不相同;

EMD 分解得到的前几个模态分量,通常集中了原始信号中最显著、最重要的信息; 

 

EMD 分解容易造成 模态混合,表现为下列现象之一:

1. 在同一个 IMF 中,尺度分布范围很宽却又各不相同的信号;

技术图片

2. 在不同 IMF 中,存在尺度相近的信号;

模态混合使得 IMF 失去单一特征尺度,形成尺度混杂的震荡,失去原有的物理意义

 

EEMD,集合经验模态分解,解决了 EMD 模态混合的现象

 

Python EMD 用法

##### 基本用法
import numpy as np
from PyEMD import EMD
import pylab as plt

s = np.random.random(100)
emd = EMD()
IMFs = emd.emd(s)


##### 示例
# Define signal
t = np.linspace(0, 1, 200)
s = np.cos(11*2*np.pi*t*t) + 6*t*t

# Execute EMD on signal
IMF = EMD().emd(s,t)
N = IMF.shape[0]+1

# Plot results
plt.subplot(N,1,1)
plt.plot(t, s, r)
plt.title("Input signal: $S(t)=cos(22\pi t^2) + 6t^2$")
plt.xlabel("Time [s]")

for n, imf in enumerate(IMF):
    plt.subplot(N,1,n+2)
    plt.plot(t, imf, g)
    plt.title("IMF "+str(n+1))
    plt.xlabel("Time [s]")

plt.tight_layout()
plt.savefig(simple_example)
plt.show()

输出

技术图片

 

Python EEMD 用法

from PyEMD import EEMD
import numpy as np
import pylab as plt


def main():
    # Define signal
    t = np.linspace(0, 1, 200)

    sin = lambda x,p: np.sin(2*np.pi*x*t+p)
    S = 3*sin(18,0.2)*(t-0.2)**2
    S += 5*sin(11,2.7)
    S += 3*sin(14,1.6)
    S += 1*np.sin(4*2*np.pi*(t-0.8)**2)
    S += t**2.1 -t

    # Assign EEMD to `eemd` variable
    eemd = EEMD()

    # Say we want detect extrema using parabolic method
    # emd = eemd.EMD
    # emd.extrema_detection="parabol"
    # eemd.trials = 50
    # eemd.noise_seed(12345)

    # Execute EEMD on S
    eIMFs = eemd.eemd(S, t, -1)
    nIMFs = eIMFs.shape[0]

    # Plot results
    plt.figure(figsize=(12,9))
    plt.subplot(nIMFs+1, 1, 1)
    plt.plot(t, S, r)

    for n in range(nIMFs):
        plt.subplot(nIMFs+1, 1, n+2)
        plt.plot(t, eIMFs[n], g)
        plt.ylabel("eIMF %i" %(n+1))
        plt.locator_params(axis=y, nbins=5)

    plt.xlabel("Time [s]")
    plt.tight_layout()
    plt.savefig(eemd_example, dpi=120)
    plt.show()


if __name__ == __main__:
    main()

 

 

 

 

参考资料:

https://www.cnblogs.com/RoseVorchid/p/12030980.html  EEMD算法python实现,原理讲得不咋地

http://www.360doc.com/content/19/0806/12/41357686_853284987.shtml#  EEMD 原理

https://www.cnblogs.com/mikawong/p/7682467.html  记录了 PyEMD 的一些资源链接

 

《Adaboost_SVM集成模型的滚动轴承早期故障诊断》  电子书

信号分解 - PyEMD

标签:5*   param   信号   using   evo   fine   sim   python实现   content   

原文地址:https://www.cnblogs.com/yanshw/p/12707121.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!