第三章平稳时间序列分析 (2)优秀课件.ppt

《第三章平稳时间序列分析 (2)优秀课件.ppt》由会员分享，可在线阅读，更多相关《第三章平稳时间序列分析 (2)优秀课件.ppt（153页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第三章平稳时间序列分析第1页，本讲稿共153页本章结构n方法性工具 nARMA模型 n平稳序列建模n序列预测 第2页，本讲稿共153页3.1 方法性工具 n差分运算n延迟算子n线性差分方程第3页，本讲稿共153页差分运算n一阶差分n 阶差分 n 步差分第4页，本讲稿共153页延迟算子n延迟算子类似于一个时间指针，当前序列值乘以一个延迟算子，就相当于把当前序列值的时间向过去拨了一个时刻 n记B为延迟算子，有 第5页，本讲稿共153页延迟算子的性质n n n n n ，其中 第6页，本讲稿共153页用延迟算子表示差分运算n 阶差分 n 步差分第7页，本讲稿共153页线性差分方程 n线性差分方程n齐

2、次线性差分方程第8页，本讲稿共153页齐次线性差分方程的解n特征方程n特征方程的根称为特征根，记作n齐次线性差分方程的通解n不相等实数根场合n有相等实根场合n复根场合第9页，本讲稿共153页非齐次线性差分方程的解 n非齐次线性差分方程的特解n使得非齐次线性差分方程成立的任意一个解n非齐次线性差分方程的通解n齐次线性差分方程的通解和非齐次线性差分方程的特解之和第10页，本讲稿共153页3.2 ARMA模型的性质 nAR模型（Auto Regression Model）nMA模型（Moving Average Model）nARMA模型（Auto Regression Moving Average

3、 model）第11页，本讲稿共153页AR模型的定义n具有如下结构的模型称为 阶自回归模型，简记为n特别当 时，称为中心化 模型第12页，本讲稿共153页 AR(P)序列中心化变换n称 为 的中心化序列，令第13页，本讲稿共153页自回归系数多项式n引进延迟算子，中心化 模型又可以简记为 n自回归系数多项式第14页，本讲稿共153页AR模型平稳性判别 n判别原因nAR模型是常用的平稳序列的拟合模型之一，但并非所有的AR模型都是平稳的 n判别方法n单位根判别法n平稳域判别法第15页，本讲稿共153页例3.1:考察如下四个模型的平稳性第16页，本讲稿共153页例3.1平稳序列时序图第17页，本讲

4、稿共153页例3.1非平稳序列时序图第18页，本讲稿共153页AR模型平稳性判别方法n特征根判别nAR(p)模型平稳的充要条件是它的p个特征根都在单位圆内n根据特征根和自回归系数多项式的根成倒数的性质，等价判别条件是该模型的自回归系数多项式的根都在单位圆外n平稳域判别 n平稳域第19页，本讲稿共153页AR(1)模型平稳条件n特征根n平稳域第20页，本讲稿共153页AR(2)模型平稳条件n特征根n平稳域第21页，本讲稿共153页例3.1平稳性判别模型特征根判别平稳域判别结论(1)平稳(2)非平稳(3)平稳(4)非平稳第22页，本讲稿共153页平稳AR模型的统计性质n均值n方差n协方差n自相关系

5、数n偏自相关系数第23页，本讲稿共153页均值 n如果AR(p)模型满足平稳性条件，则有n根据平稳序列均值为常数，且 为白噪声序列，有n推导出第24页，本讲稿共153页Green函数定义nAR模型的传递形式n其中系数 称为Green函数第25页，本讲稿共153页Green函数递推公式n原理n方法n待定系数法n递推公式第26页，本讲稿共153页方差n平稳AR模型的传递形式n两边求方差得第27页，本讲稿共153页例3.2:求平稳AR(1)模型的方差n平稳AR(1)模型的传递形式为nGreen函数为n平稳AR(1)模型的方差第28页，本讲稿共153页协方差函数n在平稳AR(p)模型两边同乘 ，再求期

