# Stability and Hopf Bifurcation Analysis of an Epidemic Model by Using the Method of Multiple Scales

Authors: Wanyong Wang; Lijuan Chen
Journal: Mathematical Problems in Engineering
 (2016)

---

## Abstract

A delayed epidemic model with nonlinear incidence rate which depends on the ratio of the numbers of susceptible and infectious individuals is considered. By analyzing the corresponding characteristic equations, the effects of time delay on the stability of the equilibria are studied. By choosing time delay as bifurcation parameter, the critical value of time delay at which a Hopf bifurcation occurs is obtained. In order to derive the normal form of the Hopf bifurcation, an extended method of multiple scales is developed and used. Then, the amplitude of bifurcating periodic solution and the conditions which determine the stability of the bifurcating periodic solution are obtained. The validity of analytical results is shown by their consistency with numerical simulations.


---

## Body

## 1. Introduction

Mathematical models describing the transmission of infectious diseases have played an important role in understanding the mechanism of disease transmission and controlling the spread of infectious diseases. In the literatures, many classical epidemic models have been proposed and studied, and many authors also attempt to develop more realistic mathematical models. In 1927, by using “compartment modeling,” Kermack and Mckendrick [1] described an epidemic model and computed the theoretical number of infectious individuals. From then on, “compartment modeling” is used until now. Compartmental models are also called S
I
R models or S
I
R
S models, where S
,
I, and R denote the number of susceptible, infectious, and recovered individuals. Some S
I
R and S
I
R
S models have been analyzed in detail [2–7]. However, many diseases incubate inside the hosts for a period of time before the hosts become infectious. Therefore, the S
E
I
R epidemic models, where E denotes the number of individuals who are infected but not yet infectious, are developed to investigate the role of the incubation period in disease transmission. Some authors have studied the S
E
I
R epidemic models [8–12]. For some diseases, it is necessary that some infectious individuals are quarantined. Then, the S
I
Q
R and S
E
I
Q
R models, where Q denotes the number of infectious individuals that have been quarantined and controlled, are studied [13–16].Incidence rate plays an important role in the transmission of disease. In earlier models [17], based on the law of mass action, the bilinear incidence rate β
S
I, where β is the average number of contacts per infectious individual per day is used. Capasso and Serio [18] introduced a saturated incidence rate g
(
I
)
S into epidemic models, where g
(
I
)
=
β
I
/
(
1
+
α
I
), where β
I measures the infection force of the disease and 1
/
(
1
+
α
I
) describes the psychological effect or inhibition effect from the behavioral change of the susceptible individuals with the increase of the infectious individuals.Based on the work of Capasso and Serio [18], in this paper, we investigate an S
I
Q
R epidemic model with a saturated incidence rate g
(
·
)
S, where g
(
·
) depends on the ratio of I and S. Therefore, we have (1)
g
=
g
I
S
=
β
I
/
S
1
+
I
/
S
=
β
I
S
+
I
.Then, the saturated incidence rate is β
S
I
/
S
+
I. Considering the effects of time delay, we obtain the following S
I
Q
R epidemic model: (2)
S
˙
=
Λ
-
β
S
t
-
τ
I
t
-
τ
S
t
-
τ
+
I
t
-
τ
-
d
S
,
I
˙
=
β
S
t
-
τ
I
t
-
τ
S
t
-
τ
+
I
t
-
τ
-
c
1
+
k
+
d
+
μ
1
I
,
Q
˙
=
k
I
-
c
2
+
d
+
μ
2
Q
,
R
˙
=
c
1
I
+
c
2
Q
-
d
R
,where S
(
t
), I
(
t
), Q
(
t
), and R
(
t
) denote the number of susceptible, infectious, quarantined, and recovered individuals at time t. Λ is the recruitment rate of the population. d is the natural death rate of the population. μ
1 and μ
2 are the extra disease-related death rate in classes I and Q, respectively. k is the rate constant for individuals leaving infectious compartment I for quarantine compartment Q. c
1 and c
2 are the removal rate constants from classes I and Q, respectively. In this model, we assume that the susceptible individuals were infected before time delay τ which is the latent period. Then, the infected individuals become infectious individuals and some of infectious individuals are quarantined. In classes I and Q, some individuals are cured and removed.The outline of this paper is as follows. In Section2, the equilibria and their stability are analyzed. The critical value of time delay at which Hopf bifurcation occurs is obtained. In Section 3, by using an extended method of multiple scales, we obtained the normal form of the Hopf bifurcation. In Section 4, The above theoretical results are validated by numerical simulations with the help of dynamical software WinPP. In Section 5, conclusions are given.


