1、甲型 H1N1 流感传播模型研究摘要本文采用了 SIR 模型对的甲型 h1n1 流感病毒的传播规律进行了研究和预测，文章收集了美国地区的甲流实验室确认病例数量的数据，对模型进行了验证，并提出了如何降低流感在人群中发病率的俩种可靠方法。一、问题重述近年来由墨西哥发端的甲型 h1n1 型流感（又称猪流感）正成为人们关注的焦点，通过相关网站获得数据，建立一个模型对甲型 h1n1 流感的走势进行预测。二、问题分析甲型 h1n1 流感的传播是一道传染病问题。在数学建模领域已经有很多关于这方面的研究，其中 SIR 模型是比较完整的模型。SIR 模型通过建立微分方程组，按照一般的传播机理建立集中模型。本文选
2、取美国地区的甲流实验室确认病例数量，建立 SIR 模型，对甲型 h1n1 流感的传播规律进行预测。美国甲型 H1N1 流感实验室确认病例数量：时间 确诊（包括死亡病例） 死亡（累计）4 月 23 日 5 04 月 24 日 8 04 月 25 日 11 04 月 26 日 20 04 月 27 日 40 04 月 28 日 64 04 月 29 日 91 04 月 30 日 109 15 月 1 日 141 15 月 2 日 160 15 月 3 日 226 15 月 4 日 279 15 月 5 日 403 15 月 6 日 642 25 月 7 日 896 25 月 8 日 1639 25
3、 月 9 日 2254 25 月 10 日 2532 35 月 11 日 2600 35 月 12 日 3009 35 月 13 日 3352 45 月 14 日 4298 45 月 15 日 4714 4三、建立模型（一） 、不考虑潜伏期的数学模型1、模型假设（1） 、在甲型 H1N1 流感传播期内，美国境内的总人数为 N 亿不变，既不考虑生死，也不考虑迁移，人群分为易感染者 S，发病人群 I 和退出人群 R(括死亡者和治愈者)四类，时刻 t 内这三类人在总人数中所占比例分别为s(t)、i(t)、r(t)。（2） 、i(t)关于时间的增长率与 s(t)成正比，比例常数为 。病人的数量减少速度
4、与当时的病人总人数成正比，比例常数为 。治愈的病人具有了免疫力，即治愈后不再会成为二次患者。（3） 、s(t)、r(t)、i(t)之和是一个常数 1。2、模型构成易感者和发病者有效接触后成为发病者者。设每个发病者平均每天有效接触的易感者数为 ， 个发病者平均每天能使 个易感者成为()StNIt ()StNI病毒潜伏者。所以有：（1）()()dttI单位时间内退出者的变化等于发病人群的减少，即（2）()dRtIt发病人群的变化等于易感人群转入的数量，即（3）()()dItStIt记初始时刻的健康者和病人的比例分别为 、 （不妨设 =0） 。0SR03、模型求解方程组（1） 、 （2） 、 （3）
5、无法求出解析解，我们定义一个新的变量 ，/于是可以求出方程的解为：（4）001()lnsisi下面分析 s(t)、i(t)、r(t)的变化情况：a、不论初始条件 、 如何，病人最终将消失，即 。0SR0ib、最终未被感染者的健康者的比例是 ，是方程s在 内的根。001()lnssi(,1/)C、若 ，则开始有： 先增加。当 时， 达到最大值，然后0/s()it01/s()it减小且趋于零， 则单调减小至 。()itstd、若 ，则 单调减小至 5， 则单调减小至 。01/s()it()sts我们发现人们的卫生水平越高，日接触率越小；医疗水平越高，日治愈率越高，于是 越小，所以提高卫生水平和医疗
6、水平有利于传染病的蔓延。结合美国的具体情况和假设条件进行分析：根据所得的数据画出美国患病人数变化曲线和治愈人数变化曲线：根据图形来看，甲型 h1n1 流感在美国呈现出蔓延的形式，即现在属于 01/s的情况，即 。由假设条件可知 的取值范围在 之间。现0/1s1.46在我们取 =1.6，则表示 ，即美国每天平均治愈的人数最多0/()1.6为 1.6 人，这与美国疾病预防与控制中心所发布的数据不同。如果美国平均每天治愈 1.6 个人的话，那么从 4 月 23 日期，治愈的总人数为 人，1.6*23.8这与实际的情况相差甚远。产生这个问题的原因有以下几个方面：第一：对每个病人每天有效接触的平均人数估
7、计值偏小。不是简单的成正比关系，应该是成多次方关系，甚至是指数关系。第二：美国疾病预防与控制中心所得到的数据具有滞后性。第三：在美国 不一定成立。可以把那些身体强壮的、注意自己个人卫生的0s人排除在外。（二） 、考虑潜伏期的数学模型1、模型假设（1） 、在甲型 H1N1 流感传播期内，美国境内的总人数为 N 亿不变，既不考虑生死，也不考虑迁移，人群分为易感染者 S，病毒潜伏人群 E,发病人群 I 和退出人群 R(括死亡者和治愈者)四类，时刻 t 内这三类人在总人数中所占比例分别为 。steitr、 、 、（2） 、每个病人每天有效接触的平均人数为 ，称为日接触率，当已感染者与易感染者有效接触时
8、，使易感染者变为病毒潜伏人群，病毒潜伏人群过一段时间再转换成发病人群，发病人群被治愈。2、模型构成易感者和发病者有效接触后成为病毒潜伏者。设每个发病者平均每天有效接触的易感者数为 ， 个发病者平均每天能使 个易感()tS()NIt ()()tSNIt者成为病毒潜伏者。所以有化简得：()()dStNtNIttStt病毒潜伏人群的变化等于易感人群转入数量减去转化为发病人群的数量，即 ()()()dEtStItE其中 表示潜伏期日发病率，即每个潜伏者平均有效发病的人数。()t单位时间内退出者的变化等于发病人群的减少，即 ()()dRtIt其中 表示日退出率，即每个病人平均有效病情结束的人数。()t发
9、病人群的变化等于潜伏人群转入的数量，即 ()()()dItEtIt()()1stirt初始时刻易感染者，已感染者与病愈免疫者的比例分别是 000(),(),sir3、模型求解由于潜伏期的人群数量不能确定，所以可视为是易感人群的一部分，因此求解过程跟忽略潜伏期的一样。四、模型的改进就如何确定日接触率 的值。就如何确定日接触率可以进行改进，根据以前的流感疫情治愈率，加权平均得到值，而不是简单的是一个正比关系。病毒在人群中的传播刚开始阶段一个有一个爆发阶段，该阶段的日接触率 很大，可设为是一个冲激变量。参考文献：1姜启源 谢金星 叶俊 数学建模（第四版）高等教育出版社2数据来源：美国疾病预防控制中心