6、望n根据n得协方差函数的递推公式第29页，本讲稿共153页例3.3:求平稳AR(1)模型的协方差n递推公式n平稳AR(1)模型的方差为n协方差函数的递推公式为第30页，本讲稿共153页例3.4:求平稳AR(2)模型的协方差n平稳AR(2)模型的协方差函数递推公式为第31页，本讲稿共153页自相关系数n自相关系数的定义n平稳AR(P)模型的自相关系数递推公式第32页，本讲稿共153页常用AR模型自相关系数递推公式nAR(1)模型nAR(2)模型第33页，本讲稿共153页AR模型自相关系数的性质n拖尾性n呈复指数衰减第34页，本讲稿共153页例3.5:考察如下AR模型的自相关图第35页，本讲稿共1

7、53页例3.5n自相关系数按复指数单调收敛到零第36页，本讲稿共153页例3.5:第37页，本讲稿共153页例3.5:n自相关系数呈现出“伪周期”性第38页，本讲稿共153页例3.5:n自相关系数不规则衰减第39页，本讲稿共153页偏自相关系数n定义对于平稳AR(p)序列，所谓滞后k偏自相关系数就是指在给定中间k-1个随机变量 的条件下，或者说，在剔除了中间k-1个随机变量的干扰之后，对 影响的相关度量。用数学语言描述就是第40页，本讲稿共153页偏自相关系数的计算n滞后k偏自相关系数实际上就等于k阶自回归模型第个k回归系数的值。第41页，本讲稿共153页偏自相关系数的截尾性nAR(p)模型偏

8、自相关系数P阶截尾第42页，本讲稿共153页例3.5续:考察如下AR模型的偏自相关图第43页，本讲稿共153页例3.5n理论偏自相关系数n样本偏自相关图第44页，本讲稿共153页例3.5:n理论偏自相关系数n样本偏自相关图第45页，本讲稿共153页例3.5:n理论偏自相关系数n样本偏自相关图第46页，本讲稿共153页例3.5:n理论偏自相关系数n样本偏自相关系数图第47页，本讲稿共153页MA模型的定义n具有如下结构的模型称为 阶自回归模型，简记为n特别当 时，称为中心化 模型第48页，本讲稿共153页移动平均系数多项式n引进延迟算子，中心化 模型又可以简记为 n 阶移动平均系数多项式第49页

9、，本讲稿共153页MA模型的统计性质n常数均值n常数方差第50页，本讲稿共153页MA模型的统计性质n自协方差函数P阶截尾n自相关系数P阶截尾第51页，本讲稿共153页常用MA模型的自相关系数nMA(1)模型nMA(2)模型第52页，本讲稿共153页MA模型的统计性质n偏自相关系数拖尾第53页，本讲稿共153页例3.6:考察如下MA模型的相关性质第54页，本讲稿共153页MA模型的自相关系数截尾n n 第55页，本讲稿共153页MA模型的自相关系数截尾n n 第56页，本讲稿共153页MA模型的偏自相关系数拖尾n n 第57页，本讲稿共153页MA模型的偏自相关系数拖尾n n 第58页，本讲稿

10、共153页MA模型的可逆性nMA模型自相关系数的不唯一性n例3.6中不同的MA模型具有完全相同的自相关系数和偏自相关系数第59页，本讲稿共153页可逆的定义n可逆MA模型定义n若一个MA模型能够表示称为收敛的AR模型形式，那么该MA模型称为可逆MA模型n可逆概念的重要性n一个自相关系数列唯一对应一个可逆MA模型。第60页，本讲稿共153页可逆MA(1)模型n n 第61页，本讲稿共153页MA模型的可逆条件nMA(q)模型的可逆条件是：nMA(q)模型的特征根都在单位圆内n等价条件是移动平滑系数多项式的根都在单位圆外第62页，本讲稿共153页逆函数的递推公式n原理n方法n待定系数法n递推公式第

