差分方程模型的理论和方法引言1、差分方程：差分方程反映的是关于离散变量的取值与变化规律。

通过建立一个或几个离散变量取值所满足的平衡关系，从而建立差分方程。

差分方程就是针对要解决的目标，引入系统或过程中的离散变量，根据实际背景的规律、性质、平衡关系，建立离散变量所满足的平衡关系等式，从而建立差分方程。

通过求出和分析方程的解，或者分析得到方程解的特别性质（平衡性、稳定性、渐近性、振动性、周期性等），从而把握这个离散变量的变化过程的规律，进一步再结合其他分析，得到原问题的解。

2、应用：差分方程模型有着广泛的应用。

实际上，连续变量可以用离散变量来近似和逼近，从而微分方程模型就可以近似于某个差分方程模型。

差分方程模型有着非常广泛的实际背景。

在经济金融保险领域、生物种群的数量结构规律分析、疾病和病虫害的控制与防治、遗传规律的研究等许许多多的方面都有着非常重要的作用。

可以这样讲，只要牵涉到关于变量的规律、性质，就可以适当地用差分方程模型来表现与分析求解。

3、差分方程建模：在实际建立差分方程模型时，往往要将变化过程进行划分，划分成若干时段，根据要解决问题的目标，对每个时段引入相应的变量或向量，然后通过适当假设，根据事物系统的实际变化规律和数量相互关系，建立每两个相邻时段或几个相邻时段或者相隔某几个时段的量之间的变化规律和运算关系（即用相应设定的变量进行四则运算或基本初等函数运算或取最运算等）等式（可以多个并且应当充分全面反映所有可能的关系），从而建立起差分方程。

或者对事物系统进行划分，划分成若干子系统，在每个子系统中引入恰当的变量或向量，然后分析建立起子过程间的这种量的关系等式，从而建立起差分方程。

在这里，过程时段或子系统的划分方式是非常非常重要的，应当结合已有的信息和分析条件，从多种可选方式中挑选易于分析、针对性强的划分，同时，对划分后的时段或子过程，引入哪些变量或向量都是至关重要的，要仔细分析、选择，尽量扩大对过程或系统的数量感知范围，包括对已有的、已知的若干量进行结合运算、取最运算等处理方式，目的是建立起简洁、深刻、易于求解分析的差分方程。

在后面我们所举的实际例子中，这方面的内容应当重点体会。

差分方程模型作为一种重要的数学模型，对它的应用也应当遵从一般的数学建模的理论与方法原则。

同时注意与其它数学模型方法结合起来使用，因为一方面建立差分方程模型所用的数量、等式关系的建立都需要其他的数学分析方式来进行；另一方面，由差分方程获得的结果有可以进一步进行优化分析、满意度分析、分类分析、相关分析等等。

第一节 差分方程的基本知识一、 基本概念1、 差分算子设数列{}n x ，定义差分算子n n n x x x -=∆∆+1:为n x 在n 处的向前差分。

而1--=∆n n n x x x 为n x 在n 处的向后差分。

以后我们都是指向前差分。

可见n x ∆是n 的函数。

