MERS传播的数学模型的建立及分析
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
2015高教社杯全国大学生数学建模竞赛
承诺书
我们仔细阅读了《全国大学生数学建模竞赛章程》和《全国大学生数学建模竞赛参赛规则》(以下简称为“竞赛章程和参赛规则”,可从全国大学生数学建模竞赛下载)。
我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括、电子、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问题。
我们知道,抄袭别人的成果是违反竞赛章程和参赛规则的,如果引用别人的成果或其他公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们重承诺,严格遵守竞赛章程和参赛规则,以保证竞赛的公正、公平性。如有违反竞赛章程和参赛规则的行为,我们将受到严肃处理。
我们授权全国大学生数学建模竞赛组委会,可将我们的论文以任何形式进行公开展示(包括进行网上公示,在书籍、期刊和其他媒体进行正式或非正式发表等)。
我们参赛选择的题号是(从A/B/C/D中选择一项填写): A
我们的参赛报名号为(如果赛区设置报名号的话):
所属学校(请填写完整的全名):工业学院
参赛队员(打印并签名) :1.王××
2.卢××
3.唐××
指导教师或指导教师组负责人(打印并签名):
(论文纸质版与电子版中的以上信息必须一致,只是电子版中无需签名。以上容请仔细核对,提交后将不再允许做任何修改。如填写错误,论文可能被取消评奖资格。)
日期:2015年7 月27 日赛区评阅编号(由赛区组委会评阅前进行编号):
2015高教社杯全国大学生数学建模竞赛
编号专用页
赛区评阅编号(由赛区组委会评阅前进行编号):
全国统一编号(由赛区组委会送交全国前编号):全国评阅编号(由全国组委会评阅前进行编号):
MERS传播的数学模型的建立与分析
摘要
本文针对MERS的传播建立了传统的SIR仓室数学模型。
针对问题一,对附件一提供的早期模型,认为“传染概率”的说法欠妥,传染期限L的确定缺乏医学上的支持,使模型的说服力降低。模型中借鉴的参数来预测的疫情走势,不失为一种方法。但在不同国家因政策,地域的不同,病毒的传播和控制呈现不同的特点,使不同国家不同城市之间的可比性降低。而且由于MERS和SARS的治病机理不同,传染率和死亡率不同,患病人数规模不同,所以对MERS的传播分析,不能简单的套用附件一所用的模型。
针对问题二,我们在WHO(World Health Orgnazation)的官方上查找到了国MERS疫情从2015年5月20日至2015年7月5日的详细数据[1]。对MERS的传播建立传统的SIR仓室模型,采用最小二乘法拟合参数,利用MATLAB编程求解,画出参数感染率λ的参数散点图和参数移出率α的散点图
对第三个问题,本文研究对MERS疫情对国入境旅游收入的影响,建立了灰色预测GM(1,1)模型。
关键字:SIR仓室模型常微分方程参数拟合灰色预测
1.问题重述
MERS(Middle East Respiratory Syndrome)病毒是一种新型的冠状病毒,这种病毒已经被命名为中东呼吸综合征冠状病毒,大多数MERS病毒感染病例发生在沙特。
2015年,MERS在国又有新一轮的爆发和蔓延,给国的经济发展和人民生活带来了较大影响。对MERS的传播建立数学模型,具体要求如下:
(1)对附件所提供的一个SARS传播的早期模型,请对其评价是否适用MERS。(2)收集MERS的国疫情数据,建立MERS传播的数学模型,说明优于附件1中的模型的原因;特别要说明怎样建立一个真正能够预测以及能为预防和控制提供可靠、足够的信息的模型。对于卫生部门所采取的措施做出评论,如:提前或延后n天采取严格的隔离措施,对疫情传播所造成的影响做出估计。
(3)收集MERS对国旅游方面影响的数据,建立相应的数学模型并进行预测。
2.模型假设
1.假设一个MERS康复者不会二度感染,他们已退出传染体系,因此将其归为“退出者”;
2.模型不考虑所研究这段时间的自然出生率和死亡率,MERS引起的死亡人数归为“退出者”;
3.假设在疾病传播期所考察地区总人数视为常数;
4.假设每个病人单位时间有效接触的人数为常数;
5. 假设国在MERS疫情流行期间和结束之后,旅游业数据的变化只与MERS疫情的影响有关,不考虑其它随机因素的影响。
3.变量说明
S:表示易感染人群(susceptible)占总人数的比例;
I;表示感染人群(infected)占总人数的比例;
R:表示移出人群占总人数的比例;
λ:表示感染者对易感染者有效感染的感染率;
α:表示移出率,即移出者的增加率
()()
x n:表示第n个单位时间的旅游业的收益;
1
()()
1
X n表示在n个单位时间时的累积旅游业的收益;
1
ϕ:称为系统发展灰数;
ψ:称为生控制变量
4.对早期模型的评价
附件1的模型主要采用“数据拟合”和“借鉴参数”的方法对疫情走势进行预测。
在数据拟合方面,该模型中有两个疑点:
1、感染期限L的确定。由于被严格隔离、治愈、死亡等原因,感染者在某
一时段后不再具有对易感人群的传染力,故对病毒的传染加上感染期限是合理的。但在对该参数的确定上,作者为了较好地拟合各阶段的数据,通过人为调试来确定L的取值,缺乏医学上的支持,使模型的说服力减弱,合理性和可靠性大大降低。
2、文中认为“K代表某种环境下一个人传染他人的平均概率”。但从模型的公式中可以看出,参数K的实际意义是一个病人平均每天传染其他人的个数。两者之间有实质的区别,文中的说法显然不妥。
从预测思想来看,该模型是借鉴先发地区——、的有关参数对的疫情进行预测的。由于、的疫情和控制都在之前,已经过了高峰期,到5月8日为止每日新增病例已降至10来例,基本处于后期控制阶段。而当时的疫情刚过了高峰期,正处于社会剧烈调整时期,数据较为凌乱,略有下降趋势,但不明显。可见在当时,采取这种借鉴是无奈之举。
但是由于城市之间的政策,风俗习惯等不同,城市之间的可比性不强,借鉴存在很大的局限性。如在,由于对传播机制认识不足,中途又出现高度感染的特殊情况。另外使用借鉴法无法对首发城市进行预测。
MERSS和SARS又有许多不同, 例如,MERS的病死率约为40.7%,传染性没有SARS强,但SARS的病死率为14%-15%,低于MERS,传染性则强于MERS。MERS和SARS治病机理不同,传染率死亡率不同,病毒潜伏时间也不同,所以不能简单的套用附件所给模型来分析和预测MERS的传播。
5.模型的建立与求解
5.1问题2的模型建立与求解
5.1.1 问题2的模型思路