# Mathematical Analysis of HIV Models with Switching Nonlinear Incidence Functions and Pulse Control

Authors: Xiying Wang; Wei Xu; Yujun Cui; Xiaomei Wang
Journal: Abstract and Applied Analysis
 (2014)

---

## Abstract

This paper aims to study the dynamics of new HIV (the human immunodeficiency virus) models with switching nonlinear incidence functions and pulse control. Nonlinear incidence functions are first assumed to be time-varying functions and switching functional forms in time, which have more realistic significance to model infectious disease models. New threshold conditions with the periodic switching term are obtained to guarantee eradication of the disease, by using the novel type of common Lyapunov function. Furthermore, pulse vaccination is applied to the above model, and new sufficient conditions for the eradication of the disease are presented in terms of the pulse effect and the switching effect. Finally, several numerical examples are given to show the effectiveness of the proposed results, and future directions are put forward.


---

## Body

## 1. Introduction

Dynamical behavior of HIV infection models has been investigated to explain different phenomena with the help of the persistence of the disease and the global stability of the disease-free equilibrium. D’Onofrio [1] studied the global asymptotic stability of the disease-free equilibrium of the HIV infection model. It has been shown that the drug efficacy functions are bang-bang type, and the stability of the infection free steady state was studied by the basic reproduction number R
=
(
1
-
e
1
)
(
1
-
e
2
)
N
β
S
/
a
b with γ
=
e
1 and η
=
e
2 [2].Assume that the parameterS is the total rate of production of healthy cells per unit time; b, a, and c are the per capita death rate of healthy cells, infected cells, and infective virus particles, respectively. Then, a basic mathematical model of HIV dynamics [1, 3–6], consisting of three state variables at time t corresponding to concentration of uninfected target cells T
(
t
), infected cells I
(
t
), and free virus particles V
(
t
), may be described by the following:
(1)
T
˙
=
S
-
b
T
-
f
(
T
,
V
)
,
I
˙
=
f
(
T
,
V
)
-
a
I
,
V
˙
=
g
(
I
)
-
c
V
,
where f
(
T
,
V
) and g
(
I
) denote incidence rate functions which are the average number of new infected cells and new virus particles per unite time, respectively. In infectious disease modeling, the incidence functions have become a crucial factor to ensure that the models have some realistic significance and may give some reasonable description. For example, in [7–9], authors assumed that f
(
T
,
V
)
=
β
T
V and g
(
I
)
=
K
I, where β is the transmission coefficient between uninfected cells and infective virus particles, and K is the average number of infective virus particles produced by an infected cell. Under the influence of HAART (highly active antiviral therapy), which is generally a combination of RTI (reverse transcriptase inhibitors) drugs and PI (protease inhibitors) drugs, it is assumed that f
(
T
,
V
)
=
(
1
-
γ
)
β
T
V and g
(
I
)
=
(
1
-
η
)
N
a
I [1, 2], where N is the average number of infective virus particles produced by an infected cell in the absence of HAART; 1
-
γ
(
0
<
γ
<
1
) is the reverse transcriptase inhibitor drug effect; and 1
-
η
(
0
<
η
<
1
) is the protease inhibitor drug effect. Here, f
(
T
,
V
) is bilinear with respect to the number of uninfected cells and virus particles, and g
(
I
) is linear with respect to the number of infected cells. In fact, both bilinear incidence rate and linear incidence rate may be modified.Recently, most of epidemic infection models with nonlinear incidence functions have been investigated by many authors. Korobeinikov [10] established global properties of two-dimension SIR and SIRS compartmental epidemic models with nonlinear transmission rate based on the method of Lyapunov functions. Liu and Stechlinski [11] analyzed infectious disease models with time-varying parameters and general nonlinear incidence rates and obtained some sufficient conditions to ensure the stability of the disease-free equilibrium. If the population is saturated with the infected individuals, the incidence rate may have nonlinear dependence on infective individuals [12]. In general, the HIV models’ parameters (e.g., the contact rate, the effect of RTI drugs, and PI drugs) were assumed to be constant in time. However, these parameters may be time-varying, due to changes in host behavior. Take the effect of RTI drugs and PI drugs; for example, the effect of RTI drugs and PI drugs is usually characterized by a quick rise to a maximum soon after drug intake, followed by a slower decay within a cycle. A more realistic approach is to assume that these parameters are time-varying, which implies that two nonlinear incidence functions are time-varying functions.A switched epidemic model is modeled by introducing switching functions into the HIV model. Moreover, according to the method of [3, 13], two switching nonlinear incidence functions are introduced into system (1) by replacing general nonlinear incidence functions.Switched systems, consisting of continuous, discrete dynamics and logic based switching rule, have gained considerable attention by authors [14–18]. One main feature of the switched system is that the included switching law may induce stability of the switched system composed of two unstable subsystems. Switched systems have been applied in various areas, such as engine control systems, neural networks, ecosystems, mechanical systems, and even biological systems. Stability results have been derived by many researchers, via the methods of Lyapunov exponents, switched Lyapunov functions, and common Lyapunov functions (see [19–23]). Until now, there are few works about the switched HIV models. This paper mainly analyzes the HIV models subjected to switching nonlinear incidence functions. Using common Lyapunov functions method, the global stability of the disease-free equilibrium is discussed and new stability criteria are established to ensure eradication of the disease.Moreover, pulse control is applied into the HIV models with switching nonlinear incidence functions. Due to switches of states and abrupt changes at the switching instants, switched systems exhibit impulsive effects, and they cannot be well described by using pure continuous or discrete switched systems [24]. Some results on the impulsive systems have been obtained [25, 26]. In this paper, new HIV models with switching nonlinear incidence functions and pulse control are developed, and the global asymptotic stability by the technique of common Lyapunov functions is analyzed.The paper is organized as follows. The HIV model with switching nonlinear incidence functions is introduced, and global asymptotic stability of the disease-free equilibrium is presented in Section2. In Section 3, pulse control is considered in the above HIV model, and sufficient conditions for the global asymptotic stability are obtained by the method of common Lyapunov functions. Numerical simulations are given in Section 4 to illustrate the threshold conditions established in the paper. Some conclusions and future directions are given in Section 5.