## 2. Stability of Equilibria and Existence of Hopf Bifurcation

In this paper, we will study system (2) which is modeled by delayed differential equations. Noticing that (
2
1
) and (
2
2
) are uncoupled with (
2
3
) and (
2
4
), we will study the system consisting of (
2
1
) and (
2
2
): (3)
S
˙
=
Λ
-
β
S
t
-
τ
I
t
-
τ
S
t
-
τ
+
I
t
-
τ
-
d
S
,
I
˙
=
β
S
t
-
τ
I
t
-
τ
S
t
-
τ
+
I
t
-
τ
-
c
1
+
k
+
d
+
μ
1
I
.System (3) always has a disease-free equilibrium E
1
=
(
Λ
/
d
,
0
). If β
>
c
1
+
k
+
d
+
μ
1, system (3) has an unique endemic equilibrium E
2
(
S
∗
,
I
∗
), where(4)
S
∗
=
Λ
β
+
d
-
c
1
+
k
+
d
+
μ
1
,
I
∗
=
β
-
c
1
+
k
+
d
+
μ
1
Λ
c
1
+
k
+
d
+
μ
1
β
+
d
-
c
1
+
k
+
d
+
μ
1
.In order to investigate the stability of the two equilibriaE
1 and E
2 simultaneously, we assume that E is one of the two equilibria. Then, the Jacobian matrix of system (3) at point E is given by (5)
-
d
-
p
e
-
λ
τ
-
λ
-
q
e
-
λ
τ
p
e
-
λ
τ
-
f
+
q
e
-
λ
τ
-
λ
=
0
;that is,(6)
λ
2
+
f
+
d
λ
+
p
-
q
λ
e
-
λ
τ
+
f
p
-
q
d
e
-
λ
τ
+
f
d
=
0
,where p
=
β
I
2
/
(
S
+
I
)
2
|
E
,
q
=
β
S
2
/
(
S
+
I
)
2
|
E and f
=
c
1
+
k
+
d
+
μ
1. Equation (6) can be rewritten as (7)
λ
2
+
A
λ
+
B
λ
e
-
λ
τ
+
C
e
-
λ
τ
+
D
=
0
,where A
=
f
+
d
,
B
=
p
-
q
,
C
=
f
p
-
q
d
,
D
=
f
d.Forτ
=
0, the characteristic equation becomes (8)
λ
2
+
A
+
B
λ
+
C
+
D
=
0
.Remark 1.
Forτ
=
0, if C
+
D
>
0 and A
+
B
>
0, then the roots of (8) are real and negative, and then the equilibrium E is asymptotically stable.
Forτ
≠
0, if λ
=
i
ω
(
ω
>
0
) is a root of (7), then we have (9)
-
ω
2
+
A
ω
i
+
B
ω
i
e
-
i
ω
τ
+
C
e
-
i
ω
τ
+
D
=
0
.Separating the real and imaginary parts, we have (10)
B
ω
sin
⁡
ω
τ
+
C
cos
⁡
ω
τ
=
ω
2
-
D
,
B
ω
cos
⁡
ω
τ
-
C
sin
⁡
ω
τ
=
-
ω
A
,which leads to the following fourth-degree polynomial equation: (11)
ω
4
+
A
2
-
B
2
-
2
D
ω
2
+
D
2
-
C
2
=
0
.Then, by discussing the distribution of the solutions of (11), we analyze the stability of the equilibrium and the existence of Hopf bifurcation. The distribution of the solutions of (11) can be divided into the following three cases.Case 1.
If(12)
A
2
-
B
2
-
2
D
>
0
,
D
2
-
C
2
>
0
,or (13)
Δ
=
A
2
-
B
2
-
2
D
2
-
4
D
2
-
C
2
<
0
,then (7) has no positive root.Case 2.
If(14)
D
2
-
C
2
<
0
,then (7) has only one positive root ω
+.Case 3.
If(15)
A
2
-
B
2
-
2
D
<
0
,
D
2
-
C
2
>
0
,
Δ
>
0
,then (7) has two positive roots: (16)
ω
±
=
2
2
-
A
2
-
B
2
-
2
D
±
Δ
,and no such solutions otherwise.Noticing thatω
+ in Case 2 is the same as ω
+ in Case 3, we suppose that (7) has two positive roots ω
±. Then, from (10), we can determine that(17)
τ
j
±
=
1
ω
±
arccos
⁡
C
-
A
B
ω
±
2
-
D
C
B
2
ω
±
2
+
C
2
+
2
j
π
ω
±
,
j
=
0,1
,
2
,
…
,at which (7) has a pair of purely imaginary roots ±
i
ω
±. Supposing that (18)
λ
τ
=
d
τ
+
i
ω
τ
,then (19)
d
τ
j
±
=
0
,
ω
τ
j
±
=
ω
±
.Substitutingλ
(
τ
) into (7) and taking the derivative with respect to τ, we have (20)
d
λ
d
τ
-
1
=
2
λ
+
A
e
λ
τ
B
λ
+
C
λ
+
B
B
λ
+
C
λ
-
τ
λ
,which, together with (10) and (16), leads to (21)
R
e
d
λ
d
τ
τ
=
τ
j
±
-
1
=
R
e
2
λ
+
A
e
λ
τ
B
λ
+
C
λ
τ
=
τ
j
±
+
R
e
B
B
λ
+
C
λ
τ
=
τ
j
±
-
R
e
τ
λ
τ
=
τ
j
±
=
1
B
2
ω
4
+
C
2
ω
2
-
B
ω
2
A
cos
⁡
ω
τ
-
2
ω
sin
⁡
ω
τ
+
C
ω
A
sin
⁡
ω
τ
+
2
ω
cos
⁡
ω
τ
-
B
2
ω
2
=
1
B
2
ω
4
+
C
2
ω
2
-
A
ω
B
ω
cos
⁡
ω
τ
-
C
sin
⁡
ω
τ
+
2
ω
2
B
ω
sin
⁡
ω
τ
+
C
cos
⁡
ω
τ
-
B
2
ω
2
=
1
B
2
ω
4
+
C
2
ω
2
-
A
ω
-
A
ω
+
2
ω
2
ω
2
-
D
-
B
2
ω
2
=
1
B
2
ω
4
+
C
2
ω
2
ω
2
2
ω
2
+
A
2
-
B
2
-
2
D
=
±
Δ
B
2
ω
2
+
C
2
.Thus, if Δ
≠
0, we have (22)
sign
⁡
R
e
d
λ
d
τ
τ
=
τ
j
+
=
sign
⁡
R
e
d
λ
d
τ
τ
=
τ
j
+
-
1
=
sign
⁡
Δ
B
2
ω
2
+
C
2
>
0
,
sign
⁡
R
e
d
λ
d
τ
τ
=
τ
j
-
=
sign
⁡
R
e
d
λ
d
τ
τ
=
τ
j
-
-
1
=
sign
⁡
-
Δ
B
2
ω
2
+
C
2
<
0
.Therefore, we can obtain the following theorem.Theorem 2.
Letω
±, τ
j
±
(
j
=
0,1
,
…
) be defined by (16) and (17), respectively. For system (3),
(i) ifA
2
-
B
2
-
2
D
>
0 and D
2
-
C
2
>
0 or Δ
<
0, then equilibrium E of system (3) is asymptotically stable for all τ
>
0;
(ii) ifD
2
-
C
2
<
0, then equilibrium E of system (3) is asymptotically stable when τ
∈
[
0
,
τ
0
+
) and unstable when τ
>
τ
0
+; system (3) undergoes a Hopf bifurcation at E when τ
=
τ
0
+;
(iii) ifA
2
-
B
2
-
2
D
<
0, D
2
-
C
2
>
0, and Δ
>
0, then there is positive integer k, such that equilibrium E switches k times from stability to instability to stability; that is, E is asymptotically stable when (23)
τ
∈
0
,
τ
0
+
∪
τ
0
-
,
τ
1
+
∪
⋯
∪
τ
k
-
1
-
,
τ
k
+
,and unstable when (24)
τ
∈
τ
0
+
,
τ
0
-
∪
τ
1
+
,
τ
1
-
∪
⋯
∪
τ
k
-
1
+
,
τ
k
-
1
-
,
τ
>
τ
k
+
.System (3) undergoes a Hopf bifurcation at E when τ
=
τ
j
±.


