MERS传播的数学模型的建立及分析

1、承诺书我们仔细阅读了全国大学生数学建模竞赛章程和全国大学生数学建模竞赛参赛规则（以下简称为“竞赛章程和参赛规则”，可从全国大学生数学建模竞赛网站下载）。我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛章程和参赛规则的，如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛章程和参赛规则，以保证竞赛的公正、公平性。如有违反竞赛章程和参赛规则的行为，我们将受到严肃处理。我们授权全

2、国大学生数学建模竞赛组委会，可将我们的论文以任何形式进行公开展示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。我们参赛选择的题号是（从A/B/C/D中选择一项填写）：A我们的参赛报名号为（如果赛区设置报名号的话）：所属学校（请填写完整的全名）：成都工业学院参赛队员（打印并签名）：1.王XX2 点XX3 .唐xx指导教师或指导教师组负责人（打印并签名）：（论文纸质版与电子版中的以上信息必须一致，只是电子版中无需签名。以上内容请仔细核对，提交后将不再允许做任何修改。如填写错误，论文可能被取消评奖资格。）日期：2015年7月27日赛区评阅编号（由赛区组委会评阅前进行编号）：编号

3、专用页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：merS专播的数学模型的建立与分析摘要本文针对MERS勺传播建立了传统的SIR仓室数学模型。针对问题一，对附件一提供的早期模型，认为“传染概率”的说法欠妥，传染期限L的确定缺乏医学上的支持，使模型的说服力降低。模型中借鉴广东香港的参数来预测北京的疫情走势，不失为一种方法。但在不同国家因政策，地域的不同，病毒的传播和控制呈现不同的特点，使不同国家不同城市之间的可比性降低。而且由于MER制SARS勺治病机理不同，

4、传染率和死亡率不同，患病人数规模不同，所以对MERS勺传播分析，不能简单的套用附件一所用的模型。针对问题二，我们在WHO(WorldHealthOrgnazation)的官方网站上查找到了韩国MERSS情从2015年5月20日至2015年7月5日的详细数据1。对MERS的传播建立传统的SIR仓室模型，采用最小二乘法拟合参数，利用MATLAB程求解，画出参数感染率九的参数散点图和参数移出率豆的散点图对第三个问题，本文研究对MERSS情对韩国入境旅游收入的影响，建立了灰色预测GM(1,1)模型。关键字：SIR仓室模型常微分方程参数拟合灰色预测1 .问题重述MERSMiddleEastRespira

5、torySyndrome)病毒是一种新型的冠状病毒，这种病毒已经被命名为中东呼吸综合征冠状病毒，大多数MERSKI毒感染病例发生在沙特。2015年，MERSE韩国又有新一轮的爆发和蔓延，给韩国的经济发展和人民生活带来了较大影响。对MERS勺传播建立数学模型，具体要求如下：(1)对附件所提供的一个SARS专播的早期模型，请又t其评价是否适用MERS(2)收集MERS勺韩国疫情数据，建立MERSt播的数学模型，说明优于附件1中的模型的原因；特别要说明怎样建立一个真正能够预测以及能为预防和控制提供可靠、足够的信息的模型。对于卫生部门所采取的措施做出评论，如：提前或延后n天采取严格的隔离措施，对疫情传

