# Stability and Hopf Bifurcation in a Delayed Predator-Prey System with Herd Behavior

Authors: Chaoqun Xu; Sanling Yuan
Journal: Abstract and Applied Analysis
 (2014)

---

## Abstract

A special predator-prey system is investigated in which the prey population exhibits herd behavior in order to provide a self-defense against predators, while the predator is intermediate and its population shows individualistic behavior. Considering the fact that there always exists a time delay in the conversion of the biomass of prey to that of predator in this system, we obtain a delayed predator-prey model with square root functional response and quadratic mortality. For this model, we mainly investigate the stability of positive equilibrium and the existence of Hopf bifurcation by choosing the time delay as a bifurcation parameter.


---

## Body

## 1. Introduction

During the last few decades, there has been great interest in the construction and study of models for the population dynamics of predator-prey systems. The classical predator-prey model can be written in the generalized form (see [1])
(1)
d
N
d
t
=
F
(
N
)
N
-
ϕ
(
N
,
P
)
P
,
d
P
d
t
=
[
k
ϕ
(
N
,
P
)
-
μ
(
P
)
]
P
,
N
(
0
)
≥
0
,
P
(
0
)
≥
0
,
where N
=
N
(
t
), P
=
P
(
t
) denote the number or density of the prey and predator species at time t, respectively, F
(
N
) is the relative growth function of prey species in the absence of predator, μ
(
P
) is the relative mortality function of predator species in the absence of prey, ϕ
(
N
,
P
) is the functional response function of the predator to the prey which can be interpreted as the change in the number of prey species attached per unit time per predator as the prey number changes, and k is the biomass conversion or consumption rate.For system (1), the dynamic behavior is determined by the three crucial components: F
(
N
), μ
(
P
), and ϕ
(
N
,
P
). Chen et al. in [1] summarized some types of F
(
N
) and μ
(
P
) which determine the species intrinsic population dynamics (primary production and mortality, resp.). For example,
(2)
F
(
N
)
=
r
(
1
-
N
K
)
,
F
(
N
)
=
r
K
-
N
K
+
ɛ
N
;
μ
(
P
)
≡
s
,
μ
(
P
)
=
a
b
+
P
.
Functional response function ϕ
(
N
,
P
) also can be classified into many different types: Holling I–IV types [2, 3], Beddington-DeAngelis type [4, 5], Ivlev type [6], ratio-dependent type [7], and the modified forms of these types [8–18]. Over the years, the mathematical models with different F
(
N
), μ
(
P
), and ϕ
(
N
,
P
) have been constructed by many researchers for describing different predator-prey systems; see [2–18].Recently, one kind of predator-prey system with herd behavior is considered by some researchers [19–23]. In this kind of systems, such as herbivores on the large savanna and their large predators, the prey species join together in herds in order to provide a self-defense against predators, and the predator interacts with the prey along the outer corridor of the herd of prey [20, 21]. Therefore, it is more appropriate to model the response functions of prey in terms of the square root of the number of prey species for this class of systems. In [20], Braza considered a predator-prey system with square root functional response and compared the dynamics of this system with the dynamics of predator-prey system that use a typical Lotka-Volterra interaction term.Moreover, there are some papers concerned with the effects of the form of mortality terms for the dynamics [23–26]. Brentnall et al. [24] pointed out that the mortality term of predator is usually described by “linear form” or “quadratic form” and the quadratic mortality is suited to intermediate predator (such as piscivore). Recently, Yuan et al. [23] studied a spatial predator-prey model with herd behavior and concluded that the Turing pattern is induced by quadratic mortality. Fulton et al. in [25] showed that, in the large and interlinked webs used in ecosystem models, model behavior is far more sensitive to the form of the grazing term than to that of the mortality terms that are close to the modeled food web.On the other hand, for most of the natural ecosystems, every species does not respond instantaneously to changes in the environment or the interactions with other species within the community. Thus models with delay are much more realistic [27–31], and many researchers paid a lot of attention to the delayed predator-prey systems; see [32–40] and the references cited therein. For example, Xiao and Ruan [34] discussed a delayed predator-prey system with nonmonotonic functional response and showed that there is a B-T singularity for any time delay value and the system can exhibit Hopf bifurcation as the time delay passes through some critical values. Xu et al. in [36] proposed a delayed periodic L-V type predator-prey system with prey dispersal in two-patch environments and gave the sufficient conditions of the existence, uniqueness, and global stability of positive periodic solution.Motivated by the above, in the present paper, we are intended to consider a predator-prey system in which the prey species exhibits herd behavior, the predator species is intermediate, and there exists a time delay in the convention of the biomass of prey to that of predator.This paper is organized as follows. We first introduce our working system in the next section. In Section3, by analyzing the characteristic equation, we discuss the local stability of a positive equilibrium and study the existence of Hopf bifurcations at the positive equilibrium. In Section 4, using the normal form theory and center manifold argument, the explicit formulae are derived to determine the direction of bifurcations and the properties of bifurcating periodic solutions. Several numerical simulations and a simple discussion are given in the last Section.