## 3. The Normal Form of Hopf Bifurcation

In this section, we suppose that system (3) undergoes a Hopf bifurcation from equilibrium E
∗
(
S
∗
,
I
∗
). Time delay τ is chosen as the bifurcation parameter and its critical value is τ
c. By using the method of multiple scales [19], we will derive the normal form of the Hopf bifurcation. By translating the equilibrium to the origin, system (3) is rewritten as the following equation: (25)
x
˙
=
F
x
,
x
τ
,where x
=
(
S
(
t
)
,
I
(
t
)
)
T
∈
R
2 is the state variable vector and x
τ
=
(
S
(
t
-
τ
)
,
I
(
t
-
τ
)
)
T.According to the MMS, a monoparametric family of solution of the type is as follows:(26)
x
t
=
ε
x
1
T
0
,
T
2
+
ε
2
x
2
T
0
,
T
2
+
ε
3
x
3
T
0
,
T
2
+
O
ε
4
,where T
k
=
ε
k
t
(
k
=
0,2
) and ε is a nondimensional parameter. The solution does not depend on slow scale T
1
=
ε
t because secular terms first appear at O
(
ε
3
). Therefore, we assume a two-scale expansion of the solution of (25).Bifurcation parameterτ is ordered as (27)
τ
=
τ
c
+
ε
2
τ
ε
.The delay term in (25) can be further expanded as (28)
x
t
-
τ
=
ε
x
1
T
0
-
τ
c
-
ε
2
τ
ε
,
T
2
-
ε
2
τ
c
-
ε
4
τ
ε
+
ε
2
x
2
T
0
-
τ
c
-
ε
2
τ
ε
,
T
2
-
ε
2
τ
c
-
ε
4
τ
ε
+
ε
3
x
3
T
0
-
τ
c
-
ε
2
τ
ε
,
T
2
-
ε
2
τ
c
-
ε
4
τ
ε
+
⋯
=
ε
x
1
T
0
-
τ
c
,
T
2
+
ε
2
x
2
T
0
-
τ
c
,
T
2
+
ε
3
-
τ
ε
D
0
x
1
T
0
-
τ
c
,
T
2
-
τ
c
D
2
x
1
T
0
-
τ
c
,
T
2
+
x
3
T
0
-
τ
c
,
T
2
+
O
ε
4
,where D
k
=
∂
/
∂
T
k.By substituting (26), (27), and (28) into (25), expanding F, and balancing the coefficients of ε
j
(
j
=
1,2
,
3
.
), a set of perturbative equations are obtained.First, forε-order terms, we have (29)
D
0
x
1
-
F
x
0
x
1
-
F
x
τ
0
x
1
τ
=
0
,where (30)
F
x
0
=
-
d
0
0
-
f
,
F
x
τ
0
=
-
p
-
q
p
q
.Equation (29) has the following general solution: (31)
x
1
=
A
T
2
u
e
i
ω
T
0
+
A
-
T
2
u
-
e
-
i
ω
T
0
,where A is complex constant and u is the right eigenvector given by (32)
F
x
0
+
F
x
τ
0
e
-
i
ω
τ
c
-
i
E
ω
u
=
0
,then (33)
u
=
1
-
d
+
p
e
-
i
ω
τ
+
i
ω
q
e
-
i
ω
τ
.Next, for ε
2-order terms, we have (34)
D
0
x
2
-
F
x
0
x
2
-
F
x
τ
0
x
2
τ
=
1
2
F
x
x
0
x
1
2
+
F
x
x
τ
0
x
1
x
1
τ
+
1
2
F
x
τ
x
τ
0
x
1
τ
2
,where F
x
x
0
,
F
x
x
τ
0
,
F
x
τ
x
τ
0 are the second-order partial derivatives of F
(
x
,
x
τ
) with respect to x and x
τ.Substituting (31) into (34), we obtain (35)
D
0
x
2
-
F
x
0
x
2
-
F
x
τ
0
x
2
τ
=
1
2
F
x
x
0
u
2
+
F
x
x
τ
0
u
2
e
-
i
ω
τ
c
+
1
2
F
x
τ
x
τ
0
u
2
e
-
2
i
ω
τ
c
·
A
2
e
2
i
ω
T
0
+
1
2
F
x
x
0
u
u
-
+
F
x
x
τ
0
u
u
-
e
i
ω
τ
c
+
1
2
F
x
τ
x
τ
0
u
u
-
·
A
A
-
+
c
.
c
..Solving (35), it yields (36)
x
2
=
A
2
z
11
e
2
i
ω
T
0
+
A
A
-
z
1
1
-
+
c
.
c
.,where vectors z
11 and z
1
1
-
∈
C
2 are obtained by solving the following equations: (37)
2
i
ω
E
-
F
x
0
-
F
x
τ
0
e
-
2
i
ω
τ
c
z
11
=
1
2
F
x
x
0
u
2
+
F
x
x
τ
0
u
2
e
-
i
ω
τ
c
+
1
2
F
x
τ
x
τ
0
u
2
e
-
2
i
ω
τ
c
,
-
F
x
0
-
F
x
τ
0
z
1
1
-
=
1
2
F
x
x
0
u
u
-
+
F
x
x
τ
0
u
u
-
e
i
ω
τ
c
+
1
2
F
x
τ
x
τ
0
u
u
-
;that is, (38)
2
i
ω
+
d
+
p
e
-
2
i
ω
τ
c
q
e
-
2
i
ω
τ
c
-
p
e
-
2
i
ω
τ
c
2
i
ω
+
f
-
q
e
-
2
i
ω
τ
c
z
111
z
112
=
1
2
·
e
-
2
i
ω
τ
c
f
20
1
u
1
2
+
f
11
1
u
1
u
2
+
f
02
1
u
2
2
f
20
2
u
1
2
+
f
11
2
u
1
u
2
+
f
02
2
u
2
2
,
d
+
p
q
-
p
f
-
q
z
1
1
-
1
z
1
1
-
2
=
1
2
·
f
20
1
u
1
u
-
1
+
f
11
1
1
2
u
1
u
-
2
+
1
2
u
2
u
-
1
+
f
02
1
u
2
u
-
2
f
20
2
u
1
u
-
1
+
f
11
2
1
2
u
1
u
-
2
+
1
2
u
2
u
-
1
+
f
02
2
u
2
u
-
2
,where f
i
j
(
1
)
=
∂
i
+
j
f
1
/
∂
S
(
t
-
τ
)
i
∂
I
(
t
-
τ
)
j
|
(
S
,
I
), f
i
j
(
2
)
=
∂
i
+
j
f
2
/
∂
S
(
t
-
τ
)
i
∂
I
(
t
-
τ
)
j
|
(
S
,
I
), f
(
1
)
=
Λ
-
β
S
(
t
-
τ
)
I
(
t
-
τ
)
/
S
(
t
-
τ
)
+
I
(
t
-
τ
)
-
d
S, f
(
2
)
=
β
S
(
t
-
τ
)
I
(
t
-
τ
)
/
S
(
t
-
τ
)
+
I
(
t
-
τ
) − (
c
1
+
k
+
d
+
μ
1
)
I. Then, we can obtain (39)
z
111
=
1
W
·
1
2
f
20
1
u
1
2
+
f
11
1
u
1
u
2
+
f
02
1
u
2
2
e
-
2
i
ω
τ
c
q
e
-
2
i
ω
τ
c
1
2
f
20
2
u
1
2
+
f
11
2
u
1
u
2
+
f
02
2
u
2
2
e
-
2
i
ω
τ
c
2
i
ω
+
f
-
q
e
-
2
i
ω
τ
c
,
z
112
=
1
W
·
2
i
ω
+
d
+
p
e
-
2
i
ω
τ
c
1
2
f
20
1
u
1
2
+
f
11
1
u
1
u
2
+
f
02
1
u
2
2
e
-
2
i
ω
τ
c
-
p
e
-
2
i
ω
τ
c
1
2
f
20
2
u
1
2
+
f
11
2
u
1
u
2
+
f
02
2
u
2
2
e
-
2
i
ω
τ
c
,
W
=
2
i
ω
+
d
+
p
e
-
2
i
ω
τ
c
q
e
-
2
i
ω
τ
c
-
p
e
-
2
i
ω
τ
c
2
i
ω
+
f
-
q
e
-
2
i
ω
τ
c
,
z
1
1
-
1
=
1
d
f
+
p
f
-
d
q
·
1
2
f
20
1
u
1
u
-
1
+
f
11
1
1
2
u
1
u
-
2
+
1
2
u
2
u
-
1
+
f
02
1
u
2
u
-
2
q
1
2
f
20
2
u
1
u
-
1
+
f
11
2
1
2
u
1
u
-
2
+
1
2
u
2
u
-
1
+
f
02
2
u
2
u
-
2
f
-
q
,
z
1
1
-
2
=
1
d
f
+
p
f
-
d
q
·
d
+
p
1
2
f
20
1
u
1
u
-
1
+
f
11
1
1
2
u
1
u
-
2
+
1
2
u
2
u
-
1
+
f
02
1
u
2
u
-
2
-
p
1
2
f
20
2
u
1
u
-
1
+
f
11
2
1
2
u
1
u
-
2
+
1
2
u
2
u
-
1
+
f
02
2
u
2
u
-
2
.Finally, for ε
3-order terms, we get (40)
D
0
x
3
-
F
x
0
x
3
-
F
x
τ
0
x
3
τ
=
1
6
F
x
x
x
0
x
1
3
+
F
x
x
0
x
1
x
2
+
1
2
F
x
x
x
τ
0
x
1
2
x
1
τ
+
1
2
F
x
x
τ
x
τ
0
x
1
x
1
τ
2
+
F
x
x
τ
0
x
1
τ
x
2
+
1
6
F
x
τ
x
τ
x
τ
0
x
1
τ
3
+
F
x
x
τ
0
x
1
x
2
τ
+
F
x
τ
x
τ
0
x
1
τ
x
2
τ
-
τ
c
F
x
τ
0
D
2
x
1
τ
-
τ
ε
F
x
τ
0
D
0
x
1
τ
-
D
2
x
1
,where F
x
x
x
0
,
F
x
x
x
τ
0
,
F
x
x
τ
x
τ
0
,
F
x
τ
x
τ
x
τ
0 are the third-order partial derivatives of F
(
x
,
x
τ
) with respect to x and x
τ.Substituting (31) and (36) into (40) and eliminating the secular terms, we can obtain the equations including D
2
A. Eliminating the coefficients of D
2
A by using the left eigenvectors and reabsorbing parameter ε [20], the normal form is determined by(41)
A
˙
=
C
τ
τ
ε
A
+
C
11
1
-
A
2
A
-
,where the expressions of coefficients C
τ
τ
ε and C
11
1
- are given by(42)
C
τ
τ
ε
=
1
M
v
H
-
i
ω
τ
ε
F
x
τ
0
u
e
-
i
ω
τ
c
,
C
11
1
-
=
1
M
v
H
1
2
F
x
x
x
0
u
2
u
-
+
F
x
x
0
u
z
-
1
1
-
+
z
11
u
-
+
z
1
1
-
u
+
1
2
F
x
x
x
τ
0
2
u
2
u
-
e
-
i
ω
τ
c
+
u
2
u
-
e
i
ω
τ
c
+
F
x
x
τ
0
u
z
-
1
1
-
e
-
i
ω
τ
c
+
z
11
u
-
e
i
ω
τ
c
+
z
1
1
-
u
e
-
i
ω
τ
c
+
1
2
F
x
x
τ
x
τ
0
2
u
2
u
-
+
u
2
u
-
e
-
2
i
ω
τ
c
+
1
2
F
x
τ
x
τ
x
τ
0
u
2
u
-
e
-
i
ω
τ
c
+
F
x
x
τ
0
u
z
-
1
1
-
+
z
11
u
-
e
-
2
i
ω
τ
c
+
z
1
1
-
u
+
F
x
τ
x
τ
0
u
z
-
1
1
-
e
-
i
ω
τ
c
+
z
11
u
-
e
-
i
ω
τ
c
+
z
1
1
-
u
e
-
i
ω
τ
c
,where M
=
v
H
(
τ
c
F
x
τ
0
e
-
i
ω
τ
c
+
E
)
u and v is the left eigenvector obtained by solving the following equations: (43)
F
x
0
+
F
x
τ
0
e
i
ω
τ
c
T
+
i
E
ω
v
=
0
,
v
H
u
=
1
,where H denotes the transpose conjugate and T denotes the transpose.To express the normal form in real form, a polar form representation for the complex amplitude is introduced:(44)
A
=
1
2
a
e
i
θ
.Substituting (44) into (41) and separating the real and imaginary parts in (41), the generalized amplitude and phase modulation equations are drawn: (45)
a
˙
=
a
R
1
+
1
4
R
11
1
-
a
3
,
θ
˙
=
I
1
+
1
4
I
11
1
-
a
2
,where R
1
=
R
e
(
C
τ
τ
ε
)
,
I
1
=
Im
⁡
(
C
τ
τ
ε
)
,
R
11
1
-
=
R
e
(
C
11
1
-
), and I
11
1
-
=
Im
⁡
(
C
11
1
-
).From (45), we can get that the amplitude of the periodic solution is (46)
a
=
2
-
R
1
R
11
1
-
.The stability of the periodic solution is determined by the sign of the eigenvalue of the Jacobian matrix of (
45
1
). The eigenvalue is (47)
λ
=
-
2
R
1
.Then, we can obtain the following theorem.Theorem 3.
The amplitude of the bifurcating periodic solution of system (3) is 2
-
R
1
/
R
11
1
-; the stability of the bifurcating periodic solution is determined by λ: the bifurcating periodic solution is stable (unstable), if λ
<
0
(
λ
>
0
).


