第十六章差分方程模型教材.doc

第十六章 差分方程模型 离散状态转移模型涉及的范围很广，可以用到各种不同的数学工具。下面我们对差分方程作一简单的介绍，下一章我们将介绍马氏链模型。 §1 差分方程 1.1 差分方程简介 规定只取非负整数。记为变量在点的取值，则称为的一阶向前差分，简称差分，称为的二阶差分。类似地，可以定义的阶差分。 由及的差分给出的方程称为的差分方程，其中含的最高阶差分的阶数称为该差分方程的阶。差分方程也可以写成不显含差分的形式。例如，二阶差分方程也可改写成。 满足一差分方程的序列称为差分方程的解。类似于微分方程情况，若解中含有的独立常数的个数等于差分方程的阶数时，称此解为该差分方程的通解。若解中不含任意常数，则称此解为满足某些初值条件的特解。 称如下形式的差分方程 （1） 为阶常系数线性差分方程，其中是常数，。其对应的齐次方程为 （2） 容易证明，若序列与均为（2）的解，则也是方程（2）的解，其中为任意常数。若是方程（2）的解，是方程（1）的解，则也是方程（1）的解。 方程（1）可用如下的代数方法求其通解： （I）先求解对应的特征方程 （3） （II）根据特征根的不同情况，求齐次方程（2）的通解。 （i）若特征方程（3）有个互不相同的实根，则齐次方程（2）的通解为 （为任意常数） （ii）若是特征方程（3）的重根，通解中对应于的项为，为任意常数。 （iii）若特征方程（3）有单重复根，通解中对应它们的项为，其中为的模，为的幅角。 （iv）若是特征方程（3）的重复根，则通解对应于它们的项为 为任意常数。 （III）求非齐次方程（1）的一个特解。若为方程（2）的通解，则非齐次方程（1）的通解为。 求非齐次方程（1）的特解一般要用到常数变易法，计算较繁。对特殊形式的也可使用待定系数法。例如，当，为的次多项式时可以证明：若不是特征根，则非齐次方程（1）有形如的特解，也是的次多项式；若是重特征根，则方程（1）有形如的特解。进而可利用待定系数法求出，从而得到方程（1）的一个特解。 例1 求解两阶差分方程。 解 对应齐次方程的特征方程为，其特征根为，对应齐次方程的通解为 原方程有形如的特解。代入原方程求得，，故原方程的通解为 例2 在信道上传输仅用三个字母且长度为的词，规定有两个连续出现的词不能传输，试确定这个信道容许传输的词的个数。 解 令表示容许传输且长度为的词的个数，，通过简单计算可求得：，。当时，若词的第一个字母是或则词可按种方式完成；若词的第一个字母是，则第二个字母是或，该词剩下的部分可按种方式完成。于是，得差分方程 ， 其特征方程为 特征根 ， 则通解为 ， 利用条件，，求得 ， 在应用差分方程研究问题时，我们常常需要讨论解的稳定性。对常系数非齐次线性差分方程（1），若不论其对应齐次方程的通解中任意常数如何取值，在时总有，则称方程（1）的解是稳定的。根据通解的结构不难看出，非齐次方程（1）稳定的充要条件为其所有特征根的模均小于1。 1.2 常系数线性差分方程的变换解法 常系数线性差分方程采用解析解法比较容易，而且对其解的意义也容易理解，但采用这种解法求解常系数线性非齐次差分方程比较繁琐，通常是采用变换，将差分方程变换为代数方程去求解。 设有离散序列，，则的变换定义为 （4） 其中是复变量。显然上式右端的级数收敛域是某个圆的外部。 的反变换记作 1.2.1 几个常用离散函数的变换 （i）单位冲激函数的变换 即单位冲激函数的变换为1。 （ii）单位阶跃函数的变换 ， 即 （iii）单边指数函数的的变换（为不等于1的正常数） 1.2.2 变换的性质 （i）线性性质 设，，则 其中为常数。收敛域为和的公共区域。 （ii）平移性 设，则 ， ， ， 例3 求齐次差分方程 ，， 的解。 解 令，对差分方程取变换，得 ， ， 对上式取反变换，便得差分方程的解为 。 §2 蛛网模型 2.1 问题提出 在自由竞争的社会中，很多领域会出现循环波动的现象。在经济领域中，可以从自由集市上某种商品的价格变化看到如下现象：在某一时期，商品的上市量大于需求，引起价格下跌，生产者觉得该商品无利可图，转而经营其它商