时间序列分析讲义第02章滞后算子

1、第二章滞后算子及其性质滞后算子是对时间序列进行动态线性运算的主要工具，利用滞后算子可以使得一些非线性运算非常简洁。 2.1基本概念时间序列是以观测值发生的时期作为标记的数据集合。一般情况下，我们是从某个特定的时间开始采集数据，直到另一个固定的时间为止，我们可以将获得的数据表示为：(yi, y2,t)如果能够从更早的时间开始观测，或者观测到更晚的时期，那么上面的数据区间可以进一步扩充。相对而言，上述数据只是一个数据的片段，整个数据序列可以表示为：(,yi,y2,，yT,)二yjH例2.1几种代表性的时间序列(1) 时间趋势本身也可以构成一个时间序列，此时：yt二t ；(2) 另一种特殊的时间序列

2、是常数时间序列，即：yt =c，c是常数，这种时间的取值不受时间的影响；(3) 在随机分析中常用的一种时间序列是高斯白噪声过程，表示为：yt二；t , ；Jt參是一个独立随机变量序列，每个随机变量都服从N (0,匚2)分布。时间序列之间也可以进行转换，类似于使用函数关系进行转换。它是将输入时间序列转换为输出时间序列。例2.2几种代表性的时间序列转换(1) 假设xt是一个时间序列，假设转换关系为：yt = 1人，这种算子是将一个时间序列 的每一个时期的值乘以常数转换为一个新的时间序列。(2) 假设Xt和wt是两个时间序列，算子转换方式为：yt二人 Wt，此算子是将两个时间序列求和。定义2.1如果

