留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

一类接种率受媒体报道影响的传染病模型分析

上一篇

下一篇

王艳妮, 刘贤宁. 一类接种率受媒体报道影响的传染病模型分析[J]. 西南师范大学学报(自然科学版), 2020, 45(5): 1-6. doi: 10.13718/j.cnki.xsxb.2020.05.001
引用本文: 王艳妮, 刘贤宁. 一类接种率受媒体报道影响的传染病模型分析[J]. 西南师范大学学报(自然科学版), 2020, 45(5): 1-6. doi: 10.13718/j.cnki.xsxb.2020.05.001
Yan-ni WANG, Xian-ning LIU. Analysis of a Kind of Epidemic Model with Vaccination Rate Related to Media Coverage[J]. Journal of Southwest China Normal University(Natural Science Edition), 2020, 45(5): 1-6. doi: 10.13718/j.cnki.xsxb.2020.05.001
Citation: Yan-ni WANG, Xian-ning LIU. Analysis of a Kind of Epidemic Model with Vaccination Rate Related to Media Coverage[J]. Journal of Southwest China Normal University(Natural Science Edition), 2020, 45(5): 1-6. doi: 10.13718/j.cnki.xsxb.2020.05.001

一类接种率受媒体报道影响的传染病模型分析

  • 基金项目: 国家自然科学基金项目(11671327)
详细信息
    作者简介:

    王艳妮(1996-), 女, 硕士研究生, 主要从事生物数学及动力系统理论及其应用研究 .

    通讯作者: 刘贤宁, 教授, 博士研究生导师
  • 中图分类号: O175

Analysis of a Kind of Epidemic Model with Vaccination Rate Related to Media Coverage

计量
  • 文章访问数:  1966
  • HTML全文浏览数:  1966
  • PDF下载数:  307
  • 施引文献:  0
出版历程
  • 收稿日期:  2019-07-18
  • 刊出日期:  2020-05-20

一类接种率受媒体报道影响的传染病模型分析

    通讯作者: 刘贤宁, 教授, 博士研究生导师
    作者简介: 王艳妮(1996-), 女, 硕士研究生, 主要从事生物数学及动力系统理论及其应用研究
  • 西南大学数学与统计学院, 重庆 400715
基金项目:  国家自然科学基金项目(11671327)

摘要: 研究了一类接种率受媒体报道影响的SIR传染病模型,同时考虑到媒体报道延迟对模型动力学性态的影响.首先计算了模型的基本再生数R0R0 < 1时,利用LaSalle不变集原理得到了无病平衡点的全局稳定性;R0>1时,研究了地方病平衡点的局部稳定性.根据媒体报道是否延迟,分别讨论了以接触率和时滞作为分支参数,系统产生Hopf分支的条件.

