差分方程模型.

1、差分方程模型一.引言数学模型按照离散的方法和连续的方法，可以分为离散模型和连续模型。1. 确定性连续模型1）微分法建模（静态优化模型），如森林救火模型、血管分支模型、最优价 格模型。2）微分方程建模（动态模型），如传染病模型、人口控制与预测模型、经济 增长模型。3）稳定性方法建模（平衡与稳定状态模型），如军备竞赛模型、种群的互相 竞争模型、种群的互相依存模型、种群弱肉强食模型。4）变分法建模（动态优化模型），如生产计划的制定模型、国民收入的增 长模型、渔业资源的开发模型。2. 确定性离散模型1）逻辑方法建模，如效益的合理分配模型、价格的指数模型。2）层次分析法建模，如旅游景点的选择模型、科研成

2、果的综合评价模型。3）图的方法建模，如循环比赛的名次模型、红绿灯的调节模型、化学制品 的存放模型。4）差分方程建模，如市场经济中的蛛网模型、交通网络控制模型、借贷模 型、养老基金设置模型、人口的预测与控制模型、生物种群的数量模型。随着科学技术的发展，人们将愈来愈多的遇到离散动态系统的问题， 差分方 程就是建立离散动态系统数学模型的有效方法。在一般情况下，动态连续模型用微分方程方法建立， 与此相适应，当时间变 量离散化以后，可以用差分方程建立动态离散模型。 有些实际问题既可以建立连 续模型，又可建立离散模型，究竟采用那种模型应视建模的目的而定。例如，人 口模型既可建立连续模型（其中有马尔萨斯模型

3、 Malthus、洛杰斯蒂克Logistic模 型），又可建立人口差分方程模型。这里讲讲差分方程在建立离散动态系统数学 模型的的具体应用。.差分方程简介在实际中，许多问题所研究的变量都是离散的形式， 所建立的数学模型也是 离散的，譬如，像政治、经济和社会等领域中的实际问题。有些时候，即使所建 立的数学模型是连续形式，例如像常见的微分方程模型、积分方程模型等。但是, 往往都需要用计算机求数值解。这就需要将连续变量在一定的条件下进行离散 化，从而将连续型模型转化为离散型模型。因此，最后都归结为求解离散形式的 差分方程解的问题。关于差分方程理论和求解方法在数学建模和解决实际问题的 过程中起着重要作用