从而可以进一步定义n x ∆的差分：n n x x 2)(∆=∆∆称之为在n 处的二阶差分，它反映的是的增量的增量。

类似可定义在n 处的k 阶差分为：))((1n k n k x x -∆∆=∆2、 差分算子 、不变算子、平移算子记n n n n x Ix x Ex ==+,1，称E 为平移算子，I 为不变算子 。

则有：n n n n x I E Ix Ex x )(-=-=∆I E -=∆∴由上述关系可得：i n ki i k i k n ik i i k i k n k n k x C x E C x I E x +=-=-∑∑-=-=-=∆00)1()1()( （1）这表明n x 在n 处的k 阶差分由n x 在k n n n ++....1,，处的取值所线性决定。

反之，由 n n n x x x -=∆+1 得 n n n x x x ∆+=+1：n n n n x x x x +-=∆++1222，得：n n n n x x x x 2122∆++-=++,这个关系表明：第n+2项可以用前两项以及相邻三项增量的增量来表现和计算。

即一个数列的任意一项都可以用其前面的k 项和包括这项在内的k+1 项增量的增量的增量……..第k 层增量所构成。

……..,)1(10k n i n k i i k i k n kx x C x ++-=-+-=∆∑得： n k i n k i i k i k k n x x C x ∆+--=+-=-+∑10)1( （2）可以看出：k n x +可以由n k n n x x x ∆∆,...,,的线性组合表示出来3、 差分方程由n x 以及它的差分所构成的方程),...,,,(1n k n n n k x x x n f x -∆∆=∆ （3）称之为k 阶差分方程。

由（1）式可知（3）式可化为),...,,,(11-+++=k n n n k n x x x n F x （4）故（4）也称为k 阶差分方程（反映的是未知数列n x 任意一项与其前，前面k 项之间的关系）。

由（1）和（2）可知，（3）和（4）是等价的。

我们经常用的差分方程的形式是（4）式。

4、 差分方程的解与有关概念（1） 如果n x 使k 阶差分方程（4）对所有的n 成立，则称n x 为方程（4）的解。

（2） 如果-=x x n （-x 为常数）是（4）的解，即),...,,(---=x x n F x则称-=x x n 为（4）的平衡解或叫平衡点。

平衡解可能不只一个。

平衡解的基本意义是：设n x 是（4）的解，考虑n x 的变化性态，其中之一是极限状况，如果x x n n =∞→lim ，则方程（4）两边取极限（x 就存在在这里面），应当有 ),...,,(---=x x n F x（3） 如果（4）的解n x 使得--x x n 既不是最终正的，也不是最终负的，则称n x 为关于平衡点-x 是振动解。

（4） 如果令：--=x x y n n ，则方程（4）会变成),...,,(1-++=k n n k n y y n G y （5）则 0=y 成为（5）的平衡点。

（5） 如果（5）的所有解是关于0=y 振动的，则称k 阶差分方程 （5）是振动方程。

如果（5）的所有解是关于0=y 非振动的，则称k 阶差分方程（5）是非振动方程。

（6） 如果（5）有解n y ，使得对任意大的y N 有 0>≥n N n y Sup y则称n y 为正则解。

（即不会从某项后全为零）（7） 如果方程（4）的解n x 使得-∞→=x x Lim n n ，则称n x 为稳定解。

5、 差分算子的若干性质（1）n n n n y x y x ∆+∆=+∆βαβα)(.)(（2）)(1)(1n n n n n n n n y x x y y y y x ∆-∆=∆+（3）n n n n n n y x x y y x ∆+∆=∆+1)(（4）∑∑==+++∆+-=∆ba k k k ab a k a b b k k y x y x y x x y 111（5）∑=∆=+∆==ni i i n nn n x C x I x E x 0000)(第二节 差分方程常用解法与性质分析1、 常系数线性差分方程的解方程)(...110n b x a x a x a n k k n k n =+++-++ （ 8）其中k a a a ,...,,10为常数，称方程（8）为常系数线性方程。

又称方程0...110=+++-++n k k n k n x a x a x a （9）为方程（8）对应的齐次方程。

如果（9）有形如n n x λ=的解，带入方程中可得：0...1110=++++--k k k k a a a a λλλ （10）称方程（10）为方程（8）、（9）的特征方程。

显然，如果能求出（10）的根，则可以得到（9）的解。

基本结果如下：（1） 若（10）有k 个不同的实根，则（9）有通解：n k k n n n c c c x λλλ+++=...2211，（2） 若（10）有m 重根λ，则通解中有构成项：n m m n c n c c λ)...(121----+++（3）若（10）有一对单复根 βαλi ±=，令：ϕρλi e ±=，αβϕβαρarctan,22=+=，则（9）的通解中有构成项： n c n c n n ϕρϕρsin cos 21--+（4） 若有m 重复根：βαλi ±=，φρλi e ±=，则（9）的通项中有构成项：n n c n c c n n c n c c n m m m m n m m ϕρϕρsin )...(cos )...(1221121---++---+++++++综上所述，由于方程（10）恰有k 个根，从而构成方程（9）的通解中必有k 个独立的任意常数。

通解可记为：-n x如果能得到方程（8）的一个特解：*n x ，则（8）必有通解：=n x -n x +*n x （11）（8） 的特解可通过待定系数法来确定。

例如：如果)(),()(n p n p b n b m m n =为n 的多项式，则当b 不是特征根时，可设成形如)(n q b m n 形式的特解，其中)(n q m 为m 次多项式；如果b是r 重根时，可设特解：r n n b )(n q m ，将其代入（8）中确定出系数即可。

例1 设差分方程1,0,0231012===++++x x x x x n n n ，求n x解：特征方程为0232=++λλ，有根：2,121-=-=λλ故：n n n c c x )2()1(21-+-=为方程的解。

由条件1,010==x x 得：n n n x )2()1(---=2、 二阶线性差分方程组设=)(n z )(n y x n，)(dc b a A =，形成向量方程组 )()1(n Az n z =+ （12）则 )1()1(z A n z n =+ （13）（13）即为（12）的解。

为了具体求出解（13），需要求出n A ，这可以用高等代数的方法计算。

常用的方法有：（1）如果A 为正规矩阵，则A 必可相似于对角矩阵，对角线上的元素就是A 的特征值，相似变换矩阵由A 的特征向量构成：)1()()1(,,111z p p n z p p A p p A n n n Λ=+∴Λ=Λ=---。

（2）将A 分解成ηξξη,,/,=A 为列向量，则有A A n n n .)(.......).(1//.//-===ηξηξηξηξηξ从而，)1(.)()1()1(1/Az z A n z n n -==+ηξ（3） 或者将A 相似于约旦标准形的形式，通过讨论A 的特征值的性态，找出n A 的内在构造规律，进而分析解)(n z 的变化规律，获得它的基本性质。