埃博拉病毒传播模型

埃博拉病毒传播模型

【范文精选】埃博拉病毒传播模型

【范文大全】埃博拉病毒传播模型

【专家解析】埃博拉病毒传播模型

【优秀范文】埃博拉病毒传播模型

范文一:埃博拉病毒如何传播

健 康 社 区 健康新知

英 国玛 丽女 王 大 学 癌 症 预 防 中心 教 授 祖

埃 博 拉 病 毒 如 何 传 播

目前 的研 究 表 明 , 埃 博 拉 病 毒 主 要 是靠 感 染 者 的  体液传播 。这些体液包 括感染者 的唾液 、 呕吐物 、 血

液、 排泄物、 l I E N、 汗水 和 母乳 等 ( 有 些专 家 对 部分 体 液

兹克在英 国   肿瘤学年报  杂志发表报告 : 中老  年人 如 果 长 期 小剂 量服 用 阿司 匹林 , 可 以显著  降低 罹患某 些癌症 的风 险或 者降低 死 亡率 。

祖 兹克领 导 的 团队通 过分析 超 过 2 0 0 组 阿

司 匹林 抗 癌效 果研 究和 临床 试验 发 现 , 如果 1 0

是 否会 传 播 有 不 同 意 见 ) 。如 果 感 染 者 身 处 晚 期 , 那

么, 他们 一丁 点体 液都 携带着 大量 病 毒 。   过往研究表 明, 如 果 医 护 人 员 不 慎 接 触 到 埃 博 拉  患者 的 血液 , 只 要 自己身 上没 有 伤 口, 这 些 血液 也 未 触

年 连 续服 用 阿 司 匹林 , 肠 癌发 病率可 降低 大约

3 5 %, 死亡 率 减少 4 0 %; 食 道 癌和 胃癌发 病率 可

降低 3 O %, 死亡率减 少3 5 % ~5 0 %。 如 果服 用

阿司 匹林 至 少 l 0 年, 这一 年龄 段 中, 男性 患癌 、   中风 和 心脏 病 的 几 率 可 降 低 9 %, 女 性发 病 率

降低 大约 7 %。

碰 到 自己 的嘴 巴 、 鼻子 和 眼睛 等 , 用 肥皂 和 清水 把 血 液

冲 洗干 净 即能 防止感 染 。   跟 中 国人谈 虎 色变 的 S A R S ( 非典) 不同, 埃 博 拉 病

毒 目前 尚不 能通 过 空 气传 播 。 不过 , 医学 专家 认 为 , 从

为达到 防癌 抗癌 效 果 , 5 O~6 5岁的人 群 需  要5 ~1 0 年 坚持 服 用 阿 司 匹林 , 每 日7 5 ~1 0 0   毫克 。祖 兹克 教授 指 出, 头3 年 看 不见效 果 , 服  用5 年后 才 能 降低 死亡 率 。但 如 果服 用 阿 司 匹

林 1 O年 , 消 化 道 出血 的 风 险 可 从 2 . 2 %增 至

今 年 的 现 实来 看 , 埃 博拉 是 高传 染 性病 毒 , 传 染性 要 强

过S A R S 、 腮腺炎、 麻疹等传染性疾病。 目 前, 紫外灯、 干  燥、 高温 、 消毒 液 ( 肥 皂水 或 酒精 ) 都可 轻 易杀 死埃 博 拉

病 毒 。病 毒离 开 动 物 体 内 、 在 阴冷 潮 湿 环 境 下 的存 活

期, 只 有短 短几 天 。

3 . 6 %, 对 少数 人 群 会 构 成致 命 危 险 。 因此 , 大  量服 用 阿 司 匹林有 风 险 , 人 们在 开 始每 天服 用

埃 博 拉病 毒 的潜 伏 期 一 般 是 5 ~ 7 天, 最短 2 天, 最

前应该向医生寻求建议 。

长可 达 2 1 天 。初期 症 状是 发烧 头痛 、 关 节 和肌 肉酸 痛 、

咽 喉疼 痛 和 肌 肉无 力 。这 跟 患 流 感 的 症状 很 相 似 , 因  此, 最 开 始 一 般 很 难 诊 断 某 人 是 否感 染 了埃 博 拉 。 目

前医护人 员只能通过询 问患者是否 到过西非“ 疫 区”

等, 进 而做 进 一步检 查 。

英 真扩大袋 - 汀类药  ’ 连更蔻图

英 国卫 生部 门2 0 1 4 年7 月更新的一项健康 指导 ,   将降胆 固醇的他汀类药物适用 范围进 一步扩大。根据

新 研 究认 为  糖 比 盐 更伤 血 压

长期 以来 , 多 项研 究 告 诫 : 吃 盐 过 多会 引发 高血

新指导意见, 未来 1 O 年 内患病风 险在 2 0 %以上的英国

人都可服用此类药物。他汀类药物可帮 助减少脂肪斑  块在血管 中的堆积, 降低血管堵塞 引发心脏病、 中风等   疾 病 的风 险。   根据这一建议, 英 国6 0 岁以上 的男性和6 5 岁 以上  的女性均可从 医生处获得他 汀类药物, 而在年轻人 中,   胆 固醇高、 肥胖 、 有心脏病家族史及吸烟人群也应服药  以降低 患病 风 险。   英国国家 卫生与临床优化研 究所教授马克 ・ 贝克  说’ 防病总 比治病要好, 疾病预防正是现代 医学的最重   要支柱之 一。不过也有反对意见认为, 关于他 汀类 药   物的副作用 目前尚无定论 , 不应 盲 目 扩大其适用 范围,   改变不 良生活习惯仍是降低心血管疾病风险最重要的

手段 。

长期 服用阿司 匹林 可降低患癌 风险

