平稳时间序列预测法PPT课件

.,1,7平稳时间序列预测法,7.1概述7.2时间序列的自相关分析7.3单位根检验和协整检验7.4ARMA模型的建模,.,2,7.1概述,时间序列取自某一个随机过程，则称：,一、平稳时间序列,过程是平稳的随机过程的随机特征不随时间变化而变化,过程是非平稳的随机过程的随机特征随时间变化而变化,.,3,宽平稳时间序列的定义：,设时间序列,，对于任意的t,k和m，满足：,则称宽平稳。,严平稳时间序列的定义：所有的统计特性不随时间的平移而变化,.,4,Box-Jenkins基本思想：用数学模型描述时间序列自身的相关性，并假定这种自相关性一直延续，用该模型预测未来的值。,ARMA模型是描述平稳随机序列的最常用的一种模型。,Box-Jenkins方法提供了对时间序列进行分析、预测，以及对ARMA模型识别、估计和诊断的系统方法。,.,5,ARMA模型的三种基本形式：自回归模型（AR：Auto-regressive）；移动平均模型（MA：Moving-Average）；混合模型（ARMA：Auto-regressiveMoving-Average）。,.,6,如果时间序列满足其中是独立同分布的随机变量序列,且满足：则称时间序列服从p阶自回归模型。,二、自回归模型,.,7,滞后算子多项式,的根均在单位圆外，即,的根大于1。,自回归模型的平稳条件：,.,8,例1AR(1)模型的平稳性条件。,对1阶自回归模型AR(1),方程两边平方再求数学期望，得到Xt的方差,由于Xt仅与t相关，因此，E(Xt-1t)=0。如果该模型稳定，则有E(Xt2)=E(Xt-12)，从而上式可变换为：,在稳定条件下，该方差是一非负的常数，从而有|1。,.,9,而AR(1)的特征方程,的根为z=1/AR(1)稳定，即|1，意味着特征根大于1。,对高阶自回模型AR(p)来说，多数情况下没有必要直接计算其特征方程的特征根，但有一些有用的规则可用来检验高阶自回归模型的稳定性：,(1)AR(p)模型稳定的必要条件是：1+2+p3时都落入置信区间，且逐渐趋于零，则该时间序列具有平稳性；若时间序列的自相关函数更多地落在置信区间外面，则该时间序列不具有平稳性。,.,27,注：在B-J方法中，只有平稳的时间序列才能建立ARMA模型，否则必须经过适当的处理使序列满足平稳性要求。例对某种趋势的时间序列进行差分处理。但很多序列不能通过差分达到平稳，而且差分虽然消除了序列的趋势易于建模，但也消除了序列的长期特征，实际的经济序列差分一般不超过两次。,.,28,3、时间序列的季节性判定准则：,月度数据，考察k=12,24,36,时的自相关系数是否与0有显著差异；季度数据，考察k=4,8,12,时的自相关系数是否与0有显著差异。,注1：实际问题中常遇到季节性和趋势性同时存在的情况，应先剔除序列趋势性，在识别季节性。注2：包含季节性的时间序列也不能直接建模，应先进行季节差分消除，季节差分一般不超过一阶。,.,29,三、ARMA模型的自相关分析,AR(p)模型的偏自相关函数是以p步截尾的，自相关函数拖尾；MA(q)模型的自相关函数具有q步截尾性，偏自相关函数拖尾；（可用以上两个性质来识别AR和MA模型的阶数）ARMA(p,q)模型的自相关函数和偏相关函数都是拖尾的。,.,30,.,31,.,32,.,33,7.4ARMA模型的建模,一、模型阶数的确定,（1）基于自相关函数和偏相关函数的定阶方法,对于ARMA(p,q)模型，可以利用其样本的自相关函数和样本偏自相关函数的截尾性判定模型的阶数。,如果样本的偏自相关函数是以p步截尾的，模型为AR(p)；如果样本的自相关函数具有q步截尾性，模型为MA(q)；如果样本的自相关函数和偏相关函数都是拖尾的，模型为ARMA(p,q)。,.,34,（1）自相关函数的截尾性统计检验：,对于每一个q,计算,.,（M取,为或者），,对于MA(q)模型，当kq时，,考察其中满足的个数是否占M个的68.3%或者95.5%以上。,.,35,（2）偏自相关函数的截尾性统计检验：,对于每一个p,计算,（M取,为或者），,对于AR(p)模型，当kp时，,考察其中满足的个数是否占M个的68.3%或者95.5%以上。,.,36,如果对于序列,和,截尾，即不存在上述的,来说，均不,和,判定平稳时间序列,，则可以,为ARMA模型。,.,37,一般地，对ARMA,模型,它们均值为0，可递推得到残量估计,现作假设检验：,是来自白噪声的样本,令,（3）残差项的白噪声检验：（Q统计量检验）,.,38,其中,取,左右。,则当,成立时，,服从,的,分布。,对给定的显著性水平,，若,，则拒绝,，即模型与原随机序列之间拟合得不好，,，则认为模型与原随机序列之间拟合,需重新考虑,得较好，模型检验被通过。,建模；若,自由度为,注：上机操作时，一般看Q统计检验的相伴概率,.,39,(1)用AR（1）拟合时间序列，考察其残差样本的自相关函数是否q1步截尾，则模型为ARMA(1,q1),否则；(2)用AR（2）拟合时间序列，考察其残差样本的自相关函数是否q2步截尾，则模型为ARMA(1,q2)，否则；(3)继续增大p，重复上述做法，直至残差序列的样本自相关函数截尾为止,（4）Tasy和TiaoARMA模型定阶法,.,40,1950年-1998年北京城乡居民定期储蓄比例,选择合适的ARMA模型拟合可以考虑拟合模型为AR(1),ARMA(1,3),.,41,连续读取70个化学反应数据,可以尝试使用AR(1)，MA(1)和ARMA(1,1)模型拟合该序列,.,42,（2）基于F检验确定阶数（3）利用信息准则法定阶（AIC准则和BIC准则）,此外，常用的方法还有：,.,43,1967年，瑞典控制论专家K.J.Astrm教授将F检验准则用于对时间序列模型的定阶。原理(模型阶数简约原则parsimonyprinciple)：设yt(1tn)是零均值平稳序列，用模型AR模型拟合检验统计量：结论若FF，则拒绝原假设，认为AR(p)合适；若FF，则拒绝原假设，模型阶数仍有上升的可能；若F1时，ARMA(1,1)预测值也是由如下差分方程决定的。,(3)向前L步预测公式(L2),.,60,三、预测误差,.,61,由于预测只能建立在到t时刻为止的可用信息的基础上，因此，根据最小均方误预测的第二个准则，以及平稳可逆序列可以表示成传递函数形式的论断，可以将预测值表示成能够估计的项t，t-1，的加权和的形式：,.,62,由上得以t为原点，向前L步的预测误差为：,由于t是白噪声，故有：,.,63,误差方差为：,.,64,注：预测误差的估算是1，p算和1，q估计都为正确的假设，实际中参数通过估计得到的，且估计量是随机变量，有均值和方差，因而实际误差大于理论估计误差。,.,65,五、预测误差的置信区间,对于正态过程，预测误差的分布为：,所以：对yt+l预测的95%的置信区间为：,因此：,.,66,设,为一AR(2)序列，,其中,求,的自协方差函数,。,例1,解答：,Yule-Walker方程为：,且：,联合上面三个方程，解出：,.,67,例2,考虑如下AR(2)序列：,若已知观测值,（1）试预报,（2）给出（1）预报的置信度为95%的预报区间。,.,68,解答：,（1）,（2）,预报的置信度为95%的预报区间分别为：,.,69,7.3单位根检验和协整检验,一、平稳性的检验,引言：前面我们讨论的是平稳时间序列的建模和预测方法，即所讨论的时间序列都是宽平稳的。一个宽平稳的时间序列的均值和方差都是常数，并且它的协方差有时间上的不变性。但是许多经济领域产生的时间序列都是非平稳的。呈现出明显得趋势性和周期性，序列不平稳，导致预测无效，产生谬误回归等问题。,.,70,1、通过时间序列的趋势图来判断,这种方法通过观察时间序列的趋势图来判断时间序列是否存在趋势性或周期性。优点：简便、直观。对于那些明显为非平稳的时间序列，可以采用这种方法。缺点：对于一般的时间序列是否平稳，不易用这种方法判断出来。,.,71,2、通过自相关函数(ACF)判断,平稳时间序列的自相关函数(ACF)要么是截尾的，要么是拖尾的。因此我们可以根据这个特性来判断时间序列是否为平稳序列。若时间序列具有上升或下降的趋势，那么对于所有短时滞来说，自相关系数大且为正，而且随着时滞k的增加而缓慢地下降。,.,72,若序列无趋势，但是具有季节性，那末对于按月采集的数据，时滞12，24，36的自相关系数达到最大(如果数据是按季度采集，则最大自相关系数出现在4，8，12，)，并且随着时滞的增加变得较小。若序列是有趋势的，且具有季节性，其自相关函数特性类似于有趋势序列，但它们是摆动的，对于按月数据，在时滞12，24，36，等处具有峰态；如果时间序列数据是按季节的，则峰出现在时滞4，8，12，等处。,.,73,3、随机游走的单位根检验(Unitroottest),随机游走是一种非平稳过程，其实随机游走一种特殊的齐次非平稳过程。检验序列是否为随机游走，通常利用DavidDickey和WayneFuller的单位根检验。单位根的含义和检验原理如下：,.,74,（1）单位根的含义,.,75,（2）单位根的检验,用Eviews进行单位根检验时给出了上述选项。如果DF检验统计量比给定显著水平临界值大，不能拒绝原假设，认为序列存在单位根，是非平稳的。,.,76,ADF检验在DF检验中，常常因为序列存在高阶滞后相关，使得随机扰动不符合白噪声假设，ADF检验修正了DF检验中的自相关问题。,此外还有：PP检验：检验具有一般形式的单位根过程DFGLS检验：DF及ADF检验对含有时间趋势的退势平稳时间序列的检验失效,古典的回归方法：只能对平稳的时间序列进行回归分析，或将非平稳的序列转化为平稳序列再做回归,.,77,非平稳时间序列分析,逐期差分后平稳，建立求和自回归移动平均模型,记为ARIMA(p,d,q),.,78,非平稳时间序列分析,逐期差分+季节差分后平稳，建立乘积季节模型，即随机季节模型与