差分方程模型(讲义).pdf
《差分方程模型(讲义).pdf》由会员分享,可在线阅读,更多相关《差分方程模型(讲义).pdf(22页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、差分方程模型差分方程模型一一.引言引言数学模型按照离散的方法和连续的方法,可以分为离散模型和连续模型。1.1.确定性连续模型确定性连续模型1)微分法建模(静态优化模型),如森林救火模型、血管分支模型、最优价格模型。2)微分方程建模(动态模型),如传染病模型、人口控制与预测模型、经济增长模型。3)稳定性方法建模(平衡与稳定状态模型),如军备竞赛模型、种群的互相竞争模型、种群的互相依存模型、种群弱肉强食模型。4)变分法建模(动态优化模型),如生产计划的制定模型、国民收入的增长模型、渔业资源的开发模型。2.2.确定性离散模型确定性离散模型1)逻辑方法建模,如效益的合理分配模型、价格的指数模型。2)层
2、次分析法建模,如旅游景点的选择模型、科研成果的综合评价模型。3)图的方法建模,如循环比赛的名次模型、红绿灯的调节模型、化学制品的存放模型。4)差分方程建模,如市场经济中的蛛网模型、交通网络控制模型、借贷模型、养老基金设置模型、人口的预测与控制模型、生物种群的数量模型。随着科学技术的发展,人们将愈来愈多的遇到离散动态系统的问题,差分方程就是建立离散动态系统数学模型的有效方法。在一般情况下,动态连续模型用微分方程方法建立,与此相适应,当时间变量离散化以后,可以用差分方程建立动态离散模型。有些实际问题既可以建立连续模型,又可建立离散模型,究竟采用那种模型应视建模的目的而定。例如,人口模型既可建立连续
3、模型(其中有马尔萨斯模型 Malthus、洛杰斯蒂克 Logistic模型),又可建立人口差分方程模型。这里讲讲差分方程在建立离散动态系统数学模型的的具体应用。二二.差分方程简介差分方程简介在实际中,许多问题所研究的变量都是离散的形式,所建立的数学模型也是离散的,譬如,像政治、经济和社会等领域中的实际问题。有些时候,即使所建立的数学模型是连续形式,例如像常见的微分方程模型、积分方程模型等。但是,往往都需要用计算机求数值解。这就需要将连续变量在一定的条件下进行离散化,从而将连续型模型转化为离散型模型。因此,最后都归结为求解离散形式的差分方程解的问题。关于差分方程理论和求解方法在数学建模和解决实际
4、问题的过程中起着重要作用。1.1.差分方程的定义差分方程的定义给定一个数列xn,把数列中的前n1项xi(i 0,1,2,n)关联起来得到的方程,则称这个方程为差分方程。差分方程。2.2.常系数线性齐次差分方程常系数线性齐次差分方程常系数线性齐次差分方程的一般形式为xn a1xn1 a2xn2 akxnk 0,(1)或者表示为F(n,xn,xn1,xnk)0 (1)其中k为差分方程的阶数为差分方程的阶数,其中a1,a2,ak为差分方程的系数,且ak 0(k n)。对应的代数方程k a1k1 a2k2 ak 0 (2)称为差分方程(1)的对应的特征方程特征方程。(2)式中的k个根1,2,k称为(1
5、)式的特征根特征根。2.12.1 差分方程的解差分方程的解常系数线性齐次差分方程的解主要是由相应的特征根的不同情况有不同的形式。下面分别就特征根为单根、重根和复根的情况给出方程解的形式。2.1.12.1.1特征根为单根特征根为单根(互不相同的根互不相同的根)设差分方程(1)有k个单特征根(互不相同的根)1,2,k,则nnnxn c11c22 ckk为该差分方程(1)的通解。其中c1,c2,ck为任意常数,且当给定初始条件xi xi(0),(i 1,2,k)(3)时,可以确定一个特解。例例 1 1 在信道上传输三个字母a,b,c且长度为n的词,规定有两个a连续出现的词不能传输,试确定这个信道允许
6、传输的词的个数。解解:令xn表示允许传输且长度为为n的词的个数,n 1,2,3,,通过简单计算可得x1 3,(a,b,c),x2 8(即 ab,ac,bc,bb,cc,ba,ca,cb)。当n 3时,若词的第一个字母是b或c,则词可按xn1种方式完成;若词的第一个字母是a,则第二个字母是b或c,该词剩下的部分可按xn2种方式完成。于是得差分方程xn 2xn1 2xn2(n 3,4,)其特征方程为2 2 2 0,特征根为113,213则通解为xn c1(13)nc2(13)n,(n 3,4,)利用条件x1 3,x2 8求参数c1,c2,即由c1(13)c2(13)3,22c1(13)c2(13)
7、8解得c1故得到原差分方程的通解为232 3,c2232 3xn232 3(13)n232 3(13)n,(n 1,2,3,4,)2.1.22.1.2特征根为重根特征根为重根设1,2,l是k阶 差 分 方 程xn a1xn1 a2xn2 akxnk 0的l(1 l k)个根,重数分别为m1,m2,ml,且mi k,则该差分方程的通解i1l为xnc1inc2in clini1lni1n1i1n2i1i1i1m1m2ml同样的,有给定的初始条件(3)可以唯一确定一个特解。例例 2 2设初始值为x01,x1 0,x21,x3 2,解差分方程xn xn13xn25xn3 2xn4 0,(n 4,5,)
8、解解:该差分方程的特征方程为43325 2 0,解得其根为1,1,1,2,故通解为xn c1(1)nc2n(1)nc3n2(1)nc42n代入初始条件x01,x1 0,x21,x3 2,得c14229710,c1,c3,c452525252故该差分方程的满足初始条件的解为xn2.1.32.1.3特征根为复根特征根为复根4229710(1)nn(1)nn2(1)n2n52525252设k阶 差 分 方 程xn a1xn1 a2xn2 akxnk 0的 一 对 共 轭 复 根1,2i和相异的k 2个单根3,4,k,则该差分方程的通解为nnxn c1ncosn c2nsin n c33 c4n4 c
9、kk其中22,arctan。同样由给定的初始条件(3)可以唯一确定一个特解。另外,对于有多个共轭复根和相异实根,或共轭复根和重根的情况,都可类似的给出差分方程解的形式。3.3.常系数线性非齐次差分方程常系数线性非齐次差分方程常系数线性非齐次差分方程的一般形式为xn a1xn1 a2xn2 akxnk f(n)(4)其中k为差分方程的阶数为差分方程的阶数,其中a1,a2,ak为差分方程的系数,且ak 0(k n),f(n)为已知函数。在差分方程(4)中,令f(n)0,所得方程xn a1xn1 a2xn2 akxnk 0 (5)称为非齐次差分方程(4)(4)对应的齐次差分方程对应的齐次差分方程,即
10、与差分方程(1)的形式相同。求解非齐次差分方程通解的一般方法:*首先求对应的齐次差分方程首先求对应的齐次差分方程(5)(5)的通解的通解xn,然后求非齐次差分方程,然后求非齐次差分方程(4)(4)的一的一(0)个特解个特解xn,则*(0)xn xn xn为非齐次差分方程(4)的通解。*关于求xn的方法同求差分方程(1)的方法相同。对于求非齐次方程(4)的特(0)解xn的方法,可以用观察法确定,也可以根据f(n)的特性用待定系数法确定,具体方法可参照常系数线性非齐次微分方程求特解的方法。4.4.差分方程的平衡点及其稳定性差分方程的平衡点及其稳定性在应用差分方程研究问题时,一般不需要求出方程的通解
11、,在给定初值后,通常可用计算机迭代求解,但常常需要讨论解的稳定性。对于差分方程F(n,xn,xn1,xnk)0,若有常数a是其解,即有F(n,a,a,a)0则称a是差分方程F(n,xn,xn1,xnk)0的平衡点平衡点,又对该差分方程的任意由初始条件确定的解xn x(n),均有lim x ann则称这个平衡点a是稳定的稳定的;否则是不稳定的不稳定的。下面给出一些特殊差分方程的平衡点和稳定性。4.14.1一阶常系数线性差分方程一阶常系数线性差分方程一阶常系数线性差分方程的一般形式为xn1 axn b,(6)其中a,b为常数,且a 1,0。它的通解为xn C(a)nba 1 (7)易知ba1是方程
12、(6)的平衡点,由(7)式知,当且仅当a 1时,ba1是方程(6)的稳定的平衡点。4.24.2二阶常系数线性差分方程二阶常系数线性差分方程二阶常系数线性差分方程的一般形式为xn2 axn1bxn r,(8)其中a,b,r为常数,当r 0时,它有一特解x*0,当r 0,且a b1 0时,它有一特解x*ra b 1,不管是哪种情形,x*是方程(8)的平衡点。设方程(8)的特征方程为2 ab 0的两个根分别为1,2,则当1,2是两个不同的实根时,方程(8)的通解为xn x*C1(1)n C2(n2);当12是两个相同实根时,方程(8)的通解为xn x*(C1 C2n)n当1,2(cos isin)是
13、一对共轭复根时,方程(8)的通解为xn x*n(C1cosn C2sin n)易知,当且仅当特征方程的任一特征根i1时,平衡点x*是稳定的。4.34.3一阶非线性差分方程一阶非线性差分方程一阶非线性差分方程的一般形式为xn1 f(xn)(9)其平衡点x*由代数方程x f(x)解出。为了分析平衡点x*的稳定性,将方程(9)的右端f(xn)在x*点作泰勒展开,只取一次项,得到xn1 f(x*)(xn x*)f(x*)(10)(10)是(9)的近似线性方程,x*是(10)的平衡点,根据一阶常系数线性差分方程(6)xn1 axn b的稳定性判定的相关结论,得:当f(x*)1时,方程(9)的平衡点是稳定
14、的;当f(x*)1时,方程(9)的平衡点是不稳定的。三三 差分方程建模实例差分方程建模实例1 1 贷款买房问题贷款买房问题某居民买房向银行贷款 6 万元,利息为月利率 1%,贷款期为 25 年,要求建立数学模型解决如下问题:1)问该居民每月应定额偿还多少钱?2)假设此居民每月可节余 700 元,是否可以去买房?1.11.1确定参变量:确定参变量:用n表示月份,An表示第n个月欠银行的钱,r表示月利率,x表示每月还钱数,A0表示贷款额。1.21.2模型的建立与求解模型的建立与求解1)1)模型的建立模型的建立时间初始一个月后二个月后三个月后欠银行款A0A1 A0(1 r)xA2 A1(1 r)xA
15、3 A2(1 r)xAn An1(1 r)xn个月后由上表可得相邻两个月的递推关系式An An1(1 r)x1.31.3模型的求解:模型的求解:(1)差分方程求解方法先求其特解。令An An1 y,则y y(1 r)x,得特解为y 再求对应齐次方程An An1(1 r)的通解。对应的特征方程为x。r(1r)0,得(1r)。齐次方程的通解为:c(1r)n因此原方程的通解为:An c(1r)nxr又因为n 0时An A0,得c A0故xrn1r1An A01r xnr(2)递推法:An A0(1 r)x11 rnn1 r11 r A01 r xn1nr令A0=60000,A300 0,n=300,
16、r=0.01得x A01 rrn1 rn1600001 0.010.013001 0.013001 632元因此,该居民每月应偿还 632 元。又 632700,所以该居民可以去买房。2 2借贷问题借贷问题中国建设银行北京市分行个人住房贷款一至二十年“月均还款金额表”(自1998 年 3 月 25 日起执行)的一部分如下:(借款额为一万元)单位:元贷款期限(年)1520年利率(%)10.20610.206还款总额(元)19569.6023488.80利息负担总和(元)9569.6013488.80月均还款额(元)108.7297.87试问他们是怎样算出来的?借贷问题的数学模型借贷问题的数学模型
17、一一.符号说明符号说明以贷款期限 20 年为例:借贷额-A010,000;贷款期限-为 N 年;月利率-r 10.206/12 0.008505;“月均还款额”-表示每月还款额是相同的,记为x;还款总额-记为S.二二.建立模型建立模型一开始借款A010,000,一个月后欠银行本利为A1 A0(1 r),但为了减少欠款,还了x元,因而A1 A0(1 r)x,第k个月情况也是这样的,即Ak Ak1(1r)x,k 1,2,N注意到了第 N 个月已经不欠银行的钱了,即AN 0,因此,我们得到以下的数学模型:Ak Ak1(1 r)xk 1,2,NA0,x,N KnownFind out such tha
18、t AN 0三三.数学模型的求解数学模型的求解首先求出用已知量表出的表达式。由A2 A1(1 r)x A0(1 r)x(1 r)x A0(1 r)2 x1(1 r)可以猜想,并用数学归纳法证明:Ak A0(1 r)k x1(1 r)(1 r)2 (1 r)k1由等比数列前k 1项的求和公式知:xAk A0(1 r)k(1 r)k1,k 1,2,Nr再由AN 0,得到:A0r(1 r)Nx(1 r)N1把已知量带入,就得到表中的x。3 3生物种群数量问题生物种群数量问题一问题的提出一问题的提出种群的数量问题是当前世界上引起普遍关注的一个问题。要预测未来种群的数量,最重要的影响因素是当前的种群数量
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 方程 模型 讲义
限制150内