第五章_微分方程模型(清華大學(xué)數(shù)學(xué)建模教程).ppt

第五章微分方程模型,5.1傳染病模型5.2經(jīng)濟(jì)增長(zhǎng)模型5.3正規(guī)戰(zhàn)與游擊戰(zhàn)5.4藥物在體內(nèi)的分布與排除5.5香煙過(guò)濾嘴的作用5.6人口預(yù)測(cè)和控制5.7煙霧的擴(kuò)散與消失5.8萬(wàn)有引力定律的發(fā)現(xiàn),動(dòng)態(tài)模型,描述對(duì)象特征隨時(shí)間(空間)的演變過(guò)程,分析對(duì)象特征的變化規(guī)律,預(yù)報(bào)對(duì)象特征的未來(lái)性態(tài),研究控制對(duì)象特征的手段,根據(jù)函數(shù)及其變化率之間的關(guān)系確定函數(shù),微分方程建模,根據(jù)建模目的和問(wèn)題分析作出簡(jiǎn)化假設(shè),按照內(nèi)在規(guī)律或用類比法建立微分方程,5.1傳染病模型,問(wèn)題,描述傳染病的傳播過(guò)程,分析受感染人數(shù)的變化規(guī)律,預(yù)報(bào)傳染病高潮到來(lái)的時(shí)刻,預(yù)防傳染病蔓延的手段,按照傳播過(guò)程的一般規(guī)律，用機(jī)理分析方法建立模型,已感染人數(shù)(病人)i(t),每個(gè)病人每天有效接觸(足以使人致病)人數(shù)為,模型1,假設(shè),若有效接觸的是病人，則不能使病人數(shù)增加,建模,?,模型2,區(qū)分已感染者(病人)和未感染者(健康人),假設(shè),1）總?cè)藬?shù)N不變，病人和健康人的比例分別為,2）每個(gè)病人每天有效接觸人數(shù)為,且使接觸的健康人致病,建模,日接觸率,SI模型,模型2,tm傳染病高潮到來(lái)時(shí)刻,(日接觸率)tm,病人可以治愈！,?,t=tm,di/dt最大,模型3,傳染病無(wú)免疫性病人治愈成為健康人，健康人可再次被感染,增加假設(shè),SIS模型,3）病人每天治愈的比例為,日治愈率,建模,日接觸率,1/感染期,一個(gè)感染期內(nèi)每個(gè)病人的有效接觸人數(shù)，稱為接觸數(shù)。,模型3,接觸數(shù)=1閾值,感染期內(nèi)有效接觸感染的健康者人數(shù)不超過(guò)病人數(shù),模型2(SI模型)如何看作模型3(SIS模型)的特例,模型4,傳染病有免疫性病人治愈后即移出感染系統(tǒng)，稱移出者,SIR模型,假設(shè),1）總?cè)藬?shù)N不變，病人、健康人和移出者的比例分別為,2）病人的日接觸率,日治愈率,接觸數(shù)=/,建模,需建立的兩個(gè)方程,模型4,SIR模型,模型4,SIR模型,相軌線的定義域,在D內(nèi)作相軌線的圖形，進(jìn)行分析,模型4,SIR模型,相軌線及其分析,s(t)單調(diào)減相軌線的方向,P1:s01/i(t)先升后降至0,P2:s0u,Q吸一支煙毒物進(jìn)入人體總量,模型建立,t=0,x=0，點(diǎn)燃香煙,q(x,t)毒物流量,w(x,t)毒物密度,1)求q(x,0)=q(x),t時(shí)刻，香煙燃至x=ut,1)求q(x,0)=q(x),2)求q(l,t),3)求w(ut,t),4)計(jì)算Q,結(jié)果分析,煙草為什么有作用?,1）Q與a,M成正比，aM是毒物集中在x=l處的吸入量,2)過(guò)濾嘴因素，,l2負(fù)指數(shù)作用,是毒物集中在x=l1處的吸入量,3）(r)煙草的吸收作用,b,l1線性作用,帶過(guò)濾嘴,不帶過(guò)濾嘴,結(jié)果分析,4)與另一支不帶過(guò)濾嘴的香煙比較，w0,b,a,v,l均相同，吸至x=l1扔掉,提高-b與加長(zhǎng)l2，效果相同,5.6人口預(yù)測(cè)和控制,年齡分布對(duì)于人口預(yù)測(cè)的重要性,只考慮自然出生與死亡，不計(jì)遷移,人口發(fā)展方程,人口發(fā)展方程,一階偏微分方程,人口發(fā)展方程,已知函數(shù)（人口調(diào)查）,生育率（控制人口手段）,生育率的分解,總和生育率,h生育模式,人口發(fā)展方程和生育率,總和生育率控制生育的多少,生育模式控制生育的早晚和疏密,正反饋系統(tǒng),滯后作用很大,人口指數(shù),1）人口總數(shù),2）平均年齡,3）平均壽命,t時(shí)刻出生的人，死亡率按(r,t)計(jì)算的平均存活時(shí)間,4）老齡化指數(shù),控制生育率,控制N(t)不過(guò)大,控制(t)不過(guò)高,5.7煙霧的擴(kuò)散與消失,現(xiàn)象和問(wèn)題,炮彈在空中爆炸，煙霧向四周擴(kuò)散，形成圓形不透光區(qū)域。,不透光區(qū)域不斷擴(kuò)大，然后區(qū)域邊界逐漸明亮，區(qū)域縮小，最后煙霧消失。,建立模型描述煙霧擴(kuò)散和消失過(guò)程，分析消失時(shí)間與各因素的關(guān)系。,問(wèn)題分析,無(wú)窮空間由瞬時(shí)點(diǎn)源導(dǎo)致的擴(kuò)散過(guò)程，用二階偏微分方程描述煙霧濃度的變化。,觀察的煙霧消失與煙霧對(duì)光線的吸收，以及儀器對(duì)明暗的靈敏程度有關(guān)。,模型假設(shè),1）煙霧在無(wú)窮空間擴(kuò)散，不受大地和風(fēng)的影響；擴(kuò)散服從熱傳導(dǎo)定律。,2）光線穿過(guò)煙霧時(shí)光強(qiáng)的減少與煙霧濃度成正比；無(wú)煙霧的大氣不影響光強(qiáng)。,3）穿過(guò)煙霧進(jìn)入儀器的光線只有明暗之分，明暗界限由儀器靈敏度決定。,模型建立,1）煙霧濃度的變化規(guī)律,熱傳導(dǎo)定律：?jiǎn)挝粫r(shí)間通過(guò)單位法向面積的流量與濃度梯度成正比,曲面積分的奧氏公式,1）煙霧濃度的變化規(guī)律,初始條件,Q炮彈釋放的煙霧總量,單位強(qiáng)度的點(diǎn)源函數(shù),對(duì)任意t,C的等值面是球面x2+y2+z2=R2;RC,僅當(dāng)t,對(duì)任意點(diǎn)(x,y,z),C0,1）煙霧濃度的變化規(guī)律,2）穿過(guò)煙霧光強(qiáng)的變化規(guī)律,光強(qiáng)的減少與煙霧濃度成正比,3）儀器靈敏度與煙霧明暗界限,煙霧濃度連續(xù)變化,煙霧中光強(qiáng)連續(xù)變化,設(shè)光源在z=-,儀器在z=,則觀測(cè)到的明暗界限為,不透光區(qū)域邊界,4）不透光區(qū)域邊界的變化規(guī)律,對(duì)任意t,不透光區(qū)域邊界是圓周,不透光區(qū)域邊界半徑,結(jié)果分析,觀測(cè)到不透光區(qū)域邊界達(dá)到最大的時(shí)刻t1，可以預(yù)報(bào)煙霧消失的時(shí)刻t2,5.8萬(wàn)有引力定律的發(fā)現(xiàn),背景,航海業(yè)發(fā)展,天文觀測(cè)精確,“地心說(shuō)”動(dòng)搖,哥白尼：“日心說(shuō)”,伽里略：落體運(yùn)動(dòng),開普勒：行星運(yùn)動(dòng)三定律,變速運(yùn)動(dòng)的計(jì)算方法,牛頓：一切運(yùn)動(dòng)有力學(xué)原因,牛頓運(yùn)動(dòng)三定律,牛頓：研究變速運(yùn)動(dòng)，發(fā)明微積分（流數(shù)法）,開普勒三定律,牛頓運(yùn)動(dòng)第二定律,萬(wàn)有引力定律,自然科學(xué)之?dāng)?shù)學(xué)原理(1687),模型假設(shè),極坐標(biāo)系(r,),太陽(yáng)(0,0),1.行星軌道,a長(zhǎng)半軸,b短半軸,e離