第六章平稳时间序列预测

1、第六章 平稳时间序列预测第一节 平稳时间序列预测概念第二节 最小均方误预测第三节 条件期望预测第四节 适时修正预测第五节 指数平滑预测与arma模型第一节 平稳时间序列预测的概念).(,)0(,1lxltlxtxxxtxttlttttt为预测值记的预测向前期步长为为原点的以这种预测称为进行预测以后的观察值对时刻用序列以前的观察值为及在时刻且零均值平稳序列设当前时刻为返回本节首页下一页上一页一、最小均方误差预测概念二、平稳arma模型最小均方误预测的推导第二节 最小均方误预测(正交投影预测)返回本节首页下一页上一页一、最小均方误差预测概念.)2(min)()( :,) 1 (:)()(:,)0(

2、,)(221值的函数预测值是过去时间序列即预测误差的方差最小足如下两条件准则步最小均方误预测要满所谓其预测误差记为进行的预测对的条件下为已知设lxxeleellxxlelxxxlxtltttlttltttt（若预测函数是线性的，则称线性最小均方误预测）（若预测函数是线性的，则称线性最小均方误预测）返回本节首页下一页上一页二、平稳arma模型最小均方误预测的推导011110111111002211001:1:)()()(:)()(:jjtjltlltlttltltlltltlttttjjtjtttttagagagagagagagagagxltgagagagagabgabbxabxbarma代入上

3、式得将下标其中如下此模型写成其传递形式模型如下设有平稳返回本节首页下一页上一页由于预测只能建立在到t时刻为止的可用信息的基础上，因此，根据最小均方误预测的第二个准则，以及平稳可逆序列可以表示成传递函数形式的论断，可以将预测值 表示成能够估计的项at，at-1，的加权和的形式：)( lxt2*21*1*0*)(tltltljjtjltagagagaglx.*达到最小的意义下确定可以在预测误差的方差系数权式中jlg由上得以t为原点，向前l步的预测误差为：0*11110)()()(jjtjljltlltlttlttaggagagaglxxle由于at是白噪声，故有：002jjoaeaajtt02*2

4、102222)()()(:jjljlaljjattltgggleelxxe预测误差的方差为所以.,:*上式达到最小值时当很容易看出jljlgg因此可得xt+l的最小均方误预测为：2211)(tltltltagagaglx预测误差为：1110)()(tlltlttlttaggaglxxle误差方差为：102221222122)1 ()(ljjalatgggglee由上推导可知，（1）最小均方误预测误差的方差和预测步长l有关，而和预测的时间原点无关。（2）预测步长l越大，预测误差的方差也越大，即预测的准确性越差。上述最小均方误预测公式中包含有无穷项求和，而在实际中我们只可能有有限的数据，因此，只能

5、用充分多项的有穷和近似，即：.)(10小于允许值即可的取值只要使ttjtjttjjtjttagtaglx第三节 条件期望预测一、条件期望预测的一般公式二、用条件期望进行预测三、arma(p,q)模型条件期望预测的一般结果四、arma(p,q)条件期望预测的置信区间返回本节首页下一页上一页一、条件期望预测的一般公式值的条件期望作为其预测取条件下指在已知条件期望预测)0(,:21lxxxxltttt用公式表示如下：),|()(21tttlttxxxxelx返回本节首页下一页上一页有关xt和at的条件期望有如下性质：)0()()0(),|() 1 (21llxlxxxxxetlttttlt)0(0)

