差分方程模型(講義).

1、差分方程模型一.引言數(shù)學模型按照離散的方法和連續(xù)的方法，可以分為離散模型和連續(xù)模型。1. 確定性連續(xù)模型1）微分法建模（靜態(tài)優(yōu)化模型），如森林救火模型、血管分支模型、最優(yōu)價 格模型。2）微分方程建模（動態(tài)模型），如傳染病模型、人口控制與預測模型、經(jīng)濟 增長模型。3）穩(wěn)定性方法建模（平衡與穩(wěn)定狀態(tài)模型），如軍備競賽模型、種群的互相 競爭模型、種群的互相依存模型、種群弱肉強食模型。4）變分法建模（動態(tài)優(yōu)化模型），如生產(chǎn)計劃的制定模型、國民收入的增 長模型、漁業(yè)資源的開發(fā)模型。2. 確定性離散模型1）邏輯方法建模，如效益的合理分配模型、價格的指數(shù)模型。2）層次分析法建模，如旅游景點的選擇模型、科研成

2、果的綜合評價模型。3）圖的方法建模，如循環(huán)比賽的名次模型、紅綠燈的調(diào)節(jié)模型、化學制品 的存放模型。4）差分方程建模，如市場經(jīng)濟中的蛛網(wǎng)模型、交通網(wǎng)絡控制模型、借貸模 型、養(yǎng)老基金設置模型、人口的預測與控制模型、生物種群的數(shù)量模型。隨著科學技術的發(fā)展，人們將愈來愈多的遇到離散動態(tài)系統(tǒng)的問題， 差分方 程就是建立離散動態(tài)系統(tǒng)數(shù)學模型的有效方法。在一般情況下，動態(tài)連續(xù)模型用微分方程方法建立， 與此相適應，當時間變 量離散化以后，可以用差分方程建立動態(tài)離散模型。 有些實際問題既可以建立連 續(xù)模型，又可建立離散模型，究竟采用那種模型應視建模的目的而定。例如，人 口模型既可建立連續(xù)模型（其中有馬爾薩斯模型

3、 Malthus、洛杰斯蒂克Logistic模 型），又可建立人口差分方程模型。這里講講差分方程在建立離散動態(tài)系統(tǒng)數(shù)學 模型的的具體應用。.差分方程簡介在實際中，許多問題所研究的變量都是離散的形式， 所建立的數(shù)學模型也是 離散的，譬如，像政治、經(jīng)濟和社會等領域中的實際問題。有些時候，即使所建 立的數(shù)學模型是連續(xù)形式，例如像常見的微分方程模型、積分方程模型等。但是, 往往都需要用計算機求數(shù)值解。這就需要將連續(xù)變量在一定的條件下進行離散 化，從而將連續(xù)型模型轉化為離散型模型。因此，最后都歸結為求解離散形式的 差分方程解的問題。關于差分方程理論和求解方法在數(shù)學建模和解決實際問題的 過程中起著重要作用

4、。1.差分方程的定義給定一個數(shù)列把數(shù)列中的前n 1項Xi （i= 0,1,2,，n）關聯(lián)起來得到的方程，則稱這個方程為差分方程。2.常系數(shù)線性齊次差分方程常系數(shù)線性齊次差分方程的一般形式為或者表示為Xn ax a2Xnd -akXy = 0，F(xiàn)(n,Xn,Xn i, ,Xn Q =0(1)(T)其中k為差分方程的階數(shù)，其中a-i,a2/ ,ak為差分方程的系數(shù)，且ak = 0（k込n）。 對應的代數(shù)方程沾+ak+a2扎心+ak = 0（2）稱為差分方程（1）的對應的特征方程。（2）式中的k個根、，2,稱為（1）式的特 征根 2.1差分方程的解常系數(shù)線性齊次差分方程的解主要是由相應的特征根的不同

5、情況有不同的形 式。下面分別就特征根為單根、重根和復根的情況給出方程解的形式。2.1.1特征根為單根（互不相同的根）設差分方程（1）有 k個單特征根（互不相同的根）1, 2 / , k，貝UnnnXn* C2 J 亠亠 Ck k為該差分方程 的通解。其中G,C2,，Ck為任意常數(shù)，且當給定初始條件Xi =x(0)，(i =1,2,k)時，可以確定一個特解。例1在信道上傳輸三個字母a,b,c且長度為n的詞，規(guī)定有兩個a連續(xù)出現(xiàn) 的詞不能傳輸，試確定這個信道允許傳輸?shù)脑~的個數(shù)。解：令Xn表示允許傳輸且長度為為n的詞的個數(shù)，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求參數(shù)C1， C2，即由3(1 +V3) +c2(1 _ J3) =3Jg(1 +V3)2 +c2(1 J3)2 =8解得2+73-2

