埃博拉病毒传播分析报告与数学建模99569

《埃博拉病毒传播分析报告与数学建模99569》由会员分享，可在线阅读，更多相关《埃博拉病毒传播分析报告与数学建模99569（15页珍藏版）》请在装配图网上搜索。

1、word*大学数学建模竞赛承诺书我们仔细阅读了*大学数学建模竞赛的参赛规如此与竞赛纪律。我们完全明白，在竞赛开始后参赛队员不能以任何方式包括、电子、网上咨询等与队外的任何人研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛纪律的, 如果引用别人的成果或其他公开的资料包括网上查到的资料，必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们重承诺，严格遵守参赛规如此和竞赛纪律，以保证竞赛的公正、公平性。如有违反竞赛纪律的行为，我们将受到严肃处理。我们授权*大学数学建模竞赛组委会，可将们的论文以任何形式进展公开展示包括进展网上公示，在书籍、期刊和其他媒体进展正式或非正式

2、发表等。参赛的题目 (从A/B中选择一项填写)B参 赛 队 员某某学号院系 日期： 2015年 05 月 04日13 / 15埃博拉病毒传播分析摘 要本文的研究对象为1976年在苏丹南部和刚果的埃博拉河地区发现的埃博拉病毒。埃博拉病毒是一种生物安全等级为4级，并且能引起人类和灵长类动物产生埃博拉出血热的烈性传染病病毒，其主要是通过病人的血液、唾液、汗水和分泌物等途径传播。其病毒的潜伏期通常只有5天至10天，感染后25天出现高热，69天死亡。面对其强大的传染力和对人类健康的巨大威胁，本文通过数学建模的方法了解埃博拉病毒的传播规律，并分析隔离措施的严格执行和药物治疗效果的提高等措施对控制疫情的作用

3、。本文中，首先我们根据已给的信息与相关假设数据，通过对条件和所给表格书记的分析，我们大致明白了猩猩从潜伏到发病再到死亡或自愈的过程，因此我们采用了excel拟合曲线，分析其发病、潜伏、自愈、死亡和隔离的相应的变化曲线，估计参数，再根据其建立数学模型，并用MATLAB求解方程组，调试参数，从而得到我们需要的结果。其次通过对已经得到的数据和曲线图的分析，可以得出人类通过严格的药物控制过后，对其发病和潜伏的影响，从而能够达到对疫情的控制的作用，并且对埃博拉病毒未来开展趋势有了更深刻的了解，以为更好的控制埃博拉病毒做出贡献。关键词：非线性曲线拟合；微分方程；MATLAB；数学模型1 问题的重述1.1

4、背景埃博拉病毒又译作伊波拉病毒于1976年在苏丹南部和刚果的埃博拉河地区被发现后，引起了医学界的广泛关注和重视。该病毒是能引起人类和灵长类动物产生埃博拉出血热的烈性传染病病毒，其生物安全等级为4级。埃博拉病毒有传染性，主要是通过病人的血液、唾液、汗水和分泌物等途径传播。各种非人类灵长类动物普遍易感，经肠道、非胃肠道或鼻途径均可造成感染，病毒的潜伏期通常只有5天至10天，感染后25天出现高热，69天死亡。发病后14天直至死亡，血液都含有病毒。埃博拉病毒感染者有很高的死亡率在50%至90%之间，致死原因主要为中风、心肌堵塞、低血容量休克或多发性器官衰竭。当前主流的认知是，埃博拉病毒主要通过接触传播

5、，而非通过空气传播；只有病人在出现埃博拉症状以后才具有传染性。在疾病的早期阶段，埃博拉病毒可能不具有高度的传染性，在此期间接触病人甚至可能不会受感染，随着疾病的进展，病人的因腹泻、呕吐和出血所排出的体液将具有高度的生物危险性；存在似乎天生就对埃博拉免疫的人，痊愈之后的人也会对入侵他们的那种埃博拉病毒有了免疫能力。埃博拉病毒很难铲除，迄今为止已有屡次疫情爆发的记录。据百度百科，最近的一次在2014年。截至2014年9月25日，此次在西非爆发的埃博拉疫情已经导致逾3000人死亡，另有6500被确诊感染。更为可怕的是，埃博拉病毒可能经过变异后可以通过呼吸传播！1.2 问题假设某地区有20万居民和30

