# A Time Delay Predator-Prey System with Three-Stage-Structure

Authors: Qiaoqin Gao; Zhen Jin
Journal: The Scientific World Journal
 (2014)

---

## Abstract

A predator-prey system was studied that has a discrete delay, stage-structure, and Beddington-DeAngelis functional response, where predator species has three stages, immature, mature, and old age stages. By using of Mawhin's continuous theorem of coincidence degree theory, a sufficient condition is obtained for the existence of a positive periodic solution.


---

## Body

## 1. Introduction

The dynamic relationship between the predator and the prey has long been and will continue to be one of the dominant themes in both ecology and mathematical ecology due to its universal existence and importance. The traditional predator-prey model has been studied extensively (see [1–3] and references cited therein). Since the pioneering and original work of Aiello and Freedman (see [4]) considered the difference between immature population and mature population, stage-structured models have attracted great attention and have been extensively researched in recent years (see, e.g., [5–9] and the references cited therein). They studied the stage-structured model with time delay, cannibalism, impulsive harvesting strategies and response functions, and many interesting results concerning stability of equilibrium state and existence of periodic solutions. At present, many authors considered species model with two stage-structures, but in real life many mammals (such as humans) growth is divided into three stages: childhood, adulthood, and old age; and almost all the insect growths are divided into eggs, larvae, and adults of three stages. Therefore, the three-stage- structure predator-prey model is studied that is more close to the reality. Recently, more and more papers, but still not too much, considered models with three-stage-structure such as [10–12].Arditi and Ginzburg [13] first put forward the following ratio-dependent predator-prey system:
(1)
x
˙
=
x
(
a
-
b
x
)
-
c
x
y
m
y
+
x
,
y
˙
=
y
(
-
d
+
f
x
m
y
+
x
)
.
Later, Fan and Kuang [14] have explored the dynamics of the nonautonomous, spatially homogeneous, and continuous time predator-prey system with the Beddington-DeAngelis functional response in a more general form:
(2)
x
˙
=
x
(
a
(
t
)
-
b
(
t
)
x
)
-
c
(
t
)
x
y
α
(
t
)
+
β
(
t
)
x
+
γ
(
t
)
y
,
y
˙
=
y
(
-
d
(
t
)
+
f
(
t
)
x
y
α
(
t
)
+
β
(
t
)
x
+
γ
(
t
)
y
)
.
Recently, Chen et al. [8] studied a delayed predator-prey system with Holling II type response and stage-structure for predator of the form
(3)
x
˙
(
t
)
=
x
(
t
)
(
r
(
t
)
-
a
(
t
)
x
(
t
-
τ
1
(
t
)
)
)
-
b
(
t
)
x
(
t
)
y
2
(
t
)
1
+
m
x
(
t
)
,
y
˙
1
(
t
)
=
c
(
t
)
x
(
t
-
τ
2
(
t
)
)
y
2
(
t
-
τ
2
(
t
)
)
1
+
m
x
(
t
-
τ
2
(
t
)
)
-
(
D
(
t
)
+
v
1
(
t
)
)
y
1
(
t
)
,
y
˙
2
(
t
)
=
D
(
t
)
y
1
(
t
)
-
v
2
(
t
)
y
2
(
t
)
.In the present paper, motivated by the above work, we investigate a time delay predator-prey system with three-stage-structure for predator and Beddington-DeAngelis functional response of the form(4)
x
˙
1
(
t
)
=
d
(
t
)
y
(
t
-
τ
2
(
t
)
)
x
2
(
t
-
τ
2
(
t
)
)
α
(
t
)
+
β
(
t
)
y
(
t
-
τ
2
(
t
)
)
+
γ
(
t
)
x
2
(
t
-
τ
2
(
t
)
)
-
(
d
1
(
t
)
+
r
1
(
t
)
)
x
1
(
t
)
,
x
˙
2
(
t
)
=
r
1
(
t
)
x
1
(
t
)
-
d
2
(
t
)
x
2
(
t
)
-
r
2
(
t
)
x
2
(
t
)
-
e
2
(
t
)
x
2
2
(
t
)
,
x
˙
3
(
t
)
=
r
2
(
t
)
x
2
(
t
)
-
d
3
(
t
)
x
3
(
t
)
,
y
˙
(
t
)
=
y
(
t
)
(
c
(
t
)
x
2
(
t
)
α
(
t
)
+
β
(
t
)
y
(
t
)
+
γ
(
t
)
x
2
(
t
)
b
(
t
)
-
a
(
t
)
x
(
t
-
τ
1
(
t
)
)
h
h
h
h
h
-
c
(
t
)
x
2
(
t
)
α
(
t
)
+
β
(
t
)
y
(
t
)
+
γ
(
t
)
x
2
(
t
)
)
,
with initial conditions
(5)
x
i
(
θ
)
=
φ
i
(
θ
)
>
0
(
i
=
1,2
,
3
)
,
y
(
θ
)
=
φ
4
(
θ
)
>
0
,
for
θ
∈
[
-
τ
,
0
]
,
x
i
(
0
)
>
0
,
i
=
1,2
,
3
,
y
(
0
)
>
0
,
where x
1
(
t
), x
2
(
t
), and x
3
(
t
) are the densities of immature, mature, and the old age predator at time t, respectively; y
(
t
) represents the density of prey at time t. τ
=
max
⁡
t
∈
[
0
,
ω
]
{
τ
1
(
t
)
,
τ
2
(
t
)
}, τ
1
(
t
)
≥
0 is delay due to prey densities and τ
2
(
t
)
≥
0 is delay due to gestation of predator. The death rates of the immature, mature, and old age predator are proportional to the fact that existing immature, mature, and old age predator population with respective proportionality d
1
(
t
), d
2
(
t
), d
3
(
t
), r
1
(
t
) denote the rate of transformation from immature predator into mature predator; r
2
(
t
) is the rate of transformation from mature predator into old age predator; k
(
t
)
=
(
d
(
t
)
/
c
(
t
)
) (0
<
k
(
t
)
≤
1) is the transformation coefficient from prey into the immature predators. Let d
1
(
t
)
+
r
1
(
t
)
=
m
(
t
) and d
2
(
t
)
+
r
2
(
t
)
=
n
(
t
).In this paper, we always make the following assumption for the system (4).(H1)
The coefficientsa
(
t
), b
(
t
), c
(
t
), d
(
t
), d
i
(
t
) (i
=
1,2
,
3), α
(
t
), β
(
t
), γ
(
t
), e
2
(
t
) and r
i
(
t
) (i
=
1,2) are all positive ω-periodic continuous function in R
=
[
0
,
∞
).The organization of this paper is as follows. In the next section, we give two lemmas for preparing the study of the existence of positive periodic solutions. In Section3, by using Gains and Mawhin’s continuation theorem of coincidence degree theory, a sufficient condition on the existence of positive periodic solutions for the model is obtained.


## 2. Preliminaries

For ecological reasons, we consider the initial function of the system (4) only in R
+
4
×
R
+
1
=
{
(
x
1
(
t
)
,
x
2
(
t
)
,
x
3
(
t
)
,
y
(
t
)
)
T
:
x
i
≥
0
,
y
≥
0
,
i
=
1,2
,
3
}
×
[
-
τ
,
0
].Supposing thatf
(
t
), t
∈
[
0
,
∞
) is a continuous function with periodic ω, we denote
(6)
f
L
=
min
⁡
t
∈
[
0
,
ω
]
{
f
(
t
)
}
,
f
M
=
max
⁡
t
∈
[
0
,
ω
]
{
f
(
t
)
}
,
f
-
=
1
ω
∫
0
ω
‍
f
(
t
)
d
t
.Using the above definitions, the coefficients of the system (4) not only are all positive and periodic function with common periodic solution but also satisfy the following two items and hold true:
(7)
min
⁡
i
=
1,2
,
3
{
a
L
,
b
L
,
c
L
,
d
L
,
m
L
,
n
L
,
d
i
L
,
r
i
L
,
e
2
L
,
α
L
,
β
L
,
γ
L
}
>
0
,
max
⁡
i
=
1,2
,
3
{
a
M
,
b
M
,
c
M
,
d
M
,
m
M
,
n
M
,
d
i
M
,
r
i
M
,
e
2
M
,
α
M
,
β
M
,
γ
M
}
<
+
∞
.In order to obtain the existence of at least one periodic solution of the system (4), the method used the applications of the continuation theorem of coincidence degree. For convenience, we introduce a few concepts and results about the coincidence degree as follows.LetX, Y be real Banach spaces, let L
:
Dom
⁡
L
⊂
X
→
Z be a linear mapping, and let N
:
X
→
Y be a continuous mapping. The mapping L will be called a Fredholm mapping of index zero if dim
⁡
⁡
Ker
⁡
L
=
co
⁡
dim
⁡
⁡
Im
⁡
L
<
∞ and Im
⁡
L is closed in Y. If L a is Fredholm mapping of index zero, then there exist continuous projectors P
:
X
→
X and Q
:
Z
→
Z such that Im
⁡
P
=
Ker
⁡
L and Ker
⁡
Q
=
Im
⁡
L
=
Im
⁡
(
I
-
Q
). It follows that the restriction L
|
Dom
⁡
L
∩
Ker
⁡
P
:
(
I
-
P
)
X
→
Im
⁡
L is invertible. Denote the inverse of L
|
Dom
⁡
L
∩
Ker
⁡
P by K
P.LetΩ be an open bounded subset of X and denote the closure of Ω by Ω
¯. A mapping N is said to be L-compact on Ω
¯. If Q
N
(
Ω
¯
) is bounded and K
P
(
I
-
Q
)
N
:
Ω
¯
→
X is compact, since Im
⁡
Q is isomorphic to KerL, then there exists an isomorphism J
:
Im
⁡
Q
→
Ker
⁡
L.Lemma 1 (Mawhin’s continuation theorem [15, page 40]).
LetX, Y be two Banach spaces, let L be a Fredholm mapping of index zero and let N
:
Ω
¯
→
Y be a continuous operator which is L-compact on Ω
¯. Assume that,(a)
for eachλ
∈
(
0,1
), x
∈
∂
Ω
∩
D
o
m
L, L
x
≠
λ
N
x;
(b)
for eachx
∈
∂
Ω
∩
K
e
r
L, Q
N
x
≠
0;
(c)
d
e
g
(
J
Q
N
,
Ω
∩
K
e
r
L
,
0
)
≠
0.Then the operator equation L
x
=
N
x has at least one solution in Ω
¯
∩
D
o
m
L.Lemma 2.
i
n
t
R
+
4 is a positively invariant region of the system (4).


## 3. Existence of Periodic Solution

Theorem 3.
Assume that the system (4) satisfies (H1) and the following:(H2)
r
1
M
d
M
>
n
L
m
L
β
L;
(H3)
b
-
>
(
c
/
r
)
¯;
(H4)
r
1
L
d
L
A
4
>
m
M
n
M
(
α
M
+
β
M
A
4
).Then the system (4) with initial condition (5) has at least one ω-periodic solution.Proof.
Let(8)
x
i
(
t
)
=
e
u
i
(
t
)
(
i
=
1,2
,
3
)
,
y
(
t
)
=
e
u
4
(
t
)
.
Then the system (4) transforms correspondingly to the following equations:
(9)
u
˙
1
(
t
)
=
d
(
t
)
e
u
4
(
t
-
τ
2
(
t
)
)
e
u
2
(
t
-
τ
2
(
t
)
)
-
u
1
(
t
)
α
(
t
)
+
β
(
t
)
e
u
4
(
t
-
τ
2
(
t
)
)
+
γ
(
t
)
e
u
2
(
t
-
τ
2
(
t
)
)
-
m
(
t
)
,
u
˙
2
(
t
)
=
r
1
(
t
)
e
u
1
(
t
)
-
u
2
(
t
)
-
n
(
t
)
-
e
2
(
t
)
e
u
2
(
t
)
,
u
˙
3
(
t
)
=
r
2
(
t
)
e
u
2
(
t
)
-
u
3
(
t
)
-
d
3
(
t
)
,
u
˙
4
(
t
)
=
b
(
t
)
-
a
(
t
)
e
u
4
(
t
-
τ
1
(
t
)
)
-
c
(
t
)
e
u
2
(
t
)
α
(
t
)
+
β
(
t
)
e
u
4
(
t
)
+
γ
(
t
)
e
u
2
(
t
)
.
It is easy to see that if the system (9) has a ω-periodic solution u
*
(
t
)
=
(
u
1
*
,
u
2
*
,
u
3
*
,
u
4
*
)
T, then the system (4) has correspondingly a positive ω-periodic solution (
x
1
*
,
x
2
*
,
x
3
*
,
y
*
)
T
=
(
e
u
1
*
,
e
u
2
*
,
e
u
3
*
,
e
u
4
*
)
T. Thus, to complete the proof, suffice it to say that the system (9) has at least one positive ω-periodic solution.
DefineX
=
Y
=
{
u
(
t
)
=
(
u
1
(
t
)
,
u
2
(
t
)
,
u
3
(
t
)
,
u
4
(
t
)
)
T
∈
C
(
R
,
R
4
)
,
u
(
t
+
ω
)
=
u
(
t
)
}, and ∥
u
∥
=
∑
i
=
1
4
max
⁡
t
∈
[
0
,
ω
]
|
u
i
(
t
)
|; then X and Y are Banach spaces with the norm ∥
·
∥.
LetN
u
=
(
N
1
(
t
)
,
N
2
(
t
)
,
N
3
(
t
)
,
N
4
(
t
)
)
T, u
∈
X, L
u
=
u
˙
(
t
)
=
d
u
(
t
)
/
d
t, P
u
=
(
1
/
ω
)
∫
0
ω
‍
u
(
t
)
d
t, u
∈
X, and Q
y
=
(
1
/
ω
)
∫
0
ω
‍
y
(
t
)
d
t, y
∈
Y, where
(10)
N
1
(
t
)
=
d
(
t
)
e
u
4
(
t
-
τ
2
(
t
)
)
e
u
2
(
t
-
τ
2
(
t
)
)
-
u
1
(
t
)
α
(
t
)
+
β
(
t
)
e
u
4
(
t
-
τ
2
(
t
)
)
+
γ
(
t
)
e
u
2
(
t
-
τ
2
(
t
)
)
-
m
(
t
)
,
N
2
(
t
)
=
r
1
(
t
)
e
u
1
(
t
)
-
u
2
(
t
)
-
n
(
t
)
-
e
2
(
t
)
e
u
2
(
t
)
,
N
3
(
t
)
=
r
2
(
t
)
e
u
2
(
t
)
-
u
3
(
t
)
-
d
3
(
t
)
,
N
4
(
t
)
=
b
(
t
)
-
a
(
t
)
e
u
4
(
t
-
τ
1
(
t
)
)
-
c
(
t
)
e
u
2
(
t
)
α
(
t
)
+
β
(
t
)
e
u
4
(
t
)
+
γ
(
t
)
e
u
2
(
t
)
.
Thus, the system (9) can be written in the form L
u
=
N
u,  u
∈
X.
Obviously,Ker
⁡
L
=
R
4, Im
⁡
L
=
{
u
(
t
)
∈
Z
,
∫
0
ω
‍
u
(
t
)
d
t
=
0
} is closed in Y and dim
⁡
⁡
Ker
⁡
L
=
4
=
Co
dim
⁡
⁡
Im
⁡
L and P, Q are continuous projectors such that Im
⁡
P
=
Ker
⁡
L, Ker
⁡
Q
=
Im
⁡
L
=
Im
⁡
(
I
-
Q
). Therefore, L is a Fredholm mapping of index zero. Furthermore, the inverse K
P
:
Im
⁡
L
→
Ket
P
∩
Dom
⁡
L exists and has the form
(11)
K
P
(
u
)
=
∫
0
t
‍
u
(
s
)
d
s
-
1
ω
∫
0
ω
‍
∫
0
t
‍
u
(
s
)
d
s
d
t
.
Thus
(12)
Q
N
u
=
1
ω
∫
0
ω
‍
N
(
t
)
d
t
,
K
P
(
I
-
Q
)
N
u
=
∫
0
t
‍
N
(
s
)
d
s
-
1
ω
∫
0
ω
‍
∫
0
t
‍
N
u
(
s
)
d
s
d
t
-
(
t
ω
-
1
2
)
∫
0
ω
‍
N
u
(
s
)
d
s
.
Obviously, Q
N and K
P
(
I
-
Q
)
N are continuous by the Lebesgue theorem; it is not difficult to show that Q
N
(
Ω
¯
) is bounded and K
P
(
I
-
Q
)
N
(
Ω
¯
)
¯ is compact for any open bounded set Ω
⊂
X by the Arzela-Ascoli theorem. Hence N is L-compact on Ω
¯ for any open bounded set Ω
⊂
X.
In order to apply Lemma1, we need to search for an appropriate open and bounded subset Ω
⊂
X. Corresponding to the operator equation, L
u
=
λ
N
u, λ
∈
(
0,1
), we have
(13)
u
˙
1
(
t
)
=
λ
(
d
(
t
)
e
u
4
(
t
-
τ
2
(
t
)
)
e
u
2
(
t
-
τ
2
(
t
)
)
-
u
1
(
t
)
α
(
t
)
+
β
(
t
)
e
u
4
(
t
-
τ
2
(
t
)
)
+
γ
(
t
)
e
u
2
(
t
-
τ
2
(
t
)
)
-
m
(
t
)
)
,
u
˙
2
(
t
)
=
λ
(
r
1
(
t
)
e
u
1
(
t
)
-
u
2
(
t
)
-
n
(
t
)
-
e
2
(
t
)
e
u
2
(
t
)
)
,
u
˙
3
(
t
)
=
λ
(
r
2
(
t
)
e
u
2
(
t
)
-
u
3
(
t
)
-
d
3
(
t
)
)
,
u
˙
4
(
t
)
=
λ
(
c
(
t
)
e
u
2
(
t
)
α
(
t
)
+
β
(
t
)
e
u
4
(
t
)
+
γ
(
t
)
e
u
2
(
t
)
b
(
t
)
-
a
(
t
)
e
u
4
(
t
-
τ
1
(
t
)
)
h
h
h
-
c
(
t
)
e
u
2
(
t
)
α
(
t
)
+
β
(
t
)
e
u
4
(
t
)
+
γ
(
t
)
e
u
2
(
t
)
)
.
Suppose thatu
=
u
(
t
)
∈
X is a solution to (13) for a certain λ
∈
(
0,1
). Integrating (13) over the interval [
0
,
ω
] leads to
(14)
∫
0
ω
‍
d
(
t
)
e
u
4
(
t
-
τ
2
(
t
)
)
e
u
2
(
t
-
τ
2
(
t
)
)
-
u
1
(
t
)
α
(
t
)
+
β
(
t
)
e
u
4
(
t
-
τ
2
(
t
)
)
+
γ
(
t
)
e
u
2
(
t
-
τ
2
(
t
)
)
d
t
=
m
-
ω
,
(15)
∫
0
ω
‍
r
1
(
t
)
e
u
1
(
t
)
-
u
2
(
t
)
d
t
=
n
-
ω
+
∫
0
ω
‍
e
2
(
t
)
e
u
2
(
t
)
d
t
,
(16)
∫
0
ω
‍
r
2
(
t
)
e
u
2
(
t
)
-
u
3
(
t
)
d
t
=
d
-
3
ω
,
(17)
∫
0
ω
‍
[
a
(
t
)
e
u
4
(
t
-
τ
1
(
t
)
)
+
c
(
t
)
e
u
2
(
t
)
α
(
t
)
+
β
(
t
)
e
u
4
(
t
)
+
γ
(
t
)
e
u
2
(
t
)
]
d
t
=
b
-
ω
.
It follows from (14) to (17) that
(18)
∫
0
ω
‍
|
u
˙
1
(
t
)
|
d
t
=
λ
∫
0
ω
‍
|
d
(
t
)
e
u
4
(
t
-
τ
2
(
t
)
)
e
u
2
(
t
-
τ
2
(
t
)
)
-
u
1
(
t
)
α
(
t
)
+
β
(
t
)
e
u
4
(
t
-
τ
2
(
t
)
)
+
γ
(
t
)
e
u
2
(
t
-
τ
2
(
t
)
)
h
h
h
h
h
-
m
(
t
)
d
(
t
)
e
u
4
(
t
-
τ
2
(
t
)
)
e
u
2
(
t
-
τ
2
(
t
)
)
-
u
1
(
t
)
α
(
t
)
+
β
(
t
)
e
u
4
(
t
-
τ
2
(
t
)
)
+
γ
(
t
)
e
u
2
(
t
-
τ
2
(
t
)
)
|
d
t
≤
2
m
-
ω
=
:
l
1
,
(19)
∫
0
ω
‍
|
u
˙
2
(
t
)
|
d
t
=
λ
∫
0
ω
‍
|
r
1
(
t
)
e
u
1
(
t
)
-
u
2
(
t
)
-
n
(
t
)
-
e
2
(
t
)
e
u
2
(
t
)
|
d
t
≤
2
n
-
ω
+
2
∫
0
ω
‍
e
2
(
t
)
e
u
2
(
t
)
d
t
,
(20)
∫
0
ω
‍
|
u
˙
3
(
t
)
|
d
t
=
λ
∫
0
ω
‍
|
r
2
(
t
)
e
u
2
(
t
)
-
u
3
(
t
)
-
d
3
(
t
)
|
d
t
≤
2
d
-
3
ω
=
:
l
3
,
(21)
∫
0
ω
‍
|
u
˙
4
(
t
)
|
d
t
=
λ
∫
0
ω
‍
|
-
c
(
t
)
e
u
2
(
t
)
α
(
t
)
+
β
(
t
)
e
u
4
(
t
)
+
γ
(
t
)
e
u
2
(
t
)
b
(
t
)
-
a
(
t
)
e
u
4
(
t
-
τ
1
(
t
)
)
h
h
h
h
h
-
c
(
t
)
e
u
2
(
t
)
α
(
t
)
+
β
(
t
)
e
u
4
(
t
)
+
γ
(
t
)
e
u
2
(
t
)
|
d
t
≤
2
b
-
ω
=
:
l
4
.
Multiplying the first equation of (13) by e
u
1
(
t
) and integrating it over [
0
,
ω
] lead to
(22)
∫
0
ω
‍
m
(
t
)
e
u
1
(
t
)
d
t
=
∫
0
ω
‍
d
(
t
)
e
u
4
(
t
-
τ
2
(
t
)
)
e
u
2
(
t
-
τ
2
(
t
)
)
α
(
t
)
+
β
(
t
)
e
u
4
(
t
-
τ
2
(
t
)
)
+
γ
(
t
)
e
u
2
(
t
-
τ
2
(
t
)
)
d
t
,
which implies
(23)
m
L
∫
0
ω
‍
e
u
1
(
t
)
d
t
≤
∫
0
ω
‍
d
(
t
)
β
(
t
)
e
u
2
(
t
-
τ
2
(
t
)
)
d
t
≤
d
M
β
L
∫
0
ω
‍
e
u
2
(
t
-
τ
2
(
t
)
)
d
t
,
∫
0
ω
‍
e
u
2
(
t
)
d
t
=
∫
τ
2
(
s
)
ω
-
τ
2
(
s
)
‍
e
u
2
(
t
)
d
t
=
∫
0
ω
‍
e
u
2
(
t
-
τ
2
(
t
)
)
d
t
≥
m
L
β
L
d
M
∫
0
ω
‍
e
u
1
(
t
)
d
t
.
Similarly, multiplying the second equation of (13) by e
u
2
(
t
) and integrating it over [
0
,
ω
] give
(24)
∫
0
ω
‍
e
2
(
t
)
e
2
u
2
(
t
)
d
t
+
∫
0
ω
‍
n
(
t
)
e
u
2
(
t
)
d
t
=
∫
0
ω
‍
r
1
(
t
)
e
u
1
(
t
)
d
t
which yields
(25)
e
2
L
∫
0
ω
‍
e
2
u
2
(
t
)
d
t
+
n
L
∫
0
ω
‍
e
u
2
(
t
)
d
t
≤
r
1
M
∫
0
ω
‍
e
u
1
(
t
)
d
t
.
By using the inequalities (
∫
0
ω
‍
e
u
2
(
t
)
d
t
)
2
≤
ω
∫
0
ω
‍
e
2
u
2
(
t
)
d
t and from (23) and (25) as well as (H1), we have
(26)
∫
0
ω
‍
e
u
2
(
t
)
d
t
≥
m
L
β
L
d
M
∫
0
ω
‍
e
u
1
(
t
)
d
t
≥
m
L
β
L
d
M
r
1
M
(
e
2
L
∫
0
ω
‍
e
2
u
2
(
t
)
d
t
+
n
L
∫
0
ω
‍
e
u
2
(
t
)
d
t
)
≥
m
L
β
L
d
M
r
1
M
[
e
2
L
ω
(
∫
0
ω
‍
e
u
2
(
t
)
d
t
)
2
+
n
L
∫
0
ω
‍
e
u
2
(
t
)
d
t
]
.
Hence
(27)
∫
0
ω
‍
e
u
2
(
t
)
d
t
≤
ω
d
M
r
1
M
-
n
L
m
L
β
L
m
L
β
L
e
2
L
=
:
ω
B
2
.
From (23) and (27), we have
(28)
∫
0
ω
‍
e
u
1
(
t
)
d
t
≤
d
M
m
L
β
L
∫
0
ω
‍
e
u
2
(
t
)
d
t
≤
ω
B
2
d
M
m
L
β
L
=
:
ω
B
1
.
Similarly, multiplying the third equation of (13) by e
u
3
(
t
) and integrating it over [
0
,
ω
] lead to
(29)
∫
0
ω
‍
d
3
(
t
)
e
u
3
(
t
)
d
t
=
∫
0
ω
‍
r
2
(
t
)
e
u
2
(
t
)
d
t
.
And, together with (27), we obtain
(30)
∫
0
ω
‍
e
u
3
(
t
)
d
t
≤
r
2
M
d
3
L
∫
0
ω
‍
e
u
2
(
t
)
d
t
≤
ω
B
2
r
2
M
d
3
L
=
:
ω
B
3
.
Integrating the fourth equation of (13) over [
0
,
ω
] leads to
(31)
∫
0
ω
‍
a
(
t
)
e
u
4
(
t
-
τ
1
(
t
)
)
d
t
=
∫
0
ω
‍
b
(
t
)
d
t
-
∫
0
ω
‍
c
(
t
)
e
u
2
(
t
)
α
(
t
)
+
β
(
t
)
e
u
4
(
t
)
+
γ
(
t
)
e
u
2
(
t
)
d
t
,
which implies
(32)
∫
0
ω
‍
e
u
4
(
t
)
d
t
≤
ω
b
M
a
L
=
:
ω
B
4
.
Since u
(
t
)
∈
X, there exist ξ
i, η
i
∈
[
0
,
ω
] such that
(33)
u
i
(
ξ
i
)
=
max
⁡
t
∈
[
0
,
ω
]
u
i
(
t
)
,
u
i
(
η
i
)
=
min
⁡
t
∈
[
0
,
ω
]
u
i
(
t
)
,
i
=
1,2
,
3
,
in which, together with (27), (28), (30), and (32), we deduce
(34)
u
i
(
ξ
i
)
≤
ln
⁡
B
i
,
i
=
1,2
,
3,4
.
In view of (19) and (27), we have
(35)
∫
0
ω
‍
|
u
˙
2
(
t
)
|
d
t
≤
2
n
-
ω
+
2
ω
e
2
M
B
2
=
:
l
2
.
Keeping in mind (18), (20), (21), and (34), we obtain that for t
∈
[
0
,
ω
]
(36)
u
i
(
t
)
≤
u
i
(
ξ
i
)
+
∫
0
ω
‍
|
u
˙
i
(
t
)
|
d
t
≤
ln
⁡
B
i
+
l
i
,
i
=
1,2
,
3,4
.
From (17), (33), and (H3), we have
(37)
a
M
ω
e
u
4
(
η
4
)
≥
a
M
∫
0
ω
‍
e
u
4
(
t
)
d
t
=
a
M
∫
0
ω
‍
e
u
4
(
t
-
τ
1
(
t
)
)
d
t
=
∫
0
ω
‍
b
(
t
)
d
t
-
∫
0
ω
‍
c
(
t
)
e
u
2
(
t
)
α
(
t
)
+
β
(
t
)
e
u
4
(
t
)
+
γ
(
t
)
e
u
2
(
t
)
d
t
≥
b
-
ω
-
(
c
γ
)
¯
ω
.
Set A
4
=
(
b
-
-
(
c
/
γ
)
¯
)
/
a
M, and hence
(38)
e
u
4
(
η
4
)
≥
A
4
.
Equation (33) leads to
(39)
d
(
η
1
)
e
u
4
(
η
1
-
τ
2
(
η
1
)
)
e
u
2
(
η
1
-
τ
2
(
η
1
)
)
-
u
1
(
η
1
)
α
(
η
1
)
+
β
(
η
1
)
e
u
4
(
η
1
-
τ
2
(
η
1
)
)
+
γ
(
η
1
)
e
u
2
(
η
1
-
τ
2
(
η
1
)
)
-
m
(
η
1
)
=
0
,
r
1
(
η
2
)
e
u
1
(
η
2
)
-
u
2
(
η
2
)
-
n
(
η
2
)
-
e
2
(
η
2
)
e
u
2
(
η
2
)
=
0
,
r
2
(
η
3
)
e
u
2
(
η
3
)
-
u
3
(
η
3
)
-
d
3
(
η
3
)
=
0
,
b
(
η
4
)
-
a
(
η
4
)
e
u
4
(
η
4
-
τ
1
(
η
4
)
)
-
c
(
η
4
)
e
u
2
(
η
4
)
α
(
η
4
)
+
β
(
η
4
)
e
u
4
(
η
4
)
+
γ
(
η
4
)
e
u
2
(
η
4
)
=
0
.
Taking into account (38) and the monotonicity of g
(
x
)
=
x
/
(
x
+
l
) leads to
(40)
m
M
e
u
1
(
η
1
)
≥
m
(
η
1
)
e
u
1
(
η
1
)
=
d
(
η
1
)
e
u
4
(
η
1
-
τ
2
(
η
1
)
)
e
u
2
(
η
1
-
τ
2
(
η
1
)
)
-
u
1
(
η
1
)
α
(
η
1
)
+
β
(
η
1
)
e
u
4
(
η
1
-
τ
2
(
η
1
)
)
+
γ
(
η
1
)
e
u
2
(
η
1
-
τ
2
(
η
1
)
)
≥
d
L
A
4
e
u
2
(
η
2
)
α
M
+
β
M
A
4
+
γ
M
e
u
2
(
η
2
)
,
(41)
e
2
M
e
2
u
2
(
η
2
)
+
n
M
e
u
2
(
η
2
)
≥
e
2
(
η
2
)
e
2
u
2
(
η
2
)
+
n
(
η
2
)
e
u
2
(
η
2
)
=
r
1
(
η
2
)
e
u
1
(
η
2
)
≥
r
1
L
e
u
1
(
η
1
)
,
(42)
d
3
M
e
u
3
(
η
3
)
≥
d
3
(
η
3
)
e
u
3
(
η
3
)
=
r
2
(
η
3
)
e
u
2
(
η
3
)
≥
r
2
L
e
u
2
(
η
2
)
.
It follows from (40) and (41) that
(43)
e
2
M
e
u
2
(
η
2
)
+
n
M
≥
r
1
L
d
L
A
4
m
M
(
α
M
+
β
M
A
4
+
γ
M
e
u
2
(
η
2
)
)
,
which yields
(44)
e
u
2
(
η
2
)
≥
(
(
e
2
M
(
α
M
+
β
M
A
4
)
+
n
M
γ
M
)
2
+
4
A
5
-
(
e
2
M
(
α
M
+
β
M
A
4
)
+
n
M
γ
M
)
h
h
+
(
e
2
M
(
α
M
+
β
M
A
4
)
+
n
M
γ
M
)
2
+
4
A
5
)
×
(
2
e
2
M
γ
M
)
-
1
=
:
A
2
,
where A
5
=
e
2
M
γ
M
(
(
r
1
L
d
L
A
4
/
m
M
)
-
n
M
(
α
M
+
β
M
A
4
)
). From (40), (42), and (44), we have
(45)
e
u
1
(
η
1
)
≥
d
L
A
4
A
2
m
M
(
α
M
+
β
M
A
4
+
γ
M
A
2
)
=
:
A
1
,
e
u
3
(
η
3
)
≥
r
2
L
A
2
d
3
M
=
:
A
3
,
which, together with (38), gives
(46)
u
i
(
η
i
)
≥
ln
⁡
⁡
A
i
,
i
=
1,2
,
3,4
.
Moreover, from (18) to (21) and (35), we obtain
(47)
u
i
(
t
)
≥
u
i
(
η
i
)
-
∫
0
ω
|
u
˙
i
(
t
)
|
d
t
≥
ln
⁡
⁡
A
i
-
l
i
,
t
∈
[
0
,
ω
]
,
i
=
1,2
,
3,4
.
It follows from (36) and (47) that
(48)
max
⁡
t
∈
[
0
,
ω
]
⁡
|
u
i
(
t
)
|
<
max
⁡
⁡
{
|
ln
⁡
B
i
+
l
i
|
,
|
ln
⁡
⁡
A
i
-
l
i
|
}
=
R
i
,
max
⁡
t
∈
[
0
,
ω
]
⁡
|
u
i
(
t
)
|
<
max
⁡
⁡
{
|
ln
⁡
B
i
+
l
i
|
,
|
-
l
i
|
}
i
=
1,2
,
3,4
.
Clearly,  R
i (i
=
1,2
,
3,4) is independent of λ. From condition of Theorem 3, it is easy to know that there exist points ζ
i
∈
[
0
,
ω
], i
=
1,2, such that the algebraic equations,
(49)
d
-
e
u
4
e
u
2
-
u
1
α
(
ζ
1
)
+
β
(
ζ
1
)
e
u
4
+
γ
(
ζ
1
)
e
u
2
-
m
¯
=
0
,
r
1
¯
e
u
1
-
u
2
-
n
-
-
e
2
¯
e
u
2
=
0
,
r
2
¯
e
u
2
-
u
3
-
d
3
¯
=
0
,
b
-
-
a
-
e
u
4
-
c
-
e
u
2
α
(
ζ
2
)
+
β
(
ζ
2
)
e
u
4
+
γ
(
ζ
2
)
e
u
2
=
0
,
have a unique solution (
v
1
*
,
v
2
*
,
v
3
*
,
v
4
*
)
T.
SetH
=
∑
i
=
1
4
R
i
+
R
0, where R
0 is taken sufficiently large such that the unique solution to (49) satisfies ∥
(
v
1
*
,
v
2
*
,
v
3
*
,
v
4
*
)
T
∥
=
|
v
1
*
|
+
|
v
2
*
|
+
|
v
3
*
|
+
|
v
4
*
|
<
R
0.
We takeΩ
=
{
u
(
t
)
=
(
u
1
(
t
)
,
u
2
(
t
)
,
u
3
(
t
)
,
u
3
(
t
)
)
T
∈
X
:
∥
u
∥
<
M
}; then Ω satisfies the condition (a) of Lemma 1. When u
∈
∂
Ω
∩
ker
⁡
⁡
L
=
∂
Ω
∩
R
4, u is a constant vector in R
4 with ∥
u
∥
=
H; then we have
(50)
Q
N
u
=
[
d
-
e
u
4
e
u
2
-
u
1
α
(
ζ
1
)
+
β
(
ζ
1
)
e
u
4
+
γ
(
ζ
1
)
e
u
2
-
m
¯
r
1
¯
e
u
1
-
u
2
-
n
-
-
e
2
¯
e
u
2
r
2
¯
e
u
2
-
u
3
-
d
3
¯
b
-
-
a
-
e
u
4
-
c
-
e
u
2
α
(
ζ
2
)
+
β
(
ζ
2
)
e
u
4
+
γ
(
ζ
2
)
e
u
2
]
=
[
0
0
0
0
]
.
This proves that condition (b) of Lemma 1 is satisfied.
We define the homotopyJ
:
Dom
⁡
L
×
[
0,1
]
→
X, and a direct calculation shows that
(51)
deg
⁡
⁡
(
J
Q
N
u
,
Ω
∩
Ker
⁡
L
,
(
0,0
,
0,0
)
T
)
=
sgn
⁡
(
a
¯
d
¯
r
1
¯
r
2
¯
e
2
u
4
+
u
2
-
u
3
(
α
(
ζ
1
)
+
β
(
ζ
1
)
e
u
4
+
γ
(
ζ
1
)
e
u
2
)
2
)
=
1
≠
0
.
Therefore Ω satisfies all the conditions of Lemma 1. Thus, by Lemma 1, we conclude that L
u
=
N
u has at least one solution in Dom
⁡
L
∩
Ω
-; that is, the system (9) has at least one ω-periodic solution in Ω. Then, there exists at least one ω-periodic solution for the system (4). The proof is now completed.


## 4. Discussion

In this paper, we have considered a predator-prey system with three-stage-structure for predator, time delay due to predator gestation and prey densities, and Beddington-DeAngelis functional response. We have derived criteria for the existence of a positive periodic solution. The method used to obtain the main result involves application of Mawhin’s continuous theorem of coincidence degree theory and integral inequalities.Chen et al. [8] considered a delayed predator-prey system with Holling II type response and structure for predator. If the predator x
2 is absent and the coefficients α
(
t
)
≡
1, β
(
t
)
≡
m(constant), and γ
(
t
)
≡
0, the corresponding subsystems of the system (4) are system (3); that is, the system of [8] belongs to the case of the system (4). The conditions for the existence of a positive periodic solution of the system (4) remain the same if delay is absent; thus the delay is “harmless” for the existence of a positive periodic solution of the system (4).In a word, despite considering the three-stage-structure, people can more protect the beneficial insect from dying out and under the preconditions; the actualization of the periodic changes is more convenient and realistic. The results accord with the change in the natural environments, so the consideration of stage-structure is necessary and important.



---

*Source: 512838-2014-03-27.xml*