Python - EMD分解

鄭重聲明：本文檔只是方便自己學習記錄

1.EMD 工具包安裝

下載地址：https://github.com/laszukdawid/PyEMD

2.解壓工具包，將文件復制到自己的python（Anaconda）的Lib的site-packages

3.cmd切換到包的目錄

4.輸入python setup.py install安裝

5.EMD分解實驗

# 導入工具庫
import numpy as np
from PyEMD import EMD, Visualisation
# 構建信號
t = np.arange(0,1, 0.01)
S = 2*np.sin(2*np.pi*15*t) +4*np.sin(2*np.pi*10*t)*np.sin(2*np.pi*t*0.1)+np.sin(2*np.pi*5*t)
# 提取imfs和剩余信號res
emd = EMD()
emd.emd(S)
imfs, res = emd.get_imfs_and_residue()
# 繪制 IMF
vis = Visualisation()
vis.plot_imfs(imfs=imfs, residue=res, t=t, include_residue=True)
# 繪制並顯示所有提供的IMF的瞬時頻率
vis.plot_instant_freq(t, imfs=imfs)
vis.show()