數(shù)學(xué)建模之微分方程建模與平衡點(diǎn)理論

微分方程列微分方程常用的方法：(1)根據(jù)規(guī)律列方程利用數(shù)學(xué)、力學(xué)、物理、化學(xué)等學(xué)科中的定理或經(jīng)過(guò)實(shí)驗(yàn)檢驗(yàn)的規(guī)律來(lái)建立微分方程模型。(2)微元分析法利用已知的定理與規(guī)律尋找微元之間的關(guān)系式，與第一種方法不同的是對(duì)微元而不是直接對(duì)函數(shù)及其導(dǎo)數(shù)應(yīng)用規(guī)律。(3)模擬近似法在生物、經(jīng)濟(jì)等學(xué)科的實(shí)際問(wèn)題中，許多現(xiàn)象的規(guī)律性不很清楚，即使有所了解也是極其復(fù)雜的，建模時(shí)在不同的假設(shè)下去模擬實(shí)際的現(xiàn)象，建立能近似反映問(wèn)題的微分方程，然后從數(shù)學(xué)上求解或分析所建方程及其解的性質(zhì)，再去同實(shí)際情況對(duì)比,檢驗(yàn)此模型能否刻畫、模擬某些實(shí)際現(xiàn)象.、模型的建立與求解1。1傳染病模型(1)基礎(chǔ)模型假設(shè)：t時(shí)刻病人人數(shù)%(/)連續(xù)可微。每天每個(gè)病人有效接觸(使病人治病的接觸)的人數(shù)為九，t=0時(shí)有%0個(gè)病人.建模：t到t+At病人人數(shù)增加TOC\o"1-5"\h\z%(t+At)一％(t)=九％(t)At (1)\o"CurrentDocument"—二九％,%(0)=% (2)\o"CurrentDocument"dt 0解得：%(t)-%e兀 (3)所以，病人人數(shù)會(huì)隨著t的增加而無(wú)限增長(zhǎng),結(jié)論不符合實(shí)際。(2)SI模型

假設(shè)：1。疾病傳播時(shí)期，總?cè)藬?shù)N保持不變.人群分為兩類，健康者占總?cè)藬?shù)的比例為s(t),病人占總?cè)藬?shù)的比例為i(t)。2。每位病人每天平均有效接觸九人，九為日接觸率.有效接觸后健康者變?yōu)椴∪恕Ｒ罁?jù)：患病人數(shù)的變化率=Ni(t)(原患病人數(shù))*九s(t)(每個(gè)病人每天使健康人變?yōu)椴∪说娜藬?shù))建模:Nd-=XNsi

dt由于設(shè)t=0時(shí)刻病人所占的比例為io,則可建立Logistic模型d八i(j)i(0)=,0(4)(6)解得:(4)(6)?i圖形如下，凰2?i圖形如下，凰2軟模犁巾:;-■■曲戲結(jié)論：在不考慮治愈情況下①當(dāng)i=1①當(dāng)i=1時(shí)d達(dá)到最大值h、dt2dt這時(shí)tm(1 、H1ln—-1u0 1②t―8時(shí)人類全被感染.未考慮治愈情況。SIS模型假設(shè)：1。疾病傳播時(shí)期，總?cè)藬?shù)N保持不變.人群分為兩類，健康者占總?cè)藬?shù)的比例為s（t），病人占總?cè)藬?shù)的比例為i（t）。.每位病人每天平均有效接觸九人,九為日接觸率。有效接觸后健康者變?yōu)椴∪?.在所有病人中，每天有比例R的人能被治愈，治愈后看作可被感染的健康者，傳染病的平均傳染期為-。R依據(jù)：患病人數(shù)的變化率二九Nsi（患病人數(shù)的變化率）一RNi（治愈率）建模:(8)(9)Nd二九Nsi-RNi

dt(8)(9)did二九3）一同"（0）30令。為整個(gè)傳染期內(nèi)每位病人有效接觸的平均人數(shù)，0=儲(chǔ)3.則有di——=一九ii—dt(10)用Matlab繪制出色?idi——=一九ii—dt(10)用Matlab繪制出色?i（圖3，圖5）和

dti~t(圖4,圖6）。N中鹿城是小的情況圖石STS模咆的曲輕圖5號(hào)塔值里的:rj曲然fY”結(jié)論：O=1為一個(gè)閾值。①。>1,i(t)極限值代)=1-1為增函數(shù)，i(t)的增減性由i的大小確定。O 0②。41,病人比例i(t)越來(lái)越小，最終趨于0。(4)SIR模型(某些疾病患者治愈后獲得了很強(qiáng)的免疫力，不會(huì)再次被感染)假設(shè):①總?cè)藬?shù)N不變，將人群分為健康者，病人，和病愈免疫的移除者，他們?cè)诳側(cè)藬?shù)中所占的比例依次為s(t),i(t)，r(t)。②九為病人的日接觸率，〃為日治愈率，O=入"為傳染期接觸數(shù).建模：由假設(shè)1得(11)(12)Ndr=從Ni