4、。1.差分方程的定义给定一个数列把数列中的前n 1项Xi （i= 0,1,2,，n）关联起来得到的 方程，则称这个方程为差分方程。2.常系数线性齐次差分方程常系数线性齐次差分方程的一般形式为或者表示为Xn ax a2Xnd -akXy = 0，F(n,Xn,Xn i, ,Xn Q =0(1)(T)其中k为差分方程的阶数，其中a-i,a2/ ,ak为差分方程的系数，且ak = 0（k込n）。 对应的代数方程沾+ak+a2扎心+ak = 0（2）称为差分方程（1）的对应的特征方程。（2）式中的k个根、，2,称为（1）式的特 征根2.1差分方程的解常系数线性齐次差分方程的解主要是由相应的特征根的不同

5、情况有不同的形 式。下面分别就特征根为单根、重根和复根的情况给出方程解的形式。2.1.1特征根为单根（互不相同的根）设差分方程（1）有 k个单特征根（互不相同的根）1, 2 / , k，贝UnnnXn* C2 J 亠亠 Ck k为该差分方程 的通解。其中G,C2,，Ck为任意常数，且当给定初始条件Xi =x(0)，(i =1,2,k)时，可以确定一个特解。例1在信道上传输三个字母a,b,c且长度为n的词，规定有两个a连续出现 的词不能传输，试确定这个信道允许传输的词的个数。解：令Xn表示允许传输且长度为为n的词的个数，n =1,2,3/，通过简单计 算可得 xi =3， (a,b,c), x2

6、 =8(即卩 ab,ac, bc, bb,cc,ba,ca,cb)当n _3时，若词的第一个字母是b或c，则词可按Xnd种方式完成；若词的第一个字母是a，则第二个字母是b或c，该词剩下的部分可按Xn种方式完成。于是得差分方程Xn = 2Xn+2Xn_2(n = 3,4,)其特征方程为丸一2扎一2 = 0，特征根为7 r = 1 + 吋 3，入 2 = 1 _ V3则通解为Xn =G(1+V3)n +C2(1V3)n，(n =3,4,)利用条件Xi =3，X2 =8求参数C1， C2，即由3(1 +V3) +c2(1 _ J3) =3Jg(1 +V3)2 +c2(1 J3)2 =8解得2+73-

7、2 + J3G 一厂,c2 -厂2*32J3故得到原差分方程的通解为2，3 n2 亠.3 nXn 二(1、3)n(1- .3)n，(n =1,2,3,4,)2 32.32.1.2特征根为重根设仆2,I是 k阶差分方程xnXn-a2Xn,akxn上=0的I1（1乞丨乞k）个根，重数分别为m1,m2/ ,ml，且mk，则该差分方程的通解 i=1为mim2m|i_!in -i_lin-ii.nxn 八 5nc2in 匕cli n li=1i 4i=1同样的，有给定的初始条件（3）可以唯一确定一个特解例2设初始值为Xo = 1, Xi = 0, X2 = 1, X3 = 2，解差分方程XnXn-3人上

8、-5Xn-2Xn =0 ,（n =4,5,）解：该差分方程的特征方程为432 35，_2二0 ，解得其根为-1,-1,-1,2，故通解为Xn =G(-1)n C2 n(-1)n C3 n2(-1)n C42n代入初始条件X0 = 1, x1 = 0, x2 = 1, x3 = 2，得422910C1 =52C152C3 _ 52，52故该差分方程的满足初始条件的解为525252522.1.3特征根为复根设k阶差分方程Xn a1Xn4 a2 Xn 2 ak Xn止=0的一对共轭复根二和相异的k-2个单根3, 4 k，则该差分方程的通解为xn =C|n cosn 丁c2： nsinn 丁Cs 3c

9、4 ；ck：, B其中_: 2，二-arctan_。ot同样由给定的初始条件(3)可以唯一确定一个特解。另外，对于有多个共轭复根和相异实根，或共轭复根和重根的情况，都可类 似的给出差分方程解的形式。3. 常系数线性非齐次差分方程常系数线性非齐次差分方程的一般形式为Xn aiXn a2Xn,akx.上二 f(n)其中k为差分方程的阶数，其中ai,a2,，ak为差分方程的系数，且ak = O(km n)， f(n)为已知函数。在差分方程中，令f(n) =0,所得方程Xn aiXn a2Xn上亠亠 akXn =0(5)称为非齐次差分方程 对应的齐次差分方程，即与差分方程(1)的形式相同。求解非齐次差

10、分方程通解的一般方法：首先求对应的齐次差分方程(5)的通解X；,然后求非齐次差分方程(4)的一个 特解x；0)，则X x； - x；0)为非齐次差分方程(4)的通解。关于求x；的方法同求差分方程(1)的方法相同。对于求非齐次方程(4)的特解 x；0)的方法，可以用观察法确定，也可以根据f( n)的特性用待定系数法确定，具 体方法可参照常系数线性非齐次微分方程求特解的方法。4. 差分方程的平衡点及其稳定性在应用差分方程研究问题时，一般不需要求出方程的通解，在给定初值后， 通常可用计算机迭代求解，但常常需要讨论解的稳定性。对于差分方程F( n,Xn，Xni,，XnQ =0，若有常数a是其解，即有F

11、( n,a,a, ,a) = 0则称a是差分方程F(n,x；,x；i,，x； Q = 0的平衡点，又对该差分方程的任意由初始条件确定的解Xn =x(n)，均有lim xn = an_.则称这个平衡点a是稳定的；否则是不稳定的。 下面给出一些特殊差分方程的平衡点和稳定性。4.1 一阶常系数线性差分方程一阶常系数线性差分方程的一般形式为Xn 1 aXn 二 b，其中a,b为常数，且a -1,0。它的通解为Xn =C(-a)na +1易知丄是方程(6)的平衡点，由 式知，当且仅当 a +1a1时，一是方程(6)的稳定的平衡点。a +14.2二阶常系数线性差分方程(8)二阶常系数线性差分方程的一般形式

12、为Xn 2 aXn 1 bXn = r，其中a,b,r为常数，当r =0时，它有一特解*X = 0， 当r = 0，且a b 0时，它有一特解* rxa+b+1不管是哪种情形，X*是方程(8)的平衡点。设方程(8)的特征方程为的两个根分别为，二、，=2，贝U 当1, 一是两个不同的实根时，方程(8)的通解为*n-nXn =XC1C1) C2C2)； 当r =，2二，是两个相同实根时，方程(8)的通解为Xn =(G 亠 C2 n) n 当,2= P(cos日+isinT)是一对共轭复根时，方程(8)的通解为Xn = x(Ci cosn ： C2 sinnR易知，当且仅当特征方程的任一特征根 R

13、1时，方程(9)的平衡点是不稳定的。三.差分方程建模实例1.贷款买房问题某居民买房向银行贷款6万元，利息为月利率1%,贷款期为25年，要求建 立数学模型解决如下问题：1) 问该居民每月应定额偿还多少钱？2) 假设此居民每月可节余700元，是否可以去买房？1.1确定参变量：用n表示月份，An表示第n个月欠银行的钱，r表示月利率，x表示每月还钱数，A0表示贷款额。1.2模型的建立与求解1)模型的建立时间欠银行款初始Ao一个月后A = A。(1 r) x二个月后A2 = A (1 r) x三个月后aA = A2 (1 r)xan个月后An 二 An_1(1r) -x由上表可得相邻两个月的递推关系式A