3、算子运算是将一个时间序列的前一期值转化为当期值，则称此算子为滞 后算子，记做L。即对任意时间序列 xt，滞后算子满足：L(xt)=人 4(1)类似地，可以定义高阶滞后算子，例如二阶滞后算子记为L2，对任意时间序列xt，二阶滞后算子满足：L2(xt)三LL(xJ =Xt,一般地，对于任意正整数 k，有：Lk (xt)二 Xt 止(3)命题2.1滞后算子运算满足线性性质：(1) L(: X：L(Xt)(2) L(Xt Wt)=L(xJ L(wJ证明：(1)利用滞后算子性质，可以得到：L(: Xt)= “斗二：L(Xt)(2) L( xt wt)=Xtnwt 1 =L(xt)L(wt)End由于滞后

4、算子具有上述运算性质和乘法的交换性质，因此可以定义滞后算子多项式，它的作用是通过它对时间序列的作用获得一个新的时间序列，并且揭示这两个时间序列之间的关系。显然，滞后算子作用到常数时间序列上，时间序列仍然保持常数，即:L(c) = c。 2.2 一阶差分方程利用滞后算子，可以将前面的一阶差分方程表示成为滞后算子形式:yt 二 ytwt 二 L yt wt也可以表示为：(1 _ L)% 訥在上述等式两边同时作用算子：(1 L J 2L2亠亠 Lt)，可以得到：(1 Lt U)(1 L)ytL 川卷过 U)wt计算得到：(1 t 1 U 1)yt =(1 Lt Lt)wt利用滞后算子性质得到：yt

5、二 t1y 二 WtWt亠 亠亠 “tw0(6)上述差分方程的解同利用叠代算法得到的解是一致的。注意到算子作用后的等式：(1 . L 亠汐 U )(1 - L)yt =yt - 4 1 yt如果时间序列yt是有界的，即存在有限的常数M，使得任意时间均有：|yt，并且：1，则上式当中的尾项随着时间增加趋于零。从而有：lim(1 L讥 U)(1 - L)yt(7)如果利用“1”表示恒等算子，则有:(8)记(需要注意的是，这里只是表示一个运算符号)：(1 一 L)二=im(1L 亠 亠Lt)因此得到了 “逆算子”的表达式，这类似于以滞后算子为变量的函数展开式。定义2.2当：1时，定义算子(1 - L

6、)的逆算子为(1 - L) J，它满足:(1) (1 一 1)(1 L)，=(1 L) 1 一 L)=l其中I表示单位算子，即对任意时间序列yt，有：I (yt) =yt(2) 在形式上逆算子可以表示为：(9)(10)这表示逆算子作为算子运算规则是：对于任意时间序列yt ,有:(1 1)4%jLj(yt)jyt_jj =0j=0当| J-1时，逆算子( L)4的定义以后讨论。如果时间序列yt是有界的，则一阶差分方程的解可以表示为:yt 二 wt wt 42 wt,j wt :j=0(11)可以验算上述表达式确实满足一阶线性差分方程。但是解并惟一，例如对于任意实数 a。，下述形式的表达式均是方程

7、的解。yt Wj %j=0上述差分方程的解中含有待定系数，这为判断解的性质留出一定的余地。 2.3二阶差分方程我们考察二阶差分方程的滞后算子表达式：yt2 yt wt将其利用滞后算子表示为：（1 一 i L 一 2 L2）yt =wt（12）对二阶滞后算子多项式进行因式分解，即寻求 1和2使得：（1 1 L 2L2）=（1 r L）（1 2 L）二1 （2）L t2 L2显然r和2是差分方程对应的特征方程的根：丸2 电丸鴨=0（13）当特征根-1和-2落在单位圆内的时候（这也是差分方程的稳定性条件），滞后算子多项式分解为：（1 - h L）=1 ：% L ：.借 L2 ： *3 L3 ，（1

8、苍 L）=1 .；”2 L 恳 L2 - 卷 L3这时二阶差分方程解可以表示为：yt =（1 一 丄）（1 一 2L） dwt注意到算子分式也可以进行分项分式分解（如此分解需要证明，参见Sargent, 1987, p.184):（1 -入 L）（1 上2 L）（人一吃）i（1 人 L） （1 -力2 L）丿(人2)l(1 人 L)(1 2L)丿1-1 丄？L2 亠 亠 1 tL2L2wt(1 - 2 )=(C C2 ) Wt (C11 C2 ) Wt 4 (C 2 C2 2 ) Wt _21 2将上述表达式带入到二阶差分方程解中：其中：C1 = ， C2 =1 _ 2 2 1利用上述公式，可

9、以得到外生扰动的动态反应乘子为：=c1j C22，j =0,1,.wt上述利用滞后算子运算得到的乘数与以前所得完全一致。例2.3对于二阶差分方程而言，其特征方程是：2 1 - 2 =0得到特征根为：A A1 匕（1 42），2 =-（ . 1 4 2）22(14)上述方程的稳定性与滞后算子多项式的根落在单位圆内是一致的。 2.4 p阶差分方程上述算子多项式的分解方法可以直接推广到p阶差分方程情形。将p阶差分方程表示成为滞后算子形式：（1 1 L 2L2i*LP）yt =Wt将上式左端的算子多项式分解为：（1 - L 2L2 丄ipLP）=：（1 、L）（1 2L） 1 这相当于寻求厂P）使得下

10、述代数多项式恒等:（1 一 Z 2Z2 -一 pZP） = （1 1 Z）（1 JZ）I定义，=z，则可以将上述多项式表示成为：（P - -P1 -2 上.p）=（，一、）（，- 这意味着算子多项式的分解，就相当于求出差分方程特征方程的根。 如果差分方程的根相异，且全部落在单位圆内，则可以进行下述分式分解:1CiC2CpCi是:(15)(1 - pL)(16)(1 -先 p z)(17)(18)C2+（1 一 1 L）（1 一 jL） （1 一，pL） （1 -L） （1 一 丄）（1_，pL）通过待定系数法，可以得到上述分式中的参数为：g 二7pr， j =1,2, P1 2-.p显然有：C

11、 c2 亠亠cp =1利用上述算子多项式分解，可以得到差分方程的解为：1y.Wt(1 _ L _ 2 L2 - PLP)C1C2cp(19)(20)(21)wtwtwt(1-、L)(1 - 2L)(1- pL)(1 L xiL2 )wt c2 (1 r2 L r2L2 )wt+cp(1 +p L +屁L2C2 Cp)Wp (C11 C2,2Cpp)Wpjg1 s2Cpp)Wp_通过上述方程通解，可以得到动态反应乘子为：+C2雪 + +CpP，j=1,2,；wt命题2.2外生变量Wt对yt现值的影响和外生变量Wt持续扰动对yt的动态影响乘子(22)(23).：Wtyt j-p：.：yt j:-W

12、t 1- Wt j证明：将差分方程的解表示为：yt仝mWtN訥厂 ，其中：:j =Ci Cpp , j =1,2,设：(L) = 丄2L2 3L3 利用算子多项式表示：yt 仝(L)wWt对y现值的影响可以表示为：Jyt.jT-j =；Xjj -C):Wt j =0j =Q； Wtj 卫注意到：(L)=%打L L2 ;：3L珂(1 一 丄)(1 一 pL)因此有：:C)=(1 1 ：) (1 ，p )=1 /： 2 ：2 一 -P-1长期乘数相当于2 T的情形，从而得到公式所示的公式。End上述命题结论是利用滞后算子多项式推导的，其结论同利用差分方程矩阵表示所得到的结论是一致的。 2.5初始条