6、0(),|()2(21llaxxxaelttttlt由于：0111101111110jjtjltlltlttltltlltltltagagagagagagagagagx利用条件期望的性质，对上式两端求条件期望，得xt+l 的条件期望预测为：221121),|()(tltltltttlttagagagxxxxelx可见， xt+l 的条件期望预测和它的最小均方误预测是一致的。二、用条件期望进行预测1.ar(1)模型的条件期望预测(参见p130)设xt适合如下ar(1)模型:tttaxx11(1)以t为原点，向前一步预测公式(l=1)ttttttttttttttttttttttxxxxaexxxx

7、exxxaxexxxxexaxx1211211211211111),|(),|(),|(),|() 1 (:得对上式两端求条件期望由返回本节首页下一页上一页(2) 向前二步预测公式(l=2) 1 (),|(),|()2(1212112122112ttttttttttttttxxxxaxexxxxexaxx(3) 向前l步预测公式(l2) 1(),|(),|()(1211121111lxxxxaxexxxxelxaxxttttltlttttlttltltlt由上推导可见，对于l0，条件期望预测值 满足如下差分方程： )( lxtltttttltttttxlxxcxxclxxxlxlx111111

8、)(:) 1 (:)(:0:)0(0) 1()(所以得由于是即得为此差分方程的特征方程2、ar(2)模型的条件期望预测n设xt适合如下ar(2)模型:ttttaxxx2211(1)以t为原点，向前一步预测公式(l=1)1212111221111211),|(),|() 1 (tttttttttttttttttxxxxxaxxexxxxexaxxx由(2) 向前二步预测公式(l=2)tttttttttttttttttxxxxxaxxexxxxexaxxx2121221121222112) 1 (),|(),|()2(3) 向前l步预测公式(l3)2() 1(),|(),|()(212122112

9、12211lxlxxxxaxxexxxxelxaxxxtttttltltlttttlttltltltlt可见，当l1时，ar(2)预测值可由如下差分方程求出：0)2() 1()(21lxlxlxttt(预测值的一般解略)3、arma(1,1)模型的条件期望预测设1111ttttaaxx(1) 向前一步预测 (l=1)tttttttttttttttttaxxxxaaxexxxxexaaxx11211112111111),|(),|() 1 (:得由(2) 向前二步预测 (l=2) 1 (),|(),|()2(12111211212112112ttttttttttttttttxxxxaaxexxx

10、xexaaxx(3) 向前l步预测公式(l2) 1(),|()(1211111lxxxxxelxaaxxttttlttltltltlt可见，当l1时，arma(1,1)预测值也是由如下差分方程决定的。0) 1()(1lxlxttltclx1)(解得：由于：tttaxcx111) 1 (所以：ttaxc11因此：ltttaxlx111)(4、ma(1)模型的条件期望预测n设11tttaax(1) 向前一步预测 (l=1)ttttttttttttttaxxxaaexxxxexaax12111211111),|(),|() 1 (:得由(2) 向前二步预测 (l=2)0),|()2(2121122t

11、tttttttxxxxexaax(3) 向前l步预测公式(l2)0),|()(21tttlttxxxxelx设：qtqttptptttaaaxxxx112211（1）向前一步预测(l=1)：qiqtitpiititaaxx111111对上式两端求条件期望得：qtqtptpttttttaaxxxxxxex11112111),|() 1 (三、arma(p,q)模型条件期望预测的一般结果返回本节首页下一页上一页(2)向前二步预测公式(l=2)qiqtitpiititaaxx122122qtqtptpttttttaaxxxxxxex2222112) 1 (),|()2(3)向前l步预测公式(lp，且

12、l q)qiqltiltpiiltiltaaxx11qltqtlpltptltlttttlttaaxxxlxlxxxxelx) 1 ()2() 1(),|()(1211(3)向前l步预测公式(lp，且l q)qiqltiltpiiltiltaaxx11)()2() 1(),|()(211plxlxlxxxxelxtpttttltt由推导可以看出，对于arma(p,q)模型的向前l 步预测(lp，且l q)，预测结果满足如下差分方程：0)()2() 1()(21plxlxlxlxtpttt（预测值解的一般形式参见课本p134）由解的一般形式可以看出，对于arma(p,q)模型，自回归部分决定了预

13、测函数的形式，而滑动平均部分则用于确定预测函数中的系数。预测举例：例1：利用对zl14所建立的模型进行预测。先对原序列零均值化，先对原序列零均值化，然后建模如下：ttttaxxx2122. 079. 0已知：78. 3,58. 4249250 xx)()3(),2(),1 (:250250250250lxxxx以及预测函数求42. 9,yyyxyttt则设原序列为251249250125022. 079. 0axxx解：787. 278. 322. 058. 479. 022. 079. 0),|() 1 (249250249250251250 xxxxxex194. 158. 422. 07

14、87. 279. 022. 0) 1 (79. 0),|()2(250250249250252250 xxxxxex同理：33. 0787. 222. 0194. 179. 0) 1 (22. 0)2(79. 0),|()3(250250249250253250 xxxxxex250250250250250250( )0.79(1)0.22(2):( )0.79(1)0.22(2)0 xlxlxlxlxlxl即2l当 时，预测值满由模型自回归部分决定的差分方程：解此差分方程即可求出预测函数。前已证明，条件期望预测与最小均方误预测是一致的，因此，预测误差和误差方差也是相同的。因此，条件期望的预测

15、误差为：01111( )( )tt ltt lt llte lxx lg ag ag a 四、arma(p,q)条件期望预测的置信区间返回本节首页下一页上一页1,)()(,) 1, 2 , 1(.0012gbgbbljgajjjjta且确定可由为格林函数的方差为白噪声预测误差的方差为：102221222122)1 ()()(ljjalattggggleeled其中：l =1时的预测误差为：1)()(ttlttalxxle于是有：) 1 () 1 (11ttttxxea可见可见arma模型中白噪声项模型中白噪声项at其实就是以其实就是以xt-1为原点，为原点，向前一步预测误差。向前一步预测误差。

16、预测误差和白噪声项的关系：再由预测误差方差的公式得：22122)()()1 ()1 (ttattaeaeeeed可见：向前一步预测误差的方差其实就是白噪声项的方差。可见：向前一步预测误差的方差其实就是白噪声项的方差。预测误差的置信区间：对于正态过程，预测误差的分布为：21212221)1 (96. 1)(latggglx所以：对所以：对xt+l预测的预测的95%的置信区间为：的置信区间为：)(, 0()(lednlett因此：)(),(),|(1ledlxnxxxttttlt)1 ()()(21222122lattgggleeled222121196. 1) 3(:,3196. 1)2(:,2

17、96. 1) 1 (:%951ggxlgxlxlatatat有时当有时当的预测区间为时的当根据预测置信区间的公式得：可见：随着预测步长的加大，预测误差的置信区间也越大，预测结果越不准确。例1：zl14磨轮剖面数据，所建模型如下：404. 022. 079. 079. 079. 01:)2(22. 0,07945. 222. 079. 021121102121ggggaraxxxatttt模型的格林函数得于是由由于802. 4787. 245. 296. 1787. 296. 1) 1 (:%951250axl的预测区间为时的当12. 6194. 179. 0145. 296. 1194. 11

18、96. 1)2(:,2221gxlat有时当42. 633. 0404. 079. 0145. 296. 133. 0196. 1) 3(:,3222221ggxlat有时当于是以t=250为原点，向前一步、二步、三步预测的95%的置信区间分别为：yxyyyyxtttt所以因为,42. 9,所以对于原序列，以t=250为原点向前一步，二步、三步的预测分别为：42. 675. 9)3(12. 6614.10)2(802. 4207.12) 1 (250250250yyy例2. 对arma21.wf1文件中的序列x建模如下：12158. 055. 097. 0tttttaaxxx58. 0,55.

19、 0,097121则已知：083. 0,336. 2,008. 1,024. 0,172. 0248247248249250axxxx的置信区间以及预测值求%95)2(),1 (:250250 xx模型的剩余平方和为260.04。（1）求预测值解：)3(),2(),1 (250250250 xxx25024925024925025025124925024925025125058. 055. 097. 0),|58. 055. 097. 0(),|() 1 (axxxxaaxxexxxexa250未知，故需先将其求出。由已知数据得：379. 0)083. 0(58. 0)34. 2(55. 0)

20、008. 1(97. 0024. 058. 055. 097. 0248247248249249axxxa578. 0)379. 0(58. 0)008. 1(55. 0)0236. 0(97. 0172. 058. 055. 097. 0249248249250250axxxa同理：因此：51. 0)578. 0(58. 0)024. 0(55. 0172. 097. 058. 055. 097. 0) 1 (250249250250axxx40. 0172. 055. 051. 097. 055. 0) 1 (97. 0),|58. 055. 097. 0(),|()2(250250249

21、250251252250251249250252250 xxxxaaxxexxxex（2）求预测值的95%的置信区间：03. 13)2250(04.260a12213021121110172. 055. 039. 097. 039. 058. 097. 01gggggggg由arma(2,1)模型的格林函数得：所以预测值的95%的置信区间为：02. 251. 003. 196. 151. 096. 1) 1 (250ax17. 240. 039. 0103. 196. 140. 0196. 1)2(221250gxa在eviews中利用arma模型进行预测。（1）eviews中进行预测时的两个

22、选项。dynamic动态预测。（含义）static一步超前预测。（含义）对于arma模型：若对序列进行拟合分析(即追溯预测)，则选static。若向前l步预测，则要选dynamic，并且要先对工作区间、样本区间进行调整如下： (1)expand first t+l (2)smpl t+1 t+l具体操作见演示。(2)通过eviews计算预测置信区间。例：根据磨轮剖面数据zl14.wf1，(1)建立模型。(2)模型追溯预测分析。(3)进行外推预测(l=3).第四节 适时修正预测一、问题的提出二、适时修正预测公式返回本节首页下一页上一页一、问题的提出.,),3(),2(, 3, 2,1,) 3()

23、,2(),1 (,111于是提出适时修正预测信息才是最新的而并未利用时刻及以前的信息因为它们只利用了就不能再用原来的等时刻的预测这时对于已成为已知的时候而当到了时刻可得到为原点进行向前预测以时刻在实际进行预测时ttttttttxxtxxttxtxxxt返回本节首页下一页上一页二、适时修正预测公式1、适时修正预测公式的推导（1）适时修正预测公式)1 () 1()(11ttlttxxglxlx返回本节首页下一页上一页2、适时修正预测公式的推导：112312112111) 1() 1 (:) 1()(:tltttltltlttltltltaglxxagagaglxagagaglx于是可得由最小均方误

24、预测公式1110) 1 () 1 (:,1,)()(:,ttttljjltjtlttaxxelaglxxleltarma向前一步预测误差为时当特别地预测误差为步的最小均方误向前为原点模型以时刻又已知综合上述推导，可得适时修正预测公式为：)1 () 1()(11ttlttxxglxlx上述公式说明：新的预测值是在旧的预测值的基础上，加上一个修正项推算出来的，而这一个修正项比例于旧的一步预测误差，比例系数随着预测超前步数的变化而变化。)2(),1 (:0 . 3252251251251xxyyx试计算假设已知例：对于 zl14磨轮剖面数据，解：)1 () 1()(11ttlttxxglxlx适时修

25、正预测公式为：所以：362. 1)787. 23(79. 0194. 1)1 ()2() 1 (2502511250251xxgxx416. 0)787. 23(404. 033. 0)1 ()3()2(2502512250251xxgxx第五节 指数平滑预测与arma模型一、一次指数平滑预测的原理二、arima(0,1,1)模型的预测返回本节首页下一页上一页一、一次指数平滑预测的原理一次指数平滑预测的基本公式为：)1 ()1 ()1 ()1 ()1 (111ttttttxxxxxx其中：01为平滑系数。将上述公式展开得：)1 ()1 ()1 () 1 (21ttttxxxx如此展开下去可得：33221)1 ()1 ()1 () 1 (tttttxxxxx返回本节首页下一页上一页设有arima(0,1,1)模型如下:ttabxb)1 ()1 (ttttxbbbxbbbxbba)1 ()1 ()1 (1 )1)(1 (1132222将其表示成逆转形式得：返回本节首页下一页上一页二、arima(0,1,1)模型的预测上式即为:3221)1 ()1 ()1 (tttttxxxax对其作向前一步预测可得：22111)1 ()1 ()1 (),|() 1 (tttttttxxxxxxex令1-=，上式可变为：221)1 ()1 () 1 (ttttxxxx其中， = 1-1由