14、n 二 An(1 r) -x1.3模型的求解：(1)差分方程求解方法先求其特解。令An二代二y，则y二y(1 r) - X，得特解为y二上r再求对应齐次方程 人=代丄(1 r)的通解。对应的特征方程为 -(1 r) =0，得=(1 - r)。齐次方程的通解为：c(V r)n因此原方程的通解为：An 二 c(1 r)n r又因为n = 0时=A，得c = A。- Xr故t (1+rJ1-1An = Ao(1 +r f -Xr(2)递推法:An =Ao(1 r)n -x,n1.,n1 r j 亠 亠1 rA0 1 r(1 + r )n -1-xrA0(1+r)nX _ 1 r n 一1,30060

15、000(1+0.01)300.(1+0.01) -1632元A0 =60000, A300 = 0 , n =300, r =0.01001因此，该居民每月应偿还632元。又632 0，加上初始条件x(0) = x。可得未来任意时刻群数量所满足的数学模型为:dxx(O) = Xoxm2.由于是利用种群繁殖周期作为时段来研究种群增长状况，则令:trt视为整数及r(x)=r_ x代入方程(1) 得:Xmrx(t 1)-x(t) =(r )x(t)Xm加上初始条件x(0) =Xo得任意时刻t种群数量所满足的离散型数学模型为x(t 1) = (r 1 - r )x(t)(Xmx(0) = X。通过这个

16、差分方程就可以很容易得到任意时刻 t种群的数量。三.模型求解1.利用Mathematica求解方程(1)，可得任意时刻t种群数量为x(t)二XmMathematic a源程序为:XmX0丿DSolvgx (t) _ r* (1 _xt/xm)*xt = 0,xt,t2.根据方程(2),只要给出初值X。就可以很容易进行递推而得到任意时刻群的数量。四.结果分析1. 上面方程（3）有时称为阻滞增长模型或Logistic模型，它有着广泛的应用。 例如传染病在封闭地区的传播，耐用消费品在有限的市场上的销售等现象，都可 以合理的、简化的用这个模型来进行描述。但它存在不足，因为随着环境的变迁， 最大种群容量

17、可能会发生变化，而且最大种群容量也不容易准确得到。2. 一方面，用离散化的时间来研究问题有时是很方便的，尤其出现了计算机 以后，人们可以很方便的对问题进行求解；另一方面，对这个种群数量问题，由 于许多种群实际上是由单一世代构成的， 在相继的世代之间几乎没有重叠，所以 种群的增长是分步进行的。这种情况下，为了准确的描述种群的数量动态就不能 用微分方程，而应利用离散的模型来描述。4. 人口的控制与预测模型一. 问题的提出常见的两个常微分方程模型（马尔萨斯（Malthus）模型和洛杰斯蒂克（Logistic） 模型）没有考虑到社会成员之间的个体差异，即不同年龄、不同体质的人在死亡、 生育方面存在的差

18、异。完全忽略了这些差异显然是不合理的。但我们不可能对每 一个人的情况逐个加以考虑，故仅考虑年龄的差异对人口的变动的影响，即假设同一年龄的人具有相同的死亡率和生育能力，这样建立的模型不但使我们能够更 细致的预测人口总数，而且能够预测老年人口、劳动力人口、学龄人口等不同年 龄组的人口信息.下面来建立离散的差分数学模型来表现人口数量的变化规律。二. 模型的建立与求解设x/t）为第t年年龄为k的人口数量，k =0,1,2,10 0,即忽略百岁以上的 人口。如果知道了第t年各年龄组的人口数，各年龄组人口的生育及死亡状态， 就可以根据人口发展变化规律推得第t 1年各年龄组的人口数。首先引入k岁人口的死亡率

19、和k岁育龄妇女的年生育率这两个概念，他们的含义和记号如下:k岁人口的年死亡率:dk一年内k岁的死亡人数这年内k岁的人口数k岁妇女的年生育率:bk一年内k岁妇女生育的婴儿数这年内k岁妇女人数第t 1年k 1岁的人口数就是第t年k岁人口数扣除它在该年的死亡人数, 即Xki(t 1) =(1-dk)Xk(t)，令pk =1 - dk称为k岁人口的存活率，故各年龄组人口随时间的变化规律可用递推公式Xk i(t 1) = PkXk(t) , (k =0,1,99)来表示。再考虑到零岁的人数100Xo(t 1)=為 bkUk(t)Xk(t)，k=0其中Uk(t)Xk (t)为第t年k岁的妇女人数，Uk (t)为第t年k岁人口的女性比(占全部k岁人口数)，bkUk(t)Xk(t)就是第t年k岁妇女所生育的婴儿数由此得到的人口模 型是：100Xo(t 1)-為 bkUk(t)Xk(t)“心(1)x“(t +1) = PkXk(t) , k =0,1，,99根据人的生理特征和人口学中的习惯，妇女的育龄区间一般取为15岁至49岁之 间，即当k % = (1)Xo