傳染病傳播的數(shù)學(xué)模型_上課

1、微分方程模型學(xué)習(xí)目的1,加深對微分方程概念的理解，掌握針對一些問題通過建立微分方程的方法及微分方程的求解過程；2,了解微分方程模型解決問題思維方法及技巧；3,領(lǐng)會建立微分方程模型的逐步改進法的核心及優(yōu)點，并掌握該方法;4,理解微分方程的解的穩(wěn)定性的意義，會用穩(wěn)定性判定模型的解是否有效；5.體會微分方程建摸的藝術(shù)性。在自然學(xué)科（如物理、化學(xué)、生物、天文）以及在工程、經(jīng)濟、軍事、社會等學(xué)科量的問題可以用微分方程來描述。正如列寧所說：“自然界的統(tǒng)一性顯示在關(guān)于各種現(xiàn)象領(lǐng)域的微分方程式的驚人的類似中."（列寧選集第二卷，人民1972年版第295頁）。要建立微分方程模型，讀者必須掌握元素法（有

2、關(guān)元素法，在高等數(shù)學(xué)中已有介紹）。所謂元素法，從某種角度上講，就是分析的方法，它是以自然規(guī)律的普遍性為根據(jù)并且以局部規(guī)律的獨立的假定為基礎(chǔ)。在解決各種實際問題時，微分方程用得極其廣泛。讀者通過下面的幾個不同領(lǐng)域中的模型介紹便有所體會，要想掌握好它，在這方面應(yīng)作大量的練習(xí)。§17.1、傳染病傳播的數(shù)學(xué)模型學(xué)習(xí)目標(biāo)1,通過學(xué)習(xí)建立傳染病傳播的數(shù)學(xué)模型的思維方法，能歸納出該類建模的關(guān)鍵性步驟及思維方法；并能指出求解傳染病傳播的數(shù)學(xué)模型的方法技巧；2,能用已知的傳染病傳播的數(shù)學(xué)模型，預(yù)報某種傳染病的傳播；3.學(xué)會從簡單到復(fù)雜的處理問題的方法。由于人體的疾病難以控制和變化莫測，因此醫(yī)學(xué)中的數(shù)學(xué)

3、模型較為復(fù)雜。生物醫(yī)學(xué)中的數(shù)學(xué)模型分為兩大類：傳染病傳播的數(shù)學(xué)模型和疾病數(shù)學(xué)模型。以下僅討論傳染病的傳播問題。人們將傳染病的統(tǒng)計數(shù)據(jù)進行處理和分析，發(fā)現(xiàn)在某一民族或地區(qū)，某種傳染病傳播時，每次所涉及的人數(shù)大體上是一常數(shù)。這一現(xiàn)象如何解釋呢？關(guān)于這個問題，醫(yī)學(xué)工作者試圖從醫(yī)學(xué)的不同角度進行解釋都得不到令人滿意的解釋。最后由于數(shù)學(xué)工作者的參與，在理論上對上述結(jié)論進行了嚴(yán)格的證明。同時又由于傳染病數(shù)學(xué)模型的建立，分析所得結(jié)果與實際過程比較吻合，這個現(xiàn)象才得到了比較滿意的解釋。傳染病傳播所涉及的因素很多，如傳染病人的多少，易受傳染者的多少，傳染率的大小，排除率的大小，人口的出生和死亡等。如果還要考慮人

4、員的遷入與遷出，潛伏期的長短以及預(yù)防疾病的傳播等因素的影響，那么傳染病的傳播就變得非常復(fù)雜。如果一開始就把所有的因素考慮在，那么將陷入多如亂麻的頭緒中不能自拔，倒不如舍去眾多的次要因素，抓住主要因素，把問題簡化，建立相應(yīng)的數(shù)學(xué)模型。將所得結(jié)果與實際比較，找出問題，修改原有假設(shè)，再建立一個與實際比較吻合的模型。下面由簡單到復(fù)雜將建模的思考過程作一示，讀者可以從中得到很好的啟發(fā)。模型一、考慮最簡單的情形：假設(shè)(1),每個病人在單位時間傳染的人數(shù)是常數(shù)K。；假設(shè)(2),一人得病后，經(jīng)久不愈，并在傳染期不會死亡。記i(t)表示t時刻病人數(shù)，K。表示每個病人單位時間傳染的人數(shù)，i(0)i。，即最初有卜個

