基于Matlab的时间序列（Time Series）（附代码）

一、模型介绍

时间序列（或称动态数列）是指将同一统计指标的数值按其发生的时间先后顺序排列而成的数列。时间序列分析的主要目的是根据已有的历史数据对未来进行预测。经济数据中大多数以时间序列的形式给出。根据观察时间的不同，时间序列中的时间可以是年份、季度、月份或其他任何时间形式。

1.1. 时间序列的不同分类

1. 按所研究的对象的多少分，有一元时间序列和多元时间序列。

2. 按时间的连续性可将时间序列分为离散时间序列和连续时间序列两种。

3. 按时间序列的分布规律来分，有高斯型时间序列和非高斯型时间序列。

4. 按序列的统计特性分，有平稳时间序列和非平稳时间序列。

   • 严平稳时间序列，也叫狭义平稳时间序列：一个时间序列的概率分布与时间$t$无关
   • 宽平稳时间序列，也叫广义平稳时间序列：一个时间序列的 一、二阶矩存在，而且对任意时刻$t$满足：
     1. 均值为常数
     2. 协方差为时间间隔$\tau$的函数。

我们以后所研究的时间序列主要是宽平稳时间序列。

1.2. 时间序列构成要素

时间序列是按照时间排序的一组随机变量，它通常是在相等间隔的时间段内依照给定的采样率对某种潜在过程进行观测的结果。时间序列数据本质上反映的是某个或者某些随机变量随时间不断变化的趋势，而时间序列预测方法的核心就是从数据中挖掘出这种规律，并利用其对将来的数据做出估计。
构成要素：长期趋势，季节变动，循环变动，不规则变动。

• 长期趋势（ T ）现象在较长时期内受某种根本性因素作用而形成的总的变动趋势。
• 季节变动（ S ）现象在一年内随着季节的变化而发生的有规律的周期性变动。
• 循环变动（ C ）现象以若干年为周期所呈现出的波浪起伏形态的有规律的变动。
• 不规则变动（I ）是一种无规律可循的变动，包括严格的随机变动和不规则的突发性影响很大的变动两种类型。

1.3. 三种时间序列模型

通常用$T_t$表示长期趋势项，$S_t$来表示季节变动趋势项，$C_t$来表示循环变动趋势，$I_t$来表示随机干扰项。常见的确定性时间序列模型有以下几种类型：
（1）加法模型$$x_t=T_t+S_t+C_t+I_t$$
（2）乘法模型$$x_t=T_t\ast S_t\ast C_t\ast I_t$$
（3）混合模型$$x_t=T_t\ast S_t+I_t$$$$x_t=S_t+T_t\ast C_t\ast I_t$$
其中$x_t$是观测目标的观测记录，$E(R_t)=0,E(R_t^2)={\sigma }^2$。

如果在预测时间范围以内，无突然变动且随机变动的方差${\sigma }^2$较小，并且有理由认为过去和现在的演变趋势将继续发展到未来，也就是假设研究数据符合宽平稳时间序列时，$S_t+C_t=0$，对$x_t=T_t+I_t$项的建模方式的不同也就构成了时间序列预测的三种不同模型：

1.3.1. AR(p)模型

AR(p)模型，表示时间$t$下的$x_t$由过去的$t-1,t-2,...,t-p$时间下的$x_{t-1},x_{t-2},...,x_{t-p}$决定，称之为p阶自回归（Auto Regression）过程：
$$x_t={\varphi }_0+{\varphi }_1{x}_{t-1}+{\varphi }_2{x}_{t-2}+\cdots +{\varphi }_p{x}_{t-p}+{\epsilon }_t(1)$$
数字特征：

1. $E(x_t)=E(x_{t-1})=\cdots =E(x_p)$
2. $D(x_t)=D(x_{t-1})=\cdots =D(x_p)$
3. $E({\epsilon }_t)=0,D({\epsilon }_t)={\sigma }_{{\epsilon }_t}^2$

结合以上的数字特征，藉由式（1）变形可得知：

1. $E(x_t)=\frac{{\varphi }_0}{1-{\varphi }_1-\cdots {\varphi }_p}$ 与p相关
2. $D(x_t)=\frac{{\sigma }_{{\epsilon }_t}^2}{1-{\varphi }_1^2}$

变形过程推导，自相关系数${\rho }_k$，协方差$\mathrm{r}(\mathrm{k})=\mathrm{c}\mathrm{o}\mathrm{v}(x_t,x_{t-k})$，偏自相关系数，统一见 3.5. 变形推导过程

