第五章 时间序列模型.docVIP

5 时间序列模型 5.1 时间序列 数字化技术的应用和发展使得随机序列的分析变得日益广泛和重要，并由平稳随机过程在时间轴上的 取样引出平稳离散随机信号或时间序列的概念。对于这类随机序列，主要采用相关函数和功率谱进行分析。对于平稳离散时间信号，还常用时间序列描述方法进行研究，由此提出时间序列模型法。它是采用各种随机差分方程表示时间序列信号的模型。在许多情况下，一个平稳离散随机信号可以视为白噪声序列通过某一离散时间线性系统所产生的。 在时间序列信号模型分析中，AR（自回归）模型、MA（滑动平均）模型和ARMA（自回归滑动平均）模型是三种最常见的标准线性模型，它们均由白噪声序列通过离散时间线性系统而产生。而实际应用中许多平稳时间序列往往可由这些模型近似表示，使得有关的分析变得更为简单，也为平稳随机序列的分析和产生提供了有效方法。另外，这些线性模型都具有连续功率谱形状，在参数谱估计方面显示出极大的优点。除非特别说明，本章只讨论具有连续谱特性的平稳时间序列。 5.2 自回归(AR)模型 设为具有零均值，方差为的平稳白噪声序列，随机序列由如下随机差分方程表示： 式中为一正整数，为实常数，不失一般性，设，并设。上式表示的信号称为阶自回归模型。显然，是它的个过去值和白噪声的线性组合。用表示上式的模型。对于上式，从统计观点讲，称以随机误差线性回归于它的个过去值。 为使分析方便，首先研究一阶和二阶AR模型，然后根据阶AR模型的分析，研究AR模型的自相关函数及功率谱密度。 一阶AR模型 根据随机序列的差分表达式，当时，可得一阶AR模型 式中为不等于零的实常数。上式为一阶随机差分方程。若设，可得： 容易得到一阶矩 如果，由上式可以看出，的均值有可能不满足平稳性，即可能不满足一阶平稳。然而，如果系数，当较大时，则有 在此情况下，是一阶渐进平稳的。 通常，，可得时间序列的自相关函数（二阶矩）为： 显然，当时，并不满足自相关平稳性，但是，当并且足够大时，有 对于实随机序列，由于对于对称分布，有 对于，不难推得，当为正数时，恒为正，且呈指数衰减。当为负数时，正负相间指数衰减。 根据可得的方差为： 说明平稳随机序列的方差比白噪声方差大。 最后讨论AR(1)模型的功率谱。对式两边取z变换，可得其传递函数为： 的功率谱为 令，有 二阶AR模型 定义随机序列的二阶AR模型为： 式中和均为实常数，。 上式二阶差分方程的特征多项式为： 定义后移算子为后移一步的运算，即 于是，二阶AR模型成为： 式中和为二阶AR模型特征多项式的根，即 所以，有特解为： 根据模型差分方程，零输入下得齐次方程 其解为： 式中和是待定系数，由初始条件确定。 模型特解和上式之和即为模型的解： 当时，上式右边齐次解随的增大而趋于零，而特解部分具有有限方差，在均方意义下收敛，随的增大而渐近收敛于特解公式的平稳结果。 实际上，二阶模型的平稳条件与其系数和是有关的，这可通过和平面表示。设，并设和，根据，在其两边同乘，有 或 其次，根据不等式，两边同乘，有 或 根据上式分析，得到以下三个条件： 以及 这就是保证二阶AR模型平稳的条件，可用系数分布图说明。图中示出了二阶系数欠阻尼、过阻尼和临界阻尼三种情况的系数区域分布，分别对应于以下三种情况： 欠阻尼：出现和一对共轭复根。 过阻尼：出现和不同的实根。 临界阻尼：出现和相同的实根。 对于平稳的情况，考察二阶AR模型的自相关函数，对模型方差方程两边同乘并作集平均，可得： 考虑到 可得： 以及 由此解得： 最后，分析AR(2)模型的功率谱密度。容易知道，其传递函数为： 于是，的功率谱为： 3. 阶AR模型 定义如下随机差分方程为阶AR模型 式中为实常数，且。 对上式两边取变换，可得： 于是，以上AP(p)模型的传递函数为： 根据它的特征多项式可解出个的极点。于是，该模型的传递函数可写为： 所以，AR模型的传递函数只有极点，除原点外没有任何零点，属于全极点模型，对应于全极点滤波器，具有无限冲激响应（IIR）。因此，模型传递函数的性质完全取决于个极点在平面上的分布情况。可以证明，如果所有个极点均满足，那么，AR模型信号满足渐近平稳性。条件意味着有界输入通过线性系统导致有界输出，系统是稳定的，这说明模型传递函数的稳定性与模型的平稳性是等价的。 根据AR模型的传递函数，阶AR模型的功率谱密度为： 可见AR模型的功率谱由各模型系数确定。 最后讨论AR(p)模型参数与相关函数的关系。根据自相关函数的定义，有 由于 于是，有： 将