代码来源于
视频讲解
做工业增速的时间序列分解
import pandas as pd
from matplotlib import pyplot as plt
from statsmodels.tsa.seasonal import seasonal_decompose
# check time series data
# air_pollution = pd.read_csv('../dataset/air_pollution.csv', parse_dates=['date'])
air_pollution = pd.read_csv('C:/Users/SMXzh/Desktop/try/Y2005M1-2022M33.csv', parse_dates=['Date'], encoding='utf-8')
# air_pollution = air_pollution.dropna()
air_pollution.set_index('Date', inplace=True)
# print(air_pollution.head())
values = air_pollution.values
# groups = [0, 1, 2, 3, 4, 5, 6, 7]
groups = [0]
# i = 1
# # plot each column
# for group i