13、件和无界序列假设给定下述线性差分方程：yt =1 yt 2yt/ 亠 亠pytT wt(24)一般情况下，求解p阶差分方程的特解，需要p个初值：y，y?-,y4，也需要外生变量的一个输入序列：wo,叫,，Wt，这样一来根据差分方程结构，便可以确定yt的时间路径。但是，在一些常见的经济或者金融时间序列当中，无法给定具体的初值或者完整的外生输入变量，那么这时差分方程解的性质如何？例2.4假设变量Pt表示股票价格，Dt表示股票派发的红利。如果一个投资者在时刻t买入股票，然后在时刻t 1卖出股票，则他将获得实际红利收入Dt d /Pt和价格收益(R+_R)/Pt，因此投资者的收益率为：rt 1 =(R

14、 ! -PJ/R 口 Pt(25)在简单的股票市场模型当中，假设收益率是常数，则上述方程可以转化为股票价格的差分方程模型：Pt -(1 * r )Pt 4、 Dt(26)如果知道红利序列D1, D2 ,，Dt和股票价格的初值 Pg，则可以得到股票价格路径为：P =(1 r)tP3 -(1 )2。1 一(1 r)tD2 -Dt(27)但是如果仅仅知道红利序列，而不知道股票价格初值，则可能有很多价格轨迹满足价格 的差分方程。为了说明这个问题，进一步假设红利为常数，则有：P =(1 r)tP -D(1 r) (1 r)t - 1丄t1(1+r)t=(1 r)tPD1 _(1 +r)=(1 r)tP0

15、 -(D/r)(D/r)(28)(1) 如果初始时期股票价格等于红利贴现，即Po =D/r，则有：R D/r , t =0,1,2,此时股票价格保持常数，股价等于红利除以收益率。这种股票价格被称为在收益率是常 数情形的股价基础成分。(2) 假设初始股价超过了 D/r，即Po D/r，这时股票价格出现了扩散现象，这与资产定价理论相符。因为为了保持资产收益率不变，股票的价格就会出现持续上升，同时假设红利是固定的，红利带来的实际收益减少将被股价的加速增长所弥补，这样就出现了股票价格膨胀的现象，即出现股票价格泡沫。(3)为了消除股价中的投资泡沫，一种方法是对股票价格路径给予有界性限制。例如， 假设对于

16、所有时期的股票价格满足：Ft P， t =0,1,2,这样一来，满足上述约束的股票价格路径便是常数的市场基础价格。 上面假设了常数红利，现在假设红利序列是有界的。将股价表示为：1R =k【R + +Dt+1 r进行向前叠代运算有：+ J*1 IH rR T 1 r j如果价格序列R满足约束条件：lim Rtt 1 - r在假设Dt和R均是有界序列，则得到股票价格水平满足：二 1R _ j =o|1 - r这是红利随时间变化时股票价格的市场基础成分。需要注意的是，对于上述情形的市场基础成分， 需要投资者对于未来红利具有完全预期。 当引入预期红利时，上述表达式仍然适用，这时可以修改为：R 川1(

17、j =o 1 r利用红利预期的股价公式，可以确定价格初值:Ro =?r7如此初值是否满足一般的股价模型，我们可以代入到具有初值的确定解中验证：R 二(1 r)t Rj -(1 r)tD1 -(1 r)tD - Dt将Po代入上式后得到：R八 1R j_1 r这正是在边界条件下所推导的向前预期解，由此可见该解与初值选择是吻合的。例2.5我们继续利用滞后算子方法讨论股票价格路径的性质。利用算子表示为：1 -(1 +r)LR =-Dt在上述表达式中，滞后算子多项式的特征根小于1,无法采用逆算子的一般表达式，为TJ1Dm 亍Dt T =0jDt jDt j)DjjDt j此我们需要采取新的定义。定义滞

18、后算子L的逆算子为L，具有性质:(1) LJ L =LL=1(2) L(yJ 二 i这样一来，滞后算子乘积就具有幕乘的性质：对于任意正整数i和j，有：Li L-j =Li-j对方程(2.12)两端乘以算子多项式：1 1 L 二一1L 之|(1 r) (1 r)21 1*(1 +r)2*(1 +r)T Df当r 0，且红利序列是有界的，则上述极限为：11(1 r)2pFDt整理得到：qQPt- Dt jtj=o(1 r)j t j根据上述运算，可以定义下述算子的逆算子：1 1 - 1 . 11-(1 r)L|_(1r)L |L (1 r)L (1 r)2L2 2.6差分方程的求解方法上面我们主要

19、论述了差分方程的表示和外生扰动的动态乘子，下面我们给出差分方程的一般求解过程。第一步：构造p阶齐次差分方程，并且寻求齐次方程的p个解：yth , i=1,2，p第二步：构造p阶非齐次差分方程的特解。第三步：齐次方程p个解的线性组合加上非齐次方程的一个特解，得到非齐次方程的通解。第四步：根据给定的边界条件，确实通解当中的未知参数，得到非齐次方程的确定解。2.6.1齐次差分方程的通解和稳定性p阶齐次差分方程的形式是：yt = yt2yta 亠.亠pyt*命题2.3对于差分方程而言，下述推断成立：(1) 如果yth是方程的解，则对任意常数A，Ayth也是解。(2) 如果y1t和y2t是方程的解，则对任意实数 A和A2，A1 y1t - A2 y2t也是方程的解。 证明：留做练习。对于p阶齐次差分方程，我们尝试地检验解的形式是：yth二A，代入差分方程为：- J p4 _ 2 PQ - =0由此可见，应该是上述特征方程的根。因此，如果差分方程具有相异实数根的时候， 可以得到p个解为：yR二f，i =1, 2 , p，此时解的稳定性要求所有根落在单位圆内。命题