11、63页，本讲稿共153页例3.6续:考察如下MA模型的可逆性第64页，本讲稿共153页(1)(2)n n n逆函数n逆转形式第65页，本讲稿共153页(3)(4)n n n逆函数n逆转形式第66页，本讲稿共153页ARMA模型的定义n具有如下结构的模型称为自回归移动平均模型，简记为n特别当 时，称为中心化 模型第67页，本讲稿共153页系数多项式n引进延迟算子，中心化 模型又可以简记为 n 阶自回归系数多项式n 阶移动平均系数多项式第68页，本讲稿共153页平稳条件与可逆条件nARMA(p,q)模型的平稳条件nP阶自回归系数多项式 的根都在单位圆外n即ARMA(p,q)模型的平稳性完全由其自回

12、归部分的平稳性决定nARMA(p,q)模型的可逆条件nq阶移动平均系数多项式 的根都在单位圆外n即ARMA(p,q)模型的可逆性完全由其移动平滑部分的可逆性决定第69页，本讲稿共153页传递形式与逆转形式n传递形式n逆转形式第70页，本讲稿共153页ARMA(p,q)模型的统计性质n均值n协方差n自相关系数第71页，本讲稿共153页ARMA模型的相关性n自相关系数拖尾n偏自相关系数拖尾第72页，本讲稿共153页例3.7:考察ARMA模型的相关性n拟合模型ARMA(1,1):并直观地考察该模型自相关系数和偏自相关系数的性质。第73页，本讲稿共153页自相关系数和偏自相关系数拖尾性n样本自相关图n

13、样本偏自相关图第74页，本讲稿共153页ARMA模型相关性特征模型自相关系数偏自相关系数AR(P)拖尾P阶截尾MA(q)q阶截尾拖尾ARMA(p,q)拖尾拖尾第75页，本讲稿共153页3.3平稳序列建模 n建模步骤n模型识别n参数估计n模型检验n模型优化n序列预测第76页，本讲稿共153页建模步骤平稳非白噪声序列计算样本相关系数模型识别参数估计模型检验模型优化序列预测YN第77页，本讲稿共153页计算样本相关系数n样本自相关系数n样本偏自相关系数第78页，本讲稿共153页模型识别n基本原则选择模型拖尾P阶截尾AR(P)q阶截尾拖尾MA(q)拖尾拖尾ARMA(p,q)第79页，本讲稿共153页模

14、型定阶的困难n因为由于样本的随机性，样本的相关系数不会呈现出理论截尾的完美情况，本应截尾的 或 仍会呈现出小值振荡的情况n由于平稳时间序列通常都具有短期相关性，随着延迟阶数 ，与 都会衰减至零值附近作小值波动？当 或 在延迟若干阶之后衰减为小值波动时，什么情况下该看作为相关系数截尾，什么情况下该看作为相关系数在延迟若干阶之后正常衰减到零值附近作拖尾波动呢？第80页，本讲稿共153页样本相关系数的近似分布nBarlettnQuenouille第81页，本讲稿共153页模型定阶经验方法n95的置信区间n模型定阶的经验方法n如果样本(偏)自相关系数在最初的d阶明显大于两倍标准差范围，而后几乎95的自

15、相关系数都落在2倍标准差的范围以内，而且通常由非零自相关系数衰减为小值波动的过程非常突然。这时，通常视为(偏)自相关系数截尾。截尾阶数为d。第82页，本讲稿共153页例2.5续n选择合适的模型ARMA拟合1950年1998年北京市城乡居民定期储蓄比例序列。第83页，本讲稿共153页序列自相关图第84页，本讲稿共153页序列偏自相关图第85页，本讲稿共153页拟合模型识别n自相关图显示延迟3阶之后，自相关系数全部衰减到2倍标准差范围内波动，这表明序列明显地短期相关。但序列由显著非零的相关系数衰减为小值波动的过程相当连续，相当缓慢，该自相关系数可视为不截尾 n偏自相关图显示除了延迟1阶的偏自相关系