## 2. The HIV Model with Switching Nonlinear Incidence Functions

In general nonlinear incidence functions [3], the HIV models’ parameters (e.g., the contact rate, the effect of RTI drugs, and PI drugs) are constants in time. In fact, the HIV models’ parameters may be time-varying, due to changes in host behavior [11, 19]. Taking the contact rate, for example, it may be time-varying and switch in time under the influence of environmental factors or changes in host behavior. According to the above approach, assume that incidence functions are time-varying functions and may change functional forms in time. Consider families of m different incidence functions f
i
(
T
,
V
) and g
i
(
I
), for i
=
1,2
,
…
,
m. Assume that the average number of new infected cells and new virus per unit time are modeled by the switching functions f
i
k
(
T
,
V
) and g
i
k
(
I
) on the interval (
t
k
-
1
,
t
k
], k
=
1,2
,
…. Suppose that these switching functions are governed by a switching rule σ
(
t
): (
t
k
-
1
,
t
k
]
→
i
k
∈
{
1,2
,
…
,
m
}, k
=
1,2
,
…, which is a piecewise constant function of time, continuous from the left. Assume that the switch time {
t
k
} satisfies t
k
>
t
k
-
1 and t
k
→
∞ as k
→
∞. Denote set of all switching rules by I. This leads to a new HIV model with switching nonlinear incidence functions:
(2)
T
˙
=
S
-
b
T
-
f
i
k
(
T
,
V
)
,
t
∈
(
t
k
-
1
,
t
k
]
,
I
˙
=
f
i
k
(
T
,
V
)
-
a
I
,
V
˙
=
g
i
k
(
I
)
-
c
V
.
The initial conditions are T
(
0
)
=
T
0, I
(
0
)
=
I
0, and V
(
0
)
=
V
0. Here T, I, and V
∈
R
+, where R
+ denotes the set of nonnegative real numbers, and S, b, a, and c are in R
+. From physical considerations, assume that f
i
(
T
,
V
) and g
i
(
I
) are continuous functions, which satisfy f
i
(
T
,
V
)
>
0 and g
i
(
I
)
>
0, for all t
≥
t
0, i
=
1,2
,
…
,
m, except along the boundaries f
i
(
T
,
0
)
=
f
i
(
0
,
V
)
=
0 and g
i
(
0
)
=
0. Then, (2) has a disease-free equilibrium E
0
=
(
S
/
b
,
0,0
) and m endemic solutions E
i
*
(
T
i
*
,
I
i
*
,
V
i
*
) for i
=
1,2
,
…
,
m, where T
i
*, I
i
*, V
i
* satisfy the following:
(3)
S
=
b
T
i
*
+
a
I
i
*
,
a
I
i
*
=
f
i
(
V
i
*
)
T
i
*
,
c
V
i
*
=
g
i
(
I
i
*
)
.
HIV drugs are most commonly prescribed to be taken on a fixed dose, fixed time interval basis. In fact, the fixed time interval between two dose can be seen as periodic and denoted by ω. Suppose that the switching rule σ satisfies t
k
-
t
k
-
1
=
τ
k with τ
k
+
m
=
τ
k, for t
∈
(
t
k
-
1
,
t
k
]. Suppose that f
i
k
(
t
)
=
f
k
(
t
) and g
i
k
(
t
)
=
g
k
(
t
). Suppose that f
k
(
t
)
=
f
k
(
t
+
ω
) and g
k
(
t
)
=
g
k
(
t
+
ω
) and β
k
(
t
)
=
β
k
(
t
+
ω
), where ω
=
τ
1
+
τ
2
+
⋯
+
τ
m. Let I
Periodic be the set of periodic switching rule with I
Periodic
⊂
I.Introducing some variables,X
1
=
T
-
S
/
b, X
2
=
I, and X
3
=
V, and substituting the corresponding variables into system (2), we get
(4)
X
˙
1
=
-
b
X
1
-
f
i
k
(
X
1
+
S
b
,
X
3
)
,
t
∈
(
t
k
-
1
,
t
k
]
,
X
˙
2
=
f
i
k
(
X
1
+
S
b
,
X
3
)
-
a
X
2
,
X
˙
3
=
g
i
k
(
X
2
)
-
c
X
3
.
For convenience of analysis, rewrite (4) in the following form:
(5)
X
˙
=
A
X
+
F
i
k
(
X
)
,
t
∈
(
t
k
-
1
,
t
k
]
,
where X
=
(
X
1
,
X
2
,
X
3
)
′, ′ denotes transposition, A is a 3
×
3 square matrix, and F
i
k is a column vector given by
(6)
A
=
(
-
b
0
0
0
-
a
0
0
0
-
c
)
,
F
i
k
=
(
-
f
i
k
(
X
1
+
S
b
,
X
3
)
f
i
k
(
X
1
+
S
b
,
X
3
)
g
i
k
(
X
2
)
.
)
.
Assume that switching nonlinear functions f
i
(
X
) and g
i
(
X
) for i
=
1,2
,
…
,
m satisfy locally Lipschitz conditions; that is, for each X
∈
R
3
+, there exist functions L
i
≥
0, R
i
≥
0, and δ
=
δ
(
X
)
>
0, such that ∥
X
-
Y
∥
<
δ implies that
(7)
∥
f
i
(
X
)
-
f
i
(
Y
)
∥
≤
L
i
∥
X
-
Y
∥
,
∥
g
i
(
X
)
-
g
i
(
Y
)
∥
≤
R
i
∥
X
-
Y
∥
,
in which ∥
X
∥
=
(
X
1
2
+
X
2
2
+
X
3
2
)
1
/
2.Obviously, due to switches, the characters of system (2) are different from most existing models (see [1, 8] and the reference therein). It is necessary to investigate the above HIV model. By transformation, the dynamics of disease-free equilibrium of system (2) are the same as the trivial solution of system (5). Via common Lyapunov functions, we consider the global asymptotic stability of the trivial solution of system (5). The idea is that we first find a common Lyapunov function for each subsystem and then impose restrictions on switching to guarantee the stability of system (5). We have the following result regarding the global asymptotic stability of the trivial solution of system (5).Theorem 1.
Assume that the switching ruleσ is periodic and that f
i
(
t
,
X
)
=
f
i
(
t
+
ω
,
X
) and g
i
(
t
,
X
)
=
g
i
(
t
+
ω
,
X
), for i
=
1,2
,
…
,
m, where f
i
(
X
) and g
i
(
X
) satisfy (7). If there exists a constant C
>
0, such that
(8)
τ
1
λ
1
+
τ
2
λ
2
+
⋯
+
τ
m
λ
m
<
-
C
,
where λ
i
=
λ
+
2
L
i
2
+
R
i
2
+
1 with λ
=
max
⁡
{
-
2
a
,
-
2
b
,
-
2
c
}, then, the trivial solution of system (5) is globally asymptotically stable, which implies that the disease-free equilibrium E
0 of system (2) is globally asymptotically stable.Proof.
Assume thati
=
i
k, for t
∈
(
t
k
-
1
,
t
k
], where i
k follows the switching rule σ
(
t
)
∈
I
Periodic. Define a common Lyapunov function U as U
=
X
′
X. Then, we have
(9)
d
U
(
t
)
d
t
=
X
˙
′
X
+
X
′
X
˙
=
X
′
(
A
′
+
A
)
X
+
2
F
i
k
′
X
.
Note that
(10)
X
′
(
A
′
+
A
)
X
≤
λ
X
′
X
,
2
F
i
k
′
X
≤
2
∥
f
i
k
(
t
,
X
)
∥
2
+
∥
g
i
k
(
t
,
X
)
∥
2
+
X
′
X
≤
(
2
L
i
k
2
+
R
i
k
2
+
1
)
X
′
X
.
Substituting (10) into (9) yields
(11)
U
˙
≤
(
λ
+
2
L
i
k
2
+
R
i
k
2
+
1
)
X
′
X
=
λ
i
k
U
.
For t
∈
(
t
k
-
1
,
t
k
], it follows that
(12)
U
(
t
)
≤
U
(
t
k
-
1
)
exp
⁡
[
∫
t
k
-
1
t
λ
i
k
(
s
)
d
s
]
.
Apply (12) on each subinterval. For t
∈
(
t
0
,
t
1
], U
(
t
)
≤
U
(
t
0
)
exp
⁡
[
∫
t
0
t
λ
i
1
(
s
)
d
s
]. For t
∈
(
t
1
,
t
2
], U
(
t
)
≤
U
(
t
1
)
exp
⁡
[
∫
t
1
t
λ
i
2
(
s
)
d
s
]
≤
U
(
t
0
)
exp
⁡
[
∫
t
0
t
1
λ
i
1
(
s
)
d
s
+
∫
t
1
t
λ
i
2
(
s
)
d
s
]. In general, for t
∈
(
t
k
-
1
,
t
k
], U
(
t
)
≤
U
(
t
k
-
1
)
exp
⁡
[
∫
t
k
-
1
t
‍
λ
i
k
(
s
)
d
s
]
≤
U
(
t
0
)
exp
⁡
[
∫
t
0
t
1
‍
λ
i
1
(
s
)
d
s
+
∫
t
1
t
2
‍
λ
i
2
(
s
)
d
s
+
⋯
+
∫
t
k
-
1
t
‍
λ
i
k
(
s
)
d
s
]. For t
∈
(
ω
-
t
m
,
ω
], U
(
t
)
≤
U
(
t
0
)
exp
⁡
[
∫
t
0
t
1
‍
λ
i
1
(
s
)
d
s
+
∫
t
1
t
2
‍
λ
i
2
(
s
)
d
s
+
⋯
+
∫
ω
-
t
m
t
‍
λ
i
m
(
s
)
d
s
]. Then, U
(
ω
)
≤
U
(
t
0
)
exp
⁡
[
τ
1
λ
1
+
τ
2
λ
2
+
⋯
+
τ
m
λ
m
]. Define ζ
=
exp
⁡
[
τ
1
λ
1
+
τ
2
λ
2
+
⋯
+
τ
m
λ
m
]. By condition (8), we have ζ
<
1; furthermore, U
(
ω
)
≤
ζ
U
(
t
0
)
<
U
(
t
0
). Similarly, it can be shown that U
(
h
ω
)
≤
ζ
U
(
(
h
-
1
)
ω
) for any integer h
=
1,2
,
… and hence
(13)
U
(
h
ω
)
≤
ζ
U
(
(
h
-
1
)
ω
)
≤
ζ
(
ζ
U
(
(
h
-
2
)
ω
)
)
≤
⋯
≤
ζ
h
U
(
t
0
)
;
it follows that lim
⁡
h
→
∞
⁡
U
(
h
ω
)
=
lim
⁡
h
→
∞
⁡
ζ
h
U
(
t
0
)
=
0. In general, for t
∈
(
t
k
-
1
,
t
k
] and h
ω
<
t
k
≤
(
h
+
1
)
ω, we have
(14)
U
(
t
)
≤
U
(
h
ω
)
exp
⁡
[
∫
t
0
t
1
λ
i
1
(
s
)
d
s
+
∫
t
1
t
2
λ
i
2
(
s
)
d
s
hhhhhhhhhhh
+
⋯
+
∫
t
k
-
1
t
λ
i
k
(
s
)
d
s
]
≤
M
U
(
h
ω
)
≤
M
U
(
t
0
)
ζ
h
=
M
U
(
t
0
)
exp
⁡
[
h
ln
⁡
ζ
]
,
where M
=
max
⁡
t
k
-
1
<
t
≤
t
k
⁡
exp
⁡
[
∫
t
0
t
1
‍
λ
i
1
(
s
)
d
s
+
∫
t
1
t
2
‍
λ
i
2
(
s
)
d
s
+
⋯
+
∫
t
k
-
1
t
‍
λ
i
k
(
s
)
d
s
]. Note that ln
⁡
ζ
<
0; then, U
(
t
) converges to zero as t
→
∞ and h
→
∞. Therefore, the trivial solution of system (5) is globally asymptotically stable, which implies that the disease-free equilibrium E
0 of system (2) is globally asymptotically stable.Remark 2.
A general Lyapunov function is not used to deal with switched differential equations. In order to display the effect of switching, we consider the dynamics of the switched HIV models via the method of the common Lyapunov function.Remark 3.
In Theorem1, general criteria are established to ensure that the disease will die out, no matter whether the subsystems are stable or unstable. Compared with [13], the stability of the disease-free equilibrium of system (2) is investigated by the common Lyapunov functions method.Remark 4.
Assume thatf
i
(
T
,
V
)
≡
(
1
-
γ
)
β
T
V and g
i
(
I
)
≡
N
a
I; system (2) reduces to the system (1.1) in [5]. Compared with the results without the switching effect of [5], the results here are closer to the reality with practical significance, characterizing the switching effect τ
i
λ
i in Theorem 1.


