SARS的预测控制模型_0.docx

《SARS的预测控制模型_0.docx》由会员分享，可在线阅读，更多相关《SARS的预测控制模型_0.docx（18页珍藏版）》请在优知文库上搜索。

1、SARS的料限制模型SARS的预料限制模型摘要摘要本问题是一个关于传染病限制的数学预料模型。首先，我们对附件1的模型进行了深化的分析，认为它具有肯定的合理性，但是对于预料而言，好用性却不强。为了能够达到精确预料的效果，我们建立了一个微分方程组的传染病限制模型来描述SARS传播的过程，此模型在探讨/SARS传播过程的基础上，采纳了差分计算的方法深化地分析了感染人数的改变规律，度量传染病扩散的程度并对制止其扩散的手段进行较深化的探讨。在模型中依据政府相关限制措施来确定日治愈率）（t,日接触率）（t的值，预料了传染病高潮的到来时刻。此外，针对SAKS对北京市接待海外旅游人数的影响，利用时间数列分析方

2、法建立了预料模型，并且得到9-12月北京地区海外旅游人数分别为：19.6、24.5、26.7、22.6（万人3一、问题的重述：一、问题的重述：SARS（Severe.AcuteRespiratorySyndrome,严峻急性呼吸道综合症，俗称：非典型肺炎）是21世纪第一个在世界范围内传播的传染病。SARS的爆发和扩散给我国的经济发展和人民生活带来了很大影响，这其中有很多重要的阅历和教训，特殊应当相识到定量地探讨传表示时间，以天数为单位。2)附件1模型的基本假设：(1)设病人在1.期限后失去传染作用，其缘由可能是被严格隔离、病愈不再传染或者死去等等：对于不同的疫区和疫情阶段，1.的值在1525之

