《时间序列-模型.ppt》由会员分享,可在线阅读,更多相关《时间序列-模型.ppt(65页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、时间序列时间序列-模型模型1现在学习的是第1页,共65页四、四、自回归模型的阶自回归模型的阶数的估计数的估计五、五、自回归模型自回归模型的参数的估计的参数的估计六、六、自回归自回归模型的检模型的检验验七、七、自回归自回归模型的预报模型的预报现在学习的是第2页,共65页 时间序列分析最重要的应用是分析和表征观察值时间序列分析最重要的应用是分析和表征观察值之间的相互依赖性与相关性,若对这种相关性进行量之间的相互依赖性与相关性,若对这种相关性进行量化处理,那么就可以方便地从系统的过去值预测将来化处理,那么就可以方便地从系统的过去值预测将来的值。的值。在数理统计中讨论的数据的线性回归模型在数理统计中讨
2、论的数据的线性回归模型,很好地表示了因变量很好地表示了因变量yt的观察值对自变量观测值的观察值对自变量观测值xt1,xt2,xtp的相关性,解决了他们之间的相关的相关性,解决了他们之间的相关性问题,但是,对一组随机观测数据,即一个性问题,但是,对一组随机观测数据,即一个时间序列内部的相关关系它却描述不出来。即时间序列内部的相关关系它却描述不出来。即它不能描述数据内部之间的相互依赖关系。它不能描述数据内部之间的相互依赖关系。现在学习的是第3页,共65页 另一方面,某些随机过程与另一些变量另一方面,某些随机过程与另一些变量取值之间的随机关系,往往根本无法用任何取值之间的随机关系,往往根本无法用任何
3、函数关系式来描述,这时就需要采用这个时函数关系式来描述,这时就需要采用这个时间序列本身的观测数据之间的依赖关系来揭间序列本身的观测数据之间的依赖关系来揭示这个时序的规律性。示这个时序的规律性。现在学习的是第4页,共65页一、自回归模型的定义一、自回归模型的定义定义定义2.1 设设xt,t=0,1,2,为时间序列,白为时间序列,白噪声序列为噪声序列为t,t=0,1,2,,且对任意的,且对任意的 s1,所以这是平稳的所以这是平稳的AR(1)模型。模型。现在学习的是第8页,共65页例例1.3 如果时间序列如果时间序列xt满足满足试问此试问此xt是否为平稳的序列模型。是否为平稳的序列模型。解:由于其自
4、回归系数多项式为解:由于其自回归系数多项式为的根为的根为u1=21与与u2=1/21,故知其根不都故知其根不都在单位圆外,所以这是非平稳的在单位圆外,所以这是非平稳的AR(2)序列模序列模型。型。现在学习的是第9页,共65页 自回归模型是描述系统内部的回归关系,故自回归模型是描述系统内部的回归关系,故称为自回归,与通常的线性回归性质是不一样称为自回归,与通常的线性回归性质是不一样的。的。现在学习的是第10页,共65页二、中心化二、中心化 AR(p)模型模型设设xt为平稳序列,且有为平稳序列,且有则对上式两端同取数学期望,即得则对上式两端同取数学期望,即得由于由于xt为平稳序列,故为平稳序列,故
5、 现在学习的是第11页,共65页即得即得则可得一个均值为则可得一个均值为0的新序列:的新序列:此时此时wt 称为称为xt 的的平稳中心化序列平稳中心化序列。现在学习的是第12页,共65页以后一般均讨论中心化的平稳模型以后一般均讨论中心化的平稳模型或序列或序列:现在学习的是第13页,共65页三、平稳模型的平稳解三、平稳模型的平稳解设平稳设平稳AR(p)模型为模型为式中式中t为白噪声序列,为白噪声序列,系数系数1,2,p 满足平稳条件:系数满足平稳条件:系数 多项式多项式(u)=0的根都在单位圆外。的根都在单位圆外。现在学习的是第14页,共65页1 后移算子后移算子若算子若算子B满足等式满足等式:
6、则称则称B为后移算子为后移算子,即即B作用作用xt 后使其转化为后使其转化为xt-1类似的类似的 于是于是,AR(p)模型可以表示为模型可以表示为现在学习的是第15页,共65页即得一差分方程即得一差分方程:其中其中(B)(B)为后移算子多项式为后移算子多项式,即称为自回归即称为自回归算子算子:现在学习的是第16页,共65页 易见,滤波器成为一个对时间序列进行变换的易见,滤波器成为一个对时间序列进行变换的实体,变换前的序列称为输入,经滤波器变换的得实体,变换前的序列称为输入,经滤波器变换的得到的序列称为输出。到的序列称为输出。自回归滤波器自回归滤波器xtt t 差分方程式可用框图表示差分方程式可
7、用框图表示:设想有一个滤波设想有一个滤波器,输入的是某种平稳序列,而输出的则是白器,输入的是某种平稳序列,而输出的则是白噪声序列,即噪声序列,即现在学习的是第17页,共65页2 AR(p)序列的平稳域与允许域序列的平稳域与允许域 定义定义2.2 AR(p)序列的平稳域为其系数取值的序列的平稳域为其系数取值的集合:集合:其允许域为其自相关函数的前其允许域为其自相关函数的前p个值的集合个值的集合:其中矩阵其中矩阵p p与与Rp,向量向量、b与与d分别为:分别为:现在学习的是第18页,共65页现在学习的是第19页,共65页例如一阶自回归模型例如一阶自回归模型AR(1):现在学习的是第20页,共65页
8、注注:实际上由平稳实际上由平稳AR(p)模型模型:现在学习的是第21页,共65页再对两端取数学期望再对两端取数学期望,并由性质并由性质:现在学习的是第22页,共65页类似的类似的,在平稳在平稳AR(p)模型两端分别同乘以模型两端分别同乘以现在学习的是第23页,共65页再对两端取数学期望再对两端取数学期望,并由上述性质可得并由上述性质可得:现在学习的是第24页,共65页其次其次,由于自相关系数等于由于自相关系数等于:现在学习的是第25页,共65页现在学习的是第26页,共65页3 AR(1)序列平稳解与自相关函数序列平稳解与自相关函数进行反复的迭代运算,则对任何自然数进行反复的迭代运算,则对任何自
9、然数n,有,有 现在学习的是第27页,共65页于是对于平稳时间序列,如果有于是对于平稳时间序列,如果有|0时时,其自相关函数:其自相关函数:现在学习的是第30页,共65页现在学习的是第31页,共65页类似的类似的,当当k0时时,其自相关函数为其自相关函数为现在学习的是第32页,共65页特别地,特别地,AR(1)序列的方差函数为序列的方差函数为其自相关系数为其自相关系数为 因为因为|10时,其自相关函数为时,其自相关函数为现在学习的是第39页,共65页综述为综述为 用上式求得用上式求得AR(2)序列的自相关函数是比较困难的,序列的自相关函数是比较困难的,在实际中常采用另一种简便有效的方法在实际中
10、常采用另一种简便有效的方法:现在学习的是第40页,共65页 设设h0,因,因t与与xt互不相关,故用互不相关,故用xt-h乘以乘以AR(2)模模型等式两端,再取数学期望,即得型等式两端,再取数学期望,即得上述递推式称为上述递推式称为AR(2)序列的序列的Yule-Walker方程。方程。利用利用Yule-Walker方程求方程求AR(2)序列的自相关序列的自相关函数方法即称为尤尔函数方法即称为尤尔-沃克法沃克法.现在学习的是第41页,共65页注注:由于此处均值函数为:由于此处均值函数为0,其自相关函数与自,其自相关函数与自协方差函数相等,为分析简单起见,我们将自相协方差函数相等,为分析简单起见
11、,我们将自相关函数作为自协方差函数,即表示为关函数作为自协方差函数,即表示为而将自相关系数而将自相关系数:称为自相关函数。称为自相关函数。因此,由因此,由Yule-Walker方程可得方程可得现在学习的是第42页,共65页其初值为其初值为 下面分三种情况讨论在上述给定初始条件下面分三种情况讨论在上述给定初始条件下下,自相关函数的差分方程的解自相关函数的差分方程的解:现在学习的是第43页,共65页再根据初值条件再根据初值条件,自回归系数方程的根与系数自回归系数方程的根与系数的关系得的关系得即得即得现在学习的是第44页,共65页 由上式可看出,当自回归系数方程的根的两实根都由上式可看出,当自回归系
12、数方程的根的两实根都在单位圆外部时,在单位圆外部时,(k)随随k的增大,向零衰减,的增大,向零衰减,若两实根中至少有一个在单位圆内部时,则若两实根中至少有一个在单位圆内部时,则(k)发散。发散。现在学习的是第45页,共65页AR(2)序列的平稳域序列的平稳域:平稳域图形见下图:平稳域图形见下图:221-20 1-1-11现在学习的是第46页,共65页再根据初值条件再根据初值条件,自回归系数方程的根与系数的关自回归系数方程的根与系数的关系得系得即得现在学习的是第47页,共65页(3)AR(2)序列的允许域序列的允许域现在学习的是第48页,共65页允许域允许域图形如下图所示图形如下图所示 21-1
13、0 1现在学习的是第49页,共65页5.一般一般AR(p)序列的平稳解与自相关函数序列的平稳解与自相关函数(1)AR(p)模型的平稳解模型的平稳解定理定理2.1.1 设设AR(p)序列的系数多项式序列的系数多项式 的所有根均在单位圆外部,即满足平稳条件,且的所有根均在单位圆外部,即满足平稳条件,且p p00。若存在实数列。若存在实数列j,j=0,1,2,满足满足p阶齐阶齐次线性差分方程:次线性差分方程:现在学习的是第50页,共65页及初值条件及初值条件:现在学习的是第51页,共65页则均方极限存在则均方极限存在 且几乎必然对一切且几乎必然对一切 t=0,1,2,,此式为,此式为AR(p)模型的
14、平稳解。模型的平稳解。(2)AR(p)模型的自相关函数模型的自相关函数定理定理2.1.2 设设xt 是是AR(p)序列,其自回归算子的序列,其自回归算子的所有根都在单位圆外部,则其自协方差函数所有根都在单位圆外部,则其自协方差函数C(k)满满足足p阶齐次差分方程阶齐次差分方程Yule-Walker方程:方程:现在学习的是第52页,共65页及初值条件及初值条件现在学习的是第53页,共65页 若在上式两端同除以若在上式两端同除以C(0)=Dx,即得自相关函即得自相关函数满足的数满足的Yule-Walker方程,即方程,即及初始条件及初始条件现在学习的是第54页,共65页 时间序列的自相关函数刻划了
15、随机序时间序列的自相关函数刻划了随机序列各时刻间的线性相关程度。实际中常用列各时刻间的线性相关程度。实际中常用自相关函数的图形自相关函数的图形相关图来分析、反映相关图来分析、反映时间序列各时刻间的线性相关性。时间序列各时刻间的线性相关性。现在学习的是第55页,共65页 6.AR(p)序列偏相关函数序列偏相关函数 对于任意平稳序列,若其自协方差函数对于任意平稳序列,若其自协方差函数C(k)满满足以下条件,即对任意足以下条件,即对任意k1,1i,jk,k k为正定为正定矩阵矩阵:现在学习的是第56页,共65页再将向量的分量记为再将向量的分量记为称为上述平稳序列的称为上述平稳序列的偏相关函数偏相关函
16、数。现在学习的是第57页,共65页 可以利用以下递推公式计算可以利用以下递推公式计算偏相关函数偏相关函数:一般来说,偏相关函数的意义由下述定理给一般来说,偏相关函数的意义由下述定理给出出:现在学习的是第58页,共65页定理定理2.1.3 对于零均值的平稳序列对于零均值的平稳序列xt而言,以下两而言,以下两条是相互等价的;条是相互等价的;(1)xt满足平稳序列满足平稳序列AR(p)模型;模型;(2)xt的偏相关函数列满足条件的偏相关函数列满足条件 此定理表明偏相关函数列的截尾性质是平此定理表明偏相关函数列的截尾性质是平稳自回归序列独有的特征。利用此定理,可以稳自回归序列独有的特征。利用此定理,可
17、以通过检查通过检查xt的偏相关函数列的的偏相关函数列的p 阶截尾性质来识阶截尾性质来识别别AR(p)模型模型。现在学习的是第59页,共65页四、四、自回归模型的阶自回归模型的阶数的估计数的估计1、利用偏相关函数的截尾性质确定阶数、利用偏相关函数的截尾性质确定阶数 因为因为AR(p)的偏相关函数在的偏相关函数在p后截尾,所以利后截尾,所以利用样本自相关函数与尤尔用样本自相关函数与尤尔-沃克方程递推求得偏沃克方程递推求得偏相关函数的值,或用软件求出:相关函数的值,或用软件求出:现在学习的是第60页,共65页递推求得偏相关函数的值:递推求得偏相关函数的值:现在学习的是第61页,共65页 从序列的自相
18、关与偏样本偏相关函数值列表观从序列的自相关与偏样本偏相关函数值列表观察察:如果样本偏相关函数值在如果样本偏相关函数值在p以后截尾以后截尾;而样本偏相关函数值明显拖尾而样本偏相关函数值明显拖尾,则可认定序列为则可认定序列为 AR(p)模型模型.现在学习的是第62页,共65页例如:序列的自相关与偏样本偏相关函数值如下例如:序列的自相关与偏样本偏相关函数值如下:k1234567k0.80.590.420.320.250.170.10kk0.8-0.15 00.08-0.03-0.06-0.02k891011121314k0.050.030.030-0.05-0.07-0.08kk0.0200.04-0.02-0.09-0.04 0.01现在学习的是第63页,共65页样本偏相关函数值折线图样本偏相关函数值折线图可认为样本偏相关函数值在可认为样本偏相关函数值在p=2以后截尾以后截尾,即在即在p=2以后样本偏相关函数值波动幅度很小以后样本偏相关函数值波动幅度很小.现在学习的是第64页,共65页样本自相关函数值折线图样本自相关函数值折线图:可认为样本自相关函数值单调下降可认为样本自相关函数值单调下降,有拖尾趋势有拖尾趋势,故故此模型可取为此模型可取为AR(2).现在学习的是第65页,共65页