信号处理经验模态分解 Empirical Mode Decomposition
经验模态分解(Empirical Mode Decomposition,简称EMD)
经验模态分解EMD是基于瞬时频率、本征模态函数(Intrinsic Mode Function,IMF)的概念,能够将复杂信号分解为若干个IMF分量,每个IMF表征信号的局部特征。依据的是数据自身的时间尺度特征来进行信号分解,无需预先设定任何基函数,因此具有自适应性。再通俗一点,EMD就像一台机器,把一堆混在一起的硬币扔进去,他会自动按照1元、5毛、1毛、5分、1分地分成几份。
经验模态分解的基本思想:将一个频率不规则的波化为多个单一频率的波 残波的形式。原波形 = ∑ IMFs 余波。
EMD方法是基于如下假设基础上的:
(1)信号至少有两个极值点,一个极大值和一个极小值;
(2)特征时间尺度通过两个极值点之间的时间定义;
(3)若数据缺乏极值点但有形变点,则可通过数据微分一次或几次获得极值点,然后再通过积分来获得分解结果。
算法流程如下所示:
EMD分解python代码:
import pandas as pd
import matplotlib.pyplot as plt
%matplotlib inline
import numpy as np
from PyEMD import EMD
import scipy.io.wavfile as wf
plt.rcParams['font.sans-serif']=['SimHei'] #用来正常显示中文标签
plt.rcParams['axes.unicode_minus']=False #用来正常显示符号
Fs = 44100
T = 1/Fs
sample_rate, noised_sigs = wf.read('noised.wav')
t = np.linspace(0,len(noised_sigs)-1,len(noised_sigs))
s = np.array(noised_sigs)
emd = EMD()
imfs = emd.emd(s,t,10)
# Plot results
plt.figure(figsize=(12,15))
plt.subplot(7,2,1)
plt.plot(t, s, 'r')
plt.title("Input signal")
plt.figure(figsize=(12,15))
plt.subplot(7,2,1)
plt.plot(t, s, 'r')
plt.title("IMFS")
for n, imf in enumerate(imfs):
plt.subplot(13,2,n 3)
plt.plot(t[:2000], imf[:2000], 'g')
plt.ylabel("IMF %i" %(n 1))
plt.tight_layout()
plt.show()
同理,集合经验模态分解(EEMD)分解展示如下:
import pandas as pd
import matplotlib.pyplot as plt
%matplotlib inline
import numpy as np
from PyEMD import EEMD
import scipy.io.wavfile as wf
plt.rcParams['font.sans-serif']=['SimHei'] #用来正常显示中文标签
plt.rcParams['axes.unicode_minus']=False #用来正常显示符号
Fs = 44100
T = 1/Fs
sample_rate, noised_sigs = wf.read('noised.wav')
t = np.linspace(0,len(noised_sigs)-1,len(noised_sigs))
s = np.array(noised_sigs)
emd = EEMD()
imfs = emd.emd(s,t,10)
# Plot results
plt.figure(figsize=(12,15))
plt.subplot(7,2,1)
plt.plot(t, s, 'r')
plt.title("Input signal")
plt.figure(figsize=(12,15))
plt.subplot(7,2,1)
plt.plot(t, s, 'r')
plt.title("IMFS")
for n, imf in enumerate(imfs):
plt.subplot(13,2,n 3)
plt.plot(t[:2000], imf[:2000], 'g')
plt.ylabel("IMF %i" %(n 1))
plt.tight_layout()
plt.show()
测试数据(noised.wav)放在这里,方便大家复现:
链接:https://pan.百度.com/s/1_nMaikrjGYDmae8u1RYIVA
提取码:z55v
这篇好文章是转载于:学新通技术网
- 版权申明: 本站部分内容来自互联网,仅供学习及演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系,请提供相关证据及您的身份证明,我们将在收到邮件后48小时内删除。
- 本站站名: 学新通技术网
- 本文地址: /boutique/detail/tanhgagbah
系列文章
更多
同类精品
更多
-
photoshop保存的图片太大微信发不了怎么办
PHP中文网 06-15 -
《学习通》视频自动暂停处理方法
HelloWorld317 07-05 -
word里面弄一个表格后上面的标题会跑到下面怎么办
PHP中文网 06-20 -
Android 11 保存文件到外部存储,并分享文件
Luke 10-12 -
photoshop扩展功能面板显示灰色怎么办
PHP中文网 06-14 -
微信公众号没有声音提示怎么办
PHP中文网 03-31 -
excel下划线不显示怎么办
PHP中文网 06-23 -
excel打印预览压线压字怎么办
PHP中文网 06-22 -
TikTok加速器哪个好免费的TK加速器推荐
TK小达人 10-01 -
怎样阻止微信小程序自动打开
PHP中文网 06-13