第六讲-时间序列分析

第一章时间序列分析简介第二章时间序列的预处理第三章平稳时间序列分析第四章非平稳序列的确定性分析第五章非平稳序列的随机分析《应用时间序列分析》第一章时间序列分析简介1.1引言按照时间的顺序把随机事件变化发展的过程记录下来就构成了一个时间序列。对时间序列进行观察、研究，找寻它变化发展的规律，预测它将来的走势就是时间序列分析。1.2时间序列的定义随机序列:按时间顺序排列的一组随机变量观察值序列:随机序列的个有序观察值，称之为序列长度为的观察值序列,,,,,21tXXXtxxx,,,21tt1.3时间序列分析方法描述性时序分析统计时序分析描述性时序分析通过直观的数据比较或绘图观测，寻找序列中蕴含的发展规律，这种分析方法就称为描述性时序分析描述性时序分析方法具有操作简单、直观有效的特点，它通常是人们进行统计时序分析的第一步。描述性时序分析案例德国业余天文学家施瓦尔发现太阳黑子的活动具有11年左右的周期统计时序分析频域分析方法时域分析方法频域分析方法原理假设任何一种无趋势的时间序列都可以分解成若干不同频率的周期波动特点非常有用的动态数据分析方法，但是由于分析方法复杂，结果抽象，有一定的使用局限性时域分析方法原理事件的发展通常都具有一定的惯性，这种惯性用统计的语言来描述就是序列值之间存在着一定的相关关系，这种相关关系通常具有某种统计规律。目的寻找出序列值之间相关关系的统计规律，并拟合出适当的数学模型来描述这种规律，进而利用这个拟合模型预测序列未来的走势特点操作步骤规范，分析结果易于解释，是时间序列分析的主流方法时域分析方法的分析步骤考察观察值序列的特征根据序列的特征选择适当的拟合模型根据序列的观察数据确定模型的口径检验模型，优化模型利用拟合好的模型来推断序列其它的统计性质或预测序列将来的发展1.4时间序列分析软件常用软件SPSS，Matlab，Eviews，SAS，S-plus/R和TSP，第二章时间序列的预处理2.1平稳性检验特征统计量平稳时间序列的定义平稳时间序列的统计性质平稳时间序列的意义平稳性的检验概率分布概率分布的意义随机变量族的统计特性完全由它们的联合分布函数或联合密度函数决定时间序列概率分布族的定义TtttmmxxxFmmtttm,,,),,,2,1()},,,({2121,,,21特征统计量均值方差自协方差自相关系数)(xxdFEXttt)()()(22xdFxXEDXttttt))((),(ssttXXEststDXDXstst),(),(平稳时间序列的定义严平稳严平稳是一种条件比较苛刻的平稳性定义，它认为只有当序列所有的统计性质都不会随着时间的推移而发生变化时，该序列才能被认为平稳。宽平稳宽平稳是使用序列的特征统计量来定义的一种平稳性。它认为序列的统计性质主要由它的低阶矩决定，所以只要保证序列低阶矩平稳（二阶），就能保证序列的主要性质近似稳定。平稳时间序列的统计定义满足如下条件的序列称为严平稳序列满足如下条件的序列称为宽平稳序列),,,(),,,(21,21,2121mtttmtttxxxFxxxFmm有，正整数，正整数Ttttmm,,,,21TtskksttskkstTtEXTtEXtt且，为常数，,,),(),()3,)2,)12严平稳与宽平稳的关系一般关系严平稳条件比宽平稳条件苛刻，通常情况下，严平稳（低阶矩存在）能推出宽平稳成立，而宽平稳序列不能反推严平稳成立特例不存在低阶矩的严平稳序列不满足宽平稳条件，例如服从柯西分布的严平稳序列就不是宽平稳序列当序列服从多元正态分布时，宽平稳可以推出严平稳平稳时间序列的统计性质常数均值自协方差函数和自相关函数只依赖于时间的平移长度而与时间的起止点无关延迟k自协方差函数延迟k自相关系数)0()(kk为整数kkttk),,()(平稳时间序列的意义时间序列数据结构的特殊性可列多个随机变量，而每个变量只有一个样本观察值平稳性的重大意义极大地减少了随机变量的个数，并增加了待估变量的样本容量极大地简化了时序分析的难度，同时也提高了对特征统计量的估计精度平稳性的检验（图检验方法）时序图检验根据平稳时间序列均值、方差为常数的性质，平稳序列的时序图应该显示出该序列始终在一个常数值附近随机波动，而且波动的范围有界、无明显趋势及周期特征自相关图检验平稳序列通常具有短期相关性。该性质用自相关系数来描述就是随着延迟期数的增加，平稳序列的自相关系数会很快地衰减向零例题例2.1检验1964年——1999年中国纱年产量序列的平稳性例2.1时序图例2.1自相关图例2.2检验1962年1月——1975年12月平均每头奶牛月产奶量序列的平稳性例2.2时序图例2.2自相关图例2.3检验1949年——1998年北京市每年最高气温序列的平稳性例2.3时序图例2.3自相关图2.2纯随机性检验纯随机序列的定义纯随机性的性质纯随机性检验纯随机序列的定义纯随机序列也称为白噪声序列，它满足如下两条性质TststststTtEXt,,,0,),()2(,)1(2标准正态白噪声序列时序图白噪声序列的性质纯随机性各序列值之间没有任何相关关系，即为“没有记忆”的序列方差齐性根据马尔可夫定理，只有方差齐性假定成立时，用最小二乘法得到的未知参数估计值才是准确的、有效的00k(k)，)0(2tDX纯随机性检验检验原理假设条件检验统计量判别原则Barlett定理如果一个时间序列是纯随机的，得到一个观察期数为的观察序列，那么该序列的延迟非零期的样本自相关系数将近似服从均值为零，方差为序列观察期数倒数的正态分布0,)1,0(~ˆknNkn假设条件原假设：延迟期数小于或等于期的序列值之间相互独立备择假设：延迟期数小于或等于期的序列值之间有相关性1,0210mHm：mkmHk，：至少存在某个1,01mm检验统计量Q统计量LB统计量)(~ˆ212mnQmkk)(~)ˆ()2(212mknnnLBmkk判别原则拒绝原假设当统计量的P值小于时，则可以以的置信水平拒绝原假设，认为该序列为非白噪声序列接受原假设当该统计量的P值大于时，则认为在的置信水平下无法拒绝原假设，即不能显著拒绝序列为纯随机序列的假定11例2.4：标准正态白噪声序列纯随机性检验样本自相关图检验结果LBQLBQ延迟统计量检验统计量值P值延迟6期2.360.8838延迟12期5.350.9454由于P值显著大于显著性水平，所以该序列不能拒绝纯随机的原假设。例2.5对1950年——1998年北京市城乡居民定期储蓄所占比例序列的平稳性与纯随机性进行检验例2.5时序图例2.5自相关图例2.5白噪声检验结果延迟阶数LB统计量检验LB检验统计量的值P值675.460.00011282.570.0001第三章平稳时间序列分析3.1方法性工具差分运算延迟算子线性差分方程差分运算一阶差分阶差分步差分pk1tttxxx111tptptpxxxkttkxx延迟算子延迟算子类似于一个时间指针，当前序列值乘以一个延迟算子，就相当于把当前序列值的时间向过去拨了一个时刻记B为延迟算子，有1,pxBxtppt延迟算子的性质，其中10B为任意常数cxcxBcxcBttt,)()(111)(ttttByBxyxBnttnxxBiniinnnBCB0)1()1()!(!!ininCin用延迟算子表示差分运算阶差分步差分pkitpiipptptpxCxBx0)1()1(tkkttkxBxx)1(线性差分方程线性差分方程齐次线性差分方程)(2211thzazazazptpttt02211ptptttzazazaz齐次线性差分方程的解特征方程特征方程的根称为特征根，记作齐次线性差分方程的通解不相等实数根场合有相等实根场合复根场合02211ppppaaap,,,21tpptttcccz2211tpptddtddtcctctccz111121)(tpptititttccececrz3321)(非齐次线性差分方程的解非齐次线性差分方程的特解使得非齐次线性差分方程成立的任意一个解非齐次线性差分方程的通解齐次线性差分方程的通解和非齐次线性差分方程的特解之和tttzzztz)(2211thzazazazptpttttz3.2ARMA模型的性质AR模型（AutoRegressionModel）MA模型（MovingAverageModel）ARMA模型（AutoRegressionMovingAveragemodel）AR模型的定义具有如下结构的模型称为阶自回归模型，简记为特别当时，称为中心化模型tsExtsEVarExxxxtsstttptptpttt,0,0)(,)(0)(0222110，p)(pAR00)(pARAR(P)序列中心化变换称为的中心化序列，令p101ttxy}{ty}{tx自回归系数多项式引进延迟算子，中心化模型又可以简记为自回归系数多项式)(pARttxB)(ppBBBB2211)(自回归模型平稳的充分必要条件是：的根在单位圆外；或者的根在单位圆内2312310pp1231230ppppp例3.1:考察如下四个模型的平稳性1(1)0.8tttxx1(2)1.1tttxx12(3)0.5ttttxxxttttxxx215.0)4(例3.1平稳序列时序图1(1)0.8tttxx12(3)0.5ttttxxx例3.1非平稳序列时序图1(2)1.1tttxxttttxxx215.0)4(例3.1平稳性判别8.011.11211i212i23112312模型特征根判别结论(1)平稳(2)非平稳(3)平稳(4)非平稳平稳AR模型的统计性质均值方差协方差自相关系数偏自相关系数均值如果AR(p)模型满足平稳性条件，则有根据平稳序列均值为常数，且为白噪声序列，有推导出p101)(110tptpttxxEExTtEExtt,0)(,}{tGreen函数定义AR模型的传递形式其中系数称为Green函数},2,1,{jGjjtjjjpijtjiipijtjiipitiittGkBkBkBx001101)(1)(ppBBBB2211)(Green函数递推公式原理方法待定系数法递推公式pkpkjGGGkkkjjkkj,0,,,2,1110其中，ttttttBGBBGxxB)()()()(方差平稳AR模型的传递形式两边求方差得函数为GreenGGxVarjjjt,)(202jtjjtGx0例3.2:求平稳AR(1)模型的方差平稳AR(