11种常见的时间序列预测方法

article/2025/8/27 6:12:57

参考内容：4大类11种常见的时间序列预测方法总结和代码示例

代码地址：
https://github.com/SeafyLiang/machine_learning_study/blob/master/time_series

11种常见的时间序列预测方法

    • 1、指数平滑Exponential Smoothing
    • 2、Holt-Winters 法
    • 3、自回归 (AR)
    • 4、移动平均模型(MA)
    • 5、自回归滑动平均模型 (ARMA)
    • 6、差分整合移动平均自回归模型 (ARIMA)
    • 7、季节性 ARIMA (SARIMA)
    • 8、包含外生变量的SARIMA (SARIMAX)
    • 9、向量自回归 (VAR)
    • 10、向量自回归滑动平均模型 (VARMA)
    • 11、包含外生变量的向量自回归滑动平均模型 (VARMAX)
    • 总结

本篇文章将总结时间序列预测方法，并将所有方法分类介绍并提供相应的python代码示例，以下是本文将要介绍的方法列表：

1、使用平滑技术进行时间序列预测

  • 指数平滑
  • Holt-Winters 法

2、单变量时间序列预测

  • 自回归 (AR)
  • 移动平均模型 (MA)
  • 自回归滑动平均模型 (ARMA)
  • 差分整合移动平均自回归模型 (ARIMA)
  • 季节性 ARIMA (SARIMA)

3、外生变量的时间序列预测

  • 包含外生变量的SARIMAX (SARIMAX)
  • 具有外生回归量的向量自回归移动平均 (VARMAX)

4、多元时间序列预测

  • 向量自回归 (VAR)
  • 向量自回归移动平均 (VARMA)

下面我们对上面的方法一一进行介绍,并给出python的代码示例
在这里插入图片描述

1、指数平滑Exponential Smoothing

指数平滑法是过去观测值的加权平均值，随着观测值变老，权重呈指数会衰减。换句话说，观察时间越近相关权重就越高。它可以快速生成可靠的预测，并且适用于广泛的时间序列。
简单指数平滑：此方法适用于预测没有明确趋势或季节性模式的单变量时间序列数据。简单指数平滑法将下一个时间步建模为先前时间步的观测值的指数加权线性函数。
它需要一个称为 alpha (a) 的参数，也称为平滑因子或平滑系数，它控制先前时间步长的观测值的影响呈指数衰减的速率，即控制权重减小的速率。a 通常设置为 0 和 1 之间的值。较大的值意味着模型主要关注最近的过去观察，而较小的值意味着在进行预测时会考虑更多的历史。简单指数平滑时间序列的简单数学解释如下所示：
在这里插入图片描述

# SES
from statsmodels.tsa.holtwinters import SimpleExpSmoothing
from random import random
# contrived dataset
data = [x + random() for x in range(1, 100)]
# fit model
model = SimpleExpSmoothing(data)
model_fit = model.fit()
# make prediction
yhat = model_fit.predict(len(data), len(data))
print(yhat)

2、Holt-Winters 法

在 1957 年初，Holt扩展了简单的指数平滑法，使它可以预测具有趋势的数据。这种被称为 Holt 线性趋势的方法包括一个预测方程和两个平滑方程（一个用于水平，一个用于趋势）以及相应的平滑参数 α 和 β。后来为了避免趋势模式无限重复，引入了阻尼趋势法，当需要预测许多序列时，它被证明是非常成功和最受欢迎的单个方法。除了两个平滑参数之外，它还包括一个称为阻尼参数 φ 的附加参数。
一旦能够捕捉到趋势，Holt-Winters 法扩展了传统的Holt法来捕捉季节性。Holt-Winters 的季节性方法包括预测方程和三个平滑方程——一个用于水平，一个用于趋势，一个用于季节性分量，并具有相应的平滑参数 α、β 和 γ。
此方法有两种变体，它们在季节性成分的性质上有所不同。当季节变化在整个系列中大致恒定时，首选加法方法，而当季节变化与系列水平成比例变化时，首选乘法方法。

# HWES
from statsmodels.tsa.holtwinters import ExponentialSmoothing
from random import random
# contrived dataset
data = [x + random() for x in range(1, 100)]
# fit model
model = ExponentialSmoothing(data)
model_fit = model.fit()
# make prediction
yhat = model_fit.predict(len(data), len(data))
print(yhat)

