差分方程模型.docx

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

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

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

4、过程中起着重要作用。1.差分方程的定义给定一个数列 把数列中的前n，1项洛(i =0,1,2,，n)关联起来得到的方程，则称这个方程为差分方程。2.常系数线性齐次差分方程常系数线性齐次差分方程的一般形式为Xn aiXn a2Xn, akXn 上=0， (1)或者表示为F(n ,Xn,Xn 1, , Xn k) (1 )其中k为差分方程的阶数，其中aa2,，ak为差分方程的系数，且a0( n)。 对应的代数方程汕+耳汕丄+a2沾二+ak = 0 (2)称为差分方程(1)的对应的特征方程。 式中的k个根1,-2/ ,5称为 式的 特征根。2.1差分方程的解常系数线性齐次差分方程的解主要是由相应的特

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

6、, x2 =8(即 ab,ac, bc, bb,cc,ba,ca,cb) 。当n 一3时，若词的第一个字母是b或c，则词可按禺种方式完成；若词的 第一个字母是a，贝U第二个字母是b或c，该词剩下的部分可按Xn,种方式完成。 于是得差分方程Xn =2xn2xn _2 ( n = 3,4)其特征方程为_ 2，_ 2 = 0，特征根为1=1.3， 2=1-3则通解为Xn =C1(1 、3)n C2(1 -、 3)n，( n =3,4/ )利用条件疋=3， X2 =8求参数C1， C2，即由3(1 +V3)+c2(1 - V3)=3Jw(1 +V3)2 +c2(1 3)2 =8解得2 3 -2 3C1

7、 ， Q2U3 2J3故得到原差分方程的通解为2 - 芒3 n 2 3 nXn (1 3)n (1- .3)n，( n =1,2,34 )2 3 2、32.1.2特征根为重根设仆2,I是 k阶差分方程xn Xn-a2Xn,akxn上=0的I1（1乞丨乞k）个根，重数分别为m1,m2/ ,ml，且mk，则该差分方程的通解 i=1为mi m2 m|i_!in - i_lin - i_li.nxn 八 5n c2in 匕 cli n li=1 i 4 i=1同样的，有给定的初始条件（3）可以唯一确定一个特解。例2 设初始值为xo = 1, xi = 0, X2 = 1, X3 = 2，解差分方程Xn

8、 Xn-3人上-5Xn-2Xn/ =0，（ n =4,5,）解：该差分方程的特征方程为4 3 2 3 5，_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，得42 / l、n 29 / |n + 7 2/ 八 n + 10Xn (T) n(-1) n (-1) 252 52 52 522.1.3特征根为复根故该差分方程的满足初始条件的解为设k阶差分方程Xn Xn4 a2Xn, ak Xn * = 0的一对共轭复根二和相异的k-2个单根3, 4 k，

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

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

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

