1、實(shí)驗(yàn)報(bào)告1 實(shí)驗(yàn)題目?jī)和檎盍餍新拥臄?shù)學(xué)模型2 實(shí)驗(yàn)問題陳述試組建一個(gè)能描述兒童麻疹流行蔓延的數(shù)學(xué)模型，我們將考慮在接種疫苗成為有效的防疫手段之前的麻疹的流行。下表一給出了英國(guó)倫敦在1647年-1660年間每年麻疹病的死亡人數(shù)。表一：倫敦每年麻疹死亡人數(shù)（1647-1660）年代4748495051525354555657585960人數(shù)5923333362852111531580674可以看出，它是以2年為周期的周期性流行。已知麻疹的潛伏期是0.5周，在這段時(shí)期內(nèi)一個(gè)被感染的孩子表面看來是正常的，但卻會(huì)傳染給別人。過了這段時(shí)間后，患病的孩子一直被隔離到痊愈為止。痊愈后孩子是免疫的。假設(shè)每個(gè)

2、感染者隨機(jī)地與他人接觸。證明你的模型有某種周期性質(zhì)。如果不然，就修改你的模型。因?yàn)槁檎畹牧餍锌隙ㄊ勤呌谥芷谑降爻霈F(xiàn)的。估計(jì)你組建的模型中的參數(shù)，以擬合0.5周的潛伏期及2年周期性流行的觀測(cè)結(jié)果。判斷估計(jì)出的參數(shù)是否實(shí)際。3 實(shí)驗(yàn)?zāi)康耐ㄟ^表中數(shù)據(jù)，建立麻疹流行蔓延模型，以擬合0.5周的潛伏期及2年周期性流行的觀測(cè)結(jié)果，判斷估計(jì)出的參數(shù)是否實(shí)際。4 實(shí)驗(yàn)內(nèi)容模型假設(shè)：（1）除感病特征外，人群中的個(gè)體間沒有差異，感病者與易感者的個(gè)體在人群中混合是均勻的。（2）人群的數(shù)量足夠大，只考慮傳染過程的平均效應(yīng)。（3）易感者感病的機(jī)會(huì)與他接觸感病者的機(jī)會(huì)成正比。（4）疾病的傳染率為常數(shù)。（5）一般的麻疹爆發(fā)在

3、幾十天，我們不考慮在一次麻疹爆發(fā)時(shí)間內(nèi)某地區(qū)的出生人口和死亡人口，以及人口的遷入和遷出。（6）感病痊愈者（即移出者）移出模型，而不再成為易感者人群中的成員。變量說明：S(t)：易感者在人群中所占的比例I(t)：感病者在人群中所占的比例R(t)：移出者在人群中所占的比例 K：疾病的傳染率 h：?jiǎn)挝粫r(shí)間內(nèi)痊愈的百分?jǐn)?shù)一個(gè)傳染期內(nèi)每個(gè)病人有效接觸易感者的平均人數(shù)，成為接觸數(shù)-初始時(shí)刻問題分析:對(duì)麻疹流行蔓延的周期性質(zhì)進(jìn)行說明。通過對(duì)SIR模型及麻疹流行的機(jī)理分析，在一次麻疹爆發(fā)以后絕大多數(shù)人體內(nèi)具有了麻疹免疫抗體，因此絕大多數(shù)新生嬰兒體內(nèi)具有抗體，考慮到引起流行周期的原因是易感人群的積累，易感人群來

4、源于新生兒因母體抗體逐漸消失而易感、既往沒有患過麻疹的兒童和成人。通常認(rèn)為在自然感染狀態(tài)下，這些易感者積累到一個(gè)以上出生隊(duì)列時(shí)，就達(dá)到爆發(fā)的“臨界”。若將同一年出生的人群組定義為一個(gè)出生隊(duì)列，出生隊(duì)列出現(xiàn)的周期性在一定程度上可以說明麻疹流行的周期性。所以我們用積累一個(gè)出生隊(duì)列的時(shí)間來表示一次爆發(fā)的臨界。模型建立通過對(duì)問題的分析，模型可以表示為：，其中考慮到初始條件，可知上述三個(gè)方程是相容的，因此可以化簡(jiǎn)為：由于方程組無法求出解析解，故可以在S-I的像平面上討論解的性質(zhì)，相軌線的定義域?yàn)椋河梢陨戏匠炭芍壘€的方程為：，其解為：。5 實(shí)驗(yàn)結(jié)果分析與討論由題目中表格給出的麻疹死亡人數(shù)與年份的對(duì)應(yīng)關(guān)系