6、播所造成的影响做出估计。(3)收集MERSt韩国旅游方面影响的数据，建立相应的数学模型并进行预测。2 .模型假设1 .假设一个MERSt复者不会二度感染，他们已退出传染体系，因此将其归为“退出者”；2 .模型不考虑所研究这段时间内的自然出生率和死亡率，MERSH起的死亡人数归为“退出者”；3 .假设在疾病传播期内所考察地区总人数视为常数；4 .假设每个病人单位时间有效接触的人数为常数；5 .假设韩国在MERSS情流行期间和结束之后，旅游业数据的变化只与MERSS情的影响有关，不考虑其它随机因素的影响。3 .变量说明S：表示易感染人群(susceptible)占总人数的比例；I;表示感染人群(i

7、nfected)占总人数的比例；R:表示移出人群占总人数的比例；九：表示感染者对易感染者有效感染的感染率；口：表示移出率，即移出者的增加率x1')(n)：表示第n个单位时间的旅游业的收益；X1°(n)表示在n个单位时间时的累积旅游业的收益；邛：称为系统发展灰数；中,称为内生控制变量4 .对早期模型的评价附件1的模型主要采用“数据拟合”和“借鉴参数”的方法对北京疫情走势进行预测。在数据拟合方面，该模型中有两个疑点：1、感染期限L的确定。由于被严格隔离、治愈、死亡等原因，感染者在某一时段后不再具有对易感人群的传染力，故对病毒的传染加上感染期限是合理的。但在对该参数的确定上，作者为

8、了较好地拟合各阶段的数据，通过人为调试来确定L的取值，缺乏医学上的支持，使模型的说服力减弱，合理性和可靠性大大降低。2、文中认为“K代表某种环境下一个人传染他人的平均概率”。但从模型的公式中可以看出，参数K的实际意义是一个病人平均每天传染其他人的个数。两者之间有实质的区别，文中的说法显然不妥。从预测思想来看，该模型是借鉴先发地区一一广东、香港的有关参数对北京的疫情进行预测的。由于广东、香港的疫情和控制都在北京之前，已经过了高峰期，到5月8日为止每日新增病例已降至10来例，基本处于后期控制阶段。而当时北京的疫情刚过了高峰期，正处于社会剧烈调整时期，数据较为凌乱，略有下降趋势，但不明显。可见在当时

9、，采取这种借鉴是无奈之举。但是由于城市之间的政策，风俗习惯等不同，城市之间的可比性不强，借鉴存在很大的局限性。如在香港，由于对传播机制认识不足，中途又出现高度感染的特殊情况。另外使用借鉴法无法对首发城市进行预测。MERSS口SAR双有许多不同,例如，MERS勺病死率约为40.7%,传染性没有SARS虽，但SARS勺病死率为14%-15%低于MERS传染性则强于MERSMERS和SARS台病机理不同，传染率死亡率不同，病毒潜伏时间也不同，所以不能简单的套用附件所给模型来分析和预测MERS勺传播。5 .模型的建立与求解5.1 问题2的模型建立与求解5.1.1 问题2的模型思路对于传染病感染区，由于

10、为了避免传染病的更大的扩散，一方面政府会对人口的流动做出限制，另一方面个人由于对传染病的警惕也不会进入感染区，所以，感染区人口流动很小。即可把它看作一个封闭区，则可以建立SIR仓室模型2,里面的总人数不变，里面的人分类为：易感染者：即正常人，但可能会被感染。感染者：已经感染这种病的人，可以传染给周围的人。移出者：包括感染者中死亡的人，感染者中自愈的人和先天对这种病毒有免疫的人，他们将不在受这种传染病的影响。5.1.2 问题2的模型建立根据传染病的感染而致的各类人口比例的变化可以建立微分方程。这里用S(t)S示易感染人群(susceptible)的比例关于时间t的函数，同样用I(t联示感染人群(

11、infected)的比例关于时间t的函数，用R(t)表示移出人群的比例关于时间t的函数。方程一：根据感染人群比例增长相等，可建立如下方程：dI(t)dt- S(t )I (t) - ： I (t)这里九表示感染者对易感染者有效感染的感染率，即单位时间内单位病人传染的人数与易感者之比值；口表示移出率，即单位时间内移出者占染病者的比率；I(t)表示感染者比例关于时间t的函数；S(t)表示易感染人者比例关于时间t的函数。方程二：根据易感染者比例减少相等，可建立如下方程：dSd?-S(t)I(t)方程三：根据移出者的增长相等，建立如下方程:dRt)dt：I(t)这里Rt)表示移出者比例关于时间t的函数