## 2. The Working System

Braza [20] proposed the basic predator-prey system with logistic growth in the prey species, linear mortality in the predator species, and a square root functional response function. It is given by
(3)
d
X
d
t
=
r
X
(
1
-
X
K
)
-
α
X
Y
1
+
t
h
α
X
,
d
Y
d
t
=
-
s
Y
+
c
α
X
Y
1
+
t
h
α
X
,
where X
=
X
(
t
), Y
=
Y
(
t
) denote the number of the prey and predator species at time t, respectively. The parameter r is the growth rate of prey species, K is its carrying capacity, s is the death rate of predator species in the absence of prey, α is the search efficiency of predator for prey, c is the biomass conversion rate, and t
h is predator's average handling time of prey.Following References [23] and [24], we choose the quadratic mortality for predator species in system (3). If we further consider the delay fact, system (3) should be modified as the following form:
(4)
d
X
d
t
=
r
X
(
1
-
X
K
)
-
α
X
Y
1
+
t
h
α
X
,
d
Y
d
t
=
-
s
Y
2
+
c
α
X
(
t
-
τ
)
Y
1
+
t
h
α
X
(
t
-
τ
)
,
where -
s
Y
2 represents the quadratic mortality of predator species and τ is the time delay which means the growth rate of predator species to depend on the number of the prey species τ units of time earlier.By the same way of [20, 23], we make some scaling and assume that the average handling time is zero. Then the working system is
(5)
d
x
d
t
=
x
(
1
-
x
)
-
x
y
,
d
y
d
t
=
-
s
y
2
+
c
x
(
t
-
τ
)
y
.
The initial conditions for system (5) take the form
(6)
x
0
(
θ
)
=
ϕ
(
θ
)
≥
0
,
y
0
(
θ
)
=
ψ
(
θ
)
≥
0
,
θ
∈
[
-
τ
,
0
)
,
m
m
m
m
m
m
m
m
m
m
m
m
m
m
m
x
0
(
0
)
>
0
,
y
0
(
0
)
>
0
,
where (
ϕ
(
θ
)
,
ψ
(
θ
)
)
∈
C
(
[
-
τ
,
0
]
,
ℝ
¯
+
2
), ℝ
¯
+
2
=
{
(
r
1
,
r
2
)
:
r
1
≥
0
,
r
2
≥
0
}.It is easy to verify that system (5) always has a trivial equilibrium E
0
(
0,0
) and a boundary equilibrium E
1
(
1,0
) for any feasible parameters. If we assume c
<
s, then system (5) has a unique positive equilibrium E
*
(
x
*
,
y
*
), where
(7)
x
*
=
1
-
c
s
,
y
*
=
c
s
x
*
.


## 3. Stability Analysis and Hopf Bifurcation

