# Periodic Solutions and Homoclinic Bifurcations of Two Predator-Prey Systems with Nonmonotonic Functional Response and Impulsive Harvesting

Authors: Mingzhan Huang; Xinyu Song
Journal: Journal of Applied Mathematics
 (2014)

---

## Abstract

Two predator-prey models with nonmonotonic functional response and state-dependent impulsive harvesting are formulated and analyzed. By using the geometry theory of semicontinuous dynamic system, we obtain the existence, uniqueness, and stability of the periodic solution and analyse the dynamic phenomenon of homoclinic bifurcation of the first system by choosing the harvesting rateβ as control parameter. Besides, we also study the homoclinic bifurcation of the second system about parameter δ on the basis of the theory of rotated vector field. Finally, numerical simulations are presented to illustrate the results.


---

## Body

## 1. Introduction

Predator-prey interaction is one of the most important relationships in the ecosystem, so it has long been a focus of study in mathematical ecology. The functional response of the predator to the prey describes how the predator density changes with the prey density. The well-known response functions, which have appeared in a lot of literatures, include linear function, Holling type-II function, Holling type-III function, and so on. The above mentioned response functions are all monotonic, and they are really accurate to describe the predator-prey interactions in many cases. But in the recent decades, scholars in different fields found through experiments that nonmonotonic functional response occurs in some predator-prey interactions. For example, if the prey exhibits group defense [1, 2], the predator growth rate will be inhibited when the prey density reaches a high level. This phenomenon is known to exist widely in nature and considerable work on it has been studied [3–8]. To study the predator-prey interaction when the prey exhibits group defense, Ruan and Xiao in [4] proposed the following model:
(1)
d
x
d
t
=
r
x
(
1
-
x
K
)
-
x
y
a
+
x
2
,
d
y
d
t
=
y
(
μ
x
a
+
x
2
-
D
)
,
where x and y represent the population densities of prey and predator, respectively; K
>
0 and r
>
0 represent the carrying capacity and the intrinsic birth rate of the prey, respectively, and D
>
0 is the death rate of the predator. The function p
(
x
)
=
x
/
(
a
+
x
2
) denotes the predator functional response which can model the phenomenon of group defense because there exists h
¯
>
0 such that p
′
(
x
)
>
0 for 0
≤
x
<
h
¯ and p
′
(
x
)
<
0 for x
>
h
¯.Harvesting strategy of biological resources is also a focus topic in mathematical bioeconomics because it relates to the optimal management of renewable resources. In many literatures, impulsive differential equations are used to model the human action of harvesting. Predator-prey systems with periodic impulsive harvesting have been studied extensively and important results have been achieved [9–12]. In recent years, more and more scholars begin to consider state dependent impulsive harvesting instead of the periodic impulsive form [13–17]. In point of ecology, we would better know some information about the amount of the species when we harvest them, thus we can avoid excessive exploitation and resource exhaustion. To this end, we introduce a reliable real-time monitoring system to estimate the number of the species. Such monitoring systems exist in many cases. According to the feedback information from the monitoring system, we can manage our resources better. For predator-prey system (1), we assume the predator has high commercial value for human beings, and people increase its production mainly through replenishing its prey. In such an ecosystem, we suppose the amount of the prey can be estimated by a monitoring system, and the monitoring data can help us to decide if we harvest the predator or not. To model this harvesting behavior, we can propose the following state-dependent impulsive differential equations:
(2)
d
x
d
t
=
r
x
(
1
-
x
K
)
-
x
y
a
+
x
2
,
d
y
d
t
=
y
(
μ
x
a
+
x
2
-
D
)
,
d
x
d
t
=
r
t
t
t
t
t
t
t
t
x
x
>
h
,
Δ
x
=
τ
,
Δ
y
=
-
β
y
,
Δ
Δ
y
x
=
h
,
where h
>
0 is a threshold. When the amount of the prey is large than h, which means the food of the predator is abundant, and the development of the system coincide with our economic interest. When the amount of the prey drops to the threshold h, which means the nutrition of the predator will be deficient, we harvest the predator at rate β
∈
(
0,1
) and replenish some prey at the same time. We denote the recruitment of the prey as τ.A lot of research has shown that the population size of the predator is influenced not only by the prey populations but also by the relative rate of prey population growth. If we take the relative growth rate effect into consideration, the above model needs some changes. To this end, we propose the following system:(3)
d
x
d
t
=
r
x
(
1
-
x
K
)
-
x
y
a
+
x
2
,
d
y
d
t
=
y
(
μ
x
a
+
x
2
-
D

)

+

δ
[
r
(

1
-
x
K
)
-
y
a
+
x
2
]
2
,
00000000000000000000000000000000000
x
>
h
,
Δ
x
=
τ
,
Δ
y
=
-
β
y
,
00,0
x
=
h
,
where the term δ
[
r
(
1
-
(
x
/
K
)
)
-
y
/
(
a
+
x
2
)
]
2 represents the effect of the prey relative growth rate on the predator.In this paper, we mainly discuss the dynamics properties of the systems (2) and (3). The paper is organized as follows. In Section 2, some notations and definitions of the geometric theory of semicontinuous dynamical systems are provided. In Section 3, by using the geometry theory of semicontinuous dynamic system, we firstly study the existence, uniqueness, and orbital stability of periodic solutions and analyse the dynamic phenomenon of homoclinic bifurcation of system (2) by choosing the harvesting rate β as control parameter. Then, we study the homoclinic bifurcation of system (3) about parameter δ on the basis of the theory of rotated vector field. The paper ends with a brief discussion and some numerical simulations.


## 2. Preliminaries

