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

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

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

3、 日期： 2021 年 7 月 27 日赛区评阅编号由赛区组委会评阅前进展编号：2021 高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号由赛区组委会评阅前进展编号：赛区评阅记录可供赛区评阅时使用：评阅人评分备注全国统一编号由赛区组委会送交全国前编号：全国评阅编号由全国组委会评阅前进展编号：MERS传播的数学模型的建立与分析摘要本文针对MERS的传播建立了传统的SIR仓室数学模型。针对问题一，对一提供的早期模型，认为“传染概率的说法欠妥，传染期限L确实定缺乏医学上的支持，使模型的说服力降低。模型中借鉴的参数来预测的疫情走势，不失为一种方法。但在不同国家因政策，地域的不同，病毒的传播

4、和控制呈现不同的特点，使不同国家不同城市之间的可比性降低。而且由于MERS和SARS的治病机理不同，传染率和死亡率不同，患病人数规模不同，所以对MERS的传播分析，不能简单的套用一所用的模型。针对问题二，我们在WHO(World Health Orgnazation)的官方上查找到了国MERS疫情从2021 年5月20日至2021 年7月5日的详细数据1。对MERS的传播建立传统的SIR仓室模型，采用最小二乘法拟合参数，利用MATLAB编程求解，画出参数感染率的参数散点图和参数移出率的散点图对第三个问题，本文研究对MERS疫情对国入境旅游收入的影响，建立了灰色预测GM1，1模型。关键字：SIR

5、仓室模型 常微分方程参数拟合 灰色预测1.问题重述MERSMiddle East Respiratory Syndrome病毒是一种新型的冠状病毒，这种病毒已经被命名为中东呼吸综合征冠状病毒，大多数MERS病毒感染病例发生在沙特。2021 年，MERS在国又有新一轮的爆发和蔓延，给国的经济开展和人民生活带来了较大影响。对MERS的传播建立数学模型，具体要求如下：1对所提供的一个SARS传播的早期模型，请对其评价是否适用MERS。2收集MERS的国疫情数据，建立MERS传播的数学模型，说明优于1中的模型的原因；特别要说明怎样建立一个真正能够预测以及能为预防和控制提供可靠、足够的信息的模型。对于卫

6、生部门所采取的措施做出评论，如：提前或延后n天采取严格的隔离措施，对疫情传播所造成的影响做出估计。3收集MERS对国旅游方面影响的数据，建立相应的数学模型并进展预测。2.模型假设1.假设一个MERS康复者不会二度感染，他们已退出传染体系，因此将其归为“退出者；2.模型不考虑所研究这段时间的自然出生率和死亡率，MERS引起的死亡人数归为“退出者；3.假设在疾病传播期所考察地区总人数视为常数；4.假设每个病人单位时间有效接触的人数为常数；5. 假设国在MERS疫情流行期间和完毕之后，旅游业数据的变化只与MERS疫情的影响有关，不考虑其它随机因素的影响。3.变量说明：表示易感染人群 (suscept

7、ible)占总人数的比例；表示感染人群(infected)占总人数的比例；：表示移出人群占总人数的比例；：表示感染者对易感染者有效感染的感染率；：表示移出率，即移出者的增加率：表示第个单位时间的旅游业的收益；表示在个单位时间时的累积旅游业的收益；：称为系统开展灰数；：称为生控制变量4.对早期模型的评价1的模型主要采用“数据拟合和“借鉴参数的方法对疫情走势进展预测。在数据拟合方面，该模型中有两个疑点：1、感染期限L确实定。由于被严格隔离、治愈、死亡等原因，感染者在*一时段后不再具有对易感人群的传染力，故对病毒的传染加上感染期限是合理的。但在对该参数确实定上，作者为了较好地拟合各阶段的数据 ,通过

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

9、，城市之间的可比性不强，借鉴存在很大的局限性。如在，由于对传播机制认识缺乏，中途又出现高度感染的特殊情况。另外使用借鉴法无法对首发城市进展预测。MERSS和SARS又有许多不同, 例如，MERS的病死率约为40.7%，传染性没有SARS强，但SARS的病死率为14%-15%，低于MERS，传染性则强于MERS。MERS和SARS治病机理不同，传染率死亡率不同，病毒潜伏时间也不同，所以不能简单的套用所给模型来分析和预测MERS的传播。5.模型的建立与求解51 问题2的模型建立与求解5.1.1 问题2的模型思路对于传染病感染区，由于为了防止传染病的更大的扩散，一方面政府会对人口的流动做出限制，另一

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