压, 使 中风危 险增加 1 / 4 , 每年 导致 全球 3 0 0 万 人 死  亡 。 日前 , ( ( 美 国 心脏 病 杂 志 刊登 一 项 新 研 究报 告 认  为, 导致 高血 压 的 罪魁 祸 首 不是 盐 而是糖 。

美 国研 究人 员认 为 , 高血 糖 水 平 影 响 大脑 的 关键

区域 下 丘脑 , 进 而导 致 心率加 快 和血 压 升 高。 另外 , 吃

糖 太多导致血糖升 高进 而刺激人体 分泌更多胰 岛素 ,   而胰岛素 同样会刺激心率加快。

美国“ 控糖行动组织” 科 学 总监 阿斯 姆 ・ 马 尔霍 特  拉 博 士表 示 , 糖 对 公 众 健 康 构 成 的风 险被 低 估 了 。控  糖 应 和控 盐一样 成 为预 防 高血 压 的重要 策略 之一 。

5 2

范文二:埃博拉病毒的传播

西安工业大学数学建模竞赛

承 诺 书

我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.

我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮

件、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问

题。

我们知道,抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他

公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正

文引用处和参考文献中明确列出。

我们郑重承诺,严格遵守竞赛规则,以保证竞赛的公正、公平性。如有违反

竞赛规则的行为,我们将受到严肃处理。

我们参赛选择的题号是(从A/B/C中选择一项填写): B

所属学校(请填写完整的全名): 西安工业大学 参赛队员 (打印并签名) :1. 陈文兴

2. 闫丽萍

3. 魏栩

指导教师或指导教师组负责人 (打印并签名):

日期: 2015 年 8 月 1 日

埃博拉病毒传播及控制分析

摘要

埃博拉病毒是能引起人类和灵长类动物产生埃博拉出血热的烈性传染病病毒,有很高的死亡率。本文根据研究人员统计所给出的前四十周人类和猩猩的发病数量和死亡数量等信息,对该病毒的传播、预测与控制进行研究并建立模型,并分析了隔离措施的严格执行和药物治疗效果的提高等措施对控制疫情的作用。 针对问题一,在了解埃博拉病毒的传播情况后,根据猩猩的发病情况建立了马尔萨斯模型:xt154.97e0.0270t。在此模型中,较好地描述病毒在“虚拟猩猩种群”中的传播情况;根据“虚拟猩猩种群”中的数据,用matlab拟合出不同状态下猩猩数量的变化曲线,并以发病状态为例建立灰色预测模型dx1

 0.0669x160.0539dt,从而较准确的预测出接下来第80、xk1(x(0)(1)b)eab947e0.0124897aa

120、200周的猩猩发病状态的数据。

针对问题二,为描述埃博拉病毒在“虚拟种群“中的相互传播规律及人和猩猩的疫情发展状况,建立SEIR模型 dQdtrI(t)11a2a1etQ(t)ttdIdt11a2a1eQ(t)11a3eI(t)

tdTdt11a3eI(t)

模型求解时,通过对模型的推导,我们发现不能给出每个函数的解析解,因此考虑利用matlab中的ode45函数进行求解。得出了患者数量随时间的变化规律。同样利用灰色预测模型预测出“虚拟人类种群”在第80、120、200周的相关数据。

针对问题三,在问题已建立的模型之上画图分析两个因素:通过某种特效药改变治愈率到80%,控制患者和健康人群的接触即控制隔离强度。对埃博拉病毒传播的影响,并通过图中控制后的患者人数,利用模型二中的关系表达式,计算出,45,50,55周的潜伏期人数,治愈人数,死亡人数。

针对问题四,对问题三进一步讨论改变隔离强度和治愈率对病毒传播的影响,分别用matlab作出患者人数随时间的变化曲线,对比分析,可得出:当降低患者与健康者的接触率和使用特殊药物提高治愈率时,随着时间的延迟,患者人数急剧下降。所以实际生活中,改变通过这两个指标可以有效的控制病情的传播。

关键词:马尔萨斯模型 SEIR模型 灰色预测 隔离强度

一、问题重述

埃博拉病毒是能引起人类和灵长类动物产生埃博拉出血热的烈性传染病病毒,埃博拉病毒感染者有很高的死亡率(在50%至90%之间),主要通过接触而非空气传播。迄今为止,已有多次疫情爆发的记录,最近的一次在2014年,截至2014年9月25日,此次在西非爆发的埃博拉疫情已经导致逾3000人死亡,另有6500被确诊为感染。

本文假设某地区有20万居民和3000只猩猩。研究人员统计了前40周人类和猩猩的发病数量和死亡数量等信息(见附件一、附件二),根据相关信息研究回答以下问题:

1、根据猩猩的发病数量和死亡数量,建立病毒传播模型,动态描述病毒在“虚拟猩猩种群”中的传播,并预测猩猩接下来的疫情变化,并给出“虚拟猩猩种群”在第80周、第120周、第200周的相关数据;

2、建立“虚拟种群”相互感染的疾病传播模型,综合描述人和猩猩疫情的发展,并预测接下来疫情在这两个群体中的发展情况,并给出 “虚拟人类种群”在第80周、第120周、第200周的相关数据;

3、由于41周外界专家的介入,及严格控制了人类与猩猩的接触,且通过某种特效药物将隔离治疗人群的治愈率提高到了80%。预测接下来疫情在“虚拟人类种群”的发展情况,并对比第2问的预测结果说明其作用和影响,给出“虚拟人类种群”在第45周、第50周、第55周的相关数据;

4、依据前述数学模型,分析各种疫情控制措施的严格执行和药物(包括防疫药物、检疫药物和治疗药物等)效果的提高等措施对控制疫情的作用。

二、问题分析

根据题意,这是一个传染性病毒随时间蔓延的过程,需要研究传染病在传播过程中各类种群的数量变化,特别是通过研究患者和疑似患者、患者和治愈者的数量变化,预测传染病的传染的高峰期和持续时间长度,从而我们可以认识到相应的疫情控制措施对控制传染病传播所能达到的效果。

针对问题一,根据附件一“虚拟猩猩种群”的数据,初步观察到发病状态、累计自愈及累计死亡的猩猩数量,并对数据做定量分析得到截至每周累计发病的数量,利用matlab编程得出病毒传播速度的散点图,针对病毒的传播过程,首先,我们用xt表示t时刻的猩猩发病个数,用表示每天每个猩猩有效接触的个数,考虑t到tt时刻发病个数的增加,建立微分方程 dxx,x0x0,dt

通过马尔萨斯模型求解得:xtx0et。利用灰色预测,预测出后期猩猩每个状态的数量,并以发病状态为例建立灰色预测模型,接着运用matlab编程假设

X0k50.2,54.4,51.2,48.4,45.4,42.4,39.2,36,再对其作一次性累加生成运算得到新的生成数列X1k50.2,104.6紧接着对,,155.8,204.2,249.6,292,331.2,367.2

X1k作紧邻均值生成得出数据阵B和数据向量Yn,再对参数列[a,b]T进行

最小二乘估计最后建立出了灰色模型(GM(1,1)模型)。我们又经过对GM(1,1)模型的残差检验,最终得出了预测结果。

针对问题二,我们把人群分为五类:健康者(易感染者)、确诊患者、疑似患者、潜伏期感染者、死亡者和治愈者,采用SIER模型,并将死亡者和治愈者都归于系统移出者统称为恢复人群,关系如下图:

在此基础上,找出单位周内这五类群体数量的变化来建立微分方程,得出模型。再利用matlab编程画出图形,改变其隔离强度药物治愈率后重新作图进行比较,对结果进行分析,并利用SEIR模型对埃博拉病毒的传播规律进行定性分析描述,对未来种群数量变化用灰色预测模型进行定量预测,并分析各种疫情控制措施对控制疫情的作用。

针对问题三,在问题已建立的模型之上画图分析两个因素:通过某种特效药改变治愈率,控制患者和健康人群的接触即控制隔离强度。对埃博拉病毒传播的影响。

针对问题四,对问题三进一步讨论改变隔离强度和治愈率对病毒传播的影响,分别用matlab作出患者人数随时间的变化曲线,对比分析,给出有效控制疫情传播的建议。

三、符号说明

符号

S(t)

E(t)

Q(t)

I(t)

T(t) 解释说明 t时刻正常者(易受感染)数量 t时刻疑似患者的数量 t时刻处于潜伏期的数量 t时刻确诊患者的数量 t时刻退出传染系统的数量(包括治愈者和死亡者)

β1

β2

a3

r 潜伏期的人数中转化为确诊患病的数量占潜伏期数量的比例 每日退出传染系统的数量比例 确诊患者的治愈时间 患者的人均周接触个体数量

因接触被感染的概率

潜伏期内的患者被隔离的强度  p

a 患者被治愈的概率

四、模型假设

1、假设单位时间内感染病毒的数量与现有的感染者成比例;

2、假设单位时间内治愈数量与现有感染者成比例;

3、假设单位时间内死亡数量与现有的感染者成比例;

4、假设患者治愈恢复后不会再被感染同种病毒,有很强的免疫能力,即被移除出此传染系统;

5、假设正常者被传染后,进入一段时间的潜伏期,处于潜伏期的群体不会表现症状,不可传染健康者,不具有传染性;

6、假设患者入院即表示患者被隔离治疗,被视为无法跟别人接触,故不会传染健康者;

7、假设实际治愈周期过后,如果患者没有治愈,则认为患者死亡,即实际治愈周期过后,患者都被移出此感染系统;

8、假设考察地区内疾病传播期间忽略个体的出生,死亡,流动等种群动力因素对总种群数量的影响。即:总种群数量不变,记为N=20万;

9、假设人能以一定的概率接触到猩猩,当接触有传播能力的猩猩后有一定的概率感染病毒,而人发病后与猩猩的接触可以忽略。

10、假设可以及时发现疑似患者并隔离治疗,并且剩下一部分未被隔离的感染者变成患者后感染正常人。

11、假设人类发病后与猩猩的接触可以忽略,即人类不作为猩猩的感染源。

五、模型的建立与求解

5.1问题一模型建立与求解

据问题一,由附件“虚拟猩猩种群”中的数据,利用Matlab作累计发病个数--时间的散点图如下:(程序参见附件1)

从图可看出:埃博拉病毒传播的速度在前40周始终呈上升的趋势,但上升的斜率有减小的趋势。

5.1.1马尔萨斯模型的建立与求解

“虚拟猩猩种群”的埃博拉病毒传播预测模型类似于人口增长的预测模型,故首先采用马尔萨斯模型(Malthusian 模型)进行建模。设时刻t的确诊患者数量It是连续、可微函数,并且每天每个患者有效接触(足以使人致病的接触)的人数为常数,考察t到t+t病人人数的增加,则有

IttItItt

再设t0时有I0个病人,即得微分方程 dII dtI0I0

解之可得:

ItI0et

其中, I0, 为常数。It154.97e0.0270t

根据“猩猩种群”疫情数据中的确诊患者的数据散点图(图1),考虑利用马尔萨斯模型xtx0et来预测埃博拉病毒的传播情况。用matlab求得

: I0154.97,0.0270。即得马尔萨斯模型如下(程序参见附件2)

I

图:

t154.97e0.0270t

结果表明,随着t的增加,猩猩感染的个数It持续增长。马尔萨斯拟合及预测图线与猩猩在前40周发病情况图线拟合程度较为符合。

由图分析知:马尔萨斯模型是关于人口或种群增长的模型,它发现人口或种群成指数增长。即在该模型中可引意为,猩猩感染病毒的个数随着时间的增长呈指数增长变化。但现实生活中,由于猩猩数量的有限,所以该模型在预测短期的数据上是可以的,但对未来长期预测情况跟实际显然是不太相符合的,因此暂不考虑用该模型进行数据预测。因此对问题一中数据的预测采用灰色预测模型。

5.1.2 利用灰色系统预测结果

我们建立的马尔萨斯较为合理的模型仅仅只是对埃博拉病毒的传播规律进了动态描述,而不能定量的进行数据预测。因此,我们先将附件中“虚拟猩猩种群”的数据用excel进行了处理:分别计算出附件中前40周猩猩的潜伏群体、处于发病状态、累计自愈、累计因病死亡的数量,然后在matlab中做出上述四个状态下猩猩数量关于周数的曲线拟合。(程序参见附件3)

1、潜伏群体

2、处于病状态

3、累计自愈

4、累计因病死亡

为了预测接下来80周,120周,200周的情况,我们分别建立灰色预测模型来预测四种状态的情况,由于模型一致,本文以猩猩发病状态为例建立模型,首先取均值将前40周分为八组(具体情况见附件excel表格)。

得出如下表(单位:只):

假设X0k50.2,54.4,51.2,48.4,45.4,42.4,39.2,36.

5.1.2.1 GM(1,1)模型的建立

为了使其成为有规律的时间序列数据,对其作一次累加生成运算,即令

X(1)(t)X(0)(k)

n1t(t,k1,27)

从而得到新的生成数列X1k50.2,104.6.对,155.8,204.2,249.6,292,331.2,367.2X1k做紧邻均值生成. 则数据阵B和数据向量Yn为

121

2121 B2

1

2

121

2

x1x211x2x311x3x411x4x511x5x611x6x711x7x81111x02104.6-77.41 0x3155.81-130.2 1x04-180 1204.21-226.9 1 Ynx05249.6 0 1x6292- 270.81-311.60331.2 1x7 1x08367.2-349.211对参数列[a,b]T进行最小二乘估计,可得

BBT10.00250.00160.0007-0.0001-0.0009-0.0016-0.0022BT 0.69200.48990.29930.1198-0.0482-0.2044-0.3483

(其中,a为发展系数,反映x的发展BTBBTYnb60.05391

a

 0.0669

趋势;b为灰色作用量,反映数据间的变化关系. )

从而可得出GM(1,1)模型:

dx1

 0.0669x160.0539dt  xk1(x(0)(1)b)eab94e70.0124897aa

bb其中,xk1(x(0)(1))ea947e0.0669897为时间响应函数形式。 aa

5.1.2.2 GM(1,1)模型的残差检验

残差大小检验,即对模型值和实际值的残差进行逐点检验. 

(1)根据预测公式,计算X

X11k,得 k50.2,104.6,155.8,204.2,249.6,292,331.2,367.2 k1,2,8

0

(2)累减生成X

k序列,k1,28

X0k50.2,54.8434,51.2963,47.9786,44.8755,41.9731,39.2584,36.7193

原始序列: X0k50.2,54.4, 51.2,48.4,45.4,42.4,39.2,36

(3)计算绝对残差和相对残差序列

绝对残差序列:0 0,0.0082,0 .0019,0.0087,0.0116,0.0101,0.0015,0.02 相对残差: 0,0.015%,0.0037%,00179%,0.0255%,0.0238%,0.0038%,0.055%GM(1,1)模型的残差检验结果:相对残差不超过0.05%,精确度高。

5.1.2.3GM(1,1)模型求解

利用Matlab进行预测,分别取预测的个数为8,16,32(程序参见附件4),得到实际值与预测值如下表(单位:只):

本文用75-80,115-120,195-200均值近似代替第80,120,200周的情况,问题一的结果预测见下表:

5.2 问题二模型的建立与求解

5.2.1 SEIR模型建立与求解

针对问题二,我们充分利用附件一、二,并注意到埃博拉病毒传播过程中每一个群体都处于动态的变化中。对S来说,一部分未被隔离的潜伏期感染者能感染正常者,使其成为潜伏期感染者流出S;对于E来说,流入者包括一部分潜伏期的感染者和一部分正常者,流出者包括一部分没有被感染的正常者和隔离后被确诊患者;对于I来说,它既有从包括隔离和未被隔离的H中确诊的流入者,也有已经治愈的流出者;对于R来说,它只有从I中治愈转化而来的流入者。以上过程在传染的每一时刻都是相同的。在某一时刻对S、E、I、R取其对时间的微分,假设患者可以随意接触和感染正常人。分析t阶段时间内,疫情的发展与变化,这样既可建立传染病传播模型的微分方程组如下:

(1)正常人-----疑似患者:

病人尚未被隔离,所以疫情发展比较迅速,此时病人人均每天接触r个正常人,假设t时刻病人人数为I(t),则新增疑似患者人数为E,

EI(t)rtrI(t)t。

(2)疑似患者-----潜伏期:

疑似患者中包括病毒携带者和非病毒携带者,病毒携带者会进入潜伏期,而非病毒携带者最终还是正常人。

设疑似患者中病毒携带者占疑似患者的比例为,假设t时刻疑似患者人数为E(t),潜伏期患者人数为Q(t),则Q(t)E(t),故新增潜伏期人数为

QE。

(3)潜伏期-----确诊患者:

因为每日潜伏期病人变为确诊患者的数量呈指数增长,用1表示这一特性。那么新增确诊患者人数为I1Q(t)t,现在要确定1,如果潜伏期天数为a1到a2,假设其变化到了一个稳定阶段,那么随着天数的增加潜伏期的病人越来越多,其概率分布呈指数稳步增长,则每天有11a2a1et概率的人变为埃博拉病毒患者,即111a2a1et。所以新增患者人数:

I11a2a1etQt。 (4)确诊患者-----治愈、死亡:

设T为退出系统人数(治愈者和死亡者),如果治愈天数设为a3,那么a3天后病人要么死亡要么被治愈,而被治愈的人产生抗体,不再会被传染,所以被治愈的人和死亡的人都算作退出系统的人。设系统退出率为a3,则有退出人数



TI(t)2t。2的求解方法与1相同,即随着天数的增加退出传染系统的人数也越来越多,则211a3et。故新退出传染系统的人数

T11a3etIt。

(1)~(4) 根据上述的式子可进一步得出:



QEt

Q(tt)Q(t)rI(t)t11a2a1eQ(t)ttt

所以得出以下:

dtrI(t)11a2a1etQ(t)tt

dIdt11a2a1eQ(t)11a3eI(t)t

I(t)dTdt11ae3 

5.2.2 SEIR模型求解

通过对模型的推导,我们发现不能给出每个函数的解析解,因此考虑利用matlab中的ode45函数进行求解。

首先,对传染病模型进行标准化,再带入参数,并由此建立微分方程组函数M文件,随后用ode函数对该文件进行调用,即可得到微分方程组的解向量,然后利用plot函数画出此解向量即可得到虚拟种群患者随周数变化的曲线图如下(程序详见附件5)

由上图分析可知:人类和猩猩患者的数量持续增多,由于人与猩猩数量差别较大,对20万人类来说,随着人类和猩猩发病者增多,感染源群体增多,易感人群发病机率增大,则人类发病数量增加较快;考虑到人发病之后与猩猩的接触可以忽略,而对于3000只猩猩来说,死亡和自愈之后立即退出传染系统,则在总量一定的情况下随着猩猩的传染源减少,猩猩的发病者数量相对于人类发病者缓慢增加。

SEIR模型综合描述了40周人类与猩猩相互传染的疫情发展状况,前20周“虚拟种群系统”患者数量增加缓慢,后期感染源增多,在缺乏特效药物及隔离强度较低的情况下,患者总量增加较快,符合实际情况,可见模型的合理性。

5.2.3 利用灰色系统进行数据预测(见附件Excel

中的灰色预测数据及结果表

格)

我们先将附件中“虚拟患者种群”的数据用excel进行了处理:分别计算出附件中前40周人类的潜伏群体、处于发病状态、隔离治疗、累计自愈、累计因病死亡的数量,然后在matlab中做出上述五个状态下人类数量关于周数的曲线拟合。(程序参见附件6) 1、潜伏群体 2、处于病状态 3、隔离治疗 4、累计自愈 5、累计因病死亡

利用灰色预测模型得到问题二的预测结果如下:(单位:人)

5.3 问题三模型的建立与求解 5.3.1 模型建立与求解

考虑到人类与猩猩的接触率和药物治愈率两个人为控制因素,在问题二中所建立的SEIR模型的基础上进行了完善。通过降低疑似患者与易感人群的接触,即虚拟种群中患者人数I(t)与易感人群S(t)的接触率r'变小,增加隔离强度P,同时依题将治愈率a提高到80%(aa180%)作图如下:

由图分析知:在未使用药物之前,40周的患者数量已达到2623人,在41周时,专家介入,有效控制易感人群(健康者)与疑似患者的接触率,即对疑似患者进行部分隔离,使得新进入潜伏期的人数在减少。因此,由于时间的延迟,患者人数短期内迅速增长,并达到峰值5436人,其后由于使用某种特效药将隔离治疗人群的治愈率提高了80%,大量的患者被治愈,且受感染的人数越来越少,导致患者人数急剧下降,最后平缓的趋于0,在80周将基本没有患者。和问题二中分别预测的第80、120、200周患者数量36、19、19人对比知,这种特效药物的使用和对患者有效隔离对埃博拉病毒的传播可以进行有效地控制。

所以,在实际对埃博拉病毒的控制中,对于特效药物的研发与使用,以及加强隔离防护措施是非常有必要的。

5.3.2 利用改进后的SIER模型进行数据预测

1、处于发病状态人数

根据上图,可以得到在45、50、55周处于发病状态的人数,即对应周数的坐标I(45)=5418, I(50)=4295, I(55)=2003。

2、退出系统人数(治愈者和死亡者)

在问题二所建立的(4)模型进行改进:退出系统人数

T(t)11a3eptI(t),其中a3=2,P=0.05,利用1中所预测的I(t)数据,即可得到T(45)、T(50)、T(55)。在此基础上,由于人类患者的治愈率a提高到a=%80,得到:累计治愈者人数=aT(t);累计死亡者人数=(1-a)T(t)。

3、潜伏期人数

因在疾病的早起阶段埃博拉病毒可能不具有:“高度的”传染性,所以对问题而中的(2)模型进行改进:取定=0.25,r=43,p=0.05,潜伏期人数Q(t)=reptI(t),利用1 中I(t)的数据,即可得出Q(t)。





4、隔离治疗人数

在此问题中我们取定隔离强度P=0.05,则得到:隔离治疗人数=PQ(t),利用3中Q(t)的预测数据,即可得到45、50、55周的隔离治疗人数。

综上,所得预测结果如下表:

5.4 问题四的模型建立与求解

本问只讨论两个因素对控制埃博拉病毒传播的作用 1.易感人群与疑似患者的隔离控制即改变隔离强度 2.通过防疫.检疫.治疗功能的药物使用改变治愈率 研究结果如下图显示(程序见附件7):

由上图分析可知:四条曲线比较可知,当隔离强度不同时,对患者人数最高峰出现的时间和病毒传播的持续时间(即患者全部痊愈没有再出现患者)有极大的影响。在隔离强度较小时,患者人数的最高峰出现时间靠后,传染病持续的传播时间较长;在隔离强度较大时,患者人数能较快的出现最高峰再较迅速的下降,因此传染病持续的时间比较短,更有利于传染病的控制。所以,在实际的传染病控制过程中,对传染病进行有效的控制,说明实行疑似患者隔离政策对控制传染病传播的效果是很明显的;

分析图可知:当通过使用防疫.检疫.治疗等功能药物使治愈率提高到90%时,峰值明显降低,且提前使患者人数降至“0点”即在短时间内有效的控制了埃博拉病毒的传播 综上,在埃博拉病毒传播过程中,我们可过通过改变隔离强度和治愈率两个指标来对埃博拉疫情进行有效控制。 目前控制埃博拉的防疫措施有:

1、研发强效的防疫药物 2、防止医护工作者感染 3、进行卫生的葬礼 4、隔离感染者 5、追踪接触者

6、设置防疫封锁线,防止疫情向城市蔓延 这些防疫的措施都是通过提高治愈率,及加强隔离制度来有效的控制疫情的发展。只有各种防护措施的结合使用,才能最大限度的控制埃博拉的传播。

六、模型评价与推广

6.1 模型评价:

在建模前期,全面分析数据相互之间的关系,及影响埃博拉病毒传播的因素,找出病毒的作用的时间段和范围,收集比较完整而准确的前期数据。在模型建立中我们采用了各种软件(如matlab,excel等)进行求解,制图精确,计算结果较为准确。但在灰色预测模型中,预测的数据为五周平均数据的近似代替,这对模型的求解结果的准确性有一定的影响。本文所建立的控制模型忽略了人口流动、

变化给该地区病毒传播带来的影响,从而模型预测结果会与实际情况有一定差距。

6.2 模型推广:

通过模型的分析可知,如果全社会的努力和投入的程度继续增加,即隔离措施的提早进行、隔离率增大、防疫药品的早日研发、公众的防御意识提高,可使得疫情周期缩短、患者人数逐步减少。实时监控病毒传播走势,采集更多的数据以验证模型和改进模型,若有预料之外的干扰因素出现,应及时修正模型,重新预测其后期走势。

七、参考文献

[1] 姜启源 谢金星 叶俊,数学模型[M],北京:高等教育出版社,2003(第三版)

[2] 尚涛 谢龙汉 杜如虚,matlab工程计算及分析,北京:清华大学出版社,2011.3

[3] 章绍辉,数学建模,北京:科学出版社,2010

八、附录

【见附件】

范文三:埃博拉病毒传播的数学模型

数 学 建 模 论 文

基于传染病模型的埃博拉疫情发展建模分析

院系名称: 地球科学学院

专业名称: 环境科学

学生姓名:崔昊天、牛星智、朱家佑

学 号: 崔昊天:2013010176

指导教师: 孟得新

完成日期2015 年 4 月 29 日

论文的分工说明:

本篇数学建模论文,是在朱家佑、崔昊天、牛星智的共同努力下完成的。其中朱家佑同学负责SEQIR模型的初步建立、Matlab程序的编译与运行,SIR模型的检验,绘制谱线图,属于首席技术人员;崔昊天同学负责SIR模型的分析与结论,提出SIR模型优化进步的建议以及论文的整体思路与撰写整合,属于文案分析人员;牛星智同学主要负责文献的查考,收集埃博拉病毒的资料以及各国疫情的数据与模型的可行性研究,对比文献资料提出模型存在的不足,并分类整理参考文献,属于认证查考人员。三位同学各司其职,有分工又有合作,环环相扣,共同学习,共同进步。

签名:

日期:2015.4.29

基于传染病模型的埃博拉疫情发展建模分析

摘要

埃博拉(Ebola virus)又译作伊波拉病毒,是一种能引起人类和灵长类动物产生埃博拉出血热的烈性传染病病毒,有很高的死亡率,在50%至90%之间。本文基于西非埃博拉疫情的发展,首先建立SIR模型,对西非埃博拉疫情进行检验,得出模型改进措施,接着建立改进后的SEQIR模型,分析预测疫情发展趋势,对目前世界各国采取的防控埃博拉疫情的措施进行简单评价,并提出埃博拉疫情对于防控的建议。

关键词:埃博拉;检验;预测;评价;建议

目录

第1章 绪论 ............................................................................................................................... 5

1.1研究背景 ............................................................................................................................. 5

1.2 研究目的及意义 ................................................................................................................ 5

1.3 国内外研究现状 ................................................................................................................ 5

1.3.1 国内外传染病模型研究现状 ......................................................................................... 5

1.3.2 国内外埃博拉病毒医疗研究现状 ................................................................................. 6

1.3.2.1 EVD治疗方法 ............................................................................................................. 6

1.3.2.2治疗 EVD 的新药研究进展 ....................................................................................... 7

1.3.2.3疫苗的研究 ................................................................................................................... 9

1.4 技术路线 .......................................................................................................................... 10

第2章 数据与方法 ................................................................................................................. 11

2.1基础知识 ........................................................................................................................... 11

2.2 数据方法 .......................................................................................................................... 11

第3章 模型建立与改进 ......................................................................................................... 12

3.1 SIR模型 ........................................................................................................................ 12

3.1.1 模型建立 ....................................................................................................................... 12

3.1.2 模型分析 ....................................................................................................................... 13

3.1.3 模型验证 ....................................................................................................................... 16

3.2 改进的SIEQR模型 ................................................................................................... 16

3.2.1 模型建立 ....................................................................................................................... 16

3.2.2 模型求解 ....................................................................................................................... 18

3.2.3 模型分析 ....................................................................................................................... 19

第4章 评价与建议 ................................................................................................................. 24

4.1 模型的评价 ...................................................................................................................... 24

4.1.1 SIR模型的评价 ............................................................................................................ 24

4.1.2 SEQIR模型的评价 ....................................................................................................... 24

4.2 现有医疗水平的评价 ...................................................................................................... 24

4.3防控建议 ........................................................................................................................... 25

致谢 ................................................................................................................................................ 26

参考文献......................................................................................................................................... 27

第1章 绪论

1.1研究背景

埃博拉(Ebola virus)又译作伊波拉病毒,是一种能引起人类和灵长类动物产生埃博拉出血热的烈性传染病病毒,有很高的死亡率,在50%至90%之间。

2014年2月西非爆发大规模病毒疫情,截至2014年12月17日,世界卫生组织(WHO)发表数据显示埃博拉出血热疫情肆虐的利比里亚、塞拉利昂和几内亚等西非三国的感染病例(包括疑似病例)已达19031人,其中死亡人数达到7373人[1]。据世界卫生组织透露,2014年西非埃博拉病毒疫情爆发的感染及死亡人数均为历史最高,并仍处于恶化状态中[17]。世界各国针对埃博拉病毒均提出了自己的防控措施,然而埃博拉病毒的致病机理尚未明确,其致死率高居不下,基于数学模型的病毒爆发特性分析因此应运而生。

1.2 研究目的及意义

国内外学者针对一般类型的传染病提出了SIR传染病模型,可用于模拟预测传染病的传播过程,评价当前的医疗卫生体系及防控措施。但是针对埃博拉这一目前无疫苗防治措施的传染病,尚未存在基于SIR模型的改进模型进行合理的模拟预测。

本文基于SIR模型,增加疑似病患、潜伏期感染者两个变量,构建SIEQR模型,分析埃博拉出血热的流行规律,并预测其未来的发展趋势,并对世界卫生组织(WHO)目前针对该疾病的防控措施进行评价,为埃博拉出血热的最佳防控提供理论基础。

1.3 国内外研究现状

1.3.1 国内外传染病模型研究现状

国内外对于传染病的数学模型研究不在少数。传染病模型的基础是Kermack和McKendrick 二人提出的仓室模型(compartmentmodel,CM)[2]。杨方廷[18]通过对SARS在某北方大城市爆发的系统仿真研究发现,恰当地建模能够有效地模拟SARS的人群传播模式。通过数学模型可推导出不同的防控策略(如限制人员流动、对疑似感染者及早隔离等)对于病毒传播的影响。

对于埃博拉病毒的数学模型研究,早在 1996年,Fauci[4]就使用SIR和 SEIR 模型,模拟扎伊尔两个时段的埃博拉爆发:1976 年 Yambuku 疫情爆发和 1995年Kikwit的疫情爆发。他们得到: 当基本再生率K满足范围1.72≤R0≤8.60时,

意味着埃博拉病毒传染性不如以前那么厉害,可以使他们减少潜在的死亡。为了预判埃博拉的发展趋势,弗吉尼亚理工大学的网络动态和模拟科学实验室(NUSSL)从2014年8月开始,采用基于个体的传染病传播工具Epi- fast来模拟埃博拉在西非的传播[22]。

国外学者在研究其发病机理以及病毒作用机理的同时,建立了新型传染病模型对疾病进行预测。由于个体在感染埃博拉病毒后经过一段时间才会出现临床症状,这段时间被称为潜伏期,对潜伏期的长度,Chowell等[9]采用最小二乘法估算出2次疫情平均潜伏期分别为5.30 d(刚果1995年)和3.35 d(乌干达2000年),Le- grand等使用临床观测结果估计的潜伏期分别为7d[12-13](刚果1995年)和12 d (乌干达,2000),与后者较为近似的潜伏期长度也被用于近期对此次西非埃博拉疫情动态的研究中[15]。美国CDC采用的潜伏期平均长度为6 d,99%的分位数为21 d。

1.3.2 国内外埃博拉病毒医疗研究现状

1.3.2.1 EVD治疗方法

由于尚无有效的经过系统研究的药物用于 EVD 患者的治疗,目前的治疗方法以积极的对症和支持治疗为主,主要包括 :维持水、电解质平衡,补充体液和电解质 ;预防和控制出血 ;维持血氧及血压平衡 ;及时控制继发感染 ;治疗肾功能衰竭和出血、DIC(disseminated intravascular coagulation,弥散性血管内凝血 ) 等并发症。临床治疗过程中可根据患者的具体情况采取下列措施进行治疗。

①一般支持、对症治疗 :卧床休息,少渣易消化半流质饮食,保证充分热量 ;

② 补液治疗 :有证据表明,早期补液,维持水电解质和酸碱平衡治疗,可明显提高存活率 ;可使用平衡盐液,维持有效血容量 ;加强胶体液补充如白蛋白、低分子右旋糖酐等,预防和治疗低血压休克 ;

③保肝抗炎治疗 :可使用甘草酸制剂 ;

④出血的治疗 :止血和输血,新鲜冰冻血浆补充凝血因子,预防 DIC ; ⑤ 控制感染 :及时发现继发感染,根据细菌培养和药敏结果应用抗菌药物 ; ⑥ 肾功能衰竭的治疗 :必要时行血液净化治疗 ;

⑦呼吸衰竭的治疗 :及时行氧疗等改善呼吸功能治疗 ;

⑧ 病原学治疗 :未经过人体学试验的三联单克隆抗体 ZMapp,在紧急状态下被批准用于 EVD 患者的治疗 ;目前已有 7 人接受此治疗,5 人获得较好疗效;另外,TKM-Ebola 也已通过 FDA 的批准,紧急情况下可用于确诊或者疑似 EVD 患者 ;

⑨ 恢复期血清治疗 :即从康复者体内提取含埃博拉病毒抗体的血液,将其

输给患者 ;该方法曾在小范围内应用,亦似有较好的效果,WHO 已出台安全使用输血疗法治疗 EVD患者的指导方针,并帮助受疫情影响国家从康复患者体内提取血浆、准备血清以治疗更多患者。但和 ZMapp 一样,目前尚无数据支持其大规模使用的效果,还有待于在应用时机、不良反应等方面做进一步观察。

1.3.2.2治疗 EVD 的新药研究进展

在埃博拉病毒被发现后的几十年里,研究者们一直致力于研发抗埃博拉病毒的新药。不过目前,用于埃博拉病毒预防和治疗的药物仍多处于实验室或前期临床研究阶段。在目前的严峻形势下,美国 FDA 已经批准ZMapp 和 TKM-Ebola 在紧急情况下用于确诊或者疑似感染埃博拉病毒的患者。WHO 委员会也给出建议,可在患者签署知情同意书的基础上,使用未经系统研究证明的药物如 ZMapp 和疫苗来控制感染,而且在伦理上也可接受。不过,即便如此,现有的可用药品也是非常有限或缺乏的。

1.ZMapp

ZMapp 是美国和加拿大共同开发的一种抗埃博拉病毒新药,由 3 种人源化的鼠单克隆抗体组成,其中的主要成分为 MB-003(又称为 :Mapp)和 ZMAb,两者优化后组成了抗埃博拉病毒感染的新药 ZMapp。MB-003是 3 种具有抗埃博拉病毒作用的单克隆抗体组成的混合抗体。2012 年,研究者用 MB-003 治疗了 6 只埃博拉病毒感染的恒河猴,给予 MB-003 的实验组中有 4 只存活;而未给药的对照组 2 只全部死亡。2013 年一项研究发现,感染埃博拉病毒 48 h 内注射 MB-003 混合抗体的恒河猴,其存活率为 43% ;如果不接受药物治疗,其存活率为零。ZMAb 是抗埃博拉病毒混合抗体联合腺病毒载体干扰素 α 的复合抗体。2013 年的另一项研究显示,使用 ZMAb 治疗感染埃博拉病毒 24 h 或 48 h 的猴,存活率分别为 100% 和 50%。2 个多月后,研究人员对存活的猴子进行再次感染试验,观察这些病愈的猴子体内是否产生了抵御埃博拉病毒的抗体,结果无 1 例猴子因再次感染而死亡。为进一步了解ZMapp的抗病毒作用,Qiu 等将致死剂量的埃博拉病毒给予三组猴子(每组6 只),分别在不同时间给予 ZMapp 治疗,结果发现所有动物都存活,更重要的是,即使是在猴感染病毒后 5 d 出现严重症状时使用,该药也能起到很好的保护作用。不过要注意的是,试验时所用毒株与本次西非暴发的毒株不同。不管怎样,ZMapp 的开发与猴试验中的良好疗效都是对抗埃博拉病毒的巨大进步。

2.TKM-Ebola

鉴于 ZMapp 的产能有限和疫情的不断扩展,国际社会也在尽可能的寻找更多的治疗 EVD 的有效药物。为此,2014 年 9 月美国 FDA 批准加拿大制药商 Tekmira公司研制的 TKM-Ebola 药物可在紧急情况下用于确诊或者疑似感染埃博拉病毒的患者。TKM-Ebola 也成为全球第二个可用于 EVD 治疗的药物。

TKM-Ebola 属于基因治疗药物,是一种靶向埃博拉病毒 RNA 多聚酶的小干扰 RNA,通过抑制埃博拉病毒相关基因 VP24 和 VP35,使致病基因 “ 沉默 ” 而起作用。2010 年发表在《柳叶刀》上的小规模的概念验证研究中,研究者将猕猴暴露于高剂量的扎伊尔型病毒后 30 min,开始给予 TKM-Ebola,随后,第一组猕猴(n=3)在暴露后 1、3、5 d 各给药 1次,而第二组猕猴(n=4)在暴露后 1 ~ 6 d 每天给药 1次,结果第一组的保护率可达 66%,第二组的保护率可达 100%。但 FDA 对其在人类的安全性上持保留态度,故仅部分准予其用于埃博拉病毒感染的患者。

3.jk-05

在应对 EVD 方面,中国的科研工作者也未曾停止前进的步伐。其中由中国军事医学科学院微生物流行病研究所针对埃博拉病毒、历时 5 年研制而成的药物 jk-05已通过总后卫生部专家评审,获得军队特需药品批件。jk-05 是一种小分子化学药物,该药能够选择性地抑制流感病毒和埃博拉病毒的 RNA 聚合酶,从而达到抑制病毒复制的目的。研究表明,jk-05 在细胞和动物水平感染实验中具有抗埃博拉病毒活性。该药物虽已完成临床前研究,并完成了用于抗流感病毒时的临床安全性评价,但对于 EVD 的治疗,目前仍仅限于紧急情况下使用。

4.法匹拉韦(favipiravir)

日本相关法规规定,在特定的条件下可采用未获批准治疗埃博拉的药物来治疗致死性的埃博拉病毒感染。根据此规定,日本富士公司和美国 MediVector 公司正在向美国 FDA 提交申请以扩大法匹拉韦(favipiravir,或称 T-705,Avigan)这一原本应用于抗流感病毒的药物,作为抗埃博拉病毒感染用药。法匹拉韦是一种新型的抗病毒化合物,主要用作抗新型流感毒株感染,今年年初,被日本政府批准上市。由于其没有经过猿猴以及人类相关研究,法匹拉韦能否抗埃博拉感染尚不明确。不过埃博拉病毒与流感病毒具有类似的结构(均为单链反义RNA 病毒),因此法匹拉韦可能对埃博拉病毒感染同样有效。但 WHO 目前并未推荐该药用于 EVD 的治疗。

5.其他用于 EVD 治疗的药物

除上述 4 种药物外,还有一些对抗埃博拉病毒的药物正在开发中,虽然很多药物均表现出良好的前景,但同样存在着相似甚至更大的不确定性。BCX-4430 是一个新的核苷类似物,可以影响病毒 RNA 多聚酶的功能而抑制病毒的复制,在啮齿类动物中具有抗埃博拉病毒活性,在非人灵长类动物中能够保护其免受马尔堡病毒感染,但并没有安全性试验的数据。AVI-7537 可通过RNA 干扰技术靶向作用于埃博拉病毒 VP24 蛋白,实验中对非人灵长类动物具有一定的保护作用。还有一些明确安全用法但用于其它疾病的药物,如氯喹、伊马替尼、雌激素受体调节剂(克罗米芬和托瑞米芬)、钙通道阻滞剂(胺碘酮、维拉帕米)等,也在

体外实验中体现出抗埃博拉病毒活性,可作为候选药物单独或联合其他药物用于治疗 EVD。

1.3.2.3疫苗的研究

寻找或研发抗埃博拉病毒的治疗药物固然重要,但管理和控制埃博拉病毒疫情暴发的最有效手段应该还是疫苗的使用,WHO 正与各方加快埃博拉病毒疫苗的临床试验。目前有数种针对埃博拉病毒的疫苗表现出良好的开发前景,但同样需要进一步的研究才能明确其应用价值。其中葛兰素史克公司开发的 “ChAd” 和美国纽琳基因公司开发的 “VSV” 疫苗最受关注,WHO 正与药物研发公司、临床专家与监管方加快这两种疫苗的临床试验,希望受疫情影响的国家能在 2015 年 1 月以后用上这些疫苗。

ChAd 疫苗是针对扎伊尔株设计研发的疫苗,它以一种良性病毒(腺病毒)为载体,仅包含了埃博拉病毒基因组的部分结构,因此不会导致接种疫苗者感染埃博拉病毒,该腺病毒也不会在人体内复制,但是如果人体处于感染状态时却会诱导免疫系统产生埃博拉病毒抗体。之前的动物实验已经证实该疫苗能够有效地防止埃博拉病毒感染,并且单剂量疫苗即可诱导产生高水平的保护效果,促进感染动物的恢复。目前,ChAd 已进入Ⅰ期临床试验阶段,这是首个进入Ⅰ期临床试验的抗击埃博拉病毒的疫苗。该试验已于 2014 年 9 月在牛津大学疫苗研究中心进行,有 60 名志愿者参与,目的是研究疫苗对健康人体是否有害,同时进一步确定最合适的疫苗剂量,但最终确定 ChAd 的疗效还要进行一系列试验。然而,在当下的情况下,人们希望把通过Ⅰ期临床试验的疫苗直接送到疫区的医护工作者手中,之后直接在疫区检测疫苗的效果。

VSV 疫苗是活病毒,该疫苗使用水泡性口炎病毒进行接种,可自我复制,感染细胞后可将埃博拉病毒的蛋白质带入宿主细胞中,之后就会激发免疫反应,防止机体出现器官衰竭和出血等症状。VSV 疫苗的价值在于,对于已经感染的猴子来说,该疫苗依然有效。美国纽琳基因公司已批准开展 VSV 疫苗的Ⅰ期临床实验。

1.4 技术路线

第2章 数据与方法

2.1基础知识

微分方程可以很好的模拟传染病的流行方式以及人感染传染病后的状态,SIR传染病模型是我们建立模型的理论基础,通过对SIR模型的改进,完善埃博拉病毒的传染机理,经过世界卫生组织数据带入验证,确立模型的准确性、科学性。

2.2 数据方法

根据已知埃博拉病毒的疾病生态学,同时考虑到埃博拉较高的病死率、短潜伏期、目前尚无疫苗防疫的特点,对传播动力学过程有较大的影响,本模型在SIR三分室的基础上增加疑似病患E,潜伏期感染者Q。设定的五分室模型如下:

各仓室之间关系如下:

dSQ(t)(1p)S(t)dtdR1I(t)dtd3dE1Q(t)(1p)[E(t)(1p)E(t)p]dtd3 dQ21I(t)Q(t)dtd1d2d3dI12Q(t)(1p)[S(t)E(t)(1p)E(t)p]I(t)d3d1d2dt

这样,就可以更加详细地描述模型变量,更好的预测埃博拉传染趋势。 本文首先运用SIR模型验证传染病传播的规律,再对模型进行参数修正,得出SEQIR模型预测某地区埃博拉疫情的发展。

第3章 模型建立与改进

3.1 SIR模型

3.1.1 模型建立

针对一般的传染病,姜启源等《数学模型》中建立了三个数学模型对其进行模拟,并逐步改进出了SIR模型,现将其作为初步分析的工具,具体过程如下。

模型假设:

1.在疾病传播期内所考察的地区范围不考虑人口的出生、死亡、流动等种群动力因素。总人口数N(t)不变,人口始终保持一个常数N。人群分为以下三类:易感染者(Susceptibles),其数量比例记为s(t),表示t时刻未染病但有可能被该类疾病传染的人数占总人数的比例;感染病者(Infectives),其数量比例记为i(t),表示t时刻已被感染成为病人而且具有传染力的人数占总人数的比例;恢复者(Recovered),其数量比例记为r(t),表示t时刻已从染病者中移出的人数(这部分人既非已感染者,也非感染病者,不具有传染性,也不会再次被感染,他们已退出该传染系统。)占总人数的比例。

2.病人的日接触率(每个病人每天有效接触的平均人数)为常数r,日治愈率(每天被治愈或死亡的病人占总病人数的比例)为常数μ,显然平均传染期为1/μ,传染期接触数为σ=r/μ。

基于以上模型假设,进行建模过程,我们看出,在假设1中显然有:

s(t) + i(t) + r(t) = 1 (3.1.1.1)

对于病愈免疫或死去的移出者的数量应为

NdrNidt (3.1.1.2)

s0>0)0我们不妨设初始时刻的易感染者,染病者,恢复者的比例分别为s(,

i0(i0>0),r0=0.

SIR基础模型用微分方程组构成表示如下:

didtsii

dssi

dt

drdti

 (3.1.1.3)

方程组中各参数及其含义如下:

日治愈率(每天被治愈或死亡的病人占总病人数的比例) μ

病人的日接触率(每个病人每天有效接触的平均人数)为常数 r

传染期有效接触数 σ

时刻t内健康人群 S(t)

时刻t内潜伏期人群 I(t)

时刻t内治愈或死亡人群 R(t)

则,在方程(3)中设r=1,μ=0.1,i(0)= 0.00001,s(0)=0.99999,用MATLAB软件编程:

运行后得出图像如图3-1所示。

图3-1 感染者比例与时间关系

3.1.2 模型分析

我们在数值计算和图形观察的基础上,利用相轨线讨论解i(t),s(t)的性质。

i ~ s平面称为相平面,相轨线在相平面上的定义域(s,i)∈D为

D = {(s,i)| s≥0,i≥0 , s + i ≤1} (3.1.2.1) 在方程(3)中消去dt并注意到σ的定义,可得

di11dssσ, i|ss0i0 (3.1.2.2)

所以,有

1di1dssσs1d1dsii0s0sσ (3.1.2.3) i

利用积分特性容易求出方程(5)的解为

i(s0i0)s1

lns

s0 (3.1.2.4)

在定义域D内,(6)式表示的曲线即为相轨线,如图3-2所示。其中箭头表示了随着时间t的增加s(t)和i(t)的变化趋向。

i

D

P2

imP1

os1/s

图3-2 s(t)和i(t)的变化曲线

下面根据(3.1.1.3)、(3.1.2.4)式和图3-2分析s(t)、i(t)和r(t)的变化情况(t→∞时它们的极限值分别记作s,i和r)。

1. 不论初始条件s0,i0如何,病人消失将消失,即

i00 (3.1.2.5)

其证明如下:

dsdr00s首先,由(3.1.1.3) dt 而 s(t)0 故 存在; 由(3.1.1.2)dt而

r(t)1 故r 存在;再由(3.1.1.1)知i存在。

dri0则由(3.1.1.1),对于充分大的t 有dt2 , 这将导致其次,若

r,与r存在相矛盾.从图形上看,不论相轨线从P1或从P2点出发,它终将

与s轴相交(t充分大)。

2.最终未被感染的健康者的比例是s,在(7)式中令i=0得到,s是方程

s0i0s1lns0s0 (3.1.2.6)

在(0,1/σ)内的根。在图形上s是相轨线与s轴在(0,1/σ)内交点的横坐标。

di1di11o1 3.若s0>1/σ,则开始有dssσ,i(t)先增加,令dssσ=0,可

得当s=1/σ时,i(t)达到最大值:

1ims0i01lns0) (3.1.2.7)

di11o

然后s

di110ss4.若01/σ,则恒有dssσ,i(t)单调减小至零,s(t)单调减小至,

如图3中由P2(s0,i0)出发的轨线。

可以看出,如果仅当病人比例i(t)有一段增长的时期才认为传染病在蔓延,那么1/σ是一个阈值,当s0>1/σ(即σ>1/s0)时传染病就会蔓延。而减小传染期接触数σ,即提高阈值1/σ使得s0≤1/σ(即σ ≤1/s0),传染病就不会蔓延(健康者比例的初始值s0是一定的,通常可认为s0接近1)。

并且,即使s0>1/σ,从(19),(20)式可以看出, σ减小时, s增加(通过作图分

i析),m降低,也控制了蔓延的程度.我们注意到在σ=λμ中,人们的卫生水平越高,日接触率λ越小;医疗水平越高,日治愈率μ越大,于是σ越小,所以提高卫生水平和医疗水平有助于控制传染病的蔓延。

从另一方面看,ss1/是传染期内一个病人传染的健康者的平均数,称为交换数,其含义是一病人被s个健康者交换。所以当 s01/ 即s01时必有σS≤1既然交换数不超过1,病人比例i(t)绝不会增加,传染病不会蔓延。

3.1.3 模型验证

几内亚是联合国公布的最不发达国家之一,人口1120万。2014年3月23日WHO宣布几内亚埃博拉病毒感染所致的埃博拉病毒病,又称埃博拉出血热暴发,并在几内亚首都科纳克里和利比里亚和塞拉利昂的边界也出现EHF疑似病例。

下表显示了2014年感染埃博拉病毒人员在几内亚的变化情况:

表3-1 2014年感染埃博拉病毒人员在几内亚的变化情况表

由数据可知在前十天内,病毒传播情况大致符合模型模拟情况,可是其增长速度并没有模拟的那么快。究其原因,是因为病毒爆发后,政府采取了相应的隔离措施,而减小传染期接触数σ,即提高阈值1/σ使得≤1/σ(即σ ≤1/),传染病就不会蔓延。下面通过改进的SIEQR模型进行模拟与预测某未知地区的埃博拉病毒爆发情况。

3.2 改进的SIEQR模型

3.2.1 模型建立

在埃博拉病毒爆发的初期,由于潜伏期的存在,社会对埃博拉病毒传播的速度和危害程度认识不够,政府和公众都没有引起足够的重视。当人们发现被感染者不断增加时,政府开始采取多种措施以控制埃博拉病毒的进一步蔓延。通过上述分析,可以把埃博拉病毒的传播过程大致分为3个阶段:

① 控制前,接近于自然传播时的传播模式;

② 过渡期,从公众开始意识到埃博拉的严重性到政府采取严格控制之的一段时间;

③ 控制后,在介入人为控制因素,严格控制之后埃博拉病毒的传播模式。 除西非等不发达地区外,美国、西班牙等国家都是在埃博拉病毒刚刚大肆传播就采取了强有力的措施,因此,这些地区的过渡期都可包括在控后期。由以上分析可知,单纯的SIR模型对于采取控制前埃博拉疫情的预测基本与实际情况较为吻合,但不能反映相关人为因素介入后的病毒传播规律和疫情的实际情况。在考虑相关人为控制因素后,针对SIR模型的不足,我们对加入了隔离等强制措施的SIR模型进行改进。

首先是人群种类的变化,我们在结合理论分析的基础上,把考察范围内的人群分为以下种类:

1、健康人群,即易感染(Susceptibles)人群。记其数量为S(t),表示t时刻未感染病但有可能感染该疾病的人数;

2、确诊患者,即被感染(Infection)该疾病的人群,记其数量为I(t) 表

示t时刻可能感染该疾病的但又不是疑似病患的人数;

3、疑似病患,记其数量为E(t)表示t时刻感染该疾病的并是疑似病患的人数;

4、潜伏期感染者,记其数量为Q(t)表示t时刻已经感染病毒但没有表现症状即处在潜伏期的人数;

5、恢复人群(Recovered),记其数量为R(t),表示t时刻已从感染病者中移出的人数,恢复者或死亡者(这部分人数既不是已感染者,也不是非感染者,不具有传染性,也不会再次被感染,他们已经推出了传染系统)。

基于以上人群分类,我们假设(埃博拉病毒传入某城市):

1. 假设易感人数的变化率与当时的易感人数和感染人数的乘积成正比;

2. 假设从感染数中移除个体的速率与当时的感染人数成正比;

3. 假设考察地区内疾病传播期间忽略人口的出生,死亡,流动等种群动力因素对总人数的影响。即:总人口数不变,记为N。

4. 假设潜伏期人群不会传染健康人,不具有传染性。

5. 假设被隔离的患者无法跟别人接触,不会传染健康人。

6. 假设治愈者已对该病毒有免疫力,不会再被该传染病传染,可以退出系统

7.假设初始时刻健康人群的总人数为S0=1.1千万(某城市),潜伏期的总人数为I0=1,疑似病患的总人数为E0=0,确诊病患的总人数为Q0=0,恢复人群的总人数为R0=0。

同样,方程组中参数及其含义如下:

病毒潜伏期(天) d1、d2

病患者治愈(死亡)时间(天) d3

病患人均每天接触人数 r

防控措施强度 p

时刻t内健康人群 S(t)

时刻t内潜伏期人群 Q(t)

时刻t内病症疑似人群 E(t)

时刻t内已患病人群 I(t)

时刻t内治愈或死亡人群 R(t)

传染病传染率 

由模型的假设得到如下关系

S(t)+I(t)+E(t)+Q(t)+R(t)=N(3.2.1.1)

根据假设在时刻t内健康人群变化有

S(tt)S(t)Q(t)(1p)S(t)t (3.2.1.2)

在时刻t内治愈或死亡人群的变化有 R(tt)R(t)11I(t)td3(d3为单位时间内患者的恢复率)(3.2.1.3)

在时刻t内病症疑似人群的变化有

E(tt)E(t)Q(t)(1p)[E(t)(1p)E(t)p1]td3(3.2.1.4)

在时刻t内已潜伏期人群的变化有(已患病人群等于潜伏期病人转为感染者减去移除人数)

Q(tt)Q(t)21I(t)Q(t)d1d2d3(3.2.1.5)

在时刻t内患病人群的变化有

12I(tt)-I(t){Q(t)(1-p)[S(t)E(t)(1-p)]-I(t)}td3d1d2 (3.2.1.6) 2

(d1d2为单位时间内潜伏期病人转为感染者的比例常数)

根据以上变化有方程组:

dSQ(t)(1p)S(t)dtdR1I(t)dtd3dE1Q(t)(1p)[E(t)(1p)E(t)p](3.2.1.7) dtd3dQ21I(t)Q(t)dtd1d2d3dI12Q(t)(1p)[S(t)E(t)(1p)E(t)p]I(t)dtddd312

3.2.2 模型求解

当d1=6, d2=12, d3=10, r=20, p=80%,我们假设患者2天后入院治疗,疑似患者2天后被隔离。有初始状态的患者人数为:

QQ(0)*(NQ(0)E(0)*r)*2N (3.2.2.1)

运用Matlab对模型进行模拟运算,得出患者人数随时间变化如图一:

图3-3 感染者比例与时间关系(P=0.8)

由上图可以得到:在当当d1=6, d2=12, d3=10, r=20, p=80%,在患者2天后入院治疗以及疑似患者2天后被隔离的条件下,当t在0~9.6607d时,患者的人数是急剧上升的,在t=9.6607达到最大值,此时患者人数为y=140822.1475在采取医疗措施,比如患者入院治疗,隔离疑似患者等后患者人数随着时间的增长呈现下降的趋势。

3.2.3 模型分析

通过上述建模过程,我们可继续分析该模型特点,其中b1,b2,b3为固定值,当P值(防控措施强度)改变时,我们得到以下几幅不同的图像。

图3-4 感染者比例与时间关系(P=0)

图3-5 感染者比例与时间关系(P=0.2)

图3-6 感染者比例与时间关系(P=0.4)

图3-7 感染者比例与时间关系(P=0.6)

图3-8 感染者比例与时间关系(P=0.8)

图3-9 感染者比例与时间关系(P=1.0)

通过Matlab拟合的曲线,我们可以很明显的看出当防控强度P为0.4时,可以达到最佳效果(患者最少,患病高峰期延长),而我们主观意识上将患者与疑似患者全部隔离是最理想状态,但曲线拟合结果显示,这反而会使患者数在初期剧增(t≈0)

,这显然是不符合常规的,况且在现实情况下不可能做到防控强

度为1,因此模型预测取P

第4章 评价与建议

4.1 模型的评价

4.1.1 SIR模型的评价

SIR模型中考虑了传染病的几个基本因素,对埃博拉疫情进行了合理的模拟,并对模型进行验证,得出模拟效果良好但仍有瑕疵,于是说该模型既是基础,又有待改进。

4.1.2 SEQIR模型的评价

模型优势:

该模型较好的弥补了SIR模型中欠缺考虑的埃博拉病毒潜伏期短、致死率高等特点,考虑到世界卫生组织针对埃博拉疫情的控制措施,增添多项参数,为埃博拉疫情在某个未知地区的出现提供了理论分析基础,可以预测病毒爆发的全过程并为医疗组织提供合理的防控力度。

模型缺陷:

1、尚无合乎模型假设的区域(美国、西班牙等地虽符合假设,但是病例较少、医疗水平对较高,无爆发迹象,因此无法纳入考察范围)可供模型验证;

2、仍未把握到埃博拉病毒致死率较高的特点,仍有改进空间;

3、病患在达到峰值时的下降速度仍较快,防控措施与模型的相关性仍相对较弱,与实际可能会产生误差。

4.2 现有医疗水平的评价

针对现有医疗水平的评价,本文还是运用SIR模型的数据进行分析,因为发达国家病例少,样本不具有代表性,而西非地区是埃博拉疫情的重灾区,考虑到其医疗卫生条件相对发达国家较差,因此较为符合不考虑隔离强度、无疑似患者仓室的系统,即SIR系统。

通过对SIR模型验证时我们发现对于西非地区采取的防控埃博拉疫情的措施效用仍未达到模型理想的水平,说明西非目前不能保证日接触率λ处于相对较低、日治愈率μ相对较高的水平,即尚不能保证日接触数σ处于较低水平,其防控体系与防控力度相对较弱。但是,西非、世界卫生组织、各国卫生部采取了一系列强制措施,如:加强出入境检疫,积极研制埃博拉病毒疫苗,隔离埃博拉患者等,这些措施起到了防控埃博拉在世界范围内传播的效果,该控制体系的优越性是我们所不能忽视的。

4.3防控建议

目前尚没有针对 EVD 的特异性疫苗,也没有特异性治疗办法,因而严格隔离控制传染源、密切接触者追踪、管理和加强个人防护是防控 EVD 的关键措施。

预防性措施 对于未发生埃博拉出血热流行的国家和地区而言,有效控制传染源的关键是及时发现和隔离控制输入性病例,为此要加强进出境检疫,加强对动物的检疫,尤其是要加强对非人灵长类和蝙蝠等野生动物的检疫工作。口岸检疫部门一旦发现可疑病例或动物,要及时通报卫生部门做好疫情调查和处理,对前往该病流行地区的旅游和进行医护工作的人员,要进行有关埃博拉出血热防病知识的普及和教育,使其提高警惕意识,做好个人防护,加强国际间信息交流与合作,密切关注该病的流行动态,尤其要高度关注曾出现过埃博拉出血热流行的地区。

疫情控制措施 一旦发现可疑病例及其接触者,应立即采取严格的隔离措施,以防止疫情的扩散及流行。与病人接触时做好个人防护,防止直接接触患者污染物。对病人的分泌物、排泄物要严格消毒,污染的针头、注射器等可用焚烧或高压蒸汽进行消毒处理,所有涉及埃博拉病毒的实验室活动应严格按照我国有关规定进行,相关的实验室检查应减少至需要的最低限度,采集病料时应注意个人防护,在病料送检过程中一定要注意防止病毒的散播,病毒分离与培养只能在生物安全 4级实验室进行,该病无特效治疗措施,主要采用对症和支持疗法,注意水,电解质平衡,预防和控制出血,控制继发感染。

致谢

在此,首先感谢孟得新老师教受我们数学建模的基本理论,让我们感受到了数学建模的独有魅力,开发了我们的思维,提升了我们的能力,拓宽了我们的视野。孟老师严谨的治学态度、一丝不苟的钻研精神使我们学习的榜样,在论文遇到瓶颈期的时候,孟老师的指导让我们找到了新的思路,有了新的发现,才得以完成最终的文章。我还要深刻感谢朱家佑、牛星智同学的大力支持与帮助。如果没有朱家佑同学强大的理科思维与熟练的软件操作能力,我们的模型不可能建立的如此迅速,如此富有创新性;如果没有牛星智同学彻夜查考文献、整理文献综述以及模型的改进意见,我们也不可能顺利完成工作,我也不可能有清晰的思路进行论文的撰写与整合。最后,还要感谢帮助和支持我的室友、同学们,谢谢你们的提点与支持。总之,是大家的努力造就了这篇论文,谢谢你们。

参考文献

[1] 曾德荣,成晓红,黄吉安.对SARS疫情预测与控制的数学模型研究——一种基于一般流传

病模型的改进方法[J].中国安全科学学报,2004,02:44-48.

[2] 徐致靖,徐展凯,刘列,许晴,祖正虎,郑涛. 对想定北京和广州发生输入性埃博拉疫情的计

算实验研究[J].生物技术通讯,2015,01:33-39+50.

[3] 吴开琛,吴开录,陈文江,林明和,李才旭. SARS传播数学模型与流行趋势预测研究[J].中国

热带医学,2003,04:421-426.

[4] Fauci AS:Ebola-underscoring the global disparities in health care resources[J].The New

England Journal of Medicine,2004,371:1084-1086.

[5] 张齐鹏.甲型H1N1流感在预防控制措施下的传播数学模型构建[J].湖南文理学院学报(自

然科学版),2009,04:16-18.

[6] 周后卿,徐幼专.埃博拉病毒感染数量的一个数学模型[J].邵阳学院学报(自然科学

版),2014,04:1-5.

[7] 徐展凯,徐致靖,祖正虎,郑涛.基于个体的传染病传播模型应用[J].生物技术通

讯,2015,01:46-50.

[8] 徐致靖,徐展凯,刘列,祖正虎,郑涛.关于西非埃博拉疫情态势预测研究的总结分析[J].生物

技术通讯,2015,01:22-32.

[9] 刘畅,丁光宏,龚剑秋,王凌程,珂张迪.SARS爆发预测和预警的数学模型研究[J].科学通

报,2004,21:2245-2251.

[10] 孙海,熊思灿,吴志强.基于改进SIR模型的甲型H1N1流感防控研究[J].东华理工大学学

报(自然科学版),2014,01:96-100.

[11] 张云辉,王姝,陈玉琪,李军.埃博拉出血热研究现状及2014年疫情进展[J].传染病信

息,2014,04:3-8.

[12] 程颖,刘军,李昱,刘翟,任翔,施一,高福,余宏杰.埃博拉病毒病:病原学、致病机制、治疗与

疫苗研究进展[J]. 科学通报,2014,30:2889-2899.

[13] 刁天喜,徐守军,赵晓宇,周巍,李鹏,刘术,王伟.埃博拉出血热药物和疫苗研究开发态势分

析[J]. 军事医学,2014,08:569-575.

[14] 孟现民,董平,卢洪洲.埃博拉病毒病的治疗及新药研究进展[J]. 上海医药,2014,21:1-5. [15] 凌春英,马秀慧.甲型H1N1流感的微分方程模型[J]. 高师理科学刊,2012,01:17-19. [16] 高秋月,肖露平,李海燕,方献平. 埃博拉病毒及其免疫研究进展[J]. 生物学教

学,2009,07:7-9.

[17] World Heal ill Organization. Ebola response roadmap situation report[R]. 2014-12-24. [18] 杨方廷.北京SARS疫情过程的仿真分析[J].系统仿真学报, 2003, 15(7): 991-998. [19] 姜启源,谢金星,叶俊.数学模型(第四版)[M].北京:高等教育出版社,2010.

[20] 宋建德,袁丽萍,朱迪国,魏荣. 近期全球埃博拉病毒病流行情况及防控进展[J]. 中国动物

检疫,2014,09:1-6.

[21] 夏生林. 埃博拉出血热的防控[J]. 畜牧与饲料科学,2010,Z2:162-163.

[22] .

[23] 唐关丽. SARS数学模型的建立及其解的存在唯一性[J].现代经济信息,2009,01:98+97.

范文四:埃博拉病毒传播问题的数学模型

摘 要

本文分别研究在单个种群、多个种群以及有防控措施、无防控措施条件下,埃博拉病毒在动物群之间的传播特性,建立模型,量化埃博拉病毒的传播规律,深刻认识该病毒的危害,并分析隔离措施的严格执行和药物治疗效果的提高等措施对控制疫情的作用。

问题一

针对描述病毒在“虚拟猩猩种群”中的传播特性并预测接下来的在猩猩中的疫情变化问题,我们认为不是单纯的预测问题,灰度预测、时间序列模型等预测模型并不能很好的适用于此。现有的SIR模型【1】是很成熟的描述一般病毒传播的数学模型,但是由于其中不包含我们需要的潜伏群,因此我们对SIR添加入潜伏群E,改进的到我们的单物种SEIR模型【2】。由SEIR模型建立微分方

问题二

针对描述埃博拉病毒在人群与猩猩群中相互传播特性以及预测人和猩猩中疫情变化的问题,本题中假设了病毒只能由猩猩传染到人。问题一中的SEIR模型虽然是对一种物种内部的传染模型,但我们认为传播模型是可以推广的,即在某一物种内部,埃博拉病毒的传播应该有着相同的规律;另外,在多个物种之间,埃博拉病毒的传播也应有类似的规律。因此,我们就以问题一的SEIR模型为基础进行改进,引入猩猩群对人群的影响因子,建立跨物种传播的SEIR模型,建立

针对此问中采取隔离与治愈感染者的措施后,要求预测疫情在人类中的发展情况、并与问题二结果作比较的问题,我们利用问题二的结论排除猩猩对人的影响后,此时病毒是在单物种内传播,适用于我们建立的单物种SEIR模型。我们利用前述数学模型,通过对比有无采取措施时的疫情变化情况,分析得出各种疫情控制措施的严格执行和药物效果的提高等对控制疫情的作用。

关键词 :单物种、跨物种SEIR模型 微分方程组 残差分析 量化分析

一 、 问题的重述与分析

1.1 问题的重述

埃博拉病毒有传染性,主要是通过病人的血液、唾液、汗水和分泌物等途径传播。病毒的潜伏期通常只有5天至10天,感染后2~5天出现高热,6~9天死亡。埃博拉病毒只有病人在出现埃博拉症状以后才具有传染性。存在似乎天生就对埃博拉免疫的人,痊愈之后的人也会对入侵他们的那种埃博拉病毒有了免疫能力。

本题希望同学们通过数学建模的方法量化埃博拉病毒的传播规律,深刻认识该病毒的危害,并分析隔离措施的严格执行和药物治疗效果的提高等措施对控制疫情的作用。

假设某地区有20万居民和3000只猩猩。人能以一定的概率接触到所有的猩猩,当接触到有传播能力的猩猩后有一定概率感染病毒,而人发病之后与猩猩的接触可以忽略。人与猩猩的潜伏期都为2周。请你根据相关信息,研究回答以下问题:

1、 建立一个病毒传播模型,动态描述病毒在“虚拟猩猩种群”中的传播,并预

测接下来的在猩猩中的疫情变化,并给出“虚拟猩猩种群”在第80周、第120周、第200周的相关数据;

2、 建立“虚拟种群”相互感染的疾病传播模型,综合描述人和猩猩疫情的发展,

并预测接下来疫情在这两个群体中的发展情况,并给出 “虚拟人类种群”在第80周、第120周、第200周的相关数据;

3、 假设在第41周,外界的专家开始介入,并立即严格控制了人类与猩猩的接

触,且通过某种特效药物将隔离治疗人群的治愈率提高到了80%。请预测接下来疫情在“虚拟人类种群”的发展情况,对比第2问的预测结果说明其作用和影响,给出“虚拟人类种群”在第45周、第50周、第55周的相关数据;

4、 请依据前述数学模型,分析各种疫情控制措施的严格执行和药物(包括防疫

药物、检疫药物和治疗药物等)效果的提高等措施对控制疫情的作用。

1.2 问题的分析

1.2.1 问题一的分析

问题一需要我们根据猩猩的发病数量和死亡数量,建立一个病毒传播模型,动态描述病毒在“虚拟猩猩种群”中的传播,并预测接下来的在猩猩中的疫情变化,这其实是一个数学上的预测问题。

针对预测问题,我们可以采用常用的时间序列模型或者灰度预测模型,但是针对病毒的传播,其自身有更加贴近实际的SIR模型,但由于SIR模型通常不考虑潜伏群的影响以及变化,而此题要求我们对潜伏群进行研究。

因此,我们改进SIR模型,加入潜伏群,建立更加适应于埃博拉病毒传播特点的SEIR模型,利用SEIR模型建立起微分方程组,利用微分方程解出的函数对题目要求的进行预测。

1.2.2 问题二的分析

第二个问题需要我们建立“虚拟种群”相互感染的疾病传播模型,描述人和猩猩疫情的发展,预测接下来疫情在这两个群体中的发展情况。这实际上是两个变量之间的关联度描述与预测问题。

本题中简化了问题,假设了病毒只能由猩猩传染到人,这就变成了单向的影响,问题一中的SEIR模型是对于一种独立的物种内部的传染模型,现在病毒从一种物种单向传播到另一物种、并在该物种内互相传播。因此,传播模型是可以推广的,即在某一物种内部,流行性疾病的传播应该有着相同的规律;另外,在两个甚至是多个物种之间,流行性疾病的传播也应该有着类似的规律。

因此,我们就以问题一中的SEIR模型为基础进行改进,初步建立两物种间病毒传播的SEIR模型。

1.2.3 问题三的分析

在第三问中,第41周后,因为外界的专家开始介入,并立即严格控制了人类与猩猩的接触,且通过某种特效药物将隔离治疗人群的治愈率提高到了80%。此时,人群与猩猩群不再有互相影响,消除了猩猩对人群的影响,所以,问题又回归到了第一问的情况,病毒在单一群体中传播的SEIR模型。 由于我们没有病毒在单一人群中传播的原始数据,这样我们无法利用SEIR模型进行预测。理论上,我们应该在第二问的基础上,计算出与猩猩接触而成为潜伏者和感染者的人数Bt,然后用已知总数减去Bt,得到病毒在人群中单一传播时的数据,以此为依据对问题进行求解。

实际上,由于易感人群数高达10E6数量级,而猩猩感染群只有10数量级,差距非常大,因此猩猩感染群对于人群的传染作用在此问题中可以忽略,我们可以直接用题目中的数据对人的SEIR群变化进行预测。

1.2.4 问题四的分析

问题四需要我们分析各种疫情控制措施的严格执行和药物(包括防疫药物、检疫药物和治疗药物等)效果的提高等措施对控制疫情的作用。那我们就需要利用没有控制措施时的疫情情况与控制时的疫情相比较,得出采取的措施对于控制

疫情的作用。

二 、 问题假设

2.1、假设附件中提供的采样数据准确无误且具有充分的代表性 2.2、假设短期内猩猩与人类的出生率死亡率对种群数量影响不大 2.3、假设死亡后的猩猩与人类不再具有感染能力

2.4、假设易感者一旦和感染者进行有效接触就会被感染成为潜伏者 2.5、假设单位时间内一个感染者能传染的数量与易感群数量成正比 2.6、假设所有自愈或治愈者不再复发

2.7、假设所有潜伏期病人都会在两周后发病

三 、 符号说明

四 、模型的建立与求解

4.1 问题一 4.1.1模型的建立

图1 猩猩中S、E、I、R关系图

在本题中,我们将未被感染过的健康猩猩作为易感群(S),易感群与感染者有效接触后,并不会立刻成为感染者,而是全都是进入潜伏期,且所有潜伏期患者都是两周后开始发病,因此,每一周的潜伏群E(t)就等于之后两周感染群的和I(t+1)+I(t+2)。同时因为死亡与自愈的猩猩不会再对别的猩猩产生影响,因此我们把死亡与自愈的猩猩同时都计入退出群。由此我们建立以下微分方程组:

dS

dt1(t)IS1(t)IS

dE(t)IS(t)E

11

dtdI

1(t)IS1(t)E1(t)I

dtdR

dt1(t)I

S(0)S0,E(0)E0,I(0)I0,R(0)R0

(0)

另外,题目中要求的是退出群中死亡与自愈分别是多少,我们的模型中的退

出群是两个之和,似乎无法得出题目要求的值,但是我们通过分析能够发现退出群中自愈量与死亡量是存在关系的,我们可以通过选取死亡群组(R1)作为参考数列,建立灰色关联度模型求解死亡群组(R1)与自愈群组(R2)的关系:

x0x0(k)|k1,2,...,n

x1x1(k)|k1,2,...,nmin|x0(t)x1(t)|max|x0(t)xs(t)|t (k)t

|x0(t)x1(t)|max|x0(t)xs(t)|t

n

r1(k)

ni1

(2)

4.1.2 模型的求解

现在我们根据建立的模型利用已知数据进行求解,我们将题目中给出的数据进行整理,得出我们建模求解需要的S、E、I、R四个群的猩猩数量(见附录1),并且利用整理出的已知数据对微分方程组中的参数进行求解: 1)潜伏者日接触率1(t):