In this section, we give some notations and definitions of the geometric theory of semicontinuous dynamical systems which will be useful for the following discussion.Definition 1 (see [9]).
Consider the state-dependent impulsive differential equations(4)
d
x
d
t
=
P
¯
(
x
,
y
)
,
d
y
d
t
=
Q
¯
(
x
,
y
)
,
(
x
,
y
)
∉
M
{
x
,
y
}
,
Δ
x
=
α
(
x
,
y
)
,
Δ
y
=
β
(
x
,
y
)
,
(
x
,
y
)
∈
M
{
x
,
y
}
.
We define the dynamic system consisting of the solution mapping of the system (4) a semicontinuous dynamical system, denoted as (
Ω
,
f
,
φ
,
M
). We require that the initial point P of the system (4) should not be in the set M
{
x
,
y
}, that is, P
∈
Ω
=
R
+
2
∖
M
{
x
,
y
}, and the function φ is a continuous mapping that satisfies φ
(
M
)
=
N. Here φ is called the impulse mapping, where M
{
x
,
y
} and N
{
x
,
y
} represent the straight lines or curves in the plane R
+
2, M
{
x
,
y
} is called the impulse set, and N
{
x
,
y
} is called the phase set.Remark 2.
For the systems (2) and (3), M
=
{
(
x
,
y
)
∣
x
=
h
,
y
≥
0
},  N
=
{
(
x
,
y
)
∣
x
=
h
+
τ
,
y
≥
0
}, and for any (
x
,
y
)
∈
M, we have φ
(
x
,
y
)
=
(
h
+
τ
,
(
1
-
β
)
y
).Definition 3 (see [9]).
For the semicontinuous dynamical system defined by the state-dependent impulsive differential equations (4), the solution mapping f
(
P
,
t
)
:
Ω
→
Ω consists of two parts.(1)
Letπ
(
P
,
t
) denote the poincaré mapping with the initial point P of the following system:
(5)
d
x
d
t
=
P
¯
(
x
,
y
)
,
d
y
d
t
=
Q
¯
(
x
,
y
)
.
If f
(
P
,
t
)
∩
M
{
x
,
y
}
=
∅, then f
(
P
,
t
)
=
π
(
P
,
t
).
(2)
If there exists a time pointT
1 such that f
(
P
,
T
1
)
=
H
∈
M
{
x
,
y
}, φ
(
H
)
=
φ
(
f
(
P
,
T
1
)
)
=
P
1
∈
N
{
x
,
y
}, and f
(
P
1
,
t
)
∩
M
{
x
,
y
}
=
∅, then f
(
P
,
t
)
=
π
(
P
,
T
1
)
+
f
(
P
1
,
t
) (see Figure 1(a)).(a) The solution mapping of the system (4). (b) The order one periodic solution of the system (4). (c) The order two periodic solution of the system (4).
(a)
(b)
(c)Remark 4.
For (2) in Definition3, if f
(
P
1
,
t
)
∩
M
{
x
,
y
}
≠
∅, and there exists a time point T
2 such that f
(
P
1
,
T
2
)
=
H
1
∈
M
{
x
,
y
}, φ
(
H
1
)
=
φ
(
f
(
P
1
,
T
2
)
)
=
P
2
∈
N
{
x
,
y
}, and f
(
P
2
,
t
)
∩
M
{
x
,
y
}
=
∅, then f
(
P
,
t
)
=
π
(
P
,
T
1
)
+
f
(
P
1
,
t
)
=
π
(
P
,
T
1
)
+
π
(
P
1
,
T
2
)
+
f
(
P
2
,
t
).
Iff
(
P
2
,
t
)
∩
M
{
x
,
y
}
≠
∅
,
…
,
f
(
P
k
-
1
,
t
)
∩
M
{
x
,
y
}
≠
∅ and f
(
P
k
,
t
)
∩
M
{
x
,
y
}
=
∅, then we can repeat the above steps and have the following form:
(6)
f
(
P
,
t
)
=
∑
i
=
1
k
‍
π
(
P
i
-
1
,
T
i
)
+
f
(
P
k
,
t
)
,
P
0
=
P
.Definition 5 (see [9]).
If there exists a pointP
∈
N
{
x
,
y
} and a time point T
1 such that f
(
P
,
T
1
)
=
H
∈
M
{
x
,
y
} and φ
(
H
)
=
φ
(
f
(
P
,
T
1
)
)
=
P
∈
N
{
x
,
y
}, then f
(
P
,
t
) is called an order one periodic solution of the system (4) whose period is T
1 (see Figure 1(b)). The orbit of the order one periodic solution is called an order one cycle. If there exists a singularity in the order one cycle, we call it an order one singular cycle. If the singularity is a saddle, we call it an order one homoclinic cycle.Definition 6 (see [9]).
If there exists a pointP
∈
N
{
x
,
y
} and a time point T
1 such that f
(
P
,
T
1
)
=
H
∈
M
{
x
,
y
}, φ
(
H
)
=
P
1
∈
N
{
x
,
y
} and P
≠
P
1, and there also exists a time point T
2 such that f
(
P
1
,
T
2
)
=
H
1
∈
M
{
x
,
y
} and φ
(
H
1
)
=
P
∈
N
{
x
,
y
}, then f
(
P
,
t
) is called an order two periodic solution of the system (4) whose period is T
1
+
T
2 (see Figure 1(c)). Analogously, we can define the order k periodic solution of the system (4).Definition 7.
SupposeΓ
=
f
(
P
,
t
) is an order one periodic solution of the system (4). If for any ε
>
0, there must exist δ
>
0 and t
0
≥
0, such that for any point P
1
∈
U
(
P
,
δ
)
∩
N
{
x
,
y
}, we have ρ
(
f
(
P
1
,
t
)
,
Γ
)
<
ε for t
>
t
0; then we call the order one periodic solution Γ orbitally asymptotically stable.Definition 8 (see [9]).
Suppose the impulse setM and phase set N of the system (4) are straight lines and a coordinate system can be defined in the phase set N. Let point A
∈
N and its coordinate is a. Assume that the trajectory from the point A intersects the impulse set M at a point A
′, and, after impulsive effect, the point A
′ is mapped to the point A
1
∈
N with the coordinate a
1; then we call point A
1 the order one successor point of point A, and the order one successor function of point A is F
1
(
A
)
=
a
1
-
a.Remark 9.
For Definition8, if the trajectory from the point A
1 intersects the impulse set M again at a point A
1
′, and, after impulsive effect, the point A
1
′ is mapped to the point A
2
∈
N with the coordinate a
2, then the point A
2 is obviously the order one successor point of point A
1; we also call point A
2 the order two successor point of point A, and the order two successor function of point A is F
2
(
A
)
=
a
2
-
a. If the process can be repeated over and over again, then we can define the order k successor point of point A (which we denote as A
k and its coordinate is a
k) and the order k successor function of A which we denote as F
k
(
A
)
=
a
k
-
a.Remark 10.
For the systems (2) and (3), for any point H
∈
N, we define the directed distance between point H and the x-axis as the coordinate of point H. In this paper, we denote the coordinate of point H as y
H.Lemma 11 (see [9]).
Successor functionF
k
(
A
) is continuous.Lemma 12.
For the systems (2) and (3), if there exists two points A
∈
N, B
∈
N such that F
1
(
A
)
F
1
(
B
)
<
0, then there must exist a point C
∈
N which is between the points A and B such that F
1
(
C
)
=
0; thus, the system must have an order one periodic solution which passes through the point C.Proof.
By Lemma11, we can easily get that there must exist a point C
∈
N which is between the points A and B such that F
1
(
C
)
=
0. According to Definition 6, we know Γ
=
f
(
C
,
t
) is an order one periodic solution. That completes the proof.Lemma 13.
For the systems (2) and (3), suppose point A
∈
N and F
k
(
A
)
≠
0, then the system does not have an order k periodic solution which passes through the point A.


## 3. Periodic Solution and Homoclinic Bifurcation

In this section, we firstly discuss the existence, uniqueness, and stability of the periodic solution of the system (2) by using differential equation geometry theory then study the dynamic phenomena of homoclinic bifurcation of both systems (2) and (3). We now list some results about the system (1) that are given by Ruan and Xiao in [4].For the system (1), there is always a hyperbolic saddle point at the origin (
0,0
) and an equilibrium (
K
,
0
) in the x-axis. When μ
2
-
4
a
D
2
>
0, the system (1) may have one or two positive equilibria. If the positive equilibria exist, we denote them as (
x
1
,
y
1
), (
x
2
,
y
2
), where
(7)
x
1
=
μ
-
μ
2
-
4
a
D
2
2
D
,
y
1
=
r
(
1
-
x
1
K
)
(
a
+
x
1
2
)
,
x
2
=
μ
+
μ
2
-
4
a
D
2
2
D
,
y
2
=
r
(
1
-
x
2
K
)
(
a
+
x
2
2
)
.
Now, we give some results appeared in [4].Lemma 14 (see [4]).
If4
a
D
2
<
μ
2
≤
(
16
/
3
)
a
D
2 and x
1
<
K
<
x
2 (or if μ
2
>
(
16
/
3
)
a
D
2 and x
1
<
K
<
-
x
1
+
2
μ
x
1
/
D), then the system (1) has three equilibria: two hyperbolic saddles (
0,0
) and (
K
,
0
) and a globally asymptotically stable equilibrium (
x
1
,
y
1
) in the interior of the first quadrant.Lemma 15 (see [4]).
If4
a
D
2
<
μ
2
<
(
16
/
3
)
a
D
2 and x
2
<
K
<
(
2
μ
-
μ
2
-
4
a
D
2
)
/
2
D (or if μ
2
>
4
a
D
2 and K
>
(
2
μ
+
μ
2
-
4
a
D
2
)
/
2
D), then the system (1) has four equilibria: two hyperbolic saddles (
0,0
) and (
x
2
,
y
2
), a hyperbolic stable node (
K
,
0
), and a stable (an unstable) equilibrium (
x
1
,
y
1
), and system (1) has no closed orbits.