5、傳染病人。則在t時間增加的病人數(shù)為t)i(t)K0i(t)t，其解為i(t)i°ek0ti(t于是得微分方程駟K"t),Koi(t)dti(0)i。結(jié)果表明：傳染病的傳播是按指數(shù)函數(shù)增加的。這個結(jié)果與傳染病傳播初期比較吻合，傳染病傳播初期，傳播快，被傳染人數(shù)按指數(shù)函數(shù)增長。但由方程(1)的解可以推出，當(dāng)t時，i(t),這顯然是不符合實際情況的。問題在于兩條假設(shè)均不合理。特別是假設(shè)(1),每個病人在單位時間傳染的人數(shù)是常數(shù)與實際不符。因為在傳播初期，傳染病人少，未被傳染者多；而在傳染病傳播中期和后期，傳染病人逐漸增多，未被傳染者逐漸減少，因而在不同時期的傳染情況是不同的。為了

6、與實際情況吻合，我們在原有基礎(chǔ)上修改假設(shè)建立新的模型。模型二、用i(t),s(t)表示t時刻傳染病人數(shù)和未被傳染人數(shù)，i(0)i0假設(shè)(1),每個病人單位時間傳染的人數(shù)與這時未被傳染的人數(shù)成正比，即K0Ks(t)；假設(shè)(2),一人得病后，經(jīng)久不愈，并在傳染期不會死亡;假設(shè)(3),總?cè)藬?shù)為n,即s(t)i(t)n.由以上假設(shè)得微分方程包9Ks(t)i(t)(2)dts(t)i(t)ni(0)io用分離變量法求得其解為i(t)1-1eKntio其圖形如圖17.1所示.模型二可以用來預(yù)報傳染較快的疾病前期傳染病高峰到來的時間。i(t)本iot圖17.1醫(yī)學(xué)上稱di/(dtt)為傳染病曲線,17.2所

7、示。由(3)式可得圖17.2它表示傳染病人增加率與時間的關(guān)系，如圖Kn2didt-1eKnti0(4)Kntd2i(t)令一2出20,得極大點為ln(n1)10Kn由此可見，當(dāng)傳染病強度K或總?cè)藬?shù)n增加時，匕都將變小，即傳染病高峰來得快，這與實際情況吻合。同時，如果知道了傳染強度K(K由統(tǒng)計數(shù)據(jù)得出)，即可預(yù)報傳染病高峰匕到來的時間，這對于防治傳染病是有益處的。模型二的缺點是：當(dāng)t時，由(3)式可知，i(t)n,即最后人人都要生病,這顯然是不符合實際情況的。造成該問題的原因是假設(shè)(2)中假設(shè)了人得病后經(jīng)久不愈。為了與實際問題更加吻合，對上面的數(shù)學(xué)模型再進一步修改，這就要考慮人得了病后有的會死亡

8、；另外不是每個人被傳染后都會傳染別人，因為其中一部分會被隔離。還要考慮人得了傳染病由于醫(yī)治和人的自身抵抗力會痊愈，并非象前面假設(shè)的那樣，人得病后經(jīng)久不愈。為此作出新的假設(shè)，建立新的模型。模型三在此模型中，雖然要考慮比前面兩個模型復(fù)雜得多的因素，但仍要把問題簡單化。設(shè)患過傳染病而完全痊愈的任何人具有長期免疫力，不考慮反復(fù)受傳染的情形，并設(shè)傳染病的潛伏期很短，可以忽略不計，即一個人患了病之后立即成為傳染者。在這種情況下，把居民分成三類：第一類是由能夠把疾病傳染給別人的那些傳染者組成的，用I表示t時刻第一類人數(shù)；第二類是由并非傳染者但能夠得病而成為傳染者的那些人組成的，用s表示t時刻第二類人數(shù)；第三