3、自回归 (AR)

在 AR 模型中，我们使用变量过去值的线性组合来预测感兴趣的变量。术语自回归表明它是变量对自身的回归。AR模型的简单数学表示如下：
在这里插入图片描述

这里，εt 是白噪声。这类似于多元回归，但是使用 yt 的滞后值作为预测变量。我们将其称为 AR§ 模型，即 p 阶的自回归模型。

#AR
from statsmodels.tsa.ar_model import AutoReg
from random import random
# contrived dataset
data = [x + random() for x in range(1, 100)]
# fit model
model = AutoReg(data, lags=1)
model_fit = model.fit()
# make prediction
yhat = model_fit.predict(len(data), len(data))
print(yhat)

4、移动平均模型(MA)

与在回归中使用预测变量的过去值的 AR 模型不同，MA 模型在类似回归的模型中关注过去的预测误差或残差。MA模型的简单数学表示如下：
在这里插入图片描述

这里，εt 是白噪声。我们将其称为 MA(q) 模型，即 q 阶移动平均模型。

# MA
from statsmodels.tsa.arima.model import ARIMA
from random import random
# contrived dataset
data = [x + random() for x in range(1, 100)]
# fit model
model = ARIMA(data, order=(0, 0, 1))
model_fit = model.fit()
# make prediction
yhat = model_fit.predict(len(data), len(data))
print(yhat)

需要说明的是不应将这里说的移动平均线方法与计算时间序列的移动平均线混淆，因为两者是不同的概念。

5、自回归滑动平均模型 (ARMA)

在 AR 模型中，我们使用变量过去值与过去预测误差或残差的线性组合来预测感兴趣的变量。它结合了自回归 (AR) 和移动平均 (MA) 模型。
AR 部分涉及对变量自身的滞后（即过去）值进行回归。MA部分涉及将误差项建模为在过去不同时间同时发生的误差项的线性组合。模型的符号涉及将 AR§ 和 MA(q) 模型的顺序指定为 ARMA 函数的参数，例如 ARMA（p，q）。ARMA 模型的简单数学表示如下所示：
在这里插入图片描述

# ARMA 
from statsmodels.tsa.arima.model import ARIMA
from random import random
# contrived dataset
data = [random() for x in range(1, 100)]
# fit model
model = ARIMA(data, order=(2, 0, 1))
model_fit = model.fit()
# make prediction
yhat = model_fit.predict(len(data), len(data))
print(yhat)

6、差分整合移动平均自回归模型 (ARIMA)

如果我们将差分与自回归和移动平均模型相结合，我们将获得 ARIMA 模型。ARIMA 是差分整合移动平均自回归模型Autoregressive Integrated Moving Average model 的首字母缩写。它结合了自回归 (AR) 和移动平均模型 (MA) 以及为了使序列平稳而对序列的差分预处理过程，这个过程称为积分(I)。ARIMA 模型的简单数学表示如下：
在这里插入图片描述

其中 y′t 是差分级数。右侧的“预测变量”包括滞后值和滞后误差。我们称之为 ARIMA(p,d,q) 模型。

这里，p 是自回归部分的阶数，d 是所涉及的一阶差分程度，q 是移动平均部分的阶数。
ACF 和 PACF 图在求 p 和 q 阶中的意义：

  • 为了找到 AR§ 模型的阶 p：我们预计 ACF 图会逐渐减小，同时 PACF 在 p 显著滞后后会急剧下降或切断。
  • 为了找到 MA(q) 模型的阶 p：我们预计 PACF 图将逐渐减小，同时 ACF 应该在某些 q 显著滞后后急剧下降或切断。
# ARIMA 
from statsmodels.tsa.arima.model import ARIMA
from random import random
# contrived dataset
data = [x + random() for x in range(1, 100)]
# fit model
model = ARIMA(data, order=(1, 1, 1))
model_fit = model.fit()
# make prediction
yhat = model_fit.predict(len(data), len(data), typ='levels')
print(yhat)

7、季节性 ARIMA (SARIMA)

ARIMA 模型还能够对广泛的季节性数据进行建模。季节性 ARIMA 模型是通过在 ARIMA 模型中包含额外的季节性项来形成的。
在这里插入图片描述

这里，m = 每个时间季节的步数。我们对模型的季节性部分使用大写符号，对模型的非季节性部分使用小写符号。