## 3. The Switched HIV Model with Pulse Vaccination

In this section, we investigate the dynamics of HIV model with switching nonlinear incidence functions and pulse control by common Lyapunov functions. Pulse control is strategy of periodically vaccinating the infectious disease in a relatively short time [11, 26]. Assume that a fraction p
(
0
<
p
<
1
) of infected cells is impulsively treated every ω
>
0 time unit, moving infected cells q
I
(
t
)
(
0
≤
q
<
p
<
1
) to the various classes. This is reasonable from a physical perspective, since some of infected cells are affected badly by pulse vaccination, and they are becoming various particles. If q
=
0 and p
>
0, then the vaccine has efficacy. If p
=
q
=
1, then vaccine is failing completely. Applying pulse vaccination to system (2), this gives a new switched HIV model
(15)
T
˙
=
S
-
b
T
-
f
i
k
(
T
,
V
)
,
t
∈
(
t
k
-
1
,
t
k
]
,
I
˙
=
f
i
k
(
T
,
V
)
-
a
I
,
V
˙
=
g
i
k
(
I
)
-
c
V
,
T
(
t
+
)
=
T
(
t
)
,
t
=
t
k
,
I
(
t
+
)
=
(
1
-
p
)
I
(
t
)
,
V
(
t
+
)
=
V
(
t
)
+
q
I
(
t
)
.
The initial conditions are T
(
t
0
+
)
=
T
0
>
0, I
(
t
0
+
)
=
I
0
≥
0, and V
(
t
0
+
)
=
V
0
≥
0. Note that system (15) has the same disease-free equilibrium E
0 as that of system (2). Under the transformation X
1
=
T
-
S
/
b, X
2
=
I, X
3
=
V, system (15) has the following form:
(16)
X
˙
1
=
-
b
X
1
-
f
i
k
(
X
1
+
S
b
,
X
3
)
,
t
∈
(
t
k
-
1
,
t
k
]
,
X
˙
2
=
f
i
k
(
X
1
+
S
b
,
X
3
)
-
a
X
2
,
X
˙
3
=
g
i
k
(
X
2
)
-
c
X
3
,
X
1
(
t
+
)
=
X
1
(
t
)
,
t
=
t
k
,
X
2
(
t
+
)
=
(
1
-
p
)
X
2
(
t
)
,
X
3
(
t
+
)
=
X
3
(
t
)
+
q
X
1
(
t
)
.
The vector system
(17)
X
˙
=
A
X
+
F
i
k
(
X
)
,
t
∈
(
t
k
-
1
,
t
k
]
,
X
(
t
+
)
=
B
X
,
t
=
t
k
,
where X
=
(
X
1
,
X
2
,
X
3
), A and B are 3
×
3 matrices, and F
i
k is a column vector given by
(18)
A
=
(
-
b
0
0
0
-
a
0
0
0
-
c
)
,
B
=
(
1
0
0
0
1
-
p
0
0
q
1
)
,
F
i
k
=
(
-
f
i
k
(
X
1
+
S
b
,
X
3
)
f
i
k
(
X
1
+
S
b
,
X
3
)
g
i
k
(
X
2
)
.
)
.
Assume that switching nonlinear functions f
i
(
X
) and g
i
(
X
) (for i
=
1,2
,
…
,
m) also satisfy (7). The following theorems give conditions for global asymptotic stability of the trivial solution of system (17) or the disease-free equilibrium of system (15).Theorem 5.
Assume that the switching ruleσ is periodic and that f
i
(
t
,
X
)
=
f
i
(
t
+
ω
,
X
) and g
i
(
t
,
X
)
=
g
i
(
t
+
ω
,
X
), for i
=
1,2
,
…
,
m, where f
i
(
X
) and g
i
(
X
) satisfy (7). If there exists a constant C
>
0, such that
(19)
m
ln
⁡
ρ
+
τ
1
λ
1
+
τ
2
λ
2
+
⋯
+
τ
m
λ
m
<
-
C
,
where ρ
=
max
⁡
{
1
,
(
(
1
-
p
)
2
+
q
2
+
1
+
(
(
1
-
p
)
2
+
q
2
+
1
)
2
-
4
(
1
-
p
)
2
)
/
2
} and λ
i
=
λ
+
2
L
i
2
+
R
i
2
+
1 with λ
=
max
⁡
{
-
2
a
,
-
2
b
,
-
2
c
}, then the trivial solution of system (17) is globally asymptotically stable, which implies that the disease-free equilibrium E
0 of system (15) is globally asymptotically stable.Proof.
Assume thati
=
i
k, for t
∈
(
t
k
-
1
,
t
k
], where i
k follows the switching rule σ
(
t
)
∈
I
Periodic. Define the common Lyapunov function U as follows: U
=
X
′
X. According to the proof of Theorem 1, we have U
˙
≤
λ
k
U. For t
∈
(
t
k
-
1
,
t
k
], it follows that
(20)
U
(
t
)
≤
U
(
t
k
-
1
+
)
exp
⁡
[
∫
t
k
-
1
t
λ
k
(
s
)
d
s
]
,
and immediately after the impulsive switch time t
=
t
k,
(21)
U
(
t
k
+
)
=
X
(
t
k
+
)
′
X
(
t
k
+
)
=
(
B
X
(
t
k
)
)
′
(
B
X
(
t
k
)
)
=
X
(
t
k
)
′
(
B
′
B
)
X
(
t
k
)
≤
ρ
X
(
t
k
)
′
X
(
t
k
)
=
ρ
U
(
t
k
)
.
Apply (20) and (21) on each subinterval. For t
∈
(
t
0
,
t
1
], U
(
t
)
≤
U
(
t
0
)
exp
⁡
[
∫
t
0
t
‍
λ
1
(
s
)
d
s
]. Further, U
(
t
1
+
)
≤
ρ
U
(
t
1
) and hence U
(
t
1
+
)
≤
ρ
U
(
t
0
)
exp
⁡
[
∫
t
0
t
1
‍
λ
1
(
s
)
d
s
]. For t
∈
(
t
1
,
t
2
], U
(
t
)
≤
U
(
t
1
+
)
exp
⁡
[
∫
t
1
t
‍
λ
2
(
s
)
d
s
]
≤
ρ
U
(
t
0
)
exp
⁡
[
∫
t
0
t
1
‍
λ
1
(
s
)
d
s
+
∫
t
1
t
‍
λ
2
(
s
)
d
s
]. Then U
(
t
2
+
)
≤
ρ
U
(
t
2
)
≤
ρ
2
U
(
t
0
)
exp
⁡
[
∫
t
0
t
1
‍
λ
1
(
s
)
d
s
+
∫
t
1
t
2
‍
λ
2
(
s
)
d
s
]. In general, for t
∈
(
t
k
-
1
,
t
k
], U
(
t
)
≤
U
(
t
k
-
1
+
)
exp
⁡
[
∫
t
k
-
1
t
‍
λ
k
(
s
)
d
s
]
≤
ρ
k
-
1
U
(
t
0
)
exp
⁡
[
∫
t
0
t
1
‍
λ
1
(
s
)
d
s
+
∫
t
1
t
2
‍
λ
2
(
s
)
d
s
+
⋯
+
∫
t
k
-
1
t
‍
λ
k
(
s
)
d
s
]. Then, for t
∈
(
ω
-
t
m
,
ω
], U
(
t
)
≤
ρ
m
-
1
U
(
t
0
)
exp
⁡
[
∫
t
0
t
1
‍
λ
i
1
(
s
)
d
s
+
⋯
+
∫
ω
-
t
m
t
‍
λ
i
m
(
s
)
d
s
] and immediately after the impulsive switch time t
=
ω,
(22)
U
(
ω
+
)
=
X
(
ω
+
)
′
X
(
ω
+
)
=
(
B
X
(
ω
)
)
′
B
X
(
ω
)
=
X
(
ω
)
′
B
′
B
X
(
ω
)
≤
ρ
X
(
ω
)
′
X
(
ω
)
=
ρ
U
(
ω
)
≤
ρ
m
U
(
t
0
)
exp
⁡
[
∫
t
0
t
1
λ
i
1
(
s
)
d
s
+
∫
t
1
t
2
λ
i
2
(
s
)
d
s
hhhhhhhhhhhhhh
+
⋯
+
∫
t
m
-
1
t
m
λ
i
m
(
s
)
d
s
]
=
U
(
t
0
)
exp
⁡
[
m
ln
⁡
ρ
+
∫
t
0
t
1
λ
i
1
(
s
)
d
s
+
∫
t
1
t
2
λ
i
2
(
s
)
d
s
hhhhhhhhhhh
+
⋯
+
∫
t
m
-
1
t
m
λ
i
m
(
s
)
d
s
]
=
U
(
t
0
)
exp
⁡
[
m
ln
⁡
ρ
+
τ
1
λ
1
+
τ
2
λ
2
+
⋯
+
τ
m
λ
m
]
=
ζ
U
(
t
0
)
,
where ζ
=
exp
⁡
[
m
ln
⁡
ρ
+
τ
1
λ
1
+
τ
2
λ
2
+
⋯
+
τ
m
λ
m
]
<
1 from (19). Similarly, it can be shown that U
(
h
ω
+
)
≤
ζ
U
(
(
h
-
1
)
ω
+
) for any integer h
=
1,2
,
… and so
(23)
U
(
h
ω
+
)
≤
ζ
U
(
(
h
-
1
)
ω
+
)
≤
ζ
(
ζ
U
(
(
h
-
2
)
ω
+
)
)
≤
⋯
≤
ζ
h
U
(
t
0
)
;
it follows that
(24)
lim
⁡
h
→
∞
U
(
h
ω
+
)
=
lim
⁡
h
→
∞
ζ
h
U
(
t
0
)
=
0
.
In general, for t
∈
(
t
k
-
1
,
t
k
] and h
ω
<
t
k
≤
(
h
+
1
)
ω, then U
(
t
)
≤
U
(
h
ω
+
)
ρ
m
-
1
exp
⁡
[
∫
t
0
t
1
‍
λ
i
1
(
s
)
d
s
+
∫
t
1
t
2
‍
λ
i
2
(
s
)
d
s
+
⋯
+
∫
t
k
-
1
t
‍
λ
i
k
(
s
)
d
s
]
≤
M
U
(
h
ω
+
)
ρ
m
-
1
≤
ρ
m
-
1
M
U
(
t
0
)
exp
⁡
[
h
ln
⁡
ζ
], where M
=
max
⁡
t
k
-
1
<
t
≤
t
k
⁡
exp
⁡
[
∫
t
0
t
1
‍
λ
i
1
(
s
)
d
s
+
∫
t
1
t
2
‍
λ
i
2
(
s
)
d
s
+
⋯
+
∫
t
k
-
1
t
‍
λ
i
k
(
s
)
d
s
]. Note that the ln
⁡
ζ
<
0; then, U
(
t
) converges to zero as t
→
∞ and h
→
∞. Thus, the trivial solution of system (17) is globally asymptotically stable, which implies that the disease-free equilibrium E
0 of system (15) is globally asymptotically stable.Remark 6.
In Theorem5, a general criterion, characterizing the switching term and the pulse term, is established to guarantee the global asymptotic stability of system (15).Remark 7.
As a special case, whenρ
=
1, the dynamics of the disease-free equilibrium are mainly determined by the switching term τ
k
λ
k, for k
=
1,2
,
….


