传染病模型.docx
《传染病模型.docx》由会员分享,可在线阅读,更多相关《传染病模型.docx(13页珍藏版)》请在三一文库上搜索。
1、病毒扩散与传播的控制模型摘 要 随着科技的发展,病毒扩散与传播越来越受到人们的关注。 本文通过建立微分方程模型,描述了病毒扩散与传播的过程,最后通过分析,得到了控制病毒扩散与传播的方法。对问题一,我们通过分析影响变量的因素,建立微分方程模型。对问题二,我们通过把增加的影响因素加入到问题一的微分方程模型中,改善后得到了新的微分方程模型。最后把变量代入,求解微分方程模型,得到结果。当t=13天时,确诊患者人数达到峰值6793000人;t=150天时,确诊人数减少到116800人。对问题三、问题四、问题五,通过把改变后的条件代入到问题二中的微分方程模型中,可以得到其对应结果。问题三的结果,当t=13
2、天时,确诊人数达到峰值6769000人,t=150天时,确诊人数减少到108400人。问题四的结果,当t=13天时,确诊人数达到峰值6795000人,t=150天时,确诊人数减少到116200人。问题五的结果,当t=12天时,确诊人数达到峰值6793000人,当t=150天时,确诊人数减少到113500人对问题六,结合前面所得到的结果,我们分析在其它因素都不变的情况下只改变一种因素,分析得到该种因素的灵敏度,最后得出各个因素的灵敏度。可以得到,尽快开始隔离、治愈时间 、 ,加强隔离强度p、减少人均日接触率r都可以改善病情。对问题七,在问题六的基础上,可以得出相应的减轻病情的方法和建议。关键词:
3、 微分方程模型 微分方程组求解(MATLAB)13一、 问题的重述已知某种不完全确知的具有传染性病毒的潜伏期为d1d2天,病患者的治愈时间为d3天。该病毒可通过直接接触、口腔飞沫进行传播、扩散,该人群的人均每天接触人数为r。为了控制病毒的扩散与传播将该人群分为五类:确诊患者、疑似患者、治愈者、死亡和正常人,可控制参数是隔离措施强度p(潜伏期内的患者被隔离的百分数)。 要求:1. 在合理的假设下试建立该病毒扩散与传播的控制模型;2. 利用你所建立的模型针对如下数据进行模拟条件1:d1=1, d2=11, d3=30, r=10,条件2:已经知道的初始发病人数为890、疑似患者为2000条件3:隔
4、离措施强度p=60%条件4:患者2天后入院治疗,疑似患者2天后被隔离,试给出患者人数随时间变化的曲线图,并明确标识图中的一些特殊点的具体数据,分析结果的合理性。3. 若将2中的条件4改为条件:患者1.5天后入院治疗,疑似患者1.5天后被隔离,模拟结果有何变化?4. 若仅将2中的条件3改为条件:隔离措施强度p=40%,模拟结果有何变化?5. 若仅将2中的条件1改为条件:d1=1, d2=11, d3=30, r=250,模拟结果有何变化?6. 分析问题中的参数对计算结果的敏感性。7. 针对如上数据给政府部门写一个不超过400字的建议报告。二、问题的分析2.1 问题一的分析问题一的解决,在于对疑似
5、患者、确诊患者、治愈者、正常人、死亡者的理解,在理解的基础上,我们分析影响它们的因素有哪些,最后通过建立微分方程模型来解释这些影响关系。2.2 问题二的分析问题二的解决基于问题一,在问题已的基础上,我们对于增加的影响因素进项分析,并且改进问题一中的微分方程模型。我们会得到考虑更加全面的微分方程模型,最后把已知的数据带入方程组中,最后用MATLAB进行求解,可以得到结果。2.3 问题三的分析问题三的解决,即把已知的条件代入问题二的微分方程组中,最后用MATLAB求解,可以得到结果。2.4 问题四的分析问题四的解决,即把已知的条件代入问题二的微分方程组中,最后用MATLAB求解,可以得到结果。2.
6、5 问题五的分析问题五的解决,即把已知的条件代入问题二的微分方程组中,最后用MATLAB求解,可以得到结果。2.6 问题六的分析通过对问题二、问题三、问题四、问题五结果的分析,通过相互之间的对比,我们会发现不同影响因素对病情的影响,最终得出结论。2.7 问题七的分析通过问题六的结果,我们会得到影响病情的因素以及影响关系,通过这些关系,提出减轻病情的建议和方法。三、模型的假设1总人数N不变,人群分为确诊患者、疑似患者、治愈者、死亡和正常人;2. 疑似患者是被病毒感染,但是没有发病的人;3确诊患者自动被隔离,不具备传染病毒的能力,所以人群中只有未被隔离的疑似患者能够传染疾病;4治愈者具有了免疫能力
7、,不会再被该病毒感染;5平均潜伏期为 5. 在问题二中,假设人群总数为1000 0000四、符号说明S疑似患者在人群中的比例I确诊患者在人群中的比例R治愈者在人群中的比例H正常人在人群中的比例D死亡人数在人群中的比例N人群的总人数(假设为1000 0000)d1d2病毒的潜伏期d3治愈时间r人群的人均日接触人数p隔离强度 平均潜伏期五、模型的建立与求解5.1 问题一的模型建立与求解易知显然有 S+I+R+H+D=1根据已知条件以及假设,分析可得:疑似患者在人群中的比例决定于未被隔离的疑似患者和疑似患者发病的比例有关。确诊患者在人群中的比例决定于疑似患者发病和患病者得到治愈的比例有关。治愈者在人
8、群中的比例决定于确诊患者被治愈的比例有关。正常人的比例决定于未被隔离的疑似患者的比例有关。死亡者的比例则是人群中每一时刻人群原来整数减去疑似患者、确诊患者、治愈者、正常人人数的比例。有以上分析以及题目中所数据和假设中的数据可将此问题满足的微分方程模型建立,如下: 到此,我们建立该病毒扩散与传播的控制模型;5.2 问题二的模型建立与求解问题二模型的建立是建立在问题一的基础上,由题中给出的额外条件患者天后得到治疗和疑似患者天后被隔离,所以我们基于问题一的模型。考虑到此事疑似患者感染人群的天数增加 天,所以此时人群中每天疑似患者的增加人数比问题一多 ;相反人群中正常人每天比问题一减少。考虑到此时患病
9、者天后开始治疗,所以这时治愈时间变为+ 天,此时人群中每天治愈者的增加人数变为;相反人群中每天确诊患者的变为 。所以得到问题二的微分方程模型,如下: 把问题二中已知的条件1、条件2、条件3、条件4可得:, , =0.000089, =0.00002 , =0由假设可得 =0.9997将这些数据代入上述方程,可得微分方程组,如下: =0.000089,=0.00002 , =0, =0.9997用MATLAB求解此微分方程模型,可以得到患者人数随时间变化的曲线图如下图所示: 问题二的结果(由于死亡者的比例几乎解决于0,甚至是一个很小的负数,所以与横坐标轴重合)在图中标出三个特殊点,(0,0.00
10、0089)(13,0.6793)(150,116800)。分析可得,起初即t=0时,人群中确诊患者的比例为0.000089,即890人;当t=13时,确诊患者人数比例达到最多0.6793,即6793000人。并且013天,确诊患者人数递增;13150天,确诊患者人数递减,最后在150天时,确诊患者比例为0.01168,即116800人。5.3 问题三的模型建立与求解由已知可得,改变问题二中的条件四, 可得微分方程组: =0.000089,=0.00002 , =0, =0.9997用MATLAB求解此微分方程模型,可以得到患者人数随时间变化的曲线图如下图所示: 问题三的结果在图中标出俩个特殊点
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 传染病 模型
链接地址:https://www.31doc.com/p-11684248.html