安装包
pyhht
pip install pyhht
from pyhht.emd import EMD
from pyhht.visualization import plot_imfs
emd = EMD(data.RUL[:10000])
imfs = emd.decompose()
PyEMD
github地址
注:作者将安装包名写错,安装后将文件名改为 PyEMD即可。
pip install EMD-signal
经验模态分解python初步实现
import math
import numpy as np
import pylab as pl
import matplotlib.pyplot as plt
import scipy.signal as signal
from scipy import fftpack
import scipy.signal as signal
from scipy import interpolate
#判定当前的时间序列是否是单调序列
def ismonotonic(x):
max_peaks=signal.argrelextrema(x,np.greater)[0]
min_peaks=signal.argrelextrema(x,np.less)[0]
all_num=len(max_peaks)+len(min_peaks)
if all_num>0:
return False
else:
return True
#寻找当前时间序列的极值点
def findpeaks(x):
return signal.argrelextrema(x,np.greater)[0]
#判断当前的序列是否为 IMF 序列
def isImf(x):
N=np.size(x)
pass_zero=np.sum(x[0:N-2]*x[1:N-1]<0)#过零点的个数
peaks_num=np.size(findpeaks(x))+np.size(findpeaks(-x))#极值点的个数
if abs(pass_zero-peaks_num)>1:
return False
else:
return True
#获取当前样条曲线
def getspline(x):
N=np.size(x)
peaks=findpeaks(x)
print '当前极值点个数:',