Python 时间序列预测：Hot-winters-轻识

1. 什么是Holt-Winters

时间序列是非常常见的数据格式，以[时间，观测值]形式表现，如下图。

现实场景中如股票走势图，国家GDP历年数据，机器cpu利用率，内存数据等都是时间序列。对未来时间的观测值进行预测是有意义的工作，提前预知未来的数据的走势，可以提前做出行动，如预测cpu使用率，如果cpu飙高，可以及早进行调整，避免机器负载过高而宕机，这个在AIOPS是很常见的一个应用场景。

今天要说到Holt-Winters是利用三次指数平滑来做时间序列预测的方法。Holt-Winters是综合了1957年Holt和1960年Winters两个人的思路的一种方法。

一次指数平滑
我们来看下，一次指数平滑如下图：

可知，si表示第i时刻的平滑估计，si可以表示为当前实际值xi和上一时刻平滑估计值得加权组合，权重由alpha来决定。那为什么称为指数平滑呢？我们来把式子展开，如下：

有点类似泰勒展开式的味道

alpha 属于[0, 1], 越大表示近期的数据影响更大

二次指数平滑：加上趋势的因素
一次指数平滑，没有考虑时间序列的趋势和季节性，二次指数平滑加上趋势因素。

从公式可知，一个时间序列的时刻值分解为baseline部分和趋势部分，t表示趋势，可以表示为连续两个时刻的差值；可知，ti也是一次的指数平滑。

Holt-Winters三次指数平滑：加上季节性因素
在二次指数平滑基础上，考虑季节性因素，就是三次指数平滑，也就是Holt-Winters。由此，一个时间序列的时刻值分解为baseline部分和趋势部分以及季节部分。由于季节性，存在周期，比如按周，按月等。pi季节性为当前季节性值和上一个周期季节性估计值的加权组合，周期在公式中以k来表示。如下：

2. Holt-Winters的实现

从第一部分可知，要实现Holt-Winters，只要知道：

•初始值：s0，t0和p0•合适的参数：alpha，beta， gamma•套入公式即可完成预测

三个重要参数：alpha，beta， gamma都属于[0, 1]之间，要么人为的搜索，要么通过数据来估计，通常采用L-BFGS优化算法来拟合数据。优化算法来自包scipy.optimize的fmin_l_bfgs_b。

from __future__ import divisionfrom sys importexitfrom math import sqrtfrom numpy import arrayfrom scipy.optimize import fmin_l_bfgs_b# 优化算法的loss function，即判断拟合效果，由RMSE MAE等def RMSE(params, *args):    Y = args[0]    type = args[1]    rmse = 0    alpha, beta, gamma = params    m = args[2]         a = [sum(Y[0:m]) / float(m)]    b = [(sum(Y[m:2* m]) - sum(Y[0:m])) / m ** 2]if type == 'additive':        s = [Y[i] - a[0] for i in range(m)]        y = [a[0] + b[0] + s[0]]for i in range(len(Y)):            a.append(alpha * (Y[i] - s[i]) + (1- alpha) * (a[i] + b[i]))            b.append(beta * (a[i + 1] - a[i]) + (1- beta) * b[i])            s.append(gamma * (Y[i] - a[i] - b[i]) + (1- gamma) * s[i])            y.append(a[i + 1] + b[i + 1] + s[i + 1])    rmse = sqrt(sum([(m - n) ** 2for m, n in zip(Y, y[:-1])]) / len(Y))return rmse# 加性的时间序列def additive(x, m, fc, alpha = None, beta = None, gamma = None):    Y = x[:]# 利用fmin_l_bfgs_b来估计参数alpha beta和gammaif(alpha == Noneor beta == Noneor gamma == None):        initial_values = array([0.3, 0.1, 0.1])        boundaries = [(0, 1), (0, 1), (0, 1)]        type = 'additive'        parameters = fmin_l_bfgs_b(RMSE, x0 = initial_values, args = (Y, type, m), bounds = boundaries, approx_grad = True)        alpha, beta, gamma = parameters[0]# 初始值 a表示baseline， b表示趋势，s表示季节性，y表示预测值， 分别取第一个周期的统计数据为初始值    a = [sum(Y[0:m]) / float(m)]    b = [(sum(Y[m:2* m]) - sum(Y[0:m])) / m ** 2]    s = [Y[i] - a[0] for i in range(m)]    y = [a[0] + b[0] + s[0]]    rmse = 0# 套用上面公式，从0开始，fc表示预测的数量，如已知前7天，预测接下来的一个小时的数据，如果数据粒度是5分钟，fc为12。for i in range(len(Y) + fc):if i == len(Y):# 预测值为            Y.append(a[-1] + b[-1] + s[-m])        a.append(alpha * (Y[i] - s[i]) + (1- alpha) * (a[i] + b[i]))        b.append(beta * (a[i + 1] - a[i]) + (1- beta) * b[i])        s.append(gamma * (Y[i] - a[i] - b[i]) + (1- gamma) * s[i])        y.append(a[i + 1] + b[i + 1] + s[i + 1])# 计算rmse值    rmse = sqrt(sum([(m - n) ** 2for m, n in zip(Y[:-fc], y[:-fc - 1])]) / len(Y[:-fc]))return y[-fc:], alpha, beta, gamma, rmse

另外，statsmodels包中也提供的实现的方法

from statsmodels.tsa.holtwinters importExponentialSmoothing

3. Holt-Winters参数

从上面实现可知，holt-winters通过预估alpha，beta和gamma来预测。算法的关键就是这三个参数和初始化值。三个参数可以通过优化算法来预估，但有可能并不是最优的。初始值的设置除了上面统计值外，还可以通过时序的分解的趋势和季节部分来初始。

import numpy as npfrom pandas import read_csvimport matplotlib.pyplot as pltfrom statsmodels.tsa.seasonal import seasonal_decomposedecomposition = seasonal_decompose(df_clean.bw, model='additive', period=288)decomposition.plot()

Holt-Winters针对波形比较稳定，没有突刺的情况下，效果会比较好。

对于存在突刺，统一的alpha，beta，gamma不能很好拟合，预测可能会滞后。

4. 总结

本文分享了时间序列预测算法Holt-Winters以及重要参数的选择，希望对你有帮助。总结如下：

•Holt-Winters是三次指数平滑，分别为baseline，趋势和季节性；•alpha、beta和gamma分别为baseline，趋势和季节性的指数加权参数，一般通过优化算法L-BFGS估计•初始化可通过平均值，也可通过时间序列分解得到•周期m或者k的选择要根据实际数据来选择•Holt-Winters针对波形比较稳定，没有突刺的情况下，效果会比较好

作者简介：wedo实验君, 数据分析师；热爱生活，热爱写作

赞赏作者

点击下方阅读原文加入社区会员