### 3.1. Periodic Solution and Homoclinic Bifurcation of System (2) about Parameter β

According to the model (2), the threshold h and the recruitment τ of the prey should satisfy the condition 0
<
h
<
h
+
τ
<
K by ecological significance. For this consideration, we have the following results.Theorem 16.
If4
a
D
2
<
μ
2
≤
(
16
/
3
)
a
D
2, x
1
<
K
<
x
2 (or if μ
2
>
(
16
/
3
)
a
D
2, x
1
<
K
<
-
x
1
+
2
μ
x
1
/
D), and x
1
<
h
<
h
+
τ
<
K, then there must exist a fixed value β
0
∈
(
0,1
) such that for every β
∈
(
β
0
,
1
), the system (2) has a unique order one periodic solution in region Ω
1, where region Ω
1 is the region enclosed by the x-axis, the impulse set x
=
h, and the unstable flow of the saddle (
K
,
0
).Proof.
According to Lemma14, the system (1) has three equilibria: two hyperbolic saddles (
0,0
) and (
K
,
0
) and a globally asymptotically stable equilibrium (
x
1
,
y
1
) in the interior of the first quadrant. For convenience, we suppose the x-axis intersects the impulse set x
=
h and the phase set x
=
h
+
τ at point A
′ and point B, respectively. The unstable flow of (
K
,
0
) also passes through the impulse set and the phase set, we denote the intersections by point A and point O, respectively. Besides, we suppose the vertical isocline d
x
/
d
t
=
0 intersects the impulse set and the phase set at point C and point D, respectively. Then the region Ω
1 is the interior of the closed curve E
0
O
A
^
∪
A
C
A
′
¯
∪
A
′
B
E
0
¯, where E
0 denotes the equilibrium (
K
,
0
) (see Figure 2).
According to the impulsive conditions of the system (2), there must exist a fixed value β
0
∈
(
0,1
), when β
=
β
0, point A is mapped to the point D after impulsive effect, that is to say, (
1
-
β
0
)
y
A
=
y
D.
Whenβ
0
<
β
<
1, after impulsive effect, point A is mapped to a point between points B and D. Select a point B
0 next to B in the phase set x
=
h
+
τ, the trajectory of the system (2) from point B
0 must intersect the impulse set x
=
h at a point B
0
′ which is next to point A, and after impulsive effect, point B
0
′ is mapped to a point B
1 which is the order one successor point of B
0; then, we have y
B
1
=
(
1
-
β
)
y
B
0
′
<
(
1
-
β
)
y
A
<
y
D. That is to say point B
1 is between points D and B
0. The trajectory from point B
1 must intersect the impulse set x
=
h again at a point B
1
′, and after impulsive effect, point B
1
′ is mapped to a point B
2. Since distinct trajectories do not intersect, we can easily have y
C
<
y
B
1
′
<
y
B
0
′ and y
B
2
=
(
1
-
β
)
y
B
1
′
<
(
1
-
β
)
y
B
0
′
=
y
B
1. Obviously, point B
2 is the order one successor point of point B
1, and then we have the following results of the order one successor function:
(8)
F
1
(
B
0
)
=
y
B
1
-
y
B
0
>
0
,
F
1
(
B
1
)
=
y
B
2
-
y
B
1
<
0
.
By Lemma 12, we know that there must exist a point M in the phase set x
=
h
+
τ which is between the points B
0 and B
1 such that F
1
(
M
)
=
0. Then we know the system (2) has an order one periodic solution which passes through the point M.
In the following, we prove the uniqueness of the order one periodic solution. Arbitrarily choose two pointsA
1 and A
2 in the phase set x
=
h
+
τ, where 0
≤
y
A
1
<
y
A
2
≤
y
D. The trajectories of the system (2) through points A
1 and A
2 must intersect the impulse set x
=
h at some points A
1
′ and A
2
′, respectively, and after impulsive effect, the points A
1
′ and A
2
′ must be mapped to the phase set x
=
h
+
τ at some points A
1
′′ and A
2
′′, respectively (see Figure 3). Because distinct trajectories do not intersect, we can easily have y
C
<
y
A
2
′
<
y
A
1
′
<
y
A and y
A
2
′′
=
(
1
-
β
)
y
A
2
′
<
y
A
1
′′
=
(
1
-
β
)
y
A
1
′; then we get the order one successor functions must satisfy
(9)
F
1
(
A
2
)
-
F
1
(
A
1
)
=
(
y
A
2
′′
-
y
A
2
)
-
(
y
A
1
′′
-
y
A
1
)
=
(
y
A
2
′′
-
y
A
1
′′
)
+
(
y
A
1
-
y
A
2
)
<
0
,
which means the order one successor function F
1 is monotone decreasing in the line segment B
D
¯, thus there exists only one point M
∈
B
D
¯ such that F
1
(
M
)
=
0.
Besides, for any pointH
∈
D
O
¯, it is easy to know that F
1
(
H
)
<
0. So the system (2) has only one order one periodic solution in the region Ω
1. That completes the proof.Figure 2
Existence of order one periodic solution of (2).Figure 3
Uniqueness of order one periodic solution of (2).Theorem 17.
Under the conditions of Theorem16, the order one periodic solution of the system (2) is orbitally asymptotically stable.Proof.
According to Theorem16, the system (2) has a unique order one periodic solution that passes through the point M which is in the phase set x
=
h
+
τ, where y
B
0
<
y
M
<
y
B
1. The trajectory of the system (2) through point B
1 intersects the impulse set x
=
h at point B
1
′, and after impulsive effect, the point B
1
′ is mapped to point B
2, where y
B
<
y
B
2
<
y
M. Besides, the trajectory of the system (2) from point B
2 must intersect the impulse set x
=
h again at a point B
2
′, and after impulsive effect, the point B
2
′ is mapped to a point B
3 in the phase set x
=
h
+
τ, where y
B
1
>
y
B
3
>
y
M (see Figure 4).
Repeat the above steps, the trajectory from pointB
0 will come across impulsive effect infinite times. Denote the phase point corresponding to the ith impulsive effect by B
i, i
=
1,2
,
…. Obviously, B
i is also the order i successor point of point B
0. Then we have
(10)
y
B
0
<
y
B
2
<
y
B
4
<
⋯
<
y
B
2
k
<
y
B
2
(
k
+
1
)
<
⋯
<
y
M
,
y
B
1
>
y
B
3
>
y
B
5
>
⋯
>
y
B
2
k
+
1
>
y
B
2
(
k
+
1
)
+
1
>
⋯
>
y
M
.
Thus {
y
B
2
k
}, k
=
0,1
,
2
,
…is a monotonically increasing sequence, and {
y
B
2
k
+
1
}, k
=
0,1
,
2
,
…is a monotonically decreasing sequence. Because the order one successor function F
1 is monotone decreasing in the line segment B
D
¯, we have
(11)
F
1
(
B
0
)
=
y
B
1
-
y
B
0
>
F
1
(
B
2
)
=
y
B
3
-
y
B
2
>
F
1
(
B
4
)
=
y
B
5
-
y
B
4
>
⋯
>
0
,
and furthermore,
(12)
y
B
2
k
⟶
y
M
,
as
k
⟶
∞
;
y
B
2
k
+
1
⟶
I
M
,
as
k
⟶
∞
.
Pick any point Q
0
∈
B
0
B
1
¯ different from the point M; without loss of generality, we assume y
B
0
<
y
Q
0
<
y
M (otherwise, y
M
<
y
Q
0
<
y
B
1; the discussions are similar), there must exist a natural number n
0 such that y
B
2
n
0
<
y
Q
0
<
y
B
2
(
n
0
+
1
). The trajectory from point Q
0 will also undergo impulsive effect infinite times, and we denote the phase point corresponding to the kth impulsive effect as Q
k, k
=
1,2
,
…. For any natural number l, we have y
B
2
(
n
0
+
l
)
<
y
Q
2
l
<
y
B
2
(
n
0
+
l
+
1
) and y
B
2
(
n
0
+
l
+
1
)
+
1
<
y
Q
2
l
+
1
<
y
B
2
(
n
0
+
l
)
+
1. Then {
y
Q
2
l
}, l
=
0,1
,
2
,
…  is also a monotonically increasing sequence, {
y
Q
2
l
+
1
}, l
=
0,1
,
2
,
… is also a monotonically decreasing sequence, and we get
(13)
y
Q
2
l
⟶
y
M
,
as
l
⟶
∞
;
y
Q
2
l
+
1
⟶
I
M
,
as
l
⟶
∞
.
This indicates that the trajectory of the system (2) from point Q
0
∈
B
0
B
1
¯ ultimately tends to the trajectory passing through point M, which means the order one periodic solution is orbitally asymptotically stable. That completes the proof.Figure 4
Illustration of the orbital asymptotic stability of the order one periodic solution of the system (2).Theorem 18.
Under the conditions of Theorem16, the system (2) has no order k periodic solution in region Ω
1, where k
≥
2.Proof.
For any pointS
∈
B
1
O
¯, y
B
1
<
y
S
<
y
O, the trajectory from point S will undergo impulsive effect infinite times. Denote the impulsive point and phase point corresponding to the k
th
(
k
=
1,2
,
…
) impulsive effect by S
k
′ and S
k, respectively. Similar to the discussions in Theorem 17, we can easily know y
C
<
y
S
k
′
<
y
B
0
′, k
=
1,2
,
…, then we have y
S
k
=
(
1
-
β
)
y
S
k
′
<
(
1
-
β
)
y
B
0
′
=
y
B
1
<
y
S. Obviously, point S
k is the order k successor point of point S, then the order k successor function F
k
(
S
)
=
y
S
k
-
y
S
≠
0. By Lemma 13, there does not exist an order k
(
k
=
1,2
,
…
) periodic solution through point S.
For any pointS
∈
B
M
¯
,
y
B
<
y
S
<
y
M, according to the proof of Theorem 17, there must exist a natural number n
0 such that y
B
2
n
0
<
y
S
<
y
B
2
(
n
0
+
1
). Denote the order k
(
k
=
1,2
,
…
) successor point of point S by point S
k, then we have y
B
2
(
n
0
+
l
)
<
y
S
2
l
<
y
B
2
(
n
0
+
l
+
1
)
<
y
M and y
M
<
y
B
2
(
n
0
+
l
+
1
)
+
1
<
y
S
2
l
+
1
<
y
B
2
(
n
0
+
l
)
+
1. Then {
y
S
2
l
}, l
=
0,1
,
2
,
…is amonotonically increasing sequence where S
0
=
S, {
y
S
2
l
+
1
}, l
=
0,1
,
2
,
…is a monotonically decreasing sequence, and
(14)
y
S
2
l
⟶
y
M
,
as
l
⟶
∞
;
y
S
2
l
+
1
⟶
I
M
,
as
l
⟶
∞
.
When k
=
2
l, we have y
M
>
y
S
k
>
y
S
0
=
y
S and the order k successor function F
k
(
S
)
=
y
S
k
-
y
S
>
0. When k
=
2
l
+
1, we have y
S
k
>
y
M
>
y
S
0
=
y
S and F
k
(
S
)
=
y
S
k
-
y
S
>
0. That is to say there does not exist an order k
(
k
=
1,2
,
…
) periodic solution passing through point S.
Analogously, for any pointS
∈
M
B
1
¯
,
y
M
<
y
S
<
y
B
1, we can prove there does not exist an order k
(
k
=
1,2
,
…
) periodic solution passing through point S.
From the above discussion, we know that the order one periodic solution is the unique periodic solution of the system (2) and there is no order k
(
k
≥
2
) periodic solution in region Ω
1. That completes the proof.Theorem 19.
If4
a
D
2
<
μ
2
<
(
16
/
3
)
a
D
2, x
2
<
K
<
(
2
μ
-
μ
2
-
4
a
D
2
)
/
2
D (or if μ
2
>
4
a
D
2, K
>
(
2
μ
+
μ
2
-
4
a
D
2
)
/
2
D), x
1
<
h
<
h
+
τ
<
x
2 and K
≤
3
a, then there must exist two fixed values β
0 and β
* which satisfy 0
<
β
0
<
β
*
<
1 such that(1)
whenβ
=
β
*, the system (2) has an order one homoclinic cycle;
(2)
ifβ
0
<
β
<
β
*, then the homoclinic cycle of system (2) disappears and bifurcates an order one periodic solution in region Ω
2, where region Ω
2 is the region enclosed by the x-axis, the impulse set x
=
h, and the unstable flow of the saddle (
x
2
,
y
2
). Furthermore, the order one periodic solution is unique and is orbitally asymptotically stable;
(3)
ifβ
*
<
β
<
1, then the homoclinic cycle of system (2) also disappears and there is no periodic solution in region Ω
2.Proof.
ForK
≤
3
a, it is easy to know that y
1
>
y
2. According to Lemma 15, the system (1) has four equilibria: two hyperbolic saddles (
0,0
) and (
x
2
,
y
2
), a hyperbolic stable node (
K
,
0
), and a node or focus (
x
1
,
y
1
). Obviously, the x-axis, the vertical isocline d
x
/
d
t
=
0, and the unstable and stable flow of (
x
2
,
y
2
) are all intersected with the impulse set x
=
h and phase set x
=
h
+
τ. For convenience, we denote the intersections by points A
′ and B
′, C and D, A and O, and A
1 and B, respectively (see Figure 5). Then the region Ω
2 is the interior of the closed curve E
0
E
2
O
A
^
∪
A
C
A
′
¯
∪
A
′
B
′
E
0
¯, where E
0 and E
2 denote the equilibria (
K
,
0
) and (
x
2
,
y
2
), respectively. According to the impulsive conditions of the system (2), there must exist a fixed value β
0
∈
(
0,1
), when β
=
β
0, point A is mapped to the point D after impulsive effect. Also, there must exist a fixed value β
*
∈
(
β
0
,
1
), when β
=
β
*, point A is mapped to the point B after impulsive effect. Thus we have (
1
-
β
0
)
y
A
=
y
D and (
1
-
β
*
)
y
A
=
y
B.
Whenβ
=
β
*, we have (
1
-
β
)
y
A
=
y
B; we know the curve B
E
2
O
A
^
∪
A
B
¯ is an order one circle which has the saddle E
2
(
x
2
,
y
2
) in it; that is to say, the system (2) has an order one homoclinic cycle.
Ifβ
∈
(
β
0
,
β
*
), after impulsive effect, point A is mapped to a point between points B and D. Similar to the discussion of Theorems 16 and 17, we can prove that the system (2) has an order one periodic solution which passes through a point M
∈
B
D
¯ and is orbitally asymptotically stable. Besides, for any point H
∈
B
′
B
¯, the trajectory of the system (2) from point H must ultimately tend to the equilibrium (
K
,
0
) and does not undergo any impulsive effect, that is to say the system (2) has no periodic solution passing through the point H; so we can know the system (2) has a unique order one periodic solution in the region Ω
2 which is orbitally asymptotically stable.
Whenβ
∈
(
β
*
,
1
), we have 0
<
(
1
-
β
)
y
A
<
y
B. Then the trajectory of the system (2) through point H
∈
O
B
¯ must tend to the equilibrium (
K
,
0
) after coming across once impulsive effect and the trajectory passing through point H
∈
B
′
B
¯ must tend to the equilibrium (
K
,
0
) without undergoing any impulsive effect. So we get that the system (2) has no periodic solution in region Ω
2. That completes the proof.Figure 5
Illustration of homoclinic bifurcation of system (2) by choosing the harvesting rate β as control parameter.


