文档详情

传染病模型.docx

发布:2017-02-05约6.52千字共9页下载文档
文本预览下载声明
传染病模型?摘要?:?当今社会,人们开始意识到通过定量地研究传染病的传播规律,建立传染病的传播模型,可以为预测和控制传染病提供可靠、足够的信息。本文利用微分方程稳定性理论对传统传染病动力学建模方式进行综述,且针对甲流,SARS等新生传染病模型进行建模和分析。??不同类型的传染病的传播过程有其各自不同的特点,我们不是从医学的角度一一分析各种传染病的传播,而是从一般的传播机理分析建立各种模型,如简单模型,SI模型,SIS模型,SIR模型等。本文中,我们应用传染病动力学模型来描述疾病发展变化的过程和传播规律,运用联立微分方程组体现疫情发展过程中各类人的内在因果联系,并在此基础上建立方程求解算法。然后,通过借助Matlab程序拟合出与实际较为符合的曲线并进行了疫情预测,评估各种控制措施的效果,从而不断完善文中的模型。??本文由简到难、全面地评价了该模型的合理性与实用性,而后对模型和数据也做了较为扼要的分析,进一步改进了模型的不妥之处。同时,在对问题进行较为全面评价的基础上又引入更为全面合理的假设,运用双线性函数模型对卫生部的措施进行了评价并给出建议,做好模型的完善与优化工作。???关键词:传染病模型,简单模型,SI,SIS,SIR,微分方程,Matlab。1问题的重述SARS(Severe Acute Respiratory Syndrome,严重急性呼吸道综合症, 俗称:非典型肺炎)是21世纪第一个在世界范围内传播的传染病。SARS的爆发和蔓延给我国的经济发展和人民生活带来了很大影响,我们从中得到了许多重要的经验和教训,认识到定量地研究传染病的传播规律、为预测和控制传染病蔓延创造条件的重要性。请你们对SARS 的传播建立数学模型,具体要求如下:1)建立传染病传播的指数模型,评价其合理性和实用性。2)建立你们自己的模型,说明为什么优于指数模型;特别要说明怎样才能建立一个真正能够预测以及能为预防和控制提供可靠、足够的信息的模型,这样做的困难在哪里?对于卫生部门所采取的措施做出评论,如:提前或延后5天采取严格的隔离措施,对疫情传播所造成的影响做出估计。附件1提供的数据供参考。3)说明建立传染病数学模型的重要性。2 定义与符号说明 N…………………………………表示为SARS病人的总数;K(感染率)……………………表示为平均每天每人的传染他人的人数;L…………………………………表示为每个病人可能传染他人的天数;N(t)………………………… 表示为每天(单位时间)发病人数;N(t)-N(t-L)………………………表示可传染他人的病人的总数减去失去传染能力的病人数;t…………………………………表示时间;R………………………………表示拟合的均方差;3 建立传染病传播的指数模型3.1模型假设1) 该疫情有很强的传播性,病人(带菌者)通过接触(空气,食物,……)将病菌传播给健康者。单位时间(一天)内一个病人能传播的人数是常数k; 2) 在 所传染的人当中不考虑已治愈的人是否被再次被传播,治愈的人数占该地区的总人数是绝对的少数,治愈者不会再被传播并不影响疫情在该时间内的感染率常数k;3) 病者在潜伏期传播可能性很小, 仍按健康人处理;4) SARS对不同的年龄组的感染率略有不同(相差不大),但我们只考虑它健康人的感染率是一样的;5) 我们所采取的隔离是非常严格的,被隔离的病人不会再感染其他人;3.2模型的分析和建立求解全国疫情从出现第一例病人起,到4月20日前后(从起点起45天左右)是疫情高峰,在此之前k值我们取k=0.16204,在此后的时间里我们取k=0.0273来计算。根据提供的数据可以建立指数模型:N(t)=n(1+K) 。在前45天我们取k=0.16204来代入,分别算出45天的病人累计数,根据45天中天病人的数量来画出图1,并与附件中所提供的数据中的日累计数来进行了比较。如图3-1所示:图3-1 根据指数模型建立的图形图3-2从两个图形中,我们可以看出,从4月20日开始计算,前45天的病人累计数和我们用k的值来代入模型画出的病人计算数基本上是吻合的。图形1中的横坐标数字表示时间的天数,如15即4月20日之后的第15天,40即4月20日之后的第40天。在45天之后的时间里,模型对k的值进行了调整,k=0.0273,我们再将k=0.0273代入模型 N(t)=n(1+K),在45天之后的时间里,我们取了30天的时间,分别算出每天的病人累计数,如图3-3所示: 图3-33.3对指数模型的验证和评价在图形3-3中的横坐标的数值表示图形1中所表示的天数之后的天数,如1即表示4月15日之后的45天之后的有第六天,也就是4月15日之后的第51天,即表示4月15日之后的第67天。首先在图形3-3结合图形3-1可以看出,图形3-1中的第45天与图形2中的第一天
显示全部
相似文档