9、類包括患病死去的人，病愈后具有長期免疫力的人，以及在病愈并出現(xiàn)長期免疫力以前被隔離起來的人，用R表示t時刻第三類人數(shù)。假設(shè)疾病傳染服從下列法則：(1)在所考慮的時期人口總數(shù)保持在固定水平N,即不考慮出生及其它原因引起的死亡，以及遷入遷出等情況；(2)易受傳染者人數(shù)s(t)的變化率正比于第一類的人數(shù)I(t)與第二類人數(shù)s(t)的乘積；(3)由第一類向第三類轉(zhuǎn)變的速率與第一類的人數(shù)成正比由(1)、(2)、(3)條得微分方程組dsdtdldtdRdtrslrslI(6)其中r、為兩個比例常數(shù)，r為傳染率，為排除率由(6)式的三個方程相加得二s(t)I(t)R(t)0dt則s(t)I(t)R(t)N(

10、人口總數(shù))故R(t)Ns(t)I(t)(7)由此可知，只要知道了s(t)和I(t),即可求出R(t)o而(6)式的第一和第二個方程與R(t)無關(guān)。因此，由dsdtdldtrslrsI(8)7日dlrslI行dsrsl1/l(s)s-Insc。r當(dāng)tt°時，l(t°)l°,s(to)s。，記sl(s)lososlnso下面討論積分曲線(9)的性質(zhì)。由(8)式知0,s(l(s)10,ss0,s所以當(dāng)s時，l(s)是s的增函數(shù)，s(9)時，l(s)是s的減函數(shù)。lso100由連續(xù)函數(shù)中間值定理及單調(diào)性知，存在唯一點s,0sso使得l(s)0。而當(dāng)sss0時，l(s)&g

11、t;0o由(7)知l=0時，ds/dt0,dl/dt0。所以(s,0)為方程組(7)的平衡點.當(dāng)tt0時，方程(9)的圖形如圖17.3當(dāng)t由t0變化到時，點(s(t),l(t)沿曲線(9)移動,并沿s減少方向移動，因為s(t)隨時間的增加而單調(diào)減少。因此，如果s0小于，則l(t)單調(diào)減小到零，s(t)單調(diào)減小到s。所以，如果為數(shù)不多的一群傳染者l0分散在居民so中，且so,則這種疾病會很快被消滅。由上分析可以得出如下結(jié)論：圖17.3如果比,則隨著s(t)減小到時，l(t)增力口，且當(dāng)s時，l(t)達到最大值。當(dāng)s(t)時，l(t)才開始減小。由以上分析可以得出如下結(jié)論：只是當(dāng)居民中的易受傳染者

12、的人數(shù)超過閾值一時傳染病才會蔓延。用一般的常識來檢r驗上面的結(jié)論也是符合的。當(dāng)人口擁擠、密度高，缺乏應(yīng)有的科學(xué)文化知識，缺乏必要的醫(yī)療條件，隔離不良而排除率低時，傳染病會很快蔓延；反之，人口密度低，社會條件好，有良好的公共衛(wèi)生設(shè)施和較好的管理而排除率高時，則疾病在有限圍出現(xiàn)卻很快被消滅。如果起初易受傳染者的人數(shù)S0大于但接近于閾值,即如果（S0）與相比是小量，則最終患病的人數(shù)近似于2（s0）,這就是著名的傳染病學(xué)中的閾值定理。生物數(shù)學(xué)家Kermack和Mekendrick在1927年首先證明了這個定理。定理（傳染病學(xué)中的閾值定理）：設(shè)Sor,且假設(shè)r/同1相比是小量,并設(shè)最初傳染者人數(shù)Io很小

13、，則最終患病的人數(shù)為2r。即易受傳染者的人數(shù)最初比閾值高多少，那最終就會比閾值低多少。證明略。根據(jù)閾值定理就可以由起初易受傳染者的人數(shù)來估計最終患病的人數(shù)。這個定理解釋了研究人員長期以來難以解釋的為什么對于某一民族或地區(qū)，某種傳染病傳播時，每次所波及的人數(shù)大體上是一常數(shù)的現(xiàn)象。在傳染病發(fā)生過程中，不可能準(zhǔn)確的調(diào)查每一天或每一星期得病的人數(shù)。因為只有那些來醫(yī)院就醫(yī)者才能被人知道他們得了病，并把他們隔離起來防止傳染。因此，統(tǒng)計的記錄是每一天或每一星期新排除者的人數(shù)，而不是新得病的人數(shù)。所以，為了把數(shù)學(xué)模型所預(yù)示的結(jié)果同疾病的實際情況進行比較，必須解出（6）式中的第三個方程：因為所以dRdtdsdR

14、s)dsdRdt出rsIIsdsdR一,，sR/s(R)s°e%eR/(10)方程（10）雖是可分離變量的，但是不能用顯式求解如果傳染病不嚴(yán)重，則R/是小量，取泰勒級數(shù)2R1Rew1R1-的前三項，取近似值得2dRdtNso包1R包R2其解為2R(t)一So其中So因此dRdt2SoSo2sNSo221一sec2tg+1So.arctg一一1,(11)方程(11)在tdR/dt平面上定義了一條對稱鐘形曲線，稱為疾病傳染曲線.疾病傳染曲線很好的說明了實際發(fā)生的傳染病。少下來。每天報告的新病案的數(shù)目逐漸上升到峰值，然后又減死亡人數(shù)300dt8oo6oo4oo2oo510152030星期數(shù)

15、r圖17.4圖17.5Kermack和Mekendrick把(11)得到的dR/dt的值,年上半年在孟買發(fā)生的瘟疫資料進行比較，他們設(shè)同取自19。5年下半年至19。6dR2890sech2(0.2t34)dt其中t按星期計，在圖17.5中，dR/dt的實際數(shù)字(圖上用”表示)同理論曲線非常一致。這就表明了模型三是在固定的居民中傳染病傳播的準(zhǔn)確而可靠的數(shù)學(xué)模型。對于同一事物，可用不同的數(shù)學(xué)工具來描述它。下面介紹一般隨機傳染病模型。模型四、般隨機傳染病的數(shù)學(xué)模型:以上建立的常微分方程描述的傳染病的傳播是確定性的模型。但人生病是隨機的，因而建立隨機的傳染病的數(shù)學(xué)模型才能更實際的反映傳染病的傳播。設(shè)X

16、(t)表示t時刻易受傳染者人數(shù)，Y(t)表示t時刻已受傳染者人數(shù)，n表示易受傳染者總數(shù)，又設(shè)t時刻有i(i>0)個易受傳染者移入已受傳染者中來。這種傳染病傳播的機制如下：(1)在群體中個體均勻的混和；(2)在區(qū)間(t,tt),一個新傳染病例出現(xiàn)的概率為xyto(t),其中0是傳染率；(3)在區(qū)間(t,tt),排除一個個體的概率為yto(t),其中(0)是排除率；(4)在區(qū)間(t,tt),有多次轉(zhuǎn)移(即多個傳染或排除)發(fā)生的概率為o(t);(5)在區(qū)間(t,tt),無變化的概率為1xytot,令Pxy(t)X(t)x,Y(t)y,x0,y0,t0。從(2)和(3)知有兩種可能的轉(zhuǎn)移(xx1

17、,yy1)及(yy1)。因此，表征隨機傳染病流行過程的差分方程為：dPxy(t)dt(x1)(y1)Px1,y1(t)y(x)Py(t)(y1)Px,y1(t)dPni(t)dt(12)i(n)Pni(t)其中0xyni,0xn,0yni。在(12)式中/為相對排除率，并已對時間變量作了變換，使方程對傳染率是無量綱的。方程(12)的初始條件為Pxy(0)1,xn,yi0,其他(6)。方程(12)所描述的隨機傳染病流行數(shù)學(xué)模型可以推得確定性模型1942年Wilson和Burker討論了潛伏期的重要性，他們用微分差分方程來描述傳染病傳播的數(shù)學(xué)模型.x(t)Arx(t)Ax(t)x(t)其中A是易受感染者乘以恢復(fù)健康的比率，是易受感染者變成傳染者的潛伏期，r是易受感染者與病人的接觸率。Cooke除了考慮潛伏期，還引進了阻尼閾值的概念，這個概念說明個別成員從易受感染者成為傳染者之前可能反復(fù)發(fā)病，Cooke模