Due to biological interpretation of the system, we only consider the positive equilibrium. In this section, by choosingτ as the bifurcation parameter and analyzing the corresponding linearized system, we investigate the stability of the positive equilibrium and the effects of the time delay on the dynamics of system (5).To study the local stability of the positive equilibriumE
*
(
x
*
,
y
*
), we first use the linear transformation x
1
=
x
-
x
*, x
2
=
y
-
y
*, for which system (5) can be written in the form
(8)
d
x
1
d
t
=
a
11
x
1
+
a
12
x
2
+
∑
i
+
j
≥
2
‍
f
i
j
(
1
)
i
!
j
!
x
1
i
x
2
j
,
d
x
2
d
t
=
a
21
x
1
(
t
-
τ
)
+
a
22
x
2
+
∑
i
+
j
≥
2
‍
f
i
j
(
2
)
i
!
j
!
x
1
i
(
t
-
τ
)
x
2
j
,
where
(9)
a
11
=
3
c
2
s
-
1
,
a
12
=
-
x
*
,
a
21
=
c
2
2
s
,
a
22
=
-
c
x
*
,
f
i
j
(
1
)
=
∂
i
+
j
f
(
1
)
(
x
1
,
x
2
)
∂
x
1
i
∂
x
2
j
|
(
0,0
)
,
f
i
j
(
2
)
=
∂
i
+
j
f
(
2
)
(
x
3
,
x
2
)
∂
x
3
i
∂
x
2
j
|
(
0,0
)
,
f
(
1
)
(
x
1
,
x
2
)
=
(
x
1
+
x
*
)
(
1
-
x
1
-
x
*
)
-
x
1
+
x
*
(
x
2
+
y
*
)
,
f
(
2
)
(
x
3
,
x
2
)
=
-
s
(
x
2
+
y
*
)
2
+
c
x
3
+
x
*
(
x
2
+
y
*
)
.
Thus, the positive equilibrium of system (5) is transformed into the zero equilibrium of system (8).The linear equations corresponding to (8) are given by
(10)
d
x
1
d
t
=
a
11
x
1
+
a
12
x
2
,
d
x
2
d
t
=
a
21
x
1
(
t
-
τ
)
+
a
22
x
2
,
and its characteristic equation is
(11)
λ
2
-
(
a
11
+
a
22
)
λ
+
a
11
a
22
-
a
12
a
21
e
-
λ
τ
=
0
.
To study the stability of the zero equilibrium and the Hopf bifurcation of system (8), we should analyze the distribution of the roots of characteristic equation (11). A simple calculation shows that a
11
a
22
-
a
12
a
21
=
c
(
x
*
)
3
/
2
>
0; thus λ
=
0 is not a root of (11).When there is no delay, that is,τ
=
0, characteristic equation (11) reduces to
(12)
λ
2
-
(
a
11
+
a
22
)
λ
+
a
11
a
22
-
a
12
a
21
=
0
.
If a
11
+
a
22
<
0, the two roots of (12) will have always negative real parts. Thus, we can obtain the following result.Lemma 1.
Assumec
<
s and a
11
+
a
22
<
0; then the two roots of characteristic equation (11) with τ
=
0 have always negative real parts; that is, the zero equilibrium of system (8) with τ
=
0 is locally asymptotically stable.Now forτ
>
0, let λ
=
i
ω (ω
>
0) be a root of (11); then ω satisfies the following equation:
(13)
-
ω
2
-
(
a
11
+
a
22
)
ω
i
+
a
11
a
22
-
a
12
a
21
(
cos
⁡
ω
τ
-
i
sin
ω
τ
)
=
0
.
Separating the real and imaginary parts, we have
(14)
-
ω
2
+
a
11
a
22
-
a
12
a
21
cos
⁡
ω
τ
=
0
,
-
(
a
11
+
a
22
)
ω
+
a
12
a
21
sin
ω
τ
=
0
,
which imply that
(15)
ω
4
+
(
a
11
2
+
a
22
2
)
ω
2
+
a
11
2
a
22
2
-
a
12
2
a
21
2
=
0
.
Notice a
11
a
22
-
a
12
a
21
>
0; we can find that (15) has no positive roots if a
11
a
22
+
a
12
a
21
≥
0 and has a positive root ω
+ if a
11
a
22
+
a
12
a
21
<
0, where
(16)
ω
+
=
-
(
a
11
2
+
a
22
2
)
+
(
a
11
2
-
a
22
2
)
2
+
4
a
12
2
a
21
2
2
.
Defining
(17)
τ
j
=
1
ω
+
(
arccos
-
(
ω
+
)
2
+
a
11
a
22
a
12
a
21
+
2
j
π
)
,
j
=
0,1
,
2
,
…
,
we have the following result.Lemma 2.
Assumec
<
s and a
11
+
a
22
<
0; then the following statements are true. (1)
Ifa
11
a
22
+
a
12
a
21
≥
0, the roots of (11) have always negative real parts for all τ
≥
0.
(2)
Ifa
11
a
22
+
a
12
a
21
<
0, (11) has a pair of purely imaginary roots ±
i
ω
+ for τ
=
τ
j.Next, under the condition of Lemma2(2), we will investigate whether the transversality condition is satisfied. Denote by λ
(
τ
)
=
α
(
τ
)
+
i
ω
(
τ
) the root of (11) near τ
=
τ
j satisfying α
(
τ
j
)
=
0, ω
(
τ
j
)
=
ω
+, j
=
0,1
,
2
,
…. Differentiating the two sides of (11) with respect to τ, we get
(18)
2
λ
d
λ
d
τ
-
(
a
11
+
a
22
)
d
λ
d
τ
+
a
12
a
21
e
-
λ
τ
(
λ
+
τ
d
λ
d
τ
)
=
0
;
that is
(19)
(
d
λ
d
τ
)
-
1
=
2
λ
-
a
11
-
a
22
-
a
12
a
21
λ
e
λ
τ
-
τ
λ
.
We can further get
(20)
sign
⁡
(
d
α
(
τ
)
d
τ
)
-
1
|
τ
=
τ
j
=
sign Re
(
d
λ
d
τ
)
-
1
|
τ
=
τ
j
=
sign Re
(
a
11
+
a
22
-
2
λ
a
12
a
21
λ
e
λ
τ
-
τ
λ
)
|
τ
=
τ
j
=
sign Re
(
a
11
+
a
22
-
2
i
ω
+
a
12
a
21
i
ω
+
(
cos
⁡
ω
+
τ
j
+
i
sin
ω
+
τ
j
)
)
=
sign Re
(
2
(
ω
+
2
-
a
11
a
22
)
a
12
2
a
21
2
+
(
a
11
+
a
22
)
2
a
12
2
a
21
2
)
=
sign Re
(
2
ω
+
2
+
a
11
2
+
a
22
2
a
12
2
a
21
2
)
>
0
.
Then we have the following conclusion.Lemma 3.
Assumec
<
s, a
11
+
a
22
<
0, and a
11
a
22
+
a
12
a
21
<
0; then the transversality conditions are satisfied; that is,
(21)
d
α
(
τ
)
d
τ
|
τ
=
τ
j
>
0
,
j
=
0,1
,
2
,
…
.Summarizing the above three lemmas, we can obtain the following theorem on stability and Hopf bifurcation of system (5).Theorem 4.
Assumec
<
s and a
11
+
a
22
<
0. For system (5), the following results are true.(1)
Ifa
11
a
22
+
a
12
a
21
≥
0, the positive equilibrium E
* is locally asymptotically stable for all τ
≥
0.
(2)
Ifa
11
a
22
+
a
12
a
21
<
0, then E
* is locally asymptotically stable when τ
∈
[
0
,
τ
0
) and unstable when τ
>
τ
0 and system (5) undergoes a Hopf bifurcation at E
* when τ
=
τ
j, j
=
0,1
,
2
,
….


