时间序列分析讲义__第01章_差分方程

时间序列分析方法讲义 第 1 章 差分方程1第一章 差分方程差分方程是连续时间情形下微分方程的特例。差分方程及其求解是时间序列方法的基础，也是分析时间序列动态属性的基本方法。经济时间序列或者金融时间序列方法主要处理具有随机项的差分方程的求解问题，因此，确定性差分方程理论是我们首先需要了解的重要内容。1.1 一阶差分方程假设利用变量 表示随着时间变量 变化的某种事件的属性或者结构，则 便是在时tyt ty间 可以观测到的数据。假设 受到前期取值 和其他外生变量 的影响，并满足下述t t 1tytw方程：(1.1)ttt wy10在上述方程当中，由于 仅线性地依赖前一个时间间隔自身的取值 ，因此称具有ty 1ty这种结构的方程为一阶线性差分方程。如果变量 是确定性变量，则此方程是确定性差分tw方程；如果变量 是随机变量，则此方程是随机差分方程。在下面的分析中，我们假设t是确定性变量。tw例 1.1 货币需求函数 假设实际货币余额、实际收入、银行储蓄利率和商业票据利率的对数变量分别表示为 、 、 和 ，则可以估计出美国货币需求函数为：tmtIbtrctctt r019.45.19.0720上述方程便是关于 的一阶线性差分方程。可以通过此方程的求解和结构分析，判断t其他外生变量变化对货币需求的动态影响。1.1.1 差分方程求解：递归替代法差分方程求解就是将方程变量表示为外生变量及其初值的函数形式，可以通过以前的数据计算出方程变量的当前值。由于方程结构对于每一个时间点都是成立的，因此可以将(1.1)表示为多个方程：0t 010wyy：1：t ttt10依次进行叠代可以得到： 1012101011 )()( wywyy 0231202)(1.2)itittit yy0110上述表达式(1.2)便是差分方程(1.1)的解，可以通过代入方程进行验证。上述通过叠代将 表示为前期变量和初始值的形式，从中可以看出 对这些变量取值的依赖性和动态变t ty化过程。1.1.2. 差分方程的动态分析：动态乘子(dynamic multiplier)在差分方程的解当中，可以分析外生变量，例如 的变化对 阶段以后的 的影响。0wtty假设初始值 和 不受到影响，则有：1ytw,时间序列分析方法讲义 第 1 章 差分方程2(1.3)ttwy10类似地，可以在解的表达式中进行计算，得到：(1.4)jtj1上述乘子仅仅依赖参数 和时间间隔 ，并不依赖观测值的具体时间阶段，这一点在1j任何差分方程中都是适用的。例 1.2 货币需求的收入乘子 在我们获得的货币需求函数当中，可以计算当期收入一个单位的变化，对两个阶段以后货币需求的影响，即： tttt IwIwmI2122利用差分方程解的具体系数，可以得到：，19.0tI7.从而可以得到二阶乘子为： 8.2tIm注意到上述变量均是对数形式，因此实际上货币需求相对于两个阶段以前收入的弹性系数，这意味着收入增长 1%，将会导致两个阶段以后货币需求增加 0.098%，其弹性是比较微弱的。定义 1.1 在一阶线性差分方程中，下述乘子系列称为 相对于外生扰动 的反应函tytw数：， (1.5)jtjjwyL1,0显然上述反应函数是一个几何级数，其收敛性依赖于参数 的取值。1(1) 当 时，反应函数是单调收敛的；01(2) 当 时，反应函数是震荡收敛的；(3) 当 时，反应函数是单调扩张的；(4) 当 时，反应函数是震荡扩张的；1可以归纳描述反应函数对于参数的依赖性：当 时，反应函数是收敛的；当1|时，反应函数是发散的。|1一个特殊情形是 的情形，这时扰动将形成持续的单一影响，即 的一个单位变1 tw化将导致其后任何时间 的一个单位变化：jty，tjjwyL,0为了分析乘子的持久作用，假设时间序列 的现值贴现系数为 ，则未来所有时间的ty流贴现到现在的总值为：ty(1.6)0jjt如果 发生一个单位的变化，而 不变，那么所产生的对于上述贴现量的影响twtsw,为边际导数：时间序列分析方法讲义 第 1 章 差分方程3， 000 1/)( jjjtj tjt wyy |上述分析的是外生变量的暂时扰动，如果 发生一个单位的变化，而且其后的t也都发生一个单位的变化，这意味着变化是持久的。这时持久扰动对于 时tsw, )(jt刻的 的影响乘数是：jy(1.7)0111 jjjttt wy当 时，对上式取极限，并将其识为扰动所产生的持久影响：|1(1.8)11)(limjtttjj yy例 1.3 货币需求的长期收入弹性 在例 1.1 中我们已经获得了货币的短期需求函数，从中可以求出货币需求的长期收入弹性为： 68.072.19tttdIwdI这说明收入增加 1%最终将导致货币需求增加 0.68%，这是收入对于货币需求反馈的持久影响效果。如果换一个角度考察扰动的影响，那么我们需要分析一个单位的外生扰动对于 以后ty路径的累积影响，这时可以将这种累积影响表示为：(1.9)10jtjwy由此可见，如果能够估计出差分方程中的系数，并且了解差分方程解的结构，则可以对经济变量进行稳定性的动态分析。另外，我们也发现，内生变量对外生变量反应函数的性质比较敏感地依赖差分方程中的系数。1.2 阶差分方程p如果在方程当中允许 依赖它的 阶前期值和输入变量，则可以得到下述 阶线性差typ p分方程( 将常数项归纳到外生变量当中 )：(1.10)ttttt wy21为了方便起见，将上述差分方程表示成为矩阵形式：(1.11)tttvF其中：， ，121pttttyy 0101321 pttv其实在方程(1.11)所表示的方程系统当中，只有第一个方程是差分方程 (1.10)，而其余方程均是定义方程：，jtjtyp,2将 阶差分方程表示成为矩阵形式的好处在于，它可以进行比较方便的叠代处理，同p时可以更方便地进行稳定性分析。另外，差分方程的系数都体现在矩阵 F 的第一行上。时间序列分析方法讲义 第 1 章 差分方程4进行向前叠代，可以得到差分方程的矩阵解为：(1.12)tttttt vFvF 1101利用 表示矩阵 中第 行、第 列元素，则方程系统(1.12) 中的第一个方程可以表)(jifij示为： ttttptttt wfwffyfyffy 1)(1)(0)(1)1(2)1()1( (1.13)需要注意，在 阶差分方程的解中需要知道 个初值： ，以及从时p ,2pyy刻 开始时的所有外生变量的当前和历史数据： 。0 ),(10t由于差分方程的解具有时间上的平移性，因此可以将上述方程(1.12)表示为：(1.14)jtjttjtjtjjt vFvF11类似地，表示成为单方程形式： (1.15)jtjttjtj pjjt wfwff yyyy 1)(1)()(12 利用上述表达式，可以得到 阶差分方程的动态反应乘子为：p，)(1jtjjfyL,0由此可见，动态反应乘子主要由矩阵 的首个元素确定。jF例 1.4 在 阶差分方程中，可以得到一次乘子为：p1)1(1FfwyLt二次乘子为： 212)(121 ft虽然可以进一步通过叠代的方法求出更高阶的反应乘子，但是利用矩阵特征根表示则更为方便，主要能够更为方便地求出矩阵 的首个位置的元素。jF根据定义，矩阵 的特征根是满足下述的 值：F(1.16)0|pI一般情况下，可以根据行列式的性质，将行列式方程转换为代数方程。例 1.5 在二阶差分方程当中，特征方程为： 01)( 212具体可以求解出两个特征根为：， (1.17)21142214上述特征根的表达式在讨论二阶线性差分方程解的稳定性时，我们还要反复用到。距阵 的特征根与 阶差分方程表达式之间的联系可以由下述命题给出：Fp命题 1.1 距阵 的特征根满足下述方程，此方程也称为 阶线性差分方程的特征方程：p0121 ppp 证明：根据特征根的定义，可知特征根满足：时间序列分析方法讲义 第 1 章 差分方程50101)(| 132 ppIF对上述行列式进行初等变化，将第 列乘以 加到第 列，然后将第 列乘p)/(1p1p以 加到第 列，依次类推，可以将上述行列式方程变化为对角方程，并求出行列)/1(2p式值为： 011 pp 这便是所求的 阶线性差分方程的特征方程。 END如果知道 阶线性差分方程的特征方程及其特征根，不仅可以分析差分方程的动态反应乘子，而且可以求解出差分方程解析解的动态形式。1.2.1 具有相异特征根的 阶线性差分方程的通解p根据线性代数的有关定理，如果一个方阵具有相异特征根，则存在非奇异矩阵 将其T化为对角矩阵，且对角线元素便是特征根：， (1.18)1TF),(1pdiag这时矩阵 的乘级或者幂方矩阵可以简单地表示为：， )(Tjjj ),(1jpjjdiag(1.19)假设变量 和 分别表示矩阵 和 的第 行、第 列元素，则可以将上述方程利ijtij ij用矩阵形式表示为： pppjpjjpppj ttttttttF 21211221212120从中可以获得：(1.19)jpjjjpj ccttt 21 121)(1 )()()其中： ， ，如此定义的序列具有下述约束条件(自行证明) ：jjt,0(1.20)21pc具有上述表达式以后，在差分方程的解：(1.15)jtjttjtj pjjt wfwff yyyy 1)(1)()(12 中可以得到动态乘子为：， jpjjjtjj ccfyL21)(1 ,0(1.21)究竟系数序列 取值如何，下述命题给出了它的具体表达式。jc命题 1.2 如果矩阵 的特征根是相异的，则系数 可以表示为：Fjc时间序列分析方法讲义 第 1 章 差分方程6(1.22)pikkiic,1)(证明：由于假设矩阵 具有相异的特征根，因此对角化的非奇异矩阵可以由特征向量F构造。令向量 为：it，1,21ipit pi,2其中 是矩阵 的第 个特征根。i经过运算可以得到： iitF由此可知 是矩阵 的对应特征根 的特征向量，利用每个 做列就可以得到矩阵 。FiitT将矩阵 的第一列表示出来：pIT1 001111,3223121 pppttt可以求解上述线性方程的解为：，)()(13121 pt 2t )()(11211 pppt 注意到： ， ，带入上述表达式即可得到结论。 iitc,END例 1.6 求解二阶差分方程： tttt wyy21.06.解：该方程的特征方程为： 02.6.2特征根为：，84.0)2.().(.124.0).(4)6.0(.212，780(21c )12c此方程的动态乘子为：，jjjjtjjwyL )24.0(.84.0(21 ,1在上述乘子的作用过程中，绝对值教大的特征根决定了乘子的收敛或者发散过程。一般情形下，如果 是绝对值最大的特征根，则有：1(1.23)(limcyjtj 时间序列分析方法讲义 第 1 章 差分方程7则动态乘子的收敛或者发散是以指数速度进行。当一些特征根出现复数的时候，差分方程解的性质出现了新的变化，扰动反应函数将出现一定的周期性质。为此，我们讨论二阶差分方程的情形。当 时，特征方程具有共扼复根，可以表示为：0421， ，biabia，/ 2/121)4)(/利用复数的三角函数或者指数表示法，可以将其写作：， ，expsinco1 iRR2baab/tn这时动态乘子可以表示为： )si()()cos)( 212121 jRcjcwyL jjjjtjj 对于实系统的扰动分析，上述反应乘子应该是实数。由于 和 也是共扼复数，因此有：，ic1ic2则有：(1.24)sin(2)os(jRjRwyLjjtjj 如果 ，即复数处于单位圆上，则上述动态乘子出现周期性变化，并且影响不会消1失；如果 ，即复数处于单位圆内，则上述动态乘子按照周期方式进行率减，其作用慢慢消失；如果 ，即复数处于单位圆外，则上述动态乘子按照周期方式进行扩散，其作用将逐渐增强。例 1.7 求解二阶差分方程： tttt wyy218.05.解：该方程的特征方程为： 08.5.2特征根为：，i86.025.).(4).(.1280502上述共扼复数的模为： 9.)6.().(2R因为 ，由此可知其动态乘子呈现收敛趋势。可以具体计算出其震荡的周期模式。1，8.0/cosa.1由此可知动态乘子的周期为： 9.42由此可知动态乘子的时间轨迹上，大于 4.9 个时间阶段便出现一次高峰。1.2.2 具有相异特征根的二阶线性差分方程的通解针对具体的二阶线性差分方程，可以讨论解的性质与参数 之间的关系。21,a. 当 时，参数取值处于抛物线 的下方。这时特征方程具有复特0421 214征根，且复数的模为： 221212 /)4()/(baR时间序列分析方法讲义 第 1 章 差分方程8因此，当 时，此时解系统是震荡收敛的；当 是震荡维持的；当102 2时是震荡发散的。12b. 当特征根为实数时，我们分析最大特征根和最小特征的性质。此时 ，0421且 212211 44 121当且仅当 时解及其动态反应乘子是稳定的。下面我们判断非稳定情形。如果： 142211即： 121求解可知，使得不等式 成立的参数解为：，或者，同理，使得不等式 成立的参数解为：2，或者，211因此当特征方程具有相异实根的时候，稳定性要求参数落入抛物线上的三角形区域内。c. 类似地可以说明，当特征方程具有相等实根的时候，即处于三角形内的抛物线上时，方程仍然具有稳定解，同时动态反应乘子也是收敛的。1.2.3 具有重复特征根的 阶线性差分方程的通解p在更为一般的情形下，矩阵 可能具有重复的特征根，即具有重根。此时可以利用FJordan 标准型表示差分方程的解及其动态反应乘子。下面以二阶差分方程为例说明。假设二阶差分方程具有重根，则可以将矩阵 表示为：F10TF计算矩阵乘积得到： 1jj于是动态反应乘子可以表示为： jjjtjj kfwyL21)(1.3 长期和现值的计算如果矩阵 的所有特征根均落在单位圆内(即所有特征根的模小于 1)，当时间间隔 逐F j渐增大时，矩阵乘积 将趋于零矩阵。如果外生变量 和 的数据均是有界的，则可以j twty利用 的所有历史数据表示差分方程的一个解：t 321tttt wwy其中 ，即矩阵 中的(1, 1) 位置元素。可以在矩阵表示下，计算 的一个暂)(iif