# Stability and Hopf Bifurcation for a Delayed SLBRS Computer Virus Model

Authors: Zizhen Zhang; Huizhong Yang
Journal: The Scientific World Journal
 (2014)

---

## Abstract

By incorporating the time delay due to the period that computers use antivirus software to clean the virus into the SLBRS model a delayed SLBRS computer virus model is proposed in this paper. The dynamical behaviors which include local stability and Hopf bifurcation are investigated by regarding the delay as bifurcating parameter. Specially, direction and stability of the Hopf bifurcation are derived by applying the normal form method and center manifold theory. Finally, an illustrative example is also presented to testify our analytical results.


---

## Body

## 1. Introduction

For the purpose of studying the prevalence of computer viruses and posing effective strategies of controlling their spread across the Internet and considering the similarity between the propagation of computer virus across the Internet and that of biological viruses across a population, many scholars have proposed computer virus models by properly modifying some epidemiological models [1–6]. In [2], Han and Tan proposed a computer virus model using an SIRS model and investigated the local stability and Hopf bifurcation for the endemic state of the model. In [4], Mishra and Pandey proposed an SEIRS model for the transmission of worms in computer network through vertical transmission and they obtained the threshold value determining whether the worms die out. All the models above assume that an infected computer in which the virus resides is in latency and cannot infect other computers.Recently, Yang et al. [7] proposed the following SLBRS computer virus model by taking into account the fact that a computer immediately possesses infection ability once it is infected. Consider (1)
d
S
(
t
)
d
t
=
μ
-
β
S
(
t
)
(
L
(
t
)
+
B
(
t
)
)
-
μ
S
(
t
)
+
α
R
(
t
)
,
d
L
(
t
)
d
t
=
β
S
(
t
)
(
L
(
t
)
+
B
(
t
)
)
-
ε
L
(
t
)
-
μ
L
(
t
)
,
d
B
(
t
)
d
t
=
ε
L
(
t
)
-
γ
B
(
t
)
-
μ
B
(
t
)
,
d
R
(
t
)
d
t
=
γ
B
(
t
)
-
α
R
(
t
)
-
μ
R
(
t
)
,
where S
(
t
) denotes the percentage of computers that have no immunity at time t; L
(
t
) denotes the percentage of infected computers that are latent at time t; B
(
t
) denotes the percentage of infected computers that are breaking out at time t; and R
(
t
) denotes the percentage of uninfected computers that have temporary immunity at time t. μ is the constant rate for the external computers connecting to the Internet and that for the internal computers disconnecting from the Internet. β is the rate of contact for the virus-free computer. α is the rate coefficient from class R to class S, ε is the rate coefficient from class L to class B, and γ is the rate coefficient from class B to class R. Yang et al. [7] studied global stability of the model by constructing a Lyapunov function for the model.It is well known that there may be a time lag when the computer uses antivirus software to clean the virus and computer virus models with time delay have been studied by many scholars [1, 2, 6, 8–15]. In [8], Feng et al. investigated the Hopf bifurcation of delayed SIRS computer virus propagation model by regarding the latent period of computer virus as a bifurcation parameter. In [10], Liu investigated the Hopf bifurcation of a delayed SEIQRS model for the transmission of malicious objects in computer network by regarding the time delay due to the temporary immunity period as a bifurcation parameter. To the best of our knowledge, few papers deal with the research of Hopf bifurcation of system (1). Motivated by the work above and considering that there may be a time lag when the computer uses antivirus software to clean the virus, we incorporate time delay into system (1) and investigate the Hopf bifurcation of the following delayed SLBRS system:
(2)
d
S
(
t
)
d
t
=
μ
-
β
S
(
t
)
(
L
(
t
)
+
B
(
t
)
)
-
μ
S
(
t
)
+
α
R
(
t
)
,
d
L
(
t
)
d
t
=
β
S
(
t
)
(
L
(
t
)
+
B
(
t
)
)
-
ε
L
(
t
)
-
μ
L
(
t
)
,
d
B
(
t
)
d
t
=
ε
L
(
t
)
-
γ
B
(
t
-
τ
)
-
μ
B
(
t
)
,
d
R
(
t
)
d
t
=
γ
B
(
t
-
τ
)
-
α
R
(
t
)
-
μ
R
(
t
)
,
where τ is the time delay due to the period that a computer uses antivirus software to clean the virus.The main purpose of this paper is to investigate the effect of delay on system (2). This paper is organized as follows. In Section 2, sufficient conditions for local stability and existence of local Hopf bifurcation are obtained by analyzing the distribution of the roots of the associated characteristic equation. In Section 3, direction and stability of the Hopf bifurcation are studied by applying the normal form theory and center manifold theorem. In Section 4, some numerical simulations are presented to verify the theoretical results.