## 4. Direction and Stability of the Hopf Bifurcation

In the previous section, we have obtained the conditions under which a family of periodic solutions bifurcate from the positive equilibriumE
* as τ crosses through the critical values τ
j. In this section, we will study the direction of Hopf bifurcation and the stability and period of the bifurcating periodic solutions by applying the normal form theory and the center manifold theorem from Hassard et al. [41]. Since the methods used are standard, we omit the detailed process and only give the main results. The readers can see [41] for more details on the derivation process.Assume that system (5) undergoes Hopf bifurcations at E
*
(
x
*
,
y
*
) for τ
=
τ
j, and then ±
i
ω
+ are the corresponding purely imaginary roots of characteristic equation (11). Let
(22)
q
(
θ
)
=
(
1
,
α
)
T
e
i
ω
+
τ
j
θ
,
q
*
(
ϑ
)
=
D
(
1
,
α
*
)
e
i
ω
+
τ
j
ϑ
m
m
m
m
m
m
m
m
m
m
m
m
m
m
m
m
for
θ
,
ϑ
∈
[
-
1,0
]
.
We can choose
(23)
α
=
i
ω
+
-
a
11
a
21
,
α
*
=
-
i
ω
+
-
a
11
a
21
e
i
ω
+
,
D
=
1
1
+
α
*
α
¯
+
τ
j
α
*
a
21
e
i
ω
+
τ
j
,
such that q
(
θ
) is the eigenvector of operator A corresponding to i
ω
+
τ
j, q
*
(
ϑ
) is the eigenvector of operator A
* corresponding to -
i
ω
+
τ
j, and
(24)
〈
q
*
(
ϑ
)
,
q
(
θ
)
〉
=
1
,
〈
q
*
(
ϑ
)
,
q
¯
(
θ
)
〉
=
0
.
Thus, we can compute
(25)
g
20
=
D
¯
τ
j
[
(
f
20
(
1
)
+
2
α
f
11
(
1
)
)
m
m
m
m
m
m
+
α
¯
*
(
e
-
2
i
ω
+
τ
j
f
20
(
2
)
+
2
α
e
-
i
ω
+
τ
j
f
11
(
2
)
+
α
2
f
02
(
2
)
)
]
,
g
11
=
D
¯
τ
j
[
(
f
20
(
1
)
+
(
α
+
α
¯
)
f
11
(
1
)
)
m
m
m
m
m
+
α
¯
*
(
f
20
(
2
)
+
(
α
e
i
ω
+
τ
j
+
α
¯
e
-
i
ω
+
τ
j
)
f
11
(
2
)
+
α
α
¯
f
02
(
2
)
)
]
,
g
02
=
D
¯
τ
j
[
(
f
20
(
1
)
+
2
α
¯
f
11
(
1
)
)
m
m
m
m
m
m
+
α
¯
*
(
e
2
i
ω
+
τ
j
f
20
(
2
)
+
2
α
¯
e
i
ω
+
τ
j
f
11
(
2
)
+
α
¯
2
f
02
(
2
)
)
]
,
g
21
=
D
¯
τ
j
(
w
1
+
α
¯
*
w
2
)
,
where
(26)
w
1
=
(
2
W
11
(
1
)
(
0
)
+
W
20
(
1
)
(
0
)
)
f
20
(
1
)
+
f
30
(
1
)
+
(
2
α
+
α
¯
)
f
21
(
1
)
w
1
+
(
2
W
11
(
2
)
(
0
)
+
W
20
(
2
)
(
0
)
+
α
¯
W
20
(
1
)
(
0
)
+
2
α
W
11
(
1
)
(
0
)
)
f
11
(
1
)
,
w
2
=
(
2
e
-
i
ω
+
τ
j
W
11
(
1
)
(
-
1
)
+
e
i
ω
+
τ
j
W
20
(
1
)
(
-
1
)
)
f
20
(
2
)
w
2
+
(
2
α
W
11
(
2
)
(
0
)
+
α
¯
W
20
(
2
)
(
0
)
)
f
02
(
2
)
w
2
+
(
2
α
W
11
(
1
)
(
-
1
)
+
α
¯
W
20
(
1
)
(
-
1
)
+
e
i
ω
+
τ
j
W
20
(
2
)
(
0
)
m
m
m
m
m
+
2
e
-
i
ω
+
τ
j
W
11
(
2
)
(
0
)
)
f
11
(
2
)
w
2
+
(
2
α
+
α
¯
e
-
2
i
ω
+
τ
j
)
f
21
(
2
)
+
f
30
(
2
)
e
-
i
ω
+
τ
j
,
W
20
(
θ
)
=
(
W
20
(
1
)
(
θ
)
W
20
(
2
)
(
θ
)
)
W
20
(
θ
)
=
i
g
20
ω
+
τ
j
q
(
0
)
e
i
ω
+
τ
j
θ
+
i
g
¯
20
3
ω
+
τ
j
q
¯
(
0
)
e
-
i
ω
+
τ
j
θ
+
E
1
e
2
i
ω
+
τ
j
θ
,
W
11
(
θ
)
=
(
W
11
(
1
)
(
θ
)
W
11
(
2
)
(
θ
)
)
W
11
(
θ
)
=
-
i
g
11
ω
+
τ
j
q
(
0
)
e
i
ω
+
τ
j
θ
+
i
g
¯
11
ω
+
τ
j
q
¯
(
0
)
e
-
i
ω
+
τ
j
θ
+
E
2
,
E
1
=
(
2
i
ω
+
-
a
11
-
a
12
-
a
21
e
-
2
i
ω
+
2
i
ω
+
-
a
22
)
-
1
E
1
×
(
f
20
(
1
)
+
2
α
f
11
(
1
)
e
-
2
i
ω
+
τ
j
f
20
(
2
)
+
2
α
e
-
i
ω
+
τ
j
f
11
(
2
)
+
α
2
f
02
(
2
)
)
,
E
2
=
-
(
a
11
a
12
a
21
a
22
)
-
1
E
2
×
(
f
20
(
1
)
+
(
α
+
α
¯
)
f
11
(
1
)
f
20
(
2
)
+
(
α
e
i
ω
+
τ
j
+
α
¯
e
-
i
ω
+
τ
j
)
f
11
(
2
)
+
α
α
¯
f
02
(
2
)
)
.
Now, let
(27)
c
1
(
0
)
=
i
2
ω
+
τ
j
(
g
20
g
11
-
2
|
g
11
|
2
-
|
g
02
|
2
3
)
+
g
21
2
,
which can be determined by the parameters in system (5). Defining
(28)
μ
2
=
-
Re
{
c
1
(
0
)
}
Re
{
λ
′
(
τ
j
)
}
,
β
2
=
2
Re
{
c
1
(
0
)
}
,
T
2
=
-
Im
⁡
{
c
1
(
0
)
}
+
μ
2
Im
⁡
{
λ
′
(
τ
j
)
}
ω
+
τ
j
,
we have the following result.Theorem 5.
The direction of the Hopf bifurcation is determined byμ
2: if μ
2
>
0
(
μ
2
<
0
), then the Hopf bifurcation is supercritical (subcritical) and the bifurcating periodic solutions exist for τ
>
τ
0
(
τ
<
τ
0
); β
2 determines the stability of the bifurcating periodic solutions: the bifurcating periodic solutions are stable (unstable) if β
2
<
0
(
β
2
>
0
); and T
2 determines the period of the bifurcating periodic solutions: the period increases (decreases) if T
2
>
0
(
T
2
<
0
).


