差分方程模型一. 引言数学模型按照离散的方法和连续的方法，可以分为离散模型和连续模型。
1. 确定性连续模型1) 微分法建模(静态优化模型)，如森林救火模型、血管分支模型、最优价格模型。
2) 微分方程建模(动态模型)，如传染病模型、人口控制与预测模型、经济增长模型。
3) 稳定性方法建模(平衡与稳定状态模型)，如军备竞赛模型、种群的互相竞争模型、种群的互相依存模型、种群弱肉强食模型。
4) 变分法建模（动态优化模型），如生产计划的制定模型、国民收入的增长模型、渔业资源的开发模型。
2. 确定性离散模型1) 逻辑方法建模，如效益的合理分配模型、价格的指数模型。
2) 层次分析法建模，如旅游景点的选择模型、科研成果的综合评价模型。
3）图的方法建模，如循环比赛的名次模型、红绿灯的调节模型、化学制品的存放模型。
4）差分方程建模，如市场经济中的蛛网模型、交通网络控制模型、借贷模型、养老基金设置模型、人口的预测与控制模型、生物种群的数量模型。
随着科学技术的发展，人们将愈来愈多的遇到离散动态系统的问题，差分方程就是建立离散动态系统数学模型的有效方法。
在一般情况下，动态连续模型用微分方程方法建立，与此相适应，当时间变量离散化以后，可以用差分方程建立动态离散模型。
有些实际问题既可以建立连续模型，又可建立离散模型，究竟采用那种模型应视建模的目的而定。
例如，人口模型既可建立连续模型(其中有马尔萨斯模型Malthus、洛杰斯蒂克Logistic模型)，又可建立人口差分方程模型。
这里讲讲差分方程在建立离散动态系统数学模型的的具体应用。
二. 差分方程简介在实际中，许多问题所研究的变量都是离散的形式，所建立的数学模型也是离散的，譬如，像政治、经济和社会等领域中的实际问题。
有些时候，即使所建立的数学模型是连续形式，例如像常见的微分方程模型、积分方程模型等。
但是，往往都需要用计算机求数值解。
这就需要将连续变量在一定的条件下进行离散化，从而将连续型模型转化为离散型模型。
因此，最后都归结为求解离散形式的差分方程解的问题。
关于差分方程理论和求解方法在数学建模和解决实际问题的过程中起着重要作用。
1. 差分方程的定义给定一个数列{}n x , 把数列中的前1+n 项i x ),,2,1,0(n i =关联起来得到的方程，则称这个方程为差分方程。
2. 常系数线性齐次差分方程常系数线性齐次差分方程的一般形式为02211=++++---k n k n n n x a x a x a x ， (1)或者表示为0),,,,(1=++k n n n x x x n F (1’)其中k 为差分方程的阶数，其中k a a a ,,,21 为差分方程的系数，且0≠k a )(n k ≤。
对应的代数方程02211=++++--k k k k a a a λλλ (2)称为差分方程(1)的对应的特征方程。
(2)式中的k 个根k λλλ,,,21 称为(1)式的特征根。
2.1 差分方程的解常系数线性齐次差分方程的解主要是由相应的特征根的不同情况有不同的形式。
下面分别就特征根为单根、重根和复根的情况给出方程解的形式。
2.1.1 特征根为单根(互不相同的根)设差分方程(1)有k 个单特征根(互不相同的根)k λλλ,,,21 ，则nk k n n n c c c x λλλ+++= 2211为该差分方程(1)的通解。
其中k c c c ,,,21 为任意常数，且当给定初始条件)0(i i x x =，),,2,1(k i = (3)时，可以确定一个特解。
例1 在信道上传输三个字母c b a ,,且长度为n 的词， 规定有两个a 连续出现的词不能传输，试确定这个信道允许传输的词的个数。
解: 令n x 表示允许传输且长度为为n 的词的个数， ,3,2,1=n ，通过简单计算可得 31=x ，(a,b,c), 82=x (即ab,ac, bc, bb,cc,ba,ca,cb)。
当3≥n 时，若词的第一个字母是b 或c ，则词可按1-n x 种方式完成; 若词的第一个字母是a ，则第二个字母是b 或c ，该词剩下的部分可按2-n x 种方式完成。
于是得差分方程2122--+=n n n x x x ( ,4,3=n )其特征方程为0222=--λλ，特征根为311+=λ， 312-=λ则通解为n n n c c x )31()31(21-++=， ( ,4,3=n )利用条件31=x ，82=x 求参数1c ，2c ，即由⎪⎩⎪⎨⎧=-++=-++8)31()31(3)31()31(222121c c c c ， 解得32321+=c ， 32322+-=c故得到原差分方程的通解为n n n x )31(3232)31(3232-+-+++=， ( ,4,3,2,1=n )2.1.2 特征根为重根设l λλλ ,,21是k 阶差分方程02211=++++---k n k n n n x a x a x a x 的l )1(k l ≤≤个根，重数分别为l m m m ,,,21 ，且k m li i =∑=1，则该差分方程的通解为n l i m i li n i m i i n i m i i n n c n c n c x lλλλ112112111121-=-=-=∑∑∑+++=同样的，有给定的初始条件(3)可以唯一确定一个特解。
例2 设初始值为2,1,0,13210====x x x x ，解差分方程02534321=---+----n n n n n x x x x x ， ( ,5,4=n )解: 该差分方程的特征方程为0253234=---+λλλλ，解得其根为2,1,1,1---，故通解为n n n n n c n c n c c x 2)1()1()1(42321+-+-+-=代入初始条件2,1,0,13210====x x x x ，得52421=c ，52291-=c ，5273=c ，52104=c 故该差分方程的满足初始条件的解为n n n n n n n x 25210)1(527)1(5229)1(52422+-+---=2.1.3 特征根为复根设k 阶差分方程02211=++++---k n k n n n x a x a x a x 的一对共轭复根βαλλi ±=21,和相异的2-k 个单根k λλλ ,,43，则该差分方程的通解为nk k n n n n n c c c n c n c x λλλθρθρ+++++= 443321sin cos其中22βαρ+=，αβθarctan=。
同样由给定的初始条件(3)可以唯一确定一个特解。
另外，对于有多个共轭复根和相异实根，或共轭复根和重根的情况，都可类似的给出差分方程解的形式。
3. 常系数线性非齐次差分方程常系数线性非齐次差分方程的一般形式为)(2211n f x a x a x a x k n k n n n =++++--- (4)其中k 为差分方程的阶数，其中k a a a ,,,21 为差分方程的系数，且0≠k a )(n k ≤，)(n f 为已知函数。
在差分方程(4)中，令0)(=n f ，所得方程02211=++++---k n k n n n x a x a x a x (5)称为非齐次差分方程(4)对应的齐次差分方程，即与差分方程(1)的形式相同。
求解非齐次差分方程通解的一般方法：首先求对应的齐次差分方程(5)的通解*n x ，然后求非齐次差分方程(4)的一个特解)0(n x ，则)0(*n n n x x x +=为非齐次差分方程(4)的通解。
关于求*n x 的方法同求差分方程(1)的方法相同。
对于求非齐次方程(4)的特解)0(n x 的方法，可以用观察法确定，也可以根据)(n f 的特性用待定系数法确定，具体方法可参照常系数线性非齐次微分方程求特解的方法。
4. 差分方程的平衡点及其稳定性在应用差分方程研究问题时，一般不需要求出方程的通解，在给定初值后，通常可用计算机迭代求解，但常常需要讨论解的稳定性。
对于差分方程0),,,,(1=++k n n n x x x n F ，若有常数a 是其解，即有0),,,,(=a a a n F则称a 是差分方程0),,,,(1=++k n n n x x x n F 的平衡点，又对该差分方程的任意由初始条件确定的解)(n x x n =，均有a x n n =∞→lim则称这个平衡点a 是稳定的；否则是不稳定的。
下面给出一些特殊差分方程的平衡点和稳定性。
4.1 一阶常系数线性差分方程一阶常系数线性差分方程的一般形式为b ax x n n =++1， (6)其中b a ,为常数，且0,1-≠a 。
它的通解为1)(++-=a ba C x n n (7) 易知1+a b是方程(6)的平衡点，由(7)式知，当且仅当 1<a时，1+a b是方程(6)的稳定的平衡点。
4.2 二阶常系数线性差分方程二阶常系数线性差分方程的一般形式为r bx ax x n n n =++++12， (8)其中r b a ,,为常数，当0=r 时，它有一特解0*=x ，当0≠r ，且01≠++b a 时，它有一特解1*++=b a rx ，不管是哪种情形，*x 是方程(8)的平衡点。
设方程(8)的特征方程为02=++b a λλ的两个根分别为1λλ=，2λλ=，则① 当21,λλ是两个不同的实根时，方程(8)的通解为n n n C C x x )()(2211*λλ++=；② 当λλλ==21是两个相同实根时，方程(8)的通解为n n n C C x x λ)(21*++=③ 当)sin (cos 2,1θθρλi +=是一对共轭复根时，方程(8)的通解为)sin cos (21*θθρn C n C x x n n ++=易知，当且仅当特征方程的任一特征根1<i λ时，平衡点*x 是稳定的。
4.3 一阶非线性差分方程一阶非线性差分方程的一般形式为)(1n n x f x =+ (9)其平衡点*x 由代数方程)(x f x =解出。
为了分析平衡点*x 的稳定性，将方程(9)的右端)(n x f 在*x 点作泰勒展开，只取一次项，得到)())((***'1x f x x x f x n n +-≈+ (10)(10)是(9)的近似线性方程，*x 是(10)的平衡点， 根据一阶常系数线性差分方程(6)b ax x n n =++1的稳定性判定的相关结论，得：① 当1)(*'<x f 时，方程(9)的平衡点是稳定的； ② 当1)(*'>x f 时，方程(9)的平衡点是不稳定的。