传染病数学建模论文

甲型H1N1流感传播模型研究
摘要
本文采用了SIR模型对的甲型h1n1流感病毒的传播规律进行了研究和预测，文章收集了美国地区的甲流实验室确认病例数量的数据，对模型进行了验证，并提出了如何降低流感在人群中发病率的俩种可靠方法。

一、问题重述
近年来由墨西哥发端的甲型h1n1型流感（又称猪流感）正成为人们关注的焦点，通过相关网站获得数据，建立一个模型对甲型h1n1流感的走势进行预测。

二、问题分析
甲型h1n1流感的传播是一道传染病问题。

在数学建模领域已经有很多关于这方面的研究，其中SIR模型是比较完整的模型。

SIR模型通过建立微分方程组，按照一般的传播机理建立集中模型。

本文选取美国地区的甲流实验室确认病例数量，建立SIR模型，对甲型h1n1流感的传播规律进行预测。

美国甲型H1N1流感实验室确认病例数量：
三、建立模型
（一）、不考虑潜伏期的数学模型 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 dt
λ=-（1） 单位时间内退出者的变化等于发病人群的减少，即
()
()dR t I t dt
ν=（2） 发病人群的变化等于易感人群转入的数量，即
()
()()()dI t S t I t I t dt
λν=-（3）
记初始时刻的健康者和病人的比例分别为0S 、0R （不妨设0R =0）。

3、模型求解
方程组（1）、（2）、（3）无法求出解析解，我们定义一个新的变量/σλν=，于是可以求出方程的解为：
000
1()ln
s
i s i s s σ
=+-+
（4） 下面分析s(t)、i(t)、r(t)的变化情况：
a 、不论初始条件0S 、0R 如何，病人最终将消失，即0i ∞=。

b 、最终未被感染者的健康者的比例是s ∞，是方程
000
1()ln
0s
s i s s σ
+-+
=在(0,1/)σ内的根。

C 、若01/s σ>，则开始有：()i t 先增加。

当01/s σ=时，()i t 达到最大值，然后()i t 减小且趋于零，()s t 则单调减小至s ∞。

d 、若01/s σ≤，则()i t 单调减小至5，()s t 则单调减小至s ∞。

我们发现人们的卫生水平越高，日接触率越小；医疗水平越高，日治愈率越高，于是σ越小，所以提高卫生水平和医疗水平有利于传染病的蔓延。

结合美国的具体情况和假设条件进行分析：
根据所得的数据画出美国患病人数变化曲线和治愈人数变化曲线：
根据图形来看，甲型h1n1流感在美国呈现出蔓延的形式，即现在属于
01/
sσ>
的情况，即
0/1/s σλν=>。

由假设条件可知λ的取值范围在1.4~1.6之间。

现在
我们取λ=1.6，则表示0/(1/) 1.6s νλ<=，即美国每天平均治愈的人数最多为1.6人，这与美国疾病预防与控制中心所发布的数据不同。

如果美国平均每天治愈1.6个人的话，那么从4月23日期，治愈的总人数为1.6*2336.8=人，这与实际的情况相差甚远。

产生这个问题的原因有以下几个方面：
第一：对每个病人每天有效接触的平均人数估计值偏小。

不是简单的成正比关系，应该是成多次方关系，甚至是指数关系。

第二：美国疾病预防与控制中心所得到的数据具有滞后性。

第三：在美国00s ≈不一定成立。

可以把那些身体强壮的、注意自己个人卫生的人排除在外。

（二）、考虑潜伏期的数学模型
1、模型假设
（1）、在甲型H1N1流感传播期内，美国境内的总人数为N 亿不变，既不考虑生死，也不考虑迁移，人群分为易感染者S ，病毒潜伏人群E,发病人群I 和退出人群R(括死亡者和治愈者)四类，时刻t 内这三类人在总人数中所占比例分别为
()()()()s t e t i t r t 、、、。

（2）、每个病人每天有效接触的平均人数为λ，称为日接触率，当已感染者与易感染者有效接触时，使易感染者变为病毒潜伏人群，病毒潜伏人群过一段时间再转换成发病人群，发病人群被治愈。

2、模型构成
易感者和发病者有效接触后成为病毒潜伏者。

设每个发病者平均每天有效接触的易感者数为()()t S t λ，()NI t 个发病者平均每天能使()()()t S t NI t λ个易感者成为病毒潜伏者。

所以有()
()()()dS t N
t S t NI t dt λ=-化简得： ()
()()()dS t t S t I t dt
λ=- 病毒潜伏人群的变化等于易感人群转入数量减去转化为发病人群的数量，即
()
()()()()()dE t t S t I t t E t dt
λα=- 其中()t α表示潜伏期日发病率，即每个潜伏者平均有效发病的人数。

单位时间内退出者的变化等于发病人群的减少，即
()
()()dR t t I t dt
ν= 其中()t ν表示日退出率，即每个病人平均有效病情结束的人数。

发病人群的变化等于潜伏人群转入的数量，即
()
()()()()dI t t E t t I t dt
αν=- ()()()1s t i t r t ++=
初始时刻易感染者，已感染者与病愈免疫者的比例分别是
00000(0),(0),0
s s i i r >>=
3、模型求解
由于潜伏期的人群数量不能确定，所以可视为是易感人群的一部分，因此求解过程跟忽略潜伏期的一样。

四、模型的改进
就如何确定日接触率λ的值。

就如何确定日接触率可以进行改进，根据以前的流感疫情治愈率，加权平均得到值，而不是简单的是一个正比关系。

病毒在人群中的传播刚开始阶段一个有一个爆发阶段，该阶段的日接触率λ很大，可设为是一个冲激变量。

参考文献：
[1]姜启源谢金星叶俊数学建模（第四版）高等教育出版社
[2]数据来源：美国疾病预防控制中心。