dt(12)令t=0時(shí)健康者與病人所占比例分別為s0(s0〉0),i0(i0>0)，則有五…i，五…i，d二一九si-,Idti(0)=i0s(0)=s0(13)利用Matlab繪制出i(t)，s(t)(圖7)，?(圖8)圖形,i?s圖形稱為相軌線。圖8L占圖?(相軌線)圖圖8L占圖?(相軌線)圖73。 40相軌線分析:利用相軌線討論解i(t),s(t)的性質(zhì)。^-i平面稱為相平面，相軌線在其上的定義域?yàn)?s,i)eD為(14)s>0,i>0,s+i<1}(14)消去方程中的dt，并由。得到(15)解得:(15)解得:(16)在定義域D內(nèi)，相軌線是上式所表示的曲線，如圖9所示，其中箭頭表示隨著時(shí)間t的增加s(t)和i(t)的變化趨勢(shì)。下面分析s(t)、i(t)和r(t)的變化情況(tf8時(shí)它們的極限值分別記做汽，二和二)圖95JR模型的相軌威①不論初始條件s0,i0如何，病人最終會(huì)消失，i8-0，證明：首先，由式(13)，4<0，而s(t)>0，所以s存在；由式(11)，蟲>0,dt d dt而r(t)<1，所以［存在；由式(11)得18存在。其次，若L―〉0,則由式(11)，對(duì)于充分大的t有dr>嗚，導(dǎo)致丫:8與7存在相矛盾。從圖形來(lái)看，無(wú)論相軌線從何點(diǎn)出發(fā)，最終都將與S軸相交.②令式(16)中i=0,則最終未被感染的健康者的比例是Sy…8為方程, 1Fs八s+i—s+—ln—y—0 (17)0 0 8oS0在(0,1/o)內(nèi)的根，在圖形上表示為相軌線與s軸在(0,1/o)內(nèi)交點(diǎn)的橫坐標(biāo)。③若S0>1/°，則i")先增加，當(dāng)s—1/o時(shí)，i(t)達(dá)到最大值i—s+i--(1+lnos) (18)8 0 0o 0然后i(t)減小旦趨于0，s(t)單調(diào)減小至“，如圖中由P1出發(fā)的相軌線。④若s0Ao，則i(t)單調(diào)減小至0，s(t)單調(diào)減小至S8，如圖中由P2出發(fā)的相軌線。結(jié)論:①若病人比例有一段時(shí)間增長(zhǎng)即認(rèn)為傳染病在蔓延，則1/。為一個(gè)閾值,s0>1/o時(shí)蔓延。可以通過(guò)減小o使s0Ao，使傳染病不蔓延。②s0>1/o,o減小時(shí)，s8增加，也能控制蔓延程度。1。2捕魚模型考察一個(gè)漁場(chǎng)，其中魚量在天然環(huán)境下按一定規(guī)律增長(zhǎng)、如果捕撈量恰好等于增長(zhǎng)量，那么漁場(chǎng)魚量將保持不變，這個(gè)捕撈量就可以持續(xù).①產(chǎn)量模型假設(shè)：X(t)為漁場(chǎng)中魚量.1。無(wú)捕撈時(shí)，魚的的增長(zhǎng)服從logistic規(guī)律，即X(t)=f(x)=rx(1-金］ (19)IN)其中：r表示固有增長(zhǎng)率，N表示環(huán)境容許的最大魚量，f(x)表示單位時(shí)間的增長(zhǎng)量。2。用E表示單位時(shí)間捕撈率，單位時(shí)間捕撈量和漁場(chǎng)魚量X(t)成正比，則有單位時(shí)間捕撈量為h(x)=Ex (2。)建模：捕撈情況下漁場(chǎng)魚量滿足,、_.、乙x、—x(t)=F(x)=rx1一n-Ex (21)其中：F(x)=f(x)一h(x)。判斷x(t)的穩(wěn)定條件，求式(21)的平衡點(diǎn)，分析其穩(wěn)定性。令式(21)為0,得兩個(gè)平衡點(diǎn)：Ex=N(1 ),x=0 (22)0 r1穩(wěn)定性判斷Ff(x)=E—r,Ff(x)=r—E當(dāng)E<r時(shí)F(x0)<0,F(xi)>0，則x0點(diǎn)穩(wěn)定，匕點(diǎn)不穩(wěn)定。當(dāng)E>r時(shí)F(x0)>0,F(xi)<0，則\點(diǎn)穩(wěn)定，x0點(diǎn)不穩(wěn)定。分析：用E表示捕撈率,丫表示固有增長(zhǎng)率.①當(dāng)E<r時(shí)，可使魚量穩(wěn)定在x0，獲得穩(wěn)定產(chǎn)量.②當(dāng)E>r時(shí)，xi穩(wěn)定，漁場(chǎng)干枯。根據(jù)(19)，(20)式分別繪制曲線尸f(x)及y=h(x)=E(x)，使用Matlab繪制圖形如下所示，得兩曲線交點(diǎn)為P，則P橫坐標(biāo)為穩(wěn)定平衡點(diǎn)x0，縱坐標(biāo)為穩(wěn)定條件下單位時(shí)間的產(chǎn)量，當(dāng)交點(diǎn)位于拋物線頂點(diǎn)時(shí)獲得最大的持續(xù)產(chǎn)量，此時(shí)的穩(wěn)定平衡點(diǎn)為x*=N,單位時(shí)間的最大持續(xù)產(chǎn)量為h=曳，捕撈率e*=r。0 2 m4 2結(jié)論：將捕撈率控制在固有增長(zhǎng)率r的一半，即使?jié)O場(chǎng)魚量保持在最大魚量的一半時(shí)，能夠獲得最大的持續(xù)產(chǎn)量。②效益模型(經(jīng)濟(jì)效益=總收入收入一成本)假設(shè)：魚銷售單價(jià)P，單位捕撈率費(fèi)用是。，單位時(shí)間收入為丁，成本為S，單位利潤(rùn)為R，則有T=ph(x)=pEx S=cER=T—S=pEx—cE (23)建模:在穩(wěn)定條件x=)下，將式(22)代入式(23)得TOC\o"1-5"\h\zR(E)=T(E)-S(E)=pNE(1-—)-cE (24)r求出使利潤(rùn)最大的捕撈強(qiáng)度為E=-[1-—] (25)R21pN)最大利潤(rùn)下的漁場(chǎng)穩(wěn)定魚量、和單位時(shí)間的持續(xù)產(chǎn)量hRNcxR=2*2P (26)h=rx(1-鼠)=且[1-。] (27)RRN4Ip2N2J結(jié)論：當(dāng)有最大效益時(shí)，捕撈率和持續(xù)產(chǎn)量都減小，漁場(chǎng)應(yīng)保持的穩(wěn)定魚量增加，捕撈成本越大或銷售價(jià)格越低所需減少增大的部分越大.③捕撈過(guò)度：封閉式捕撈追求利益最大，開放式捕撈只追求利潤(rùn)。令式(24)中R(E)=0，解Ej則E=/1-—] (28)SIPNJ當(dāng)E<Es時(shí)，利潤(rùn)r(e)〉0經(jīng)營(yíng)者加大捕撈強(qiáng)度,當(dāng)E>E/r(e)<0經(jīng)營(yíng)者減小捕撈強(qiáng)度，E為盲目捕撈下的臨界強(qiáng)度.S或利用Matlab繪制—?T(E),S(E)曲線如圖(12)，則t(E),S(E)交點(diǎn)橫坐標(biāo)即為ES。二、微分方程與平衡點(diǎn)理論2。1一階微分方程設(shè)一階微分方程為x(t)=f(x) (1)求解方程f(X)=0即可出平衡點(diǎn)x=x0。再判斷平衡點(diǎn)x0是否穩(wěn)定。判斷平衡點(diǎn)的常用方法有以下兩種(1)直接法將f(x)在x0點(diǎn)作泰勒展開，僅取一次項(xiàng)，則得方程(1)的近似線性方程為