11、一：根据感染人群比例增长相等，可建立如下方程：这里表示感染者对易感染者有效感染的感染率，即单位时间单位病人传染的人数与易感者之比值；表示移出率，即单位时间移出者占染病者的比率；表示感染者比例关于时间的函数；表示易感染人者比例关于时间的函数。方程二：根据易感染者比例减少相等，可建立如下方程：方程三：根据移出者的增长相等，建立如下方程：这里表示移出者比例关于时间的函数。方程四：易感染者比例，感染者比例和移出者比例之和恒为1：综上：5.1.3问题2的模型求解5.2.1问题2的模型微分方程初始值确定 我们统计的数据是从2021 年5月20日到2021 年7月7日的数据见表一，共计49天，对于移除者比例

12、，在病情开场没有死亡人数，也没有治愈人数，所以移除者比例的初始条件，对于感染者比例，根据数据看出开场只有1人，所以感染者比例的初始条件，对于易感染者的初始条件可有算出，即。表一 对MERS疫情每天数据的统计日期当天确诊病例人数累计确诊病例人数当天死亡人数累计死亡人数20/05/2021 330021/05/2021 030022/05/2021 030023/05/2021 030024/05/2021 030025/05/2021 140026/05/2021 150027/05/2021 050028/05/2021 270029/05/2021 6130030/05/2021 31600

13、31/05/2021 8240001/06/2021 4281102/06/2021 7350103/06/2021 5401204/06/2021 13531305/06/2021 7601406/06/2021 16760407/06/2021 16920408/06/2021 6983709/06/2021 161141810/06/2021 1212631111/06/2021 613221312/06/2021 1314511413/06/2021 314831714/06/2021 515342115/06/2021 415712216/06/2021 516222417/06/2

14、021 216452918/06/2021 316702919/06/2021 016713020/06/2021 216903021/06/2021 317203022/06/2021 317503023/06/2021 417903024/06/2021 017913125/06/2021 118013226/06/2021 118123427/06/2021 018103428/06/2021 018103429/06/2021 018103430/06/2021 018103401/07/2021 118203402/07/2021 218413503/07/2021 01840350

15、4/07/2021 118503505/07/2021 018503506/07/2021 018503507/07/2021 01851365.2.2问题2的模型参数拟合 对于模型上的参数感染率和移出率受很多因素影响，很显然不是一个常数。但是对于不是常数的参数很难拟合，所以，我们把所采集的数据见表一分为适当个组，把每个组的参数看作常数，这样就很方便拟合出每个组的参数，最后再综合各个组拟合的参数，把参数拟合成关于自变量时间的函数当然这里的组数分得越多则参数拟合得更加准确。这样就完整拟合出了每个时刻的参数值。对于每个组的参数拟合，根据参数感染率的定义单位时间单位病人传染的人数与易感者之比值而易感

16、染人数远大于单位病人传染人数，所以。移出率单位时间移出者占染病者的比率同样可以得知，所以可以定两个参数的围都在0到1之间，这里可以用尝试参数值去逼近最正确参数，具体步骤如下：Step1：因为模型的参数感染率和移出率围都在0到1之间，这里我们先用0.1的梯度把参数分成0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9九个值，两个参数共组成81组值。Step2：把分成的每组参数值代入模型的微分方程组，用matlab中的函数可算出一个数值解。Step3：用上面求出来的函数计算出的值和实际值作差计算平方和即最小二乘法，找出所有组参数值中最小的二乘数，即那组的感染率和移出率就是在梯度

17、0.1下的最正确参数值。Step4：画出已经计算出来的函数图，与实际数据的函数图比照，如果相差太远这进展下一步，反之转到Step7。Step5：为了更加靠近最终的参数最正确值，这里再在上次拟合出来的参数值左右一个梯度的围再以上次梯度的0.1倍进展参数再分组。例如：假设上次算出来的最正确感染率=0.5移出率=0.2，这次的参数围就是，这次就用0.01的梯度把参数分为20组参数。Step6：回到Step2。Step7：到了这里我们已经算出了一组数据的最正确参数值。然后重复以上步骤，算出每组数据的最正确参数值。Step8：根据各组数据的最正确参数值画出参数值得散点图，拟合成适应函数。到这里，全部参数