## 2. Stability and Existence of Local Hopf Bifurcation

By a simple computation, we can conclude that ifR
0
>
1, then system (2) has a unique positive equilibrium D
*
(
S
*
,
L
*
,
B
*
,
R
*
), where
(3)
S
*
=
(
ε
+
μ
)
(
γ
+
μ
)
β
(
ε
+
γ
+
μ
)
,
L
*
=
μ
(
α
+
μ
)
(
γ
+
μ
)
(
1
-
1
/
R
0
)
(
α
+
μ
)
(
ε
+
μ
)
(
γ
+
μ
)
-
α
ε
γ
,
B
*
=
ε
L
*
γ
+
μ
,
R
*
=
γ
B
*
α
+
μ
and R
0
=
β
(
ε
+
γ
+
μ
)
/
(
ε
+
μ
)
(
γ
+
μ
) is called the basic reproduction number.The linearized system of system (2) is
(4)
d
S
(
t
)
d
t
=
a
1
S
(
t
)
+
a
2
L
(
t
)
+
a
3
B
(
t
)
+
a
4
R
(
t
)
,
d
L
(
t
)
d
t
=
a
5
S
(
t
)
+
a
6
L
(
t
)
+
a
7
B
(
t
)
,
d
B
(
t
)
d
t
=
a
8
L
(
t
)
+
a
9
B
(
t
)
+
b
1
B
(
t
-
τ
)
,
d
R
(
t
)
d
t
=
a
10
R
(
t
)
+
b
2
B
(
t
-
τ
)
,
where
(5)
a
1
=
-
μ
-
β
(
L
*
+
B
*
)
,
a
2
=
-
β
S
*
,
a
3
=
-
β
S
*
,
a
4
=
α
,
a
5
=
β
(
L
*
+
B
*
)
,
a
6
=
β
S
*
-
ε
-
μ
,
a
7
=
β
S
*
,
a
8
=
ε
,
a
9
=
-
μ
,
a
10
=
-
(
α
+
μ
)
,
b
1
=
-
γ
,
b
2
=
γ
.
The characteristic equation of system (4) is
(6)
λ
4
+
A
3
λ
3
+
A
2
λ
2
+
A
1
λ
+
A
0
+
(
B
3
λ
3
+
B
2
λ
2
+
B
1
λ
+
B
0
)
e
-
λ
τ
=
0
,
where
(7)
A
0
=
a
1
a
6
a
9
a
10
+
a
8
a
10
(
a
3
a
5
-
a
1
a
7
)
,
A
1
=
a
7
a
8
a
10
+
a
8
(
a
1
a
7
-
a
3
a
5
)
-
a
1
a
6
(
a
9
+
a
10
)
-
a
9
a
10
(
a
1
+
a
6
)
,
A
2
=
a
1
a
6
+
a
9
a
10
-
a
7
a
8
+
(
a
1
+
a
6
)
(
a
9
+
a
10
)
,
A
3
=
-
(
a
1
+
a
6
+
a
9
+
a
10
)
,
B
0
=
a
10
b
1
(
a
1
a
6
-
a
2
a
5
)
,
B
1
=
(
a
2
a
5
-
a
1
a
6
-
a
1
a
10
-
a
6
a
10
)
b
1
,
B
2
=
(
a
1
+
a
6
+
a
10
)
b
1
,
B
3
=
-
b
1
.
For τ
=
0, (6) reduces to
(8)
λ
4
+
A
13
λ
3
+
A
12
λ
2
+
A
11
λ
+
A
10
=
0
,
where
(9)
A
10
=
A
0
+
B
0
,
A
11
=
A
1
+
B
1
,
A
12
=
A
2
+
B
2
,
A
13
=
A
3
+
B
3
.Obviously, if condition (H
1) (10) holds, D
* is locally asymptotically stable in the absence of delay according to the Routh-Hurwitz criterion. One has
(10)
Det
1
=
A
13
>
0
,
Det
2
=
|
A
13
1
A
11
A
12
|
>
0
,
Det
3
=
|
A
13
1
0
A
11
A
12
A
13
0
A
10
A
11
|
>
0
,
Det
4
=
|
A
13
1
0
0
A
11
A
12
A
13
1
0
A
10
A
11
A
12
0
0
0
A
10
|
>
0
.
For τ
>
0, let λ
=
i
ω
(
ω
>
0
) be the root of (6). Then, we get
(11)
(
B
1
ω
-
B
3
ω
3
)
sin
τ
ω
+
(
B
0
-
B
2
ω
2
)
cos
⁡
τ
ω
=
A
2
ω
2
-
ω
4
-
A
0
,
(
B
1
ω
-
B
3
ω
3
)
cos
⁡
τ
ω
-
(
B
0
-
B
2
ω
2
)
sin
τ
ω
=
A
3
ω
3
-
A
1
ω
,
from which one can get
(12)
ω
8
+
c
3
ω
6
+
c
2
ω
4
+
c
1
ω
2
+
c
0
=
0
,
with
(13)
c
0
=
A
0
2
-
B
0
2
,
c
1
=
A
1
2
-
B
1
2
-
2
A
0
A
2
+
2
B
0
B
2
,
c
2
=
A
2
2
-
B
2
2
+
2
A
0
-
2
A
1
A
3
+
2
B
1
B
3
,
c
3
=
A
3
2
-
B
3
2
-
2
A
2
.
Let ω
2
=
v; then, (12) becomes
(14)
v
4
+
c
3
v
3
+
c
2
v
2
+
c
1
v
+
c
0
=
0
.
In order to give the main results in this paper, we make the following assumption.(H
2)
Equation (14) has at least one positive root.If condition (H
2) holds, we know that (12) has at least a positive root ω
0 such that (6) has a pair of purely imaginary roots ±
i
ω
0. For ω
0, the corresponding critical value of the delay is
(15)
τ
0
=
1
ω
0
arccos
p
6
ω
0
6
+
p
4
ω
0
4
+
p
2
ω
0
2
+
p
0
q
6
ω
0
6
+
q
4
ω
0
4
+
q
2
ω
0
2
+
q
0
,
where
(16)
p
0
=
-
A
0
B
0
,
p
2
=
A
0
B
2
-
A
1
B
1
+
A
2
B
0
,
p
4
=
A
1
B
3
-
A
2
B
2
+
A
3
B
1
-
B
0
,
p
6
=
B
2
-
A
3
B
3
,
q
0
=
B
0
2
,
q
2
=
B
1
2
-
2
B
0
B
2
,
q
4
=
B
2
2
-
2
B
1
B
3
,
q
6
=
B
3
2
.
Taking the derivative with respect to τ on both sides of (6), we obtain
(17)
[
d
λ
d
τ
]
-
1
=
-
4
λ
3
+
3
A
3
λ
2
+
2
A
2
λ
+
A
1
λ
(
λ
4
+
A
3
λ
3
+
A
2
λ
2
+
A
1
λ
+
A
0
)
+
3
B
3
λ
2
+
2
B
2
λ
+
B
1
λ
(
B
3
λ
3
+
B
2
λ
2
+
B
1
λ
+
B
0
)
.
Thus,
(18)
Re
[
d
λ
d
τ
]
τ
=
τ
0
-
1
=
f
′
(
v
*
)
B
3
2
ω
0
6
+
(
B
2
2
-
2
B
1
B
3
)
ω
0
4
+
(
B
1
2
-
2
B
0
B
2
)
ω
0
2
+
B
0
2
,
where v
*
=
ω
0
2 and f
(
v
)
=
v
4
+
c
3
v
3
+
c
2
v
2
+
c
1
v
+
c
0.Thus, if condition (H
3) f
′
(
v
*
)
≠
0, then Re
[
d
λ
/
d
τ
]
τ
=
τ
0
-
1
≠
0. According to the Hopf bifurcation theorem in [16], we have the following results.Theorem 1.
If conditions (H
1)–(H
3) hold, then the positive equilibrium D
*
(
S
*
,
L
*
,
B
*
,
R
*
) of system (2) is asymptotically stable for τ
∈
[
0
,
τ
0
) and system (2) undergoes a Hopf bifurcation at the positive equilibrium D
*
(
S
*
,
L
*
,
B
*
,
R
*
) when τ
=
τ
0.


