设为首页
加入收藏
期刊导航
网站地图
首页
期刊
数学与物理
地球与环境
信息通讯
经济与管理
生命科学
工程技术
医药卫生
人文社科
化学与材料
会议
合作
新闻
我们
招聘
千人智库
我要投稿
办刊
期刊菜单
●领域
●编委
●投稿须知
●最新文章
●检索
●投稿
文章导航
●Abstract
●Full-Text PDF
●Full-Text HTML
●Full-Text ePUB
●Linked References
●How to Cite this Article
AdvancesinAppliedMathematics
应用数学进展
,2021,10(6),2083-2094
PublishedOnlineJune2021inHans.http://www.hanspub.org/journal/aam
https://doi.org/10.12677/aam.2021.106218
具有潜伏期和分布时滞的传染病模型及动力学
分析
忽丹丹,杨雨
中国地质大学数学与物理学院,湖北武汉
收稿日期:
2021
年
5
月
21
日;录用日期:
2021
年
6
月
9
日;发布日期:
2021
年
6
月
24
日
摘要
本文研究了一类潜伏期具有感染性的疾病传播模型
,
其中引入分布时滞刻画了潜伏期患者转症为感
染者的概率
.
通过分析相应的特征方程
,
讨论了无病平衡点和地方病平衡点的局部稳定性
.
再通过
构造合适的
Lyapunov
泛函
,
应用
LaSalle
不变性原理
,
得到当疾病传播阈值即基本再生数小于
1
时
,
无病平衡点是全局渐近稳定的
,
疾病将会消失
;
当基本再生数大于
1
时
,
地方病平衡点是全
局渐近稳定的
,
疾病将持续存在
.
最后用数值模拟验证了理论结果。
关键词
SEIR
模型,分布时滞,基本再生数,
Lyapunov
泛函,全局稳定性
DynamicAnalysisonLatentPeriodand
DistributedDelayEpidemiologicalModel
DandanHu,YuYang
SchoolofMathematicandPhysics,ChinaUniversityofGeosciences,WuhanHubei
Received:May21
st
,2021;accepted:Jun.9
th
,2021;published:Jun.24
th
,2021
文章引用
:
忽丹丹
,
杨雨
.
具有潜伏期和分布时滞的传染病模型及动力学分析
[J].
应用数学进展
,2021,10(6):
2083-2094.DOI:10.12677/aam.2021.106218
忽丹丹,杨雨
Abstract
Inthispaper,amodelofdiseasetransmissionwithinfectiousincubationperiodis
studied,in which adistributed timedelayis introducedto characterize theprobability
ofpatientsinlatencybecominginfected.Byanalyzingthecorrespondingcharacter-
isticequations,thelocalstabilityofdisease-freeequilibriumandendemicequilibrium
isdiscussed.Then,byconstructinganappropriateLyapunovfunctionalandapplying
LaSalleinvarianceprinciple,weobtainedthatthedisease-freeequilibriumpointis
globallyasymptoticallystableandthediseasewilldisappearwhenthediseasetrans-
missionthreshold,namelythebasicregenerationnumber,islessthan1.Whenthe
basicreproductionnumberisgreaterthan 1,theendemicequilibriumpoint isglobally
asymptoticallystableandthediseasewillpersist.Finally,thetheoreticalresultsare
verifiedbynumericalsimulation.
Keywords
SEIRModel,DistributedDelay,BasicReproductionNumber,LyapunovFunctional,
GlobalStability
Copyright© 2021byauthor(s)andHansPublishersInc.
This work is licensed under theCreative Commons Attribution InternationalLicense (CCBY4.0).
http://creativecommons.org/licenses/by/4.0/
1.
引言
传染病的传播阻碍了社会经济的发展
,
给人类的生存带来了长期而严峻的威胁
,
对传染病发病
机理、流行规律和防治策略的研究
,
其重要性日益突出
,
且已成为当今世界需要迫切解决的一个重
大问题
[1].
自
1760
年荷兰的数学家
,
医学家
DanielBernoulli
开始构建数学模型研究牛痘接种对
天花的预防和治疗以来
,
人们就开始从数学上对传染病进行研究
,
确定性传染病模型的研究应该说
始于
20
世纪初
.1906
年
,Hamer
为了研究荨麻疹的反复流行
,
构造并分析了一类离散的时间模
型
[2].1911
年
,
英国微生物学家
RonaldRoss
博士利用连续的微分方程模型对疟疾在蚊子与人群
之间传播的动力学行为进行了研究
.
自
1927
年美国数学家
Kermack
和苏格兰医学家、流行病学家
McKendrick
构造了著名的
SIR
“仓室”模型以来
,
数学模型成为研究疾病传播、评估感染风险、
DOI:10.12677/aam.2021.1062182084
应用数学进展
忽丹丹,杨雨
优化控制策略的重要工具
[3].
2020
年
,
新型冠状病毒肺炎(
CoronaVirusDisease2019,COVID-19
)疫情在全球多个国家暴
发,随后呈现出大规模蔓延趋势
[4].
据世界卫生组织实时统计数据显示,截至
2021
年
4
月
22
日,
全球确诊
COVID-19
人数超
1
亿
4439
万例,累计死亡病例超
307
万例,目前全球已有
24
个国家
COVID-19
确诊病例超过
100
万例
[5].
此类病毒存在潜伏期,属于聚集性疫情,且存在无症状感
染者
[6].
近年来
,
一些研究者也提出了有关病毒在潜伏期具有感染性的传染病模型
[7–10].
在此模
型的基础上并结合新冠病毒潜伏期具有感染力的特征
,
我们考虑了一个潜伏期具有感染性的
SEIR
流行病模型
.
该模型包含对
COVID-19
疾病的易感群体
(S),
潜伏期群体
(E),
感染群体
(I)
和康复
群体
(R).
这里
,
被感染的群体已经表现疾病症状
,
并可以传染给其他人
.
暴露类个体处于潜伏期
,
虽
然他们没有表现出发病症状
,
但仍然有能力感染他人
.
为了研究新型冠状病毒在潜伏期对疾病传播动态的影响
,
同时考虑到潜伏期群体转症到有症状
感染群体需要一定的时间
,
一般为
1-14d,
多为
3-7d[11].
我们假设这个延迟在
[0
,τ
]
的区间内为
一个分布式参数
[12],
其中
τ
是易感者从暴露于病毒到出现发病症状之间的延迟最大值
.
p
(
s
)
是关
于时滞
s
的非负连续分布函数
,
满足
τ
0
p
(
s
)
d
s
=1
.
在此,我们考虑了以下由分布时滞微分方程描
述的
SEIR
流行病模型
d
S
(
t
)
d
t
=Λ
−
β
1
SE
−
β
2
SI
−
µS,
d
E
(
t
)
d
t
=
β
1
SE
+
β
2
SI
−
(
α
+
µ
)
E,
d
I
(
t
)
d
t
=
τ
0
αE
(
t
−
s
)
p
(
s
)
d
s
−
(
γ
+
µ
)
I,
d
R
(
t
)
d
t
=
γI
−
µR.
(1.1)
在模型
(1.1)
中
,
Λ
为新增易感者的数量
,
µ
是自然死亡率
,
γ
为恢复率
.
α
是常数
,
在这个常数
速率下
,
具有传染性的无症状感染者开始出现发病症状
.
标准发生率的模型更符合实际也最常用
,
因
此我们将无症状感染者和感染者的疾病传播率分别表示为
β
1
和
β
2
.
模型
(1.1)
中所有参数都是非负
常数
.
由于
R
(
t
)
没有出现在系统
(1.1)
的前三个方程中,故可将系统
(1.1)
简化为以下子系统
d
S
(
t
)
d
t
=Λ
−
β
1
SE
−
β
2
SI
−
µS,
d
E
(
t
)
d
t
=
β
1
SE
+
β
2
SI
−
(
α
+
µ
)
E,
d
I
(
t
)
d
t
=
τ
0
αE
(
t
−
s
)
p
(
s
)
d
s
−
(
γ
+
µ
)
I.
(1.2)
系统
(1.1)
满足如下初始条件
S
(
θ
)=
ϕ
1
(
θ
)
,E
(
θ
)=
ϕ
2
(
θ
)
,I
(
θ
)=
ϕ
3
(
θ
)
,R
(
θ
)=
ϕ
4
(
θ
)
,ϕ
i
(
θ
)
>
0
,θ
∈
[
−
τ,
0]
,
(1.3)
其中
Φ=(
ϕ
1
(
θ
)
,ϕ
2
(
θ
)
,ϕ
3
(
θ
)
,ϕ
4
(
θ
))
∈
C
[
−
τ,
0]
,
R
4
+0
,
C
表示从区间
[
−
τ,
0]
到
R
4
+0
的全体连续
函数所构成的集合
,
对一致收敛拓扑构成
Banach
空间
.
DOI:10.12677/aam.2021.1062182085
应用数学进展
忽丹丹,杨雨
2.
解的非负性和有界性
引理
2.1
(1.1)
(1.3)
,
S
(
t
)
,E
(
t
)
,I
(
t
)
,R
(
t
)
,
.
证明
利用参数变分法可得系统
(1.1)
有如下的解
S
(
t
)=
S
(0)
e
−
∫
t
0
(
β
1
E
(
ξ
)+
β
2
I
(
ξ
)+
µ
)
d
ξ
+Λ
t
0
e
−
∫
t
η
(
β
1
E
(
ξ
)+
β
2
I
(
ξ
)+
µ
)
d
ξ
d
η,
E
(
t
)=
E
(0)
e
−
∫
t
0
(
α
+
µ
−
β
1
S
(
ξ
)
d
ξ
+
β
2
SI
t
0
e
−
∫
t
η
(
α
+
µ
−
β
1
S
(
ξ
))
d
s
d
η,
I
(
t
)=
I
(0)
e
−
(
γ
+
µ
)
t
+
α
τ
0
E
(
t
−
s
)
p
(
s
)
d
s
t
0
e
−
(
γ
+
µ
)(
t
−
η
)
d
η,
R
(
t
)=
R
(0)
e
−
µt
+
γI
t
0
e
−
µ
(
t
−
η
)
d
η.
这意味着该系统的解
S
(
t
)
,E
(
t
)
,I
(
t
)
,R
(
t
)
是非负的
.
接下来证明解的有界性,我们定义
D
(
t
)=
τ
0
p
(
s
)
S
(
t
−
s
)
d
s
+
τ
0
p
(
s
)
E
(
t
−
s
)
d
s
+
I
(
t
)+
R
(
t
)
.
(2.1)
将
D
(
t
)
沿着系统
(1.1)
的解的方向关于时间
t
求导
,
可得
d
D
(
t
)
d
t
=
τ
0
p
(
s
)Λ
−
β
1
S
(
t
−
s
)
E
(
t
−
s
)
−
β
2
S
(
t
−
s
)
I
(
t
−
s
)
−
µS
(
t
−
s
)
d
s
+
τ
0
p
(
s
)
β
1
S
(
t
−
s
)
E
(
t
−
s
)+
β
2
S
(
t
−
s
)
I
(
t
−
s
)
−
(
α
+
µ
)
E
(
t
−
s
)
d
s
+
τ
0
p
(
s
)
E
(
t
−
s
)
d
s
−
µI
(
t
)
−
µR
(
t
)
=Λ
−
µ
τ
0
p
(
s
)
S
(
t
−
s
)
d
s
−
µ
τ
0
p
(
s
)
E
(
t
−
s
)
d
s
−
µI
(
t
)
−
µR
(
t
)
=Λ
−
µD
(
t
)
<
0
,D
(
t
)
>
Λ
µ
,
>
0
,D
(
t
)
<
Λ
µ
.
这意味着
D
(
t
)
是有界的,因此系统
(1.1)
所有的解
S
(
t
)
,E
(
t
)
,I
(
t
)
,R
(
t
)
都是有界的
,
得证
.
3.
基本再生数与平衡点的存在性
确定性流行病模型的一个主要关注点是找到一个引入地区的疾病能够发展成大规模暴发的条
件
,
如果疾病大规模暴发,该疾病就可能成为流行病
.
在传染病动力学中
,
基本再生数是一个很重要
的概念
.
它表示在发病初期
,
当所有人均为易感者时
,
一个病人在其平均患病期内所传染的人数
.
在
这方面,确定性模型的一个有用的阈值被称为基本再生数
ℜ
0
.
接下来我们将给出模型
(1.2)
的基本
再生数表达式
[13].
DOI:10.12677/aam.2021.1062182086
应用数学进展
忽丹丹,杨雨
令
x
=(
E,I
)
T
,
模型
(1.2)
可表示为
d
x
d
t
=
F
(
t,x
(
t
))
−
V
(
t,x
(
t
))
,
(3.1)
此处
F
=
F
1
F
2
=
β
1
SE
+
β
2
SI
0
,
V
=
V
1
V
2
=
(
α
+
µ
)
E
τ
0
αE
(
t
−
s
)
p
(
s
)
d
s
+(
γ
+
µ
)
I
.
接下来,我们设定
F
(
t
)
和
V
(
t
)
是
2×2
的矩阵,定义为
F
(
t
)=(
∂
F
i
(
t,x
(
t
))
∂x
j
)
1
≤
i,j
≤
2
与
V
(
t
)=
(
∂
V
i
(
t,x
(
t
))
∂x
j
)
1
≤
i,j
≤
2
.
然后,由
(3.1)
可以得出
F
(
E
0
)=
β
1
S
0
β
2
S
0
00
,V
(
E
0
)=
α
+
µ
0
−
αγ
+
µ
.
模型
(1.2)
的再生矩阵为
FV
−
1
=
β
1
S
0
α
+
µ
+
αβ
2
S
0
(
α
+
µ
)(
γ
+
µ
)
β
2
S
0
γ
+
µ
00
,
以及其再生矩阵的谱半径
ρ
FV
−
1
=
β
1
Λ(
γ
+
µ
)+
β
2
α
Λ
(
γ
+
µ
)(
α
+
µ
)
µ
.
对于模型
(1.2),
从以上方程得到基本再生数为
ℜ
0
=
β
1
Λ(
γ
+
µ
)+
β
2
α
Λ
(
γ
+
µ
)(
α
+
µ
)
µ
=
Λ
µ
(
α
+
µ
)
(
β
1
+
αβ
2
γ
+
µ
)
.
(3.2)
该流行病模型
(1.2)
有两个稳态点
,
一个是无病平衡点
(DFE)
E
0
=(
S
0
,E
0
,I
0
)=(Λ/
µ,
0
,
0)
.
另一个是地方性平衡点
E
∗
=(
S
∗
,E
∗
,I
∗
)
,
它是以下系统的正解
Λ
−
β
1
S
∗
E
∗
−
β
2
S
∗
I
∗
−
µS
∗
=0
,
β
1
S
∗
E
∗
+
β
2
S
∗
I
∗
−
(
α
+
µ
)
E
∗
=0
,
α
τ
0
E
∗
p
(
s
)
d
s
−
(
γ
+
µ
)
I
∗
=0
.
(3.3)
由
(1.1),(3.2)
和
(3.3),
我们可以得到
S
∗
=
Λ
µ
ℜ
0
,E
∗
=
µ
(
γ
+
µ
)
β
1
(
γ
+
µ
)+
β
2
α
(
ℜ
0
−
1)
,I
∗
=
αµ
β
1
(
γ
+
µ
)+
β
2
α
(
ℜ
0
−
1)
.
显然,模型
(1.2)
存在一个正的地方病平衡点
E
∗
=(
S
∗
,E
∗
,I
∗
)
,
当且仅当
ℜ
0
>
1
.
接下来我们
将分析以上两个平衡点的稳定性
.
DOI:10.12677/aam.2021.1062182087
应用数学进展
忽丹丹,杨雨
4.
动力学分析
定理
4.1
1)
当
ℜ
0
<
1
时,系统
(1.2)
的无病平衡点
E
0
是局部渐近稳定的
;
2)
当
ℜ
0
>
1
时,系统
(1.2)
的地方性平衡点
E
∗
是局部渐近稳定的
.
证明
1)
系统
(1.2)
在
DFE
处的
Jacobi
矩阵形式如下
J
(
S,E,I
)=
−
µ
−
β
1
S
0
−
β
2
S
0
0
β
1
S
0
−
(
α
+
µ
)
β
2
S
0
0
α
τ
0
p
(
s
)
e
−
λs
d
s
−
(
γ
+
µ
)
.
(4.1)
计算
E
0
处的特征方程,得到
(
λ
+
µ
)
λ
2
+
Aλ
+
B
−
αβ
2
S
0
τ
0
p
(
s
)
e
−
λs
d
s
=0
.
(4.2)
这里
A
=
γ
+
µ
+
αβ
2
S
0
γ
+
µ
−
(
ℜ
0
−
1)(
α
+
µ
)
,
B
=
αβ
2
S
0
−
(
ℜ
0
−
1)(
α
+
µ
)(
γ
+
µ
)
.
显然
,
λ
0
=
−
µ
是方程
(4.2)
的根
.
研究方程
(4.2)
根的性态我们只需要考虑方程
λ
2
+
Aλ
+
B
−
αβ
2
S
0
τ
0
p
(
s
)
e
−
λs
d
s
=0
.
(4.3)
当
τ
=0
时,由
A
与
B
都是正数,我们知道方程
(4.3)
的根都具有负实部且
λ
=0
不是方
程
(4.3)
的根
.
当
τ>
0
时,我们考虑如果方程
(4.3)
有成对的纯虚根存在
,
不失一般性
,
我们假设
λ
=
±
iω
(
ω>
0
)
是它的根当且仅当
ω
满足
−
ω
2
+
γ
+
µ
+
αβ
2
S
0
γ
+
µ
−
(
ℜ
0
−
1)(
α
+
µ
)
iω
+
αβ
2
S
0
−
(
ℜ
0
−
1)(
α
+
µ
)(
γ
+
µ
)
−
αβ
2
S
0
τ
0
p
(
s
)(
cos
(
ωs
)
−
isin
(
ωs
))
d
s
=0
.
(4.4)
把实部和虚部分开
,
可以得到
αβ
2
S
0
−
(
ℜ
0
−
1)(
α
+
µ
)(
γ
+
µ
)
−
ω
2
=
αβ
2
S
0
τ
0
p
(
s
)(
cos
(
ωs
))
d
s,
γ
+
µ
+
αβ
2
S
0
γ
+
µ
−
(
ℜ
0
−
1)(
α
+
µ
)
ω
=
−
αβ
2
S
0
τ
0
p
(
s
)(
sin
(
ωs
))
d
s.
(4.5)
由
τ
0
p
(
s
)
d
s
=1
,
可知
τ
0
p
(
s
)(
−
ω
2
+
αβ
2
S
0
−
(
ℜ
0
−
1)(
α
+
µ
)(
γ
+
µ
))
d
s
=
αβ
2
S
0
−
(
ℜ
0
−
DOI:10.12677/aam.2021.1062182088
应用数学进展
忽丹丹,杨雨
1)(
α
+
µ
)(
γ
+
µ
)
−
ω
2
.
那么
(4.5)
就等价于
−
ω
2
+
αβ
2
S
0
−
(
R
0
−
1)(
α
+
µ
)(
γ
+
µ
)=
αβ
2
S
0
cos
(
ωτ
)
,
γ
+
µ
+
αβ
2
S
0
γ
+
µ
−
(
R
0
−
1)(
α
+
µ
)
ω
=
−
αβ
2
S
0
sin
(
ωτ
)
.
(4.6)
对
(4.6)
的两个方程各自平方后相加可得
ω
4
+
pω
2
+
q
=0
.
(4.7)
令
z
=
ω
2
,
方程
(4.7)
可记为
f
(
z
)=
z
2
+
pz
+
q
=0
,
(4.8)
这里
p
=2(
ℜ
0
−
1)(
α
+
µ
)(
γ
+
µ
)+
γ
+
µ
+
αβ
2
S
0
γ
+
µ
2
−
2
αβ
2
S
0
−
2(
ℜ
0
−
1)(
α
+
µ
)
γ
+
µ
+
αβ
2
S
0
γ
+
µ
=
−
2(
ℜ
0
−
1)(
α
+
µ
)
αβ
2
S
0
γ
+
µ
+(
γ
+
µ
)
2
+
αβ
2
S
0
γ
+
µ
2
,
q
=((
ℜ
0
−
1)(
α
+
µ
)(
γ
+
µ
))
2
−
2
αβ
2
S
0
(
ℜ
0
−
1)(
α
+
µ
)(
γ
+
µ
)
.
如果
ℜ
0
<
1
,
可得
p>
0
和
q>
0
,
对于任意的
τ>
0
都会使得
f
(
z
)
>
0
,
这与
f
(
z
)=0
形成
矛盾
.
这说明方程
(4.2)
所有特征值都具有负实部,即无病平衡点
E
0
是局部渐近稳定的
.
2)
系统
(1.2)
在平衡点
E
∗
(
S
∗
,I
∗
,E
∗
)
处的
Jacobi
矩阵为
J
(
S
∗
,I
∗
,E
∗
)=
−
β
1
E
∗
−
β
2
I
∗
−
µ
−
β
1
S
∗
−
β
2
S
∗
β
1
E
∗
+
β
2
I
∗
β
1
S
∗
−
(
α
+
µ
)
β
2
S
∗
0
α
τ
0
p
(
s
)
e
−
λs
d
s
−
(
γ
+
µ
)
.
(4.9)
因为
E
∗
(
S
∗
,I
∗
,E
∗
)
是系统
(1.2)
的正平衡点,我们记
β
1
S
∗
−
(
α
+
µ
)=
−
αβ
2
S
∗
µ
+
γ
,
µ
+
β
1
E
∗
+
β
2
I
∗
=
µ
ℜ
0
,
β
1
S
∗
(
β
1
E
∗
+
β
2
I
∗
)=
β
1
Λ
ℜ
0
(
ℜ
0
−
1)
.
(4.10)
系统
(1.2)
在
E
∗
(
S
∗
,I
∗
,E
∗
)
处的特征方程可以写成
DOI:10.12677/aam.2021.1062182089
应用数学进展
忽丹丹,杨雨
(
λ
+
µ
ℜ
0
)
λ
2
+(
γ
+
µ
+
αβ
2
S
∗
γ
+
µ
)
λ
+
αβ
2
S
∗
+
β
1
Λ
ℜ
0
(
ℜ
0
−
1)(
λ
+
γ
+
µ
)
=(
λ
+
µ
)
αβ
2
S
∗
τ
0
p
(
s
)
e
−
λs
d
s.
我们采用反证法,假如存在一个特征值
˜
λ
,
并且其实部非负
.
通过
(4),
显然有
˜
λ
2
+(
γ
+
µ
+
αβ
2
S
∗
γ
+
µ
)
˜
λ
+
αβ
2
S
∗
>
αβ
2
S
∗
τ
0
p
(
s
)
e
−
˜
λs
d
s
.
当
ℜ
0
>
1
,
我们知道
˜
λ
+
µ
ℜ
0
>
˜
λ
+
µ
以及
β
1
Λ
ℜ
0
(
ℜ
0
−
1)(
˜
λ
+
γ
+
µ
)
>
0
.
系统
(1.2)
的特征
方程方程
(4)
等号左端显然是大于该方程的右端
.
这是矛盾的,也就是说特征方程方程
(4)
的根都具
有负实部
.
因此,系统
(1.2)
的地方性平衡点
E
∗
是局部渐近稳定的
.
接下来我们主要讨论系统
(1.2)
无病平衡点和地方性平衡点的全局渐近稳定性
.
定理
4.2
1)
当
ℜ
0
<
1
时,系统
(1.2)
的无病平衡点
E
0
是全局渐近稳定的
;
2)
当
ℜ
0
>
1
时,系统
(1.2)
的地方性平衡点
E
∗
是全局渐近稳定的
.
证明
1)
构造
Lyapunov
泛函如下
V
1
(
t
)=
ℜ
0
E
(
t
)+
β
2
S
0
γ
+
µ
I
(
t
)+
αβ
2
S
0
γ
+
µ
τ
0
s
0
p
(
s
)
E
(
t
−
σ
)
d
σ
d
s.
(4.11)
V
1
(
t
)
沿着系统
(1.2)
的解的方向关于时间
t
求导
,
得到
d
V
1
(
t
)
d
t
=
ℜ
0
d
E
(
t
)
d
t
+
β
2
S
0
γ
+
µ
d
I
(
t
)
d
t
+
d
d
t
αβ
2
S
0
γ
+
µ
τ
0
s
0
p
(
s
)
E
(
t
−
σ
)
d
σ
d
s
=
ℜ
0
(
β
1
SE
+
β
2
SI
−
αE
−
µE
)+
αβ
2
S
0
γ
+
µ
τ
0
d
d
t
s
0
p
(
s
)
E
(
t
−
σ
)
d
σ
d
s
+
β
2
S
0
γ
+
µ
α
τ
0
E
(
t
−
s
)
p
(
s
)
d
s
−
(
γ
+
µ
)
I
=
ℜ
0
(
β
1
SE
+
β
2
SI
)
−
β
1
S
0
E
−
β
2
S
0
I
−
αβ
2
S
0
γ
+
µ
E
+
αβ
2
S
0
γ
+
µ
τ
0
E
(
t
−
s
)
p
(
s
)
d
s
+
αβ
2
S
0
γ
+
µ
τ
0
(
E
(
t
)
−
E
(
t
−
s
))
p
(
s
)
d
s
=
ℜ
0
(
β
1
SE
+
β
2
SI
)
−
β
1
S
0
E
−
β
2
S
0
I
≤
(
ℜ
0
−
1)(
β
1
S
0
E
+
β
2
S
0
I
)
.
显然
,
如果
ℜ
0
<
1
,
就有
d
V
1
(
t
)
d
t
≤
0
.
其中
d
V
1
(
t
)
d
t
=0
,
当且仅当
S
=
S
0
,
E
=
I
=0
.
也就是说
E
0
=
(
λ
µ
,
0
,
0)
是
Γ
1
=
(
S,E,I
)
∈
Ω
|
d
V
1
d
t
=0
的最大不变集
.
由
LaSalle’s
不变性原理
,
我们可以
得出,
ℜ
0
<
1
时
,
无病平衡点
E
0
是全局渐近稳定的
.
2)
定义如下
Lyapunov
泛函
V
2
(
t
)=
V
S
+
V
E
+
V
I
+
V
+
,
(4.12)
DOI:10.12677/aam.2021.1062182090
应用数学进展
忽丹丹,杨雨
这里
V
S
=
S
−
S
∗
−
S
∗
ln
S
S
∗
,V
E
=
E
−
E
∗
−
E
∗
ln
E
E
∗
,V
I
=
β
2
S
∗
I
∗
αE
∗
I
−
I
∗
−
I
∗
ln
I
I
∗
,
V
+
=
β
2
S
∗
I
∗
τ
0
s
0
p
(
s
)
E
(
t
−
σ
)
E
∗
−
1
−
ln
E
(
t
−
σ
)
E
∗
d
σ
d
s.
将这四个函数
V
S
,V
E
,V
I
,V
+
沿着系统
(1.2)
的解的方向关于时间
t
求导可得
d
V
S
(
t
)
d
t
=
S
−
S
∗
S
(Λ
−
β
1
SE
−
β
2
SI
−
µS
)
=
S
−
S
∗
S
(
β
1
S
∗
E
∗
−
β
1
SE
+
β
2
S
∗
I
∗
−
β
2
SI
+
µS
∗
−
µS
)
=
µS
∗
2
−
S
S
∗
−
S
∗
S
+
β
1
S
∗
E
∗
1
−
S
∗
S
−
SE
S
∗
E
∗
+
E
E
∗
+
β
2
S
∗
I
∗
1
−
S
∗
S
−
SI
S
∗
I
∗
+
I
I
∗
.
d
V
E
(
t
)
d
t
=
E
−
E
∗
E
(
β
1
SE
+
β
2
SI
−
αE
−
µE
)
=
E
−
E
∗
E
β
1
SE
+
β
2
SI
−
β
1
S
∗
E
−
β
2
S
∗
I
∗
E
E
∗
=
β
2
S
∗
I
∗
1+
SI
S
∗
I
∗
−
E
E
∗
−
SIE
∗
S
∗
I
∗
E
+
β
1
S
∗
E
∗
1+
SE
S
∗
E
∗
−
S
S
∗
−
E
E
∗
.
同样的
d
V
I
(
t
)
d
t
=
I
−
I
∗
I
β
2
S
∗
I
∗
αE
∗
α
τ
0
E
(
t
−
s
)
p
(
s
)
d
s
−
(
γ
+
µ
)
I
=
I
−
I
∗
I
β
2
S
∗
I
∗
αE
∗
α
τ
0
E
(
t
−
s
)
p
(
s
)
d
s
−
αE
∗
I
∗
I
=
β
2
S
∗
I
∗
1+
τ
0
E
(
t
−
τ
)
E
∗
p
(
s
)
d
s
−
τ
0
I
∗
E
(
t
−
τ
)
IE
∗
p
(
s
)
d
s
−
I
I
∗
.
d
V
+
(
t
)
d
t
=
β
2
S
∗
I
∗
d
d
t
τ
0
s
0
p
(
s
)
E
(
t
−
σ
)
E
∗
−
1
−
ln
E
(
t
−
σ
)
E
∗
d
σ
d
s
=
β
2
S
∗
I
∗
τ
0
p
(
s
)
d
d
t
s
0
E
(
t
−
σ
)
E
∗
−
1
−
ln
E
(
t
−
σ
)
E
∗
d
σ
d
s
=
β
2
S
∗
I
∗
τ
0
p
(
s
)
−
E
(
t
−
s
)
E
∗
+
E
E
∗
+
ln
E
(
t
−
s
)
E
d
s
=
β
2
S
∗
I
∗
τ
0
p
(
s
)
E
E
∗
d
s
−
β
2
S
∗
I
∗
τ
0
p
(
s
)
E
(
t
−
s
)
E
∗
d
s
+
β
2
S
∗
I
∗
τ
0
p
(
s
)
ln
E
(
t
−
s
)
E
d
s.
DOI:10.12677/aam.2021.1062182091
应用数学进展
忽丹丹,杨雨
结合
d
V
S
(
t
)
d
t
,
d
V
E
(
t
)
d
t
,
d
V
I
(
t
)
d
t
,
d
V
+
(
t
)
d
t
,我们可以得到
d
V
2
(
t
)
d
t
=
d
V
S
(
t
)
d
t
+
d
V
E
(
t
)
d
t
+
d
V
I
(
t
)
d
t
+
d
V
+
(
t
)
d
t
=(
µS
∗
+
β
1
S
∗
E
∗
)
2
−
S
S
∗
−
S
∗
S
+
β
2
S
∗
I
∗
τ
0
p
(
s
)
3
−
S
∗
S
−
SIE
∗
S
∗
I
∗
E
−
I
∗
I
E
(
t
−
s
)
E
∗
+
ln
E
(
t
−
s
)
E
∗
d
s
=(
µS
∗
+
β
1
S
∗
E
∗
)
2
−
S
S
∗
−
S
∗
S
+
β
2
S
∗
I
∗
1
−
S
∗
S
+
ln
S
∗
S
+
β
2
S
∗
I
∗
1
−
SIE
∗
S
∗
I
∗
E
+
ln
SIE
∗
S
∗
I
∗
E
(3
.
13)
+
β
2
S
∗
I
∗
τ
0
p
(
s
)
1
−
I
∗
E
(
t
−
s
)
IE
∗
+
ln
I
∗
E
(
t
−
s
)
IE
∗
d
s.
(3
.
14)
由算术平均数大于或等于几何平均数可知
2
−
S
S
∗
−
S
∗
S
≤
0
.
又由于函数
H
(
x
)=1
−
h
(
x
)+
ln
h
(
x
)
,x
∈
R
+
,
对任意
h
(
x
)
≥
0
都有
H
(
x
)
≤
0
,
当且仅当
h
(
x
)=1
时有最小值
H
(1)=0
.
因此
(3.13)
和
(3.14)
是非正的
.
由上述不等式可知
,
正定函数
V
2
(
t
)
具有非正导数
d
V
2
(
t
)
d
t
.
其中
,
d
V
2
(
t
)
d
t
=0
当且仅当
S
(
t
)=
S
∗
,E
(
t
−
s
)=
E
(
t
)=
E
∗
,I
(
t
)=
I
∗
.
Γ
2
=
(
S
(
t
)
,E
(
t
)
,I
(
t
))
∈
Ω
|
d
V
2
d
t
=0
的最大不变集为
E
∗
.
因此,由
LaSalle’s
不变性原理
,
我们可以得出,
ℜ
0
>
1
时
,
地方性平衡点
E
∗
是全局渐近稳定的
.
5.
数值模拟
在模型
(1.2)
中
,
当取
p
(
s
)=
δ
(
s
−
τ
)
,
这里
δ
为狄拉克函数,模型
(1.2)
可简化为离散时滞的
SEI
模型
d
S
(
t
)
d
t
=Λ
−
β
1
SE
−
β
2
SI
−
µS,
d
E
(
t
)
d
t
=
β
1
SE
+
β
2
SI
−
(
α
+
µ
)
E,
d
I
(
t
)
d
t
=
αE
(
t
−
τ
)
−
(
γ
+
µ
)
I.
(5.1)
这里我们选取
[14]
中新冠肺炎潜伏期的中位数
τ
=7
.
设模型
(5.1)
中各参数值为
Λ=0
.
1
,β
1
=
0
.
4
,β
2
=0
.
9
,µ
=0
.
3
,α
=0
.
15
,γ
=0
.
35
.
初值
ϕ
1
=0
.
5
,ϕ
2
=0
,ϕ
3
=0
.
3
.
经计算
ℜ
0
=0
.
4501
<
1
.
由定理
(4.2)
可知,系统
(5.1)
的无病平衡点是全局渐近稳定的,如图
1
所示
.
从图
1
中我们可以看到,
随着时间的推移,无症状感染者与发病者的数量都趋于零
.
也就是说
,
当
ℜ
0
<
1
时,疾病终将消
亡
.
数值模拟结果与定理结论相符
.
再设系统
(5.1)
中各参数值为
Λ=0
.
22
,β
1
=0
.
6
,β
2
=1
.
3
,µ
=0
.
3
,α
=0
.
15
,γ
=0
.
15
.
初值为
ϕ
1
=0
.
45
,ϕ
2
=0
.
05
,ϕ
3
=0
.
35
.
经计算
ℜ
0
=1
.
684
>
1
.
如图
2
所示
,
该数据验证了定理
(4.2),
系
统
(5.1)
是持久的
,
且
t
趋于无穷的时候
,
疾病将持久存在
.
DOI:10.12677/aam.2021.1062182092
应用数学进展
忽丹丹,杨雨
Figure1.
ℜ
0
<
1
,thediseasewilldie
图
1.
当
ℜ
0
<
1
时
,
疾病终将消亡
Figure2.
ℜ
0
>
1
,thediseasewillpersist
图
2.
当
ℜ
0
>
1
时
,
疾病将会持续
理论分析结果和数值拟合表明
,
系统
(1.1)
的全局动力学性质仅依赖于基本再生数
,
即当基本再
生数
ℜ
0
<
1
时
,
无病平衡点全局渐近稳定
,
疾病将会消亡
;
当基本再生数
ℜ
0
>
1
时
,
地方病平衡点
全局渐近稳定
,
疾病将会持续
.
从数值模拟中我们可以看出
,
当减小基本再生数时
,
疾病可以从持续
流行变成消亡
.
由基本再生数的表达式
(3.2)
可以看出
,
降低
β
1
,β
2
会减缓疫情的暴发
,
潜伏者与感
染者人数也会随之降低
.
同样的
,
当
Λ
减小时
,
也即易感者人数减少时
,
基本再生数也会随之减小
,
使疾病消亡
.
接种疫苗可以降低易感染者的人数
,
因此政府可以采取核酸检测、隔离以及疫苗接种
等措施降低潜伏期患者以及感染者与易感染者的接触传播
,
从而降低阈值
ℜ
0
,
延缓疾病暴发
.
参考文献
[1]
陆征一
,
周义仓
.
数学生物学进展
[J].
北京
:
科学出版社
,2006.
DOI:10.12677/aam.2021.1062182093
应用数学进展
忽丹丹,杨雨
[2]Hamer,W.H.(1906)EpidemicDiseaseinEngland.
TheLancet
,
1
,733-739.
[3]Kermack, W.O.andMcKendrick,A.G. (1927)AContribution totheMathematicalTheory of
Epidemics.
ProceedingsoftheRoyalSocietyofLondon.SeriesA
,
115
,700-721.
https://doi.org/10.1098/rspa.1927.0118
[4]Wang,C.,Horby,P.W.,Hayden,F.G.,
etal.
(2020)ANovelCoronavirusOutbreakofGlobal
HealthConcern.
TheLancet
,
395
,470-473.https://doi.org/10.1016/S0140-6736(20)30185-9
[5]Worldometer(2021)COVID-19CoronavirusPandemic.
https://www.worldometers.info/coronavirus/
[6]
张龙浩
,
李柏宏
,
贾鹏
,
蒲剑
,
白蓓
,
李音
,
朱培嘉
,
李雷
,
曾国军
,
赵欣
,
董珊珊
,
刘梦菡
,
张楠
.
新
型冠状病毒
(SARS-CoV-2)
全球研究现状分析
[J].
生物医学工程学杂志
,2020, 37(2):236-249.
[7]Yang,C.Y.andWang,J.(2020)AMathematicalModelfortheNovelCoronavirusEpidemic
inWuhan,China.
MathematicalBiosciencesandEngineering
,
17
,2708-2724.
https://doi.org/10.3934/mbe.2020148
[8]
王冰杰
.
基于潜伏期有传染力的
SEIR
传染病模型的控制策略
[J].
东北师大学报
(
自然科学版
),
2014,46(1):28-32.
[9]Li,G.H.andJin,Z.(2005)GlobalStabilityofaSEIREpidemicModelwithInfectiousForce
inLatent,InfectedandImmunePeriod.
Chaos,SolitonsandFractals
,
25
,1177-1184.
https://doi.org/10.1016/j.chaos.2004.11.062
[10]Çakan,S.(2020)DynamicAnalysisofaMathematicalModelwithHealthCareCapacityfor
COVID-19Pandemic.
Chaos,SolitonsandFractals
,
139
,ArticleID:110033.
https://doi.org/10.1016/j.chaos.2020.110033
[11]Liu,M.,Liu,X.andLiang,W.N.(2020)HandbookofBasicEpidemiologicalKnowledge
onCoronavirusDisease2019(COVD-19)-110QuestionsAnswers.PekingUniversityMedical
Press,Beijing.
[12]Huang,G. and Liu,A.P. (2013)ANote on GlobalStability for aHeroin Epidemic Modelwith
DistributedDelay.
AppliedMathematicsLetters
,
26
,687-691.
https://doi.org/10.1016/j.aml.2013.01.010
[13]VandenDriessche,P.andWatmough,J.(2002)ReproductionNumbersandSub-Threshold
EndemicEquilibriaforCompartmentalModelsofDiseaseTransmission.
MathematicalBio-
sciences
,
180
,29-48.https://doi.org/10.1016/S0025-5564(02)00108-6
[14]
刘小惠
,
何阳
,
麻先思
,
罗良清
.
有关新冠肺炎潜伏期和疑似期的统计数据分析
:
基于湖北省外
2172
条确诊数据
[J].
应用数学学报
,2020,43(2):278-294.
DOI:10.12677/aam.2021.1062182094
应用数学进展