18、拟合完毕。算法流程图如图一所示：开场把参数分组将每一组参数值代入微分方程组，求方程组的数值解由每组参数算出函数的值与实际值比拟，记录最小二乘数反响的两个参数画出函数图，与实际函数图比照再次将参数分组，代入微分方程组，求数值解记录这组数据的额参数值每组数是否找到了参数相差不大相差很大否完毕图1 算法流程图根据以上算法用matlab编程程序见一，可以画出参数感染率的参数散点图见图2和参数移出率的散点图见图3。程序没有写出来，没有结果5.3问题3的模型建立与求解5.3.1问题3的模型建立对于问题3，MERS疫情对国入境旅游游客人数的影响，这里建立了一阶一元灰色预测模型GM(1，1)，对旅游入境旅游游

19、客人数进展预测。以单位是表示的时间序列就为：这里表示第个单位时间的国入境旅游游客人数令：即通过上面公式7累积形成新的列：这里表示表示在个单位时间时的国累积入境旅游的游客人数则可以建立GM(1,1)模型：这里称为系统开展灰数，称为生控制变量。5.3.2问题3的模型求解 确定参数： 可利用一阶微分方程的解法求得其解，然后将7式代入，使用最小二乘法，便可以得出：其中： 生成列可以按以下公式计算出：其中，。计算预测：根据采集的数据见表2，可以通过以上公式11可以用matlab算出参数，然后用matlab解出微分方程。表 2 2021 年国旅游各月份的游客数月份游客万人月份游客万人1月524月672月5

20、65月303月616月20通过matlab程序见2，可以画出出在第个月的国入境旅游的游客人数的函数图表3。表3 预计游客人数结果月份12345678游客人数万人79.6960.1754.9845.6937.9631.5426.2121.786、模型的分析、推广与改良要真正建立能够预测以及能为预防和控制提供可靠，足够的信息的模型。要怎样做，以及困难。事实上，真正情况要比我们的模型复杂的多，比方，由于初期传播机理不明，医护人员为收到好的保护，很多受到传染。而传播率事实上也是按照*种趋势变化的，他与人们的警觉心理等因素有关。要建立这样的模型，我们认为：1、要不断的研究传染病理，完善模型。2、建立实时

21、反响机制，不断修改模型参数，以求更准确的预测以后的数据。3、尽量以实际的研究数据建模，力求准确反映事件的实质。这样，当人为修改参数进展模拟预测时，可以为采取什么措施效果最好做出较准确的建议。4、研究各种措施实际对疾病传播的影响。要做到这些，也恰恰是建模的难点，因为当涉及到新的传染病的时候，以上数据或参数都必须通过长期的研究或拟合才能得到的。比方措施对疾病传播的影响，是要统计大量实际数据才能获得的。以上也可以看作我们模型的改良方向。7、参考文献1世界卫生组织. ./csr/disease/coronavirus_infections/MERS-CoV-cases-rok-21Jul

22、15.*ls*?ua=12 周宁 林，基于灰色微分方程的传染病预测，.ki./Article/CJFDTotal-LLYY202104016.htm，2021 .07.298、附录1：2：y= dsolve('Dy=81.9185-0.1852*y','y(1)=52','t') y = 819185/1852 - (722881*e*p(-(463*t)/2500)*e*p(463/2500)/1852 >> *=81.9185-0.5*0.1852*(

23、819185/1852 - (722881*e*p(-(463*t)/2500)*e*p(463/2500)/1852)+(819185/1852 - (722881*e*p(-(463*(t-1)/2500)*e*p(463/2500)/1852) >> * >> t=1 2 3 4 5 6 7 8t =     1    

24、; 2     3     4     5     6     7     8>> *=81.9185-0.5*0.1852*(819185/1852 - (722881*e*p(-(463*t)/2500)*e*p(463/2500)/1852)

25、+(819185/1852 - (722881*e*p(-(463*(t-1)/2500)*e*p(463/2500)/1852)* =   79.6419   66.1775   54.9894   45.6928   37.9679   31.5490   26.2152   21.7832y= dsolve(&#

26、39;Dy=81.9185-0.1852*y','y(1)=52','t') y = 819185/1852 - (722881*e*p(-(463*t)/2500)*e*p(463/2500)/1852t=1 2 3 4 5 6 7 8t =     1     2     3

27、     4     5     6     7     8>> *=81.9185-0.5*0.1852*(819185/1852 - (722881*e*p(-(463*t)/2500)*e*p(463/2500)/1852)+(819185/1852 - (722881*e*p(-(463*(t-1)/2500)*e*p(463/2500)/1852)* =   79.6419   66.1775   54.9894   45.6928   37.9679   3