English Abstract

  • 传染病作为人类的一项重大威胁,数学模型的建立始终是从理论上预测和控制传染病传播的重要方法[1-4].近年来,媒体作为当今时代的重要信息来源,影响着公众对疾病的认识和防范.

    媒体对传染病的报道会使公众提高意识、加强防备,因此在传染病模型中考虑信息反馈更符合实际[5-6].其中媒体报道的信息主要分为两种反馈方式:一种是对公众行为的直接影响,例如减少外出,戴口罩等,由此减小有效接触率[5, 7];另一种是对接种的影响,父母会根据媒体对传染病的报道判断是否对孩子进行接种[6, 8],从而影响接种率.事实上,在自愿接种的前提下,当媒体对传染病的报道较少时,人类预测到的感染风险也较低,导致接种的可能性降低[6, 9].文献[9]建立了一类接种率受媒体报道影响的传染病模型,其中报道的信息只与当前的疾病流行率相关,讨论了系统的平衡点全局稳定的条件.由于人类获取到的信息会流失,同时考虑到媒体对传染病报道的延迟,本文建立如下模型

    初值条件为

    其中:SIR代表易感染人群、已感染人群和免疫人群的密度;M为媒体报道的信息量;Λ为易感染人群的常数输入率;β为有效接触率;p0为不受媒体报道影响的接种率,p(M)则表示与媒体报道相关的接种率;d为人群的自然死亡率;r为感染者的恢复率;α0I为媒体依据传染病的流行情况而报道的的信息增长量;τ为媒体对病情报道的延迟时滞;λ0为信息的流失率.在本文的研究中,设p(M)满足下列条件

    1) p(0)=0;

    2) p′(M)>0,M>0;

    3) $\left(\frac{p(M)}{M}\right)^{\prime} \leqslant 0$.

    由于系统(1)中的第一、二、四个方程与R无关,因此,系统(1)的动力学性质等价于下列系统

    X=C([-τ,0],$\mathbb{R}$+3),由解的存在唯一性定理,对任意的ϕ∈X,系统(2)存在唯一解.

  • 本节研究系统(2)在条件(1)下解的非负有界性,记

    则由下面的定理可知Γ为系统(2)的正向不变集.

    定理1  系统(2)在条件(1)下的解(S(t),I(t),M(t))始终非负且有界.

      首先证明解的非负性.由系统(2)的第二个方程得

    I(0)≥0得对任意的t>0,有I(t)≥0.令t1=inf{tt>0,S(t)=0},由系统(2)的第一个方程得

    则存在充分小的ε1>0,在(t1-ε1t1)上有S(t) < 0,这与(0,t1)上有S(t)>0相矛盾,因此对任意的t>0,S(t)>0.同理可证对任意的t>0,有M(t)>0.

    下面证明解的有界性.令N=S+I,则根据解的非负性有

    由比较定理得

    由系统(2)的第三个方程得

    因此

    证毕.

  • 系统(2)始终存在一个无病平衡点E0=$\left(\frac{\mathit{\Lambda}}{p_{0}+d}, 0,0\right)$.计算基本再生数可得

    定理2  当R0>1时,系统(2)存在唯一地方病平衡点E*=(S*I*M*),其中$S^{*}=\frac{d+r}{\beta}$,M*=$\frac{\alpha_{0} I^{*}}{\lambda_{0}}$.

      令系统(2)的右端等于0,若I≠0,计算得$S=\frac{d+r}{\beta}$,$M=\frac{\alpha_{0} I^{*}}{\lambda_{0}}$,整理得

    由于

    因此,当R0>1时,g(0)>0,方程(3)存在唯一正实根,则系统(2)存在唯一的地方病平衡点E*=(S*I*M*).证毕.

  • 定理3  对一切τ≥0,当R0 < 1时,系统(2)的无病平衡点E0局部渐近稳定,当R0>1时不稳定.

      系统(2)关于无病平衡点E0的线性化矩阵为

    可得系统(2)关于无病平衡点E0的特征方程

    因此无病平衡点的稳定性与时滞τ无关,当且仅当R0 < 1时,所有特征值有负实部,由Routh-Hurwitz判据得,R0 < 1时,E0局部渐近稳定,R0>1时,E0不稳定.证毕.

    定理4  对一切τ≥0,当R0 < 1时,系统(2)的无病平衡点E0Γ内全局渐近稳定.

      构造Lyapunov函数$V=\frac{1}{2} I^{2}$,V沿着系统(2)轨线的全导数为

    R0≤1时,$\frac{\mathrm{d} V}{\mathrm{d} t} \leqslant 0$,且$\left\{x \in \mathit{\Gamma} | \frac{\mathrm{d} V}{\mathrm{d} t}=0\right\}=\left\{E_{0}\right\}$,由LaSalle不变集原理得,无病平衡点E0Γ内全局渐近稳定.证毕.

  • 下面研究地方病平衡点的局部稳定性,将系统(2)关于E*线性化得

    其中Y(t)=(SIM)T

    系统(2)关于E*的特征方程为

    其中

    定理5  τ=0时,若1 < R0 < R*,则系统(2)的地方病平衡点E*局部渐近稳定.

      τ=0时,系统(2)关于E*的特征方程为

    显然,Bi>0(i=1,2,3,4).由于$I^{*}=\frac{\lambda_{0}}{\alpha_{0}} M^{*}$,$\left(\frac{p(M)}{M}\right)^{\prime} \leqslant 0$,且由方程(3)得

    因此

    若$\left(\lambda_{0}+\frac{\mathit{\Lambda} }{S^{*}}\right)$$\left(\frac{\mathit{\Lambda}}{S^{*}} \mathit{\lambda} _{0}\right)$$-\beta \mathit{\Lambda} \lambda_{0}+(d+r)\left(p_{0}+d\right) \lambda_{0}>0$,即R0 < R*时,有B1B2-(B3+B4)>0,由Routh-Hurwitz判据得,地方病平衡点E*局部渐近稳定.证毕.

    接下来研究τ=0时,固定除β外的其他参数,系统(2)以β为参数产生Hopf分支的情况.假设

    定理6  τ=0时,设β*为方程B1B2-(B3+B4)=0的根.若假设(L1)成立,则β < β*时,E*局部渐近稳定,β>β*时,E*变为不稳定,且系统(2)在β=β*处产生Hopf分支.

      设B1B2-(B3+B4)=0时,β=β*,此时方程(6)变为

    因此,方程(6)必有一对纯虚根$\pm \sqrt{B_{2}}$i,下面讨论系统(2)在β=β*处产生Hopf分支的横截条件.将方程(6)关于β求导得

    将纯虚根$\sqrt{B_{2}}$i代入得

    当假设(L1)成立时,$\frac{\mathrm{d}({Re}(\mathit{\Lambda}))}{\mathrm{d} \beta}>0$,则β < β*时,E*局部渐近稳定,β>β*时,E*不稳定,且系统(L1)在β=β*处产生Hopf分支.

    接下来讨论B1B2-(B3+B4)>0时,时滞τ对系统(2)的影响.若τ>0,地方病平衡点E*稳定性改变的必要条件是方程(4)存在纯虚根.假设方程(4)存在纯虚根iθ(θ>0).由方程(4)得

    整理得

    其中

    l=θ2,则方程(8)变为

    接下来讨论方程(9)存在正实根的条件.令

    根据文献[10],有如下结论

    引理1[10]1)若ω3 < 0,则方程(9)至少存在一个正根;

    2) 若ω3>0,Δ>0,则方程(9)存在正根当且仅当z1>0且h(z1) < 0;

    3) 若ω3>0,Δ<0,则方程(9)存在正根当且仅当存在z*∈{z1z2z3}使得z*>0且h(z*)≤0.

    假设引理1中保证方程(9)存在正根的任意一个条件成立,不妨设lkk=1,…,k0(1≤k0≤3)为方程(9)的正根,则方程(8)存在相应的正根$\theta_{k}=\sqrt{l_{k}}$.令

    定义

    Λ(τ)=ξ(τ)+iθ(τ)为方程(4)的根,满足ξ(τ*)=0,θ(τ*)=θ0.假设

    定理7  τ>0时,若B1B2-(B3+B4)>0,且条件(L2)成立,则地方病平衡点E*对于τ < τ*局部渐近稳定,对于τ>τ*变为不稳定,且系统(2)在τ=τ*处产生Hopf分支.

      根据以上讨论,现在研究系统(2)在τ=τ*处产生Hopf分支的横截条件.将方程(4)关于时滞τ求导得

    利用式(7)可得

    当条件(L2)成立时,$\frac{\mathrm{d}(\operatorname{Re}(\lambda))}{\mathrm{d} \tau}>0$,因此地方病平衡点E*对于τ < τ*局部渐近稳定,对于τ>τ*变为不稳定,且系统(2)在τ=τ*处产生Hopf分支.证毕.

  • 研究了一类接种率受媒体报道影响的传染病模型,且考虑了媒体对传染病的报道存在延迟的情况.讨论了无病平衡点和地方病平衡点的稳定性,发现无病平衡点的全局稳定性不受时滞影响,当R0 < 1时,无病平衡点全局渐近稳定.而对于地方病平衡点,一方面,若时滞不存在,即媒体根据当前的疾病流行情况进行报道时,系统会以接触率为参数产生Hopf分支;另一方面,若存在时滞,随着时滞的增加,在地方病平衡点附近也会产生周期解.时滞越长,周期解的振荡幅度越大,这意味着传染病越难控制.因此,媒体对传染病病情的及时报道将有利于对传染病的控制.

参考文献 (10)

目录

/

返回文章
返回