1.3.2. MA(q)模型

MA(q)模型，表示时间$t$下的$x_t$由过去的$t-1,t-2,...,t-q$时间下的历史扰动项${\epsilon }_{t-1},{\epsilon }_{t-2},...,{\epsilon }_{t-q}$决定，称之为q阶的移动平均（Moving Average）过程：

$$x_t=\mu +{\epsilon }_t-{\theta }_1{\epsilon }_{t-1}-{\theta }_2{\epsilon }_{t-2}-\cdots -{\theta }_q{\epsilon }_{t-q}(2)$$

数字特征： $E({\epsilon }_t)=0,D({\epsilon }_t)={\sigma }_{{\epsilon }_t}^2$

结合以上的数字特征，藉由式（2）变形可得知：

1. $E(x_t)=\mu$
2. $D(x_t)=(1+{\theta }_1^2+{\theta }_2^2+\cdots +{\theta }_q^2)\cdot {\sigma }_{\epsilon }^2$ 与q相关

变形过程推导，自相关系数${\rho }_k$，协方差$\mathrm{r}(\mathrm{k})=\mathrm{c}\mathrm{o}\mathrm{v}(x_t,x_{t-k})$，偏自相关系数，统一见 3.5. 变形推导过程

1.3.3. ARMA(p,q)模型

将纯AR(p)与纯MA(q)结合，得到一个一般的自回归移动平均（Auto Regressive Moving Average）过程ARMA(p,q)：
$$x_t={\varphi }_0+{\varphi }_1{x}_{t-1}+{\varphi }_2{x}_{t-2}+\cdots +{\varphi }_p{x}_{t-p}+{\epsilon }_t-{\theta }_1{\epsilon }_{t-1}-{\theta }_2{\epsilon }_{t-2}-\cdots -{\theta }_q{\epsilon }_{t-q}(3)$$

把前半部分看出随着时间变化的$T_t$，后面半部分看出历史扰动项$I_t$，可以得出数字特征：

1. $E(x_t)=\frac{{\varphi }_0}{1-{\varphi }_1-\cdots -{\varphi }_p}$
2. $D(x_t)=D(T_t)+D(I_t)$ 同时与p,q相关

变形过程推导，自相关系数${\rho }_k$，协方差$\mathrm{r}(\mathrm{k})=\mathrm{c}\mathrm{o}\mathrm{v}(x_t,x_{t-k})$，偏自相关系数，模型变量太多了，这里我就懒得算了，谁考试让你用手推这个那完全就是变态了

1.3.4. ARIMA(p,d,q)模型

ARIMA模型（Auto Regressive Integrated Moving Average model），差分整合移动平均自回归模型，又称整合移动平均自回归模型（移动也可称作滑动）。ARIMA(p，d，q)中，AR是“自回归”，p为自回归项数；MA为“滑动平均”，q为滑动平均项数，d为使之成为平稳序列所做的差分次数（阶数）

这里需要证明一下，当时间序列本身不是平稳的时候，如果它的增量，即的一次差分，稳定在零点附近，可以将看成是平稳序列。所以在实际的问题中，所遇到的多数非平稳序列可以通过一次或多次差分后成为平稳时间序列，再藉由差分后的平稳序列带入之前的ARMA模型，这就是ARIMA的建模构成方法：

1.3.5. 变形推导过程

二、基于Matlab进行计算

Step1：指定非季节性 ARIMA 模型

建立一个 ARIMA(2,1,2) 模型：
$$(1-{\varphi }_{1}L-{\varphi }_2{L}^2)(1-L)y_t=(1+{\theta }_{1}L+{\theta }_2{L}^2){\epsilon }_t$$

Mdl = arima(2,1,2)

Step2：参数评估，并查看生成模型的

[EstMdl,EstParamCov,logL,info] = estimate(Mdl,y,Name,Value)

其中生成info中由参数详情，EstParamCov有结果分析

Step3：预测结果

[Y,YMSE] = forecast(Mdl,numperiods,Y0,Name,Value)

numperiods为预测范围，或预测周期内的时间点数，指定为正整数

Y0用于初始化预测模型的预采样响应数据，指定为长度为 numpreobs 的数值列向量或 numpreobs-by-numpaths 数值矩阵

生成Y代表预测结果，YMSE代表准确性MSE得分