#-*-?coding:?utf-8?-*-

import?math
import?numpy?as?np
import?pylab?as?pl
import?matplotlib.pyplot?as?plt
import?scipy.signal?as?signal
from?scipy?import?fftpack
from?scipy.fftpack?import?fft?ifft
import?scipy.signal?as?signal
from?scipy?import?interpolate


#?判定當(dāng)前的時(shí)間序列是否是單調(diào)序列
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


#?尋找當(dāng)前時(shí)間序列的極值點(diǎn)
def?findpeaks（x）:
????return?signal.argrelextrema（x?np.greater）[0]

#?判斷當(dāng)前的序列是否為?IMF?序列
def?isImf（x）:
????N?=?np.size（x）
????pass_zero?=?np.sum（x[0:N?-?2]?*?x[1:N?-?1]?????peaks_num?=?np.size（findpeaks（x））?+?np.size（findpeaks（-x））??#?極值點(diǎn)的個(gè)數(shù)
????if?abs（pass_zero?-?peaks_num）?>?1:
????????return?False
????else:
????????return?True


#?獲取當(dāng)前樣條曲線
def?getspline（x）:
????N?=?np.size（x）
????peaks?=?findpeaks（x）
????print?（‘當(dāng)前極值點(diǎn)個(gè)數(shù)：‘?len（peaks））
????if?（len（peaks）?<=?3）:
????????if?（len（peaks）?????????????peaks?=?np.concatenate（（[0]?peaks））
????????????peaks?=?np.concatenate（（peaks?[N?-?1]））??#?這里是為了防止樣條次數(shù)不夠，無(wú)法插值的情況
????????t?=?interpolate.splrep（peaks?y=x[peaks]?w=None?xb=None?xe=None?k=len（peaks）?-?1）
????????return?interpolate.splev（np.arange（N）?t）
????t?=?interpolate.splr