16、数显著大于2倍标准差之外，其它的偏自相关系数都在2倍标准差范围内作小值随机波动，而且由非零相关系数衰减为小值波动的过程非常突然，所以该偏自相关系数可视为一阶截尾 n所以可以考虑拟合模型为AR(1)第86页，本讲稿共153页例3.8美国科罗拉多州某一加油站连续57天的OVERSHORT序列 第87页，本讲稿共153页序列自相关图第88页，本讲稿共153页序列偏自相关图第89页，本讲稿共153页拟合模型识别n自相关图显示除了延迟1阶的自相关系数在2倍标准差范围之外，其它阶数的自相关系数都在2倍标准差范围内波动。根据这个特点可以判断该序列具有短期相关性，进一步确定序列平稳。同时，可以认为该序列自相关

17、系数1阶截尾n偏自相关系数显示出典型非截尾的性质。n综合该序列自相关系数和偏自相关系数的性质，为拟合模型定阶为MA(1)第90页，本讲稿共153页例3.9n1880-1985全球气表平均温度改变值差分序列 第91页，本讲稿共153页序列自相关图第92页，本讲稿共153页序列偏自相关图第93页，本讲稿共153页拟合模型识别n自相关系数显示出不截尾的性质n偏自相关系数也显示出不截尾的性质n综合该序列自相关系数和偏自相关系数的性质，可以尝试使用ARMA(1,1)模型拟合该序列第94页，本讲稿共153页参数估计n待估参数n 个未知参数n常用估计方法n矩估计n极大似然估计n最小二乘估计第95页，本讲稿共

18、153页矩估计n原理n样本自相关系数估计总体自相关系数n样本一阶均值估计总体均值，样本方差估计总体方差第96页，本讲稿共153页例3.10:求AR(2)模型系数的矩估计nAR(2)模型nYule-Walker方程n矩估计（Yule-Walker方程的解）第97页，本讲稿共153页例3.11:求MA(1)模型系数的矩估计nMA(1)模型n方程n矩估计第98页，本讲稿共153页例3.12:求ARMA(1,1)模型系数的矩估计nARMA(1,1)模型n方程n矩估计第99页，本讲稿共153页对矩估计的评价n优点n估计思想简单直观n不需要假设总体分布n计算量小（低阶模型场合）n缺点n信息浪费严重n只用到

19、了p+q个样本自相关系数信息，其他信息都被忽略n估计精度差n通常矩估计方法被用作极大似然估计和最小二乘估计迭代计算的初始值 第100页，本讲稿共153页极大似然估计n原理n在极大似然准则下，认为样本来自使该样本出现概率最大的总体。因此未知参数的极大似然估计就是使得似然函数（即联合密度函数）达到最大的参数值 第101页，本讲稿共153页似然方程n由于 和 都不是 的显式表达式。因而似然方程组实际上是由p+q+1个超越方程构成，通常需要经过复杂的迭代算法才能求出未知参数的极大似然估计值 第102页，本讲稿共153页对极大似然估计的评价n优点n极大似然估计充分应用了每一个观察值所提供的信息，因而它的

20、估计精度高n同时还具有估计的一致性、渐近正态性和渐近有效性等许多优良的统计性质n缺点n需要假定总体分布第103页，本讲稿共153页最小二乘估计n原理n使残差平方和达到最小的那组参数值即为最小二乘估计值 第104页，本讲稿共153页条件最小二乘估计n实际中最常用的参数估计方法n假设条件n残差平方和方程n解法n迭代法第105页，本讲稿共153页对最小二乘估计的评价n优点n最小二乘估计充分应用了每一个观察值所提供的信息，因而它的估计精度高n条件最小二乘估计方法使用率最高n缺点n需要假定总体分布第106页，本讲稿共153页例2.5续n确定1950年1998年北京市城乡居民定期储蓄比例序列拟合模型的口径