### 3.2. Homoclinic Bifurcation of System (3) about Parameter δ

By above analysis, we know there exists aβ
*
∈
(
0,1
) such that system (2) has an order one homoclinic cycle and when the parameter β is appropriately changed, the homoclinic cycle disappears and bifurcates a unique stable order one periodic solution. In the following, we will choose δ as a control parameter and study homoclinic bifurcation of system (3) by using the theory of rotated vector fields. For the sake of convenience, we give some properties about rotated vector fields.The perturbed system of system (4) with parameter δ is as follows:
(15)
d
x
d
t
=
P
¯
(
x
,
y
,
δ
)
,
d
y
d
t
=
Q
¯
(
x
,
y
,
δ
)
,
d
x
d
t
=
P
¯
(
x
,
y
,
δ
)
,
(
x
,
y
)
∉
M
{
x
,
y
}
,
Δ
x
=
α
(
x
,
y
)
,
Δ
y
=
β
(
x
,
y
)
,
(
x
,
y
)
∈
M
{
x
,
y
}
.
Let
(16)
Δ
=
|
P
¯
Q
¯
∂
P
¯
∂
δ
∂
Q
¯
∂
δ
|
,
then we have the following definitions.Definition 20 (see [18]).
For any point on the trajectory of system (15), if Δ
>
0, then system (15) constitutes positive rotated vector fields concerning the parameter δ; otherwise, if Δ
<
0, system (15) constitutes negative rotated vector fields.Lemma 21 (see [18]).
In the positive (negative) rotated vector fields of system (15), the rotated direction of vector fields is counterclockwise (clockwise) when parameter δ changes from δ
=
0 to δ
>
0.Theorem 22.
If4
a
D
2
<
μ
2
<
(
16
/
3
)
a
D
2, x
2
<
K
<
(
2
μ
-
μ
2
-
4
a
D
2
)
/
2
D (or if μ
2
>
4
a
D
2, K
>
(
2
μ
+
μ
2
-
4
a
D
2
)
/
2
D), x
1
<
h
<
h
+
τ
<
x
2, K
≤
3
a, β
=
β
*, δ
>
0, and δ
≪
1, then the system (3) has a unique order one periodic solution. Furthermore, it is orbitally asymptotically stable.Proof.
From the discussion in Theorem19, we know that when β
=
β
* and δ
=
0, the system (3) has an order one homoclinic circle B
E
2
O
A
^
∪
A
B
¯.
For the system (3), we denote P
(
x
,
y
)
=
r
x
(
1
-
(
x
/
K
)
)
-
x
y
/
(
a
+
x
2
) and  Q
(
x
,
y
)
=
y
(
μ
x
/
(
a
+
x
2
)
-
D
)
+
δ
[
r
(
1
-
(
x
/
K
)
)
-
y
/
(
a
+
x
2
)
]
2. By simple calculation, we have
(17)
Δ
=
|
P
Q
∂
P
∂
δ
∂
Q
∂
δ
|
=
x
[
r
(
1
-
x
K
)
-
y
a
+
x
2
]
3
.
Divided by the vertical isocline [
r
(
1
-
(
x
/
K
)
)
-
y
/
(
a
+
x
2
)
]
=
0, the rotated direction of vector fields below the vertical isocline is counterclockwise, but above the vertical isocline, the rotated direction of vector fields of system (3) is clockwise. When the parameter δ changes from 0 (the value in which heteroclinic cycle exists) to 0
<
δ
1
≪
1, both the unstable and stable flow of (
x
2
,
y
2
) move away from their original positions. We denote the new intersections of the unstable and stable flow of (
x
2
,
y
2
) and the impulse set x
=
h and phase set x
=
h
+
τ are A
¯ and O
¯, and A
¯
1 and B
¯, respectively (see Figure 6). According to Theorem 19, we have (
1
-
β
*
)
y
A
=
y
B; then there must exist a point F
¯ which is in the impulse set x
=
h and below the point A such that (
1
-
β
*
)
y
F
¯
=
y
B
¯. Obviously, system (3) must have a trajectory passing through the point F
¯, and we suppose it intersects with the phase set x
=
h
+
τ at point F.
Obviously, pointB
¯ is the order one successor point of point F, then the order one successor function F
1
(
F
)
=
y
B
¯
-
y
F
<
0. Besides, select a point B
¯
0 next to B
¯ in the phase set x
=
h
+
τ; the trajectory of the system (3) from point B
¯
0 must intersect the impulse set x
=
h at a point B
¯
0
′ which is next to point A
¯ and above point A, and after impulsive effect, point B
¯
0
′ is mapped to a point B
¯
1 which is the order one successor point of B
¯
0, then we have y
B
¯
1
=
(
1
-
β
*
)
y
B
¯
0
′
>
(
1
-
β
*
)
y
A
=
y
B. That is to say point B
¯
1 is above point B. Then we have the successor function F
1
(
B
¯
0
)
=
y
B
¯
1
-
y
B
0
>
0. By Lemma 12, we know that there must exist a point M in the phase set x
=
h
+
τ which is between the points B
¯ and F such that F
1
(
M
)
=
0. Then we know the system (3) has an order one periodic solution which passes through the point M.
Similar to the discussion of Theorems16 and 17, we can prove the uniqueness and stability of the order one periodic solution. That completes the proof.Figure 6
Illustration of homoclinic bifurcation of system (3) about δ.