1(t)

E/t

I(t)*S(t)

(3)

2)感染者日接触率1(t):

1(t)

I/t

I*S

(4)

3)退出率1(t):

1(t)

1(t):

1(t)=

R(t)

I(t)

(5)

4)潜伏群发病率

每日新确认病例数每日疑似累计数

(6)

由我们整理出来的数据带入以上公式,估算得:

11.085E-4,10,10.3018,10.5121

求出微分方程的参数后,我们根据整理出来的数据特点进行判断,选取第11周时的S、E、I、R值作为初值,即:

S02752,E030,I052,R0166

利用MATLAB编写程序(程序源代码见附录2),将参数以及初始值导入微分方程组(1)进行求解,预测出的S、E、I、R四个群的数量变化分别如下:

猩猩数/个

周数/周

图2 猩猩易感群(S)数量的实测值与预测值

猩猩数/个

周数/周

图3 猩猩潜伏群(E)数量的实测值与预测值

猩猩数/个

周数/周

图4 猩猩感染群(I)数量的实测值与预测值

猩猩数/个

周数/周

图5 猩猩退出群(R)数量的实测值与预测值

以上我们分别预测出了处于S、E、I、R四个群中的猩猩数量,现在我们需要将退出群(R)分解成两部分即分别还需要求解出死亡猩猩(R1)和自愈猩猩(R2)的数量。根据附件一,我们作出死亡群(R1)与自愈群(R2)的图像如下:

图6 猩猩死亡群(R1)与自愈群(R2)数量的实测值

将附件一中的死亡群与自愈群带入方程组(2)可以求得每周死亡数量与每周自愈数量的相关度:

r=0.869

说明死亡群组与自愈群组的相关度很高,即我们可以视为每周死亡数量与自愈数量的变化是同步的,因此他们各自占退出群的比例可以视为保持不变。

R1k1*R(t),R2k2*R(t) k10.6668,k21k1

通过计算我们得知第一问题的结果 :

4.1.3 结果的分析

根据模型求解结果我们发现,易感群、感染群与潜伏群的变化趋势基本一样。通常在疫情爆发时,感染群数量和潜伏群数量都是呈指数型增长,因此本题中给出的数据应该是在疫情已经爆发了一段时间以后,这时疫情已经由于生态系统自身的调节能力而得到控制,因此感染者与潜伏者数量会逐渐下降最终趋于平缓,而退出群虽然会逐渐增加,但是由于最后疫情得到控制,感染者与潜伏者趋于平缓,那么退出群也会逐渐趋于平缓直到达到上限。

从图像我们可以分析出,在开始统计的第11周-100周的时间内,每一周存在的感染者与潜伏者在逐渐减小,但是同时,因病死亡与自愈康复的猩猩数量都在较快的增长,说明这段时间疫情仍然较为严重。