## 3. Properties of the Hopf Bifurcation

Letu
1
(
t
)
=
S
(
t
)
-
S
*, u
2
(
t
)
=
L
(
t
)
-
L
*, u
3
(
t
)
=
B
(
t
)
-
B
*, u
4
(
t
)
=
R
(
t
)
-
R
*, and τ
=
τ
0
+
μ, μ
∈
R, and normalize the delay by t
→
(
t
/
τ
). Then, system (2) can be transformed into a functional differential equation (PDF) as
(19)
u
˙
(
t
)
=
L
μ
u
t
+
F
(
μ
,
u
t
)
,
where
(20)
u
t
=
(
u
1
(
t
)
,
u
2
(
t
)
,
u
3
(
t
)
,
u
4
(
t
)
)
T
∈
C
=
C
(
[
-
1,0
]
,
R
4
)
,
L
μ
ϕ
=
(
τ
0
+
μ
)
(
A
′
ϕ
(
0
)
+
B
′
ϕ
(
-
1
)
)
,
F
(
μ
,
u
t
)
=
(
τ
0
+
μ
)
(
-
β
(
ϕ
1
(
0
)
ϕ
2
(
0
)
+
ϕ
1
(
0
)
ϕ
3
(
0
)
)
β
(
ϕ
1
(
0
)
ϕ
2
(
0
)
+
ϕ
1
(
0
)
ϕ
3
(
0
)
)
0
0
)
,
where
(21)
A
′
=
(
a
1
a
2
a
3
a
4
a
5
a
6
a
7
0
0
a
8
a
9
0
0
0
0
a
10
)
,
B
′
=
(
0
0
0
0
0
0
0
0
0
0
b
1
0
0
0
b
2
0
)
.By the Riesz representation theorem, there is a4
×
4 matrix function with bounded variation components η
(
θ
,
μ
), θ
∈
[
-
1,0
] such that
(22)
L
μ
ϕ
=
∫
-
1
0
d
η
(
θ
,
μ
)
ϕ
(
θ
)
,
ϕ
∈
C
.
In fact, we choose
(23)
η
(
θ
,
μ
)
=
(
τ
0
+
μ
)
(
A
′
δ
(
θ
)
+
B
′
δ
(
θ
+
1
)
)
.
For ϕ
∈
C
(
[
-
1,0
]
)
,
R
4
), we define
(24)
A
(
μ
)
ϕ
=
{
d
ϕ
(
θ
)
d
θ
,
-
1
≤
θ
<
0
,
∫
-
1
0
d
η
(
θ
,
μ
)
ϕ
(
θ
)
,
θ
=
0
,
R
(
μ
)
ϕ
=
{
0
,
-
1
≤
θ
<
0
,
F
(
μ
,
ϕ
)
,
θ
=
0
.
Then, system (19) is equivalent to the following form:
(25)
u
˙
(
t
)
=
A
(
μ
)
u
t
+
R
(
μ
)
u
t
.
Next, we define the adjoint operator A
* of A as
(26)
A
*
(
μ
)
s
=
{
-
d
φ
(
s
)
d
s
,
0
<
s
≤
1
,
∫
-
1
0
d
η
T
(
s
,
μ
)
φ
(
-
s
)
,
s
=
0
and a bilinear inner product as
(27)
〈
φ
,
ϕ
〉
=
φ
-
(
0
)
ϕ
(
0
)
-
∫
θ
=
-
1
0
∫
ξ
=
0
θ
φ
-
(
ξ
-
θ
)
d
η
(
θ
)
ϕ
(
ξ
)
d
ξ
,
where η
(
θ
)
=
η
(
θ
,
0
).Letρ
(
θ
)
=
(
1
,
ρ
2
,
ρ
3
,
ρ
4
)
T
e
i
τ
0
ω
0
θ be the eigenvector of A
(
0
) corresponding to the eigenvalue +
i
τ
0
ω
0 and let ρ
*
(
θ
)
=
D
(
1
,
ρ
2
*
,
ρ
3
*
,
ρ
4
*
)
e
i
τ
0
ω
0
s be the eigenvector of A
*
(
0
) corresponding to the eigenvalue -
i
τ
0
ω
0. Then, from the definition of A
(
0
) and A
*
(
0
), we obtain
(28)
ρ
2
=
a
5
+
a
7
ρ
3
τ
0
ω
0
-
a
6
,
ρ
3
=
a
5
a
8
(
i
ω
0
-
a
6
)
(
i
ω
0
-
a
9
-
b
1
e
-
i
τ
0
ω
0
)
,
ρ
4
=
i
ω
0
-
a
1
-
a
2
ρ
2
-
a
3
ρ
3
a
4
,
ρ
2
*
=
-
i
ω
0
+
a
1
a
5
,
ρ
3
*
=
-
a
2
+
(
i
ω
0
+
a
6
)
ρ
2
*
a
8
,
ρ
4
*
=
-
a
3
+
a
7
ρ
2
*
+
(
i
ω
0
+
a
9
+
b
1
e
i
τ
0
ω
0
)
ρ
2
*
b
2
e
i
τ
0
ω
0
.
From (27), we can get
(29)
〈
q
*
(
s
)
,
q
(
θ
)
〉
=
ρ
-
[
1
+
ρ
2
ρ
-
2
*
+
ρ
3
ρ
-
3
*
+
ρ
4
ρ
-
4
*
e
-
i
τ
0
ω
0
iiiiiiiiiiiiiiiiiiiiiiiiiiii
+
τ
0
ρ
3
e
-
i
τ
0
ω
0
(
b
1
ρ
-
3
*
+
b
2
ρ
-
4
*
)
]
.
Then, we choose
(30)
D
-
=
[
1
+
ρ
2
ρ
-
2
*
+
ρ
3
ρ
-
3
*
+
ρ
4
ρ
-
4
*
e
-
i
τ
0
ω
0
iiiiii
+
τ
0
ρ
3
e
-
i
τ
0
ω
0
(
b
1
ρ
-
3
*
+
b
2
ρ
-
4
*
)
]
-
1
such that 〈
ρ
*
,
ρ
〉
=
1 and 〈
ρ
*
,
ρ
-
〉
=
0.Following the algorithms given in [16] and using similar computation process in [17], we can get the following important coefficients:
(31)
g
20
=
2
τ
0
D
-
β
(
ρ
-
2
*
ρ
3
-
ρ
2
)
,
g
11
=
τ
0
D
-
β
(
ρ
-
2
*
(
ρ
3
+
ρ
-
3
)
-
(
ρ
2
+
ρ
-
2
)
)
,
g
02
=
2
τ
0
D
-
β
(
ρ
-
2
*
ρ
-
3
-
ρ
-
2
)
,
g
21
=
2
τ
0
D
-
β
(
ρ
-
2
*
(
W
11
(
1
)
(
0
)
ρ
3
+
1
2
W
20
(
1
)
(
0
)
ρ
-
3
iiiiiiiiiiiiiiiiiiiiii
+
W
11
(
3
)
(
0
)
+
1
2
W
20
(
3
)
)
iiiiiiiiiiiiiiii
-
(
W
11
(
1
)
(
0
)
ρ
2
+
1
2
W
20
(
1
)
(
0
)
ρ
-
2
iiiiiiiiiiiiiiiiiiiiiiiii
+
W
11
(
2
)
(
0
)
+
1
2
W
20
(
2
)
)
)
,
with
(32)
W
20
(
θ
)
=
i
g
20
ρ
(
0
)
τ
0
ω
0
e
i
τ
0
ω
0
θ
+
i
g
-
02
ρ
-
(
0
)
3
τ
0
ω
0
e
-
i
τ
0
ω
0
θ
+
E
1
e
2
i
τ
0
ω
0
θ
,
W
11
(
θ
)
=
-
i
g
11
ρ
(
0
)
τ
0
ω
0
e
i
τ
0
ω
0
θ
+
i
g
-
11
ρ
-
(
0
)
τ
0
ω
0
e
-
i
τ
0
ω
0
θ
+
E
2
,
where E
1 and E
2 can be determined by the following equations, respectively,
(33)
(
2
i
ω
0
-
a
1
-
a
2
-
a
3
-
a
4
-
a
5
2
i
ω
0
-
a
6
-
a
7
0
0
-
a
8
2
i
ω
0
-
a
9
-
b
1
e
-
2
i
τ
0
ω
0
0
0
0
-
b
2
e
-
2
i
τ
0
ω
0
2
i
ω
0
-
a
10
)
=
2
(
E
1
(
1
)
E
1
(
2
)
0
0
)
,
(
a
1
a
2
a
3
a
4
a
5
a
6
a
7
0
0
a
8
a
9
+
b
1
0
0
0
b
2
a
10
)
E
2
=
-
(
E
2
(
1
)
E
2
(
2
)
0
0
)
,
with
(34)
E
1
(
1
)
=
-
β
ρ
2
,
E
1
(
2
)
=
β
ρ
2
,
E
2
(
1
)
=
-
β
(
ρ
2
+
ρ
-
2
)
,
E
2
(
2
)
=
β
(
ρ
3
+
ρ
-
3
)
.
Then, we can get the following coefficients:
(35)
C
1
(
0
)
=
i
2
τ
0
ω
0
(
g
11
g
20
-
2
|
g
11
|
2
-
|
g
02
|
2
3
)
+
g
21
2
,
μ
2
=
-
Re
{
C
1
(
0
)
}
Re
{
λ
′
(
τ
0
)
}
,
β
2
=
2
Re
{
C
1
(
0
)
}
,
T
2
=
-
Im
⁡
{
C
1
(
0
)
}
+
μ
2
Im
⁡
{
λ
′
(
τ
0
)
}
τ
0
ω
0
.
In conclusion, we have the following results.Theorem 2.
For system (2), If μ
2
>
0
(
μ
2
<
0
), then the Hopf bifurcation is supercritical (subcritical). If β
2
<
0
(
β
2
>
0
), then the bifurcating periodic solutions are stable (unstable). If T
2
>
0
(
T
2
<
0
), then the period of the bifurcating periodic solutions increases (decreases).


## 4. Numerical Simulation

In this section, we present a numerical example to justify the theoretical results above. We use the same coefficients as those used in [7] and consider the following system:
(36)
d
S
(
t
)
d
t
=
0.001
-
0.15
S
(
t
)
(
L
(
t
)
+
B
(
t
)
)
-
0.001
S
(
t
)
+
0.01
R
(
t
)
,
d
L
(
t
)
d
t
=
0.15
S
(
t
)
(
L
(
t
)
+
B
(
t
)
)
-
0.05
L
(
t
)
-
0.001
L
(
t
)
,
d
B
(
t
)
d
t
=
0.05
L
(
t
)
-
0.02
B
(
t
-
τ
)
-
0.001
B
(
t
)
,
d
R
(
t
)
d
t
=
0.02
B
(
t
-
τ
)
-
0.01
R
(
t
)
-
0.001
R
(
t
)
,
from which we obtain the basic reproduction number R
0
=
9.944
>
1 and the unique positive equilibrium D
*
(
0.1006,0.1160,0.2762,0.5022
) of system (36). Further, we obtain ω
0
=
2.3091, τ
0
=
58.5472. According to Theorem 1 in Section 2, we know that, when τ
∈
[
0,58.5472
), the positive equilibrium D
*
(
0.1006,0.1160,0.2762,0.5022
) of system (36) is locally asymptotically stable. As can be seen from Figures 1 and 2, when τ
=
55.2386
<
τ
0, the positive equilibrium D
*
(
0.1006,0.1160,0.2762,0.5022
) of system (36) is locally asymptotically stable. However, when the delay passes through the critical value of τ
0
=
58.5472, the positive equilibrium D
* will lose stability and a Hopf bifurcation occurs. This property can be illustrated by Figures 3 and 4.Figure 1
The phase plot of statesS, L, and B for τ
=
55.2386
<
τ
0
=
58.5472.Figure 2
The phase plot of statesL, B, and R for τ
=
55.2386
<
τ
0
=
58.5472.Figure 3
The phase plot of statesS, L, and B for τ
=
60.5386
>
τ
0
=
58.5472.Figure 4
The phase plot of statesL, B, and R for τ
=
60.5386
>
τ
0
=
58.5472.In addition, by some complex computations and from (35), we obtain μ
2
=
-
0.8143
<
0, β
2
=
-
0.9001
<
0, and T
2
=
0.0466
>
0. Therefore, we can conclude that the Hopf bifurcation of system (36) is subcritical, the bifurcating periodic solutions are stable, and the period of the periodic solutions decreases according to Theorem 2 in Section 3.


## 5. Conclusions

A delayed SLBRS computer virus model has been studied in this paper based on the model proposed in [7]. Compared with literature [7], we mainly consider the effect of the time delay due to the period that the computers use antivirus software to clean the virus on the model. The main results are presented in terms of local stability and local Hopf bifurcation. Sufficient conditions for local stability and existence of local Hopf bifurcation are obtained by analyzing the distribution of the roots of the associated characteristic equation. It has been shown that when the delay is suitably small (τ
<
τ
0), the computer virus model is locally asymptotically stable. In such case, the propagation of the computer virus can be easily controlled. However, once the delay passes though the critical value τ
0, a Hopf bifurcation occurs and a family of periodic solutions bifurcates from the positive equilibrium of the model. This is not welcome in networks because the computer virus will be out of control in this case. Further, direction and stability of the Hopf bifurcation have also been determined by using the normal form method and center manifold theory. Numerical simulations are also included to verify the analytical predictions.



---

*Source: 373171-2014-08-17.xml*