12、。方程四：易感染者比例，感染者比例和移出者比例之和恒为1:S(t)I(t)Rt)=1dI (t)dt S(t )I (t ) - ： I(t )(5)j孚=-九S(t)I(t)4dtdRt)aI(t)dtSt)I(t)Rt)=15.1.3 问题2的模型求解5.2.1问题2的模型微分方程初始值确定我们统计的数据口是从2015年5月20日到2015年7月7日的数据(见表一)，共计49天，对于移除者比例Rt),在病情开始没有死亡人数，也没有治愈人数，所以移除者比例Rt)的初始条件R0)=0,对于感染者比例I(t),根据数据看出开始只有1人，所以感染者比例I(t)的初始条件I(0)=1/N,对于易感染

13、者S(t)的初始条件可有S(t)+I(t)+R(t)=1算出，即S(0)=11/N。对MERSS情每天数据的统计日期当天确诊病例人数累计确诊病例人数当大夕匕亡人数累计死亡人数20/05/2015330021/05/2015030022/05/2015030023/05/2015030024/05/2015030025/05/2015140026/05/2015150027/05/2015050028/05/2015270029/05/20156130030/05/20153160031/05/20158240001/06/20154281102/06/20157350103/06/201554

14、01204/06/201513531305/06/20157601406/06/201516760407/06/201516920408/06/20156983709/06/2015161141810/06/20151212631111/06/2015613221312/06/20151314511413/06/2015314831714/06/2015515342115/06/2015415712216/06/2015516222417/06/2015216452918/06/2015316702919/06/2015016713020/06/2015216903021/06/2015317

15、203022/06/2015317503023/06/2015417903024/06/2015017913125/06/2015118013226/06/2015118123427/06/2015018103428/06/2015018103429/06/2015018103430/06/2015018103401/07/2015118203402/07/2015218413503/07/2015018403504/07/2015118503505/07/2015018503506/07/2015018503507/07/201501851365.2.2问题2的模型参数拟合对于模型上的参数感

16、染率九和移出率口受很多因素影响，很显然不是一个常数。但是对于不是常数的参数很难拟合，所以，我们把所采集的数据（见表一）分为适当个组，把每个组的参数看作常数，这样就很方便拟合出每个组的参数，最后再综合各个组拟合的参数，把参数拟合成关于自变量时间t的函数（当然这里的组数分得越多则参数拟合得更加精确）。这样就完整拟合出了每个时刻的参数值。对于每个组的参数拟合，根据参数感染率人的定义（单位时间内单位病人传染的人数与易感者之比值）而易感染人数远大于单位病人传染人数，所以九<1。移出率口（单位时间内移出者占染病者的比率）同样可以得知a<1,所以可以定两个参数的范围都在0至M之间，这里可以用尝试

17、参数值去逼近最佳参数，具体步骤如下：Stepl:因为模型的参数感染率儿和移出率a范围都在0到1之间，这里我们先用0.1的梯度把参数分成（）九个值，两个参数共组成81组值。Step2：把分成的每组参数值代入模型的微分方程组，用matlab中的函数可算出一个数值解。Step3:用上面求出来的函数计算出的值和实际值作差计算平方和（即最小二乘法），找出所有组参数值中最小的二乘数，即那组的感染率九和移出率a就是在梯度0.1下的最佳参数值。Step4：画出已经计算出来的函数图，与实际数据的函数图对比，如果相差太远这进行下一步，反之转到Step7。Ste

18、p5:为了更加靠近最终的参数最佳值，这里再在上次拟合出来的参数值左右一个梯度的范围再以上次梯度的0.1倍进行参数再分组。例如：假设上次算出来的最佳感染率=0.5移出率a=0.2,这次的参数范围就是0.4九E0.60.1<a<0.3,这次就用0.01的梯度把参数分为20组参数。Step6:回至UStep2。Step7:到了这里我们已经算出了一组数据的最佳参数值。然后重复以上步骤，算出每组数据的最佳参数值。Step8:根据各组数据的最佳参数值画出参数值得散点图，拟合成适应函数。到这里，全部参数拟合结束。算法流程图如图一所示：开始相差不大图1算法流程图根据以上算法用matlab编程（程序