在大约第120周以后,感染者、潜伏者、死亡与自愈者都已经基本达到自己的稳定值,这是说明疫情已经基本上得到控制,猩猩已经基本不再感染埃博拉病

毒。

4.1.4 模型的检验

我们从图1-图4中的实测值散点图和预测值曲线图可以看出,我们利用SEIR模型建立的微分方程组解出后得到的S(t)、E(t)、I(t)、R(t)与实际的群体数量变化非常接近,图像基本吻合,因此说明我们的预测符合实际变化情况,我们的模型可以认为是可信的。

为了更加精确的描述我们模型的准确程度,我们引入残差分析,分析实际观察值与我们建模预测的值之间的残差,用以评判我们的模型是否精准:

图7 猩猩易感群(S)的残差

图8 猩猩潜伏群(E)的残差

图9 猩猩感染群(I)的残差

图10 猩猩退出群(R)的残差

求出残差后我们发现,我们的每一组群体的残差置信区间都基本包含零点,每一组仅有少数残差的置信区间不包含零点在内,因此我们验证了我们的模型是可靠的。

4.1.5 模型的优点与缺点

一、 模型的优点

1. 首先SEIR模型进行了较为详细的分类,建立的微分方程组较为紧凑,相关性较强。因为各参数的设置合理,所以参数的设计在详细的数据的支持下和实际情况较为接近,为图形的拟合和对最终控制期的预测奠定了基础。

2.SEIR模型适用范围较广,只要数据足够地精确详细,则求出估计参数便可求解。

3.该模型虽然是控后模型,但只要人的警惕性增强,把控制期提前到疫情初期,那么只要做适当的参数修改便可以作为疫情发展的全过程的预测模型。

二、 模型的缺点

1.首先是SEIR模型没有考虑年龄结构层次对疫情的影响。因为根据医学研究表明,年龄较小和较大者极易感染病毒,而青壮年由于身体机能强免疫力高不易感染病毒。

2.SEIR模型是基于我们将α,β等描述感染能力的参数当做定值,即感染者潜伏者感染其他人的能力不会发生变化而成立,而实际上病毒的感染能力会随着接触水平等因素变化而变化。

4.2 问题二 4.2.1模型的建立

事实上病毒在两种物种之间相互传播时传播特点非常复杂,但本题已经限定了埃博拉病毒从物种A(猩猩)向物种B(人)的传播,而忽略B对A的影响(虚线部分)。经过这样简化,我们可以得到下图(图11)中的传播关系:

图11 猩猩与人的S、E、I、R关系图

由第一问中的SEIR模型,我们经过可以改进可以得到病毒在两物种之间传播时的SEIR模型,并由此建立新的微分方程组:

dAs

dt1(t)AiAs1(t)AiAs

dAe(t)AiAs(t)Ae

11

dt

dAi(t)AiAs(t)Ae(t)Ai

111

dtdAr1(t)Aidt

(7)

dBt

dt(t)AiBsp(t)Btq(t)Bt

dBs(t)BiBs(t)BiBs(t)

22

dtdBe2(t)BiBsp(t)Bt2(t)Be

dt

dBi

2(t)BiBs2(t)Beq(t)BtdtdBr

2(t)Bidt(t)p(t)q(t)

(8)

同样,还需要求解出死亡人群(Br1)、治愈人群(Br2)及隔离治疗(Br3)的

数量,我们选取死亡群组(Br1)作为参考数列建立灰色关联度模型求解死亡群组

与治愈群组(Br2)和隔离治疗(Br3)的关系:

x0x0(k)|k1,2,...,n

x1x1(k)|k1,2,...,nmin|x0(t)x1(t)|max|x0(t)xs(t)|t (k)t

|x0(t)x1(t)|max|x0(t)xs(t)|t

n

r1(k)

ini1

(9)

4.2.2 模型的求解

现在我们根据建立的模型利用已知数据进行求解,我们将题目中给出的数

据进行整理,得出我们建模求解需要的S、E、I、R四个群的人数(见附录3),并且利用整理出的已知数据对微分方程组中的参数进行求解: 1) 潜伏者日接触率2(t):

2(t)

Be/t

Bi(t)*Bs(t)

(10)

2)感染者日接触率2(t):

2(t)

Bi/t

Bi*Bs

(11)

3)退出率2(t):

2(t)

Br(t)

Bi(t)

(12)

4)潜伏群发病率2(t):

2(t)=每日新确认病例数每日疑似累计数

(13)

5)B物种与A物种接触群日增长率λ(t):

(t)

p(t)

B物种接触群增长数量

Bs∗Ai

(14)

6)B物种接触群到潜伏群的比率p(t):

B物种由接触者成为潜伏者的日新增数量

B物种接触者日累计数量

(15)

7)B物种接触到感染群的发病率

q(t)

B物种由接触者成为感染者的日新增数量

B物种接触者日累计数量

(16)

由我们整理出来的数据带入以上公式,估算得:

2p2-1.23907,27.51951,21.49931,q1.23907

求出微分方程的参数后,我们根据整理出来的数据特点进行判断,选取第11周时的Bs、Be、Bi、Br值作为初值,即:

Bs0,Be0,Bi0,Br0

利用MATLAB编写程序(程序源代码见附录4),将参数以及初始值导入微分方程组(8)进行求解,预测出人的S、E、I、R四个群体的数量变化分别如下:

图12 人群中易感群(Bs)的实际值与预测值

图13 人群中潜伏群(Be)的实际值与预测值

图14 人群中感染群(Bi)的实际值与预测值

图15 人群中退出群(Br)的实际值与预测值

将附件3中的死亡群与自愈群带入方程组(2)可以求得每周死亡人数与每周自治愈人数、每周隔离治疗人数的相关度:

r1=0.832;r2=0.256

说明死亡群组与自愈群组的相关度很高,即我们可以视为每周死亡人数与治愈人数的变化是同步的,但是死亡人数与隔离人数之间并无明显关系。

而我们通过计算发现,隔离治疗人数与感染群人数有非常高的相关度,因此我们可以用感染群人数估计隔离治疗的人数,再利用退出群中减去隔离治疗的部分计算死亡人数与治愈人数:

Br3k*I

Br1k1*(R(t)Br3),Br2k2*(R(t)Br3)

k0.710547032948828

k10.7216,k21k1

通过计算我们得知第二个问题的结果 :

表2 “虚拟人类种群”群体数量预测结果(单位:个)

4.2.3 结果的分析

根据模型求解结果我们发现,在猩猩对人群的单向传染中,人群的易感群、感染群与潜伏群数量与在单个物种中传播时的变化趋势基本一样。且利用模型求解出的结果和实际数据仍然非常吻合,说明我们的SEIR模型经过修改,仍然适用于简化情况下病毒在两个物种之间的传播特点。

从图像我们可以分析出,在开始统计的第11周-100周的时间内,每一周存在的感染者与潜伏者在逐渐减小,但是同时,因病死亡、治愈后康复以及隔离治疗的总人数在较快的增长,说明这段时间疫情仍然较为严重。

在大约第120周以后,感染者、潜伏者、死亡、治愈以及隔离的人数都已经基本达到自己的稳定值,这是说明疫情已经基本上得到控制,人群已经基本不再大规模感染埃博拉病毒。

4.2.4 模型的检验

我们从图12-图15中的实测值散点图和预测值曲线图可以看出,我们利用SEIR模型建立的微分方程组解出后得到的S(t)、E(t)、I(t)、R(t)与实际的群体数量变化非常接近,图像基本吻合,因此说明我们的预测符合实际变化情况,我们的模型可以认为是可信的。

为了更加精确的描述我们模型的准确程度,我们引入残差分析,分析实际观察值与我们建模预测的值之间的残差,用以评判我们的模型是否精准:

图16 人群中易感群(Bs)的残差

图17 人群中潜伏群(Be)的残差

图18 人群中感染群(Bi)的残差

图19 人群中退出群(Br)的残差

求出残差后我们发现,我们的每一组群体的残差置信区间都基本包含零点,每一组仅有少数残差的置信区间不包含零点在内,因此我们验证了我们的模型是可靠的。

4.2.5 模型的优点与缺点

一、模型的优点

1. 首先SEIR模型对于模拟埃博拉病毒在单种生物中传播时具有非常好的效果,预测疫情发展时非常精确

2.我们对单种生物中的SEIR模型进行改进,添加入A种群对B种群的影响因子,贴近实际,使得我们的模型对于更为复杂的跨物种传播仍然适用。 二、模型的缺点

1. SEIR模型自身的缺点在此次改进中不能得到修正,仍然忽略了很多诸如年龄等对于疫情传播有影响的因素。

2.改进后的SEIR模型虽然预测后仍然很精确,残差也非常小,但是我们发现在预测潜伏群体时出现了与常识相违背的结果,这可能是由于我们的参数计算时数据太少而且偶然因素较多导致,这也说明了我们模型的抗干扰能力不够强,实用性也仍然不够广。

4.3 问题三 4.3.1模型的建立

我们在分析中已经得出,我们可以忽略猩猩感染群对于人类的感染作用,因而我们直接利用题目中的数据带入我们一问中建立的单一物种间的病毒传播SEIR模型,建立微分方程组进行求解:

dBs

dt3(t)BiBs2(t)BiBs

dBe(t)BiBs(t)Be

33

dt

dBi(t)BiBs(t)Be(t)Bi

333

dtdBr3(t)Bidt

(17)

在第三问中,题目中要求解出死亡人群(Br1)、治愈人群(Br2)及隔离治

疗(Br3)的数量,我们在第二问中已经通过灰色相关度分析得出死亡群组与自愈群组的相关度很高,即我们可以视为每周死亡人数与治愈人数的变化是同步的;同时,隔离治疗人数与感染群人数有非常高的相关度,因此我们可以用感染群人数估计隔离治疗的人数,再求解死亡群组与感染群组人数。

4.3.2 模型的求解

根据模型解得在41周以后通过采取隔离和治疗后,疫情在“虚拟人类种群”的发展情况:

图20 人群中易感群(Bs)的实际值与预测值

图21 人群中潜伏群(Be)的实际值与预测值

图22 人群中感染群(Bi)的实际值与预测值

图23 人群中退出群(Br)的实际值与预测值

根据建模解出的结果分别计算隔离治疗群(Br3)、累计治愈群(Br2)以及累计死亡群(Br1)人数如下:

Br3k*I

Br1k1*(R(t)Br3),Br2k2*(R(t)Br3)

k0.710547032948828

k10.7216,k21k1

通过计算我们得知第三个问题的结果 :

4.3.3 结果的分析

根据模型求解结果我们发现,在去除猩猩种群对人群的影响以后,病毒单独在人群中进行传播时的易感群、感染群与潜伏群数量与在猩猩中传播时的变化趋势总体上仍然一直。

从图像和数据我们可以分析出, 在第41周开始隔离开猩猩与人群并且提高治疗药效后,潜伏群与感染群下降的速度明显快于在第二问没有采取任何措施时

的下降速度。

这说明采取隔离措施时,减少了人的接触被感染率,抑制潜伏群和感染群的增长,直接体现在加快了潜伏群与感染群趋于最终稳定的速度。

而治疗药物效果的提高,能及时检测人群中的潜伏群,若能采取及时的隔离措施,则能够大大降低其与易感群的有效接触,使传染得到控制,减少了潜伏人群和处于发病状态的数量以及死亡数。

4.3.4 模型的检验

我们采用与第一问一样的模型求解,可以发现最终的预测曲线与实际值之间同样非常吻合,我们同样作出残差图对模型的精确度进行检验分析:

图24 人群中易感群(Bs)的残差

图25 人群中潜伏群(Be)的残差