## 3.1. Periodic Solution and Homoclinic Bifurcation of System (2) about Parameter β

According to the model (2), the threshold h and the recruitment τ of the prey should satisfy the condition 0
<
h
<
h
+
τ
<
K by ecological significance. For this consideration, we have the following results.Theorem 16.
If4
a
D
2
<
μ
2
≤
(
16
/
3
)
a
D
2, x
1
<
K
<
x
2 (or if μ
2
>
(
16
/
3
)
a
D
2, x
1
<
K
<
-
x
1
+
2
μ
x
1
/
D), and x
1
<
h
<
h
+
τ
<
K, then there must exist a fixed value β
0
∈
(
0,1
) such that for every β
∈
(
β
0
,
1
), the system (2) has a unique order one periodic solution in region Ω
1, where region Ω
1 is the region enclosed by the x-axis, the impulse set x
=
h, and the unstable flow of the saddle (
K
,
0
).Proof.
According to Lemma14, the system (1) has three equilibria: two hyperbolic saddles (
0,0
) and (
K
,
0
) and a globally asymptotically stable equilibrium (
x
1
,
y
1
) in the interior of the first quadrant. For convenience, we suppose the x-axis intersects the impulse set x
=
h and the phase set x
=
h
+
τ at point A
′ and point B, respectively. The unstable flow of (
K
,
0
) also passes through the impulse set and the phase set, we denote the intersections by point A and point O, respectively. Besides, we suppose the vertical isocline d
x
/
d
t
=
0 intersects the impulse set and the phase set at point C and point D, respectively. Then the region Ω
1 is the interior of the closed curve E
0
O
A
^
∪
A
C
A
′
¯
∪
A
′
B
E
0
¯, where E
0 denotes the equilibrium (
K
,
0
) (see Figure 2).
According to the impulsive conditions of the system (2), there must exist a fixed value β
0
∈
(
0,1
), when β
=
β
0, point A is mapped to the point D after impulsive effect, that is to say, (
1
-
β
0
)
y
A
=
y
D.
Whenβ
0
<
β
<
1, after impulsive effect, point A is mapped to a point between points B and D. Select a point B
0 next to B in the phase set x
=
h
+
τ, the trajectory of the system (2) from point B
0 must intersect the impulse set x
=
h at a point B
0
′ which is next to point A, and after impulsive effect, point B
0
′ is mapped to a point B
1 which is the order one successor point of B
0; then, we have y
B
1
=
(
1
-
β
)
y
B
0
′
<
(
1
-
β
)
y
A
<
y
D. That is to say point B
1 is between points D and B
0. The trajectory from point B
1 must intersect the impulse set x
=
h again at a point B
1
′, and after impulsive effect, point B
1
′ is mapped to a point B
2. Since distinct trajectories do not intersect, we can easily have y
C
<
y
B
1
′
<
y
B
0
′ and y
B
2
=
(
1
-
β
)
y
B
1
′
<
(
1
-
β
)
y
B
0
′
=
y
B
1. Obviously, point B
2 is the order one successor point of point B
1, and then we have the following results of the order one successor function:
(8)
F
1
(
B
0
)
=
y
B
1
-
y
B
0
>
0
,
F
1
(
B
1
)
=
y
B
2
-
y
B
1
<
0
.
By Lemma 12, we know that there must exist a point M in the phase set x
=
h
+
τ which is between the points B
0 and B
1 such that F
1
(
M
)
=
0. Then we know the system (2) has an order one periodic solution which passes through the point M.
In the following, we prove the uniqueness of the order one periodic solution. Arbitrarily choose two pointsA
1 and A
2 in the phase set x
=
h
+
τ, where 0
≤
y
A
1
<
y
A
2
≤
y
D. The trajectories of the system (2) through points A
1 and A
2 must intersect the impulse set x
=
h at some points A
1
′ and A
2
′, respectively, and after impulsive effect, the points A
1
′ and A
2
′ must be mapped to the phase set x
=
h
+
τ at some points A
1
′′ and A
2
′′, respectively (see Figure 3). Because distinct trajectories do not intersect, we can easily have y
C
<
y
A
2
′
<
y
A
1
′
<
y
A and y
A
2
′′
=
(
1
-
β
)
y
A
2
′
<
y
A
1
′′
=
(
1
-
β
)
y
A
1
′; then we get the order one successor functions must satisfy
(9)
F
1
(
A
2
)
-
F
1
(
A
1
)
=
(
y
A
2
′′
-
y
A
2
)
-
(
y
A
1
′′
-
y
A
1
)
=
(
y
A
2
′′
-
y
A
1
′′
)
+
(
y
A
1
-
y
A
2
)
<
0
,
which means the order one successor function F
1 is monotone decreasing in the line segment B
D
¯, thus there exists only one point M
∈
B
D
¯ such that F
1
(
M
)
=
0.
Besides, for any pointH
∈
D
O
¯, it is easy to know that F
1
(
H
)
<
0. So the system (2) has only one order one periodic solution in the region Ω
1. That completes the proof.Figure 2
Existence of order one periodic solution of (2).Figure 3
Uniqueness of order one periodic solution of (2).Theorem 17.
Under the conditions of Theorem16, the order one periodic solution of the system (2) is orbitally asymptotically stable.Proof.
According to Theorem16, the system (2) has a unique order one periodic solution that passes through the point M which is in the phase set x
=
h
+
τ, where y
B
0
<
y
M
<
y
B
1. The trajectory of the system (2) through point B
1 intersects the impulse set x
=
h at point B
1
′, and after impulsive effect, the point B
1
′ is mapped to point B
2, where y
B
<
y
B
2
<
y
M. Besides, the trajectory of the system (2) from point B
2 must intersect the impulse set x
=
h again at a point B
2
′, and after impulsive effect, the point B
2
′ is mapped to a point B
3 in the phase set x
=
h
+
τ, where y
B
1
>
y
B
3
>
y
M (see Figure 4).
Repeat the above steps, the trajectory from pointB
0 will come across impulsive effect infinite times. Denote the phase point corresponding to the ith impulsive effect by B
i, i
=
1,2
,
…. Obviously, B
i is also the order i successor point of point B
0. Then we have
(10)
y
B
0
<
y
B
2
<
y
B
4
<
⋯
<
y
B
2
k
<
y
B
2
(
k
+
1
)
<
⋯
<
y
M
,
y
B
1
>
y
B
3
>
y
B
5
>
⋯
>
y
B
2
k
+
1
>
y
B
2
(
k
+
1
)
+
1
>
⋯
>
y
M
.
Thus {
y
B
2
k
}, k
=
0,1
,
2
,
…is a monotonically increasing sequence, and {
y
B
2
k
+
1
}, k
=
0,1
,
2
,
…is a monotonically decreasing sequence. Because the order one successor function F
1 is monotone decreasing in the line segment B
D
¯, we have
(11)
F
1
(
B
0
)
=
y
B
1
-
y
B
0
>
F
1
(
B
2
)
=
y
B
3
-
y
B
2
>
F
1
(
B
4
)
=
y
B
5
-
y
B
4
>
⋯
>
0
,
and furthermore,
(12)
y
B
2
k
⟶
y
M
,
as
k
⟶
∞
;
y
B
2
k
+
1
⟶
I
M
,
as
k
⟶
∞
.
Pick any point Q
0
∈
B
0
B
1
¯ different from the point M; without loss of generality, we assume y
B
0
<
y
Q
0
<
y
M (otherwise, y
M
<
y
Q
0
<
y
B
1; the discussions are similar), there must exist a natural number n
0 such that y
B
2
n
0
<
y
Q
0
<
y
B
2
(
n
0
+
1
). The trajectory from point Q
0 will also undergo impulsive effect infinite times, and we denote the phase point corresponding to the kth impulsive effect as Q
k, k
=
1,2
,
…. For any natural number l, we have y
B
2
(
n
0
+
l
)
<
y
Q
2
l
<
y
B
2
(
n
0
+
l
+
1
) and y
B
2
(
n
0
+
l
+
1
)
+
1
<
y
Q
2
l
+
1
<
y
B
2
(
n
0
+
l
)
+
1. Then {
y
Q
2
l
}, l
=
0,1
,
2
,
…  is also a monotonically increasing sequence, {
y
Q
2
l
+
1
}, l
=
0,1
,
2
,
… is also a monotonically decreasing sequence, and we get
(13)
y
Q
2
l
⟶
y
M
,
as
l
⟶
∞
;
y
Q
2
l
+
1
⟶
I
M
,
as
l
⟶
∞
.
This indicates that the trajectory of the system (2) from point Q
0
∈
B
0
B
1
¯ ultimately tends to the trajectory passing through point M, which means the order one periodic solution is orbitally asymptotically stable. That completes the proof.Figure 4
Illustration of the orbital asymptotic stability of the order one periodic solution of the system (2).Theorem 18.
Under the conditions of Theorem16, the system (2) has no order k periodic solution in region Ω
1, where k
≥
2.Proof.
For any pointS
∈
B
1
O
¯, y
B
1
<
y
S
<
y
O, the trajectory from point S will undergo impulsive effect infinite times. Denote the impulsive point and phase point corresponding to the k
th
(
k
=
1,2
,
…
) impulsive effect by S
k
′ and S
k, respectively. Similar to the discussions in Theorem 17, we can easily know y
C
<
y
S
k
′
<
y
B
0
′, k
=
1,2
,
…, then we have y
S
k
=
(
1
-
β
)
y
S
k
′
<
(
1
-
β
)
y
B
0
′
=
y
B
1
<
y
S. Obviously, point S
k is the order k successor point of point S, then the order k successor function F
k
(
S
)
=
y
S
k
-
y
S
≠
0. By Lemma 13, there does not exist an order k
(
k
=
1,2
,
…
) periodic solution through point S.
For any pointS
∈
B
M
¯
,
y
B
<
y
S
<
y
M, according to the proof of Theorem 17, there must exist a natural number n
0 such that y
B
2
n
0
<
y
S
<
y
B
2
(
n
0
+
1
). Denote the order k
(
k
=
1,2
,
…
) successor point of point S by point S
k, then we have y
B
2
(
n
0
+
l
)
<
y
S
2
l
<
y
B
2
(
n
0
+
l
+
1
)
<
y
M and y
M
<
y
B
2
(
n
0
+
l
+
1
)
+
1
<
y
S
2
l
+
1
<
y
B
2
(
n
0
+
l
)
+
1. Then {
y
S
2
l
}, l
=
0,1
,
2
,
…is amonotonically increasing sequence where S
0
=
S, {
y
S
2
l
+
1
}, l
=
0,1
,
2
,
…is a monotonically decreasing sequence, and
(14)
y
S
2
l
⟶
y
M
,
as
l
⟶
∞
;
y
S
2
l
+
1
⟶
I
M
,
as
l
⟶
∞
.
When k
=
2
l, we have y
M
>
y
S
k
>
y
S
0
=
y
S and the order k successor function F
k
(
S
)
=
y
S
k
-
y
S
>
0. When k
=
2
l
+
1, we have y
S
k
>
y
M
>
y
S
0
=
y
S and F
k
(
S
)
=
y
S
k
-
y
S
>
0. That is to say there does not exist an order k
(
k
=
1,2
,
…
) periodic solution passing through point S.
Analogously, for any pointS
∈
M
B
1
¯
,
y
M
<
y
S
<
y
B
1, we can prove there does not exist an order k
(
k
=
1,2
,
…
) periodic solution passing through point S.
From the above discussion, we know that the order one periodic solution is the unique periodic solution of the system (2) and there is no order k
(
k
≥
2
) periodic solution in region Ω
1. That completes the proof.Theorem 19.
If4
a
D
2
<
μ
2
<
(
16
/
3
)
a
D
2, x
2
<
K
<
(
2
μ
-
μ
2
-
4
a
D
2
)
/
2
D (or if μ
2
>
4
a
D
2, K
>
(
2
μ
+
μ
2
-
4
a
D
2
)
/
2
D), x
1
<
h
<
h
+
τ
<
x
2 and K
≤
3
a, then there must exist two fixed values β
0 and β
* which satisfy 0
<
β
0
<
β
*
<
1 such that(1)
whenβ
=
β
*, the system (2) has an order one homoclinic cycle;
(2)
ifβ
0
<
β
<
β
*, then the homoclinic cycle of system (2) disappears and bifurcates an order one periodic solution in region Ω
2, where region Ω
2 is the region enclosed by the x-axis, the impulse set x
=
h, and the unstable flow of the saddle (
x
2
,
y
2
). Furthermore, the order one periodic solution is unique and is orbitally asymptotically stable;
(3)
ifβ
*
<
β
<
1, then the homoclinic cycle of system (2) also disappears and there is no periodic solution in region Ω
2.Proof.
ForK
≤
3
a, it is easy to know that y
1
>
y
2. According to Lemma 15, the system (1) has four equilibria: two hyperbolic saddles (
0,0
) and (
x
2
,
y
2
), a hyperbolic stable node (
K
,
0
), and a node or focus (
x
1
,
y
1
). Obviously, the x-axis, the vertical isocline d
x
/
d
t
=
0, and the unstable and stable flow of (
x
2
,
y
2
) are all intersected with the impulse set x
=
h and phase set x
=
h
+
τ. For convenience, we denote the intersections by points A
′ and B
′, C and D, A and O, and A
1 and B, respectively (see Figure 5). Then the region Ω
2 is the interior of the closed curve E
0
E
2
O
A
^
∪
A
C
A
′
¯
∪
A
′
B
′
E
0
¯, where E
0 and E
2 denote the equilibria (
K
,
0
) and (
x
2
,
y
2
), respectively. According to the impulsive conditions of the system (2), there must exist a fixed value β
0
∈
(
0,1
), when β
=
β
0, point A is mapped to the point D after impulsive effect. Also, there must exist a fixed value β
*
∈
(
β
0
,
1
), when β
=
β
*, point A is mapped to the point B after impulsive effect. Thus we have (
1
-
β
0
)
y
A
=
y
D and (
1
-
β
*
)
y
A
=
y
B.
Whenβ
=
β
*, we have (
1
-
β
)
y
A
=
y
B; we know the curve B
E
2
O
A
^
∪
A
B
¯ is an order one circle which has the saddle E
2
(
x
2
,
y
2
) in it; that is to say, the system (2) has an order one homoclinic cycle.
Ifβ
∈
(
β
0
,
β
*
), after impulsive effect, point A is mapped to a point between points B and D. Similar to the discussion of Theorems 16 and 17, we can prove that the system (2) has an order one periodic solution which passes through a point M
∈
B
D
¯ and is orbitally asymptotically stable. Besides, for any point H
∈
B
′
B
¯, the trajectory of the system (2) from point H must ultimately tend to the equilibrium (
K
,
0
) and does not undergo any impulsive effect, that is to say the system (2) has no periodic solution passing through the point H; so we can know the system (2) has a unique order one periodic solution in the region Ω
2 which is orbitally asymptotically stable.
Whenβ
∈
(
β
*
,
1
), we have 0
<
(
1
-
β
)
y
A
<
y
B. Then the trajectory of the system (2) through point H
∈
O
B
¯ must tend to the equilibrium (
K
,
0
) after coming across once impulsive effect and the trajectory passing through point H
∈
B
′
B
¯ must tend to the equilibrium (
K
,
0
) without undergoing any impulsive effect. So we get that the system (2) has no periodic solution in region Ω
2. That completes the proof.Figure 5
Illustration of homoclinic bifurcation of system (2) by choosing the harvesting rate β as control parameter.