21、 n拟合模型：AR(1)n估计方法：极大似然估计n模型口径第107页，本讲稿共153页例3.8续n确定美国科罗拉多州某一加油站连续57天的OVERSHORTS序列拟合模型的口径 n拟合模型：MA(1)n估计方法：条件最小二乘估计n模型口径第108页，本讲稿共153页例3.9续n确定1880-1985全球气表平均温度改变值差分序列拟合模型的口径 n拟合模型：ARMA(1,1)n估计方法：条件最小二乘估计n模型口径第109页，本讲稿共153页模型检验n模型的显著性检验n整个模型对信息的提取是否充分n参数的显著性检验n模型结构是否最简第110页，本讲稿共153页模型的显著性检验n目的n检验模型的有效

22、性（对信息的提取是否充分）n检验对象n残差序列n判定原则n一个好的拟合模型应该能够提取观察值序列中几乎所有的样本相关信息，即残差序列应该为白噪声序列 n反之，如果残差序列为非白噪声序列，那就意味着残差序列中还残留着相关信息未被提取，这就说明拟合模型不够有效第111页，本讲稿共153页假设条件n原假设：残差序列为白噪声序列n备择假设：残差序列为非白噪声序列第112页，本讲稿共153页检验统计量nLB统计量第113页，本讲稿共153页例2.5续n检验1950年1998年北京市城乡居民定期储蓄比例序列拟合模型的显著性 n残差白噪声序列检验结果延迟阶数LB统计量P值检验结论65.830.3229拟合模

23、型显著有效1210.280.50501811.380.8361第114页，本讲稿共153页参数显著性检验n目的n检验每一个未知参数是否显著非零。删除不显著参数使模型结构最精简 n假设条件n检验统计量第115页，本讲稿共153页例2.5续n检验1950年1998年北京市城乡居民定期储蓄比例序列极大似然估计模型的参数是否显著 n参数检验结果检验参数t统计量P值结论均值46.120.0001显著6.720.0001显著第116页，本讲稿共153页例3.8续:对OVERSHORTS序列的拟合模型进行检验 n残差白噪声检验n参数显著性检验检验参数t统计量P值结论均值3.750.0004显著10.600.

24、0001显著延迟阶数LB统计量P值结论63.150.6772模型显著有效129.050.6171第117页，本讲稿共153页例3.9续:对1880-1985全球气表平均温度改变值差分序列拟合模型进行检验 n残差白噪声检验n参数显著性检验检验参数t统计量P值结论16.340.0001显著3.50.0007显著延迟阶数LB统计量P值结论65.280.2595模型显著有效1210.300.4247第118页，本讲稿共153页模型优化n问题提出n当一个拟合模型通过了检验，说明在一定的置信水平下，该模型能有效地拟合观察值序列的波动，但这种有效模型并不是唯一的。n优化的目的n选择相对最优模型 第119页，

25、本讲稿共153页例3.13:拟合某一化学序列第120页，本讲稿共153页序列自相关图第121页，本讲稿共153页序列偏自相关图第122页，本讲稿共153页拟合模型一n根据自相关系数2阶截尾，拟合MA(2)模型n参数估计n模型检验n模型显著有效 n三参数均显著 第123页，本讲稿共153页拟合模型二n根据偏自相关系数1阶截尾，拟合MA(1)模型n参数估计n模型检验n模型显著有效 n两参数均显著 第124页，本讲稿共153页问题n同一个序列可以构造两个拟合模型，两个模型都显著有效，那么到底该选择哪个模型用于统计推断呢？n解决办法n确定适当的比较准则，构造适当的统计量，确定相对最优第125页，本讲稿