图26 人群中感染群(Bi)的残差

图27 人群中退出群(Br)的残差

4.3.5 模型的优点与缺点

由于我们在此问中采用了与第一次完全一致的SEIR模型,因此模型具有SEIR模型自身具有的诸多优缺点,在此不再赘述。

需要在此说明的是,在第三问的处理中,我们根据分析得出猩猩对于人的传染影响很小,因此忽略了猩猩群体对人群体的直接影响是对问题进行了简化,但是实际上在41周将猩猩与人群隔离开以前,猩猩的存在会对人群的S、E两个群体数量产生直接影响,且通过S、E两个群体对I、R两个群体产生间接影响,虽然影响很小,但是仍然可能造成我们的简化模型产生较大误差,这是我们模型最大的缺点。

4.4 问题四

通过将第一二问没有采取任何控制与防治措施时的疫情与第三问采取隔离和投放药物的疫情相比较,我们可以得出如下结论:

投放防疫药物、采取隔离措施时,能减少人的接触被感染率,即减少有效接

触,即抑制潜伏群的增长,在我们模型中的直接体现是降低了有效接触的数量,从而降低感染群的峰值。

检疫药物效果的提高,能及时检测人群中的潜伏群,若能采取及时的隔离措施,则能够大大降低其与易感群的有效接触,使传染得到控制,减少了潜伏人群和处于发病状态的数量以及死亡数。

而治疗药物药效的提高,直接作用于患病人群,能够降低其死亡率。

故防疫与检疫药物的效果提升和相关措施的严格执行,能够对疫情的控制起到决定性的作用,能够减少疫情的扩散;而治疗药物的效果的提升则降低了死亡率,控制疫情的死亡人数。

五、参考文献

[1]百度百科,http://baike.baidu.com/link?url=TXifg9YS2W-h23d9kE9BemlthlBguG7KaSJwdRSBjeRXsxmAt-9v7dfT926bUSPO3htwKpaJ-ONYARiSaVaBHa,2015年5月1日

[2]周广发、陈月华,SARS传播模型的探讨,2015年5月1日

六、附录

6.1 附录1

6.2 附录2

问题一微分方程组(1)MATLAB求解程序代码

function dy=seir(t,y) dy = zeros(4,1);

α=0.000107593196449557; β=0;

w=0.301185431268702; =0.512238823647188;

dy(1)= -y(3)*y(1)* α-β*y(1)*y(3); dy(2)= α*y(1)*y(3)-y(2)*;

dy(3)= β*y(1)*y(3)+*y(2)-y(3)*w; dy(4)=y(3)*w; end

6.3 附录3

6.4 附录4

问题二微分方程组(7)(8)MATLAB求解程序代码

function dz=qq(t,z) dz=zeros(8,1);

a=0.000107593196449557; b=0;

w=0.301185431268702;

q=0.506343491;

b2=7.51950955852228E-06; q2=-1.23907471853569E-07; w2=1.49931331589255;

dz(1)= -q2*z(7)*z(1)-b2*z(1)*z(3); dz(2)= q2*z(1)*z(7)-z(2)*q2;

dz(3)= b2*z(3)*z(1)+q2*z(2)-z(3)*w2; dz(4)=z(2)*w2;

dz(5)= -z(7)*z(5)*a-b*z(5)*z(7); dz(6)= a*z(5)*z(7)-z(6)*q;

dz(7)= b*z(5)*z(7)+q*z(6)-z(7)*w; dz(8)=z(7)*w; end

范文五:埃博拉病毒传播模型及规律预测

摘要:埃博拉病毒是一种能引起人类和灵长类动物产生埃博拉出血热的烈性传染病病毒。文章以2014年西非疫区为参照,建立虚拟环境的常微分方程组,利用四阶龙格―库塔法求解其数值解,具体通过C语言程序设计实现,并据此研究埃博拉病毒的传播规律,分析隔离措施的严格执行和药物治疗效果的提高等措施对控制疫情的作用。

关键词:数学建模;埃博拉病毒;常微分方程组数值解;四阶龙格―库塔法;西非疫区 文献标识码:A

中图分类号:O175 文章编号:1009-2374(2016)04-0194-03 DOI:10.13535/j.cnki.11-4406/n.2016.04.096

1 模拟真实环境

埃博拉病毒的自然宿主虽尚未最后确定,但已有多方证据表明猴子及猩猩等野生非人灵长类动物有埃博拉感染现象。该病毒的传播途径分为人畜传播、人人传播两种。2014年,在几内亚、塞拉利昂和利比里亚等国,许多受埃博拉病毒影响的人口都以丛林肉为重要的蛋白质和营养物质来源,与丛林中动物接触频繁。这为人畜之间的病毒传播创造了条件。

我们现假设两个感染埃博拉病毒的虚拟种群:即某地区内的20万居民和3000只猩猩。人能以一定的概率接触到所有的猩猩,当接触到有传播能力的猩猩后有一定概率感染病毒,而人发病之后与猩猩的接触可以忽略。人与猩猩的潜伏期都为2周。并在出现疫情41周后模拟外界医疗力量的介入,使得人类与猩猩不再发生接触,且隔离治疗人群的治愈率提高到80%。模拟数据详见附录。

2 建立数学模型

2.1 模型假设

(1)依据人或猩猩的健康状态,将人或猩猩划分为健康者、埃博拉感染者(也称患病者)、退出者(含自愈者、死亡者);(2)自然封闭条件下,猩猩无自然迁移,故无病源的流入、流出,种群数量不变。人类数量庞大,在无大规模迁移的情况下,认为人类数目为一定值,保持不变;(3)健康者中不包含退出者;(4)人和猩猩自愈后二度感染的概率均为0,人被治愈后二度感染的几率为0;(5)不存在有效免疫药物可使人对埃博拉病毒产生免疫,同时猩猩对病毒也不免疫;(6)人的传染途径有人传染人、猩猩传染人两条。两条途径的传染率并不相同,分别假设为传染率C1和传染率C2。C1猩猩与猩猩之间传染途径只有猩猩传染猩猩一条,假设猩猩之间的传染率为C0;(7)患病人无法传染患病猩猩;(8)41周外界介入后,猩猩与人的传播途径切断,隔离患者的治愈率提高到80%,同时未被隔离的患者治愈率不变。

2.2 符号说明

符号说明如表1所示:

3 模型的建立与求解

3.1 数据处理

根据累计死亡个体数,求得每周死亡个体数。同理,根据累计自愈个体数,求得每周自愈个体数。由每周仍处于发病状态的个体数加本周自愈个体数和本周死亡个体数得到每周总患病个体数。

由每周总患病个体数比总体数目得到每周患病个体在总体中的比例A(t);由相邻两周每周患病数相减得到每周新增患病数。由总体个数减去新增病例累计和获得健康个体数,并由此得到健康个体在总体中的比例J(t);由累加自愈治愈人数与累加死亡人数得到退出者总数量,并由此得到退出者在总体中的比例T(t)。

由于埃博拉病毒的潜伏期是两周,所以任意一周的新增病例是两周前处于患病状态的个体传染的。由新增病例数比两周前处于患病状态的个体数得到该周埃博拉病毒传染率,由此计算出每一周的传染率。通过MATLAB绘图,我们得到其近似曲线为一条平行于X轴的曲线,所以通过加权平均求得平均传染率。

因为人类最初患病个体不可能为人类传染所致,所以两周内出现的患病者必为由猩猩传播而来的。最初两周,人每周新增的患病个数除以处于两周前患病状态的猩猩个数得到猩猩与人之间平均传染率C1。

两周之后人患病可由猩猩和人传播两种途径导致,猩猩每周处于患病状态的数量和C1已知,所以每周由猩猩传播导致人患病的数量可求出,用每周新增患病人数减去每周猩猩传播导致人患病的数量,即每周由人传染导致的患病数量。由每周人传染导致的患病数量除以两周前未被隔离处于患病状态的人数可得到每周埃博拉人与人之间传染率C1,通过MATLAB画出C2的图像,可以发现其图像为平行于x轴的曲线,可通过加权平均求出人与人之间平均传染率C2。

死亡率是由本周死亡个数比本周总病例数得到。我们由此求得每一周的死亡率。通过MATLAB绘图我们得到一条同样近似平行于X轴的曲线,所以通过加权平均,求得疫情稳定后平均死亡率。

每周处于未隔离的患者人数处于每周的处于患病状态的总人数可得到每周的未隔离率G,其图像为平行于x轴的曲线,通过加权平均求出平均未隔离率G。

同理,我们还得到疫情稳定后的平均自愈及治愈率。我们近似地认为,在病情爆发后不久,即疫情稳定后,周感染率、周死亡率、周治愈自愈率、周未隔离率都是常值。

3.2 埃博拉病毒的传播模型

由假设知,猩猩患病只能由猩猩传播。每个患病猩猩每周可使得只的健康猩猩变为患病猩猩,由患病猩猩数量得每周共有只健康猩猩被感染。即患病猩猩的增加率,又因为每周自愈的猩猩数目为,死亡的猩猩数目为,所以猩猩患病数随时间变化满足:

同理,我们得到人类的传播模型(由前述,所有脚标为2的符号均为人类相关数据,脚标为1的符号为与猩猩相关数据):

3.3 模型求解

通过联立方程组和数据处理,我们使用四阶龙格―库塔法分别求出人和猩猩群体中健康者和患病者的比例的数值解。

通过C程序设计编译程序求解模型所用常微分方程组的数值解。该程序在前四十组解得检验中拟合程度极高,故由此得到较为可靠的预测数据。C语言程序代码详见附录。

数据拟合如下:

3.4 建立使用免疫药物后的模型

未被隔离患者的治愈率和被隔离患者的治愈率加权平均后得到患者的平均治愈率Zh。1-Zh为死亡率与未被治愈率之和。默认死亡率与未被治愈率权重不变。在(1-Zh)中,可以算出死亡率的权重和未被治愈率的权重。在严格控制人类与黑猩猩接触并使用特效药后的数学模型如下:

将之与隔离前模型对照后发现,特效防疫药物的使用极大地提高了治愈率Zh,快速地降低了患病数。而隔离黑猩猩的措施将黑猩猩传染致病人数降为0。

综上,两种措施都有效地预防了疾病的进一步扩散,抑制了疾病的传播,使患病人数的增长率由正变负,从而导致患病数在短期内大量且持续减少。

4 模型的评价与推广

4.1 模型优势

(1)种群数量较小时,通过求解比例的变化得出结果较为精确;(2)可以动态地描述种群发病率、死亡率、自愈率、治愈率、传染率等多种特征量;(3)四阶龙格―库塔法通过数值求解常微分方程组,很好地拟合了给定的原始数据。

4.2 模型缺陷

(1)由于通过比例而非数量求解,所以当种群数量较大且比例变化不明显时误差较大;(2)认为两个“虚拟种群”内部个体总数在随时间变化时基本不变,未考虑个体总数的时间变化率。

5 附录

5.1 模拟数据

模拟数据如表2所示:

5.2 C语言程序代码

参考文献

[1] 李信真.计算方法[M].西安:西北工业大学出版社,2013.

[2] H.Nishiura,G.Chowell.EARLY TRANSMISSION DYNAMICS OF EBOLA VIRUS DISEASE(EVD)[J].WEST AFRICA,2014,(8).

[3] 甄西丰.实用数值计算方法[M].北京:清华大学出版社,2006.

[4] 刘来福,何青.用Maple和MATLAB解决科学计算问题[M].北京:高等教育出版社,1999.

(责任编辑:王 波)

范文六:埃博拉病毒传播问题

一、 问题重述

埃博拉病毒至今已多次爆发,很难根除。本题希望我们通过建模来量化病毒

传播规律以达到了解疫情和实施有效措施控制疫情的目的。

1、仅针对猩猩群体建立病毒传播模型,根据已知数据描述并预测疫情变化,完成第80、120、200周的数据。

2、对人和猩猩的群体建立模型,综合描述并预测人和猩猩两个群体疫情的发展,完成第80、120、200周的数据。

3、从第41周开始,在实施严格控制人与猩猩的接触和将人群治愈率提高至80%两措施的情况下,对比第2问的预测结果说明措施的作用和影响。完成第45、50、55周数据。.

4、据前分析各种疫情控制措施的严格执行和药物效果的提高等措施对控制疫情的作用。

二、 问题分析

题目按照研究猩猩单独群体、人和猩猩混合群体、隔离措施和药物治疗对病

毒传播的影响的顺序层层深入,逐步展开对控制疫情的探索研究。

1、问题一要求建模描述并预测病毒在猩猩种群中的传播。实际需要找到几个已知数据类即发病个体数、自愈个体数和死亡个体数之间的关系。观察到已知数据中发病猩猩个数、自愈猩猩个数、死亡猩猩个数随时间单调变化,所以考察各变量从时刻t到时刻t+t的增量,然后得到变量关于t的微分方程,编程进行数值计算。通过已知数据验证模型的可靠性,再以此为基础进行预测。其中,要合理处理潜伏群体和处于发病状态个体数量的关系。

2、问题二建立在问题一的基础上,将人群进行分类。考虑到猩猩感染人而人不感染猩猩,猩猩群体的传播规律不受影响,人类患病者的感染源为患病猩猩和患病的人,这些患病人类中一部分会被隔离治疗,隔离群体则有相对应的治愈率和死亡率,借鉴问题一的想法,找到问题二中给定量和假设量他们之间的关系。然后进一步建立被分类人群(易患病的人、潜伏期的人、已患病的人、患病被隔离的人、隔离被治愈的人、自愈的人、隔离后死亡的人、未被隔离死亡的人)的数量关系。

3、问题三中,采取隔离措施后,人的群体和猩猩的群体相互独立,即猩猩不再将病毒传染给人类。同时,被隔离人群的治愈率提高到了80%。所以对比于第二问,忽略患病猩猩对人的传染,同时被隔离群体的治愈率为一给定值0.8 。可将问题二的模型修改直接利用。对措施产生的作用进行探究应从处于发病状态个体的个数和死亡个体的个数的变化入手,考虑干预手段的大小对疾病预防和控制的影响。

4、注意到前面模型中对患病人群的隔离率没有相应的措施,所以,问题四中应重点分析隔离率对疫情控制的影响,可具体考察与隔离率相关的式子并进行分析。

三、 模型假设

1、埃博拉病毒变异可通过呼吸传播的情况不会发生; 2、该地区人和猩猩均无先天对埃博拉病毒免疫的个体; 3、根据埃博拉病毒的特点,原发病体与其所感染的患病体数量相比非常小,因此可假设患病体都是由传染而得病的。

4、处于潜伏期的个体必定发病;

四、符号说明

t:单位时间(周);

S1(t): 易患病猩猩的数量; E1(t): 潜伏期猩猩的数量; I1(t): 已患病猩猩的数量; R1(t): 自愈猩猩的数量; D1(t): 死亡猩猩的数量; a:猩猩的自愈率; b:猩猩的死亡率

c:猩猩的有效接触率(易患病猩猩一接触患病猩猩则患病称为有效接触);

S2(t): 易患病的人的数量;

E2(t): 潜伏期的人的数量; I2(t): 已患病的人的数量; G (t) : 患病被隔离的人; R*(t): 隔离被治愈的人的数量; R2(t): 自愈的人的数量; D2(t): 死亡的人的数量; d:人的自愈率;; e:人的死亡率;

f:人的有效接触率;

. k:人的隔离率;

g:人的隔离治愈率; 1-g:人的隔离死亡率;

五、模型的建立与求解

问题一:

1、动态描述:

图1是对前40周猩猩群体中个体发病、自愈、死亡数量变化趋势的折线图形描述,先纵向分析三者各自的变化趋势,可以看出:前40周,死亡个体数和自愈个体数均随周数增加,发病个体数在前6周波动,后随周数减少,前40周发病个体数和自愈个体数增加的趋势有所减缓;后横向分析,可以看出:前40周,死亡个体数和自愈个体数同时增加,具有一致性,但不确定这种一致性会保持多久,而他们与发病个体数有相反的变化趋势。 从以上分析可以得出下面结论:

在不具备集体隔离和治疗能力的猩猩群体内,发病和死亡的趋势不会无限进行下去,猩猩们的自愈能力会扭转颓势。

图1 2、预测:

猩猩的群体可分为5类,分别为易患病猩猩、潜伏期猩猩、患病猩猩、自愈类猩猩、死亡类猩猩。由假设得:潜伏期猩猩2周后均变为患病猩猩。

自愈率=自愈猩猩数/(自愈猩猩数+死亡猩猩数)

a= R1 /( R1+D1 ), 已知数据中R1与 形表现为图2:

( R1+D1 )关系的图

图2

容易得到 a=0.33;

对于易患病猩猩,t时刻的数量

S1(t)是连续、可微函数,并且每周每

c,由于易患病猩猩随时间是减少

个患病猩猩有效接触的易患病猩猩的数量为的,所以考察t

到t+t时间内易患病猩猩数量的增加量,就有

S1(t+t)- S1(t)= - c I1S1t

再设t=0时有S0 个易患病猩猩,即得微分方程

dS1

c I1S1 (1)dt

依照(1)式的获得方式,即得如下微分方程

dE1

cI1S1E1t(dt

2 ) (2)

dI1

E1(t2)aI1bI1 (3) dtdR1

aI1 (4) dtdD1

bI1 (5) dt

由于方程(1)~(5)无法求出解析解,所以利用已知附件数据,通过MATLAB

软件进行数值计算,计算所得数据的图表形式如下:

图3

4

图5

然后将指定周数的数据填入下表:

问题二:

混合群体中类似于猩猩群体,人的群体分为8类:易患病的人、潜伏期的人、已患病的人、患病被隔离的人、隔离被治愈的人、自愈的人、隔离后死亡的人、未被隔离死亡的人。由于人发病后与猩猩的接触可以忽略,所以只有猩猩会感染人而人不会感染猩猩。则猩猩群体病毒传播的规律不发生改变,人的感染源有患病猩猩和患病的人。则人的群体的对应微分方程为:

dS2

 (If1+I2)S1 (6) dt

dE2

f(I2+I1)S2E2(t2) (7) dt

dI2

E2(t2)(1k)(dI2eI2)kI2dt

(8)

dR2

(1k)I2d (9) dtdG

I2k (10) dt

dR*

kI2g (11)

dt

dD2

(1k)I2ekI2(1g) (12) dt

由于在问题二中,模型未知量较多和前40周数据有所波动造成了模型求解的复杂,所以在使用已知数据时,将处于发病状态人数数据的前8周看做第一阶段,以后数据为第二阶段;将被隔离人数数据的前5周看做第一阶段,以后看做第二阶段;另外两组数据完整使用。这样在编程求解时初值做相应改变。 通过MATLAB软件进行数值计算,计算所得数据的图表形式如下:

图6

图7

图8