## 4. Numerical Simulations

In order to illustrate the effectiveness of the proposed results above, the stability of HIV models with switching nonlinear incidence functions and pulse control is presented. Moreover, the comparison between results in HIV models with and without the switching effect is presented. Here, we assume thatt
0
=
0.Example 1.
Consider a HIV model with switching nonlinear incidence functions(25)
T
˙
=
S
-
b
T
-
10
(
1
-
γ
i
k
)
β
i
k
sin
(
T
V
)
V
,
t
∈
(
t
k
-
1
,
t
k
]
,
I
˙
=
10
(
1
-
γ
i
k
)
β
i
k
sin
(
T
V
)
V
-
a
I
,
V
˙
=
(
1
-
η
i
k
)
N
a
sin
(
I
)
I
-
c
V
with i
k
∈
{
1,2
} following a periodic switching rule
(26)
σ
(
t
)
=
{
1
,
if
t
∈
(
k
,
k
+
0.75
]
,
2
,
if
t
∈
(
k
+
0.75
,
k
+
1
]
,
k
=
0,1
,
2
,
…
,
and a HIV model without switching nonlinear incidence functions
(27)
T
˙
=
S
-
b
T
-
10
(
1
-
γ
)
β
sin
(
T
V
)
V
,
I
˙
=
10
(
1
-
γ
)
β
sin
(
T
V
)
V
-
a
I
,
V
˙
=
(
1
-
η
)
N
a
sin
(
I
)
I
-
c
V
.
By (26), we have τ
1
=
0.75 and τ
2
=
0.25. In this example, we first consider the global asymptotic stability of system (25) with initial values T
0
=
1000, I
0
=
500, and V
0
=
500. Take parameters S
=
1
0
3, b
=
0.7, a
=
0.7, c
=
0.7, and N
=
2 and the switching parameters γ
1
=
0.9, γ
2
=
0.7, η
1
=
0.8, η
2
=
0.6, β
1
=
0.32, and β
2
=
0.15. For system (25) with these parameters specifications, we have ∥
10
(
1
-
γ
i
k
)
β
i
k
sin
(
T
V
)
V
∥
≤
10
(
1
-
γ
i
k
)
β
i
k
∥
V
∥ and ∥
(
1
-
η
i
k
)
N
a
sin
(
I
)
I
∥
≤
(
1
-
η
i
k
)
N
a
∥
I
∥ in (7) and τ
1
λ
1
+
τ
2
λ
2
=
-
0.0079 with λ
1
=
-
0.1168 and λ
2
=
0.3186 in (8). Choosing C
=
0.0080, the disease-free equilibrium (
T
,
I
,
V
)
=
(
S
/
b
,
0,0
) is globally asymptotically stable by Theorem 1. The trajectory (solid lines) of system (25) is plotted in Figure 1, which indicates that this system is globally asymptotically stable. Furthermore, taking the parameters γ
=
γ
2, η
=
η
2, and β
=
β
2 of the system (27), and other parameters being the same as that of the system (25), we present the solutions (dashed lines) of the system (27) in Figure 1 for comparison purposes. In the first graph, we see that the total population of uninfected cells corresponding to the switching effect was much higher than the one without the switching effect, while, for infected cells and virus with the switching effect, the total populations were lower than the one without the switching effect (see the second graph and the third graph). Thus, computer simulations of these models agree well with mathematical theory.Numerical solutions of system (25) and system (27). The solutions with the switching effect (solid lines) were compared with the solutions without the switching effect (dashed lines).
(a)
(b)
(c)Example 2.
Consider a HIV model with switching nonlinear incidence functions and pulse control,(28)
T
˙
=
S
-
b
T
-
(
1
-
γ
i
k
)
β
i
k
sin
(
T
V
)
V
,
t
∈
(
t
k
-
1
,
t
k
]
,
I
˙
=
(
1
-
γ
i
k
)
β
i
k
sin
(
T
V
)
V
-
a
I
,
V
˙
=
(
1
-
η
i
k
)
N
a
sin
(
I
)
I
-
c
V
,
T
(
t
+
)
=
T
(
t
)
,
t
=
t
k
,
I
(
t
+
)
=
(
1
-
p
)
I
(
t
)
,
V
(
t
+
)
=
V
(
t
)
+
q
I
(
t
)
,
with i
k
∈
{
1,2
} following the periodic switching rule σ
(
t
): τ
1
=
0.5, τ
2
=
0.5, m
=
2. According to Example 1, take L
i
k
=
(
1
-
γ
i
k
)
β
i
k and R
i
k
=
(
1
-
η
i
k
)
N
a in (7). Take parameters S
=
1
0
4, b
=
0.6, a
=
0.6, c
=
2, and N
=
2 and the switching parameters γ
1
=
0.4, γ
2
=
0.3, η
1
=
0.7, η
2
=
0.6, β
1
=
1
0
-
4, and β
2
=
1
0
-
5. Take the pulse control parameters p
=
0.3 and q
=
0.1. We can derive λ
1
=
-
0.0704, λ
2
=
0.0304, and ρ
=
1.0193. Thus, (19) is satisfied with C
=
1
0
-
5. By Theorem 5, we can draw the conclusion that the disease-free equilibrium (
T
,
I
,
V
)
=
(
S
/
b
,
0,0
) is globally asymptotically stable. Figure 2 shows a simulation result of this system. Moreover, it can be seen from Figure 2 that pulse vaccination is successful to cure the HIV infection.Figure 2
Numerical solutions of system (28).


## 5. Conclusions

In this paper, new HIV models with switching nonlinear incidence functions and pulse control are investigated. It is reasonable from a physical perspective that nonlinear incidence functions are assumed to be switching nonlinear incidence to incorporate into HIV models, since nonlinear incidence functions are changing in time, which may change functional form in time, due to changes in host behavior. For the periodic switching rule, some new sufficient conditions are established to ensure the global asymptotic stability of the disease-free equilibrium by constructing common Lyapunov functions. The obtained results have more advantages than those in [2, 7] and are very useful for a large class of infection disease models. The results indicated that the HIV model with the switching effect plays an important role in understanding the dynamics of the disease. Furthermore, taking pulse vaccination into the above model, a new HIV model with switching nonlinear incidence functions and pulse control is developed. Some sufficient conditions characterizing the pulse term and the switching term are derived to determine whether the pulse vaccination succeeded in preventing disease. Numerical examples are carried out to verify the proposed results. One future direction is to study multicity HIV infections models with switching parameters.



---

*Source: 853960-2014-10-16.xml*