7、 + 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）個根，重數(shù)分別為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 3c4

9、 ；ck：, B其中_: 2，二-arctan_。ot同樣由給定的初始條件(3)可以唯一確定一個特解。另外，對于有多個共軛復根和相異實根，或共軛復根和重根的情況，都可類 似的給出差分方程解的形式。3. 常系數(shù)線性非齊次差分方程常系數(shù)線性非齊次差分方程的一般形式為Xn aiXn a2Xn,akx.上二 f(n)其中k為差分方程的階數(shù)，其中ai,a2,，ak為差分方程的系數(shù)，且ak = O(km n)， f(n)為已知函數(shù)。在差分方程中，令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)的方法，可以用觀察法確定，也可以根據(jù)f( n)的特性用待定系數(shù)法確定，具 體方法可參照常系數(shù)線性非齊次微分方程求特解的方法。4. 差分方程的平衡點及其穩(wěn)定性在應用差分方程研究問題時，一般不需要求出方程的通解，在給定初值后， 通常可用計算機迭代求解，但常常需要討論解的穩(wěn)定性。對于差分方程F( n,Xn，Xni,，XnQ =0，若有常數(shù)a是其解，即有F(

11、 n,a,a, ,a) = 0則稱a是差分方程F(n,x；,x；i,，x； Q = 0的平衡點，又對該差分方程的任意由初始條件確定的解Xn =x(n)，均有l(wèi)im xn = an_.則稱這個平衡點a是穩(wěn)定的；否則是不穩(wěn)定的。 下面給出一些特殊差分方程的平衡點和穩(wěn)定性。4.1 一階常系數(shù)線性差分方程一階常系數(shù)線性差分方程的一般形式為Xn 1 aXn 二 b，其中a,b為常數(shù)，且a -1,0。它的通解為Xn =C(-a)na +1易知丄是方程(6)的平衡點，由 式知，當且僅當 a +1a1時，一是方程(6)的穩(wěn)定的平衡點。a +14.2二階常系數(shù)線性差分方程(8)二階常系數(shù)線性差分方程的一般形式為

12、Xn 2 aXn 1 bXn = r，其中a,b,r為常數(shù)，當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 1

13、時，方程(9)的平衡點是不穩(wěn)定的。三.差分方程建模實例1.貸款買房問題某居民買房向銀行貸款6萬元，利息為月利率1%,貸款期為25年，要求建 立數(shù)學模型解決如下問題：1) 問該居民每月應定額償還多少錢？2) 假設此居民每月可節(jié)余700元，是否可以去買房？1.1確定參變量：用n表示月份，An表示第n個月欠銀行的錢，r表示月 利率，x表示每月還錢數(shù)，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?？傻梦磥砣我鈺r刻群數(shù)量所滿足的數(shù)學模型為:dxxmx(O) = Xo2.由于是利用種群繁殖周期作為時段來研究種群增長狀況，則令:trt視為整數(shù)及r(x)=r_ x代入方程(1) 得:Xmrx(t 1)-x(t) =(r )x(t)Xm加上初始條件x(0) =Xo得任意時刻t種群數(shù)量所滿足的離散型數(shù)學模型為x(t 1) = (r 1 - r )x(t)(Xmx(0) = X。通過這個

16、差分方程就可以很容易得到任意時刻 t種群的數(shù)量。三.模型求解1.利用Mathematica求解方程(1)，可得任意時刻t種群數(shù)量為x(t)二XmXmX0丿Mathematic a源程序為:DSolvgx (t) _ r* (1 _xt/xm)*xt = 0,xt,t2.根據(jù)方程(2),只要給出初值X。就可以很容易進行遞推而得到任意時刻群的數(shù)量。四.結果分析281. 上面方程（3）有時稱為阻滯增長模型或Logistic模型，它有著廣泛的應用。 例如傳染病在封閉地區(qū)的傳播，耐用消費品在有限的市場上的銷售等現(xiàn)象，都可 以合理的、簡化的用這個模型來進行描述。但它存在不足，因為隨著環(huán)境的變遷， 最大種群

17、容量可能會發(fā)生變化，而且最大種群容量也不容易準確得到。2. 一方面，用離散化的時間來研究問題有時是很方便的，尤其出現(xiàn)了計算機 以后，人們可以很方便的對問題進行求解；另一方面，對這個種群數(shù)量問題，由 于許多種群實際上是由單一世代構成的， 在相繼的世代之間幾乎沒有重疊，所以 種群的增長是分步進行的。這種情況下，為了準確的描述種群的數(shù)量動態(tài)就不能 用微分方程，而應利用離散的模型來描述。4. 人口的控制與預測模型一. 問題的提出常見的兩個常微分方程模型（馬爾薩斯（Malthus）模型和洛杰斯蒂克（Logistic） 模型）沒有考慮到社會成員之間的個體差異，即不同年齡、不同體質(zhì)的人在死亡、 生育方面存在

18、的差異。完全忽略了這些差異顯然是不合理的。但我們不可能對每 一個人的情況逐個加以考慮，故僅考慮年齡的差異對人口的變動的影響，即假設同一年齡的人具有相同的死亡率和生育能力，這樣建立的模型不但使我們能夠更 細致的預測人口總數(shù)，而且能夠預測老年人口、勞動力人口、學齡人口等不同年 齡組的人口信息.下面來建立離散的差分數(shù)學模型來表現(xiàn)人口數(shù)量的變化規(guī)律。二. 模型的建立與求解設x/t）為第t年年齡為k的人口數(shù)量，k =0,1,2,10 0,即忽略百歲以上的 人口。如果知道了第t年各年齡組的人口數(shù)，各年齡組人口的生育及死亡狀態(tài)， 就可以根據(jù)人口發(fā)展變化規(guī)律推得第t 1年各年齡組的人口數(shù)。首先引入k歲人口的死

19、亡率和k歲育齡婦女的年生育率這兩個概念，他們的含義和記號如下:k歲人口的年死亡率:dk一年內(nèi)k歲的死亡人數(shù)這年內(nèi)k歲的人口數(shù)k歲婦女的年生育率:bk一年內(nèi)k歲婦女生育的嬰兒數(shù)這年內(nèi)k歲婦女人數(shù)第t 1年k 1歲的人口數(shù)就是第t年k歲人口數(shù)扣除它在該年的死亡人數(shù), 即Xki(t 1) =(1-dk)Xk(t)，令pk =1 - dk稱為k歲人口的存活率，故各年齡組人口隨時間的變化規(guī)律可用遞推公式Xk i(t 1) = PkXk(t) , (k =0,1,99)來表示。再考慮到零歲的人數(shù)100Xo(t 1)=為 bkUk(t)Xk(t)，k=0其中Uk(t)Xk (t)為第t年k歲的婦女人數(shù)，Uk (t)為第t年k歲人口的女性比(占全部k歲人口數(shù))，bkUk(t)Xk(t)就是第t年k歲婦女所生育的嬰兒數(shù)由此得到的人口模型是：100Xo(t 1)-為 bkUk(t)Xk(t)“心(1)x“(t +1) = PkXk(t) , k =0,1，,99根據(jù)人的生理特征和人口學中的習慣，婦女的育齡區(qū)間一般取為15歲至49歲之 間，即當k % = (1)Xo