6、00只猩猩。人能以一定的概率接触到所有的猩猩，当接触到有传播能力的猩猩后有一定概率感染病毒，而人发病之后与猩猩的接触可以忽略。研究人员统计了前40周人类和猩猩的发病数量和死亡数量等信息，请你根据相关信息，研究回答以下问题：1、 根据猩猩的发病数量和死亡数量，建立一个病毒传播模型，动态描述病毒在“虚拟猩猩种群中的传播，并预测接下来的在猩猩中的疫情变化，并以下述格式给出“虚拟猩猩种群在第80周、第120周、第200周的相关数据；“虚拟猩猩种群群体数量预测结果单位：只潜伏群体处于发病状态累计自愈累计因病死亡第80周第120周第200周2、 建立“虚拟种群相互感染的疾病传播模型，综合描述人和猩猩疫情的

7、开展，并预测接下来疫情在这两个群体中的开展情况，并以下述格式给出 “虚拟人类种群在第80周、第120周、第200周的相关数据；“虚拟人类种群群体数量预测结果单位：个潜伏人群处于发病状态隔离治疗累计治愈累计因病死亡第80周第120周第200周3、假设在第41周，外界的专家开始介入，并立即严格控制了人类与猩猩的接触，且通过某种特效药物将隔离治疗人群的治愈率提高到了80%。请预测接下来疫情在“虚拟人类种群的开展情况，比照第2问的预测结果说明其作用和影响，给出“虚拟人类种群在第45周、第50周、第55周的相关数据，数据格式同问题2；4、请依据前述数学模型，分析各种疫情控制措施的严格执行和药物包括防疫药

8、物、检疫药物和治疗药物等效果的提高等措施对控制疫情的作用。2 问题分析通过对条件的分析，并通过给出的表格数据，大致明白猩猩从潜伏到发病再到死亡或自愈。我们通过excel作出发病随时间的变化曲线，潜伏随时间变化曲线，估计参数。然后通过建立数学模型用MATLAB解出方程组，调试参数使其死亡，自愈等曲线与给出表格大致一样，然后通过建立的模型求出问题一。 同问题一分析，我们通过excel作出相应处于发病状态的曲线，自愈以与死亡和隔离的曲线，估计模型相应的参数。然后通过建立的数学模型用MATLAB解出方程组，调试参数使其自愈，处于发病等曲线和表格给出的数据大致一致。 同问题二分析，我们通过excel作出

9、治愈率提高80%后相应处于发病状态的曲线，自愈以与死亡和隔离的曲线，估计模型相应的参数。然后通过建立的数学模型用MATLAB解出方程组，调试参数使其自愈，处于发病等曲线和表格给出的数据大致一致。 通过上术数据和曲线图的分析，可以很清楚的看出当有人类干预后即就是严格的通过药物后，发病和潜伏等都有很明显的改善。3 假设与符号3.1模型的假设：n 由于埃博拉病毒的传播期限不是很长，故假设不考虑这段时间的人口出生率和自然死亡率；n 平均潜伏期限为6天；n 处于潜伏期的埃博拉病人不具有传染性。3.2符号说明：t0 表示从最初发现埃博拉患者到卫生部门采取预防措施的时间间隔；N 表示疫区总人口数；S(t)

10、表示t时刻健康人数占总人口数的比例；I(t) 表示t时刻感染人数占总人口数的比例；E(t) 表示t时刻潜伏期的人口数占总人口数的比例；Q(t) 表示t时刻退出类的人数占总人数的比例；(t) 表示日接触率，即表示每个病人平均每天有效接触的人数；N 表示疫区总猩猩口数；S(t) 表示t时刻健康猩猩数占总猩猩数的比例；I(t) 表示t时刻感染猩猩数占总猩猩数的比例；E(t) 表示t时刻潜伏期的猩猩数占总猩猩数的比例；Q(t) 表示t时刻退出类的猩猩数占总猩猩的比例；(t) 表示日接触率，即表示每个病猩猩平均每天有效接触的猩猩数；(t) 表示日接触率，即表示每个病猩猩平均每天有效接触的人数；g(t)

11、表示政府控制力度；f(t) 表示疫情指标。4 模型的建立与求解 由问题的分析，将猩猩群分为易感猩猩群S，病毒潜伏猩猩群E，发病猩猩群I，退出者Q四类：l 易感人群S与病毒潜伏人群E之间的转化易感者和发病者有效接触后成为病毒潜伏者，设每个病人平均每天有效接触的健康人数为(t)S，NI个病人平均每天能使(t)SNI个易感者成为病毒潜伏者。故,即l 病毒潜伏人群E与发病人群I间的转化潜伏人群的变化等于易感人群转入的数量减去转为发病人群的数量，即。l 发病人群I与退出者Q间的转化 单位时间退出者的变化等于发病人群的减少，即 很明显从我们建立的模型是无法得到E,S,I,Q的解析解的。为了解决这个问题，我