## 4. An Example

In this section, we perform some numerical simulations to verify our analysis. For the following parameter valuesΛ
=
10
,
β
=
0.6
,
d
=
0.02
,
c
1
=
0.08
,
k
=
0.01
,
μ
1
=
0.03, we can see that endemic equilibria E
∗
(
20.8333,68.4524
) of system (3) exist. Then, a periodic solution with frequency ω
=
0.32395 bifurcates from the endemic equilibria at critical value τ
=
τ
c
=
4.94223. By Theorem 2, endemic equilibria E
∗ are asymptotically stable when τ
∈
[
0
,
τ
c
) and unstable when τ
>
τ
c.By using the dynamical software WinPP, some numerical simulations are given. First, by Theorem2, endemic equilibria E
∗ are asymptotically stable when τ
∈
[
0
,
τ
c
) as shown in Figure 1 and the time histories of the state variables Q and R are also given in Figure 1.Figure 1
The response of state variables of system (2) with Λ
=
10
,
β
=
0.6
,
d
=
0.02
,
c
1
=
0.08
,
k
=
0.01
,
μ
1
=
0.03
,
c
2
=
0.02
,
μ
2
=
0.04, and τ
=
4.8
<
τ
c
=
4.94223.Endemic equilibriaE
∗ are unstable when τ
>
τ
c. When τ crosses critical value τ
c, a stable periodic solution yields as shown in Figure 2.Figure 2
The response of state variables of system (2) with Λ
=
10
,
β
=
0.6
,
d
=
0.02
,
c
1
=
0.08
,
k
=
0.01
,
μ
1
=
0.03
,
c
2
=
0.02
,
μ
2
=
0.04, and τ
=
5.1
>
τ
c
=
4.94223.By using the method developed in this paper, the amplitude of the periodic solution can be obtained. The relationship between the amplitude and bifurcation parameterτ is given in (45). Then, the analysis results are compared with that of numerical simulations as shown in Figure 3.Figure 3
Comparison between the theoretical solution (solid) and numerical simulation (dot) for the amplitude of periodic solution.