## 5. Numerical Simulations and Discussions

In this section, we present some numerical simulations to verify our theoretical results proved in previous sections by using MATLAB DDE solver. We simulate the system (5) with the initial value (
x
0
,
y
0
)
=
(
0.4,0.39
). In Figures 1–3 the initial point (
x
0
,
y
0
) and the positive equilibrium E
* are represented by green star and red star, respectively.The phase portraits of system (5). Initial values x
0
=
0.4, y
0
=
0.39. Parameters s
=
1, c
=
0.4. Delay values (a) τ
=
0, (b) τ
=
5, (c) τ
=
10, and (d) τ
=
15.
(a)
(b)
(c)
(d)Figure 2
The trajectories and phase portrait of system (5). Initial values x
0
=
0.4, y
0
=
0.39. Parameters s
=
1, c
=
0.6. Delay value τ
=
5.248.Figure 3
The trajectories and phase portrait of system (5). Initial values x
0
=
0.4, y
0
=
0.39. Parameters s
=
1, c
=
0.6. Delay value τ
=
6.049.First, we choose parameterss
=
1, c
=
0.4. In this case, system (5) has only one positive equilibrium E
*
=
(
0.6,0.3098
). Simple calculations show that a
11
+
a
22
=
-
0.7098
<
0, a
11
a
22
+
a
12
a
21
=
0.062
>
0, which satisfy the conditions of Theorem 4(1). By Theorem 4(1), we expect that E
* is stable for any delay value. We further choose four different delay values: τ
a
=
0, τ
b
=
5, τ
c
=
10, and τ
d
=
15. The computation simulations are depicted in Figure 1. Obviously, the computation simulations can support the result of Theorem 4(1).Next, we choose parameterss
=
1, c
=
0.6. The positive equilibrium of system (5) is E
*
=
(
0.4,0.3795
). It is easy to show that a
11
+
a
22
=
-
0.4795
<
0, a
11
a
22
+
a
12
a
21
=
-
0.0759
<
0, which satisfy the conditions of Theorem 4(2). From the formulae in the previous section, we can compute the values of τ
0, μ
2, β
2, and T
2 as
(29)
τ
0
=
6.048
,
μ
2
=
130888
,
β
2
=
-
503.3349
,
T
2
=
681.8543
,
from which we conclude that E
* is asymptotically stable for τ
<
τ
0 and the Hopf bifurcation of system (5) occurring at the critical value τ
0 is supercritical and the bifurcating periodic solution exists when τ crosses τ
0 to the right; also the bifurcating periodic solution is stable. By taking τ
=
5.248
<
τ
0 and τ
=
6.049
>
τ
0, we can show the computation simulations in Figures 2 and 3, respectively. It is easy to see that the computation simulations support the results of Theorem 4(2) and Theorem 5.In the present paper, we have considered a delayed predator-prey system in which the prey species exhibits herd behavior and the predator species with quadratic mortality. Our research shows that, for system (5), the positive equilibrium E
* is always asymptotically stable under certain conditions (see Theorem 4(1)). In other words, under such conditions, the dynamics behavior of the corresponding ODE system is robust with respect to the time delay. But, under some other conditions (see Theorem 4(2)), E
* will lose its stability and Hopf bifurcations occur when the delay τ passes through some critical values. In this case, the time delay has an important effect for the system dynamics.Our results may enrich the dynamics in the predator-prey system and help us to better understand the interaction of predator with prey in a real ecosystem. Further studies are necessary to analyze the dynamics of more realistic but complex systems, such as delayed diffusive predator-prey system with herd behavior.



---

*Source: 568943-2014-04-08.xml*