12、们求助于计算机软件MATLAB来求出它们的数值解。我们先通过附件中给的数据算出每一天的E,S,I,Q,做出它们与时间的函数图象，然后画出我们通过模型解出的数值解随时间变化的图象。比照这两组图，可以发现实际和理论存在着一定的差异。这必然是因为我们的参数估计不合理造成的。所以，我们必须通过不断调整那些非计算得到的参数，来使实际图象和理论图象趋于一致。经过屡次调试,我们发现，当=0.680人,=0.9,=0.58时,实际图象和理论图象有最好的符合。而这三个值均在我们估计的围，所以我们认为这三个值的得到是合理的。一旦参数确定，就可以通过MATLAB软件求出该方程组在某个区间段的数值解，从而可推算出我们

13、所需的数值如下表所示。周数SEQ第80周第120周第200周0在根据逻辑关系式计算可得下表的预测值 表1 “虚拟猩猩种群群体数量预测结果单位：只周数潜伏人群处于发病状态累计治愈累计因病死亡第80周30283596第120周00299598第200周00300600结果分析 根据上表可知，在第80周以后，处于潜伏状态的猩猩接近于0 ，处于发病状态的猩猩也趋近与0，且猩猩的治愈数和因病死亡数变化不大，由该模型预测出的结果与附件中的数据的得出的发病率和累计死亡率趋势一样。 健康人数占总数比例 比对图1.1 健康人数占总数比例图参考数据图1.2 健康人数占总数的比例图模拟数据潜伏人数占总数比例比对图2

14、.1 潜伏人数占总数的比例图参考数据图2.2 潜伏人数占总数的比例图模拟数据退出人数占总数比例比对图3.1 退出人数占总数的比例图参考数据图3.2 退出人数占总数的比例图模拟数据MATLAB主要程序functiondx=rossler(t,x,flag,a,b,c)dx=-a*x(1)+a*x(1)*x(3)+a*x(1)*x(2)+a*x(1)*x(1);a*x(1)-a*x(1)*x(3)-a*x(1)*x(2)-a*x(1)*x(1)-b*x(2);c-c*x(3)-c*x(2)-c*x(1);a=0.680;b=0.90;c=0.580;x0=0.9950.0050;t,y=ode45

15、(rossler,080,x0,a,b,c);flot(t,y); 由问题的分析，将人群分为易感人群S，病毒潜伏人群E，发病人群I，退出者Q四类：l 易感人群S与病毒潜伏人群E之间的转化易感者和发病者有效接触后成为病毒潜伏者，设每个病人平均每天有效接触的健康人数为(t)S，NI个病人平均每天能使(t)SNI个易感者成为病毒潜伏者。故,即l 病毒潜伏人群E与发病人群I间的转化潜伏人群的变化等于易感人群转入的数量减去转为发病人群的数量，即。l 发病人群I与退出者Q间的转化 单位时间退出者的变化等于发病人群的减少，即 很明显从我们建立的模型是无法得到E,S,I,Q的解析解的。为了解决这个问题，我们求

16、助于计算机软件MATLAB来求出它们的数值解。 我们先通过附件中给的数据算出每一天的E,S,I,Q,做出它们与时间的函数图象，然后画出我们通过模型解出的数值解随时间变化的图象。比照这两组图，可以发现实际和理论存在着一定的差异。这必然是因为我们的参数估计不合理造成的。所以，我们必须通过不断调整那些非计算得到的参数，来使实际图象和理论图象趋于一致。隔离治疗人数占总人数的比例图4.1 隔离治疗人数占总人数的比例图模拟数据图4.2 隔离治疗人数占总人数的比例图参考数据死亡人数占总人数的比例图5.1 死亡人数占总数的饿比例图模拟数据图5.2 死亡人数占总数的比例图参考数据自愈人数占总人数的比例图6.1

17、自愈人数占总数的比例图模拟数据图6.2 自愈人数占总人数的比例参考数据发病人数占总数的比例图图7.1 发病人数占总人数的比例参考数据图7.2 发病人数占总数的比例图模拟数据表2 “虚拟人类种群群体数量预测结果 单位：个潜伏人群处于发病状态隔离治疗累计治愈累计因病死亡第80周6547610002370第120周7250814504980第200周5946521007650结果分析 ：由上表可知，在第80周以后，处于潜伏状态的人群变化幅度不大，处于发病状态的人群也变化幅度不大，且人群的治愈数和因病死亡数持续增长，由该模型预测出的结果与附件中的数据的得出的发病率和累计死亡率趋势一样。外界的专家开始介