## 3.2. Homoclinic Bifurcation of System (3) about Parameter δ

By above analysis, we know there exists aβ
*
∈
(
0,1
) such that system (2) has an order one homoclinic cycle and when the parameter β is appropriately changed, the homoclinic cycle disappears and bifurcates a unique stable order one periodic solution. In the following, we will choose δ as a control parameter and study homoclinic bifurcation of system (3) by using the theory of rotated vector fields. For the sake of convenience, we give some properties about rotated vector fields.The perturbed system of system (4) with parameter δ is as follows:
(15)
d
x
d
t
=
P
¯
(
x
,
y
,
δ
)
,
d
y
d
t
=
Q
¯
(
x
,
y
,
δ
)
,
d
x
d
t
=
P
¯
(
x
,
y
,
δ
)
,
(
x
,
y
)
∉
M
{
x
,
y
}
,
Δ
x
=
α
(
x
,
y
)
,
Δ
y
=
β
(
x
,
y
)
,
(
x
,
y
)
∈
M
{
x
,
y
}
.
Let
(16)
Δ
=
|
P
¯
Q
¯
∂
P
¯
∂
δ
∂
Q
¯
∂
δ
|
,
then we have the following definitions.Definition 20 (see [18]).
For any point on the trajectory of system (15), if Δ
>
0, then system (15) constitutes positive rotated vector fields concerning the parameter δ; otherwise, if Δ
<
0, system (15) constitutes negative rotated vector fields.Lemma 21 (see [18]).
In the positive (negative) rotated vector fields of system (15), the rotated direction of vector fields is counterclockwise (clockwise) when parameter δ changes from δ
=
0 to δ
>
0.Theorem 22.
If4
a
D
2
<
μ
2
<
(
16
/
3
)
a
D
2, x
2
<
K
<
(
2
μ
-
μ
2
-
4
a
D
2
)
/
2
D (or if μ
2
>
4
a
D
2, K
>
(
2
μ
+
μ
2
-
4
a
D
2
)
/
2
D), x
1
<
h
<
h
+
τ
<
x
2, K
≤
3
a, β
=
β
*, δ
>
0, and δ
≪
1, then the system (3) has a unique order one periodic solution. Furthermore, it is orbitally asymptotically stable.Proof.
From the discussion in Theorem19, we know that when β
=
β
* and δ
=
0, the system (3) has an order one homoclinic circle B
E
2
O
A
^
∪
A
B
¯.
For the system (3), we denote P
(
x
,
y
)
=
r
x
(
1
-
(
x
/
K
)
)
-
x
y
/
(
a
+
x
2
) and  Q
(
x
,
y
)
=
y
(
μ
x
/
(
a
+
x
2
)
-
D
)
+
δ
[
r
(
1
-
(
x
/
K
)
)
-
y
/
(
a
+
x
2
)
]
2. By simple calculation, we have
(17)
Δ
=
|
P
Q
∂
P
∂
δ
∂
Q
∂
δ
|
=
x
[
r
(
1
-
x
K
)
-
y
a
+
x
2
]
3
.
Divided by the vertical isocline [
r
(
1
-
(
x
/
K
)
)
-
y
/
(
a
+
x
2
)
]
=
0, the rotated direction of vector fields below the vertical isocline is counterclockwise, but above the vertical isocline, the rotated direction of vector fields of system (3) is clockwise. When the parameter δ changes from 0 (the value in which heteroclinic cycle exists) to 0
<
δ
1
≪
1, both the unstable and stable flow of (
x
2
,
y
2
) move away from their original positions. We denote the new intersections of the unstable and stable flow of (
x
2
,
y
2
) and the impulse set x
=
h and phase set x
=
h
+
τ are A
¯ and O
¯, and A
¯
1 and B
¯, respectively (see Figure 6). According to Theorem 19, we have (
1
-
β
*
)
y
A
=
y
B; then there must exist a point F
¯ which is in the impulse set x
=
h and below the point A such that (
1
-
β
*
)
y
F
¯
=
y
B
¯. Obviously, system (3) must have a trajectory passing through the point F
¯, and we suppose it intersects with the phase set x
=
h
+
τ at point F.
Obviously, pointB
¯ is the order one successor point of point F, then the order one successor function F
1
(
F
)
=
y
B
¯
-
y
F
<
0. Besides, select a point B
¯
0 next to B
¯ in the phase set x
=
h
+
τ; the trajectory of the system (3) from point B
¯
0 must intersect the impulse set x
=
h at a point B
¯
0
′ which is next to point A
¯ and above point A, and after impulsive effect, point B
¯
0
′ is mapped to a point B
¯
1 which is the order one successor point of B
¯
0, then we have y
B
¯
1
=
(
1
-
β
*
)
y
B
¯
0
′
>
(
1
-
β
*
)
y
A
=
y
B. That is to say point B
¯
1 is above point B. Then we have the successor function F
1
(
B
¯
0
)
=
y
B
¯
1
-
y
B
0
>
0. By Lemma 12, we know that there must exist a point M in the phase set x
=
h
+
τ which is between the points B
¯ and F such that F
1
(
M
)
=
0. Then we know the system (3) has an order one periodic solution which passes through the point M.
Similar to the discussion of Theorems16 and 17, we can prove the uniqueness and stability of the order one periodic solution. That completes the proof.Figure 6
Illustration of homoclinic bifurcation of system (3) about δ.


