《(完整版)差分方程模型(讲义)》由会员分享，可在线阅读，更多相关《(完整版)差分方程模型(讲义)（25页珍藏版）》请在金锄头文库上搜索。
1、差分方程模型一.引言数学模型按照离散的方法和连续的方法，可以分为离散模型和连续模型。1. 确定性连续模型1）微分法建模（静态优化模型），如森林救火模型、血管分支模型、最优价格模型。2）微分方程建模（动态模型），如传染病模型、人口控制与预测模型、经济增长模型。3）稳定性方法建模（平衡与稳定状态模型），如军备竞赛模型、种群的互相竞争模型、种群的互相依存模型、种群弱肉强食模型。4）变分法建模（动态优化模型），如生产计划的制定模型、国民收入的增长模型、渔业资源的开发模型。2. 确定性离散模型1）逻辑方法建模，如效益的合理分配模型、价格的指数模型。2）层次分析法建模，如旅游景点的选择模型、科研成果的综合
2、评价模型。3）图的方法建模，如循环比赛的名次模型、红绿灯的调节模型、化学制品的存放模型。4）差分方程建模，如市场经济中的蛛网模型、交通网络控制模型、借贷模型、养老基金设置模型、人口的预测与控制模型、生物种群的数量模型。随着科学技术的发展，人们将愈来愈多的遇到离散动态系统的问题，差分方程就是建立离散动态系统数学模型的有效方法。在一般情况下，动态连续模型用微分方程方法建立，与此相适应，当时间变量离散化以后，可以用差分方程建立动态离散模型。有些实际问题既可以建立连续模型，又可建立离散模型，究竟采用那种模型应视建模的目的而定。例如，人口模型既可建立连续模型（其中有马尔萨斯模型Malthus、洛杰斯蒂克
3、Logistic模型），又可建立人口差分方程模型。这里讲讲差分方程在建立离散动态系统数学模型的的具体应用。二.差分方程简介在实际中，许多问题所研究的变量都是离散的形式，所建立的数学模型也是离散的，譬如，像政治、经济和社会等领域中的实际问题。有些时候，即使所建立的数学模型是连续形式，例如像常见的微分方程模型、积分方程模型等。但是，往往都需要用计算机求数值解。这就需要将连续变量在一定的条件下进行离散化，从而将连续型模型转化为离散型模型。因此，最后都归结为求解离散形式的差分方程解的问题。关于差分方程理论和求解方法在数学建模和解决实际问题的过程中起着重要作用。1.差分方程的定义给定一个数列x,把数列中
4、的前n+1项x(i=0,1,2,n)关联起来得到的ni方程，则称这个方程为差分方程。2.常系数线性齐次差分方程常系数线性齐次差分方程的一般形式为x+ax+axHFax=0,(1)n1n-12n-2kn-k或者表示为F(n,x,x,x)=0(1nnF1nFk其中k为差分方程的阶数，其中a,a，,a为差分方程的系数，且a丰0(k3时，若词的第一个字母是b或c，则词可按x种方式完成；若词的n-1第一个字母是a，则第二个字母是b或c，该词剩下的部分可按x种方式完成。n-2于是得差分方程x=2x+2x(n=3,4,)nn-1n-2其特征方程为九2-2九-2=0，特征根为九=1+“3，九=1312则通解为
5、x=c(1+3)n+c(13)n，(n=3,4，)n12利用条件x=3，x=8求参数c，c，即由1212c(1+J3)+c(1-V3)=312，c(1+间2+c(173)2=8V12解得2+.32+3c=，c=123223故得到原差分方程的通解为xn=寻(1)n(1-n，(n=1,2,3,4，)2.1.2特征根为重根设九，九，九是k阶差分方程x12l+ax+ax+ax=0的1n-12n-2kn-kl(1lk)个根，重数分别为m,m，,m，12l且m=k，则该差分方程的通解ii=1x=cni-1九n+cni-1九n+cni-1九nn1i12i2lili=1i=1i=1同样的，有给定的初始条件(3
6、)可以唯一确定一个特解。例2设初始值为x二1,x二0,x二1,x二2，解差分方程0123x+x一3x一5x一2x=0,(n=4,5，)nn-1n-2n-3n-4解:该差分方程的特征方程为九4+九3-3九2-5九一2=0,解得其根为-1,-1,-1,2，故通解为x=c(-1)n+cn(-1)n+cn2(-1)n+c2nn12代入初始条件x=1,x=0,x=1,x=2，01234229c=，c=-，152152故该差分方程的满足初始条件的解为4229710x=(-1)nn(-1)n+n2(-1)n+2nn525252522.1.3特征根为复根7c=-35210c=-452设k阶差分方程x+ax+a
7、x+ax=0的一对共轭复根n1n-12n-2kn-k九，九=ai和相异的k-2个单根九，九，九，则该差分方程的通解为1234kx=cpnCOSn0+cpnsinn0+c九n+c九n+c九nn123344kk其中p=Ja2+p2,9=arctan。a同样由给定的初始条件(3)可以唯一确定一个特解。另外，对于有多个共轭复根和相异实根，或共轭复根和重根的情况，都可类似的给出差分方程解的形式。3. 常系数线性非齐次差分方程常系数线性非齐次差分方程的一般形式为x+ax+axHFax=f(n)(4)n1n-12n-2kn-k其中k为差分方程的阶数，其中a,a，,a为差分方程的系数，且a丰0(kn),12k
8、kf(n)为已知函数。在差分方程(4)中，令f(n)=0，所得方程x+ax+axHFax=0(5)n1n-12n-2kn-k称为非齐次差分方程(4)对应的齐次差分方程，即与差分方程(1)的形式相同。求解非齐次差分方程通解的一般方法：首先求对应的齐次差分方程(5)的通解x*，然后求非齐次差分方程的一个n特解x(0)，则nx=x*Fx(0)nnn为非齐次差分方程(4)的通解。关于求x*的方法同求差分方程(1)的方法相同。对于求非齐次方程(4)的特解nx(0)的方法，可以用观察法确定，也可以根据f(n)的特性用待定系数法确定，具n体方法可参照常系数线性非齐次微分方程求特解的方法。4. 差分方程的平衡
9、点及其稳定性在应用差分方程研究问题时，一般不需要求出方程的通解，在给定初值后，通常可用计算机迭代求解，但常常需要讨论解的稳定性。对于差分方程F(n,x,x,x)=0，若有常数a是其解，即有nnF1nFkF(n,a,a,a)=0则称a是差分方程F(n,x,x，,x)=0的平衡点，又对该差分方程的任意由nnF1nFk#初始条件确定的解x=x(n)，均有nlimx=an则称这个平衡点a是稳定的；否则是不稳定的下面给出一些特殊差分方程的平衡点和稳定性4.1一阶常系数线性差分方程一阶常系数线性差分方程的一般形式为x+ax=b,n+1n(6)其中a,b为常数，且aH-1,0。它的通解为=C(-a)n+(7
10、)b易知丄是方程(6)的平衡点，由式知，当且仅当a+1b时，丄是方程(6)的稳定的平衡点。a+14.2二阶常系数线性差分方程二阶常系数线性差分方程的一般形式为x+ax+bx=r，n+2n+1n其中a,b,r为常数，当r=0时，它有一特解(8)当r丰0，且a+b+1丰0时，它有一特解rx*=a+b+1不管是哪种情形，x*是方程(8)的平衡点。设方程(8)的特征方程为的两个根分别为九=九，九=九，则当九，九是两个不同的实根时，方程(8)的通解为x=x*+C(九)n+C(九)n；n11221212# 当X=X=X是两个相同实根时，方程(8)的通解为12x=x*+(C+Cn)Xnn12 当X二p(co
11、s0+isin0)是一对共轭复根时，方程的通解为1,2x=x*+pn(Ccosn0+Csinn0)n12易知，当且仅当特征方程的任一特征根|XJ1时，平衡点x*是稳定的。4.3一阶非线性差分方程一阶非线性差分方程的一般形式为x二f(x)(9)n+1n其平衡点x*由代数方程x=f(x)解出为了分析平衡点x*的稳定性，将方程(9)的右端f(x)在x*点作泰勒展开，n只取一次项，得到x-f(x*)(x一x*)+f(x*)(10)n+1n(10)是(9)的近似线性方程，x*是(10)的平衡点，根据一阶常系数线性差分方程(6)x+ax=b的稳定性判定的相关结论，得：n+1n 当|f(x*)|1时，方程(
12、9)的平衡点是不稳定的。三差分方程建模实例1贷款买房问题某居民买房向银行贷款6万元，利息为月利率1%，贷款期为25年，要求建立数学模型解决如下问题：1) 问该居民每月应定额偿还多少钱？2) 假设此居民每月可节余700元，是否可以去买房?11确定参变量：用n表示月份，A表示第n个月欠银行的钱，r表示月n利率，X表示每月还钱数，A表示贷款额。01.2模型的建立与求解1)模型的建立时间欠银彳丁款初始A0一个月后A=A(1+r)-x10二个月后A=A(1+r)-x21三个月后A=A(1+r)-x32n个月后A=A(1+r)-xnn-1由上表可得相邻两个月的递推关系式A=A(1+r)一xnn-11.3模型的求解：(1)差分方程求解方法先求其特解。令A=A=y，则y=y(1+r)-x，得特解为y=-nn一1r再求对应齐次方程A二A(1+r)的通解。对应的特征方程为nn-1九-(1+r)=0，得X=(1+r)。齐次方程的通解为：c(1+r)n因此原方程的通解为：xA=c(1+r)n+nrx又因为n=0时A=A，得c=A一一n00r故(2)递推法：A=A(1+r)n一x1+U+rn0A=60000,A=0,n=3