它将 ARIMA 模型与在季节性数据级别执行相同的自回归、差分和移动平均建模的能力相结合。

# SARIMA 
from statsmodels.tsa.statespace.sarimax import SARIMAX
from random import random
# contrived dataset
data = [x + random() for x in range(1, 100)]
# fit model
model = SARIMAX(data, order=(1, 1, 1), seasonal_order=(0, 0, 0, 0))
model_fit = model.fit(disp=False)
# make prediction
yhat = model_fit.predict(len(data), len(data))
print(yhat)

8、包含外生变量的SARIMA (SARIMAX)

SARIMAX 模型是传统 SARIMA 模型的扩展，包括外生变量的建模，是Seasonal Autoregressive Integrated Moving-Average with Exogenous Regressors 的缩写
外生变量是其值在模型之外确定并施加在模型上的变量。它们也被称为协变量。外生变量的观测值在每个时间步直接包含在模型中，并且与主要内生序列的使用不同的建模方式。
SARIMAX 方法也可用于通过包含外生变量来模拟具有外生变量的其他变化，例如 ARX、MAX、ARMAX 和 ARIMAX。

# SARIMAX 
from statsmodels.tsa.statespace.sarimax import SARIMAX
from random import random
# contrived dataset
data1 = [x + random() for x in range(1, 100)]
data2 = [x + random() for x in range(101, 200)]
# fit model
model = SARIMAX(data1, exog=data2, order=(1, 1, 1), seasonal_order=(0, 0, 0, 0))
model_fit = model.fit(disp=False)
# make prediction
exog2 = [200 + random()]
yhat = model_fit.predict(len(data1), len(data1), exog=[exog2])
print(yhat)

9、向量自回归 (VAR)

VAR 模型是单变量自回归模型的推广，用于预测时间序列向量或多个并行时间序列，例如 多元时间序列。它是关于系统中每个变量的一个方程。
如果序列是平稳的，可以通过将 VAR 直接拟合到数据来预测它们（称为“VAR in levels”）。如果序列是非平稳的，我们会取数据的差异以使其平稳，然后拟合 VAR 模型（称为“VAR in differences”）。
我们将其称为 VAR§ 模型，即 p 阶向量自回归模型。

# VAR 
from statsmodels.tsa.vector_ar.var_model import VAR
from random import random
# contrived dataset with dependency
data = list()
for i in range(100):v1 = i + random()v2 = v1 + random()row = [v1, v2]data.append(row)
# fit model
model = VAR(data)
model_fit = model.fit()
# make prediction
yhat = model_fit.forecast(model_fit.y, steps=1)
print(yhat)

10、向量自回归滑动平均模型 (VARMA)

VARMA 方法是 ARMA 对多个并行时间序列的推广，例如 多元时间序列。具有有限阶 MA 误差项的有限阶 VAR 过程称为 VARMA。
模型的公式将 AR§ 和 MA(q) 模型的阶数指定为 VARMA 函数的参数，例如 VARMA（p，q）。VARMA 模型也可用于VAR 或 VMA 模型。

# VARMA 
from statsmodels.tsa.statespace.varmax import VARMAX
from random import random
# contrived dataset with dependency
data = list()
for i in range(100):v1 = random()v2 = v1 + random()row = [v1, v2]data.append(row)
# fit model
model = VARMAX(data, order=(1, 1))
model_fit = model.fit(disp=False)
# make prediction
yhat = model_fit.forecast()
print(yhat)

11、包含外生变量的向量自回归滑动平均模型 (VARMAX)

Vector Autoregression Moving-Average with Exogenous Regressors (VARMAX) 是 VARMA 模型的扩展，模型中还包含使用外生变量的建模。它是 ARMAX 方法对多个并行时间序列的推广，即 ARMAX 方法的多变量版本。
VARMAX 方法也可用于对包含外生变量的包含模型进行建模，例如 VARX 和 VMAX。

# VARMAX 
from statsmodels.tsa.statespace.varmax import VARMAX
from random import random
# contrived dataset with dependency
data = list()
for i in range(100):v1 = random()v2 = v1 + random()row = [v1, v2]data.append(row)
data_exog = [x + random() for x in range(100)]
# fit model
model = VARMAX(data, exog=data_exog, order=(1, 1))
model_fit = model.fit(disp=False)
# make prediction
data_exog2 = [[100]]
yhat = model_fit.forecast(exog=data_exog2)
print(yhat)