12、a +1a1时，b是方程(6)的稳定的平衡点。a +14.2二阶常系数线性差分方程二阶常系数线性差分方程的一般形式为Xn 2 aXn 1 bXn 二 r， (8)其中a,b,r为常数，当r =0时，它有一特解*X = 0，当r = 0，且a b 1 = 0时，它有一特解* rxa+b +1不管是哪种情形，x*是方程(8)的平衡点。设方程(8)的特征方程为2 a儿-b = 0的两个根分别为，二、，=2，贝U1当1, 2是两个不同的实根时，方程(8)的通解为Xn =X* G(1)n。2(-2)；2当r =，2二，是两个相同实根时，方程(8)的通解为Xn = (G 亠 C2 n) n3当,2= P(

13、cos日+i sin日)是一对共轭复根时，方程(8)的通解为Xn = x (Ci cosn ： C2 sinnR易知，当且仅当特征方程的任一特征根 R 1时，平衡点x*是稳定的。4.3阶非线性差分方程一阶非线性差分方程的一般形式为Xni 二 f(Xn) (9)其平衡点X由代数方程X = f (x)解出。为了分析平衡点X*的稳定性，将方程(9)的右端f(Xn)在X*点作泰勒展开, 只取一次项，得到 * * *Xn i ： f (X )(Xn - X ) f(X ) (10)(10)是(9)的近似线性方程，X*是(10)的平衡点， 根据一阶常系数线性差分方程xn 1 axb的稳定性判定的相关结论，

14、得：1当f(x*)| 1时，方程(9)的平衡点是稳定的；2当f(x*)| 时，方程(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个月

15、后An = An_1(1 r) X由上表可得相邻两个月的递推关系式An 二 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故代讥1 rn-x丄r递推法:, , n 二 1 . , n1 r j 亠 亠1 r A0 1 rA0 =60000, A300 = o , n =300 , r

16、 =0.01001因此，该居民每月应偿还632元。又632 0 ,加上初始条件x(0) =x0可得未来任意时刻t种群数量所满足的数学模型为：dx ； x =r 1 _ !x* dt i Xm丿 x( 0) = xo2.由于是利用种群繁殖周期作为时 段来研究种群增长状 况，则令r氏=1，t视为整数及r(x) x代入方程 得:Xmrx(t 1)-x(t) =(r )x(t)Xm加上初始条件x(0) =Xo得任意时刻t种群数量所满足的离散型数学模型为x(t 1) = (r 1 - r )x(t) Xmx(0) = Xo通过这个差分方程就可以很容易得到任意时刻 t种群的数量。三.模型求解Xm1 .利用

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

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

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

20、如下: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岁妇女所生育的婴儿数.由此得到的人口 模型是：100X

21、o(t 1)-為 bkUk(t)Xk (t)* 心 (1)和住+1) = PkXk(t) , k=0,1，,99根据人的生理特征和人口学中的习惯，妇女的育龄区间一般取为15岁至49岁之间，即当k 15和k 49时，bk =0， 令x(t) =(X0(t),X1(t), ,Xk(t), ,X100(t)Tu0(t)boui (t)bi比佻 U99 (t )b99ui00 (t)b100Po00 00L =0Pi0 00- 000 p990则人口模型(1)的矩阵形式为x(t 1)=Lx(t) 其中L称为莱斯利(Lwslie )矩阵.当第t0年的人口状况已知时，从式就可以推得第t年的人口为x(t 1

22、) = LtJx(to).5.市场经济中的蛛网模型在自由竞争的市场经济中，商品的价格是由市场上该商品的供应量决定的， 供应量越大，价格就越低。另一方面，生产者提供的商品数量又是由该商品的价 格决定的，价格上升将刺激生产者的生产积极性，导致商品生产量的增加。反之, 价格降低会影响生产者的积极性，导致商品生产量的下降。在没有外界干扰的情 况下，这种现象将如此反复下去。这样的需求和供应关系决定了市场经济中商品 的价格和数量必然是振荡的。这种振荡越小越好，如果振荡太大就会影响人民群 众的正常生活。(1)商品数量与价格的振荡在什么条件下趋向稳定 ？(2)当不稳定时政府能采取什么干预手段使之稳定 ？下面用

23、差分方程理论建模，讨论市场经济趋于稳定的条件，再用图形方法建 立“蛛网模型”对上述现象进行分析，对结果进行解释，然后作适当推广。3.1模型的假设和符号说明1记第n时段商品数量为xn，价格为yn， n= 1,2/ 0这里我们把时间离散化为时段，1个时段相当于商品的1个生产周期，如蔬 菜、水果可以是1年，肉类可以是一个饲养周期。2在门时段商品的价格yn取决于数量Xn。设yn二f (Xn)。它反映消费者对 这种商品的需求关系，称为需求函数。因为商品的数量越多，价格越低。需求函数在图1中用一条下降的曲线f表 示，f称为需求曲线。3在n 1时段商品的数量Xn 1由上一时段的价格yn决定，用Xn1=g(y

24、n) 表示。它反映生产者的供应关系，称为 供应函数。因为价格越高，生产量越大。供应函数在图 1中用一条上升的曲线g表示， g称为供应曲线xo图1 商品供求关系曲线3.2 模型的建立与求解设需求曲线f和供应曲线g相交于点P0(x0, y0)，在F0附近取函数f和g的线性近似，即需求曲线f:yn - y0 -八(xn - X0) , ? 0 (11)供应曲线g :Xn 1-X。(yn-y。), 1 0 (12)由式(11)(12)消去yn，得到一阶线性差分方程Xn .1 - % (1； 一:片)x0， n =1,2, (13)因此Xo是其平衡点，即P0是平衡点。对式(13)进行递推，得Xn 1 -

25、 ) X1 1 ( - ) x。， n - 1,2,由此可得，平衡点稳定的条件是： 厂1 ;不稳定的条件是：门.1。下面用图形解释此模型。若对某一个k有xk =x0，则由(11)式得，当n k时xn =x0，从而yn = y。， 即商品的数量和价格将永远保持在 Po(Xo，yo)点。但是实际生活中的种种干扰使得Xn, yn不可能停止在Po(Xo,y。)上。不妨设X1偏离Xo (见图2，图3)，我们来分 析随着n的增加，Xn,yn的变化情况。图2 Po点是稳定的数量X1给定后，价格y1由曲线f上的P点决定，下一时段的数量X2由曲线g上的P2点决定，这样得到一序列的点R(xi,yj , P2(x2,y2) , P3(xa,ya),P4(X4,y4)，,在图2上，这些点将按照箭头所示方向趋向 Po(xo,yo)，表明Po