传染病模型
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
病毒扩散与传播的控制模型
摘 要
随着科技的发展,病毒扩散与传播越来越受到人们的关注。
本文通过建立微分方程模型,描述了病毒扩散与传播的过程,最后通过分析,得到了控制病毒扩散与传播的方法。
对问题一,我们通过分析影响变量的因素,建立微分方程模型。
对问题二,我们通过把增加的影响因素加入到问题一的微分方程模型中,改善后得到了新的微分方程模型。最后把变量代入,求解微分方程模型,得到结果。当t=13天时,确诊患者人数达到峰值6793000人;t=150天时,确诊人数减少到116800人。
对问题三、问题四、问题五,通过把改变后的条件代入到问题二中的微分方程模型中,可以得到其对应结果。问题三的结果,当t=13天时,确诊人数达到峰值6769000人,t=150天时,确诊人数减少到108400人。问题四的结果,当t=13天时,确诊人数达到峰值6795000人,t=150天时,确诊人数减少到116200人。问题五的结果,当t=12天时,确诊人数达到峰值6793000人,当t=150天时,确诊人数减少到113500人
对问题六,结合前面所得到的结果,我们分析在其它因素都不变的情况下只改变一种因素,分析得到该种因素的灵敏度,最后得出各个因素的灵敏度。可以得到,尽快开始隔离、治愈时间1t 、2t ,加强隔离强度p 、减少人均日接触率r 都可以改善病情。
对问题七,在问题六的基础上,可以得出相应的减轻病情的方法和建议。
关键词: 微分方程模型 微分方程组求解(MATLAB )
一、问题的重述
已知某种不完全确知的具有传染性病毒的潜伏期为d1~d2天,病患者的治愈时间为d3天。该病毒可通过直接接触、口腔飞沫进行传播、扩散,该人群的人均每天接触人数为r。为了控制病毒的扩散与传播将该人群分为五类:确诊患者、疑似患者、治愈者、死亡和正常人,可控制参数是隔离措施强度p(潜伏期内的患者被隔离的百分数)。
要求:
1.在合理的假设下试建立该病毒扩散与传播的控制模型;
2.利用你所建立的模型针对如下数据进行模拟
条件1:d1=1, d2=11, d3=30, r=10,
条件2:已经知道的初始发病人数为890、疑似患者为2000
条件3:隔离措施强度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 问题一的分析
问题一的解决,在于对疑似患者、确诊患者、治愈者、正常人、死亡者的理解,在理解的基础上,我们分析影响它们的因素有哪些,最后通过建立微分方程模型来解释这些影响关系。
2.2 问题二的分析
问题二的解决基于问题一,在问题已的基础上,我们对于增加的影响因素进
项分析,并且改进问题一中的微分方程模型。我们会得到考虑更加全面的微分方程模型,最后把已知的数据带入方程组中,最后用MATLAB 进行求解,可以得到结果。 2.3 问题三的分析
问题三的解决,即把已知的条件代入问题二的微分方程组中,最后用MATLAB 求解,可以得到结果。 2.4 问题四的分析
问题四的解决,即把已知的条件代入问题二的微分方程组中,最后用MATLAB 求解,可以得到结果。 2.5 问题五的分析
问题五的解决,即把已知的条件代入问题二的微分方程组中,最后用MATLAB 求解,可以得到结果。 2.6 问题六的分析
通过对问题二、问题三、问题四、问题五结果的分析,通过相互之间的对比,我们会发现不同影响因素对病情的影响,最终得出结论。 2.7 问题七的分析
通过问题六的结果,我们会得到影响病情的因素以及影响关系,通过这些关系,提出减轻病情的建议和方法。
三、模型的假设
1.总人数N 不变,人群分为确诊患者、疑似患者、治愈者、死亡和正常人; 2. 疑似患者是被病毒感染,但是没有发病的人;
3.确诊患者自动被隔离,不具备传染病毒的能力,所以人群中只有未被隔离的疑似患者能够传染疾病;
4.治愈者具有了免疫能力,不会再被该病毒感染; 5.平均潜伏期为
12
1
d d 5. 在问题二中,假设人群总数为1000 0000
四、符号说明
五、模型的建立与求解
5.1 问题一的模型建立与求解
易知显然有
S+I+R+H+D=1
根据已知条件以及假设,分析可得:
疑似患者在人群中的比例决定于未被隔离的疑似患者和疑似患者发病的比例有关。
确诊患者在人群中的比例决定于疑似患者发病和患病者得到治愈的比例有关。
治愈者在人群中的比例决定于确诊患者被治愈的比例有关。
正常人的比例决定于未被隔离的疑似患者的比例有关。
死亡者的比例则是人群中每一时刻人群原来整数减去疑似患者、确诊患者、治愈者、正常人人数的比例。
有以上分析以及题目中所数据和假设中的数据可将此问题满足的微分方程模型建立,如下:
4(1)dS S N
Nr p SH N dt d =-- 43
dI S I
N
N N dt d d =- 3
dR I
N N dt d = (1)dH
N
Nr p SH dt
=--
1D S I R H =----
到此,我们建立该病毒扩散与传播的控制模型;
5.2 问题二的模型建立与求解
问题二模型的建立是建立在问题一的基础上,由题中给出的额外条件患者1
t 天后得到治疗和疑似患者2t 天后被隔离,所以我们基于问题一的模型。
考虑到此事疑似患者感染人群的天数增加 2t 天,所以此时人群中每天疑似患者的增加人数比问题一多2
rSH
N
t ;相反人群中正常人每天比问题一减少2
rSH
N
t 。 考虑到此时患病者1t 天后开始治疗,所以这时治愈时间变为1t +3d 天,此时
人群中每天治愈者的增加人数变为31
I
N
d t +;相反人群中每天确诊患者的变为431
S I N
N d d t -+ 。 所以得到问题二的微分方程模型,如下:
24(1)dS rSH S
N
N Nr p SH N dt t d =+-- 431
dI S I
N
N N dt d d t =-+ 31
dR I
N N dt d t =+
2
(1)dH NrSH N
Nr p SH dt t -=-- 1D S I R H =----
把问题二中已知的条件1、条件2、条件3、条件4可得: