# Time Delayed Stage-Structured Predator-Prey Model with Birth Pulse and Pest Control Tactics

Authors: Mei Yan; Yongfeng Li; Zhongyi Xiang
Journal: Abstract and Applied Analysis
 (2014)

---

## Abstract

Normally, chemical pesticides kill not only pests but also their natural enemies. In order to better control the pests, two-time delayed stage-structured predator-prey models with birth pulse and pest control tactics are proposed and analyzed by using impulsive differential equations in present work. The stability threshold conditions for the mature prey-eradication periodic solutions of two models are derived, respectively. The effects of key parameters including killing efficiency rate, pulse period, the maximum birth effort per unit of time of natural enemy, and maturation time of prey on the threshold values are discussed in more detail. By comparing the two threshold values of mature prey-extinction, we provide the fact that the second control tactic is more effective than the first control method.


---

## Body

## 1. Introduction

The outbreak of pest often triggers serious ecological and economic problems. In recent years, the management of pest has increasingly become the focus of attention. How to effectively control pest is one primarily concern problem. In practice, lots of factors can affect the efficiency of pest control, for instance, the time of impulsive effects, the number of prey stocked or naturally released, and the proportion of killing or catching pests. Mathematical modelling is one of the main ways for estimating and predicting the range of ecological interactions between pest and predator. Lately, many papers have been devoted to propose and analyze the predator-prey systems [1–4].In one aspect, many species have the life history that goes through two stages, immature and mature, which has significant morphological and behavioral differences between them. Therefore, it is necessary to account for these differences, and the dynamics of mathematical models with stage-structured prey-predator model has been widely studied [5–11] in recent years. In other aspects, in the natural world, the immature predator always undergo a certain time τ (which is called maturation time delay) to be mature, so time delays play an important role in biological meanings in age-structured models. Time delay and stage-structured systems are introduced into predator-prey models, which greatly enriches biologic background. Many age-structured models with time delay were extensively studied [12–14]. A single species model with stage structure and time delay was invested by Aiello and Freedman [15].As far as the population dynamics is concerned, most models often considered that the population reproduces throughout all year. However, many species give birth seasonally or in regular pules. In this regard, the continuous reproduction of mature species should be removed from the model and termed this growth form as birth pulse. For instance, Tang and Chen study an age-structured model with density-dependent birth pulse in [16] and Xiang et al. [17] have considered a delayed Lotka-Volterra model with birth pulse and impulsive effect at different moment on the prey. Theories of impulsive differential equations have been studied lately [18–20]. Impulsive systems are found in many domains of applied sciences. The application of impulsive differential equations to population dynamics is an interesting topic since it is reasonable and corrects in modeling the evolution of populations. Virtually, the pulses occur at different moments more realistic, which is more meaningful than the pulses that occur at the same time. Recently, the qualitative analysis of impulsive differential equations at different moments has been widely investigated in many works [21–23].Although, many authors have devoted to study the effects of pesticide on pest and its natural enemies and lots of instructive control strategies also have been given. An optimal time of pesticide applications still seems to be a novel interesting area. Following the practical pest management, we firstly propose the predator-prey model with pulse at the same time. Further, we assumed that the birth pulse and pest control tactics occur at different time. Discussing and comparing the mature prey-extinction of the two models, we get some new effective pest management.


## 2. Predator-Prey Model with Impulsive Effects at the Same Time

The purpose of this paper is to address how the time of impulsive effects influences the pest control. On the basis of the above discussion and motivated by [8, 16], we first will extend the following impulsive model which introduces the impulsive control tactics at the same time:
(1)
x
˙
1
(
t
)
=
r
x
2
(
t
)
-
r
e
-
ω
τ
x
2
(
t
-
τ
)
-
ω
x
1
(
t
)
,
x
˙
2
(
t
)
=
r
e
-
ω
τ
x
2
(
t
-
τ
)
-
β
x
2
(
t
)
y
(
t
)
-
d
1
x
2
(
t
)
-
d
2
x
2
2
(
t
)
,
y
˙
(
t
)
=
k
β
x
2
(
t
)
y
(
t
)
-
d
3
y
(
t
)
,
t
≠
n
T
,
x
1
(
t
+
)
=
(
1
-
p
0
)
x
1
(
t
)
+
μ
,
x
2
(
t
+
)
=
x
2
(
t
)
,
y
(
t
+
)
=
(
1
-
p
)
y
(
t
)
+
b
y
(
t
)
a
+
y
(
t
)
,
t
=
n
T
,
(
ω
1
(
s
)
,
ω
2
(
s
)
,
ω
3
(
s
)
)
∈
C
+
=
C
(
[
-
τ
,
0
]
,
R
+
3
)
,
ω
i
(
0
)
>
0
,
i
=
1,2
,
3
,
where x
1
(
t
) and x
2
(
t
) are densities of the immature prey and mature prey, respectively. y
(
t
) denotes the density of predator. τ represents a constant time to maturity and ω
,
d
1
,
d
2
,
d
3
,
k
,
a
,
b
,
β
,
r are positive constants. r is the birth rate of immature prey, ω
(
ω
>
d
3
), d
1 are the mortality rates of the immature prey and mature prey, and d
3
(
d
3
>
d
1
) is the death coefficient of y
(
t
). d
2 is the intraspecific competition rate of the mature prey. β is the transmission coefficient. The term r
e
-
ω
τ
x
2
(
t
-
τ
) represents the immatures who were born at time t
-
τ and survive at time t (with the immature death rate ω) and therefore represents the transformation of immatures to matures (r
e
-
ω
τ
>
d
1). k is the rate of conversing prey into predator. p
0
,
p
(
0
≤
p
0
,
p
<
1
) represent the killing (or poisoning) rate of the prey and the predator, respectively. μ
≥
0 is the stocking amount of the immature prey at t
=
n
T. b
y
(
t
)
/
(
a
+
y
(
t
)
) denotes the birth effort of predator population at t
=
n
T, where b means the maximum birth effort per unit of time and a the predator size at which 50% saturation (
b
y
(
t
)
/
(
a
+
y
(
t
)
)
|
y
(
t
)
=
a
=
b
/
2
) measures how soon saturation occurs.Lemma 1 (see [24]).
Consider the following equation:(2)
x
˙
(
t
)
=
a
x
(
t
-
τ
)
-
b
x
(
t
)
-
c
x
2
(
t
)
,
where a
,
b
,
c
,
τ
>
0 for -
τ
≤
t
≤
0. We have(i)
ifa
<
b, then lim
⁡
t
→
∞
x
(
t
)
=
0;
(ii)
ifa
>
b, then lim
⁡
t
→
∞
x
(
t
)
=
(
a
-
b
/
c
).Lemma 2.
Consider the following system:(3)
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
u
˙
(
t
)
=
-
d
3
u
(
t
)
,
t
≠
n
T
,
u
(
t
+
)
=
(
1
-
p
)
u
(
t
)
+
b
u
(
t
)
a
+
u
(
t
)
,
t
=
n
T
,
(i)
ifp
>
1
+
b
/
a
-
e
d
3
T, the fixed point u
0
*
=
0 is globally asymptotically stable;
(ii)
ifp
<
1
+
b
/
a
-
e
d
3
T, the fixed point u
*
=
(
a
(
1
-
p
)
+
b
-
a
e
d
3
T
)
/
(
1
-
(
1
-
p
)
e
-
d
3
T
) is globally asymptotically stable.Proof.
Integrating the first equation of (3) on n
T
<
t
<
(
n
+
1
)
T, we have
(4)
u
(
t
)
=
u
(
n
T
+
)
e
-
d
3
(
t
-
n
T
)
,
n
T
≤
t
≤
(
n
+
1
)
T
.
With the successive pulse, we can obtain the following stroboscopic map of system (3):
(5)
u
(
(
n
+
1
)
T
+
)
=
(
1
-
p
)
u
(
n
T
+
)
e
-
d
3
T
+
b
u
(
n
T
+
)
e
-
d
3
T
a
+
u
(
n
T
+
)
e
-
d
3
T
.
Let F
(
u
)
=
(
1
-
p
)
u
e
-
d
3
T
+
(
b
u
e
-
d
3
T
)
/
(
a
+
u
e
-
d
3
T
). If p
>
1
+
(
b
/
a
)
-
e
d
3
T, then the fixed point u
0
*
=
0, since
(6)
|
d
F
(
u
)
d
u
|
u
=
0
=
(
1
-
p
)
e
-
d
3
T
+
b
e
-
d
3
T
a
<
1
;
and unstable if p
<
1
+
(
b
/
a
)
-
e
d
3
T, since
(7)
|
d
F
(
u
)
d
u
|
u
=
0
=
(
1
-
p
)
e
-
d
3
T
+
b
e
-
d
3
T
a
>
1
.
If p
<
1
+
(
b
/
a
)
-
e
d
3
T, (5) has a unique positive fixed point:
(8)
u
*
=
a
(
1
-
p
)
+
b
-
a
e
d
3
T
1
-
(
1
-
p
)
e
-
d
3
T
.
By the same way, we can obtain the fact that the fixed point u
*
=
(
a
(
1
-
p
)
+
b
-
a
e
d
3
T
)
/
(
1
-
(
1
-
p
)
e
-
d
3
T
) is globally asymptotically stable if p
<
1
+
(
b
/
a
)
-
e
d
3
T. This completes the proof.which implies that corresponding periodic solution of (3) on (
n
T
,
(
n
+
1
)
T
] is
(9)
u
~
(
t
)
=
a
(
1
-
p
)
+
b
-
a
e
d
3
T
1
-
(
1
-
p
)
e
-
d
3
T
e
-
d
3
(
t
-
n
T
)
.
Thus (1) has a mature prey-extinction periodic solution (
x
~
1
(
t
)
,
0
,
y
~
(
t
)
), where x
~
1
(
t
)
=
(
μ
e
-
ω
(
t
-
n
T
)
)
/
(
1
-
(
1
-
p
0
)
e
-
ω
T
)
,
y
~
(
t
)
=
(
a
(
1
-
p
)
+
b
-
a
e
d
3
T
)
/
(
1
-
(
1
-
p
)
e
-
d
3
T
)
e
-
d
3
(
t
-
n
T
).Theorem 3.
The mature prey-extinction periodic solution(
x
~
1
(
t
)
,
0
,
y
~
(
t
)
) is globally asymptotically stable if
(10)
R
1
<
1
holds, where R
1
=
(
r
e
-
ω
τ
-
d
1
)
[
1
-
(
1
-
p
)
e
-
d
3
T
]
/
β
[
a
(
1
-
p
)
e
-
d
3
T
+
b
e
-
d
3
T
-
a
].Proof.
Let(
x
1
(
t
)
,
x
2
(
t
)
,
y
(
t
)
) be any solution of system (1). Following from the third equation of system (1), we notice that y
˙
(
t
)
≥
-
d
3
y
(
t
). Consider the following impulsive differential system:
(11)
y
˙
1
(
t
)
=
-
d
3
y
1
(
t
)
,
t
≠
n
T
,
y
1
(
t
+
)
=
(
1
-
p
)
y
1
(
t
)
+
b
y
1
(
t
)
a
+
y
1
(
t
)
,
t
=
n
T
.
According to Lemma 2, we obtain the fact that y
~
1
(
t
)
=
(
(
a
(
1
-
p
)
+
b
-
a
e
d
3
T
)
/
(
1
-
(
1
-
p
)
e
-
d
3
T
)
)
e
-
d
3
(
t
-
n
T
), for n
T
<
t
≤
(
n
+
1
)
T. By using the comparison differential theorem of impulsive equation, we have y
(
t
)
≥
y
1
(
t
) and y
1
(
t
)
→
y
~
1
(
t
) as t
→
∞. Therefore, for an arbitrarily small positive constant ε and all t large enough:
(12)
y
(
t
)
≥
y
1
(
t
)
≥
y
~
1
(
t
)
-
ε
holds true.
For simplicity, we assume that (12) holds for all t
>
0. From (12), we have
(13)
d
x
2
(
t
)
d
t
≤
r
e
-
ω
τ
x
2
(
t
-
τ
)
-
(
β
(
y
~
1
(
t
)
-
ε
)
+
d
1
)
x
2
(
t
)
-
d
2
x
2
2
(
t
)
,
t
>
n
T
+
τ
.
Consider the following comparison differential equation:
(14)
d
x
(
t
)
d
t
=
r
e
-
ω
τ
x
(
t
-
τ
)
-
(
β
(
y
~
1
(
t
)
-
ε
)
+
d
1
)
x
(
t
)
-
d
2
x
2
(
t
)
,
ccccccccccccccccccccccccc
t
>
n
T
+
τ
.
From (10) and Lemma 1, we have lim
⁡
t
→
∞
x
(
t
)
=
0. By the comparison differential theorem, we get lim
⁡
t
→
∞
x
2
(
t
)
<
lim
⁡
t
→
∞
x
(
t
)
=
0. Noticing the positivity of x
2
(
t
), we know that lim
⁡
t
→
∞
x
2
(
t
)
=
0. Without loss of generality, we assume that for any ε
1 (sufficiently small)
(15)
x
2
(
t
)
≤
ε
1
.
Following from (1) and (15), we can obtain that
(16)
-
r
e
-
ω
τ
ε
1
-
ω
x
1
(
t
)
≤
x
˙
1
(
t
)
≤
r
ε
1
-
ω
x
1
(
t
)
,
which yields z
1
(
t
)
≤
x
1
(
t
)
≤
z
2
(
t
) and z
1
(
t
)
→
x
~
1
(
t
)
,
z
2
(
t
)
→
x
~
1
(
t
) as t
→
∞
→
∞, while z
1
(
t
)
,
z
2
(
t
) are the solutions of
(17)
z
˙
1
(
t
)
=
-
r
e
-
ω
τ
ε
1
-
ω
z
1
(
t
)
,
t
≠
n
T
,
Δ
z
1
(
t
)
=
-
p
0
z
1
(
t
)
+
μ
,
t
=
n
T
,
n
∈
Z
+
,
(18)
z
˙
2
(
t
)
=
r
ε
1
-
ω
z
2
(
t
)
,
t
≠
n
T
,
Δ
z
2
(
t
)
=
-
p
0
z
2
(
t
)
+
μ
,
t
=
n
T
,
n
∈
Z
+
,
respectively. Here, z
~
1
(
t
)
=
-
r
e
-
ω
τ
ε
1
/
ω
+
(
(
μ
+
(
(
r
e
-
ω
τ
ε
1
)
/
ω
)
p
0
)
/
(
1
-
(
1
-
p
0
)
e
-
ω
T
)
)
e
-
ω
(
t
-
n
T
)and  z
~
2
(
t
)
=
(
r
ε
1
/
ω
)
+
(
(
μ
+
(
r
ε
1
/
ω
)
p
0
)
/
(
1
-
(
1
-
p
0
)
e
-
ω
T
)
)
e
-
ω
(
t
-
n
T
), for t
∈
(
n
T
,
(
n
+
1
)
T
].
In view of the comparison differential theorem, for anyε
2, let ε
1
→
0; then z
~
1
(
t
)
-
ε
2
<
x
1
(
t
)
<
z
~
2
(
t
)
+
ε
2, for t large enough, which means x
1
(
t
)
→
x
~
1
(
t
) as t
→
∞.
Next, we will prove thaty
(
t
)
→
y
~
(
t
) as t
→
∞. From (1) and (15), we get
(19)
-
d
3
y
(
t
)
≤
y
˙
(
t
)
≤
-
(
d
3
-
k
β
ε
1
)
y
(
t
)
.
For the left hand inequality, it follows from impulsive differential equation (5) that y
(
t
)
≥
y
1
(
t
) and y
1
(
t
)
→
y
~
1
(
t
) as t
→
∞. For the right hand inequality, we consider the following impulsive differential equation:
(20)
z
˙
3
(
t
)
=
-
(
d
3
-
k
β
ε
1
)
z
3
(
t
)
,
t
≠
n
T
,
z
3
(
t
+
)
=
(
1
-
p
)
z
3
(
t
)
+
b
z
3
(
t
)
a
+
z
3
(
t
)
,
t
=
n
T
,
n
∈
Z
+
,
and z
~
3
(
t
)
=
(
(
a
(
1
-
p
)
+
b
-
a
e
d
3
T
)
/
(
1
-
(
1
-
p
)
e
-
d
3
T
)
)
e
-
(
d
3
-
λ
β
ε
1
)
(
t
-
n
T
).
Therefore, for anyε
3, y
~
1
(
t
)
-
ε
3
<
y
(
t
)
<
z
~
3
(
t
)
+
ε
3. Let ε
1
→
0; then we get
(21)
y
~
(
t
)
-
ε
3
<
y
(
t
)
<
y
~
(
t
)
+
ε
3
,
for t large enough, which means y
(
t
)
→
y
~
(
t
) as t
→
∞. The proof is completed.Corollary 4.
(i) Ifp
<
1
+
(
b
β
/
(
(
r
e
-
ω
τ
-
d
1
)
+
a
β
)
)
-
e
d
3
T holds, then the mature prey-extinction periodic solution (
x
~
1
(
t
)
,
0
,
y
~
(
t
)
) is globally asymptotically stable.
(ii) IfT
<
(
1
/
d
3
)
ln
⁡
(
1
+
(
b
β
/
(
(
r
e
-
ω
τ
-
d
1
)
+
a
β
)
)
-
p
) holds, then the mature prey-extinction periodic solution (
x
~
1
(
t
)
,
0
,
y
~
(
t
)
) is globally asymptotically stable.
In biological terms, sincep
<
1
+
(
b
β
/
(
(
r
e
-
ω
τ
-
d
1
)
+
a
β
)
)
-
e
d
3
T or   T
<
(
1
/
d
3
)
ln
⁡
(
1
+
(
b
β
/
(
(
r
e
-
ω
τ
-
d
1
)
+
a
β
)
)
-
p
), the mature prey will extinct and immature prey and predators will coexist.


## 3. Predator-Prey Model with Impulsive Effects at the Different Time

In particular, in order to avoid the adverse effects of pesticides on the newly released natural enemies, we consider the following method implemented in practice to avoid such antagonism. That is, we assume that the pulse occurs att
=
(
n
+
l
)
T
(
0
<
l
<
1
) and t
=
(
n
+
1
)
T. The detailed changes are shown as follows:
(22)
x
˙
1
(
t
)
=
r
x
2
(
t
)
-
r
e
-
ω
τ
x
2
(
t
-
τ
)
-
ω
x
1
(
t
)
,
x
˙
2
(
t
)
=
r
e
-
ω
τ
x
2
(
t
-
τ
)
-
β
x
2
(
t
)
y
(
t
)
-
d
1
x
2
(
t
)
-
d
2
x
2
2
(
t
)
,
y
˙
(
t
)
=
k
β
x
2
(
t
)
y
(
t
)
-
d
3
y
(
t
)
,
t
≠
(
n
+
l
)
T
,
t
≠
(
n
+
1
)
T
,
x
1
(
t
+
)
=
(
1
-
p
0
)
x
1
(
t
)
,
x
2
(
t
+
)
=
x
2
(
t
)
,
y
(
t
+
)
=
(
1
-
p
)
y
(
t
)
,
t
=
(
n
+
l
)
T
,
x
1
(
t
+
)
=
x
1
(
t
)
+
μ
,
x
2
(
t
+
)
=
x
2
(
t
)
,
y
(
t
+
)
=
y
(
t
)
+
b
y
(
t
)
a
+
y
(
t
)
,
t
=
(
n
+
1
)
T
,
n
∈
Z
+
.
(
ω
1
(
s
)
,
ω
2
(
s
)
,
ω
3
(
s
)
)
∈
C
+
=
C
(
[
-
τ
,
0
]
,
R
+
3
)
,
ω
i
(
0
)
>
0
,
i
=
1,2
,
3
.
If the mature prey x
2
(
t
) is absent, then system (22) reduces to
(23)
x
˙
1
(
t
)
=
-
ω
x
1
(
t
)
,
y
˙
(
t
)
=
-
d
3
y
(
t
)
,
t
≠
(
n
+
l
)
T
,
t
≠
(
n
+
1
)
T
,
x
1
(
t
+
)
=
(
1
-
p
0
)
x
1
(
t
)
,
y
(
t
+
)
=
(
1
-
p
)
y
(
t
)
,
t
=
(
n
+
l
)
T
,
x
1
(
t
+
)
=
x
1
(
t
)
+
μ
,
y
(
t
+
)
=
y
(
t
)
+
b
y
(
t
)
a
+
y
(
t
)
,
t
=
(
n
+
1
)
T
,
n
∈
Z
+
.
We can easily obtain the analytic solution of system (23) at the interval [
n
T
,
(
n
+
1
)
T
):
(24)
x
1
(
t
)
=
{
x
1
(
n
T
+
)
e
-
ω
(
t
-
n
T
)
,
t
∈
(
n
T
,
(
n
+
l
)
T
]
,
x
1
(
(
n
+
l
)
T
+
)
e
-
ω
(
t
-
(
n
+
l
)
T
)
,
t
∈
(
(
n
+
l
)
T
,
(
n
+
1
)
T
]
,
y
(
t
)
=
{
y
(
n
T
+
)
e
-
d
3
(
t
-
n
T
)
,
t
∈
(
n
T
,
(
n
+
l
)
T
]
,
y
(
(
n
+
l
)
T
+
)
e
-
d
3
(
t
-
(
n
+
l
)
T
)
,
t
∈
(
(
n
+
l
)
T
,
(
n
+
1
)
T
]
.Lemma 5.
(i) Ifp
>
1
-
a
e
d
3
T
/
(
a
+
b
), the fixed point y
0
*
=
0 is globally asymptotically stable.
(ii) Ifp
<
1
-
a
e
d
3
T
/
(
a
+
b
), the fixed point y
*
=
(
(
a
+
b
)
(
1
-
p
)
-
a
e
d
3
T
)
/
(
1
-
p
)
[
1
-
(
1
-
p
)
e
-
d
3
T
] is globally asymptotically stable.Proof.
Consider system (24), which yields the following stroboscopic map:
(25)
y
(
(
n
+
1
)
T
+
)
=
(
1
-
p
)
y
(
n
T
+
)
e
-
d
3
T
+
b
(
1
-
p
)
y
(
n
T
+
)
e
-
d
3
T
a
+
(
1
-
p
)
y
(
n
T
+
)
e
-
d
3
T
.
Let F
(
y
)
=
(
1
-
p
)
y
e
-
d
3
T
+
b
(
1
-
p
)
y
e
-
d
3
T
/
(
a
+
(
1
-
p
)
y
e
-
d
3
T
). If p
>
1
-
a
e
d
3
T
/
(
a
+
b
), there exists a unique trivial fixed point y
0
*
=
0. The trivial fixed point y
0
* is locally stable if p
>
1
-
a
e
d
3
T
/
(
a
+
b
), since
(26)
|
d
F
(
y
)
d
y
|
y
=
0
=
a
(
1
-
p
)
e
-
d
3
T
+
b
(
1
-
p
)
e
-
d
3
T
a
<
1
,
and unstable if p
<
1
-
a
e
d
3
T
/
(
a
+
b
), since
(27)
|
d
F
(
y
)
d
y
|
y
=
0
=
a
(
1
-
p
)
e
-
d
3
T
+
b
(
1
-
p
)
e
-
d
3
T
a
>
1
.
If p
<
1
-
a
e
d
3
T
/
(
a
+
b
), (25) has a positive fixed point denoted by y
*, where y
*
=
(
(
a
+
b
)
(
1
-
p
)
-
a
e
d
3
T
)
/
(
1
-
p
)
[
1
-
(
1
-
p
)
e
-
d
3
T
], when p
<
1
-
a
e
d
3
T
/
(
a
+
b
). The positive fixed point y
* is locally stable since d
F
(
y
)
/
(
d
y
)
|
y
=
y
*
=
|
(
1
-
p
)
e
-
d
3
T
+
a
[
1
-
(
1
-
p
)
e
-
d
3
T
]
2
/
b
(
1
-
p
)
e
-
d
3
T
|
<
1, if p
<
1
-
a
e
d
3
T
/
(
a
+
b
).
Further, we can show thaty
* is globally asymptotically stable, if the following statements are satisfied:(i)
ify
*
>
y
>
0, then y
*
>
F
(
y
)
>
y;
(ii)
ify
>
y
*, then y
>
F
(
y
)
>
y
*.
By calculation, we getd
F
/
d
y
=
(
1
-
p
)
e
-
d
3
T+a
b
(
1
-
p
)
e
-
d
3
T
/
[
a
+
(
1
-
p
)
e
-
d
3
T
y
]
2
>
0, which shows that y
*
>
F
(
y
), when y
*
>
y
>
0; we know (28)
F
(
y
)
-
y
=
(
1
-
p
)
y
e
-
d
3
T
+
b
(
1
-
p
)
e
-
d
3
T
a
+
(
1
-
p
)
y
e
-
d
3
T
-
y
>
y
(
(
1
-
p
)
e
-
d
3
T
+
b
(
1
-
p
)
e
-
d
3
T
a
+
(
1
-
p
)
y
*
e
-
d
3
T
-
1
)
=
0
.
This yields F
(
y
)
>
y. So we have y
*
>
F
(
y
)
>
y under the assumption y
*
>
y
>
0; otherwise y
>
F
(
y
)
>
y
* if y
>
y
*. Thus the statements (i) and (ii) are satisfied. This completes the proof.Now we can deduce that the positive equilibriumy
* of system (25) is globally asymptotically stable. So we have the corresponding positive periodic solution (
x
~
1
(
t
)
,
y
~
(
t
)
) in the following, which is globally stable:
(29)
x
~
1
(
t
)
=
{
μ
1
-
(
1
-
p
0
)
e
-
ω
T
e
-
ω
(
t
-
n
T
)
,
t
∈
(
n
T
,
(
n
+
l
)
T
]
,
μ
(
1
-
p
0
)
1
-
(
1
-
p
0
)
e
-
ω
T
e
-
ω
(
t
-
n
T
)
,
t
∈
(
(
n
+
l
)
T
,
(
n
+
1
)
T
]
,
y
~
(
t
)
=
{
(
a
+
b
)
(
1
-
p
)
-
a
e
d
3
T
(
1
-
p
)
[
1
-
(
1
-
p
)
e
-
d
3
T
]
e
-
d
3
(
t
-
n
T
)
,
t
∈
(
n
T
,
(
n
+
l
)
T
]
,
(
a
+
b
)
(
1
-
p
)
-
a
e
d
3
T
1
-
(
1
-
p
)
e
-
d
3
T
e
-
d
3
(
t
-
n
T
)
,
t
∈
(
(
n
+
l
)
T
,
(
n
+
1
)
T
]
.Theorem 6.
Assume that(30)
R
1
*
<
1
holds, where R
1
*
=
(
r
e
-
ω
τ
-
d
1
)
[
1
-
(
1
-
p
)
e
-
d
3
T
]
/
β
[
(
1
-
p
)
(
a
+
b
)
e
-
d
3
T
-
a
]. Then the mature prey-extinction periodic solution (
x
~
1
(
t
)
,
0
,
y
~
(
t
)
) is globally asymptotically stable.Proof.
Following from the second equation of system (22), we notice that y
˙
(
t
)
≥
-
d
3
y
(
t
); thus consider the following impulsive differential system:
(31)
z
˙
4
(
t
)
=
-
d
3
z
4
(
t
)
,
t
≠
(
n
+
l
)
T
,
t
≠
(
n
+
1
)
T
,
Δ
z
4
(
t
)
=
-
p
z
4
(
t
)
,
t
=
(
n
+
l
)
T
,
Δ
z
4
(
t
)
=
b
z
4
(
t
)
a
+
z
4
(
t
)
,
t
=
(
n
+
1
)
T
,
n
∈
Z
+
.
According to Lemma 5 and the comparison theorem on impulsive differential equations, we have y
(
t
)
≥
z
4
(
t
) and z
4
(
t
)
→
y
~
(
t
) as t
→
∞. Then there exists an integer k
2
>
k
1, n
>
k
2, such that
(32)
y
(
t
)
≥
z
4
(
t
)
≥
y
~
(
t
)
-
ε
holds for allt large enough.
This is(33)
y
(
t
)
≥
y
~
(
t
)
-
ε
≥
{
(
a
+
b
)
(
1
-
p
)
-
a
e
d
3
T
(
1
-
p
)
[
1
-
(
1
-
p
)
e
-
d
3
T
]
e
-
d
3
l
T
,
t
∈
(
n
T
,
(
n
+
l
)
T
]
,
(
a
+
b
)
(
1
-
p
)
-
a
e
d
3
T
1
-
(
1
-
p
)
e
-
d
3
T
e
-
d
3
T
,
t
∈
(
(
n
+
l
)
T
,
(
n
+
1
)
T
]
.
Since(
(
(
a
+
b
)
(
1
-
p
)
-
a
e
d
3
T
)
/
(
1
-
p
)
[
1
-
(
1
-
p
)
e
-
d
3
T
]
)
e
-
d
3
l
T > (
(
(
a
+
b
)
(
1
-
p
)
-
a
e
d
3
T
)
/
(
1
-
(
1
-
p
)
e
-
d
3
T
)
)
e
-
d
3
T, we have
(34)
y
(
t
)
≥
(
a
+
b
)
(
1
-
p
)
-
a
e
d
3
T
1
-
(
1
-
p
)
e
-
d
3
T
e
-
d
3
T
-
ε
≜
ψ
,
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
t
∈
(
n
T
,
(
n
+
1
)
T
]
,
n
>
k
2
.
For simplicity, we assume that (34) holds for all t
≥
0. From (22) we have
(35)
d
x
2
(
t
)
d
t
≤
r
e
-
ω
τ
x
2
(
t
-
τ
)
-
(
β
(
y
~
(
t
)
-
ε
)
+
d
1
)
x
2
(
t
)
-
d
2
x
2
2
(
t
)
,
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
t
>
n
T
+
τ
,
n
>
k
2
.
Consider the following comparison differential system:
(36)
d
z
5
(
t
)
d
t
=
r
e
-
ω
τ
z
5
(
t
-
τ
)
-
(
β
(
y
~
(
t
)
-
ε
)
+
d
1
)
z
5
(
t
)
-
d
2
z
5
2
(
t
)
,
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
t
>
n
T
+
τ
,
n
>
k
2
.
From (23) and Lemma 1, we have lim
⁡
t
→
∞
z
5
(
t
)
=
0. By the comparison differential theorem, we have lim
⁡
t
→
∞
x
2
(
t
)
<
lim
⁡
t
→
∞
z
5
(
t
)
=
0. Noticing the positivity of x
2
(
t
), we know that lim
⁡
t
→
∞
x
2
(
t
)
=
0. Without loss of generality, we assume that there exists an integer k
3
,
k
3
>
k
2, for any ε
1 (sufficiently small), such that
(37)
x
2
(
t
)
≤
ε
1
.
In view of (22) and (37), it is easy to obtain that
(38)
-
r
e
-
ω
τ
ε
1
-
ω
x
1
(
t
)
≤
x
˙
1
(
t
)
≤
r
ε
1
-
ω
x
1
(
t
)
,
we get z
6
(
t
)
≤
x
1
(
t
)
≤
z
7
(
t
), and z
6
(
t
)
→
x
~
1
(
t
)
,
z
7
(
t
)
→
x
~
1
(
t
) as t
→
∞, where
(39)
z
˙
6
(
t
)
=
-
r
e
-
ω
τ
ε
1
-
ω
z
6
(
t
)
,
t
≠
(
n
+
l
)
T
,
t
≠
(
n
+
1
)
T
,
z
6
(
t
+
)
=
(
1
-
p
0
)
z
6
(
t
)
,
t
=
(
n
+
l
)
T
,
z
6
(
t
+
)
=
z
6
(
t
)
+
b
z
6
(
t
)
a
+
z
6
(
t
)
,
t
=
(
n
+
1
)
T
,
n
∈
Z
+
,
z
˙
7
(
t
)
=
r
ε
1
-
ω
z
7
(
t
)
,
t
≠
(
n
+
l
)
T
,
t
≠
(
n
+
1
)
T
,
z
7
(
t
+
)
=
(
1
-
p
0
)
z
7
(
t
)
,
t
=
(
n
+
l
)
T
,
z
7
(
t
+
)
=
z
7
(
t
)
+
b
z
7
(
t
)
a
+
z
7
(
t
)
,
t
=
(
n
+
1
)
T
,
n
∈
Z
+
.
From the comparison differential theorem, for any ε
2, there exists an integer k
4
,
n
>
k
4; let ε
1
→
0; then z
~
6
(
t
)
-
ε
2
<
x
1
(
t
)
<
z
~
7
(
t
)
+
ε
2, for t large enough, which means x
1
(
t
)
→
x
~
1
(
t
) as t
→
∞, where
(40)
z
~
6
(
t
)
=
{
-
r
e
-
ω
τ
ε
1
ω
+
(
μ
1
-
(
1
-
p
0
)
e
-
ω
T
-
-
r
e
-
ω
τ
ε
1
ω
)
e
-
ω
(
t
-
n
T
)
,
t
∈
(
n
T
,
(
n
+
l
)
T
]
,
-
r
e
-
ω
τ
ε
1
ω
+
(
μ
(
1
-
p
0
)
e
-
ω
l
T
1
-
(
1
-
p
0
)
e
-
ω
T
-
-
r
e
-
ω
τ
ε
1
ω
)
e
-
ω
(
t
-
(
n
+
l
)
T
)
,
t
∈
(
(
n
+
l
)
T
,
(
n
+
1
)
T
]
,
z
~
7
(
t
)
=
{
r
ε
1
ω
+
(
μ
1
-
(
1
-
p
0
)
e
-
ω
T
-
r
ε
1
ω
)
e
-
ω
(
t
-
n
T
)
,
t
∈
(
n
T
,
(
n
+
l
)
T
]
,
r
ε
1
ω
+
(
μ
(
1
-
p
0
)
e
-
ω
l
T
1
-
(
1
-
p
0
)
e
-
ω
T
-
r
ε
1
ω
)
e
-
ω
(
t
-
(
n
+
l
)
T
)
,
t
∈
(
(
n
+
l
)
T
,
(
n
+
1
)
T
]
.
Next, we will prove that y
(
t
)
→
y
~
(
t
) as t
→
∞. From (22) and (37), we get
(41)
-
d
3
y
(
t
)
≤
y
˙
(
t
)
≤
-
(
d
3
-
λ
β
ε
1
)
y
(
t
)
.
For the left hand inequality, it follows from impulsive differential equation (31) that y
(
t
)
≥
z
4
(
t
) and z
4
(
t
)
→
y
~
1
(
t
) as t
→
∞. For the right hand inequality, we consider the following impulsive differential equation:
(42)
z
˙
8
(
t
)
=
-
(
d
3
-
λ
β
ε
1
)
z
8
(
t
)
,
t
≠
(
n
+
l
)
T
,
t
≠
(
n
+
1
)
T
,
Δ
z
8
(
t
)
=
-
p
z
8
(
t
)
,
t
=
(
n
+
l
)
T
,
Δ
z
8
(
t
)
=
b
z
8
(
t
)
a
+
z
8
(
t
)
,
t
=
(
n
+
1
)
T
,
n
∈
Z
+
,
z
~
8
(
t
)
=
{
(
a
+
b
)
(
1
-
p
)
-
a
e
(
d
3
-
λ
β
ε
1
)
T
(
1
-
p
)
[
1
-
(
1
-
p
)
e
-
(
d
3
-
λ
β
ε
1
)
T
]
e
-
(
d
3
-
λ
β
ε
1
)
(
t
-
n
T
)
,
t
∈
(
n
T
,
(
n
+
l
)
T
]
,
(
a
+
b
)
(
1
-
p
)
-
a
e
(
d
3
-
λ
β
ε
1
)
T
1
-
(
1
-
p
)
e
-
(
d
3
-
λ
β
ε
1
)
T
e
-
(
d
3
-
λ
β
ε
1
)
(
t
-
n
T
)
,
t
∈
(
(
n
+
l
)
T
,
(
n
+
1
)
T
]
.
Therefore, for any ε
3, there exists an integer k
5
,
n
>
k
5 such that
(43)
z
~
4
(
t
)
-
ε
3
<
y
(
t
)
<
z
~
8
(
t
)
+
ε
3
.
Let ε
1
→
0; then we get
(44)
y
~
(
t
)
-
ε
3
<
y
(
t
)
<
y
~
(
t
)
+
ε
3
,
for t large enough, which means y
(
t
)
→
y
~
(
t
) as t
→
∞. The proof is completed.Corollary 7.
(i) Ifp
<
1
-
[
(
r
e
-
ω
τ
-
d
1
)
+
a
β
]
e
d
3
T
/
(
(
r
e
-
ω
τ
-
d
1
)
+
(
a
+
b
)
β
) holds, then the mature prey-extinction periodic solution (
x
~
1
(
t
)
,
0
,
y
~
(
t
)
) is globally asymptotically stable.
(ii) IfT
<
(
1
/
d
3
)
ln
⁡
⁡
(
[
(
r
e
-
ω
τ
-
d
1
)
+
(
a
+
b
)
β
]
(
1
-
p
)
/
[
(
r
e
-
ω
τ
-
d
1
)
+
a
β
]
e
d
3
T
) holds, then the mature prey-extinction periodic solution (
x
~
1
(
t
)
,
0
,
y
~
(
t
)
) is globally asymptotically stable.
The biological significance of Corollary7 is the same as Corollary 4, so we omit it.
In order to investigate the permanence of system (22), we should give the following Definition and Lemma.Lemma 8.
There exists a constantM
>
0 such that x
1
(
t
)
≤
M
/
k
,
x
2
(
t
)
≤
M
/
k
,
y
(
t
)
≤
M for each solution (
x
1
(
t
)
,
x
2
(
t
)
,
y
(
t
)
) of system (22) with all t large enough.Proof.
DefineV
(
t
)
=
k
x
1
(
t
)
+
k
x
2
(
t
)
+
y
(
t
), note that ω
>
d
3, then t
≠
(
n
+
l
)
T
,
t
≠
(
n
+
1
)
T and we have
(45)
D
+
V
(
t
)
+
d
3
V
(
t
)
=
-
k
(
ω
-
d
3
)
x
1
(
t
)
+
k
(
r
+
d
3
-
d
1
)
x
2
(
t
)
-
k
d
2
x
2
2
(
t
)
≤
-
k
d
2
(
x
2
2
(
t
)
-
r
+
d
3
-
d
1
2
d
2
)
2
+
k
(
r
+
d
3
-
d
1
)
2
4
d
2
≤
k
(
r
+
d
3
-
d
1
)
2
4
d
2
≜
M
1
.
When t
=
(
n
+
l
)
T,
(46)
V
(
(
n
+
l
)
T
+
)
=
k
x
1
(
(
n
+
l
)
T
+
)
+
k
x
2
(
(
n
+
l
)
T
+
)
+
y
(
(
n
+
l
)
T
+
)
=
k
(
(
1
-
p
0
)
)
x
1
(
(
n
+
l
)
T
)
+
k
x
2
(
(
n
+
l
)
T
)
+
(
1
-
p
)
y
(
(
n
+
l
)
T
)
≤
V
(
(
n
+
l
)
T
)
.
When t
=
(
n
+
1
)
T,
(47)
V
(
(
n
+
1
)
T
+
)
=
k
x
1
(
(
n
+
1
)
T
+
)
+
k
x
2
(
(
n
+
1
)
T
+
)
+
y
(
(
n
+
1
)
T
+
)
=
V
(
(
n
+
1
)
T
)
+
k
μ
+
b
y
(
(
n
+
1
)
T
)
a
+
y
(
(
n
+
1
)
T
)
≤
V
(
(
n
+
1
)
T
)
+
k
μ
+
b
.
Let us denote M
2
=
k
μ
+
b.
It follows from the comparison theorem of impulsive differential equations (see lemma 2.2 [25], page 23) that
(48)
V
(
t
)
≤
V
(
0
+
)
e
-
d
3
t
+
M
1
1
-
e
-
d
3
t
d
3
+
M
2
e
-
d
3
(
t
-
T
)
1
-
e
d
3
T
+
M
2
e
d
3
T
e
d
3
T
-
1
⟶
M
1
d
3
+
M
2
e
d
3
T
e
d
3
T
-
1
,
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
as
t
⟶
∞
.
So V
(
t
) is ultimately bounded. Hence, by the definition of V
(
t
), there exists a constant M
>
0, such that x
1
(
t
)
≤
M
/
k
,
x
2
(
t
)
≤
M
/
k
,
y
(
t
)
≤
M for all t large enough. This completes the proof.Remark 9.
According to Lemma8, it is clear that limsup
⁡
t
→
∞
V
(
t
)
→
M
1
/
d
3
+
M
2
(
e
d
3
T
/
(
e
d
3
T
-
1
)
). For convenience, we note that M
=
(
k
(
r
+
d
3
-
d
1
)
2
/
4
d
2
d
3
)
+
(
k
μ
+
b
)
(
e
d
3
T
/
(
e
d
3
T
-
1
)
).Definition 10.
The system (22) is said to be permanent if there are constants m
,
M
>
0 (independent of all initial values) and a finite time T
0, such that for all solutions (
x
1
(
t
)
,
x
2
(
t
)
,
y
(
t
)
) with all initial values x
1
(
0
+
)
>
0
,
x
2
(
0
+
)
>
0
,
y
(
0
+
)
>
0
,
m
≤
x
1
(
t
)
≤
M
/
k
,
m
≤
x
2
(
t
)
≤
M
/
k, and m
≤
y
(
t
)
≤
M holds for all t
≥
T
0.In biological terms, the permanence of (22) implies that prey (both immature and mature) and predators will coexist, none of them facing extinction or growing indefinitely. Denote
(49)
R
*
=
(
a
+
b
)
(
1
-
p
)
-
a
e
(
d
3
-
k
β
θ
)
T
(
1
-
p
)
[
1
-
(
1
-
p
)
e
-
(
d
3
-
k
β
θ
)
T
]
,
θ
=
(
a
+
b
)
(
1
-
p
)
a
e
d
3
T
(
1
-
1
R
2
*
)
+
1
R
2
*
.Theorem 11.
If(50)
R
2
*
>
1
holds, where R
2
*
=
(
(
r
e
-
ω
τ
-
d
1
-
d
2
(
M
/
k
)
)
/
β
R
*
), system (22) is permanent.Proof.
It is seen that the second equation of system (22) can be rewritten as
(51)
d
x
2
(
t
)
d
t
=
(
r
e
-
ω
τ
-
β
y
(
t
)
-
d
1
-
d
2
x
2
(
t
)
)
x
2
(
t
)
-
r
e
-
ω
τ
d
d
t
∫
t
-
τ
t
‍
x
2
(
u
)
d
u
.
Let us consider any positive solution (
x
1
(
t
)
,
x
2
(
t
)
,
y
(
t
)
) of system (22). By (51), V
(
t
) is defined as
(52)
V
(
t
)
=
x
2
(
t
)
+
r
e
-
ω
τ
d
d
t
∫
t
-
τ
t
‍
x
2
(
u
)
d
u
.
Taking the first derivative of V
(
t
) with respect to t, we have
(53)
d
V
(
t
)
d
t
=
(
r
e
-
ω
τ
-
β
y
(
t
)
-
d
1
-
d
2
x
2
(
t
)
)
x
2
(
t
)
.
As seen in Lemma 8, (53) can be written as
(54)
d
V
(
t
)
d
t
>
(
r
e
-
ω
τ
-
β
y
(
t
)
-
d
1
-
d
2
M
k
)
x
2
(
t
)
.
From the definition of R
2
*, we have
(55)
e
k
β
T
θ
=
(
a
+
b
)
(
1
-
p
)
a
e
d
3
T
(
1
-
1
R
2
*
)
+
1
R
2
*
>
1
,
and then θ
>
0; it is easy to know that there exists a sufficiently small ε such that
(56)
r
e
-
ω
τ
>
β
(
R
*
+
ε
)
+
d
1
+
d
2
M
k
.
We claim that, for any t
0
>
0, it is impossible that x
2
(
t
)
<
θ for all t
>
t
0. Otherwise, there is a t
0
>
0 such that θ for all t
>
t
0. It follows from the third equation of (22) that
(57)
y
˙
(
t
)
<
(
k
β
θ
-
d
3
)
y
(
t
)
,
for all t
>
t
0. Consider the following comparison impulsive system of (57):
(3)
z
˙
9
(
t
)
=
-
(
d
3
-
k
β
θ
)
z
9
(
t
)
,
t
≠
(
n
+
l
)
T
,
t
≠
(
n
+
1
)
T
,
Δ
z
9
(
t
)
=
-
p
z
9
(
t
)
,
t
=
(
n
+
l
)
T
,
Δ
z
9
(
t
)
=
b
z
9
(
t
)
a
+
z
9
(
t
)
,
t
=
(
n
+
1
)
T
,
n
∈
Z
+
,
and we have
(59)
z
~
9
(
t
)
=
{
(
a
+
b
)
(
1
-
p
)
-
a
e
(
d
3
-
k
β
θ
)
T
(
1
-
p
)
[
1
-
(
1
-
p
)
e
-
(
d
3
-
k
β
θ
)
T
]
e
-
(
d
3
-
k
β
θ
)
(
t
-
n
T
)
,
t
∈
(
n
T
,
(
n
+
l
)
T
]
,
(
a
+
b
)
(
1
-
p
)
-
a
e
(
d
3
-
k
β
θ
)
T
[
1
-
(
1
-
p
)
e
-
(
d
3
-
k
β
θ
)
]
e
-
(
d
3
-
k
β
θ
)
(
t
-
n
T
)
,
t
∈
(
(
n
+
l
)
T
,
(
n
+
1
)
T
]
.
In view of the comparison differential theorem for equation, there exists a T
1
>
0, such that
(60)
y
(
t
)
<
z
9
(
t
)
<
z
~
9
(
t
)
+
ε
1
≤
{
(
a
+
b
)
(
1
-
p
)
-
a
e
(
d
3
-
k
β
θ
)
T
(
1
-
p
)
[
1
-
(
1
-
p
)
e
-
(
d
3
-
k
β
θ
)
T
]
+
ε
1
,
t
∈
(
n
T
,
(
n
+
l
)
T
]
,
e
-
(
d
3
-
k
β
θ
)
l
T
(
a
+
b
)
(
1
-
p
)
-
a
e
(
d
3
-
k
β
θ
)
T
[
1
-
(
1
-
p
)
e
-
(
d
3
-
k
β
θ
)
]
+
ε
1
,
t
∈
(
(
n
+
l
)
T
,
(
n
+
1
)
T
]
,
for t
>
T
1, which implies that y
(
t
)
≤
z
~
9
(
t
)
+
ε
1, where
(61)
z
~
9
(
t
)
≤
(
a
+
b
)
(
1
-
p
)
-
a
e
(
d
3
-
k
β
θ
)
T
(
1
-
p
)
[
1
-
(
1
-
p
)
e
-
(
d
3
-
k
β
θ
)
T
]
=
R
*
.
Thus
(62)
y
(
t
)
≤
R
*
+
ε
1
≜
η
,
t
>
T
1
.
From (54), we get
(63)
r
e
-
ω
τ
>
β
η
+
d
1
+
d
2
M
k
.
By (54) and (63), we have
(64)
V
˙
(
t
)
>
(
r
e
-
ω
τ
-
β
η
-
d
1
-
d
2
M
k
)
x
2
(
t
)
,
t
≥
T
1
.
Let
(65)
x
2
m
=
min
⁡
t
∈
[
T
1
,
T
1
+
τ
]
x
2
(
t
)
.
In the following we show thatx
2
(
t
)
≥
x
2
m for all t
≥
T
1. Otherwise, there is a T
2
>
0 such that x
2
(
t
)
≥
x
2
m for T
1
≤
t
≤
T
1
+
τ
+
T
2
,
x
2
(
T
1
+
τ
+
T
2
)
=
x
2
m and x
˙
2
(
T
1
+
τ
+
T
2
)
<
0, which implies that,
(66)
x
˙
2
(
T
1
+
T
2
+
τ
)
=
r
e
-
ω
τ
x
2
(
T
1
+
T
2
+
τ
)
-
β
x
2
(
T
1
+
T
2
+
τ
)
y
(
T
1
+
T
2
+
τ
)
-
d
1
x
2
(
T
1
+
T
2
+
τ
)
-
d
2
x
2
2
(
T
1
+
T
2
+
τ
)
≥
(
r
e
-
ω
τ
-
β
η
-
d
1
-
d
2
M
)
x
2
m
>
0
,
it is a contradiction. Thus, x
2
≥
x
2
m for all t
≥
T
1. Hence, for t
≥
T
1, we get
(67)
V
˙
(
t
)
>
x
2
m
(
r
1
e
-
ω
τ
-
β
η
-
d
1
-
d
2
M
)
>
0
,
which means V
(
t
)
→
∞ as t
→
∞. This is contrary to the boundedness of V
(
t
). The claim is proved.
By the claim, consider the following two possibilities.
Case  1.
x
2
(
t
)
≥
x
2
* for all t large enough.
Case  2.
x
2 oscillates about x
2
* for t large enough. Set
(68)
q
=
min
⁡
{
x
2
*
2
,
x
2
*
e
-
(
β
M
+
d
1
+
d
2
x
2
*
)
τ
}
.
We will show thatx
2
(
t
)
≥
q for all t large enough. The conclusion is evident in Case  1. For Case  2, suppose that there exists the positive constant t
¯
,
ξ satisfying x
2
(
t
¯
)
=
x
2
(
t
¯
+
ξ
)
=
x
2
* and x
2
(
t
)
<
x
2
* for all t
¯
<
t
<
t
¯
+
ξ, where t
¯ is enough large such that
(69)
x
2
(
t
)
>
q
for
t
¯
<
t
<
t
¯
+
ξ
,
since x
2
(
t
) is continuous and bounded and is not affected by impulses. Thus x
2
(
t
) is uniformly continuous. Consequently, there exists a constant T
3
(
0
<
T
3
<
τ and T
3 is independent of the choice of  t
¯) such that x
2
(
t
)
>
(
x
2
*
/
2
) for t
¯
≤
t
≤
t
¯
+
T
3. If ξ
≤
T
3, our aim is obtained. Then we consider the case T
3
<
ξ
<
τ. From the second equation of (22) and the above assumption, we have that x
˙
2
(
t
)
≥
-
(
β
M
+
d
1
+
d
2
x
2
*
)
x
2
(
t
) for t
¯
<
t
≤
t
¯
+
ξ
≤
t
¯
+
τ. Hence one obtains x
2
(
t
)
≥
x
2
*
e
-
(
β
M
+
d
1
+
d
2
x
2
*
)
τ. It is clear that x
2
(
t
)
≥
q for t
¯
<
t
≤
t
¯
+
ξ. If ξ
≥
τ, then we have x
2
(
t
)
≥
q for t
¯
<
t
≤
t
¯
+
τ. The same arguments can be continued; we can obtain x
2
(
t
)
≥
q for t
¯
+
τ
<
t
≤
t
¯
+
ξ. Since the kind of interval [
t
¯
,
t
¯
+
ξ
] is chosen in arbitrary way, we have x
2
(
t
)
≥
q for t large enough. In view of the above discussion, the selection of q is independent of the positive solution of (22) which satisfies x
2
(
t
)
≥
q for all t large enough.
Next, from the first, the fourth, and the seventh equation of system (22), we have
(70)
x
˙
1
(
t
)
≥
-
ω
x
1
(
t
)
,
t
≠
(
n
+
l
)
T
,
t
≠
(
n
+
1
)
T
,
Δ
x
1
(
t
)
=
-
p
0
x
1
(
t
)
,
t
=
(
n
+
l
)
T
,
Δ
x
1
(
t
)
=
μ
,
t
=
(
n
+
1
)
T
,
n
∈
Z
+
.
Considering the comparison system,
(71)
z
˙
10
(
t
)
=
-
ω
z
10
(
t
)
,
t
≠
(
n
+
l
)
T
,
t
≠
(
n
+
1
)
T
,
Δ
z
10
(
t
)
=
-
p
0
z
10
(
t
)
,
t
=
(
n
+
l
)
T
,
Δ
z
10
(
t
)
=
μ
,
t
=
(
n
+
1
)
T
,
n
∈
Z
+
.
Therefore for any ε
4
>
0 small enough such that x
1
(
t
)
≥
z
~
10
(
t
)
-
ε
4 for t sufficiently large,
(72)
z
~
10
(
t
)
=
{
μ
1
-
(
1
-
p
0
)
e
-
ω
T
e
-
ω
(
t
-
n
T
)
,
t
∈
(
n
T
,
(
n
+
l
)
T
]
,
μ
(
1
-
p
0
)
1
-
(
1
-
p
0
)
e
-
ω
T
e
-
ω
(
t
-
n
T
)
,
t
∈
(
(
n
+
l
)
T
,
(
n
+
1
)
T
]
.
In view of the comparison theorem of impulsive differential equation and the Theorem 6, we get
(73)
x
1
(
t
)
≥
μ
(
1
-
p
0
)
1
-
(
1
-
p
0
)
e
-
ω
T
e
-
ω
T
-
ε
2
≜
m
1
.
From (34) and (68), set ψ
≜
m
3 and q
≜
m
2.
Definem
=
min
⁡
{
m
1
,
m
2
,
m
3
}, and thus we have x
1
(
t
)
≥
m
,
x
2
(
t
)
≥
m
,
y
(
t
)
≥
m. By Lemma 8 and the above discussion, the system is permanent. The proof is completed.


## 4.  Numerical Analysis

In the previous sections, we introduced the analytical tools and used them for a qualitative analysis of the system obtaining some results about the dynamics of the system. In this section, we perform a numerical analysis of the model based on the previous results. What we are interested in is how the key factors affect the thresholdsR
1 and R
1
*. Let R
1
(
p
,
T
,
τ
) = (
(
r
e
-
ω
τ
-
d
1
)
[
1
-
(
1
-
p
)
e
-
d
3
T
]
/
β
[
a
(
1
-
p
)
e
-
d
3
T
+
b
e
-
d
3
T
-
a
]
). Taking the first derivatives R
1
(
p
,
T
,
τ
) with respect to p
,
T and τ, respectively, one obtains
(74)
∂
R
1
∂
p
=
(
r
e
-
ω
τ
-
d
1
)
β
b
e
-
2
d
3
T
[
a
(
1
-
p
)
e
-
d
3
T
+
b
e
-
d
3
T
-
a
]
2
>
0
,
∂
R
1
∂
T
=
(
r
e
-
ω
τ
-
d
1
)
β
d
3
b
e
-
d
3
T
[
a
(
1
-
p
)
e
-
d
3
T
+
b
e
-
d
3
T
-
a
]
2
>
0
,
∂
R
1
∂
τ
=
-
ω
r
e
-
ω
τ
[
1
-
(
1
-
p
)
e
-
d
3
T
]
β
[
a
(
1
-
p
)
e
-
d
3
T
+
b
e
-
d
3
T
-
a
]
<
0
.
Those inequalities demonstrate that threshold value R
1 is a monotonic increasing function with respect to p and T and is a monotonic decreasing function with τ. These results indicate that the smaller killing (or poisoning) rate p or the shorter pulse period T, the smaller the threshold value R
1 which follows, so, the result shows that it is more effective for pest control. Similarly, when pests have very short maturity τ, the bigger the threshold value R
1 which follows, and hence make it more difficult for pest control. Therefore, we must carefully select the impulsive catching (or poisoning) rate or the time of pesticide applications or reduce the dosages of the pesticide.So far, we have considered the global attractive mature prey-eradication solution of systems (1) and (22). We first focus on the system (1); when threshold value R
1
<
1, the mature prey is eradicated totally and the immature prey and the predator population will tend to a stable level; see Figure 1, where r
=
1.2
,
τ
=
0.5, ω
=
0.2
,
β
=
2
,
μ
=
2
,
k
=
0.2
,
d
1
=
0.1
,
d
2
=
0.2, d
3
=
0.15
,
p
0
=
0.1
,
p
=
0.5, a
=
5
,
b
=
6
,
T
=
3, and R
1
≈
0.7997
<
1, while if we chose the parameters set as those in Figure 2, then the all populations oscillate periodically, which indicates that the system (1) is permanent. Furthermore, we can carry out two-parameter bifurcation analysis for the threshold value R
1, as those shown in Figures 3(a)–3(d). In each subplot, we let two key parameters vary simultaneously and see how those parameters affect the threshold value R
1. All simulation results shown in Figure 3 indicate that the R
1 appears to be quite sensitive to small changes in parameters, p
,
T
,
τ, and b. However, we note that the dynamical behaviour of system (1) is dominated by several types of periodic solutions and their coexistence, even with very high intrinsic growth rates. In Figure 4, we studied the influence of killing rate p on the complexity of system (1); numerical results show that when killing rate p varies from 0.45 to 0.95, the dynamical behaviour of system (1) is very complex. A typical chaos oscillation is captured when p
=
0.7 (see Figure 5). Figure 6 shows the bifurcation diagrams of system (1); the system analyzed here could take on many forms of complexity, including period doubling bifurcation, symmetry-breaking pitchfork bifurcation, and chaotic solutions. Meanwhile, system (1) also has a certain type of coexistence if the intrinsic growth rate of the pest population is large enough.Dynamical behavior of system (1); the parameter values are as follows: r
=
1.2
,
τ
=
0.5
,
ω
=
0.2, β
=
2
,
μ
=
2
,
k
=
0.2, d
1
=
0.1
,
d
2
=
0.2
,
d
3
=
0.15, p
0
=
0.1
,
p
=
0.5, a
=
5
,
b
=
6, and T
=
3.
(a)
(b)
(c)Dynamical behavior of system (1); the parameter values are as follows: r
=
1.2
,
τ
=
0.5
,
ω
=
0.2
,
β
=
1.2, k
=
0.2
,
d
1
=
0.1
,
d
2
=
0.2, d
3
=
0.15
,
p
0
=
0.1
,
p
=
0.7, a
=
5
,
b
=
6
,
μ
=
2, and T
=
3.
(a)
(b)
(c)The influence of some key parameters on the threshold levelR
1; the other parameter values are as follows: r
=
3, ω
=
0.2
,
k
=
0.2
,
d
1
=
0.1, d
2
=
0.2
,
d
3
=
0.15, and a
=
5
.(a) b
=
7
,
τ
=
0.2, T
:
1
→
3.8
,
p
:
0.2
→
0.6; (b) b
=
7
,
p
=
0.25, T
:
3.5
→
4.5, τ
:
0.1
→
0.5; (c) b
=
7
,
T
=
4, p
:
0.1
→
0.5
,
τ
:
0.1
→
0.5; (d) τ
=
0.2
,
T
=
3
,
b
:
5
→
8, p
:
0.1
→
0.4.
(a)
(b)
(c)
(d)Complex dynamic properties of system (1) with r
=
3
,
τ
=
0.5
,
ω
=
0.2, μ
=
4
,
β
=
4
,
k
=
0.2
,
p
0
=
0.1, d
1
=
0.1
,
d
2
=
0.2, d
3
=
0.15
,
a
=
5
,
b
=
7, and T
=
4, where killing rate p varies from 0.45 to 0.95.
(a)
(b)
(c)Time-series of the system (1) with pulse, where killing rate p
=
0.7; other parameters are the same as of those in Figure 4.
(a)
(b)
(c)Bifurcation diagrams of system (1) with ω
=
0.2
,
β
=
4
,
μ
=
4
,
k
=
0.2, p
0
=
0.1
,
p
=
0.6
,
d
1
=
0.1, d
2
=
0.2, d
3
=
0.15
,
τ
=
0.05, a
=
5
,
b
=
7, and T
=
4, where the birth rate r varies from 8 to 17.
(a)
(b)
(c)For system (22), according to Theorem 6, we know that, if R
1
*
<
1, the mature prey-extinction periodic solution is globally attractive. This can be seen clearly from Figure 7, where r
=
5
,
τ
=
0.5, ω
=
0.2
,
β
=
0.6
,
k
=
0.2, d
1
=
0.15
,
d
2
=
0.2
,
d
3
=
0.1, p
0
=
0.2
,
p
=
0.5, a
=
3
,
b
=
6
,
T
=
2, l
=
0.7, and μ
=
2; then R
1
*
≈
0.8733
<
1, while if we chose the parameters set as those in Figure 8, it is clear that the system (22) is permanent. By using the same methods as those for system (1), we can investigate the effects of killing rate p, pulse periodic T, birth rate r of pest populations, and the biggest birth rate b of predator on the threshold value R
1
* and discuss their biological implications, so we omitted them here.Dynamical behavior of system (22); the parameter values are as follows: r
=
5
,
τ
=
0.5
,
ω
=
0.2, β
=
0.6
,
k
=
0.2
,
d
1
=
0.15, d
2
=
0.2
,
d
3
=
0.1, p
0
=
0.2
,
p
=
0.5, a
=
3
,
b
=
6
,
T
=
2
,
l
=
0.7, and μ
=
2. Here, R
1
*
≈
0.8733, p
<
1
-
(
(
a
e
d
3
T
)
/
(
a
+
b
)
)
≈
0.7181.
(a)
(b)
(c)Dynamical behavior of system (22); the parameter values are as follows: r
=
5
,
τ
=
0.5
,
ω
=
0.2, β
=
0.6
,
k
=
0.2
,
d
1
=
0.15, d
2
=
0.2
,
d
3
=
0.1, p
0
=
0.2
,
p
=
0.5, a
=
1
,
b
=
2
,
T
=
2
,
l
=
0.5, and μ
=
3.
(a)
(b)
(c)We are more interested in how different patterns of insecticide applications affect the two threshold values (R
1
,
R
1
*) and which strategy is more conducive to pest control. To do this, we investigate the superiority of the two strategies. Assume that the R
1 and R
1
* are with the same parameters; Figure 9 provides the details of how different control strategies affect the threshold values R
1 and R
1
*. If p
=
0, that is, do not spray insecticide, R
1
=
R
1
*, the effect of two control strategies is; if p
≠
0, then R
1
>
R
1
*; this result indicates that if natural enemies were born at the same time as spraying insecticide, the side-effects of pesticides on the natural enemy population are very strong. Obviously, part of the newborn natural enemies will be killed and the number of natural enemies will be sharply reduced; the contribution of chemical control for the threshold value R
1 is stronger than R
1
*. From the practical point of view, we should avoid to spray insecticide when the predators are born.Figure 9
Comparing two kinds of control strategy; the parameter values are as follows:r
=
4.5
,
τ
=
0.2
,
ω
=
0.2,β
=
2
,
d
1
=
0.1
,
d
3
=
0.15, a
=
2
,
b
=
4, and T
=
3; the killing rate p varies from 0 to 0.45.


## 5.  Discussion

When using integrated pest management as an approach to control insect pests, one must be committed to a long-term strategy. It is well known that pesticides usually act not only on the pest species but also, with even stronger impact, on their natural enemies; then, despite the fact that pesticides kill the target pest, they are simultaneously reducing the population density of natural enemies; the influence of pesticides on natural enemies is even greater than on insect pests. As a result of this indirect effect, treatments can counterintuitively lead to an effective increase of the pest species, even pest breakout again.In this paper, our idea is to contribute to pest control programs with pesticide applications and with birth pulse in natural enemies and to provide some strategies; that is, spraying pesticides and natural enemies were born at the same time and at a different time. The threshold conditions which guarantee the existence and stability of the mature prey-extinction periodic solution are provided. If the integrated control methods cannot completely eradicate the mature prey, the pest population can have outbreaks at different scales. By numerical bifurcation investigations, we found that, when choosing different parameter spaces, two attractors from which the pest population oscillates with different amplitudes can coexist for a wide range of parameters. Moreover, the effects of times of spraying pesticides (or the enemy was born) and control tactics on the threshold conditions were carefully investigated. In particular, the effects of the killing rates of pesticides for pest and natural enemy populations and spraying period on the stability of the mature prey-extinction periodic solutions were discussed. The results imply that the modeling methods described can help the pest control specialist to decide appropriate control strategies and assist management decision-making.



---

*Source: 935923-2014-04-07.xml*