具有时滞的SIR计算机病毒模型的稳定性分析_免费论文全文下载

申明:本网站内容仅用于学术交流，如有侵犯您的权益，请及时告知我们，本站将立即删除有关内容。

【摘 要】研究了一类具有非线性发生率的时滞SIR传染病模型. 确定了决定计算机病毒消失或继续存在的基本再生数， 通过分析系统对应的特征方程，得到无病平衡点与地方平衡点的局部稳定性。通过构造适当的Lyapunov函数，利用LaSalle不变原理，证明当基本再生数小于1时，无病平衡点是全局渐近稳定的； 最后，通过MATLAB进行数值模拟验证了所得理论分析结果的正确性。

【关键词】时滞； 传染病模型；全局稳定性； Lyapunov函数；基本再生数
Stability Analysis of a Delayed SIR Computer Virus Model
LEI Xue-hong XU Yun-xia
（Department of mathematical science， Kaili College， Kaili Guizhou 556011，China）
【Abstract】In this paper， A delayed SIR epidemic model with nonlinear incidence rate is investigated. the basic reproduction number is determining whether the disease dies is found， and the existence of the model is discussed. By analyzing the corresponding characteristic equation.the local stability of a disease-free equilibrium and endemic equilibrium are discussed. According to the suitable Lyapunov function and LaSalle invariance principle， it is proved that the disease-free equilibrium is globally asymptotically stable as the basic reproduction number for viral infection is less than unity， Finally，the theoretical results obtained are verified by numerical simulations for the numerical model with a specific incidence by MATLAB.
【Key words】Time delay； Epidemic model； Global stability； Lyapunov Function；Basic Reproduction Number
0 引言
在�^去的几十年里，互联网的迅速普及。然而，随着计算机迅速普及也极大提高了计算机病毒的传播能力。由于计算机病毒具有极大的破坏性、不可预测性，多态性等特点， 早已成为现代信息社会的重要威胁之一。随着计算机与通信技术的快速发展，计算机病毒程序也变得越来越复杂，人们开始运用传染病动力学中仓室模型来了解计算机病毒传播的一般规律，并得到了大量的研究成果[1-2]。最典型的仓室模型是Kephart借鉴经典的传染病SIS模型建立了最早的计算机病毒传播仓室模型，宣告计算机病毒传播动力学的诞生。SIS模型则用来描述计算机康复后不具有免疫力再次被病毒感染的情况。大量的临床研究发现，由于个体的免疫系统完善水平不同，对于相同的病毒，并不是每个个体感染后都能产生抵抗被该病毒再次感染的免疫力。因此 ，更合理的发生率应该是非线性的。随着研究的不断深入，人们开始引入不同的非线性传染率；如带有非线性传染率为βSqI的SIS传染病毒模型；带有非线性传染率的的SIR或SIRS模型，带有非线性传染病毒βI（1+αIk-1）S的SIR模型，带有非线性传染率的的SEIV模型，带有非线性传染病毒的模型和带有非线性传染率的SIR模型；本文研究如下
S'（t）=？撰--μSI'（t）=-（μ+r+ε）IR'（t）=rI-μR（1）
其中S（t），I（t），R（t）分别表示在t时刻因特网中易感节点、感染节点，免疫节点所占的比例。设？撰表示单位时间内新增加易感节点的数量，μ表示因故障等引起的宕机率，β表示平均有效接触率，ε表示因病毒感染病引起的宕机率，r表示为感染节点到免疫节点的恢复系数。τ为时滞，即表示病毒的潜伏期，0μ（μ+r+ε）时，还有一个病毒平衡点P*（S*，I*，R*）.I*=；S*=；R*=.
令R0=为系统（1）的基本再生数。
在本节，通过讨论系统（1）对应的特征方程来讨论平衡点的局部渐近稳定性。 定理1：当R01时，系统（1）的无病平衡点P0是不稳定的。
证明：系统（1）的无病平衡点P0的线性特征方程为（λ+μ）2（λ-βe-（δ+λ）+μ+r+ε）=0（3）
λ1=λ2=-μ，
λ3β？撰e-（δ+λ）τ-（μ+r+ε），
令方程f（λ）=λ-β？撰e-（δ+λ）τ+（μ+r+ε）（4）
由文献[5]定理1当R01时，
f（0）=-e-δτ+（μ+r+ε）1时，系统（1）不稳定的。
当R0>1时，系统（1）在病毒平衡点P*处的特征方程为
λ3+a1λ2+a2λ+a3+（b1λ2+b2λ+b3）e-λτ=0（5）
其中：a1=3μ+ε+γ+；
b1=-
a2=3μ2+2με+2μγ+（2μ+ε+γ）
b2=-
a3=μ（μ+ε+γ）（μ+）；
b3=-；
当τ=0时，有λ3+c1λ2+c2λ+c3=0.（6）
c1=（3μ+ε+γ）αI*+，
c2=
c3=I*；
此时有c1>0，c2>0，c3>0，
容易计算c1c2-c3>0.
因此，由Routh-Hurwitz判断定理可知，平衡点P*局部渐近稳定性的。
当τ>0时，设λ=iω（ω>0）是方程（4）的一个根，分离实部与虚部，得
b2ωcos（ωτ）+（b1ω2-b3）sin（ωτ）=ω3-a2ω（b3-b1ω2）cos（ωτ）+b2ωsin（ωτ）=a1ω2-a3（7）
（7）式两个平方和，得ω6+r1ω4+r2ω2+r3=0
其中：r1=a21-2a2-b21；
r2=a22-2a1a3+2b1b3-b22；r3=a23-b23；
令z=ω2，Δ=r21-3r2，
g（z）=z3+r1z2+r2z+r3（8）
定理2 若R0>1成立；
（1）当τ=0时；
（2）当τ>0时，r3>0且Δ≤0；
（3）当τ>0时，r31时解是局部渐近的。通过构造Lyapunov函数，利用LaSalle不变原理，证明当R01时，得到了地方平衡点稳定的充分条件。最后，选取适当的数值，通过MATLAB进行数值模拟，进一步验证了所得的主要结果。
【参考文献】
[1]Liu W M， Levin S A， Iw asa Y.Influence of nonlinear incidence rates upon the behavior of SIRS epidemiological models[J].Joumal of Mathematical Biology，1986， 23（2）：187-204.
[2]胡宝安，李兵，李亚玲.具有时滞的SIR计算机病毒传播模型[J].计算机工程.2016.42（5）：168-172
[3]杨亚莉，李建全，刘万萌，等.一类具有分布时滞和非线性发生率的媒介传染病模型的全局稳定性[J].应用数学和力学Vol.34，No.12.Dec，2013，1291-1299.
[4]陈辉，徐瑞.一类具有饱和感染率和胞内时滞的病毒感染模型的稳定性和Hopf分支[J].应用数学，2016，29（2）：398-408.
[5]杨俊仙，闫萍.一类具有饱和率的时滞SEIR传染病模型的分析[J].中山大学学报（自然科学版）.Vol.54，No.3.May，2015，51-55.
[6]徐昌进，姚凌云.具有时滞的计算机网络病毒传染模型分支分析[J].河南科技大学学报（自然科学版）Vol.34，No.1.Feb.，2013，55-59.
[7]叶志勇，刘原，赵彦勇.一类SIQR传染病模型在无尺度网络上的传播行为分析[J].计算机工程与科学.Vol.36，No.8.Aug. 2014，1524-1527.
[8]邓长松，刘启明.无标度网络上具有时滞的计算机病毒传播模型研究[J].河北师范大学学报（自然科学版）Vol.39，No.2.Mar.，2015，111-115.
[9]李盈科，李浓溟，杨文成，等.具有时滞的生态―传染病模型的定性分析[J].数学的实践与认识.Vol.45，No.13.Jun.，2015，273-277.
[责任编辑：朱丽娜]

本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系，我们将在24小时内删除。QQ:522-52-5970