西安電子科技大學(xué)數(shù)學(xué)建模講義第三講.ppt
《西安電子科技大學(xué)數(shù)學(xué)建模講義第三講.ppt》由會員分享，可在線閱讀，更多相關(guān)《西安電子科技大學(xué)數(shù)學(xué)建模講義第三講.ppt（75頁珍藏版）》請?jiān)谘b配圖網(wǎng)上搜索。
數(shù)學(xué)建模講義,主講人:穆學(xué)文,西安電子科技大學(xué)數(shù)學(xué)系Email:mxw1334@,第三講微分方程模型,,,動態(tài)模型,描述對象特征隨時間(空間)的演變過程,分析對象特征的變化規(guī)律,預(yù)報對象特征的未來性態(tài),研究控制對象特征的手段,根據(jù)函數(shù)及其變化率之間的關(guān)系確定函數(shù),微分方程建模,根據(jù)建模目的和問題分析作出簡化假設(shè),按照內(nèi)在規(guī)律或用類比法建立微分方程,主要內(nèi)容,生物單種群增長模型3.1人口增長模型3.2傳染病模型生物多種群增長模型3.3正規(guī)戰(zhàn)與游擊戰(zhàn)3.4捕食系統(tǒng)的Volterra方程,,,為了保持自然資料的合理開發(fā)與利用，人類必須保持并控制生態(tài)平衡，甚至必須控制人類自身的增長。本節(jié)將建立幾個簡單的單種群增長模型，以簡略分析一下這方面的問題。一般生態(tài)系統(tǒng)的分析可以通過一些簡單模型的復(fù)合來研究，大家若有興趣可以根據(jù)生態(tài)系統(tǒng)的特征自行建立相應(yīng)的模型。,美麗的大自然,種群的數(shù)量本應(yīng)取離散值，但由于種群數(shù)量一般較大，為建立微分方程模型，可將種群數(shù)量看作連續(xù)變量，甚至允許它為可微變量，由此引起的誤差將是十分微小的。,離散化為連續(xù)，方便研究,3.1如何預(yù)報人口的增長--Malthus模型與Logistic模型,背景,世界人口增長概況,中國人口增長概況,研究人口變化規(guī)律,控制人口過快增長,指數(shù)增長模型——馬爾薩斯提出(1798),常用的計(jì)算公式,x(t):時刻t的人口,基本假設(shè):人口(相對)增長率r是常數(shù),不考慮移民,今年人口x0,年增長率r,k年后人口,隨著時間增加，人口按指數(shù)規(guī)律無限增長,Malthus模型實(shí)際上只有在群體總數(shù)不太大時才合理，到總數(shù)增大時，生物群體的各成員之間由于有限的生存空間，有限的自然資源及食物等原因，就可能發(fā)生生存競爭等現(xiàn)象。,所以Malthus模型假設(shè)的人口凈增長率不可能始終保持常數(shù)，它應(yīng)當(dāng)與人口數(shù)量有關(guān)。,指數(shù)增長模型的應(yīng)用及局限性,與19世紀(jì)以前歐洲一些地區(qū)人口統(tǒng)計(jì)數(shù)據(jù)吻合,適用于19世紀(jì)后遷往加拿大的歐洲移民后代,可用于短期人口增長預(yù)測,不符合19世紀(jì)后多數(shù)地區(qū)人口增長規(guī)律,不能預(yù)測較長期的人口增長過程,19世紀(jì)后人口數(shù)據(jù),阻滯增長模型(Logistic模型),人口增長到一定數(shù)量后，增長率下降的原因：,資源、環(huán)境等因素對人口增長的阻滯作用,且阻滯作用隨人口數(shù)量增加而變大,假設(shè),r~固有增長率(x很小時),xm~人口容量（資源、環(huán)境能容納的最大數(shù)量）,,,,x(t)~S形曲線,x增加先快后慢,,阻滯增長模型(Logistic模型),參數(shù)估計(jì),用指數(shù)增長模型或阻滯增長模型作人口預(yù)報，必須先估計(jì)模型參數(shù)r或r,xm,利用統(tǒng)計(jì)數(shù)據(jù)用最小二乘法作擬合,例：美國人口數(shù)據(jù)（單位~百萬）,專家估計(jì),阻滯增長模型(Logistic模型),繼續(xù),最小二乘法,設(shè)經(jīng)實(shí)際測量已得到n組數(shù)據(jù)（xi,yi），i=1,…,n。將數(shù)據(jù)畫在平面直角坐標(biāo)系中，見圖。如果建模者判斷這n個點(diǎn)很象是分布在某條直線附近，令該直線方程為y=ax+b，進(jìn)而利用數(shù)據(jù)來求參數(shù)a和b。由于該直線只是數(shù)據(jù)近似滿足的關(guān)系式，故yi-(axi+b)=0一般不成立，但我們希望,,最小,,此式對a和b的偏導(dǎo)數(shù)均為0，解相應(yīng)方程組，求得：,,,,用MATLAB作線性最小二乘擬合,1.作多項(xiàng)式f(x)=a1xm+…+amx+am+1擬合,可利用已有程序:,a=polyfit(x,y,m),1.lsqcurvefit已知數(shù)據(jù)點(diǎn)：xdata=（xdata1，xdata2，…，xdatan），ydata=（ydata1，ydata2，…，ydatan）,用MATLAB作非線性最小二乘擬合,Matlab的提供了兩個求非線性最小二乘擬合的函數(shù)：lsqcurvefit和lsqnonlin。兩個命令都要先建立M-文件fun.m，在其中定義函數(shù)f(x)，但兩者定義f(x)的方式是不同的,可參考例題.,lsqcurvefit用以求含參量x（向量）的向量值函數(shù)F(x,xdata)=（F（x，xdata1），…，F(xiàn)（x，xdatan））T中的參變量x(向量),使得,,輸入格式為:(1)x=lsqcurvefit(‘fun’,x0,xdata,ydata);(2)x=lsqcurvefit(‘fun’,x0,xdata,ydata,options);(3)x=lsqcurvefit(‘fun’,x0,xdata,ydata,options,’grad’);(4)[x,options]=lsqcurvefit(‘fun’,x0,xdata,ydata,…);(5)[x,options,funval]=lsqcurvefit(‘fun’,x0,xdata,ydata,…);(6)[x,options,funval,Jacob]=lsqcurvefit(‘fun’,x0,xdata,ydata,…);,說明：x=lsqcurvefit(‘fun’,x0,xdata,ydata,options);,lsqnonlin用以求含參量x（向量）的向量值函數(shù)f(x)=(f1(x),f2(x),…,fn(x))T中的參量x，使得最小。其中fi（x）=f（x，xdatai，ydatai）=F(x,xdatai)-ydatai,2.lsqnonlin,已知數(shù)據(jù)點(diǎn)：xdata=（xdata1，xdata2，…，xdatan）ydata=（ydata1，ydata2，…，ydatan）,輸入格式為：1）x=lsqnonlin（‘fun’，x0）；2）x=lsqnonlin（‘fun’，x0，options）；3）x=lsqnonlin（‘fun’，x0，options，‘grad’）；4）[x，options]=lsqnonlin（‘fun’，x0，…）；5）[x，options，funval]=lsqnonlin（‘fun’，x0，…）；,說明：x=lsqnonlin（‘fun’，x0，options）；,,例用下面一組數(shù)據(jù)擬合中的參數(shù)a，b，k,該問題即解最優(yōu)化問題：,1）編寫M-文件curvefun1.mfunctionf=curvefun1(x,tdata)f=x(1)+x(2)*exp(-0.02*x(3)*tdata)%其中x(1)=a;x(2)=b；x(3)=k;,2）輸入命令tdata=100:100:1000cdata=1e-03*[4.54,4.99,5.35,5.65,5.90,6.10,6.26,6.39,6.50,6.59];x0=[0.2,0.05,0.05];x=lsqcurvefit(curvefun1,x0,tdata,cdata)f=curvefun1(x,tdata),F(x，tdata)=，x=(a，b，k),解法1.用命令lsqcurvefit,3）運(yùn)算結(jié)果為：f=0.00430.00510.00560.00590.00610.00620.00620.00630.00630.0063x=0.0063-0.00340.2542,4）結(jié)論：a=0.0063,b=-0.0034,k=0.2542,返回,模型檢驗(yàn),用模型計(jì)算2000年美國人口，與實(shí)際數(shù)據(jù)比較,實(shí)際為281.4(百萬),,模型應(yīng)用——預(yù)報美國2010年的人口,加入2000年人口數(shù)據(jù)后重新估計(jì)模型參數(shù),Logistic模型在經(jīng)濟(jì)領(lǐng)域中的應(yīng)用(如耐用消費(fèi)品的售量),阻滯增長模型(Logistic模型),大量實(shí)驗(yàn)資料表明用Logistic模型來描述種群的增長，效果還是相當(dāng)不錯的。例如，高斯把5只草履蟲放進(jìn)一個盛有0.5cm3營養(yǎng)液的小試管，他發(fā)現(xiàn)，開始時草履蟲以每天230.9%的速率增長，此后增長速度不斷減慢，到第五天達(dá)到最大量375個，實(shí)驗(yàn)數(shù)據(jù)與r=2.309，a=0.006157，N(0)=5的Logistic曲線：幾乎完全吻合，見圖3.6。,圖3-6,Malthus模型和Logistic模型的總結(jié),Malthus模型和Logistic模型均為對微分方程（3.7）所作的模擬近似方程。前一模型假設(shè)了種群增長率r為一常數(shù)，（r被稱為該種群的內(nèi)稟增長率）。后一模型則假設(shè)環(huán)境只能供養(yǎng)一定數(shù)量的種群，從而引入了一個競爭項(xiàng)。,用模擬近似法建立微分方程來研究實(shí)際問題時必須對求得的解進(jìn)行檢驗(yàn)，看其是否與實(shí)際情況相符或基本相符。相符性越好則模擬得越好，否則就得找出不相符的主要原因，對模型進(jìn)行修改。,Malthus模型與Logistic模型雖然都是為了研究種群數(shù)量的增長情況而建立的，但它們也可用來研究其他實(shí)際問題，只要這些實(shí)際問題的數(shù)學(xué)模型有相同的微分方程即可，下面我們來看兩個較為有趣的實(shí)例。,年齡分布對于人口預(yù)測的重要性,只考慮自然出生與死亡，不計(jì)遷移,人口發(fā)展方程,人口發(fā)展方程,,,一階偏微分方程,,,例2:傳染病模型,問題,描述傳染病的傳播過程,分析受感染人數(shù)的變化規(guī)律,預(yù)報傳染病高潮到來的時刻,預(yù)防傳染病蔓延的手段,按照傳播過程的一般規(guī)律，用機(jī)理分析方法建立模型,已感染人數(shù)(病人)i(t),每個病人每天有效接觸(足以使人致病)人數(shù)為?,模型1,假設(shè),若有效接觸的是病人，則不能使病人數(shù)增加,建模,,,,?,模型2,區(qū)分已感染者(病人)和未感染者(健康人),假設(shè),1）總?cè)藬?shù)N不變，病人和健康人的比例分別為,2）每個病人每天有效接觸人數(shù)為?,且使接觸的健康人致病,建模,,?~日接觸率,SI模型,模型2,tm~傳染病高潮到來時刻,?(日接觸率)??tm?,病人可以治愈！,?,t=tm,di/dt最大,模型3,傳染病無免疫性——病人治愈成為健康人，健康人可再次被感染,增加假設(shè),SIS模型,3）病人每天治愈的比例為?,?~日治愈率,建模,?~日接觸率,1/?~感染期,?~一個感染期內(nèi)每個病人的有效接觸人數(shù)，稱為接觸數(shù)。,,接觸數(shù)?=1~閾值,如果感染期內(nèi)有效接觸感染的健康者人數(shù)不超過病人數(shù),患者就會全部治愈。,模型4,傳染病有免疫性——病人治愈后即移出感染系統(tǒng)，稱移出者,SIR模型,假設(shè),1）總?cè)藬?shù)N不變，病人、健康人和移出者的比例分別為,2）病人的日接觸率?,日治愈率?,接觸數(shù)?=?/?,建模,需建立的兩個方程,模型4,SIR模型,,模型4,SIR模型,相軌線的定義域,在D內(nèi)作相軌線的圖形，進(jìn)行分析,SIR模型,相軌線及其分析,s(t)單調(diào)減?相軌線的方向,P1:s0>1/??i(t)先升后降至0,P2:s0<1/??i(t)單調(diào)降至0,1/?~閾值,模型4,SIR模型,預(yù)防傳染病蔓延的手段,?(日接觸率)??衛(wèi)生水平?,?(日治愈率)??醫(yī)療水平?,傳染病不蔓延的條件——s00，積分曲線在N軸上的投影曲線（稱為軌線）將趨于K。這說明，平衡點(diǎn)N=0和N=K有著極大的區(qū)別。,圖3-17,定義1自治系統(tǒng)的相空間是指以（x1,…,xn）為坐標(biāo)的空間Rn。,特別，當(dāng)n=2時，稱相空間為相平面。,空間Rn的點(diǎn)集{(x1,…,xn)}|xi=xi(t)滿足(3.28)，i=1,…,n}稱為系統(tǒng)的軌線，所有軌線在相空間的分布圖稱為相圖。,定義2設(shè)x0是（3.28）的平衡點(diǎn)，稱：,（1）x0是穩(wěn)定的，如果對于任意的ε>0，存在一個δ>0，只要|x(0)-x0|<δ，就有|x(t)-x0|xo時，又有f(x)0時（3.30）有周期解，零點(diǎn)是穩(wěn)定的中心（非漸近穩(wěn)定）；在其他情況下，零點(diǎn)均為不穩(wěn)定的。,非線性方程組（3.29）平衡點(diǎn)穩(wěn)定性討論可以證明有下面定理成立:,定理2若（3.30）的零點(diǎn)是漸近穩(wěn)定的，則（3.29）的平衡點(diǎn)也是漸近穩(wěn)定的；若（3.30）的零點(diǎn)是不穩(wěn)定的，則（3.29）的平衡點(diǎn)也是不穩(wěn)定的。,捕食系統(tǒng)的Volterra方程,問題背景：,意大利生物學(xué)家D’Ancona曾致力于魚類種群相互制約關(guān)系的研究，在研究過程中他無意中發(fā)現(xiàn)了一些第一次世界大戰(zhàn)期間地中海沿岸港口捕獲的幾種魚類占捕獲總量百分比的資料，從這些資料中他發(fā)現(xiàn)各種軟骨掠肉魚，如鯊魚、鰩魚等我們稱之為捕食者（或食肉魚）的一些不是很理想的魚類占總漁獲量的百分比。在1914~1923年期間，意大利阜姆港收購的魚中食肉魚所占的比例有明顯的增加：,他知道，捕獲的各種魚的比例近似地反映了地中海里各種魚類的比例。戰(zhàn)爭期間捕魚量大幅下降，但捕獲量的下降為什么會導(dǎo)致鯊魚、鰩魚等食肉魚比例的上升，即對捕食者有利而不是對食餌有利呢？他百思不得其解，無法解釋這一現(xiàn)象，就去求教當(dāng)時著名的意大利數(shù)學(xué)家V.Volterra，希望他能建立一個數(shù)學(xué)模型研究這一問題。,Volterra將魚劃分為兩類。一類為食用魚（食餌），數(shù)量記為x1(t)，另一類為食肉魚（捕食者），數(shù)量記為x2(t)，并建立雙房室系統(tǒng)模型。,1、模型建立,大海中有食用魚生存的足夠資源，可假設(shè)食用魚獨(dú)立生存將按增長率為r1的指數(shù)律增長（Malthus模型），既設(shè)：,由于捕食者的存在，食用魚數(shù)量因而減少，設(shè)減少的速率與兩者數(shù)量的乘積成正比（競爭項(xiàng)的統(tǒng)計(jì)籌算律），即：,對于食餌（Prey）系統(tǒng)：,對于捕食者（Predator）系統(tǒng)：,捕食者設(shè)其離開食餌獨(dú)立存在時的死亡率為r2，即：,但食餌提供了食物，使生命得以延續(xù)。這一結(jié)果也要通過競爭來實(shí)現(xiàn)，再次利用統(tǒng)計(jì)籌算律，得到：,方程組（3.31）反映了在沒有人工捕獲的自然環(huán)境中食餌與捕食者之間的相互制約關(guān)系。下面我們來分析該方程組。,2、模型分析,方程組（3.31）是非線性的，不易直接求解。容易看出，該方程組共有兩個平衡點(diǎn)，即：,方程組還有兩組平凡解：,和,和,當(dāng)x1(0)、x2(0)均不為零時，，應(yīng)有x1(t)>0且x2(t)>0，相應(yīng)的相軌線應(yīng)保持在第一象限中。,求（3.31）的相軌線,將兩方程相除消去時間t，得：,令,用微積分知識容易證明：,有：,,與的圖形見圖3-20,易知僅當(dāng)時（3.32）才有解,當(dāng)時，軌線退化為平衡點(diǎn)。,當(dāng)時，軌線為一封閉曲線（圖3-21），即周期解。,證明具有周期解。,只需證明：存在兩點(diǎn)及，時，方程無解。,由的性質(zhì)，，而，使得：,。同樣根據(jù)的性質(zhì)知，當(dāng)下載提示(請認(rèn)真閱讀)
- 1.請仔細(xì)閱讀文檔，確保文檔完整性，對于不預(yù)覽、不比對內(nèi)容而直接下載帶來的問題本站不予受理。
- 2.下載的文檔，不會出現(xiàn)我們的網(wǎng)址水印。
- 3、該文檔所得收入（下載+內(nèi)容+預(yù)覽）歸上傳者、原創(chuàng)作者；如果您是本文檔原作者，請點(diǎn)此認(rèn)領(lǐng)！既往收益都?xì)w您。
文檔包含非法信息？點(diǎn)此舉報后獲取現(xiàn)金獎勵！
下載文檔到電腦，查找使用更方便
14.9 積分
| 下載 |
還剩頁未讀，繼續(xù)閱讀
- 配套講稿：
如PPT文件的首頁顯示word圖標(biāo)，表示該P(yáng)PT已包含配套word講稿。雙擊word圖標(biāo)可打開word文檔。
- 特殊限制：
部分文檔作品中含有的國旗、國徽等圖片，僅作為作品整體效果示例展示，禁止商用。設(shè)計(jì)者僅對作品中獨(dú)創(chuàng)性部分享有著作權(quán)。
- 關(guān) 鍵 詞：
- 西安電子科技大學(xué) 數(shù)學(xué) 建模 講義 三講
裝配圖網(wǎng)所有資源均是用戶自行上傳分享，僅供網(wǎng)友學(xué)習(xí)交流，未經(jīng)上傳用戶書面授權(quán)，請勿作他用。
相關(guān)資源
更多
正為您匹配相似的精品文檔
相關(guān)搜索
鏈接地址：http://m.appdesigncorp.com/p-3586285.html