传染病模型数学建模论文(2)

1、甲型H1N1流感传播模型研究摘要本文采用了 SIR模型对的甲型hlnl流感病毒的传播规律进行了研究和预测，文章收集了美国地区的甲流实验室确认病例数量的数据，对模型进行了验证，并提出了如何降低流感在人群中发病率的俩种可靠方法。一、问题重述近年来由墨西哥发端的甲型hlnl型流感（又称猪流感）正成为人们关注的 焦点，通过相关网站获得数据，建立一个模型对甲型hlnl流感的走势进行预测。二、问题分析甲型hlnl流感的传播是一道传染病问题。在数学建模领域已经有很多关于 这方面的研究，其中SIR模型是比较完整的模型。SIR模型通过建立微分方程组, 按照一般的传播机理建立集中模型。本文选取美国地区的甲流实验室

2、确认病例数 量，建立SIR模型，对甲型hlnl流感的传播规律进行预测。美国甲型H1N1流感实验室确认病例数量:时间确诊（包括死亡病例）死亡（累计）4月23日504月24日804月25日1104月26日20104月27日4004月28日6404月29日91104月30日10915月1日14115月2日160r 15月3日22615月4日27915月5日403r 15月6日64225月7日89625月8日1639回5月9日225425月10日253235月11日2600回5月12日300935月13日335245月14日429845月15日47144三、建立模型(一)、不考虑潜伏期的数学模型1、模

3、型假设(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)个发病者平均每

4、天能使S(t)NI (t)个易感者成为病毒潜伏者。所以有：(1)(2)(3)等 S(t)I(t)dt单位时间内退出者的变化等于发病人群的减少，即飒 I(t)dt发病人群的变化等于易感人群转入的数量，即个 S(t)I(t) I(t)dt记初始时刻的健康者和病人的比例分别为 So、R0 (不妨设R=0)3、模型求解方程组(1)、(2)、(3)无法求出解析解，我们定义一个新的变量/于是可以求出方程的解为：1 s/ ,、i (So io) s ln (4)So下面分析s、i(t) 、r(t)的变化情况：a、不论初始条件&、Ro如何，病人最终将消失，即i 0。b、最终未被感染者的健康者的比例是s ,是方

5、程1 s(So io) s -ln 0 在(0,1/ )内的根。 soG若So 1/ ，则开始有：i先增加。当So 1/时，i(t)达到最大值，然后i(t)减小且趋于零，s(t)则单调减小至s od、若& 1/ ,则i(t)单调减小至5, s(t)则单调减小至s 。我们发现人们的卫生水平越高，日接触率越小；医疗水平越高，日治愈率越高，于是越小，所以提高卫生水平和医疗水平有利于传染病的蔓延。结合美国的具体情况和假设条件进行分析：根据所得的数据画出美国患病人数变化曲线和治愈人数变化曲线：根据图形来看，甲型h1n1流感在美国呈现出蔓延的形式，即现在属于 与1/的情况，即/1/soo由假设条件可知的取

6、值范围在1.41.6之间。现在我们取 =1.6,则表示 /(1/so) 1.6,即美国每天平均治愈的人数最多为1.6人，这与美国疾病预防与控制中心所发布的数据不同。如果美国平均每天治愈1.6个人的话，那么从4月23日期，治愈的总人数为1.6*23 36.8人，这与 实际的情况相差甚远。产生这个问题的原因有以下几个方面：第一：对每个病人每天有效接触的平均人数估计值偏小。 不是简单的成正比关系， 应该是成多次方关系，甚至是指数关系。第二：美国疾病预防与控制中心所得到的数据具有滞后性。第三：在美国so o不一定成立。可以把那些身体强壮的、注意自己个人卫生的人排除在外。(二)、考虑潜伏期的数学模型1、

7、模型假设(1)、在甲型H1N1流感传播期内，美国境内的总人数为 N亿不变，既不考虑生 死，也不考虑迁移，人群分为易感染者S,病毒潜伏人群E,发病人群I和退出人 群R(括死亡者和治愈者)四类，时刻t内这三类人在总人数中所占比例分别为st、et、it、rt。(2)、每个病人每天有效接触的平均人数为 ，称为日接触率，当已感染者与易 感染者有效接触时，使易感染者变为病毒潜伏人群，病毒潜伏人群过一段时间再 转换成发病人群，发病人群被治愈。2、模型构成易感者和发病者有效接触后成为病毒潜伏者。设每个发病者平均每天有效接 触的易感者数为(t)S(t), NI (t)个发病者平均每天能使 S(t)NI个易感者成 为病毒潜伏者。所以有N吧D(t)S(t)NI(t)化简得：dt病毒潜伏人群的变化等于易感人群转入数量减去转化为发病人群的数量，即 其中(t)表示潜伏期日发病率，即每个潜伏者平均有效发病的人数。单位时间内退出者的变化等于发病人群的减少，即其中(t)表示日退出率，即每个病人平均有效病情结束的人数。发病人群的变化等于潜伏人群转入的数量，即初始时刻易感染者，已感染者与病愈免疫者的比例分别是3、模型求解由于潜伏期的人群数量不能确定，所以可视为是易感人群的一部分，因此求 解过程跟忽略潜伏期的一样。四、模型的改进就如何确定日接触率 的值。就如何确定日接触率可以进行改进， 根据以前 的流感疫情治