26、共153页AIC准则n最小信息量准则（An Information Criterion）n指导思想n似然函数值越大越好 n未知参数的个数越少越好 nAIC统计量第126页，本讲稿共153页SBC准则nAIC准则的缺陷n在样本容量趋于无穷大时，由AIC准则选择的模型不收敛于真实模型，它通常比真实模型所含的未知参数个数要多 nSBC统计量第127页，本讲稿共153页例3.13续n用AIC准则和SBC准则评判例3.13中两个拟合模型的相对优劣 n结果nAR(1)优于MA(2)模型AICSBCMA(2)536.4556543.2011AR(1)535.7896540.2866第128页，本讲稿共153

27、页序列预测n线性预测函数n预测方差最小原则第129页，本讲稿共153页序列分解预测误差预测值第130页，本讲稿共153页误差分析n估计误差n期望n方差第131页，本讲稿共153页AR(p)序列的预测n预测值n预测方差n95置信区间第132页，本讲稿共153页例3.14n已知某超市月销售额近似服从AR(2)模型（单位：万元/每月）n今年第一季度该超市月销售额分别为：101，96，97.2万元n请确定该超市第二季度每月销售额的95的置信区间 第133页，本讲稿共153页例3.14解：预测值计算n四月份n五月份n六月份第134页，本讲稿共153页例3.14解：预测方差的计算nGREEN函数n方差第1

28、35页，本讲稿共153页例3.14解：置信区间n公式n估计结果预测时期95置信区间四月份（85.36，108.88）五月份（83.72，111.15）六月份（81.84，113.35）第136页，本讲稿共153页例2.5:北京市城乡居民定期储蓄比例序列拟合与预测图 第137页，本讲稿共153页MA(q)序列的预测n预测值n预测方差第138页，本讲稿共153页例3.15n已知某地区每年常驻人口数量近似服从MA(3)模型（单位：万）：n最近3年的常驻人口数量及一步预测数量如下：n预测未来5年该地区常住人口的95置信区间年份统计人数预测人数200210411020031081002004105109

29、第139页，本讲稿共153页例3.15解：随机扰动项的计算第140页，本讲稿共153页例3.15解：估计值的计算第141页，本讲稿共153页例3.15解：预测方差的计算第142页，本讲稿共153页例3.15解：置信区间的计算预测年份95置信区间2005（99，119）2006（83，109）2007（87，115）2008（86，114）2009（86，114）第143页，本讲稿共153页ARMA(p,q)序列预测n预测值n预测方差第144页，本讲稿共153页例3.16n已知模型为：n且 n预测未来3期序列值的95的置信区间。第145页，本讲稿共153页例3.16解：估计值的计算第146页，本

30、讲稿共153页例3.16解：预测方差的计算nGreen函数n方差第147页，本讲稿共153页例3.16解：置信区间的计算时期95置信区间101（0.136，0.332）102（0.087，0.287）103（0.049，0.251）第148页，本讲稿共153页修正预测n定义n所谓的修正预测就是研究如何利用新的信息去获得精度更高的预测值 n方法n在新的信息量比较大时把新信息加入到旧的信息中，重新拟合模型 n在新的信息量很小时不重新拟合模型，只是将新的信息加入以修正预测值，提高预测精度第149页，本讲稿共153页修正预测原理n在旧信息的基础上，的预测值为n假设新获得一个观察值 ，则n 的修正预测值为n修正预测误差为n预测方差为第150页，本讲稿共153页一般情况n假设新获得p个观察值 ，则n 的修正预测值为n修正预测误差为n预测方差为第151页，本讲稿共153页例3.14续:假如四月份的真实销售额为100万元，求二季度后两个月销售额的修正预测值 n计算四月份的预测误差n计算修正预测值n计算修正方差第152页，本讲稿共153页修正置信区间预测时期修正前置信区间修正后置信区间四月份（85.36，108.88）五月份（83.72，111.15）（87.40，110.92）六月份（81.84，113.35）（85.79，113.21）第153页，本讲稿共153页