为了描述淋病和其他传染病的传播动态,文献[1, 2, 3, 4, 5] 建立和研究了一类非线性时滞微分系统
文献[6, 7, 8, 9, 10]的作者们在系统(1.1)中考虑了连续分布时滞,得到如下系统
系统(1.2)比系统(1.1)更符合实际,其中ci(t)是第i个区域无病易感者数量 (如果不被传染就不会变成染病者), xi(t)表示第i个区域在t时刻染病者数量,βij(t) 是第j个区域的染病者对 第i个区域易感者在t时刻的传染率,ai(t)表示第i个区域染病者的恢复率, τij(t)≥0表示病毒 在体内的潜伏期,即染病者发病到感染其他人的时段.此外,此模型染病者不会死亡并且不具有免疫力.
收敛性在动力系统行为特征中扮演着关键的角色,种群动力学模型解的收敛性研究已经取 得了广泛的结果(见文献[11, 12, 13]). 特别地,恢复率函数ai(t)为不震动时,即
系统(1.1)和(1.2)已经建立了非震动条件下确保解收敛的若干充分性条件, 其结论可见文献[1, 2, 3, 4, 5, 14, 16] 和 [6, 7, 8, 9, 10, 15].
另一方面,文献[17]指出,随着季节波动,震动系数出现在线性种群动力学模型方程中, 在某些季节死亡率(收获率) 可能会大于出生率,此外由于缺乏药物和及时治疗,某些突发传染病会出现负恢复率. 那么,一个问题自然产生了: 寻求确立带有震动恢复率的系统(1.1)和(1.2)的全局指数稳定性的条件.基于上述讨论, 本文中我们将研究如下 带有时变时滞和连续分布时滞传染病模型的指数稳定性
其中 ai:R→R,ci, βij,˜ci, ˜βij,τij:R→[0,+∞) 和 Kij:(−∞,0]→[0,+∞)是有界连续函数,且 Kij(t)eκt 在[0,+∞)可积, κ为一个确定正的常数,i, j∈J. 显然,方程(1.1)和(1.2)是模型(1.4) 的特殊情形.
基于模型(1.4)的生物学意义,我们只考虑如下初始条件
这里
C+=CB((−∞,0],[0,+∞))×CB((−∞,0],[0,+∞))×⋯×CB((−∞,0],[0,+∞))⏟n,
且 CB((−∞,0],[0,+∞)) 是有界连续函数的集合.
本文主要目的是在舍弃假设(1.3)的前提下,给出初值问题(1.4)和(1.5)的全局指数稳定性的充分 条件. 据我们所知,具有震动恢复率的模型(1.4)的收敛性问题未见他人研究,本 文的结果是全新的. 本文第三部分数值模拟例子很好地说明了我们结果的有效性.
本文假设以下条件成立.
(H1) 存在一个正的常数Fi 和一个连续函数 a∗i:R→(0, +∞) 使得
e−∫tsai(u)du≤Fie−∫tsa∗i(u)du t,s∈R, t−s≥0, i∈J;
(H2) 存在正的常数 ξ1,ξ2,⋯,ξn 使得
supt∈R{−a∗i(t)+ξ−1iFici(t)n∑j=1βij(t)ξj+ξ−1iFi˜ci(t)n∑j=1ξj˜βij(t)∫0−∞Kij(v)dv}<0, i∈J.
为了方便,我们介绍一些符号.
|x|=(|x1|, |x2|, ⋯,|xn|)T,‖x‖=max1≤i≤n|xi|,g+=supt∈R|g(t)|,g−=inft∈R|g(t)|.
记初值问题(1.4)和(1.5)的解为x(t;0,φ)=(x1(t;0,φ),x2(t;0,φ),⋯,xn(t;0,φ))T, 并记其最大右存在区间为[0,η(φ)).
定理 2.1 假设 (H1) 和(H2) 成立. 则,当xi(t;0,φ)>0,t∈[0,η(φ)), i∈J, 集合 {xi(t;0,φ):t∈[0,η(φ)),i∈J} 有界,且 η(φ)=+∞. 此外,当 t→+∞时,x(t;0,φ)=(x1(t;0,φ),x2(t;0,φ),⋯,xn(t;0,φ))T 收敛于(0,0,⋯,0)T.
证 令
(x1(t),x2(t),⋯,xn(t))T=(x1(t;t0,φ),x2(t;t0,φ),⋯,xn(t;t0,φ))T.
首先证明
利用反证法,假设(2.1)式不成立, 那么必存在 i∈J 和 T∈[0,η(φ)) 使得
由(1.4)式,有
x′i(s)=−{ai(s)+n∑j=1[βij(s)xj(s−τij(s))+˜βij(s)∫0−∞Kij(v)xj(s+v)dv]}xi(s)+ci(s)n∑j=1βij(s)xj(s−τij(s))+˜ci(s)n∑j=1˜βij(s)∫0−∞Kij(v)xj(s+v)dv, s∈[0,T]
和
0=xi(T)=xi(0)e−∫T0{ai(u)+n∑j=1[βij(u)xj(u−τij(u))+˜βij(u)∫0−∞Kij(v)xj(u+v)dv]}du+∫T0e−∫Ts{ai(u)+n∑j=1[βij(u)xj(u−τij(u))+˜βij(u)∫0−∞Kij(v)xj(u+v)dv]}du[ci(s)n∑j=1βij(s)xj(s−τij(s))+˜ci(s)n∑j=1˜βij(s)∫0−∞Kij(v)xj(s+v)dv]ds≥xi(0)e−∫T0{ai(u)+n∑j=1[βij(u)xj(u−τij(u))+˜βij(u)∫0−∞Kij(v)xj(u+v)dv]}du>0,
这是一个矛盾,因此(2.1) 式成立.
其次证明 {xi(t):t∈[0,η(φ)),i∈J} 有界且 η(φ)=+∞.
令y(t)=(y1(t),y2(t),⋯,yn(t))T=(ξ−11x1(t),ξ−12x2(t),⋯,ξ−1nxn(t))T. 由 (1.4)式 有
y′i(t)=−ai(t)yi(t)+ξ−1i(ci(t)−xi(t))n∑j=1βij(t)xj(t−τij(t)) +ξ−1i(˜ci(t)−xi(t))n∑j=1˜βij(t)∫0−∞Kij(s)xj(t+s)ds, i∈J.
由 (H2),选择一个常数λ∈(0,min{κ,mini∈Ja∗i −}) 使得
设
对任意 ε>0,有‖y(t)‖<(‖φ‖ξ+ε)e−λt<M(‖φ‖ξ+ε)e−λt, t∈(−∞,0], 这里 M=n∑i=1Fi+1..
下面证明
否则,必存在 i∈J 和 θ∈(0,η(φ)) 使得
注意到
在(2.7)式的两边乘以 e∫s0ai(u)du,并在[0,t]上积分 ,得
那么,由(2.3),(2.6)和(2.8)式,有
|yi(θ)|=yi(θ)=yi(0)e−∫θ0ai(u)du+∫θ0e−∫θsai(u)du[ξ−1i(ci(s)−xi(s))n∑j=1βij(s)xj(s−τij(s)) +ξ−1i(˜ci(s)−xi(s))n∑j=1˜βij(s)∫0−∞Kij(v)xj(s+v)dv]ds≤yi(0)e−∫θ0ai(u)du+∫θ0e−∫θsai(u)du[ξ−1ici(s)n∑j=1βij(s)xj(s−τij(s))+ξ−1i˜ci(s)n∑j=1˜βij(s)∫0−∞Kij(v)xj(s+v)dv]ds≤yi(0)Fie−∫θ0a∗i(u)du+Fi∫θ0e−∫θsa∗i(u)du[ξ−1ici(s)n∑j=1βij(s)ξjyj(s−τij(s)) +ξ−1i˜ci(s)n∑j=1˜βij(s)∫0−∞Kij(v)ξjy(s+v)dv]ds≤(‖φ‖ξ+ε)Fie−∫θ0a∗i(u)du +Fi∫θ0e−∫θsa∗i(u)du[ξ−1ici(s)n∑j=1βij(s)ξjM(‖φ‖ξ+ε)e−λ(s−τij(s)) +ξ−1i˜ci(s)n∑j=1˜βij(s)∫0−∞Kij(v)ξjM(‖φ‖ξ+ε)e−λ(s+v)dv]ds≤M(‖φ‖ξ+ε)e−λθ{FiMe−∫θ0(a∗i(u)−λ)du +∫θ0e−∫θs(a∗i(u)−λ)duFi[ξ−1ici(s)n∑j=1βij(s)ξjeλτ+ij +ξ−1i˜ci(s)n∑j=1ξj˜βij(s)∫0−∞Kij(v)e−λvdv]ds}≤M(‖φ‖ξ+ε)e−λθ[FiMe−∫θ0(a∗i(u)−λ)du+∫θ0e−∫θs(a∗i(u)−λ)du(a∗i(s)−λ)ds]≤M(‖φ‖ξ+ε)e−λθ[1−(1−FiM)e−∫θ0(a∗i(u)−λ)du]<M(‖φ‖ξ+ε)e−λθ,
这与(2.6)式的第一个方程矛盾,所以(2.5) 式成立.令 ε⟶0+,由(2.5)式知
‖y(t)‖≤M‖φ‖ξe−λt≤M‖φ‖ξ, t∈[0,η(φ)),
说明集合{xi(t):t∈[0,η(φ)),i∈J} 有界. 根据文献[18]中的定理2.3和2.4,易得η(φ)=+∞. 因此,
‖x(t)‖≤maxi∈Jξi‖y(t)‖≤maxi∈JξiM‖φ‖ξe−λt, t∈[0,+∞),
这就完成了定理2.1的证明.
下面处理与模型(1.4)类似的如下系统
其中 ai:R→R,ci, βij,˜ci, ˜βij,τij,ˉτij:R→[0,+∞) 和 Kij:[−ˉτ+ij,0]→[0,+∞) 都是有界连续函数,i,j∈J.
令 τ=maxi,j∈J{τ+ij,ˉτ+ij}, 和
˜Cn+=C([−τ,0],[0,+∞))×C([−τ,0],[0,+∞))×⋯×C([−τ,0],[0,+∞))⏟n.
与系统(2.9)相关的初始条件形式如下
与定理2.1的论证类似,我们得到
定理 2.2 假设(H1) 成立, 存在正的常数 ξ1,ξ2,⋯,ξn 使得
supt∈R{−a∗i(t)+ξ−1ici(t)n∑j=1βij(t)ξj+ξ−1i˜ci(t)n∑j=1˜βij(t)∫0−ˉτij(t)Kij(v)dvξj}<0, i∈J.
设初值问题(2.9)和(2.10)的解为 x(t;0,φ)=(x1(t),x2(t),⋯,xn(t))T ,并记其最大右存在区间为 [0,η(φ)). 则当 xi(t)≥0, t∈[0, η(φ)),i∈J, 集合 {xi(t):t∈[0,η(φ)),i∈J} 有界,且 η(φ)=+∞. 此外,当 t→+∞时,x(t;0,φ) 指数收敛于 (0,0,⋯,0)T.
注 2.1 当 ai(t)≡a∗i(t)(i∈J) 不震动时,在
的假设下,系统(1.1)和(2.9)的指数收敛性分别在文献[14, 16]和[15]已经得证. 显然,(2.11)式是条件(H2) 在 a(t)≡a∗(t)和 ˜βij(t)≡0的特殊情形, (2.12)式是条件(H2) 的当a(t)≡a∗(t) 和 βij(t)≡0的特殊情形. 所以文献[14, 16] 和[15]的结果分别是本文定理2.1和2.2的推论.因此, 本文的结论改进和扩展了先前文献的结果.
例 3.1 考虑如下带有时滞和震动恢复率的传染病模型
显然,
a1(t)=18+20sin2000t, a2(t)=18+20cos2000t, a∗1(t)=a∗2(t)=18,ξi=1.
e−∫tsai(u)du≤e1100e−18(t−s), t,s∈R, t−s≥0,i=1,2.
supt∈R{−a∗i(t)+ξ−1iFici(t)n∑j=1βij(t)ξj+ξ−1iFi˜ci(t)n∑j=1˜βij(t)∫0−∞Kij(v)dvξj}≤−18+e1100(2+2+2+2)<−1, ξi=1, i=1,2.
即系统(3.1)满足了定理2.1的所有条件.故当t→+∞时, 系统(3.1)具有初值φ∈C+条件下的所有解指数收敛于(0,0).数值模拟验证了此结果(见图 1).
注 3.1 系统(3.1)是一个带有震动恢复率a1(t)=18+20sin2000t 和 a2(t)=18+20cos2000t的简单传染病模型,很显然,文献[1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 14, 15, 16]和其中的参考文献的所有结果不能 直接获得系统(3.1)的所有解的指数收敛性, 这说明本文结论是全新的.