第三章平稳时间序列分析本章结构方法性工具ARMA模型平稳序列建模序列预测方法性工具ARMA模型平稳序列建模序列预测3.1方法性工具差分运算延迟算子线性差分方程方法性工具ARMA模型平稳序列建模序列预测差分运算一阶差分阶差分步差分pk1tttxxx111tptptpxxxkttkxx例如:原数列12510172637一阶差分1357911二阶差分22222方法性工具ARMA模型平稳序列建模序列预测延迟算子延迟算子类似于一个时间指针,当前序列值乘以一个延迟算子,就相当于把当前序列值的时间向过去拨了一个时刻记B为延迟算子,有1,pxBxtppt方法性工具ARMA模型平稳序列建模序列预测延迟算子的性质,其中10B为任意常数cxcxBcxcBttt,)()(111)(ttttyxyxBnttnxxBiniinnnBCB0)1()1()!(!!ininCin方法性工具ARMA模型平稳序列建模序列预测用延迟算子表示差分例如,三阶差分3221112121122312301233313233()()=(2)=2()()=33=CCCCttttttttttttttttttttttttxxxxxxxxxxxxxxxxxxxxxxxx333300(1)(1)iiiiititiiCxCBx3(1)tBx方法性工具ARMA模型平稳序列建模序列预测用延迟算子表示差分运算阶差分步差分pkitpiipptptpxCxBx0)1()1(tkkttkxBxx)1(方法性工具ARMA模型平稳序列建模序列预测线性差分方程线性差分方程齐次线性差分方程)(2211thzazazazptpttt02211ptptttzazazaz齐次线性差分方程的解特征方程特征方程的根称为特征根,记作齐次线性差分方程的通解不相等实数根场合有相等实根场合复根场合02211ppppaaap,,,21tpptttcccz2211tpptddtddtcctctccz111121)(tpptititttccececrz3321)(02211ptptttzazazaz1234,,,,,iipaibreaibre为互不相等的实根方法性工具ARMA模型平稳序列建模序列预测非齐次线性差分方程的解非齐次线性差分方程的特解使得非齐次线性差分方程成立的任意一个解非齐次线性差分方程的通解齐次线性差分方程的通解和非齐次线性差分方程的特解之和tttzzztz)(2211thzazazazptpttttz方法性工具ARMA模型平稳序列建模序列预测本章结构方法性工具ARMA模型平稳序列建模序列预测方法性工具ARMA模型平稳序列建模序列预测3.2ARMA模型的性质AR模型(AutoRegressionModel)MA模型(MovingAverageModel)ARMA模型(AutoRegressionMovingAveragemodel)方法性工具ARMA模型平稳序列建模序列预测AR模型的定义具有如下结构的模型称为阶自回归模型,简记为特别当时,称为中心化模型tsExtsEVarExxxxtsstttptptpttt,0,0)(,)(0)(0222110,p)(pAR00)(pAR方法性工具ARMA模型平稳序列建模序列预测AR(P)序列中心化变换称为的中心化序列,令p101ttxy}{ty}{tx方法性工具ARMA模型平稳序列建模序列预测自回归系数多项式引进延迟算子,中心化模型又可以简记为自回归系数多项式)(pARttxB)(ppBBBB2211)(1122----=tttptptxxxx方法性工具ARMA模型平稳序列建模序列预测AR模型平稳性判别判别原因AR模型是常用的平稳序列的拟合模型之一,但并非所有的AR模型都是平稳的方法性工具ARMA模型平稳序列建模序列预测例3.1:考察如下四个模型的平稳性1(1)0.8tttxx1(2)1.1tttxx12(3)0.5ttttxxxttttxxx115.0)4(例3.1平稳序列时序图1(1)0.8tttxx12(3)0.5ttttxxx例3.1非平稳序列时序图1(2)1.1tttxxttttxxx115.0)4(方法性工具ARMA模型平稳序列建模序列预测AR模型平稳性判别方法特征根判别法平稳域判别法方法性工具ARMA模型平稳序列建模序列预测1122121212()0,----=0----=0ttttptppppppBxxxxx即的特征方程为特征根为,,11221212121212122----=0----=0===,,2jjtttptppppppdiwiwjjjjddpmxxxxdrerempdm设齐次线性差分方程的特征方程有p个特征根,,其中为个相等实根为对共轭复根为个互不相等的实根11221212121212122----=0----=0===,,2jjtttptppppppdiwiwjjjjddpmxxxxdrerempdm设齐次线性差分方程的特征方程有p个特征根,,其中为个相等实根为对共轭复根为个互不相等的实根11222'1112111----=0++()jjtjtttptppmdmitwitwjtttjjjjjjjdjxxxxxctcrcece那么齐次线性差分方程的通解为12()()0(B)0111(B)0tpARpBx可以证明,模型的自回归系数多项式的根是齐次差分方程的特征根的倒数,即,,是的根2121221212()11111()11=()0pppppppppBBBB而自回归系数多项式那么11221212()0,----=0----=0ttttptpppppBxxxxx事实上,即的特征方程为1122----=()=tttptptttxxxxBxx设非齐次线性差分方程的一个特解为1()()=(1)piiBBB根据以上性质可以因子分解成111()=(1)===()(1)(1)ptittipttittpiiiiBxBxkxBBB即,那么1122'211121111----=++()(1)jjjtttptpttttpmpdmitwitwjtttijjjjjtjjdjiixxxxxxxkctcrceceB那么齐次线性差分方程的通解为1212(),,,,(1,2,,),lim0,||1,1,2,,2||1,1,2,,()pmjjttiiARpccccjmxipmrimARpp要使得模型平稳,即对任意的实数有该式成立的条件是即要求模型的个特征根都在单位圆内。方法性工具ARMA模型平稳序列建模序列预测AR模型平稳性判别方法特征根判别AR(p)模型平稳的充要条件是它的p个特征根都在单位圆内根据的自回归特征根和自回归系数多项式的根成倒数的性质,等价判别条件是该模型系数多项式φ(B)=0的根都在单位圆外平稳域判别平稳域},,,{21单位根都在单位圆内p方法性工具ARMA模型平稳序列建模序列预测11平稳域〈1特征根111AR(1),-=0tttxx模型特征方程为平稳域21121211224242特征根}11,{12221,且1122212AR(2)+,--=0ttttxxx模型特征方程为1211平稳性条件方法性工具ARMA模型平稳序列建模序列预测例3.1平稳性判别8.010.81.111.1211i212i221210.5,0.5,1.523112312221210.5,1.5,0.5模型特征根判别平稳域判别结论(1)平稳(2)非平稳(3)平稳(4)非平稳方法性工具ARMA模型平稳序列建模序列预测平稳AR模型的统计性质均值方差协方差自相关系数偏自相关系数方法性工具ARMA模型平稳序列建模序列预测均值如果AR(p)模型满足平稳性条件,则有根据平稳序列均值为常数,且为白噪声序列,有推导出p101)(110tptpttxxEExTtEExtt,0)(,}{t方法性工具ARMA模型平稳序列建模序列预测21211221212()1----=(B)0,,,()()0111,(B)0pptttptpttppBBBBxxxxxARpB齐次差分方程即的特征根是模型的自回归系数多项式的根的倒数,即,,是的根方差:借助格林函数计算1122----=()=tttptpttxxxxBx非齐次线性差分方程1()()=(1)piiBBB根据以上性质可以因子分解成111()=(1)===()(1)(1)ptittipttittpiiiiBxBxkxBBB即,那么方法性工具ARMA模型平稳序列建模序列预测Green函数定义AR模型的传递形式其中系数称为Green函数01{1,,1,2,}pjjiiiGGkj1100100()()1ppjtittiitiijipjjiitjjtjjtjijjkxkBBBkGGB0(),()jjttjGBGBxGB记=那么方法性工具ARMA模型平稳序列建模序列预测方差平稳AR模型的传递形式两边求方差得函数为GreenGGxVarjjjt,)(202jtjjtGx0Green函数递推公式原理待定系数法得递推公式pkpkjGGGkkkjjkkj,0,,,2,1110其中,()()()()ttttttBxBGBxGB2121()11ppkpkkBBBBB101)()pkjkjttkjBGB上式写成('111+-))jjjkjkttjkGGB整理得(('1-0,1,2,jjkjkkGGj方法性工具ARMA模型平稳序列建模序列预测平稳AR(1)模型的传递形式平稳AR(1)模型的传递形式为Green函数为itiitiittBBx01011)(1,1,0,1jGjj方法性工具ARMA模型平稳序列建模序列预测求平稳AR(1)模型的方差平稳AR(1)模型的传递形式为Green函数为平稳AR(1)模型的方差itiitiittBBx0