第五章 差分方模型.docVIP

第五章 差分方程模型 在第四章中，我们利用微分方程方法研究了一些连续变化的变量。如果将变量离散化，即可得到相应的差分方程模型，为了方便不熟悉差分方程的读者，先对本章用到的差分方程的知识作一简略介绍。 5.1差分方程简介 一、差分方程及其通解 以表示时间，规定只取非负整数。表示第一周期初，表示第二周期初等。记为变量在时刻时的取值，则称为的一阶差分，称为yt的二阶差分。类似地，可以定义yt的阶差分。 由、及的差分给出的方程称为差分方程，其中含的最高阶差分的阶数称为该差分方程的阶。差分方程也可以写成不显含差分的形式。例如，二阶差分方程也可改写成。 满足差分方程的序列称为此差分方程的解。类似于微分方程情况，若解中含有的独立常数的个数等于差分方程的阶数时，则称此解为该差分方程的通解。若解中不含任意常数，则称此解为满足某些初值条件的特解，例如，考察两阶差分方程： ??? 易见与均是它的特解，而则为它的通解，其中，为两个任意常数。 类似于微分方程，称差分方程 （1） 为阶线性差分方程，当时称其为阶非齐次线性差分方程，而 （2） 称为方程（1）对应的齐次线性差分方程。 若（1）中所有的均为与无关的常数，则称其为常系数差分方程，即阶常系数线性差分方程可写成 （3） 其对应的齐次方程为 （4） 容易证明，若序列与均为方程（4）的解，则也是方程（4）的解，其中，为任意常数，这说明，齐次方程的解构成一个线性空间（解空间）。若是方程（4）的解，是方程（3）的解，则也是方程（3）的解。 方程（3）可用如下的代数方法求其通解： （5） ?（步二）根据特征根的不同情况，求齐次方程（4）的通解。 情况1 若特征方程（5）有个互不相同的实根，，…，，则齐次方程（4）的通解为 ????????? (，，…，为任意常数)???情况2 若是特征方程（5）的重根，通解中对应于的项为， （）为任意常数。 情况3 若特征方程（5）有单重复根，通解中对应于的项为，，其中，为的模，为的幅角。 情况4 若为特征方程（）的k重复根，则通解对应于的项为 （）为任意常数。 （步三） 求非齐次方程（）的一个特解。若为方程（4）的通解，则非齐次方程（）的通解为。((((((((3）的特解一般要用到常数变易法，计算较繁。对特殊形式的b(t)也可使用待定系数法。例如，当，为的次多项式时可以证明：若不是特征根，则非齐次方程（）有形如的特解， 也是的次多项式；若是重特征根，则方程（）有形如的特解。进而可利用待定系数法求出，从而得到方程（）的一个特解。例求解两阶差分方程。?解 对应齐次方程的特征方程为，其特征根为，对应齐次方程的通解为 原方程有形如的特解。代入原方程求得，，故原方程的通解为 二、差分方程在应用差分方程研究问题时，一般不需要求出方程的通解，在给定初值后，通常可用计算机迭代求解，但我们常常需要讨论解的稳定性。 （6） 在（6）式中令得到的代数方程的根称为差分方程（6）的平衡点。如果时，则称平衡点是稳定的，否则是不稳定的。由（6）式得: 　 （7） 即方程（6）的解可表为 　　　　　　　　 　　　　　　　　　　（8） 其中由初始值确定。 显然，由（8）式可知差分方程（6）的平衡点稳定的充要条件是 | a | 1 （9） 顺便指出，对于维向量x(k)和常数矩阵A构成的方程组 x( k+1) + Ax( k ) = 0 （10） 其平衡点稳定的条件是A的特征根（）均有 （11） 即均在复平面上的单位圆内。这个结果可由将化为对角阵（或Jordan）阵得到。 2、二阶线性常系数差分方程， （12） 考察方程（12）的平衡点（＝0）的稳定性，设（12）的特征方程的根为，，则不难验证，（12）的通解可表为 （13） 其中常数，