## 4. Numerical Simulations and Discussions

In this paper, we have proposed two predator-prey models with nonmonotonic functional response and state dependent impulsive harvesting. If we take no account of the predator harvesting, Ruan and Xiao (2001) showed that the system (1) has no periodic solution under the conditions listing in Theorem 16. Our numerical simulations also confirm such conclusions (see Figure 7). However, under the same conditions, if we further consider an impulsive predator harvesting, we get that the system (2) has a periodic solution in some cases. We also showed that the existence is mainly dependent on the harvesting rate β. We found that there exists a threshold β
0
∈
(
0,1
), and the system must have a unique orbitally asymptotically stable order one periodic solution when the harvesting rate β is higher than β
0. When the periodic solution exists, our numerical simulations show that, with a bigger harvesting rate, the solution has a longer period and a larger amplitude (see Figure 8).Figure 7
The time series and the portrait phase of the system (1) when r
=
0.5, K
=
15, a
=
100, μ
=
0.5, D
=
0.01, and (
x
(
0
)
,
y
(
0
)
)
=
(
12,20
).Figure 8
The time series and the portrait phase of the system (2) when r
=
0.5, K
=
15, a
=
100, μ
=
0.5, D
=
0.01, h
=
6, τ
=
4, and (
x
(
0
)
,
y
(
0
)
)
=
(
12,20
).Besides, we also prove that the system (2) does not have order k, (
k
=
2,3
,
…
) periodic solutions in Theorem 18; then the order one periodic solution is the unique periodic solution. These results illustrate that if we can satisfy the conditions listing in Theorem 16, we can not only form a sustainable ecological system but also obtain large amounts of valuable biological resources.Under the conditions listing in Theorem19, the system (1) also does not have any periodic solution (see Figure 9) and we also give the conditions under which the system (2) has an order one periodic solution. Besides, according to the conclusions of Theorem 19, we can choose the parameter β as a bifurcation parameter such that the impulsive differential equation system (2) exhibits the phenomenon of homoclinic bifurcation. Similar to many ordinary differential equation systems, we can see that there exists a bifurcation point β
=
β
* for the system (2). When β
=
β
*, the system (2) has an order one homoclinic cycle which is the unique order one circle in region Ω
2. When β gradually changes from β
=
β
* to β
0
<
β
<
β
*, the order one homoclinic cycle is broken and a new order one periodic solution which is orbitally asymptotically stable is generated (see Figure 10). When β gradually changes from β
=
β
* to β
>
β
*, the order one homoclinic cycle is also broken but no new order one periodic solution is generated at the same time; that is to say, the system (2) will have no periodic solution (see Figure 11). These results illustrate that it demands reasonable control of the harvest yield to form a good ecological environment.Figure 9
The time series and the portrait phase of the system (1) when r
=
0.5, K
=
15, a
=
100, μ
=
0.8, D
=
0.038, and (
x
(
0
)
,
y
(
0
)
)
=
(
12,20
).Figure 10
Time series and portrait phase of the system (2) when r
=
0.5, K
=
15, a
=
100, μ
=
0.8, D
=
0.038, h
=
9, τ
=
2, β
=
0.6, and (
x
(
0
)
,
y
(
0
)
)
=
(
12,20
).Figure 11
Time series and portrait phase of the system (2) when r
=
0.5, K
=
15, a
=
100, μ
=
0.8, D
=
0.038, h
=
9, τ
=
2, β
=
0.8, and (
x
(
0
)
,
y
(
0
)
)
=
(
12,20
).



---

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