第五章 差分方程模型.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）的通解，则非齐次方程（）的通解为。求非齐次方程（）的特解一般要用到常数变易法，计算较繁。对特殊形式的b(t)也可使用待定系数法。例如，当，为的次多项式时可以证明：若不是特征根，则非齐次方程（）有形如的特解，也是的次多项式；若是重特征根，则方程（）有形如的特解。进而可利用待定系数法求出，从而得到方程（）的一个特解。例求解两阶差分方程。?解 对应齐次方程的特征方程为，其特征根为，对应齐次方程的通解为原方程有形如的特解。代入原方程求得故原方程的通解为差分方程在应用差分方程研究问题时，一般不需要求出方程的通解，在给定初值后，通常可用计算机迭代求解，但我们常常需要讨论解的稳定性。非齐次线性差分方程线性方程 (16) 方程（16）的平衡点由代数方程解出。为分析的稳定性，将方程（16）的右端在点作Taylor展开，只取一次项，（16）近似为 (17) （17）是（16）的近似线性方程，也是（17）的平衡点。关于线性方程（17）的稳定平衡点的讨论已由（6）～（9）给出，而当时方程（16）与（17）平衡点的稳定性相同。于是得到当 （18） 时，非线性方程（16）的平衡点是稳定的；当 （19） 时，非线性方程（16）的平衡点是不稳定的。 5.2市场经济中的蛛网模型 时段的上市数量为，价格为。这里我们把时间离散化为时段，1个时段相当于商品的1个生产周期，如蔬菜、水果是一个种植周期，肉类是牲畜的饲养周期。 同一时段商品的价格取决于数量，设 （1） 它反映消费者对这种商品的需求关系，称需求函数。因为商品的数量越多价格越低，所以在图1中用一条下降曲线表示它，称需求函数。 下一时段商品的数量