19、见附件一），可以画出参数感染率九的参数散点图（见图2）和参数移出率a的散点图（见图3）。程序没有写出来，没有结果5.3问题3的模型建立与求解5.3.1 问题3的模型建立对于问题3,MER盟情对韩国入境旅游游客人数的影响，这里建立了一阶一元灰色预测模型LJGM（1,1）,对旅游入境旅游游客人数进行预测。以单位是t表示的时间序列就为：Xi0（1）,Xi0（2）,|hXifnj（6）这里x1（）（n）表示第n个单位时间的韩国入境旅游游客人数令：Xi（t）=zxQ（i）1 1即通过上面公式（7）累积形成新的列：汉。,X”2）,川,X“n>（8）这里X11）（n）表示表示在n个单位时间时的韩国累积

20、入境旅游的游客人数则可以建立GM（1,1）模型：反/）。）十邛丫（）（t）=中-Xi二,di1这里丫丁）口）=2甘）口）+X，（t-1q（10）邛称为系统发展灰数，中称为内生控制变量。5.3.2问题3的模型求解确定参数2:可利用一阶微分方程的解法求得其解，然后将（7）式代入，使用最小二乘法，便可以得出：=ATA'ATM（11）其中:12冲+"）（2）1a-2W）（2）+»）（3）1*+-2/i'）（n-n+xi0（n）VLu.X，（2）XQMn=：*X10（n）LJ生成列xj）（ty可以按以下公式计算出：（12）（13）X10（t尸X，-。尸口十（14）其中

21、，t=2,3,III,n。计算预测：根据采集的数据（见表2）,可以通过以上公式（11）可以用matlab算出参数，然后用matlab解出微分方程。月份游客（万人）月份游客（万人）1月;524月672月1565月303月616月20表22015年韩国旅游各月份的游客数通过matlab程序(见附件2),可以画出出X1()(t)(在第t个月的韩国入境旅游的游客人数)的函数图(表3)。表3预计游客人数结果月份12345678游客人数(万人)79.6960.1754.9845.6937.9631.5426.2121.786、模型的分析、推广与改进要真正建立能够预测以及能为预防和控制提供可靠，足够的信息的

22、模型。要怎样做，以及困难。事实上，真正情况要比我们的模型复杂的多，比如，由于初期传播机理不明，医护人员为收到好的保护，很多受到传染。而传播率事实上也是按照某种趋势变化的，他与人们的警惕心理等因素有关。要建立这样的模型，我们认为：1、要不断的研究传染病理，完善模型。2、建立实时反馈机制，不断修改模型参数，以求更准确的预测以后的数据。3、尽量以实际的研究数据建模，力求准确反映事件的实质。这样，当人为修改参数进行模拟预测时，可以为采取什么措施效果最好做出较准确的建议。4、研究各种措施实际对疾病传播的影响。要做到这些，也恰恰是建模的难点，因为当涉及到新的传染病的时候，以上数据或参数都必须通过长期的研究

23、或拟合才能得到的。比如措施对疾病传播的影响，是要统计大量实际数据才能获得的。以上也可以看作我们模型的改进方向。7、参考文献1 世界卫生组织./csr/disease/coronavirus_infections/MERS-CoV-cases-rok-21Jul15.xlsx?ua=12 周宁李林，基于灰色微分方程的传染病预测，8、附录附件1:附件2:y=dsolve('Dy=81.9185-0.1852*y','y(1)=52','t')y=819185/1852-(722881*exp(-(463*t)/250

24、0)*exp(463/2500)/1852>>x=81.9185-0.5*0.1852*(819185/1852-(722881*exp(-(463*t)/2500)*exp(463/2500)/1852)+(819185/1852-(722881*exp(-(463*(t-1)/2500)*exp(463/2500)/1852)>>x>>t=123t=4567812378456>>x=81.9185-0.5*0.1852*(819185/1852*exp(463/2500)/1852)+(819185/1852*exp(463/2500)/1852)-(722881*exp