总结

在这篇文章中，基本上覆盖了所有主要时间序列预测的问题。我们可以把上面提到的方法整理成以下几个重要的方向：

  • AR：自回归
  • MA：平均移动
  • I：差分整合
  • S：季节性
  • V：向量（多维输入）
  • X：外生变量

本文中提到的每种算法基本上都是这几种方法的组合，本文中已将每种的算法都进行了重点的描述和代码的演示，如果你想深入了解其中的知识请查看相关的论文。


http://chatgpt.dhexx.cn/article/wpet6FUT.shtml

相关文章

时间序列预测方法最全总结！

时间序列预测方法最全总结！

时间序列预测就是利用过去一段时间的数据来预测未来一段时间内的信息，包括连续型预测（数值预测，范围估计）与离散型预测（事件预测）等，具有非常高的商业价值。 需要明确一点的是，与回归…
阅读更多...
时间序列（一）：时间序列数据与时间序列预测模型

时间序列（一）：时间序列数据与时间序列预测模型

时间序列系列文章： 时间序列（一）：时间序列数据与时间序列预测模型 时间序列（二）：时间序列平稳性检测 时间序列（三）：ARIMA模型实战 时间序列及其预测是日常工…
阅读更多...
组合预测模型 | ARIMA-LSTM时间序列预测（Python）

组合预测模型 | ARIMA-LSTM时间序列预测（Python）

组合预测模型 | ARIMA-LSTM时间序列预测（Python） 目录 组合预测模型 | ARIMA-LSTM时间序列预测（Python）预测结果基本介绍程序设计参考资料 预测结果 基本介绍 ARIMA-LSTM时间序列预测（Python完整源码和数据&#xff09…
阅读更多...
时间序列预测系列文章总结（代码使用方法）

时间序列预测系列文章总结（代码使用方法）

前言 这篇文章是对前面所写的LSTM时序预测文章的代码使用方法的总结。强烈建议使用代码前先阅读本文，而不是直接询问！ 此外，代码数据中除了负荷列其他列都已经归一化了，所以代码中只对负荷列进行了归一化，如果使用自己…
阅读更多...
4大类11种常见的时间序列预测方法总结和代码示例

4大类11种常见的时间序列预测方法总结和代码示例

本篇文章将总结时间序列预测方法，并将所有方法分类介绍并提供相应的python代码示例，以下是本文将要介绍的方法列表： 1、使用平滑技术进行时间序列预测 指数平滑Holt-Winters 法 2、单变量时间序列预测 自回归 (AR)移动平均模型 (MA)自回归…
阅读更多...
【深度学习时间序列预测案例】零基础入门经典深度学习时间序列预测项目实战（附代码+数据集+原理介绍）

【深度学习时间序列预测案例】零基础入门经典深度学习时间序列预测项目实战（附代码+数据集+原理介绍）

🚨注意🚨：最近经粉丝反馈，发现有些订阅者将此专栏内容进行二次售卖，特在此声明，本专栏内容仅供学习，不得以任何方式进行售卖，未经作者许可不得对本专栏内容行使发表权、署名权、修改…
阅读更多...
如何理解vcc,vdd,vss

如何理解vcc,vdd,vss

常见解释 VCC： Ccircuit 表示电路的意思, 即接入电路的电压 VDD：Ddevice 表示器件的意思, 即器件内部的工作电压 VSS： Sseries 表示公共连接的意思，通常指电路公共接地端电压 个人理解 VCC ：双极器件的正，一…
阅读更多...
VCC、VDD、VSS、VDDA、VSSA、VREF+等标识

VCC、VDD、VSS、VDDA、VSSA、VREF+等标识

阅读更多...
VCC,GND,VSS,VDD的理解

VCC,GND,VSS,VDD的理解

Definition： VCC：Ccircuit 表示电路的意思, 即接入电路的电压。 GND：在电路里常被定为电压参考基点。 VDD：Ddevice 表示器件的意思, 即器件内部的工作电压 VSS：Sseries 表示公共连接的意思，通常指电路公共接…
阅读更多...
VCS简介

VCS简介

1.2.1 关于VCS VCS是Verilog Compiled Simulator的缩写。VCS MX是一个编译型的代码仿真器。它使你能够分析，编译和仿真Verilog，VHDL，混合HDL，SystemVerilog，OpenVera和SystemC描述的设计。 它还为您提供了一系列仿真和…
阅读更多...
VCC、VDD、VSS以及VBAT的区别