5、，用MATLAB編程畫出的曲線圖如圖一所示：【圖一】 由于對(duì)S(t)和I(t)的求1解非常困難，所以先用數(shù)值計(jì)算的方法來預(yù)估計(jì)S(t)和I(t)的一般變化規(guī)律。在方程（1）中設(shè)k=1，h=0.3，I(0)=0.02，S(0)=0.98。編寫MATLAB程序并運(yùn)行得到如圖二， 【圖二】【注】圖中藍(lán)色曲線為I(t)，即病人比例；綠色曲線為S(t),即健康人比例。從時(shí)間流程圖中可以看出，隨著時(shí)間的增加，S(t)單調(diào)遞減，I(t)在時(shí)達(dá)到峰值以后會(huì)隨時(shí)間減小，當(dāng)時(shí)，S(t)值很小，而I(t)=0。說明在一次麻疹疫情爆發(fā)以后絕大多數(shù)的人體內(nèi)已經(jīng)具有麻疹免疫抗體，被移出除傳染系統(tǒng)。c. 對(duì)結(jié)果的分析參數(shù)中

6、取h為0.3，則潛伏期為1/h=3.33天，約等于0.5周。同時(shí)我們可以看到把h取為0.3得到的曲線符合實(shí)際情況，說明潛伏期為麻疹病毒的潛伏期是0.5周是正確的觀點(diǎn)。另外，當(dāng)時(shí)，說明在一次麻疹疫情爆發(fā)以后絕大多數(shù)的人體內(nèi)已經(jīng)具有麻疹免疫抗體，被移出除傳染系統(tǒng)。查資料知，絕大多數(shù)的嬰兒在9個(gè)月時(shí)血內(nèi)的母親抗體已測(cè)不出，有些嬰兒體內(nèi)的抗體存在時(shí)間可以長(zhǎng)達(dá)15個(gè)月，所以可以取1年為一個(gè)出生隊(duì)列產(chǎn)生的時(shí)間，用時(shí)間坐標(biāo)來表示出生隊(duì)列與麻疹流行周期的關(guān)系如圖三： 【圖三】【注】：表示第一年的年初；表示第一年年末；表示第二年年初；表示第二年年末；在時(shí)刻出生的嬰兒到時(shí)刻抗體消失，時(shí)刻出生的嬰兒到時(shí)刻抗體消失。

7、易感人群從時(shí)刻開始積累，在時(shí)刻易感人群剛好積累一個(gè)出生隊(duì)列，因此易感人群積累一個(gè)出生隊(duì)列的時(shí)間為2年。當(dāng)易感人群積累到一個(gè)出生隊(duì)列時(shí)，就是第二次麻疹爆發(fā)的“臨界”，因此可以說麻疹流行的周期為2年。綜上所述，兒童麻疹流行蔓延的模型具有周期性，且以0.5周的潛伏期和2年的周期性流行。故模型所得結(jié)果與題目要求是一致的。6 實(shí)驗(yàn)程序（Matlab或者其它軟件語言陳述）用MATLAB編寫程序如下：（1） 畫圖一：x=47:60;y=5 92 3 33 33 62 8 52 11 153 15 80 6 74;plot(x,y,rp-.)xlabel(年份); ylabel(倫敦每年麻疹病死亡人數(shù));（2） S(t)和I(t)的變化規(guī)律及畫圖三M文件為：(chuanran.m)function y=chuanran(t,x)a=1;b=0.3;y=a*x(1)*x(2)-b*x(1);-a*x(1)*x(2);命令框中輸入： ts=0