ARMA模型的时域特性.ppt

第三章 ARMA模型的特性,ARMA模型，一方面，它基于观测时间序列 建立起来的随机微分方程，因而它解释了动态数据的统计特性；另一方面，由于 可视为某一系统的输出，因而，它又揭示了产生此动态数据的系统的动态特性。 同时，不论是数据的统计特性，还是系统的动态特性，均可在时域和频域中得到描述，所有这些特性，构成了ARMA模型的基本特性。,本章重点讨论ARMA模型的最主要的时域特性系统的单位脉冲响应函数 和动态数据的自协方差函数 。前者表征系统特性，在时序方法中又称为Green函数，后者表征数据的统计特性。 同时，还将介绍ARMA模型的另外两个时域特性逆函数和偏自相关函数。,讨论模型特性的目的在于，一方面，它是实际应用的理论基础，很多实际问题的解决往往就是模型特性直接应用的结果；另一方面，它又是建立模型的必要准备。,线性常系数差分方程及其解的一般形式,在时间序列的时域分析中，线性差分方程是非常重要，也是极为有效的工具。 任何一个ARMA模型都是一个线性差分方程；因此，ARMA模型的性质往往取决于差分方程根的性质。 为了更好地讨论ARMA模型的特性，先简单介绍线性差分方程的一般知识。,时间序列模型与线性差分方程,线性差分方程在时间序列分析中有着重要的应用，常用的时间序列模型和某些模型的自协方差函数和自相关函数都可以视为线性差分方程，而线性差分方程对应的特征根的性质对判断模型的平稳性有着非常重要的意义。,是普通的n阶差分方程，其中 为系统参数的函数，当 为常数时，就是常系数n阶差分方程， 是个离散序列，也叫驱动函数； 是系统的响应。当 时，上式变为 称为n阶齐次差分方程。,线性差分方程,线性差分方程 齐次线性差分方程,设,AR(1)模型的Green函数,1、AR(1)模型的Green函数,首先，将最简单的AR(1)模型作为一个例子。 AR(1)模型： 反复进行迭代,即：,Green函数的定义,当一个相关的平稳时间序列可以用一个无关的平稳时间序列的现在值和过去值的线性组合表示时，其“权”定义为Green函数，即 式中， 称为Green函数， ，,（1）式可以记为,其中,式（1）表明具有传递形式的平稳序列可以由现在时刻以前的白噪声通过系统“ ”的作用而生成， 是j个单位时间以前加入系统的干扰项 对现实响应的权，亦即系统对 的“记忆”。,格林函数的意义,格林函数的含义：格林函数是描述系统记忆扰动程度的函数。,Green函数刻画了系统动态响应衰减的快慢程度。 Green函数所描述的动态性完全取决于系统参数 。,则AR(1)模型的格林函数可以表示为： AR(1)模型可表示为 同时，可用一个无限阶MA来逼近。,例：下面是参数分别为0.9、0.1的AR（1）系统对,扰动的记忆情况。（P46）,AR(1)系统的平稳性,系统稳定性的概念以及稳定性与平稳性的关系,一阶系统的稳定性,Green函数的另一个重要作用是， 可表明系统的稳定性这一重要的动态特性。所谓一个系统是不稳定的，是指它在任意瞬间受到一个一瞬即逝的干扰(即脉冲)后，其运动状态偏离平衡位置越来越远，这相当于 ， 是发散的；反之，如果其运动状态最终能回到平衡位置上，这相当于 ，则称系统是渐进稳定的；,线性系统的稳定性仅由系统本身的固有特性所决定，而与外界无关，即，ARMA模型所描述的线性系统，其稳定性只与AR部分有关，而与MA部分无关，因此，AR(1)，ARMA(1,1)，ARMA(1,m)系统的稳定性问题实质上是一致的，从而可根据Green函数的取值情况判断它们所对应的不同的一阶系统的稳定性。,2、 AR(1)系统的平稳性条件,平稳性的涵义就是干扰项对系统的影响逐渐减弱，直到消失，对于一个AR（1）系统，将其写成格林函数的表示形式:,如果系统是平稳的，则预示随着j，扰动的权数,对于AR(1)系统,即,这要求,上述条件等价于AR(1)系统的特征方程,的根在单位圆内(或方程,的根在单位圆外).,AR(1)的结论可以推广到AR(n),ARMA(2,1)模型的Green函数,AR(2)和ARMA(1,1)模型的Green函数,AR(2)和ARMA(1,1)模型是ARMA(2,1)模型的特殊形式； 描述动态性的Green函数也有上述关系；,ARMA(1,1)模型的Green函数,ARMA(2,1)系统的平稳性,1、用特征根表示的平稳性条件 这个推论在AR(1)中平稳性的条件，同样对ARMA(2,1)模型也依然适应；此时， ARMA(2,1)系统的平稳性条件为： 即，特征方程的特征根的模在单位圆内,ARMA(n,n-1)系统的平稳性,2、用自回归系数表示的平稳性条件,AR(n)模型的Green函数,AR(n)模型Green函数的递推公式为：,第二节 逆函数和可逆性 （Invertibility）,是零均值平稳序列，如果白噪声序列,能够表示为,一、逆函数的定义,设,则称上式为平稳序列,式中的加权系数,称为逆函数。,的”逆转形式“。,1、逆函数 类似Green函数，逆函数定义为：当一个无关的平稳时间序列 可以用一个相关的平稳时间序列 的现在值和过去值的线性组合来表示时，其负“权”定义为逆函数.,可逆的定义,可逆定义 若一个模型能够表示成为收敛的AR模型形式，那么该模型具有可逆性，也就是可逆的。 可逆概念的重要性 一个自相关系数列唯一对应一个可逆MA模型。,AR(1)模型的逆函数,逆函数,Green函数,MA(1)模型的逆函数,逆函数,Green函数,格林函数与逆函数间关系,格林函数与逆函数间的这种对偶性不只是一阶模型所有，对于任意阶模型都成立。 例如：ARMA(2,1)与ARMA(1,2),MA(m)模型逆函数的递推公式,如果一个MA(m)模型满足可逆性条件，它就可以写成如下两种等价形式：,MA(m)模型逆函数的递推公式,MA模型的可逆条件,MA(m)模型的可逆条件是： MA(m)模型的特征根都在单位圆内,ARMA(1,2)模型的可逆性条件,例3.6续:考察如下MA模型的可逆性,(1)(2),逆函数 逆转形式,(3)(4),逆函数 逆转形式,ARMA模型,一、ARMA(n,m)模型可分别表示为：,其中：,平稳条件与可逆条件,ARMA(n,m)模型的平稳条件 n阶自回归系数多项式 的根都在单位圆外 即ARMA(n,m)模型的平稳性完全由其自回归部分的平稳性决定 ARMA(n,m)模型的可逆条件 m阶移动平均系数多项式 的根都在单位圆外 即ARMA(n,m)模型的可逆性完全由其移动平滑部分的可逆性决定,理论自协方差函数和自相关函数 对于ARMA系统来说，设序列的均值为零，则自协方差函数,第三节 自相关函数与偏自相关函数,自相关函数,样本自相关函数的计算 在拟合模型之前，我们所有的只是序列的一个有限样本数据，无法求得理论自相关函数，只能求样本的自协方差函数和自相关函数。样本自协方差有两种形式：,一、自相关函数,则相应的样本自相关函数为：,1、AR(n)过程自相关函数ACF,1阶自回归模型AR(1) Xt=Xt-1+ at 的k阶滞后自协方差函数为：,0,1,1,),(,(,g,j,jg,a,j,g,k,k,t,t,k,t,k,X,X,E,=,=,+,=,-,-,-,k=1,2,因此，AR(1)模型的自相关函数为,k=1,2,若AR(1) 稳定，则| 1，因此，k时，呈指数形衰减，直到零。这种现象称为拖尾或称AR(1)有无穷记忆（infinite memory）。 注意， 0时，呈振荡衰减状。,一般地，n阶自回归模型AR(n) Xt=1Xt-1+ 2Xt-2 + nXt-n + at,k期滞后协方差为:,n,k,n,k,k,t,n,t,n,t,t,K,t,k,X,X,X,X,E,-,-,-,-,-,-,-,+,+,+,=,+,+,+,+,=,g,j,g,j,g,j,a,j,j,j,g,L,L,2,2,1,1,2,2,1,1,),(,(,从而有自相关函数 :,可见，无论k有多大， k 的计算均与其到n阶滞后的自相关函数有关，因此呈拖尾状。 如果AR(n)是平稳的，则|k|递减且趋于零。,其中：zi是AR(n)特征方程(z)=0的特征根，由AR(n)平稳的条件知，|zi|1; 因此，当zi均为实数根时，k呈几何型衰减（单调或振荡）； 当存在虚数根时，则一对共扼复根构成通解中的一个阻尼正弦波项， k呈正弦波衰减。,事实上，自相关函数,是一n阶差分方程，其通解为,对MA(1)过程,2、MA(m)过程,1,-,-,=,t,t,t,X,qa,a,可容易地写出它的自协方差函数：,0,),1,(,3,2,2,1,2,2,0,=,=,=,-,=,+,=,L,g,g,qs,g,s,q,g,a,a,于是，MA(1)过程的自相关函数为：,可见，当k1时，k=0，即Xt与Xt-k不相关，MA(1)自相关函数是截尾的。,其自协方差系数为,一般地，m阶移动平均过程MA(m),相应的自相关函数为,可见，当km时， Xt与Xt-k不相关，即存在截尾现象，因此，当km时， k=0是MA(m)的一个特征。 于是：可以根据自相关系数是否从某一点开始一直为0来判断MA(m)模型的阶。,二、偏自相关函数,自相关函数ACF(k)给出了Xt与Xt-1的总体相关性，但总体 相关性可能掩盖了变量间完全不同的隐含关系。 例如，在AR(1)随机过程中，Xt与Xt-2间有相关性可能主要是由于它们各自与Xt-1间的相关性带来的:,即自相关函数中包含了这种所有的“间接”相关。 与之相反，Xt与Xt-k间的偏自相关函数(partial autocorrelation，简记为PACF)则是消除了中间变量Xt-1，Xt-k+1 带来的间接相关后的直接相关性，它是在已知序列值Xt-1，Xt-k+1的条件下，Xt与Xt-k间关系的度量。,从Xt中去掉Xt-1的影响，则只剩下随机扰动项at，显然它与Xt-2无关，因此我们说Xt与Xt-2的偏自相关函数 为零，记为,在AR(1)中，,0,),(,2,*,2,=,=,-,t,t,X,Corr,a,r,对于AR(1) 过程，当k = 1时, 1 0，当k 1时, k* =0，所以AR(1) 过程的偏自相关函数特征是在k = 1出现峰值（ 1 = 1*）然后截尾。,AR(n)模型,自相关函数： 平滑地指数衰减,偏自相关函数： k=1时有正峰值然后截尾,AR(1)模型相关函数与偏自相关函数对比,同样地，在AR(n)过程中，对所有的kn，Xt与Xt-k间的偏自相关函数为零。 AR(n)的一个主要特征是: kn时，k*=Corr(Xt,Xt-k)=0 即k*在n以后是截尾的。,一随机时间序列的识别原则： 若Xt的偏自相关函数在n以后截尾，即kn时，k*=0，而它的自相关函数k是拖尾的，则此序列是自回归AR(n)序列。,对于一个k阶AR模型，有：,由此得到Yule-Walker 方程，记为：,已知时，由该方程组可以解出,利用克莱姆法则，对k=1，2，3， 依次求解得,上述,序列为AR模型的偏自相关函数。,偏自相关性是条件相关，是在给定,的条件下，,和,的条件相关。换名话说，偏自相关,函数是对,和,所解释的相关的度量。,之间未被,易得，在AR(k)模型中，第k个偏自相关系数就是AR(k)模型中Xt-k的回归系数。,如果自回归过程的阶数为n，则对于kn应该有 =0。,MA(m)模型 可以验证MA(m)过程的偏自相关函数是非截尾但趋于零的。,MA(m)模型的识别规则：若随机序列的自相关函数截尾，即自m以后，k=0（ km）；而它的偏自相关函数是拖尾的