2019年7月1日出版
文章编号:1000⁃0887(2019)07⁃0713⁃15
AppliedMathematicsandMechanics Vol.40,No.7,Jul.1,2019
ⓒ应用数学和力学编委会,ISSN1000⁃0887
一类具有非线性发生率与时滞的非局部扩散
SIR模型的临界波的存在性
∗
张 秋1, 陈广生1,2
2.广西科技师范学院数学与计算机科学学院,广西来宾546199)
摘要: 研究了一类具有时滞的非局部扩散SIR传染病模型的行波解.首先,利用反证法证明了I是有界的,并根据I的有界性研究了波速c>c∗时行波解(波速大于最小波速的行波)的存在性.其次,利用c>c∗的行波的存在性结果证明了临界波(波速等于最小波速的行波)的存在性.最后,讨论了R0对临界波存在性的影响.
关 键 词: 行波解; 临界波速; 非局部扩散; 基本再生数中图分类号: O175.14 文献标志码: A
DOI:10.21656/1000⁃0887.390208
(1.西安电子科技大学数学与统计学院,西安710071;
引 言
有大量的结果[1⁃10].例如,Wang和Wu[11]研究了一般的Kermack⁃McKendrick模型行波解的存在性.Li等[12]研究了人口数量变化的非局部时滞SIR传染病模型.文献[13]研究了具有非线性发生率的时滞SIR模型等等.
ìï∂S(x,t)=d1(J∗S(x,t)-S(x,t))-f(S(x,t))g(I(x,t-τ)),ï∂tïï∂I(x,t)
=d2(J∗I(x,t)-I(x,t))+f(S(x,t))g(I(x,t-τ))-γI(x,t), í(1)∂tï
ï∂R(x,t)
=d3(J∗R(x,t)-R(x,t))+γI(x,t),ïï∂t
î
其中di>0(i=1,2,3)表示扩散系数,γ表示恢复率,τ≥0表示疾病的潜伏时间,J∗u(x,t)是关于空间x的标准卷积.J∗S(x,t)-S(x,t),J∗I(x,t)-I(x,t),J∗R(x,t)-R(x,t)分别表示由易感者、感染者、恢复者的扩散导致的净增长率.假设J,f和g满足条件:
(H1)f,g∈C((0,+∞),(0,+∞)),f(0)=g(0)=0,f′(S)>0,∀S≥0,且g′(I)>0,g″(I)≤0,∀I≥0;
∗
传染病模型是生物数学中一类典型的数学模型,关于传染病模型行波解的存在性研究已
最近,邹霞和吴事良[14]研究了以下具有时滞的非局部扩散SIR传染病模型:
收稿日期: 2018⁃07⁃26;修订日期: 2018⁃11⁃09
基金项目: 国家自然科学基金(面上项目)(11671315)
作者简介: 张秋(1989—),女,硕士生(通讯作者.E⁃mail:1204142234@qq.com);
陈广生(1979—),男,博士生(E⁃mail:cgswavelets@126.com).
713
714张 秋 陈 广 生
(H2)J∈C1(R),J(x)=J(-x)≥0,
(f(S-∞)g′(0))/γ>1时,系统(1)存在满足边界条件S(-∞)=S-∞,I(±∞)=0的行波解,并得到了当R0≤1或者R0>1,0<c<c∗时,系统(1)不存在行波解的结论,其中σ0=S-∞
在条件(H1)和(H2)下,邹霞和吴事良[14]证明了当c>max{c∗,(3/2)d2σ0},R0=
∫J(x)dx=1且J具有紧支撑.
R
Schauder不动点定理和双边Laplace变换建立了系统(1)的行波解的存在性与不存在性.文献[14]是在大传播速度(即c>max{c∗,(3/2)d2σ0})条件下得到系统(1)的行波解的存在性,文献[15]是在小传播速度(即c>c∗)条件下得到行波解的存在性.然而,文献[14⁃15]没有讨论当波速c=c∗时,系统(1)的行波解的存在性问题.众所周知,反应扩散系统在c>c∗参阅文献[9,16⁃18],其中Wu[18]在Fu等[6]研究的基础上,讨论了以下模型临界波的存在性:条件下行波解的存在性已有很多研究,但当c=c∗时行波解的存在性讨论很少,相关结果可ìï∂Sj(t)=d[S(t)+S(t)-2S(t)]-βS(t)I(t),
j+1j-1jjj
ï∂t í(2)
∂I(t)ïj
ï∂t=[Ij+1(t)+Ij-1(t)-2Ij(t)]+βSj(t)Ij(t)-γIj(t).î
Wu[18]通过讨论系统(2)中I的有界性,证得了在波速c>c∗条件下行波解的存在性,从而得到当c=c∗时行波解的存在性.最近,Yang等[17]利用文献[18]中的方法建立了以下模型临界波的存在性:
ìï∂S(x,t)=d1(J∗S(x,t)-S(x,t))-βS(x,t)I(x,t),ï∂tí
ï∂I(x,t)=d(J∗I(x,t)-I(x,t))+βS(x,t)I(x,t)-γI(x,t).
2ï∂t
î
∫
+∞-∞
J(y)|y|dy.随后,Zhang等[15]也研究了系统(1)的行波解的存在性问题,他们利用
(3)
注意到,系统(3)是系统(1)的一个特例,那自然要问:系统(1)是否存在具有临界波速的行波解?回答是肯定的.
为了研究系统(1)的临界波,仅需考虑系统(1)的子系统(由于系统(1)的第三个方程与前两个方程不是耦合的):
ìï∂S(x,t)=d1(J∗S(x,t)-S(x,t))-f(S(x,t))g(I(x,t-τ)),ï∂t í(4)
∂I(x,t)ï=d2(J∗I(x,t)-I(x,t))+f(S(x,t))g(I(x,t-τ))-γI(x,t).ï∂tî
令(S(x+ct),I(x+ct))=(S(z),I(z)),其中z=x+ct.将其代入系统(4)可得如下相应的行波系统:
{
cS′(z)=d1(J∗S(z)-S(z))-f(S(z))g(I(z-cτ)),
其中边界条件为
cI′(z)=d2(J∗I(z)-I(z))+f(S(z))g(I(z-cτ))-γI(z),
(5)(6)
(S,I)(-∞)=(S-∞,0),I(+∞)=0,
且S在+∞处的取值是任意的.
现在,给定f(S-∞)g′(0)>γ,定义最小传播速度:
一类具有非线性发生率与时滞的非局部扩散SIR模型的临界波的存在性
c∗
λ>0
inf
d2
∫
715
+∞-∞
J(y)e-λydy-d2+f(S-∞)g′(0)e-λcτ-γ
λ
.
本文受文献[17⁃18]的启发,主要研究系统(4)当c=c∗时行波解的存在性.因此,首先利用反证法证明I是有界的,并根据I的有界性研究波速c>c∗时行波解的存在性.其次,利用c>c∗时行波的存在性结果证明临界波的存在性.最后,讨论R0对临界波存在性的影响.
这里需要指出的是,文献[14]在证明行波解的存在性时,对建立I的有界性附加了一定
的条件,而本文使用文献[17⁃18]中的技巧,在没有限制的条件下,利用反证法证明了I的有界性,进而建立了当c>c∗时行波解的存在性.由于时滞的存在和函数g的性质的影响,在证明I的有界性时,很多证明细节与文献[17⁃18]是不同的.
本文的安排如下:第1节研究当c>c∗时行波解的有界性;第2节证明当c=c∗时行波解的存在性;第3节给出具体实例,验证了所得的临界波结果;最后,对全文内容进行总结.
1 行波解的有界性
(4)存在满足0<S<S-∞,I>0的解(S,I),且S(-∞)=S-∞,I(-∞)=0.I(z
∫J(y)I(Iz(z)y)dy,f(S(z))gI((z)
-
R
本节主要研究系统(4)行波解的有界性.由文献[14]可得,当R0>1和c>c∗时,系统引理1 若函数(S,I)是系统(5)波速为c>c∗的行波解且0<S<S-∞,I>0,则函数
-cτ))
,
I′(z)
在R上有界.I(z)
证明 首先,对方程(5)的第二个方程的两边同时除以I,可得
c
I′(z)I(z-y)f(S(z))g(I(z-cτ))
=d2J(y)-γ.dy-d2+I(z)I(z)I(z)R
-
R
因为0<S<S-∞,可得I′/I的有界性与有关.
为了方便起见,记u(z)
u(z)=
d2cd2
d2c
∫
(7)
-cτ))
的有界性
I(z
∫J(y)I(Iz(z)y)dy,f(S(z))gI((z)
记κ
d2c
和G(z)
所以G(z)是非减的且limz→-∞G(z)=0.由条件(H2)可知,r>0,其中r是suppJ的半径.因此,取0<2r0<r.对方程(8)从-∞到z积分,可得
G(z)≥κ
∫u(s)ds).因此,对G(z)关于z求导,有
G′(z)=(μ+u(z))G(z)≥κ∫J(y)e∫dyG(z)≥0.
expμz+
c∫
∫
R
J(y)e∫zJ(y)e∫z
∫
R
J(y)
I(z-y)f(S(z))g(I(z-cτ))
=dy-μ+
I(z)cI(z)
u(s)ds
I′/I和μ=(d2+γ)/c.则式(7)为
z-y
dy-μ+dy-μ.
z0
f(S(z))g(I(z-cτ))
≥
cI(z)
z-y
u(s)ds
R
(
z-yz
u(s)ds
R
(8)
∫∫
z-∞
R
J(y)e∫x
x-y
u(s)ds
dyG(x)dx=
716
κJ(y)eμy
R
则
然后,对方程(8)关于x从z-r0到z积分,可得
∫
∫∫G(x-y)dxdy≥κ∫J(y)e∫G(x-y)dxdy≥κr∫J(y)eG(z-r-y)dy,
z-∞μy
z
R
z-r0
0
μy
R
0
张 秋 陈 广 生
0-∞
J(y)eμy
G(z-r0-y)
G(z)
dy≤(κr0)-1.
(9)
G(z)∫
-G(z-r0-y)≥κ
κμy
∫z
ξ-r0∫
R
J(y)e∫x-y
x
u(s)ds
dyG(x)dx=
R
J(y)e
≥
κr-2r∫
zz-r0
G(x-y)dxdy0
0μy-∞J(y)eG(z-r0-y)dy≥ κr∫-2r0
0-∞
J(y)eμydyG(z+r0).因为r>2r-∫2r0
和∫
0-∞
J(y)eμy1
dy>0且
ω0
κr-2r0
0-∞
所以,∫
J(y)eμydy
,
对任意的 G(z+ξr∈R,有由于
0)≤ω0G(z).
∫
R
J(y)
I(Iz(-z)
y)
dy=
∫
R
J(y)eμy
G(Gz(-z)
y)
dy,
结合式(9)和(10)可得
J(y)I(∫Iz(-
R
z)y)dy=
∫ 0J(y)eμy
G(z-y)
+∞
J(y)eμy
G(Gz(-z)
y)
dy≤
∫
-∞
G(z)dy+
0J(y)e
G(0μy
∫∞
-∞
Gz(-z)
y)
dy+
(y)eμydy≤
ω00(y)e
μy
G(z-r0-y∫
+)
0
J-∞
G(z)
y+
+∞
0
J(y)eμydy≤
ωκr0
∫
J+
又因为0
∫
d+∞
0
J(y)eμydy.
∫
整理得
G(Gz)(z是非减的)≥G(z,-故
cτ).(10)
一类具有非线性发生率与时滞的非局部扩散SIR模型的临界波的存在性
则
因此
expμz+
((
∫u(s)ds)≥exp(μz-μcτ+∫
z0
717
z-cτ
0
u(s)ds,
)
expμcτ-exp
由条件(H1)知,对任意的t∈[0,I(z-cτ)],有
f(S(z))g(I(z-cτ))f(S(z))g(I(z-cτ))I(z-cτ)
= ≤
I(z)I(z-cτ)I(z)
zg(I(z-cτ))
f(S-∞)expu(s)ds≤
I(z-cτ)z-cτ
f(S-∞)g′(0)exp(μcτ),
(∫
∫
z-cτz
u(s)ds≥1.
)
z-cτz
u(s)ds≤exp(μcτ).
)
(11)
(∫
)
则
所以
f(S(z))g(I(z-cτ))
≤exp(μcτ)f(S-∞)g′(0).
I(z)|u(z)|≤κJ(y)
d2+γexp(cμτ)I(z-y)
+dy+f(S-∞)g′(0).
I(z)ccR
I′I(z-y)
因此,J(y)dy和是有界的.
II(z)R
下面介绍几个引理用于证明I的有界性.
引理2 任取ck∈(c∗,c∗+1)且(ck,Sk,Ik)是系统(5)波速为c=ck的一序列行波解,且0<Sk<S-∞,Ik>0.如果存在一个序列{zk}使得当k→+∞时Ik(zk)→+∞,则Sk(zk)→0.
∗
中,k∈N和ε>0.由式(5)的第一个方程,可得S′.因此,Sk(z)≥ε/2,∀zk(z)≤2d1S-∞/c∈[zk-δ,zk].由引理1得,存在实数L1>0满足|I′k/Ik|≤L1.因此,对所有的k∈N,有
∫
∫
□
证明 反证法.假设不成立,则存在一个序列{zk}(仍用zk表示)使得Sk(zk)≥ε,其
Ik(zk)Ik(z)
Ik(s)当k→+∞时,Ik(zk)→+∞,则
Ik(z)≥e-δL1Ik(zk)→+∞.又因为
z
=exp
{∫
zk
I′k(s)
Lδ
ds≤e1, ∀z∈[zk-δ,zk].
}
因此,当k→+∞,有
minIk(z-cτ)≥e-cτL1Ik(z)≥e-cτL1e-δL1Ik(zk)→+∞.注意到
z∈[zk-δ,zk]
Ik(z-cτ)
Ik(z)
=exp
{∫
zz-cτ
I′k(s)Ik(s)
cτL
ds≤e1.
}
g(Ik(z-cτ))-g(0)
Ik(z-cτ)-0
=g′(t0), t0∈[0,Ik(z-cτ)].
718张 秋 陈 广 生
根据条件(H1)可知0<g′(Ik(z-cτ))<g′(t0)<g′(0).则由式(5)的第一个方程得这与Sk(z)<S-∞矛盾.证毕.
z∈[zk-δ,zk]
maxS′k(z)≤
2d1S-∞
c∗
εö
-fæç÷g′(t0)Ik(z-cτ)→-∞.
2èø
□
引理3 若(S,I)是系统(5)波速为c>c∗时的行波解且0<S(ξ)<S-∞,I>0,如果
limsupz→+∞I(z)=+∞,那么limz→+∞I(z)=+∞.
∞
=+∞和liminfz→+∞I(z)<+∞.则存在点列{zk}k+=1满足limk→+∞zk
I(zk)→m.因此,对任意的k∈N,有I(zk)≤m+1.
再者,对每一个k∈N,假设存在一个点yk∈[zk,zk+1]使得I(yk)=maxz∈[zk,zk+1]I(z).
证明 假设m
又因为limsupz→+∞I(z)=+∞,则当k→+∞时,I(yk)→+∞.因为L1是引理1中给定的常数,所以,对所有的k∈N,有I(yk)>(m+1)eL1r,其中r是suppJ的半径.那么,当|y-yk|≤r时,有
ykI′(z)æö=expçdz÷≤eL1r, ∀k∈N,
I(y)èyI(z)ø
可得[yk-r,yk+r]⊂(zk,zk+1).因此,由式(5)的第二个方程可得出
0=cI′(yk)=
I(yk)
∫
f(S(yk))g(I(yk-cτ))-γI(yk).(12)由limk→+∞I(yk)=+∞和引理2,所以limk→+∞S(yk)=0.因此,limz→+∞I(z)=+∞.证毕.□一个可测函数且满足
的特征根.
cZ(z)=cψ=
引理4[19] 假设c>0,B(·)是一个连续函数满足B(±∞)
limz→±∞B(z).令Z(z)为
d2J(y)(I(yk-y)-I(y))dy+f(S(yk))g(I(yk-cτ))-γI(yk)≤
R
∫
因此,Z是一致连续和一致有界.再者,ψ±
定理1 若函数(S,I)是系统(5)波速为c>c∗时的行波解且S(-∞)=S-∞,I>0,则I在R上有界.
证明 反证法.假设I在R上无界,则limsupz→+∞I(z)=+∞.因此由引理2和引理3,可得I(+∞)=+∞和S(+∞)=0.
记u(x)=I′(x)/I(x).对式(5)的第二个方程两边同除以I,则有
d2d2+γf(S(z))g(I(z-cτ))I(z-y)
+ u(z)=J(y)dy-.
cRI(z)cI(z)
因为S(+∞)=0和S(-∞)=S-∞,所以由引理4可得limz→+∞u(z)存在,且满足
d2
定义连续函数:
∫J(y)e
R
j
∫
R
Jj(y)e∫z
-ψy
z-y
Z(s)ds
dy+B(z).
dy+B(±∞) (j=1,2)
limz→±∞Z(z)存在且是方程
(∫J(y)e
R
∫
-λy
dy-1-cλ-γ=0.
-λy
)
(13)
h(c,λ)=d2
(∫J(y)e
R
dy-1-cλ-γ.
)
一类具有非线性发生率与时滞的非局部扩散SIR模型的临界波的存在性719
因此,h满足h(0,c)=-γ<0且
∂h∂λ
∂2h
=-c<0,=d2J(y)y2e-λydy>0,limh(λ,c)=+∞. 2λ→+∞R∂λλ=0
所以方程(13)存在唯一的正根λ0.因为I(z)>0和limz→+∞I(z)=+∞,因此limz→+∞u(z)≥0,则limz→+∞u(z)=λ0.由文献[14],可得I(z)≤eλ1z,其中0<λ1<λ2是方程
d2d2
的两个特征根.将λ2代入方程(14),整理得
∫
∫∫
+∞-∞
J(y)e-λydy-d2-cλ+f(S-∞)g′(0)e-λcτ-γ=0
(14)
+∞-∞
则0<λ1<λ2<λ0.因为limz→+∞θ(z)=λ0,从而存在z∗和L满足 I(z)≥Le((λ0+λ2)/2)z∗.
这与λ1<(λ2+λ0)/2矛盾.因此,I有界.证毕.
J(y)e-λ2ydy-d2-cλ2-γ=-f(S-∞)g′(0)e-λ2cτ<0,
定理2 如果R0>1和c>c∗,那么式(5)存在行波解(S(z),I(z))满足0<S(z)<S-∞,0<I(z),S(+∞)=σ<S(-∞)=S-∞和I(±∞)=0,其中σ>0.同时
证明 由文献[14]中定理2.1,可得式(5)存在一个解(S(z),I(z))且满足0<S(z)<S-∞,0<I(z),S(+∞)=σ,S(-∞)=S-∞和I(±∞)=0.对任意的φ,υ∈R,有
因此
因此
正如文献[14]中所描述的,只要I有界,下面的结果都成立.因此在这里只给出部分证明.
□
∫f(S(z))g(I(z-cτ))dz<+∞,∫I(z)dz<+∞.
R
R
φυ
φυ
+∞-∞φυ
∫(J∗S(z)-S(z))dz=
∫∫J(y)[S(υ-y)-S(υ)]dydξ= -∫∫J(y)y∫S′(υ-ty)dtdydξ=
∫J(y)y∫[S(υ-ty)-S(φ-ty)]dtdy.
+∞-∞
10
+∞-∞
10
再对式(5)的第一个方程从υ到φ积分,可得
φυ
∫f(S(z))g(I(z-cτ))dz=
d∫(J∗S(z)-S(z))dz+cS(υ)-cS(φ)≤dC
φυ1
1
∫
φυ
(J∗S(z)-S(z))dz<2S-∞
∫
+∞-∞
J(y)|y|dy
C0.
0
+2cS-∞.
又因I(z)有界,所以存在实数M使得I(z)≤M成立.对式(5)的第二个方程从ζ到η积分,有
γ
∫f(S(z))g(I(z-cτ))dz<+∞.
R
∫
ηζ
I(z)dz=d2J(y)(-y)
R
ηζ
∫f(S(z))g(I(z-cτ))dz-c(I(η)-I(ζ)).
∫∫[I(η-θy)-I(ζ-θy)]dθdy+
10
(15)
720张 秋 陈 广 生
因此
所以
γ
2 临界波行波解的存在性
因为I(z)在R上有界,所以,由式(5)的第二个方程可得|I′(z)|<+∞且I(+∞)=0.最后,类似于文献[14],可得S(+∞)=σ,其中,σ>0.因此,定理2成立.证毕.□
本节考察当c=c∗和R0>1时,系统(5)的行波解的存在性.为了证明当c=c∗时的行波解
∫I(z)dξ
R
∫I(z)dz≤2dM∫J(y)|-y|dy+∫f(S(z))g(I(z-cτ))dz+2cM.
ηζ
2
R
R
<+∞.
的存在性,取一列严格递减序列{ck}使得limk→+∞ck=c∗.不失一般性,假设ck∈[c∗,c∗+1]且(Sk,Ik)是式(5)的行波解.
引理5 {Ik}在R上一致有界.
证明 反证法.假设不成立,则存在一个序列{zk}使得当k→+∞时,有Ik(zk)→+∞.由引理2可得,Sk(zk)→0.又因为对每一个k,均有Ik(±∞)=0和Ik(z)有界.因此,假设Ik(zk)=maxRIk(z).则I′=0.因此k(zk)
= 0=ckI′k(zk) f(Sk(zk))g(Ik(zk-cτ))-γIk(zk).
这与足够大的k矛盾,故{Ik}是一致有界的.证毕.
d2J(y)(Ik(zk-y)-Ik(zk))dy+f(Sk(zk))g(Ik(zk-cτ))-γIk(zk)≤
R
∫
□
证明 因为c∗<ck<c∗+1,所以0<Sk(z)<S-∞和0<Ik(z)<+∞.并且,Sk(-∞)=S-∞,Ik(±∞)=0和Sk(+∞)=σ<S-∞.又因为(Sk,Ik)满足
ì=d1J(y)(Sk(z-y)-Sk(z))dy-f(Sk(z))g(Ik(z-ckτ)),ïckS′k(z)
Rï
ï í=d2J(y)(Ik(z-y)-Ik(z))dy+k(z)ïckI′R
ïï f(S(z))g(I(z-cτ))-γI(z).îkkkk
对式(16)的第二个方程在R上积分,可得
则
定理3 假设R0>1和c=c∗.则方程(5)存在行波解(S∗(z),I∗(z))满足0<S∗(z)<S-∞,I∗(z)>0,S∗(-∞)=S-∞和I∗(±∞)=0,其中σ≥0.
∫∫
(16)
∫
+∞-∞
(f(Sk(z))g(Ik(z-ckτ))-γIk(z))dz=0,
γIk(z)
, ∀k∈Z.
g(Ik(z-ckτ))
因为|I′k/Ik|≤K和条件(H1),因此
γIk(z)
= inff(Sk(z))<
z∈Rg(Ik(z-ckτ))
Ik(z)Ik(z-ckτ)1
γ≤γeckKτ, ∀t∈[0,Ik(z-ckτ)].
Ik(z-ckτ)g(Ik(z-ckτ))g′(t)
z∈R
inff(Sk(z))<
一类具有非线性发生率与时滞的非局部扩散SIR模型的临界波的存在性721
又因为Ik是一致有界的,所以存在一个常数M使得Ik(z-cτ)<M.所以
由条件(H1)可得f(x)是单调递增且f(S-∞)g′(0)>γ,则存在zk满足f(Sk(zk))=γeckKτ(1/g′(M)).经过平移,可得f(Sk(0))=γeckKτ(1/g′(M)).又因为0<Sk<S-∞且{Ik}Ik→I∗在C1loc(R)上成立且满足
以及
}有界.因此,根据Arzela⁃Ascoli定理,当k→+∞时,有Sk→S∗和是一致有界的,所以{I′k
ì=d1J(y)(S∗(z-y)-S∗(z))dy-f(S∗(z))g(I∗(z-cτ)),ïc∗S′∗(z)
Rï
ï
ícI′(z)=dJ(y)(I(z-y)-I(z))dy+
2∗∗ï∗∗R
ïï f(S(z))g(I(z-cτ))-γI(z),î∗∗∗
∗Kτ
z∈R
inff(Sk(z))<γeckKτ
1
.
g′(M)
∫∫
(17)
类似于c>c∗,可得I∗(±∞)=0.
f(S∗(0))=γec
1
, 0≤S∗≤S-∞,0≤I∗≤+∞.g′(M)
因此,式(17)的第一个方程为J(y)S∗(y0-y)dy=0.所以,S∗≡0且f(S∗)≡0.这与f(Sk(0))矛盾.因此,S∗>0.
R
=0.首先,证明S∗>0.反证法,假设存在常数y0∈R使得S∗(y0)=0成立,则S′∗(y0)
∫
明,可得I∗≡0.则式(17)的第一个方程为
因为
记
H(v)
d1
其次,证明I∗>0.反证法,假设存在y∈R使得I∗(y)=0成立.类似于S∗(z)>0的证
=d1J(y)(S∗(z-y)-S∗(z))dy.c∗S′∗(z)
R
(∫J(y)e
R
v=0
∫
(18)
-vy
dy-1-c∗v.
∗
)
因此,H(v)=0有两个实根0和v0.由式(18)、文献[19]中性质3.6和S∗(z)在R上有界,可得S∗(z)≡C,其中C∈[0,S-∞].因为f(S∗(0))=γecγec
∗Kτ
∗Kτ
∂H
H(0)=0,
∂v
=-c
∂2H
<0,2=d1J(y)y2e-vydy>0,limH(v)=+∞.
v→+∞R∂v
∫
(1/g′(M)).再者,因为f(S-∞)g′(M)>γec
f(Sk(z))>f(Sk(ρk))=
f(S-∞)+γec
∗Kτ
∗Kτ
>γ,所以,选择ρk<0使得
(1/g′(M)),因此f(C)≡
成立.
f(S-∞)+γec
∗Kτ
2
(1/g′(M))
2
(1/g′(M))
, ∀ξ<ρk,
义两个函数;
接下来,证明存在一个正常数T0使得ρk≥-T0.反证法,当k→+∞时,有ρk→+∞.定
722张 秋 陈 广 生
Ik(ρk)
则对每一个k∈N,有(Uk(z),Vk(z))满足
Uk(z)=Sk(ρk+z),Vk(z)=
Ik(ρk+z)
, ∀z∈R.
因为
ì=d1J(y)(Uk(z-y)-Uk(z))dy-f(Uk(z))g(Ik(ρk+z-cτ)),ïckU′k(z)
Rï
ïïcV′(z)=dJ(y)(V(z-y)-V(z))dy-
2kk(19)íkk
R
ï
g(Ik(z+ρk-cτ))ï
-γVk(z).ï f(Uk(z))
Ik(ρk)î
Ik(ρk)
=
g(Ik(z+ρk-cτ))Ik(z+ρk-cτ)
Uk(z-cτ)=g′(t1)Uk(z-cτ),
∫∫
g(Ik(z+ρk-cτ))
则式(19)的第二个方程为
又因为
∀t1∈[0,Ik(z+ρk-cτ)],
=d2J(y)(Vk(z-y)-Vk(z))dy-f(Uk(z))g′(t1)Vk(z-cτ)-γVk(z).ckV′k(z)
R
æz+ρkI′ök(x)= expçdx÷,Ik(ρk)èρkIk(x)ø
则Vk(z)在R上局部一致有界.又因为对所有的k∈N,均有Uk(z)<S-∞.因此,当k→+∞
Vk(z)=
Ik(z+ρk)
∫
∫
1
时,有Uk(z)→U+∞(z)和Vk(z)→V+∞(z),其中U+∞(z)∈C1loc(R),V+∞(z)∈Cloc(R).因为limk→+∞Ik(ρk)=0,所以(U+∞(z),V+∞(z))满足
c∗U′+∞(z)=d1J(y)(U+∞(z-y)-U+∞(z))dy,
R
C1是任意常数.因为
所以
f(U+∞(z))g′(0)V+∞(z-cτ)-γV+∞(z).(21)类似地,由文献[19]中性质3.6、式(19)和U+∞(z)有界,可得U+∞(z)≡C1成立,其中
f(Uk(0))=f(Sk(ρk))=f(U+∞(z))=
f(S-∞)+ec
∗Kτ
∗Kτ
c∗V′+∞
∫(z)=d∫J(y)(V
2R
(20)
+∞
(z-y)-V+∞(z))dy+
2
(γ/g′(M))
.
2
将f(U+∞(z))代入到方程(21),整理得
f(S-∞)+ec
f(S-∞)+ec(γ/g′(M))
.
c∗V′+∞(z)=d2J(y)(V+∞(z-y)-V+∞(z))dy+
R
∗Kτ
令
∫
2
(γ/g′(M))
g′(0)V+∞(z-c∗τ)-γV+∞(z).
Q(χ)d2
(∫J(y)e
R
-χy
dy-1-c∗χ+
)
一类具有非线性发生率与时滞的非局部扩散SIR模型的临界波的存在性723
因为Sk(-∞)=S-∞和引理4,因此,由式(5)的第二个方程可得
I′k(z)
ψklim
z→-∞I(z)k存在且是方程
d2J(y)e-ψydy-c∗ψ+f(S-∞)g′(0)e-c
R
f(S-∞)+ec
∗Kτ
2
(γ/g′(M))
g′(0)e-
χc∗τ
-γ.
的一个特征根{ψ±},则当k→+∞时,ψ±→ψ∗且满足
d2J(y)e-ψydy-c∗ψ+f(S-∞)g′(0)e-c
R
∫∫
∗ψτ
-γ=0
∗ψτ
-γ=0.(22)
因此
f(S-∞)g′(0)e-ψcτ-f(S-∞)g′(0)e-ψcτ=0.
所以Q(χ)=0有两个根χ1,χ2,且0<χ1<χ2.类似地,再由文献[19]中性质3.6和式(21),
∗∗
∗∗
γ-(f(S-∞)g′(0)e-ψ
Q(ψ)=
∗
f(S-∞)+ec
∗Kτ
2
(γ/g′(M))
∗c∗τ
g′(0)e-ψ
∗c∗τ
-
-γ)≤
可得
V+∞(z)=T3e1+T4e2,
其中Tj(j=3,4)是常数.
另一方面,当k→+∞时,有
χ
χ
(23)
Vk(z)=
∗z
Ik(ρk+z)Ik(ρk)
k
=e∫ρk
ρ+z
(I′k(s)/Ik(s))ds
则V+∞(z)=eψ成立.因此
.因为χ1<ψ∗<χ2,这与式(22)矛盾.所以,存在一个T0>0使得ρk≥-T0
f(S-∞)+ec
⇒
2∗
f(S-∞)+ecKτγ(1/g′(M))
2
∗Kτ
→e
ψ∗z
,
f(Sk(-T0))>
γ(1/g′(M))
这与f(Sk(z))≡ec=0且S∗()满足
f(S∗(-T0))≥
∗Kτ
.
使得S∗()=S-∞.从而,S′)∗(
-cτ))≤
接下来,可证S∗(z)<S-∞在R上成立.假设存在实数 0=c∗S′)=d1J(y)(S∗(∗(
-f(S∗())g(I∗(
R
γ(1/g′(M))矛盾.因此,可得I∗(ξ)>0.
∫
这与I∗>0矛盾.
-cτ)).
-y)-S∗())dy-f(S∗())g(I∗(
最后,证明S∗(-∞)=S-∞.尽管证明的思想与文献[17]中给定的相似,但为了文章的
完整性还是给出详细的证明.设=liminfx→-∞S∗(x)=S-∞.反证法,假设<S-∞,则存在点
+∞
列{ηn}n成立.定义函数:=1使得limn→+∞ηn=-∞和limn→+∞S∗(ηn)=
ϕn(x)S∗(x+θn),ψn(x)I∗(x+θn), x∈R.
724张 秋 陈 广 生
因此,当n→+∞时,ψn(x)→0.提取子序列(仍然用{θn}表示),当n→+∞时,limn→+∞θn=-∞且ϕn→ϕ+∞,其中ϕ+∞∈C1
loc(R).因此ϕ+∞满足所以,ϕ+∞≡.因此
limS∗(x+θn)=.
n→+∞
c∗ϕ′+∞(z)=d1J(y)(ϕ+∞(z-y)-ϕ+∞(z))dy.
R
∫
又因为(ck,Sk,Ik)是系统(16)的行波解,所以
(24)(25)
对式(19)的第一个方程从-∞到θn积分,可得
=d1J(y)(Sk(z-y)-Sk(z))dy-f(Sk(z))g(Ik(z-cτ)).ckS′k(z)
R
∫
cθn
k[S∫
k(θn)-σ]=d1
-
θ∞
R
k
z-y)-Sk
(z))dydz-
n
-∞
f(Sk(z))g(Ik(∫z-∫J(y)(S(cτ))dz,
其中n∈N∫.又因为θn-∞
∫SkR
J(y)(S(k
(-z∞-)y)=S-∞且J(x)=J(-x),可得
-Sk
θ(z))dydz=
n
R
-∞
k
k
∫J(y)∫(S(θz-y)-S(z))dzdy= n
1 ∫∫R
J(y)(-y)J(y)(-y)∫∫S′(z-ty)dtdzdy=
1∫1-∞
[S′0
k
R
0k
(θn
-ty)-S-∞
]dtdy= 因此,对k取极限∫R
J(y)y,有
∫0
Sk
(θn
+ty)dtdy.
c∗
[S∗(θn)-S-∞]=d1J(y)y1∫
R
S∗θ(θn+ty)dtdy-
n
∫
0
-∞
f(S∗因为S∗,I∗均有界且∫
(z))g(I∗(z∫
-cτ))dz.
f(S∗(z))g(因此 I∗(z-cτ))dz<+∞.则对n取极限,有,S 3 应∗(0->c∗
[-R
∞)=S-S-∞]≥0.∞.证毕.
这一节主要通过具体的生物模型进一步验证临界波的存在性用实例
.在系统(5)中,令f(S)=S,g(I)=βI/(1+αI),则
f(S(x,t))g(I(x,t-τ))=βS(x,t)I(x,t-τ)
那么,系统(5)就变为如下具有非线性发生率的时滞反应扩散方程1+αI(x,t-τ)
,
ìï ï
cS′(z)=dí
1(J∗S(z)-S(z))-βS(z)I(z-cτ):
ï1+αI(,ïî
cI′(z)=d2(J∗I(z)-I(z))+βS(z)1+I(αIz(-zcτz--)cτ)cτ)-γI(z),□
(26)
一类具有非线性发生率与时滞的非局部扩散SIR模型的临界波的存在性725
其中β>0表示感染率,γ表示恢复率,对于系统(26),条件(H1)和(H2)显然成立.边界条件:(S(-∞),I(-∞))=(S-∞,0)是初始无病平衡点,(S(+∞),I(+∞))=(S+∞,0)是某一个无病平衡点,且S(+∞)为任意实数.显然,系统(26)是在文献[20]的基础上考虑到霍乱具有潜伏期且易感者与感染者的活动可能是非局部的因素下建立的模型.
对系统(26)的第二个方程在初始无病平衡点线性化,则特征方程为
cλ=d2
显然,给定βS-∞>γ,则系统(26)的最小传播速度为
c∗
infd2
∫
+∞-∞
J(y)e-λydy-d2+βS-∞e-λcτ-γ.
+∞-∞
λ>0
∫
J(y)e-λydy-d2+βS-∞e-λcτ-γ
λ
.
因为g′(I)=β/(1+αI)2,所以g′(I)>0且g′(0)=β.又S(z)≤S-∞,因此f(S(z))≤S-∞,则由系统的(26)的第二个方程可得
d2I′≤Ic
所以引理1成立.由此可知引理2、引理3和引理4都成立,则I(z)是有界的,所以定理1成立.因此,类似于文献[14]的方法,易得当波速c>c∗时,系统(26)存在满足边界条件的行波解,所以定理2成立.
存在满足边界条件的行波解.不失一般性,假设ck∈[c∗,c∗+1]且(Sk,Ik)是系统(26)的行波解.则(Sk,Ik)满足下面的模型:
其中
Ik(z-ckτ)ìïcS′(z)=d(J∗S(z)-S(z))-βS(z),kk1kkkï1+αIk(z-ckτ)ïí
Ik(z-ckτ)ï
=-+-γIk(z),cI′(z)d(J∗I(z)I(z))βS(z)2kkkïkk
+-1αI(zcτ)îkk
任意选取一列严格递减序列{ck}使得limk→+∞ck=c∗.由于当波速c>c∗时,系统(26)
∫
d2+γexp(τ(d2+γ))I(z-y)
++J(y)dyβS-∞,
I(z)ccR
(27)
成立,因而{Ik}是一致有界的,即总可以找到一个M>0使得Ik≤M.又因为Ik有界,则由系统(27)的第二个方程可得,|I′k/Ik|有界,则存在一个常数K>0使得|I′k/Ik|≤K.的行波解,即定理3对于模型(26)依然成立.
因此由理论结果可知,当R0>(1+αM)2ec
∗Kτ
又因为引理2和定理1成立,所以对每一个k,都有Ik有界且Ik(±∞)=0,则证得引理5
(Sk(-∞),Ik(-∞))=(S-∞,0),(Sk,Ik)(+∞)=(S+∞,0).
且c=c∗时,系统(26)存在满足边界条件
4 总 结
本文研究了一类具有非线性发生率与时滞的非局部扩散SIR模型的临界波的存在性,文
∗
献[14⁃15]中没有讨论当c=c∗时,系统(5)的行波解是否存在.受到文献[17⁃18]的启发,本文证明了系统(5)在R0>g′(0)ecKτ/g′(M)和c=c∗时,存在行波解.而对于1<R0≤g′(0)ecKτ/g′(M)和c=c∗,同样认为其是临界波存在的条件,这将是笔者下一步需要研究
∗
的问题.
726张 秋 陈 广 生
参考文献(References):
[1] [2] [3] [4] [5] [6] [7] [8] [9]
VANDENDRIESSCHEP,WATMOUGHJ.Reproductionnumbersandsub⁃thresholdendemic2002,180(1/2):29⁃48.1969⁃1993.
equilibriaforcompartmentalmodelsofdiseasetransmission[J].MathematicalBiosciences,YANGFY,LIY,LIWT,etal.TravelingwavesinanonlocaldispersalKermack⁃McKendrick
epidemicmodel[J].DiscreteandContinuousDynamicalSystems(SeriesB),2013,18(7):
WANGHY.Spreadingspeedsandtravelingwavesfornon⁃cooperativereaction⁃diffusionsys⁃tems[J].JournalofNonlinearScience,2011,21(5):747⁃783.
WANGHY,WANGXS.TravelingwavephenomenainaKermack⁃MckendrickSIRmodel[J].JournalofDynamics&DifferentialEquations,2016,28(1):143⁃166.
WANGJB,LIWT,YANGFY.TravelingwavesinanonlocaldispersalSIRmodelwithnon⁃tions,2015,27(1/3):136⁃152.
FUSC,GUOJS,WUCC.Travelingwavesolutionsforadiscretediffusiveepidemicmodel[J].JournalofNonlinearandConvexAnalysis,2016,17(9):1739⁃1751.
ZHANGTR,WANGWD,WANGKF.Minimalwavespeedforaclassofnon⁃cooperativedif⁃fusion⁃reactionsystem[J].JournalofDifferentialEquations,2016,260(3):2763⁃2791.WANGXS,WANGHY,WUJH.Travelingwavesofdiffusivepredator⁃preysystems:diseaseoutbreakpropagation[J].Discrete&ContinuousDynamicalSystems,2017,32(9):3303⁃3324.
LIY,LIWT,YANGFY.TravelingwavesforanonlocaldispersalSIRmodelwithdelayandexternalsupplies[J].AppliedMathematicsandComputation,2014,247:723⁃740.nicationsonPure&AppliedAnalysis,2015,14(3):1001⁃1022.
localdelayedtransmission[J].CommunicationsinNonlinearScience&NumericalSimula⁃
[10] LIY,LIWT,LING.TravelingwavesofadelayeddiffusiveSIRepidemicmodel[J].Commu⁃[11] WANGZC,WUJH.TravellingwavesofadiffusiveKermack⁃McKendrickepidemicmodel
Sciences,2010,466(2113):237⁃261.
[12] LIWT,YANGFY,MAC,etal.TravelingwavesolutionsofanonlocaldelayedSIRmodel
2014,19(2):467⁃484.
withoutoutbreakthreshold[J].DiscreteandContinuousDynamicalSystems(SeriesB),
withnon⁃localdelayedtransmission[J].ProceedingsMathematicalPhysical&Engineering
[13] BAIZG,WUSL.TravelingwavesinadelayedSIRepidemicmodelwithnonlinearincidence[14] 邹霞,吴事良.一类具有非线性发生率与时滞的非局部扩散SIR模型的行波解[J].数学物理学
2018,38(3):496⁃513.(inChinese))011513.DOI:10.1063/1.5021761.
[J].AppliedMathematicsandComputation,2015,263:221⁃232.
报,2018,38(3):496⁃513.(ZOUXia,WUShiliang.TravelingwavesinanonlocaldispersalSIRepidemicmodelwithdelayandnonlinearincidence[J].ActaMathematicaScientia,
[15] ZHANGSP,YANGYR,ZHOUYH.TravelingwavesinadelayedSIRmodelwithnonlocal
dispersalandnonlinearincidence[J].JournalofMathematicalPhysics,2018,59(1):
[16] CHENYY,GUOJS,HAMELF.Travelingwavesforalatticedynamicalsystemarisingina
diffusiveendemicmodel[J].Nonlinearity,2016,30(6).DOI:10.1088/1361⁃6544/aa6b0a.
[17] YANGFY,LIWT.TravelingwavesinanonlocaldispersalSIRmodelwithcriticalwavespeed
一类具有非线性发生率与时滞的非局部扩散SIR模型的临界波的存在性
[J].JournalofMathematicalAnalysis&Applications,2017,458(2):1131⁃1146.model[J].JournalofDifferentialEquations,2017,262(1):272⁃282.
727
[18] WUCC.Existenceoftravelingwaveswiththecriticalspeedforadiscretediffusiveepidemic[19] ZHANGGB,LIWT,WANGZC.Spreadingspeedsandtravelingwavesfornonlocaldispersal
2012,252(9):5096⁃5124.
equationswithdegeneratemonostablenonlinearity[J].JournalofDifferentialEquations,
[20] CAPASSOV,SERIOG.AgeneralizationoftheKermack⁃McKendrickdeterministicepidemic
model[J].MathematicalBiosciences,1978,42(1/2):43⁃61.
ExistenceofCriticalTravelingWavesforNonlocal
DispersalSIRModelsWithDelayand
NonlinearIncidence
(1.SchoolofMathematicsandStatistics,XidianUniversity,
Xi’an710071,P.R.China;
2.CollegeofMathematicsandComputerScience,GuangxiScience&
ZHANGQiu1, CHENGuangsheng1,2
TechnologyNormalUniversity,Laibin,Guangxi546199,P.R.China)
Abstract:TheexistenceoftravelingwavesolutionsfornonlocaldispersalSIRepidemicmodelswithdelaywasstudied.Firstly,theboundednessofIwasprovedbycontradiction.Thenac⁃cordingtotheboundednessofI,theexistenceoftravelingwaveswithc>c∗wasestablished.Secondly,throughfurtheranalysisoftravelingwaveswithsuper⁃criticalspeeds,theexistenceoftravelingwaveswiththecriticalspeedwasderived.Finally,theinfluenceofbasicreproduc⁃tionnumberR0ontheexistenceofc>c∗wasdiscussed.
tionnumber
Foundationitem:TheNationalNaturalScienceFoundationofChina(GeneralProgram)
(11671315)
Keywords:travelingwavesolution;criticalwavespeed;nonlocaldispersal;basicreproduc⁃
张秋,陈广生.一类具有非线性发生率与时滞的非局部扩散SIR模型的临界波的存在性[J].应用数学
和力学,2019,40(7):713⁃727.
ZHANGQiu,CHENGuangsheng.ExistenceofcriticaltravelingwavesfornonlocaldispersalSIR713⁃727.
modelswithdelayandnonlinearincidence[J].AppliedMathematicsandMechanics,2019,40(7):
引用本文/Citethispaper:
因篇幅问题不能全部显示,请点此查看更多更全内容