3、间，为了简洁把1.固定在20天这个值。(2)不考虑疫情出现失控或反复的状态。(3)将整个SARS疫情的过程分为初期、过渡期、稳定期。初期：指从疫情起先到疫情的高峰期，此阶段，整个社会的防范程度都比较低，K值相对高：过渡期：指初期过后的10天，此时由于社会加强了宣扬力度，提高了人们的防范意识，使得K值逐步下降到很小。稳定期：指疫情得到基本限制，K降低到一个很小的稳定值，直到没有病例。3)附件1模型的建立：假定初始时刻的病例数为ON,平均每病人每天可传染K个人(K一般为小数)，平均每个病人可以干脆感染他人的时间为1.天。则在1.天之内，病例数目的增长随时间1(单位天)的关系是：tKNtN)1()(

4、0假如不考虑对传染期的限制，则病例数将依据指数规律增长。考虑传染期1.的作用后，改变将显著偏离指数率，增长速度会所以得到1(1.N的病例总数为：)0()l)(0()()1(NKN1.N1.N同理，可以这样求出：)1()1)(1()1()2(NKN1.N1.N因为在高峰期前，由假设K值是一个定值，所以，用数学归纳法可以得到：)1()1)(1()1()(iNKiNi1.Ni1.N(*)北京在3月1号发觉了第一例患者，那么可以认为，对于北京来N(0)=1,且北京疫情的起先是3月1号。将1)0(N代入(*),(*)式，用C语言对半模拟循环计算方法进行编程(源程序见附录1),可以很简洁的得到北京从疫情起

5、先到高峰期的这59天里每天的患者总数(见下表)患者数(计算值)患者数(实际值)口期患者数(计算值)患者数(实际值)日期患者数(计算值)患者数(实际值3.21-3.22154.111643.313.23174.121843.413.24194.132073.523.25224.142333.623.26274.152623.723.27284.162943.823.28314.173313.933.29354.183723.1033.30404.194183.1143.31454.204703393.1244.1504.215284823.1354.2574.225945883.1454.364

6、4.236686933.1564.4724.247517743.1674.5814.258448773.1784.6914.269499883.1894.71024.27106711143.19104.81154.28119911993.20124.91294.29134713473.2114-4.10146-从上表中可以看出，由该模型计算出的数值和实际的数值相差不大，将得到数值的误差和实际数值相比，不会超过5%可见，这个模型还是有肯定的合理性，比较正确的反映了在疫情发生的初期，患者总数和天数之间的关系。5）附件1模型的不足：1）附件1建立的早期模型事实上是个指数增长模型，虽然考虑到了传染期的

7、限制，但是这个模型的误差随着时间的推移将会增大。所以，用这个模型去预料过渡期以及稳定期的状况就会产生较大的误差。2）附件1的模型中由于将1.表示的平均每个病人在被发觉前后可以造成干脆传染的期限定在了20天这个固定的值，而没有考虑到高峰期时政府为了有效限制该病的传播而加大了宣扬的力度，可能会使得病人的有效传播期减小。3）附件1模型中对广东和香港的数据拟合得到的K值只是适用于当地，假如用于北京的状况则表现出较大的误差。这就说明白由公布的SARS病人数据得到的相对固定的K值对其它地方下次预料没有太大的实际作用。三、对SARS问题建立新的模型三、对SARS问题建立新的模型1.问题的分析与假设社会、经济

8、、文化、风俗习惯等因素都会影响SARS的传播及最终的结果，但是，最干脆的因素是：自由传染者的数量及其在健康人群中的分布，被传染者的数量,传播形式及病毒本身的传播实力等。在建立模型时不行能也没有必要考虑全部的因素，只能抓主要的因素进行合理的假设和建模。由此,我们做如下的假设：1）国家卫生部供应的全国疫情统计真实牢靠。2）将SARS全部可能的传播途径都视为与病源的干脆接触。3）在疾病传播期内所考察的地区的总人数N视为常数，不考虑人口的流淌。4）依据目前的医学调查资料，SARS康复者尚未复发。因此，我们可以假设一个SARS康复者再次感染SARS的概率为0.这些人势必会留意个人卫生远离传播源，所以他们

9、既不是易感染者，也不是已感染者，可视为他们已经退出SARS传染体系。5）相对于传统的传染病，SARS的传播时间不是很长，故假设不考虑这段时间内的人口诞生率和自然死亡率，而对于由SARS引起的死亡人数,也将其视为退出者。2.模型的参数说明:）（tI：第t天患病的人数。）（tR：第t天退出传染系统的人数。）（tS：第t天易感染人群的人数，即健康者的人数。）（ItR：第t天死亡的人数。）（2tR：第t天治愈出院的人数。）（ti：第t天病人占该地区总人口的比例。）（tr：第I天退出传染系统的人数占该地区总人口的比例。）（ts：第t天易感染人群的人数占该地区总人口的比例。日接触率，即每个病人在传染期内每

10、天有效接触的平均人数。日治愈率，即每天被治愈的病人数占病人总数的比例。传染期接触数，每个病人在传染期内有效接触人数。3.模型的建立：3.1 微分方程组模型：由假设3)明显有:NtStRtI)0()(BIJ：1 )OO(tstrti设)(II是连续、可微的函数，考察从t到tt病人人数的增加，就有：)()()(/)()(tItstIttIttI当0tBj,得到微分方程：)0()0/00OIIIIlIsdtCll为记初始时刻的患病人数除以总人数N即为：0)0(iiiisdtdi(2)对于病愈或死亡的退出者而言同理应有：)0(0)0(/统的人数为记初始时刻退出传染系RldtdR即为：0)0(rridt

11、dr(3)再记初始时刻的易感者人数是)0(00SS其比例是)0(00Ss,那么可以将易感人群的函数)(tS表示成为：0)0(SSSidtdS除以总人数N即为：0)0(sSSidtds(4) 3.2差分方程组的模型:由于微分方程组(1)(2)(3)(4)无法给出解析解,为了得到模型的数值解,特建立以下差分方程组：NtStsNtStRtItItsttStStRtIttRtIttsttItItI/)()O()()O1()1()1()1()O1()1()1()O1()1()1()1()1()1()(其中)(t,)(t是关于t的函数。四、模型的求解为了得至J(5)式中的)(t与)(t的表达式，由附件(二

12、)的数据得到)(t与日接触率)(t与时间t的关系。(5) (6)(7)(8)(9)从图中可知，在第25天前，)(t基本为一条直线，由(6)经过拟合得258119.0487.810274.310442.510422.32400.0108）（223446ttttttt在第25天后,）（t也基本为一条直线，经过拟合得25002.02404412.007335.0004559.00.0009419-）（23tttttt解得以下一组解：日期累计病例数（计弊值）累计病例数（实际值）日期累计病例数（计算值）累计病例数（实际值）4.214334825.23247025124.225515885.2424742

13、5144.236786935.25247725174.248107745.26248125204.259418775.27248425214.2610669885.28248625224.27118211145.29249025224.28128911995.30249325224.29138413475.31249525224.30147214406.1249825225.1155415536.2250025225.2163116366.3250225235.3170517416.4250425225.4178118036.5250625225.5185818976.6250825225.6

14、193819606.7250925235.7202220496.8251125235.8210821366.9251325225.9219521776.10251425225.10228022276.11251525225.11235722656.12251625215.12242023046.13251725215.13246123706.14251925215.14247123886.15251925215.15244024056.16251925215.16244324206.17252025215.17244824346.18252025215.18245224446.19252125

15、215.19245524656.20252125215.20246024906.21252125215.21246224996.22252125215.22246725046.2325212521该组解与北京实际值的比较见图（1）：五、模型的评价该模型与附件1的模型相比，具有以下优点：1 .预料性强，与可以依据当时的状况改变而改变，并且应用此模型可以有效预料以后的状况。2 .考虑了免疫人群SARS患者经治愈后，至今为止尚未发觉一例熨发的状况，因此可以将SARS患者以及由于患SARS而残废的人群视为免疫人群，从而保证了模型的精确性。模型的不足：1 .未给出）（t,）（t与政府制定的相关政策的解析关系式；2.由于影响的因素多、收集信息困难、传播途径不明确，模型得出的计算值与实际值仍有肯定偏差。对卫生部门实行的措施的评价：由图可知，在卫生部门实行严格措施（4月20日）以后，累计病例仍有一个为期25天左右的上升过程，因此