1、时间序列分析方法讲义 第 2 章 滞后算子及其性质1第二章 滞后算子及其性质滞后算子是对时间序列进行动态线性运算的主要工具,利用滞后算子可以使得一些非线性运算非常简洁。2.1 基本概念时间序列是以观测值发生的时期作为标记的数据集合。一般情况下,我们是从某个特定的时间开始采集数据,直到另一个固定的时间为止,我们可以将获得的数据表示为: ),(21Ty如果能够从更早的时间开始观测,或者观测到更晚的时期,那么上面的数据区间可以进一步扩充。相对而言,上述数据只是一个数据的片段,整个数据序列可以表示为: tTy),(21例 2.1 几种代表性的时间序列(1) 时间趋势本身也可以构成一个时间序列,此时:
2、;ty(2) 另一种特殊的时间序列是常数时间序列,即: , 是常数,这种时间的取值ct不受时间的影响;(3) 在随机分析中常用的一种时间序列是高斯白噪声过程,表示为: ,ty是一个独立随机变量序列,每个随机变量都服从 分布。t ),0(2N时间序列之间也可以进行转换,类似于使用函数关系进行转换。它是将输入时间序列转换为输出时间序列。例 2.2 几种代表性的时间序列转换(1) 假设 是一个时间序列,假设转换关系为: ,这种算子是将一个时间序tx ttxy列的每一个时期的值乘以常数转换为一个新的时间序列。(2) 假设 和 是两个时间序列,算子转换方式为: ,此算子是将两个时ttwttw间序列求和。
3、定义 2.1 如果算子运算是将一个时间序列的前一期值转化为当期值,则称此算子为滞后算子,记做 。即对任意时间序列 ,滞后算子满足:Ltx(1)1)(ttx类似地,可以定义高阶滞后算子,例如二阶滞后算子记为 ,对任意时间序列 ,二2Ltx阶滞后算子满足:(2)22)()(ttt xL一般地,对于任意正整数 ,有:k(3)kttkx命题 2.1 滞后算子运算满足线性性质:(1) )()(ttxL(2) )(tw证明:(1) 利用滞后算子性质,可以得到:)(1tttxL(2) End)(1ttLxx由于滞后算子具有上述运算性质和乘法的交换性质,因此可以定义滞后算子多项式,它的作用是通过它对时间序列的
4、作用获得一个新的时间序列,并且揭示这两个时间序列之时间序列分析方法讲义 第 2 章 滞后算子及其性质2间的关系。显然,滞后算子作用到常数时间序列上,时间序列仍然保持常数,即: 。cL)(2.2 一阶差分方程利用滞后算子,可以将前面的一阶差分方程表示成为滞后算子形式:(4)ttttt wyLy1也可以表示为:(5)tt)(在上述等式两边同时作用算子: ,可以得到:)1(2tLL tttt wyL)(1计算得到: tttty)1)(利用滞后算子性质得到:(6)01wy ttttt 上述差分方程的解同利用叠代算法得到的解是一致的。注意到算子作用后的等式: tttt yyLL1)( 如果时间序列 是有
5、界的,即存在有限的常数 ,使得任意时间均有: ,tyMMyt|并且 ,则上式当中的尾项随着时间增加趋于零。从而有:1|(7)tttt )1(lim如果利用“1”表示恒等算子,则有:(8)(li LLtt记(需要注意的是,这里只是表示一个运算符号) :(9)1litt 因此得到了“逆算子”的表达式,这类似于以滞后算子为变量的函数展开式。定义 2.2 当 时,定义算子 的逆算子为 ,它满足:|)1(L1)(L(1) (10)IL)()(1其中 表示单位算子,即对任意时间序列 ,有:I tytty(2) 在形式上逆算子可以表示为:(11)01)(jj这表示逆算子作为算子运算规则是:对于任意时间序列
6、,有:ty001)()( jjtjtjt yLy当 时,逆算子 的定义以后讨论。|1如果时间序列 是有界的,则一阶差分方程的解可以表示为:t 021jjttttt wwy可以验算上述表达式确实满足一阶线性差分方程。但是解并惟一,例如对于任意实数,下述形式的表达式均是方程的解。0a时间序列分析方法讲义 第 2 章 滞后算子及其性质30jjttt way上述差分方程的解中含有待定系数,这为判断解的性质留出一定的余地。2.3 二阶差分方程我们考察二阶差分方程的滞后算子表达式: tttt wyy21将其利用滞后算子表示为:(12)ttL)(对二阶滞后算子多项式进行因式分解,即寻求 和 使得:12)()
7、1)(1 2122 LL 显然 和 是差分方程对应的特征方程的根:(13)022当特征根 和 落在单位圆内的时候( 这也是差分方程的稳定性条件 ),滞后算子多项1式分解为:,312)( LL212这时二阶差分方程解可以表示为: tt wy11)()(注意到算子分式也可以进行分项分式分解(如此分解需要证明,参见 Sargent,1987,p. 184): )1()()()1)( 2212 LL将上述表达式带入到二阶差分方程解中: 2212121 2121 )()()( )()()( ttt ttt wcwcwc LLy 其中: , 利用上述公式,可以得到外生扰动的动态反应乘子为:, (14)jj
8、tjcwy21,10上述利用滞后算子运算得到的乘数与以前所得完全一致。例 2.3 对于二阶差分方程而言,其特征方程是: 0212得到特征根为:,)4(211)4(121上述方程的稳定性与滞后算子多项式的根落在单位圆内是一致的。时间序列分析方法讲义 第 2 章 滞后算子及其性质42.4 p 阶差分方程上述算子多项式的分解方法可以直接推广到 p 阶差分方程情形。将 p 阶差分方程表示成为滞后算子形式:(15)ttpwyLL)1(2将上式左端的算子多项式分解为:(16)1()(1221 Lpp 这相当于寻求 使得下述代数多项式恒等:),(1(17)( 1zzzz定义 ,则可以将上述多项式表示成为:(
9、18)()(221 pppp 这意味着算子多项式的分解,就相当于求出差分方程特征方程的根。如果差分方程的根相异,且全部落在单位圆内,则可以进行下述分式分解:(19)1()1()()1()(1( 22 LcLcLL pp 通过待定系数法,可以得到上述分式中的参数为:, (20)1121ppiic pj,显然有:(21)pc利用上述算子多项式分解,可以得到差分方程的解为:(22) jpjjj pp tptt tptt tpt wcc wcLLcccwy)( )()1(1( )1()(21 1212221通过上述方程通解,可以得到动态反应乘子为:, (23)jpjjtjwy21 ,命题 2.2 外生
10、变量 对 现值的影响和外生变量 持续扰动对 的动态影响乘子ttytwty是: pjjtty210 pjttjtjj wyw 21lim证明:将差分方程的解表示为:,321tttty其中:时间序列分析方法讲义 第 2 章 滞后算子及其性质5,1jpjj c ,21设: 320)( LL利用算子多项式表示: ttwy对 现值的影响可以表示为: 000 )(jjjtjtt y注意到: 113210 )()()( LLL p因此有: 21)()p 长期乘数相当于 的情形,从而得到公式所示的公式。 End上述命题结论是利用滞后算子多项式推导的,其结论同利用差分方程矩阵表示所得到的结论是一致的。2.5 初
11、始条件和无界序列假设给定下述线性差分方程:(24)tptttt wyyy21一般情况下,求解 p 阶差分方程的特解,需要 p 个初值: ,也需要外pyy,21生变量的一个输入序列: ,这样一来根据差分方程结构,便可以确定 的时t,10 t间路径。但是,在一些常见的经济或者金融时间序列当中,无法给定具体的初值或者完整的外生输入变量,那么这时差分方程解的性质如何?例 2.4 假设变量 表示股票价格, 表示股票派发的红利。如果一个投资者在时刻tPtD买入股票,然后在时刻 卖出股票,则他将获得实际红利收入 和价格收益t 1 ttPD/1,因此投资者的收益率为:tP/)(1(25)ttttr/)(1在简
12、单的股票市场模型当中,假设收益率是常数,则上述方程可以转化为股票价格的差分方程模型:(26)ttt Dr1)(如果知道红利序列 和股票价格的初值 ,则可以得到股票价格路径为:,2t 0P(27)ttttt DrrPr 210)()()(但是如果仅仅知道红利序列,而不知道股票价格初值,则可能有很多价格轨迹满足价格的差分方程。为了说明这个问题,进一步假设红利为常数,则有:(28)/()/()1(11)()()(020rDPrt tt ttt (1) 如果初始时期股票价格等于红利贴现,即 ,则有:rDP/0时间序列分析方法讲义 第 2 章 滞后算子及其性质6,rDPt/,210t此时股票价格保持常数
13、,股价等于红利除以收益率。这种股票价格被称为在收益率是常数情形的股价基础成分。(2) 假设初始股价超过了 ,即 ,这时股票价格出现了扩散现象,这与资rD/rP/0产定价理论相符。因为为了保持资产收益率不变,股票的价格就会出现持续上升,同时假设红利是固定的,红利带来的实际收益减少将被股价的加速增长所弥补,这样就出现了股票价格膨胀的现象,即出现股票价格泡沫。(3) 为了消除股价中的投资泡沫,一种方法是对股票价格路径给予有界性限制。例如,假设对于所有时期的股票价格满足:,Pt,210t这样一来,满足上述约束的股票价格路径便是常数的市场基础价格。上面假设了常数红利,现在假设红利序列是有界的。将股价表示
14、为: 11ttt Dr进行向前叠代运算有: 111 tTtTtTtt DrDrrPr 如果价格序列 满足约束条件:t01limTtTr在假设 和 均是有界序列,则得到股票价格水平满足:tDtP0j jtjtr这是红利随时间变化时股票价格的市场基础成分。需要注意的是,对于上述情形的市场基础成分,需要投资者对于未来红利具有完全预期。当引入预期红利时,上述表达式仍然适用,这时可以修改为: 0)(1j jtjt DErP利用红利预期的股价公式,可以确定价格初值 :0P0jjr如此初值是否满足一般的股价模型,我们可以代入到具有初值的确定解中验证: ttttt DrDP21)()()1(将 代入上式后得到
15、:0j jtjtr这正是在边界条件下所推导的向前预期解,由此可见该解与初值选择是吻合的。例 2.5 我们继续利用滞后算子方法讨论股票价格路径的性质。利用算子表示为: ttDPLr)1(在上述表达式中,滞后算子多项式的特征根小于 1,无法采用逆算子的一般表达式,为此我们需要采取新的定义。时间序列分析方法讲义 第 2 章 滞后算子及其性质7定义滞后算子 的逆算子为 ,具有性质:L1(1) I11(2) )(tty这样一来,滞后算子乘积就具有幂乘的性质:对于任意正整数 和 ,有:ijjijiL对方程(2.12)两端乘以算子多项式: )1(21 )1()()( TrrLr整理得到: TttttT Dr
16、DrrP)1()()1()( 12当 ,且红利序列是有界的,则上述极限为:0r)1(j jtjtDP根据上述运算,可以定义下述算子的逆算子: 2)1()(1)()( LrrLrr2.6 差分方程的求解方法上面我们主要论述了差分方程的表示和外生扰动的动态乘子,下面我们给出差分方程的一般求解过程。第一步:构造 p 阶齐次差分方程,并且寻求齐次方程的 p 个解: ,htiyp,21第二步:构造 p 阶非齐次差分方程的特解。第三步:齐次方程 p 个解的线性组合加上非齐次方程的一个特解,得到非齐次方程的通解。第四步:根据给定的边界条件,确实通解当中的未知参数,得到非齐次方程的确定解。2.6.1 齐次差分
17、方程的通解和稳定性p 阶齐次差分方程的形式是: ptttt yyy21命题 2.3 对于差分方程而言,下述推断成立:(1) 如果 是方程的解,则对任意常数 , 也是解。ht Ahty(2) 如果 和 是方程的解,则对任意实数 和 , 也是方程的解。y1t2 12hthtyA21证明:留做练习。对于 p 阶齐次差分方程,我们尝试地检验解的形式是: ,代入差分方程为:thty0121 pp由此可见, 应该是上述特征方程的根。因此,如果差分方程具有相异实数根的时候,时间序列分析方法讲义 第 2 章 滞后算子及其性质8可以得到 p 个解为: , ,此时解的稳定性要求所有根落在单位圆内。tihtiyp,
18、21命题 2.4 对于齐次差分方程而言:(1) 齐次方程所有特征根落在单位圆内的必要条件是: ;pi1(2) 齐次方程所有特征根落在单位圆内的充分条件是: ;i|(3) 齐次方程至少具有一个单位根的充要条件是: pi1如果齐次方程的特征根出现重根,则应该寻求多项式与指数函数乘机形式的解。例如,如果二阶齐次差分方程具有重根,则两个解应该分别是 , 。thty1tht122.6.2 非齐次差分方程的特解如何寻求非齐次线性差分方程的特解,需要根据非齐次项的具体性质判断。(1) 指数形式的非齐次项此时方程形式是: rttt dbyy10可以尝试特解形式为: ,可以求解出特解为:rtptdc10rpt
19、)/()/(如果 ,尝试解的形式为:1rttdcy如果 ,可以选取其他形式的尝试解。(2) 确定性时间趋势此时方程形式是: dpiitt byy10此时尝试解的形式选为: dpt tctc02.6.3 非齐次差分方程的通解获得齐次方程的 p 个解以后,它们的线性组合构成了齐次方程的通解;如果再获得一个非齐次方程的一个特解,则可以得到非齐次方程的通解,它是齐次方程 p 个解的线性组合加上非齐次方程的特解。这里需要注意的是,上述通解中的线性组合中有 p 个未知系数,这是需要 p 个条件来确定这些系数,这些系数确定以后,就可以得到具有边界或者初值条件的差分方程的具体解,或者唯一解。这个唯一解的系数则具有相应的经济含义。