甲型 H1N1 流感传播模型研究.docxVIP

----宋停云与您分享---- ----宋停云与您分享---- 甲型 H1N1 流感传播模型研究 摘要 本文采用了 SIR 模型对的甲型 h1n1 流感病毒的传播规律进行了 研究和预测，文章收集了地区的甲流实验室确认病例数量的数据， 对模型进行了验证，并提出了如何降低流感在人群中发病率的俩种可 靠方法。 一、问题重述 近年来由墨西哥发端的甲型 h1n1 型流感（又称猪流感）正成为人们关注的焦点，通过相关网站获得数据，建立一个模型对甲型 h1n1 流感的走势进行预测。 二、问题分析 甲型 h1n1 流感的传播是一道传染病问题。在数学建模领域已经有很多关于这方面的研究，其中 SIR 模型是比较完整的模型。SIR 模型通过建立微分方程组， 按照一般的传播机理建立集中模型。本文选取地区的甲流实验室确认病例数量，建立 SIR 模型，对甲型 h1n1 流感的传播规律进行预测。 甲型 H1N1 流感实验室确认病例数量： 时间 确诊（包括死亡病例） 死亡（累计） 4 月 23 日 5 0 4 月 24 日 8 0 4 月 25 日 11 0 4 月 26 日 20 0 4 月 27 日 40 0 4 月 28 日 64 0 4 月 29 日 91 0 4 月 30 日 109 1 5 月 1 日 141 1 5 月 2 日 160 1 5 月 3 日 226 1 5 月 4 日 279 1 5 月 5 日 403 1 ----宋停云与您分享---- ----宋停云与您分享---- 5 月 6 日 642 2 5 月 7 日 896 2 5 月 8 日 1639 2 5 月 9 日 2254 2 5 月 10 日 2532 3 5 月 11 日 2600 3 5 月 12 日 3009 3 5 月 13 日 3352 4 5 月 14 日 4298 4 5 月 15 日 4714 4 三、建立模型 （一）、不考虑潜伏期的数学模型 1、模型假设 （1）、在甲型 H1N1 流感传播期内，境内的总人数为 N 亿不变，既不考虑生死，也不考虑迁移，人群分为易感染者 S，发病人群 I 和退出人群 R(括死亡者和治愈者)四类，时刻 t 内这三类人在总人数中所占比例分别为 s(t)、 i(t)、r(t)。 （2）、i(t)关于时间的增长率与 s(t)成正比，比例常数为 ? 。 病人的数量减少速度与当时的病人总人数成正比，比例常数为? 。治愈的病人具有了免疫力，即治愈后不再会成为二次患者。 （3）、s(t)、r(t)、i(t)之和是一个常数 1。 2、模型构成 易感者和发病者有效接触后成为发病者者。设每个发病者平均每天有效接触的易感者数为? S (t) ，NI (t) 个发病者平均每天能使? S (t)NI (t) 个易感者成为病毒 潜伏者。所以有： dS (t) ? ?? S (t)I (t) （1） dt 单位时间内退出者的变化等于发病人群的减少，即 dR(t) ?? I (t) （2） dt 发病人群的变化等于易感人群转入的数量，即 dI (t) ? ?S (t)I (t) ?? I (t) （3） dt ----宋停云与您分享---- ----宋停云与您分享---- 记初始时刻的健康者和病人的比例分别为 S 、 R 0 0 3、模型求解 （不妨设 R 0 =0）。 ----宋停云与您分享---- ----宋停云与您分享---- 方程组（1）、（2）、（3）无法求出解析解，我们定义一个新的变量 ? ? ? /? ，于是可以求出方程的解为： ----宋停云与您分享---- ----宋停云与您分享---- i ? (s 0 ? i ) ? s ? 1 0 ? s ln s （4） 0 ----宋停云与您分享---- ----宋停云与您分享---- 下面分析 s(t)、i(t)、r(t)的变化情况： ----宋停云与您分享---- ----宋停云与您分享---- a、不论初始条件 S 、 R 0 0 如何，病人最终将消失，即i ? ? 0 。 ----宋停云与您分享---- ----宋停云与您分享---- ----宋停云与您分享---- ----宋停云与您分享---- b、最终未被感染者的健康者的比例是 s ? ，是方程 ----宋停云与您分享---- ----宋停云与您分享---- (s ? i ) ? s ? 1 ln s ? 0 在(0,1/ ? ) 内的根。 0 0 ? s 0 ----宋停云与您分享---- ----宋停云与您分享---- C、若 s 0 ? 1/ ? ，则开始有：i(t) 先增加。当 s 0 ? 1/ ? 时，i(t) 达到最大值，然后i(t) ----宋停云与您分享---- ----宋停云与您分享---- 减小