VCC、VDD、VSS以及VBAT的区别

在STM32 的学习中，发现有几种看起来相关的名称，分别是VCC、VDD、VSS、VBAT，在经过搜索查找之后，总结如下： 1.VCC的C是Circuit的意思，是指整个供电回路的电压, 也有人说VCC是双极器件的正极 2.VDD的D是Dev…
阅读更多...
Vcc、Vee、Vdd、Vss傻傻分不清楚？

Vcc、Vee、Vdd、Vss傻傻分不清楚？

Vcc、Vee、Vdd、Vss傻傻分不清楚？ 以下内容源于网络。 电源 左边两个符号均表示电池类直流电源，左边第一个表示电池组，第二个表示单电池。这些符号不一定总是带有 和 - 号。左边第三个（带 和 - 的圆圈符号）表示非电…
阅读更多...
semi-supervised classification

semi-supervised classification

半监督学习任务主要分为半监督聚类、半监督分类、半监督回归等问题，我们主要针对半监督分类问题。 半监督学习的假设 基于低密度假设 模型的决策边界不应该将该密度区域划分开，而应该处在低密度区域内。基于平滑化假设 输入空间中&#xf…
阅读更多...
Weakly-Supervised Semantic Segmentation via Sub-category Exploration

Weakly-Supervised Semantic Segmentation via Sub-category Exploration

Weakly-Supervised Semantic Segmentation via Sub-category Exploration Abstract问题 1.Introduction2. Related Work2.1weakly-supervised semantic segmentation(WSSS)弱监督语义分割的初始预测。---Initial Prediction for WSSS.2.2 、Response Refinement for WSSS&#…
阅读更多...
Scribble-Supervised Medical Image Segmentation

Scribble-Supervised Medical Image Segmentation

Scribble-Supervised Medical Image Segmentation via Dual-Branch Network and Dynamically Mixed Pseudo Labels Supervision 摘要 采用草率的分割标注心脏分割双分枝网络：一个编码器，两个解码器动态结合两个解码器的输出伪标签 方法 模型结构 Lpce…
阅读更多...
Learning Affinity from Attention End-to-End Weakly-Supervised Semantic Segmentation withTransformers

Learning Affinity from Attention End-to-End Weakly-Supervised Semantic Segmentation withTransformers

Learning Affinity from Attention: End-to-End Weakly-Supervised Semantic Segmentation with Transformers 来源：CVPR 2022，武汉大学、京东、悉尼大学 导言 本文是一篇做自然图像弱监督语义分割的论文，利用图像级的类别标签来实现像素级…
阅读更多...
Token Contrast for Weakly-Supervised Semantic Segmentation

Token Contrast for Weakly-Supervised Semantic Segmentation

文章来源：[CVPR2023] Keywords：Weakly-Supervised Semantic Segmentation（WSSS）；over-smoothing; ViT 一、本文提出的问题以及解决方案: 本文解决了over-smoothing问题，该问题其实是在之前的GCN网络中提出…
阅读更多...
Self-supervised Learning整理

Self-supervised Learning整理

Self-supervised Learning Pre-train Fine-tune Pre-train Fine-tune算是一种Transfer Learning。 首先，假设按照SimCLR中的设定，将一般的模型分为两部分，分别为Encoder和Projection Head。Encoder结构类似于AutoEncoder中的Encoder部分&a…
阅读更多...
Supervised Discrete Hashing

Supervised Discrete Hashing

Supervised Discrete Hashing 2015 CVPR 问题： 处理施加在追踪的哈希码上的离散约束，使哈希优化具有挑战性(通常是NP- hard)。 解决： 提出了一个新的监督哈希框架，其中的学习目标是生成最优的二进制哈希码用于线性分类。 通过引…
阅读更多...
NetVLAD: CNN architecture for weakly supervised place recognition

NetVLAD: CNN architecture for weakly supervised place recognition

背景知识： Vector of Locally Aggregated Descriptors（VLAD）image retrieval. 【CC】是广泛使用的图像提取方式，本文是在在这个提取器上做改进；具体是啥下面有介绍 weakly supervised ranking loss 【CC】本文的另外…
阅读更多...
推荐文章

Copyright @ 2022~2023