1. 前言
自从Kermack与McKendrick在1927年文献 [1] 中建立易感(Susceptible)-染病(Infective)-康复(Reco- vered)模型(简称SIR模型)以来,传染病动力学发展迅速,大量的传染病数学模型被应用于分析各种各样的传染病,这些传染病模型大多数是用于研究各种传染病的一般传播规律,对人们研究传染病的传播机制,流行规律和防控理论具有重要的理论意义。近年来在众多学者的探索研究下,涌现出一大批关于传染病模型动力学性质研究的优秀成果 [2] [3] [4] [5] [6] [8] - [16] ,对于连续的SIR模型动力学性质研究中,也涌现出许多优秀的研究成果 [2] - [6] [14] [15] [16] 。对于连续传染病模型而言,不仅仅要研究模型的稳定性,更要研究其分岔性质。因为稳定性仅仅说明了模型所代表的系统受外界干扰时模型解的性质长时间的变化趋势 [7] [10] [14] ,而研究模型分岔性质 [19] 可以得到模型所代表的系统的拓扑结构的变化与模型系数参数的关系,从而得到传染病模型所代表的传染病的相关性质 [11] [12] [13] [15] [16] 。所以研究传染病模型的分岔性质对研究传染病的传播机制与流行规律是非常必要的。
在2008年,Xinzhu Meng和Lansun Chen在文献 [10] 中,根据传统的免疫性连续传染病SIR模型:
(1)
建立并研究了一类具有水平和垂直传播的新的免疫性连续传染病SIR模型
(2)
无病平衡点处的局部稳定性和地方病平衡点处的全局渐进稳定性,并且得出系统不存在极限环 [10] 。在模型(2)中分别代表易感者、染病者、康复者。是有效接触率,是出生率,是死亡率,模型(2)假设在种群中出生率等于死亡率,是康复率,分别代表易感者群体和染病者群体的后代比例,且满足,表示在易感群体中成功免疫的比例,包括成年人。在本文我们首先给出平衡点类型与模型(2)系数参数的关系,进而利用中心流形定理得到了SIR模型跨临界分岔的性质与正规形,最后给出模型在无病平衡点处发生跨临界分岔的生物学解释。
2. 平衡点的类型
在本节,我们在参数空间,其中上讨论平衡点类型与传染病模型系数参数的关系,确定平衡点的双曲和非双曲性质和平衡点的类型。在文献 [10] 中假设种群数量恒为常数,即,实施尺度变换,有
(3)
将(3)代入(2)有
(4)
系统(4)有两个平衡点,第一个平衡点是没有感染者的群体,称为无病平衡点,其中。第二个平衡点是有感染者的群体,称为地方病平衡点,其中:。其中为基本再生率,详见文献 [10] 。从上面实施的尺度变换和传染病模型(4)的生物学意义可知,,故系统(4)的解总是有界的。
定理2.1无病平衡点是非双曲的当且仅当位于直线,否则无病平衡点将是如下几种类型;
i) 当时,无病平衡点是不稳定的鞍点;
ii) 当时,无病平衡点是稳定的结点。
证明 系统(4)在的雅可比矩阵为
其特征值为
由模型假设可知,故,当时,,所以无病平衡点是不稳定的鞍点;当时,,所以无病平衡点是稳定的结点。
定理2.2 地方病平衡点不存在非双曲的情形,地方病平衡点将是如下几种类型;
i) 当时,地方病平衡点是稳定焦点;
ii) 当,时,地方病平衡点是稳定结点;
iii) 当,时,地方病平衡点是不稳定的鞍点。
证明系统(4)在处的雅可比矩阵为
其特征方程为
整理得
当时,,,所以地方病平衡点是稳定焦点;当时,记,并且与同号,由韦达定理知,,故当时,即,,所以平衡点是稳定结点,当时,即,与一正一负,所以地方病平衡点是不稳定的鞍点。当地方病平衡点是非双曲情形时有,,,所以与不可能同时为0,故当且仅当时地方病平衡点是非双曲的,此时地方病平衡点退化成无病平衡点,两个平衡点重合,此时系统(4)只有一个平衡点,故地方病平衡点不存在非双曲情形,证毕。
3. 系统(4)的跨临界分岔
在本节,我们将用中心流形定理来研究系统(4)在无病平衡点处跨临界分岔的性质与正规形,下面的两个引理来自于文献 [18] 。
引理3.1 [18] 若系统
(5)
满足条件:矩阵的特征值为0,矩阵的特征值小于0,且
则系统(5)存在一中心流形,在局部充分小范围内可表示为
进一步,限制在中心流形上的系统(5)在局部充分小范围可以表示为
引理3.2 [18] 设参数族是一维系统
(6)
满足条件:,。若
则系统(6)在点出现跨临界分岔。
定理3.1当,在平衡点处将发生跨临界分岔.具体地讲,对较小的扰动使得时,系统(4)出现两个平衡点:一个是稳定的无病平衡点,另一个是不稳定的负的地方病平衡点,它们在时发生粘结合成一点;当时,系统(4)也出现两个平衡点:一个是不稳定的无病平衡点,另一个是稳定的正地方病平衡点。因此,当在附近发生微小变化时,系统(4)平衡点的个数发生了改变,并且平衡点与的稳定性发生了交替,即系统(4)在无病平衡点处发生跨临界分岔,且系统(4)在无病平衡点处的正规形为 [17] 。
证明 对于位于直线,有,为了体现系统(4)在无病平衡点处的跨临界分岔对参数的依赖,我们把写成;施行如下坐标变换,易知系统(4)在无病平衡点处的雅可比矩阵为,
及特征值和对应的特征向量为
(7)
记,特征值和对应的特征向量可以写成如下形式
(8)
由特征向量组(8)我们可以得到如下变换,
(9)
由变换(9)可以得到如下系统
(10)
把系统(10)写成悬挂式
(11)
由引理3.1可知在附近可以通过研究中心流形上一参数系统来确定平衡点处的稳定性与分岔性质,由引理3.1可知中心流形形式如下;
由引理3.1假设中心流形有如下形式;
(12)
由引理3.1知可以通过(13)来计算中心流形(12),其中表示次数大于或等于3的项;
(13)
比较系数得
解得;故系统(11)在附近中心流形为,代入系统(11)得
(14)
可以验算如下式子成立
(15)
由引理3.2可以知道系统(4)在发生跨临界分岔,式子(14)表明在附近所有高阶项不影响其分岔性质,同时表明在附近系统(12)的拓扑结构与(16)的拓扑结构是局部拓扑等价的;
(16)
式子(16)可以看作系统(4)在无病平衡点处跨临界分岔的正规形。
4. 跨临界分岔的生物学解释
由定理3.1可以知道,在处系统(4)的拓扑结构容易发生质的变化,在时,地方病平衡点是不稳定的且是负的,表明传染病并没有在种群中流行;当时,地方病平衡点是稳定的,即传染病会在种群中稳定流行.所以系统(4)的拓扑结构与传染病的流行与否有着非常紧密的联系,人们可以通过研究传染病模型的拓扑结构并且通过控制传染病模型的有关参数来达到控制流行病的目的,对传染病的传播机理,流行规律的研究与流行病的预防和控制具有重要理论意义。
基金项目
国家自然科学基金项目(11561019);广东省创新强校科技重大项目(2014KZDXM065);广东省大学生科技创新重点项目(pdjh2016a0301);广东省大学生科技创新项目(pdjh2017b0314);岭南师范学院攀峰计划项目。
参考文献