图9

问题三: 1、预测:

问题三中,采取隔离措施后,人的群体和猩猩的群体相互独立,即猩猩不

再将病毒传染给人类。同时被隔离群体的治愈率为一给定值

0.8

g=0.8,1-g=0.2 。(7)~(12)式修改为:

dS2

 fI2S1 dt

dE2

fI2S2E2(t2) dt

dI2

E2(t2)(1k)(dI2eI2)kI2 dtdR2

(1k)I2d dt

dG

I2kdt

dR*

0.8k2I dt

dD2

(1k)I2e0.2kI2 dt

可利用修改过的问题二的模型进行预测。即可完成下表:

2、作用和影响说明:

严格控制患病猩猩与人类的接触后,人的感染源减少了近40%,使人类处于潜伏期和发病状态的数量大幅减少;药物治疗提高治愈率至80%后,隔离治疗人群的死亡率下降了62%,总的患病人群的死亡率下降了41%,对疫情的控制

效果明显。说明干预手段的程度的大小对疾病控制的影响显著。

问题四:

前述数学模型中疫情控制措施有:隔离措施(切断患病猩猩与人的接触、将患病人隔离)和提高治愈率的措施。

前面实现了严格控制患病猩猩与人的接触,并将隔离人群的治愈率提高至80%,对疫情的控制效果比较明显。但患病人群的隔离率k未经改善,观察问题二中(8)式并进行变形得下式:

dI2E2(t2)I2(de)kI2(de1)dt

当提高隔离率时,处于发病状态的个体数量会减少,因病死亡人数也会减少,对控制疫情有积极效果。

模型的建立帮助我们搞清楚了现有干预手段究竟有多大效果,预测表明,隔离措施的严格执行和药物治疗是控制疫情的最好手段。

六、模型的优缺点评价及改进

模型优点:

1、模型变量全面,描述和联系性强,可以将问题描述清楚并就问题给出可参考答案。

2、模型具有灵活性和适应性,如解决问题二时,取用部分数据,模型仍然成。

模型缺点:

1、变量多,模型实现时遇到了较多困难。

2、在进行函数拟合时,有较多点在函数曲线外,误差较大。

3、细节处理不到位。

模型改进:

可以先只找S和I的关系,列微分方程求解,再列E、R、D的微分方程。这样就将一个大问题转化为两个较容易的问题,相当于减少了变量个数,而且更能考虑他们之间的细节关系。

参考文献:

【1】姜启源,谢金星,叶俊 . 数学模型(第三版)【M】. 北京:高等教育出版社,2003.

【2】周后卿,徐幼专:埃博拉病毒感染数量的一个数学模型 .

【3】杨玉华:SARS传染模型的研究及实证分析 .

11

范文七:狗会感染和传播埃博拉病毒吗?

狗会感染和传播埃博拉病毒吗?

狗会感染病毒,但是否传播尚无定论!

一项由美国疾控中心(CDC)开展的研究,调查狗是不是真的能感染埃博拉病毒。研究检查了2001年到2002年加蓬暴发埃博拉疫情期间,当地狗体内的埃博拉病毒抗体的普遍程度,结果显示狗能够感染病毒,但是不会表现出任何症状,而感染最终也会自行消失。

研究者最终推测:“狗可能是人类埃博拉病毒爆发以及扩散的源头”,即动物体内的病毒被清除之前的一小段时间内,通过唾液、尿液和粪便中携带的病毒颗粒传播疾病;被狗舔、咬或者给狗梳毛等行为将病毒传给人类。

对于人是否会遭狗传染,医学人员看法有分歧。美国兽医学会发言人莎伦·柯蒂斯·格兰斯科格如是说:“我认为有这种可能。但是,如果狗不靠近或啃食病死的动物或人的尸体,这种可能性不大。”

美国北卡罗来纳州立大学兽医学者彼得·考恩博士则持相反看法。考恩曾在各地就动物传播疾病的风险提供咨询。他说:“从来没有记载显示狗能传播埃博拉病毒。这显然不是非洲本次疫情中(病毒)传播的主要途径。”

美国疾病和预防控制中心主任托马斯·弗里登说:“显然,我们将研究所有可能性。我们尚未把这认定为传播途径。”

相关阅读---一条狗的命运!

尽管这一问题没有定论,西班牙一名染病女护士的宠物狗则面临“死刑”。马德里地方政府7日接到法庭命令,将对这条狗处以安乐死。政府方面的理由是,“没有科学证据”显示狗不会传播埃博拉病毒。

这名女护士9月30日出现发病症状,是首个在非洲以外感染埃博拉病毒的病例。她在马德里的卡洛斯三世医院工作,护理过两名在非洲染病的神职人员。两名患者分别于8月12日和9月25日死亡。这名护士的丈夫和两名同事也被隔

离观察。

西班牙动物主题脱口秀节目主持人、兽医卡洛斯·罗德里格斯说,女护士和丈夫不忍心狗遭处死。男主人从医院给他发信息,拜托他暂时收养这条狗。但是因为法院有命令,罗德里格斯表示爱莫能助。“他又哭着拜托我,至少要确保这条狗没有痛苦地死去。”

西班牙动物权益组织不满当局不给狗做体检或隔离就直接处死。“这明显是反应过度。”兽医学者考恩说。

这名护士的丈夫与宠物狗Excalibur的合影

哪些其他生物会感染埃博拉病毒?

埃博拉病毒是人畜共患病,这意味着它可以在不同物种间传播。感染该病毒后,受害最大的是灵长类动物,包括大猩猩、黑猩猩、猴及人类。其他已知被感染的动物包括非洲果蝠、羚羊、豪猪、兔子、猪等。现今,还没有猫会被感染的报道。

埃博拉病毒如何入侵人类社会?

美国疾控中心的主任汤姆·弗里登(Tom Frieden)医生在周二的新闻发布会上表示:“我们知道在非洲的农村地区,埃博拉可以感染哺乳动物。事实上,病毒就是靠这种方式传播的,从大概是蝙蝠传到了生活在丛林中的其他动物身上,而当地人会捕猎这些动物。”

到目前为止,已经在许多丛林动物体内发现了埃博拉病毒,包括丛林野猪、啮齿动物、豪猪以及羚羊。任何感染了埃博拉并死去的动物都有可能将病毒传播给狩猎者,或是任何以丛林肉为食的人们。

利比里亚当地的告示牌,警告人们不要猎捕和食用野生动物

远离恐怖------现阶段我们能做些什么?

科学家还没能确定埃博拉病毒的天然宿主,即自然状态下就会被病毒感染并且成为病毒扩散的主要源头的动物,但是传播病毒则是另一码事。

英国兽医协会主席约翰·布莱克威尔(John Blackwell)指出:“既然这种病毒目前仅在灵长类动物(包括人类)和少数几种野生动物体内被发现,那么传播的主要途径应该还是人传人。”

他补充道,狗身上的发病过程和它们在传播病毒过程中的作用还未确定,但是对与确认感染或疑似感染的人接触过的动物实施严格的隔离措施,还是很合理的预防方法。

感谢您阅读本篇文章,如果感觉文章不错,可分享至朋友圈;想得到相关文章,请关注“宠儿香微信服务号”,我们会定期把文章推送给您!

范文八:埃博拉病毒传播问题

第一部分 问题重述

第二部分

问题分析

第三部分

模型的假设

第四部分

定义与符号说明

第五部分

模型的建立与求解

第六部分

模型的评价与推广

第七部分

参考文献

第八部分

附录

(www.wenku1.com)一、 问题重述

假设某地区有 20 万居民和 3000 只猩猩。 人能以一定的概率接触到所有的 猩猩, 当接触到有传播能力的猩猩后有一定概率感染病毒,而人发病之后与猩猩 的接触可以忽略。研究人员统计了前 40 周人类和猩猩的发病数量和死亡数量等 信息(见附件一,附件二) ,相关问题如下; (1)根据猩猩的发病数量和死亡数量,建立一个病毒传播模型,动态描述 病毒在“虚拟猩猩种群”中的传播,并预测接下来的在猩猩中的疫情变化,并给 出“虚拟猩猩种群”在第 80 周、第 120 周、第 200 周的相关数据; (2)建立“虚拟种群”相互感染的疾病传播模型,综合描述人和猩猩疫情 的发展,并预测接下来疫情在这两个群体中的发展情况,并给出 “虚拟人类种 群”在第 80 周、第 120 周、第 200 周的相关数据; (3)假设在第 41 周,外界的专家开始介入,并立即严格控制了人类与猩猩 的接触,且通过某种特效药物将隔离治疗人群的治愈率提高到了 80%。预测接下 来疫情在“虚拟人类种群”的发展情况, 对比第 2 问的预测结果说明其作用和影 响,给出“虚拟人类种群”在第 45 周、第 50 周、第 55 周的相关数据; (4)依据前述数学模型,分析各种疫情控制措施的严格执行和药物(包括 防疫药物、检疫药物和治疗药物等)效果的提高等措施对控制疫情的作用。

(www.wenku1.com)二、 问题分析

(一)问题 1 的分析 建立一个病毒传播模型, 动态描述病毒在“虚拟猩猩种群”中的传播,有助 于政府和相关医疗机构预测接下来在猩猩中的疫情变化情况,科学预测疫情在 “虚拟猩猩种群”中的发展以及可能对“虚拟人类种群”带来的影响。 题中给出了“虚拟猩猩种群”的有关数据 (截至本周末处于发病状态的数量、 截至本周末总共死亡个数、截至本周末总共自愈个数) ,要求我们预测接下来在 猩猩种群中的疫情变化,分别给出在第 80 周、第 120 周、第 200 周“虚拟猩猩 种群”中潜伏群体、处于发病状态、累计自愈、累计因病死亡的数据。 由以上原因,我们可以建立一个 SLIRD 模型,利用 matlab 做出图像,对所 要求的结果分别进行预测,并将结果进行比较。 (二)问题 2 的分析 建立“虚拟种群”相互感染的疾病传播模型, 综合描述人和猩猩疫情的发展, 有助于指导政府通过控制人类与猩猩种群间接触来达到控 制疫情的目的, 可以 帮助政府分析控制两种群接

触所花的成本和疫情控制达到的效果之间的关联程 度。 题中给出了“虚拟人类种群”的有关数据(截至本周末处于发病未隔离状态 的数量、截至本周末总共累计治愈个数、截至本周末总共累计死亡人数、截至本 周末正在被隔离治疗的人数) ,要求我们预测接下来疫情在这两个群体中的发展 情况,并给出 “虚拟人类种群”在第 80 周、第 120 周、第 200 周的相关数据。 由以上原因, 我们可以建立一个考虑因素较问题 1 更多的 SLIRD 模型,利用 matlab 做出图像,对所要求的结果分别进行预测,并将结果进行比较。 (三)问题 3 的分析 在外界专家介入控制人类与猩猩接触的情况下以及某种特效药将被隔离治 疗人群的治愈率提高的条件下。 题中给出了明确的相关条件,即“假设在第 41 周,外界的专家开始介入, 并立即严格控制了人类与猩猩的接触, 且通过某种特效药物将隔离治疗人群的治 愈率提高到了 80%”,即相关参数的值直接由题设给建立病毒的传播模型,有助 于政府分析各种疫情控制措施的严格执行和药物效果的提高等措施对控制疫情 的作用,有助于政府合理规划投入医学研究成本,有重要指导意义 定,其余数据也可使用问题 1、问题 2 中已有数据。要求我们预测接下来疫 情在“虚拟人类种群”的发展情况,对比第 2 问的预测结果说明其作用和影响, 给出“虚拟人类种群”在第 45 周、第 50 周、第 55 周的相关数据。 由问题 3 的以上特点,我们可以直接利用问题 1、问题 2 中所建立的模型, 将相关参数设定成符合问题 3 条件的值,利用 matlab 进行做出图像,对所要求 的结果分别进行预测,并将结果进行比较。

(www.wenku1.com)(四)问题 4 的分析 我们可以依据前述数学模型,分析各种疫情控制措施的严格执行和药物(包 括防疫药物、检疫药物和治疗药物等)效果的提高等措施对控制疫情的作用,写 出分析结果。

三、 模型假设

(1)假设题目所给的数据真实可靠; (2)假设模型体系中没有猩猩和人类的迁入迁出,以及出生和其它原因的死亡; (3)鉴于被隔离的人数一直处于较少且比较稳定的状态,我们假设被隔离的病人 在本周即会被治愈或不幸死亡,不会长久地占用医疗设施; (4)鉴于被隔离的人数与所有发病人数(包括隔离和未隔离的)比值稳定在 0.42 左右,所以我们有理由相信,由于各种实际情况,该地区只能够将发病人群的百 分之四十二隔离治疗,故我们假设,隔离人数=0.72*发病且未隔离人数; (5)由题目可知,发病后的人与猩猩之间的接触不考虑,故假设猩猩的发病不受 人类影响。

四、定义与符号说明

(1)S1—截止到本周末健康猩猩占总猩

猩数的比例。 (2)S2—截止到本周末健康人占总人数的比例。 (3)L1—截止到本周末潜伏者占总猩猩数的比例。 (4)L2—截止到本周末潜伏者占总人数的比例。 (5)I1—截止到本周末处于发病状态的猩猩占总猩猩数的比例。 (6)I2—截止到本周末处于发病状态且未被隔离的人占总人数的比例。 (7)M—截止到本周末处于被隔离状态的人占总人数的比例,M=0.72*I2。

(www.wenku1.com)(8)R1—截止到本周末累计自愈的猩猩占总猩猩数的比例。 (9)R2—截止到本周末累计治愈的人占总人数的比例。 (10)D1—截止到本周末累计死亡的猩猩占总猩猩数的比例。 (11)D2—截止到本周末累计死亡的人占总人数的比例。 (12) ε 1—本周新增发病的猩猩数与截止本周潜伏的猩猩数的比值, 经过计算, ε 1 取平均值 1.05。 (13)ε 2—本周新增发病的人数与截止本周潜伏的人数的比值,经过计算,ε 2 取平均值 0.99。 (14)μ 1—本周自愈的猩猩数与截止本周处于发病的猩猩数的比值,经过计算 μ 1 取平均值 0.1。 (15) Ω —本周新增死亡的猩猩数与截止本周处于发病的猩猩数的比值,经过计 算,取平均值 0.2。 (16) μ 2—人的自愈率, 即自愈的总人数占发病且未被隔离的人的总数的比例, 经过计算取平均值 0.215。 (17)λ 1—每个发病猩猩每周有效接触(足以致病)的猩猩的数目,经过查阅 相关资料,并由题目所给数据测试得,当 λ 1 取 0.28 时,最为符合,故 λ 1 取 0.28。 (18)λ 2—每个病人每周有效接触(足以致病)的人的数目,经过查阅相关资 料, 并由题目所给数据测试得, 当 λ 2 取 0.7228 时, 最为符合, 故 λ 1 取 0.7228。 (19)λ 3—每个发病猩猩每周有效接触(足以致病)的人的数目,经过查阅相 关资料,并由题目所给数据测试得,当 λ 1 取 0.001 时,最为符合,故 λ 1 取 0.001。 (20)β 1—第二问中被隔离的病人的治愈率,即总共治愈的人数与总共隔离治 疗的人数。计算得 β 1=0.59。 (21)β 2—第三问中被隔离的病人的治愈率,由条件知 β 2=0.8。 (22)N1—“虚拟猩猩种群”总数,即 N1=3000。 (23)N2—“虚拟人类种群”总数,即 N2=200000。

(www.wenku1.com)五、模型的建立与求解 第一部分:问题 1 的 SLIRD 模型

(1) 猩猩总数 3000 不变,猩猩种群分为健康猩猩、处于潜伏期猩猩、发病状 态猩猩、已自愈猩猩和已死亡猩猩,称 SLIRD 模型。时刻 t 五类猩猩在猩 猩总数 N 中占的比例分别记作 S1(t) 、L1(t) 、I1(t) 、R1(t)和 D1 (t) 。病猩猩的周接触率为 λ 1,周治愈率为 μ 1,周死亡率为 Ω ,处于 潜伏期猩猩发病比例为 ε 1。 (2) 由模型假设显然有: S1(t)+L1(t

)+I1(t)+R1(t)+D1(t)=1(1.01) 对于健康猩猩、处于潜伏期猩猩、处于发病状态猩猩、已自愈猩猩和已 死亡猩猩有: N*(dS1/dt)=-λ 1*S1*I1*N(1.02) N*(dI1/dt)=N*(ε *L1-μ 1*I1-Ω 1*I1) (1.03) N*(dR1/dt)=N*μ 1*I1(1.04) N*(dD1/dt)=N*Ω *I1(1.05) 则由(1.01) 、 (1.02) 、 (1.03) 、 (1.04) 、 (1.05)式,SLIRD 模型的方 程可以写作: dS1/dt=-λ 1*S1*I1 dL1/dt=λ 1*S1*I1-ε *L1 dI1/dt=ε *L1-μ 1*I1-Ω *I1 dR1/dt=μ 1*I1 dD1/dt=Ω *I1 上述方程无法求出 L1(t) 、I1(t) 、R1(t) 、D1(t)的解析解,我们对其进行

(www.wenku1.com)matlab 数值计算(下图为根据上述模型做出的前 40 周“虚拟猩猩种群”各类状 态猩猩数量变化以及部分放大图) ,并与附件中的数据进行比较,对误差进行数 据分析:

前 40 周“虚拟猩猩种群”各类状态猩猩数量变化图

前 40 周“虚拟猩猩种群”各类状态猩猩数量变化部分放大图 前四周数据与附件所给数据拟合得比较吻合,根据上述模型,利用 matlab 进行数值计算,将时间轴延长(下图为根据上述模型做出的前 200

(www.wenku1.com)周 “虚拟猩猩种群” 各类状态猩猩数量变化以及预测结果) 我们能够得到“虚 拟猩猩种群”在第 80 周、第 120 周、第 200 周的相关数据:

前 200 周“虚拟猩猩种群”各类状态猩猩数量变化图 预测结果: 潜伏群体 80 周 120 周 200 周 0 0 0 处于发病 2 0 0 累计自愈 212 215 216 累计死亡 424 431 431

第二部分:问题 2 的 SLIMRD 模型

(1)猩猩的各种条件不变,人的总数 200000 不变,人群分为健康人、处于潜伏 期的人、发病且被隔离的的人、发病且未被隔离的人,隔离且被治愈的人,死亡 的人,称 SLIMRD 模型。时刻 t 六类人在总数 N2 中占的比例分别记作 S2(t) 、 L2(t) 、I1(t) 、M2(t),R1(t)和 D1(t) 。

(www.wenku1.com)(2)与问题一类似,对于健康人、处于潜伏期的人、处于发病状态且未隔离的 人、隔离且治愈的人和累计死亡的人有:

dS2/dt=-λ 3*S2*I1-I2*λ 2*S2 dL2/dt=λ 3*S2*I1+I2*λ 2*S2-ε 2*L2 dI2/dt=ε 2*L2-M dR2/dt=Mβ 1 dD2/dt=I2(1-μ 2)+M(1-β 1) 此外有 M=0.72*I2

与问题 1 方程组联立: dS1/dt=-λ 1*S1*I1 dL1/dt=λ 1*S1*I1-ε *L1 dI1/dt=ε *L1-μ 1*I1-Ω *I1 dR1/dt=μ 1*I1 dD1/dt=Ω *I1 上述两个方程组无法求出 S2(t) 、L2(t) 、I2(t) 、R2(t) 、D2(t)的解析解, 我们对其进行 matlab 数值计算 (下图为根据上述模型做出的前 40 周 “虚拟人群” 各类状态人的数量变化以及部分放大图) ,并与附件中的数据进行比较,对误差 进行数据分析:

(www.wenku1.com)前 40 周“虚拟人群”各类状态人数量变化图

前 40 周“虚拟人群”各类状态人数量变化部分放大图 前四周数据与附件所给数据拟合得比较吻合,根据上述

模型,利用 matlab 进行数值计算,将时间轴延长(下图为根据上述模型做出的前 200 周“虚拟人群”各类状态人数量变化以及预测结果)我们能够得到“虚拟人 群”在第 80 周、第 120 周、第 200 周的相关数据:

(www.wenku1.com)前 200 周“虚拟人群”各类状态人数量变化图 预测结果 潜伏人群 第 80 周 第 120 周 第 200 周 62 30 12 处于发病状态 46 46 16 隔离治疗 34 34 12 累计治愈 1850 2782 3768 累计因病死亡 4716 7086 9611

第三部分:问题 2 的 SLIMRD 模型

(1) 在问题 3 的情境中,虚拟猩猩种群各类猩猩数量比例的变化同问题 1、问 题 2; (2) 在问题 3 的情境中,前 40 周虚拟人类种群病情发展情况同问题 2 中虚拟 人类种群模型。从第 41 周开始,被隔离人群的治愈率由 β 1 变为了 β 2 (β 2=0.8)

(www.wenku1.com)(3) 故因此我们有从第 41 周开始的两种群常系数微分方程组: 虚拟猩猩种群: dS1/dt=-λ 1*S1*I1 dL1/dt=λ 1*S1*I1-ε *L1 dI1/dt=ε *L1-μ 1*I1-Ω *I1 dR1/dt=μ 1*I1 dD1/dt=Ω *I1 虚拟人类种群: dS2/dt=-I2*λ 2*S2 dL2/dt=I2*λ 2*S2-ε 2*L2 dI2/dt=ε 2*L2-M dR2/dt=Mβ 2 dD2/dt=I2(1-μ 2)+M(1-β 2) 上述两个方程组无法求出 S2(t) 、L2(t) 、I2(t) 、R2(t) 、D2(t)的解析解, 我们对其进行 matlab 数值计算(下图为根据上述模型做出的第 40 到 55 周“虚 拟人群”各类状态人的数量变化以及部分放大图,并附有第 40 到 200 周各类状 态人的数量变化图) ,与附件中的数据进行比较,对误差进行数据分析:

(www.wenku1.com)第 40 到 55 周“虚拟人群”各类状态人数量变化图

第 40 到 55 周“虚拟人群”各类状态人数量变化部分放大图

(www.wenku1.com)第 40 到 200 周“虚拟人群”各类状态人数量变化部分放大图 预测结果: 潜伏群体 45 周 50 周 55 周 22 19 18 处于发病 36 31 27 累计自愈 722 760 794 累计死亡 1851 1922 1989

(www.wenku1.com)四、 模型评价与推广

1.本文将所给的问题归纳为传染病问题,所建立的传染病模型科学合理 2.文中所列的方程符合实际情况和数学法则,合理的设置了一些变量和参数, 使问题清晰明了。 3.文中列出的微分方程无法求出解析值,故使用 MATLAB 进行数值求解,降低了 运算量,并科学地预测出了 40 周以后猩猩种群和人群里疾病的发展趋势 4.方程中涉及到的大部分参数,如 β ,ε 等,均由所给的数据计算取平均值得 到,结果科学而具有说服力。 5.本文中也存在一些无法避免的缺陷,例如有些参数(人或猩猩的有效接触数 等)无法通过计算得到,只能由不断变换数值,观察曲线与所给数值的吻合 度来估算,这可能会导致在预测过程中产生难免的误差。 6.从预测的数据可以看出,如果在第四十周外界专家介入,并立即严格控制了 人类与猩猩

的接触,且通过特效药物将隔离人群的治愈率提高到了 0.8,将 使得第 45 周, 第 50 周和第 55 周以致第 200 周的发病人数,潜伏期人数明显 且持续降低,累计死亡人数也会相应减缓增加。因此,我们有理由相信:控 制人类与猩猩的接触,使用特效药物可以有效遏制埃博拉的发展,同时我们 可以做出合理的推测:加强对人群的监控,改善医疗条件,提高隔离人数占 发病总人数的比例,也可以有效地抵抗埃博拉的传播。

(www.wenku1.com)五、 参考文献

[1]姜启源,数学模型(第二版) ,高等教育出版社,1991 [2]任善强,数学模型,重庆大学出版社,1987 [3]叶其孝主编,大学生数学建模竞赛辅导教材(四) ,湖南教育出版社,2001 [4]张志涌等编著,精通 MATLAB,北京航空航天大学出版社,2011 [5]何万生等,数学模型与建模,甘肃教育出版社,2001 [6]杨学桢,数学建模方法,河北大学出版社,2000 [7]贾晓峰等,微积分与数学模型,高等教育出版社,1999 [8]孙祥编著,清华大学出版社,matlab7.0 基础教程 [9]于润伟编著,matlab 基础及应用 [10]薛定宇,高等应用数学问题求解 matlab 求解 [11]薛定宇,数学物理方程的 matlab 解法与可视化

(www.wenku1.com)六、 附件 第一部分:问题 1 的 matlab 程序代码

t 属于 0 到 40 之间时: function dy=rigid(t,y) dy=zeros(5,1); dy(1)=-0.28*y(3)*y(1); dy(2)=0.28*y(3)*y(1)-1.05*y(2); dy(3)=1.05*y(2)-0.3*y(3); dy(4)=0.1*y(3); dy(5)=0.2*y(3); end

[T,Y]=ode45('rigid',[0 40],[1 1/30 0 0 0]); plot(T,3000*Y(:,1),'-',T,3000*Y(:,2),'*',T,3000*Y(:,3),'+',T,3000*Y(: ,4),'') t 属于 0 到 200 之间时: function dy=rigid(t,y) dy=zeros(5,1); dy(1)=-0.28*y(3)*y(1); dy(2)=0.28*y(3)*y(1)-1.05*y(2); dy(3)=1.05*y(2)-0.3*y(3); dy(4)=0.1*y(3);

(www.wenku1.com)dy(5)=0.2*y(3); end

[T,Y]=ode45('rigid',[0 200],[1 1/30 0 0 0]); plot(T,3000*Y(:,1),'-',T,3000*Y(:,2),'*',T,3000*Y(:,3),'+',T,3000*Y(: ,4),'')

第二部分:问题 2 的 matlab 程序代码

t 属于 0 到 40 之间时: function dy=rigid(t,y) dy=zeros(10,1); dy(1)=-0.28*y(3)*y(1); dy(2)=0.28*y(3)*y(1)-1.05*y(2); dy(3)=1.05*y(2)-0.3*y(3); dy(4)=0.1*y(3); dy(5)=0.2*y(3); dy(6)=-0.001*y(3)*y(6)-0.7228*y(8)*y(6); dy(7)=0.001*y(3)*y(6)+0.7228*y(8)*y(6)-0.99*y(7); dy(8)=0.99*y(7)-0.724*y(8); dy(9)=0.427*y(8); dy(10)=1.089*y(8); end

(www.wenku1.com)[T,Y]=ode45('rigid',[0 40],[1 1/30 0 0 0 1 0 0 0 0]);..... plot(T,200000*Y(:,6),'*',T,200000*Y(:,7),'+',T,200000*Y(:,8),'-',T,20 0000*Y(:,9),'') t 属于 0 到 200 之间时: function dy=rigid(t,y) dy=zeros(10,1); dy(1)=-0.28*y(3)*y(1); dy(2)=0.28*y(3)*y(1)-1.05*y(2); dy(3)=1.05*y(2)-0.3*y(3); dy(4)=0.1*y(3); dy(5)=0.2*y(3); dy(6)=-0.001*y(3)*y(6)-0.7228*y(8)*y(6); dy(7)=0.001*y(3)*y(6)+0.7228*y(8)*y(6)-0.99*y(7); dy(8)=0.99*y(7)-0.724*y(8); dy(9)=0.427*y(8); dy(10)=1.089*y(8); end

[T,Y]=ode45('rigid',[0 200],[1 1

/30 0 0 0 1 0 0 0 0]);..... plot(T,200000*Y(:,6),'*',T,200000*Y(:,7),'+',T,200000*Y(:,8),'-',T,20 0000*Y(:,9),'')

(www.wenku1.com)第三部分:问题 3 的 matlab 程序代码

t 属于 40 到 55 之间时: function dy=rigid(t,y) dy=zeros(5,1); dy(1)=-0.7228*y(3)*y(1); dy(2)=0.7228*y(3)*y(1)-0.99*y(2); dy(3)=0.99*y(2)-0.724*y(3); dy(4)=0.5792*y(3); dy(5)=0.9298*y(3); end [T,Y]=ode45('rigid',[0 15],[0.98772 0.00035 0.000205 0.003415 0.008825]);..... plot(T,200000*Y(:,1),'*',T,200000*Y(:,2),'+',T,200000*Y(:,3),'-',T,20 0000*Y(:,4),'') t 属于 40 到 200 之间时: function dy=rigid(t,y) dy=zeros(5,1); dy(1)=-0.7228*y(3)*y(1); dy(2)=0.7228*y(3)*y(1)-0.99*y(2); dy(3)=0.99*y(2)-0.724*y(3); dy(4)=0.5792*y(3); dy(5)=0.9298*y(3);

(www.wenku1.com)end [T,Y]=ode45('rigid',[0 160],[0.98772 0.00035 0.000205 0.003415 0.008825]);..... plot(T,200000*Y(:,1),'*',T,200000*Y(:,2),'+',T,200000*Y(:,3),'-',T,20 0000*Y(:,4),'')

(www.wenku1.com)

范文九:埃博拉病毒不太可能经由空气传播

埃博拉 病毒不太可 能经 由空气传播

T i n a   He s ma n   S a e y

猪 身上 携 带埃 博 拉 病 毒 。在 菲律 宾 的 猪 身上 曾发现 过 雷 斯  顿 ( R e s t o n ) 埃博拉病毒 , 但这种病毒不会对人引发疾病。   目前所 知 ,还 没 有 猪 受到 扎 伊 尔埃 博 拉 病毒 的感 染 ,而 这

种病毒正导致 了当前的疾病大暴发。   即使猪能够 通过空 气传播病 毒,它 们也有可能是 唯

具 备 此 能 力 的生 物 。科 宾 格和 另一 个 团队 对 此进 行 了进  步 研 究 ,他 们 将 最 新 研 究 成 果 发 表 在 7 月2 5日出 版 的  ( ( S c i e n t i f i c   R e p o r t s ) ) 上 ,他 们 发现 没 有 证 据表 明患 病 的

在一项 实验 中感染埃博 拉病 毒的恒河猴 ( 左 ),不能将病毒 经空气传播 给

猕猴 ( 右 )

短尾 猿可 将病 毒通 过空 气粒 子传 播给 健康 的猴 子 。   研 究人 员将装 有 2只受N? L 伊尔埃博拉 病毒 感染 的恒 河  猴 的笼子放 在未 受感 染的猕 猴旁 。猴子 间不能接 触 ,但 也没  有特 殊 的屏 蔽对 未受 感染 的猴 子作 出保护 。6 天后 ,受 到感

染 的猴 子 死亡 。而 那 些猕 猴 则 保持 未 被 传染 的状 态 达 到 2 8   天之久 ,这个时 间远长 于埃博拉病 毒 6~ 1 6天 的传染期 。

研 究 人 员报 告 说 ,猴 子 相 互 间 不能 通 过 空 气 传播 埃 博

拉 病 毒 。这 个结 果确 定 了对 人 类疾 病 暴 发 的 观 测结 论 ,即  感 染 这种 致 命的 病毒 需要 进 行体 液接 触 。   这项 研 究是 对 2 0 1 2年一 项 研究 的继 续 。在那项 研 究 中 ,   以 来 自干 加 拿 大 公 共 卫 生机 构 的 感 染 病学 专 家 加 里 ・ 科 宾  格 ( G a r y   K o b i n g e r )为 首 的研 究 团 队提 出埃 博 拉 病 毒 有  可 能 会 经 由空 气 传 播 。他 们 将 放 置短 尾 猿 的 笼 子 放 在一 个  猪圈 旁 ,猪 圈 中的 猪 均 感 染 了埃 博 拉 病毒 ,不 久 他 们发 现  那 些 短 尾 猿 也 受 到感 染 。这 期 间短 尾 猿 与 猪从 未 接 触 过 。   由此 ,研 究 人 员认 为 ,他们 的 发 现 意 味着 病 毒 会 像 完 美 的  空 气 粒 子 飞沫 飘 到 短 尾 猿 的 笼子 中 ,而那 些 粒 子 飞 沫 正是  由猪释 放 出来 的 。

猪 看起 来 会 比其 他 生 物 释 放 出更 多的 雾 化病 毒 颗 粒 ,

自从 2 0 1 3 年 l 2月以来 ,这次有史以来最大规模的埃  博 拉 病毒 暴 发 已经在 几 内亚 、利 比 里亚 、塞拉 利 昂 和尼 日

利 亚 至 少 夺去 了 9 3 2 人 的 生

命 ,死 亡率 高于 5 0 % 。 在那 些

被 感 染 的 人 中 ,包括 2名 来 自于 美 国卫 生 保 健机 构 的 护理

人 员 ,他们 均 对 受到 致 命 病毒 感 染 的患 者 进 行过 护 理 。 目   前这 2名美 国患者肯特与南希已经被转送到美国亚特兰大

E mo r y大 学 医院进 行 隔离治 疗 。

来 自 于 柏 林 罗 伯 特 ・科 赫 研 究 院 ( R o b e r t   K o c h

I n s t i t u t e ) 的 流 行 病 学 及 疾 病 生 态 学 家 法 比 安 ・林 德 兹  ( F a b i a n   L e e n d e r t z )说 道 :疾 病正 在通 过 人 与 人 的 接 触

来 自于 英 国 兰 开斯 特 大 学 的病 毒 进 化 学 家 德 里克 ・ 盖泽 尔  ( D e r e k   G a t h e r e r )说道 。 “ 如果 病毒 真 能经 浮 质传播 的 话 ,

那 么 猪绝 对就 是 帮 凶了 。”他表 示 道 。

进 行 传 播 。其 中女性 的死 亡 率 最 高 ,因 为她 们需 要 照顾 患  者、 她 们的 孩子 , 以及那 些 参加过 葬 礼而触 碰到 尸体 的人 。

卫 生保 健 人 员也是 高危 人群 。

但 是 他 并 不 认 为研 究 人 员 已最 终 确 定埃 博拉 病 毒 能 够  经 由空 气 传 播 。 含有 病 毒 的水 滴 也 有 可 能 在研 究 人 员清 洗  猪圈 时溅 到短 尾 猿的 笼子 中 ,他 说道 。   猪有 可 能 并 非是 埃 博 拉 病 毒 暴发 的 源 头 。 只有 极 少 数

但 是埃 博拉 病 毒也并 非 像许 多 人想象 的那 样 易于 传播 ,   林 德兹 表 示 道 。即 使一 名感 染者 乘坐 飞 机 飞往 美 国 、欧 洲  或 者 其 他 地 方 ,密 切 的卫 生 保健 措 施 也 将确 保 埃博 拉 病 毒  “ 不 会走 得很 远 ” ,林 德兹 说道 。 ■

王麒 译自 科学新闻 ( S c i e n c e Ne ws) 2 01 4年 9月 6日

本 栏 目 文字 与 图 片 由

( 美 国科 学新 闻 )版 权 所 有 ,授 权 我 刊 使 用 , 不得 以任 何 形 式转 载 、摘 编 ,特 此 声 明。

2 0 1 4 . 1 0 中国科技教 育  6 7

范文十:埃博拉病毒的传播分析

摘 要 由于猩猩的潜伏期无法忽略,因此我们建立了经典的传染病传播SEIR模型。但由于在发病的不同时期患病者的感染能力不同,因此我们在此模型的基础上对感染人数增加了一个扰动项。我们未考虑猩猩的出生率和死亡率,从而得到适合此问题的改进的SEIR模型。通过对模型中参数的求解,我们可以对埃博拉病毒在猩猩群体的传播进行预测。   关键词 埃博拉病毒 改进的SEIR模型 最小二乘法   1问题分析   1.1 问题一的分析   问题一是对埃博拉病毒在猩猩种群中传播感染规律的研究。此问题属于传染病传播问题,解决此类问题一般建立经典的SIR仓室模型进行求解。考虑到本问题中埃博拉病毒的潜伏期对本问题的结果有不可忽略的影响,因此我们在经典SIR模型的基础上增添潜伏期这一仓室建立SEIR模型。考虑到猩猩种群的总数较小(3000只),并且猩猩的繁殖周期�^长(3-4年),因此在本模型中忽略猩猩自然出生率及死亡率对问题结果的影响。   2模型假设   (1)假设猩猩的自然出生率及死亡率对预测结果无影响;   (2)假设猩猩和人都处在除埃博拉疫情外正常的生活状态,不受其它因素影响;   (3)假设感染者在自愈的同时获得免疫力,不会再感染病毒;   3模型的建立   3.1 问题一模型的建立   问题一是考虑在猩猩种群中埃博拉病毒的传播情况,本问题中埃博拉病毒的潜伏期对本问题的结果有不可忽略的影响。猩猩种群的总个体数较小(3000只),并且猩猩的繁殖周期较长(3-4年),因此在本模型中忽略猩猩自然出生率及死亡率对问题结果的影响。同时考虑到有效传染率�%[、潜伏期到患病的转换率�%e、死亡率�%]、治愈率r等因素的影响。   参考文献   [1] 罗建军,杨琪.MATLAB教程[M].北京:电子工业出版社,2005.   [2] 张发.传染病传播模型综述[J].系统工程理论与实践,2011,9,31(9).   [3] 赵一楠.一类具有无穷时滞的SEIR传染病模型及其边界问题[D].扬州大学,2012.