基于媒体报道下的一类SIRS传染病模型研究

张 林,李存林,郭文娟

(1.北方民族大学数学与信息科学学院,宁夏银川 750021)

(2.北方民族大学管理学院,宁夏银川 750021)

1 引言

随着社会经济的快速发展,传染病和流行病给人们的生活、社会发展等带来了很大的危害.越来越引起了许多生物数学和疾病预防等工作者的重视,并取得了一些研究成果[1–5].考虑如下经典的SIRS传染病模型[6]:

其中符号皆为正的常数,意义表示如表1.

表1:符号及其含义

众所周知,通过媒体报道可以降低人与人之间的接触率,这在2003年的非典[4,7,8]和2013年的H7N9[9]传染病中已经得到了证实.作为一种新型的传染病H7N9于2013年首次出现在上海,它和非典一样,短时间内在人类迅速传播开.后来人们通过电视的报道了解到H7N9是通过人与人之间的接触传播的,于是大部分人采取尽量少出门,少参加一些社会活动来降低人与人之间的接触率,这确实在一定程度上减少了疾病的传播.

然而,模型(1.1)中并没有考虑媒体报道对传染病的影响.事实上,由于媒体报道的因素,传染率β是会减小的[4].因此将传染率β表示为媒体报道的函数,即β=β1−β2f(I),那么模型(1.1)转换为

其中β1是不考虑感染者的一般接触率,β2是因感染者的存在而减少的最大接触率.因为每个人与他人的接触不可避免,故而假设β1＞β2.函数f(I)满足

模型(1.2)是在模型(1.1)的基础上考虑了媒体报道对疾病的影响,并且证明了R0=1的情形,即当R0=1时,无病平衡点仍然是全局渐进稳定的.本文第二节首先给出了模型(1.2)的基本再生数,并讨论了模型(1.2)平衡点的存在性;第三节在平衡点存在的情况下讨论其无病平衡点和地方病平衡点的全局渐进稳定性,从而得知疾病的灭绝与持久是由基本再生数控制的.最后通过数值模拟对得出的结果进行了验证.

2 基本再生数和平衡点的存在性

首先给出模型(1.2)的一个控制疾病持久与灭绝的临界值[10]–基本再生数[11]

设Nt=St+It+Rt,将模型(1.2)中的三个方程相加可得

两边同时对t积分

因此

根据上述推理现在定义一个有界集Γ:

下面继续讨论地方病平衡点 E∗=(S∗,I∗,R∗)的存在性.设E∗=(S∗,I∗,R∗)是下列方程组的一个解

对方程组(2.3)求解可得

且

令

根据假设(H1)可知F(I)是一个减函数,因此

若R0＞1,F(I)=0存在一个正解I∗,则模型(1.2)存在唯一的地方病平衡点E∗=(S∗,I∗,R∗),其中

3 平横点的全局渐进稳定性

下面讨论模型(1.2)的无病平衡点E0和地方病平衡点E∗的全局渐进稳定性.

定理3.1 若R0≤1,则模型(1.2)的无病平衡点是全局渐进稳定的.

证 定义一个Lyapunov函数

∈足够小.则有

令

因此

分析 θ1Λ(β1− β2f(I))− µ(1+αI2)(θ1(µ + ν+ δ)− θ2ν),有

将上面的不等式代入(3.2)式,可以得到

(1)当R0＜1时,

(2)当R0=1时,θ2=0,从(3.3)式可以得到当且仅当根据LaSalle不变原理,模型(1.2)的任意解都收敛到B,其中B⊂{(S,I,R):,I=0,R=0}是模型(1.2)的最大不变子集,即B={E0}是一个单点集.因此,当R0≤1时,在有界集Γ中,E0是全局渐进稳定的.证毕.

定理3.2 当R0＞1时,模型(1.2)存在唯一的地方病平衡点E∗=(S∗,I∗,R∗)是全局渐进稳定的,且E0是不稳定的.

证 将E∗代入模型(1.2)的Jacobian矩阵

J(E∗)的特征多项式为

为了计算b1,b2,b3的值,令

其中

进一步计算可得

因为

由于b1＞0,b2＞0且b3＞0,则J(E∗)的三个特征值的实部都非负,由Routh-Hurwitz准则可知,E∗是局部渐进稳定的.

令Nt=St+It+Rt,把(1.2)式中的三个方程相加得则方程(1.2)变为

且N∗=S∗+I∗+R∗.接下来继续证明模型(1.2)的地方病平衡点E∗是全局渐进稳定的,则只需要证明模型(3.6)的解(N∗,I∗,R∗)是全局渐进稳定的.考虑下面的Lyapunov函数

其中k1和k2是正常数.V的倒数为

根据LaSalle渐进稳定定理[12,13],可知模型(1.2)的地方病平衡点E∗是全局渐进稳定的,又

4 数值分析

下面给出例子对以上结论进行验证.选取一个满足条件(H1)的函数模型(1.2)可化为

下面对(4.2)式中的参数赋值计算如下表2.

表2:各参数取值计算

通过MATLAB软件运算作图如图1.

图1:对于三组不同的参数值,模型(4.1)的解S(t),I(t)和R(t)的时间序列图

分析上图可知在图1(a)中,因为R0＞1,计算得

选取不同的初值时,S,I,R最后都将趋于唯一的地方病平衡点E∗,这说明疾病是持久的;在图1(b)中,因为R0=1,计算得E0=(2.0000,0,0),此时E∗是不存在的,通过选取不同的初值,S最后都将趋于=2.0000,而I和R逐渐趋于0,即疾病灭绝;在图1(c)中,因为R0=0.5333＜1,计算得E0=(1.3333,0,0),此时,E∗也是不存在的,且选取不同的初值,S最后都将趋于=1.3333,而I和R逐渐趋于0,即疾病灭绝.

综上所述,无病平衡点E0=(,0,0)是恒存在的,当R0＞1时,存在一个地方病平衡点E∗是全局渐进稳定的,而当R0≤1时,疾病灭绝.

5 结论

流行病的传播给人们的生活带来了巨大的损失和伤害.因此,本文讨论了一类确定性SIRS传染病模型的持久与灭绝:当基本再生数R0≤1时,疾病是灭绝的;当R0＞1时,疾病是持久的.然而在现实生活中,流行病不可避免地受到随机因素的影响,所以下一步将要讨论不确定环境下随机扰动对模型(1.2)的影响.