xG）=/,（x）（x-x） （2）0所以，X也是方程（2）的平衡點(diǎn).令廣（x）=〃,則方程（2）的一般解為0 0xG）=+x c為常數(shù)o對(duì)于x點(diǎn)的穩(wěn)定性有如下結(jié)論:0TOC\o"1-5"\h\z如果廣（x）＜0,則X對(duì)于方程（2）和（1）都是穩(wěn)定的；0 0如果尸G）＞0,則X對(duì)于方程（2）和（1）都是不穩(wěn)定的；0 0（2）間接法如果存在x某個(gè)鄰域內(nèi)的任意值，使方程（1）的解式方）滿足0limxG）=x （3）…0那么X是穩(wěn)定的，否則X是不穩(wěn)定的.0 02.2二階微分方程設(shè)二階微分方程為(4)f（xx）—0求出方程/J2 的解,即為二階微分方程的平衡點(diǎn)X=X0,X=X。記作TOC\o"1-5"\h\zgkx,x）=0 1122l12PGo,xo）0 1 2利用直接法判斷平衡點(diǎn)的穩(wěn)定性，由線性常系數(shù)微分方程組xax+ax＜1zX11 22 （5）x\t）=bx+bxl2 11 22得系數(shù)矩陣記(6)(7)(6)(7)(8)A=i2bbL1 2」為求出方程（5）的惟一平衡點(diǎn)P（0,0）的穩(wěn)定性，令A(yù)的行列式為odetA^OP（0,0）的穩(wěn)定性可由方程（5）的特征方程的根九決定。即odet（A-1/）=0方程（8）可以寫為入2+p入+q=0TOC\o"1-5"\h