## 5. Conclusions

In this paper, the stability and Hopf bifurcation of an epidemic model with time delay are studied. The stability of disease-free equilibrium and endemic equilibrium is studied and the conditions of the stability and unstability of the two equilibria are obtained. Then, the critical value ofτ at which a Hopf bifurcation occurs is obtained. By using an extended method of multiple scales, the normal form of Hopf bifurcation is obtained. In order to verify the theoretical results, the Runge-Kutta scheme is adopted to produce the numerical results. The comparison between analytical predictions and numerical results of the amplitude of Hopf bifurcation reveals a qualitatively and quantificationally excellent agreement.In order to obtain the normal form of Hopf bifurcation occurring in delayed differential equations, we develop an extended method of multiple scales. Choosing time delay as a bifurcation parameter and using the method, one can get the normal form of Hopf bifurcation without the tedious computation encountered in the center manifold reduction. Moreover, when we develop the method, we choosen-dimensional delayed differential equation. Therefore, the method can be used for studying the Hopf bifurcation occurring in a general n-dimensional delayed differential equation.Through observing the results, there are two contributions in this paper. First, we improve an epidemic model and introduce a nonlinear incidence rate which depends on the ratio of the numbers of susceptible and infectious individuals. Second, we extend the method of multiple scales. By using this method, the normal form of Hopf bifurcation occurring inn-dimensional delayed differential equations can be obtained.



---

*Source: 2034136-2016-08-30.xml*