python時間序列分解STL

參考鏈接：參考鏈接1，參考鏈接2.

輸入數據

 代碼：

import statsmodels.api as sm
import matplotlib.pyplot as plt
import pandas as pd
#from dateutil import get_gran, format_timestamp

data = pd.read_csv('data.csv',usecols=['time', 'value'],index_col='time')#將時間列作爲索引
print(data)
rd = sm.tsa.seasonal_decompose(data['value'].values, freq=72)#該時間序列週期爲72
rd.plot()
plt.show()

結果：