18、入，并立即严格控制了人类与猩猩的接触，且通过某种特效药物将隔离治疗人群的治愈率提高到了80%。专家的预防措施力度g(t)在控制疫情的过程中起到了重要的作用，与如下因素有关：l 专家关注的疫情来自于最近几天的疫情，不妨取近三天的平均值；l 当t=t0时，g(t)有一个初始值，即为潜在的政府力度K0；综上所述，可以给出g(t)随疫情变化的曲线，形态如下列图，横坐标为疫情，纵坐标为g(t)，其表达式为其中k0+k1=1,。根据有关数据，令k0=0.2，k1=0.8，当=0.58时，取g(t0)=0.7，得参数估计=0.1803.政府控制力度g(t)与日传染率(t)的关系：(1)当政府控制力度为0的时

19、候(t)取最大值；(2)随着g(t)的增大，(t)减小；(3)当g(t)不强时，对(t)的变化所起的作用较小；(4)当g(t)超过一定的数值时对(t)的影响效果明显；(5)当g(t)趋近于1的时候(不可能为1)，如此(t)趋近0。由以上几点可以确定(t)随g(t)的变化关系曲线，采用函数刻画此形态，其中为常数。图8 g(t)和(t)的关系图 表3 “虚拟人类种群群体数量预测结果 单位：个潜伏人群处于发病状态隔离治疗累计治愈累计因病死亡第45周6334127151786第50周5827118241825第55周4613149581876结果分析 由上表可知，在专家介入后，埃博拉病毒的预防控制力度

20、加大，累计治愈的人数在增多，因病死亡人数虽然在增加，但是其增加幅度不大，说明埃博拉病毒已经得到了良好的控制，与预期估测结果相吻合。在发病初期，由于人们对埃博拉病毒的认识不够，重视不足，防措施较差，没有有效的防疫药物、检疫药物和治疗药物治疗，也没有相应的政府控制措施，随着时间t的增长，病情不断恶化，感染病情所占比例I呈现不断增加的趋势，健康人数占总人数的比例S不断下降，退出率Q也呈现持续增长的趋势，造成了巨大的经济损失和人员伤亡。在发病中后期，随着相关政府的介入和对该病毒的相关知识的普与，提高了人们对埃博拉病毒的预防意识，同时，在科研人员的不断努力下，防疫药物、检疫药物和治疗药物的种类增多、疗效

21、增强，随着时间的增长，感染患者的比例I呈下降趋势，健康人数所占比例S的下降趋势由急变缓，治愈率不断提高，死亡人数得到控制，一场殃与全人类的疫情风波得到较好的控制。5 模型的评价本模型中，我们根据已给的信息与相关假设数据，通过对条件和所给表格书记的分析，我们大致明白了猩猩从潜伏到发病再到死亡或自愈的过程，因此我们采用了excel拟合曲线，分析其发病、潜伏、自愈、死亡和隔离的相应的变化曲线，估计参数，再根据其建立数学模型，并用MATLAB求解方程组，调试参数，从而得到我们需要的结果。其次通过对已经得到的数据和曲线图的分析，可以得出人类通过严格的药物控制过后，对其发病和潜伏的影响，从而能够达到对疫情

22、的控制的作用，并且对埃博拉病毒未来开展趋势有了更深刻的了解，以为更好的控制埃博拉病毒做出贡献。 本模型重点是分析规律和进展预测。因为数据受很多随机因素的影响，规律性受到干扰，所以其变化情况不能较好地表达总体的规律性，进而不能对疫情进展较准确的预测；针对这个问题，我们对数据进展了统计平均，从总体的平均规律入手，没有局限于仅对现有数据的模拟。但是也要根据现有的数据对模型进展检验。从前面求解方程得到的图形结果来看，模拟的曲线确实较好地代表了现有数据的总体变化规律。不论是本论文模型还是概率模型，进一步的工作和更准确的结果给出将有待于收集传染病学实际资料。相信随着人们对埃博拉的进一步认识，随着社会各界的深入研究，从数学角度看，其传播模型将更加完善，预测结果将更准确，从医学角度看，埃博拉将有更好的治疗方案和防控措施，疫期将进一步缩短。参考文献1 秀兰 林峰. 数学建模于实验化学工业. 20132 贺超英 王少喻. MATLAB应用与实验教程. 20133 德丰. MATLAB实用数值分析. 200124 SARS数学建模优秀论文.5 学文 炳照 王宏洲 数学建模优秀论文精选与点评.2011