# Geometric Analysis of an Integrated Pest Management Model Including Two State Impulses

Authors: Wencai Zhao; Yulin Liu; Tongqian Zhang; Xinzhu Meng
Journal: Abstract and Applied Analysis
 (2014)

---

## Abstract

According to integrated pest management strategies, we construct and investigate the dynamics of a Holling-Tanner predator-prey system with state dependent impulsive effects by releasing natural enemies and spraying pesticide at different thresholds. Applying the Dulacs criterion, the global stability of the positive equilibrium in the system without impulsive effect is discussed. By using impulsive differential equation geometry theory and the method of successor functions, we prove the existence of periodic solution of the system with state dependent impulsive effects. Furthermore, the stability conditions of periodic solutions are obtained. Some simulations are exerted to illustrate the feasibility of our main results.


---

## Body

## 1. Introduction

In 2012, the main maize area in Northern China embraced the large-scale plant diseases and insect pests [1]. With the good weather condition,Mythimna separata Walker,Ostrinia furnacalis, aphid, and other pests had a mass propagation which brought a large damage to the production of corn. How human beings effectively control the pests has been a significant task. Due to its simple operation and quick effect, spraying pesticide has always been the major way to kill pests for a long time. However, there may be pesticide residues in vegetables and crops threatening people’s good health and damaging the environment. With the increasing awareness of the environment, people are paying more attention to developing green agriculture. Releasing natural enemies and artificial capture are becoming significant means of controlling pests.Different control strategies will be applied due to different behavior features of pests and their damages for crops in different stages. For instance, in the spawning period ofOstrinia furnacalis, artificial releasing ofTrichogramma (the natural enemy ofOstrinia furnacalis) will be applied with the control effect of 70%–80%; if the hatchability is over 30%, pesticide control will be used instantly [2]. In consideration of the rapidity of chemical control and nonpollution of biological control, people control pests integrating with biological, physical, and chemical means under the EIL (economic injury level) to realize environmental, economic, and social profits together. Spraying pesticides and releasing natural enemies are instantaneous; these phenomena can be described as the impulsive differential equations. In recent decades, the theoretical research on the impulsive differential equation has represented a significant development and has been widely used in various mathematical ecological models [3–11] and many scholars made a deep analysis of the impulsive differential ecological system at a fixed time and have got some important products [12–26]. However, in the actual process of pest control, relevant measures will be used according to pest quantity and its damage to crops, which is the state impulsive differential system. Tang et al. [27, 28], Zeng et al. [29], Zhao et al. [30], Nie et al. [31, 32], and literatures [33–37] had a further exploration for the system and made great progress. Based on the study of a Holling-Tanner system, an integrated pest control model with two impulses is established aiming at the specific pest conditions in different threshold.Letx
=
x
(
t
), y
=
y
(
t
) denote the population densities of pest (prey) and natural enemy (predator) at time t, respectively; then the predator-prey system usually can be expressed as [38]
(1)
d
x
d
t
=
x
g
(
x
)
-
y
Φ
(
x
)
,
d
y
d
t
=
y
[
-
q
(
x
)
+
c
Φ
(
x
)
]
,
in which g
(
x
) denotes the relative growth rate of the prey; Φ
(
x
) is the functional response function of the predator; q
(
x
) is the mortality rate of the prey. The literature [39] studied a class of Holling-Tanner system with functional response function Φ
(
x
)
=
m
x
/
(
A
+
x
). This system was given by
(2)
d
x
d
t
=
r
x
(
1
-
x
K
)
-
m
x
A
+
x
y
,
d
y
d
t
=
y
s
(
1
-
h
y
x
)
,
h
h
h
h
h
x
(
0
)
>
0
,
y
(
0
)
>
0
.
Here, let the prey population be growth in logistic and the environmental capacity is K; the intrinsic growth rate of predator is s and the carrying capacity is proportional to the number of prey. Introducing transformation t
~
=
r
t, x
~
(
t
~
)
=
x
(
t
)
/
K, y
~
(
t
~
)
=
m
y
(
t
)
/
r
K, and letting δ
=
s
/
r, β
=
s
h
/
m, a
=
A
/
K, the system (2) is changed into a dimensionless form:
(3)
d
x
d
t
=
x
(
1
-
x
)
-
x
a
+
x
y
,
d
y
d
t
=
y
(
δ
-
β
y
x
)
,
h
h
h
x
(
0
)
>
0
,
y
(
0
)
>
0
.
In order to carry out integrated control of pests, we adopt strategies as follows.(
1
) When the pest density x
(
t
) reached a lower level x
=
h
1, we release natural enemies to control pests for low damage of insect pests on crops; that is,
(4)
Δ
x
(
t
)
=
0
,
Δ
y
(
t
)
=
λ
,
x
=
h
1
,
where λ is amount of natural enemies y
(
t
) released one time.(
2
) If the pest density x
(
t
) reached a higher level x
=
h
2, due to the fact that the damage of insect pests on crops is severe at this time, we effectively combine spraying pesticides with releasing natural enemies to control pests; that is,
(5)
Δ
x
(
t
)
=
-
p
x
(
t
)
,
Δ
y
(
t
)
=
-
q
y
(
t
)
+
τ
,
x
=
h
2
.
Here, τ is the amount of natural enemies y
(
t
) released one time, p, q are mortality rates of pests and natural enemies which die from spraying pesticides, and p
,
q
∈
(
0,1
).Synthesizing systems (3), (4), and (5), the following integrated pest management model is obtained:
(6)
d
x
(
t
)
d
t
=
x
(
t
)
(
1
-
x
(
t
)
)
-
x
(
t
)
a
+
x
(
t
)
y
(
t
)
,
d
y
(
t
)
d
t
=
y
(
t
)
(
δ
-
β
y
(
t
)
x
(
t
)
)
,
x
≠
h
1
,
h
2
,
or
x
=
h
1
,
y
>
y
~
,
Δ
x
(
t
)
=
0
,
Δ
y
(
t
)
=
λ
,
x
=
h
1
,
y
≤
y
~
,
Δ
x
(
t
)
=
-
p
x
(
t
)
,
Δ
y
(
t
)
=
-
q
y
(
t
)
+
τ
,
x
=
h
2
,
x
(
0
)
>
0
,
y
(
0
)
>
0
,
where x
(
t
) denotes the population density of pests at time t and y
(
t
) denotes the population density of natural enemies at time t; a, δ, β are positive numbers; p, q, λ, τ are control parameters and positive numbers; the point (
h
1
,
y
~
) is the intersection of the isoclinic line y
=
(
1
-
x
)
(
x
+
a
) and the straight line x
=
h
1.


## 2. Preliminaries

In order to analyze the dynamics of the system (6), we introduce the basic knowledge of the state impulsive differential equations.Consider the state impulsive differential equation:(7)
d
x
(
t
)
d
t
=
P
(
x
,
y
)
,
d
y
(
t
)
d
t
=
Q
(
x
,
y
)
,
(
x
,
y
)
∉
M
(
x
,
y
)
,
Δ
x
(
t
)
=
α
(
x
,
y
)
,
Δ
y
(
t
)
=
β
(
x
,
y
)
,
(
x
,
y
)
∈
M
(
x
,
y
)
,
where P
(
x
,
y
) and Q
(
x
,
y
) have order-one continuous partial derivatives.Definition 1.
The dynamic system which is formed by solution mapping of system (7) is called semicontinuous dynamic system, denoted by (
Ω
,
f
,
φ
,
M
), where f is semicontinuous dynamical system mapping and f
:
Ω
→
Ω, φ
(
M
)
=
N, and in which φ is pulse mapping. Here, M
(
x
,
y
) and N
(
x
,
y
) are straight or curved line in the plane. M
(
x
,
y
) is called impulsive set, and N
(
x
,
y
) is called corresponding image set.In the system (6), let M
1
=
{
(
x
,
y
)
∣
x
=
h
1
,
0
<
y
≤
y
~
}, and the image set corresponding to impulsive mapping (4) is N
1
=
φ
1
(
M
1
)
=
{
(
x
,
y
)
∣
x
=
h
1
,
0
<
y
≤
y
~
+
λ
}. Let M
2
=
{
(
x
,
y
)
∣
x
=
h
2
,
y
>
0
}, and the image set corresponding to impulsive mapping (5) is N
2
=
φ
2
(
M
2
)
=
{
(
x
,
y
)
∣
x
(
t
+
)
=
(
1
-
p
)
h
2
,
y
(
t
+
)
=
(
1
-
q
)
y
(
t
)
+
τ
}.Definition 2.
Assume that impulsive setM and image set N are straight lines, the orbit Π
(
A
,
t
) of system (7) starting from point A on N hits M at point A
1 and then jumps onto point A
1
+, then the function f
(
A
)
=
y
A
1
+
-
y
A is defined as a successor function about point A, and then point A
1
+ is called successor point of A.Lemma 3 (see [10, 11]).
Successor function is continuous.Definition 4.
If there exists a pointP
0 on the image set N, and a constant T
>
0 such that Π
(
P
0
,
T
)
=
P
∈
M, φ
(
P
)
=
P
0
∈
N, then the orbit Π
(
P
0
,
t
) starting from P
0 is called an order-one periodic solution of the system (7).Lemma 5 (Bendixson theorem of impulsive differential equations [10, 11]).
Assume thatG is a Bendixson region of system (7); if G does not contain critical points of system (7), then system (7) contains a closed orbit in G.For the system (6), from Lemma 5, the following conclusion is obtained.Lemma 6.
In system (6), if there exist points A and B on the image set N, such that the successor function satisfies f
(
A
)
f
(
B
)
<
0, then there must exist an order-one periodic solution in system (6).Lemma 7 (Analogue of Poincaré Criterion [3, 4]).
Assume thatx
=
ξ
(
t
), y
=
η
(
t
) is the T-periodic solution of the following impulsive differential equations:
(8)
d
x
(
t
)
d
t
=
P
(
x
,
y
)
,
d
y
(
t
)
d
t
=
Q
(
x
,
y
)
,
Φ
(
x
,
y
)
≠
0
,
Δ
x
(
t
)
=
α
(
x
,
y
)
,
Δ
y
(
t
)
=
β
(
x
,
y
)
,
Φ
(
x
,
y
)
=
0
,
where P
(
x
,
y
) and Q
(
x
,
y
) contain order-one continuous partial derivatives and Φ
(
x
,
y
) is a sufficiently smooth function with g
r
a
d
Φ
(
x
,
y
)
≠
0.
If the multiplierμ satisfies the condition |
μ
|
<
1, then the periodic solution (
ξ
(
t
)
,
η
(
t
)
) of the system (8) is orbitally asymptotically stable, where
(9)
μ
=
∏
j
=
1
n
‍
κ
j
exp
⁡
[
∫
0
T
‍
(
∂
P
(
ξ
(
t
)
,
η
(
t
)
)
∂
x
h
μ
h
h
h
h
h
h
h
h
h
h
h
h
h
h
+
∂
Q
(
ξ
(
t
)
,
η
(
t
)
)
∂
y
)
d
t
]
,
with
(10)
κ
j
=
(
(
∂
β
∂
y
·
∂
Φ
∂
x
-
∂
β
∂
x
·
∂
Φ
∂
y
+
∂
Φ
∂
x
)
P
+
h
κ
j
h
h
h
+
(
∂
α
∂
x
·
∂
Φ
∂
y
-
∂
α
∂
y
·
∂
Φ
∂
x
+
∂
Φ
∂
y
)
Q
+
)
×
(
∂
Φ
∂
x
P
+
∂
Φ
∂
y
Q
)
-
1
,
and P, Q, ∂
α
/
∂
x, ∂
α
/
∂
y, ∂
β
/
∂
x, ∂
β
/
∂
y, ∂
Φ
/
∂
x, ∂
Φ
/
∂
y are calculated at the point (
ξ
(
τ
j
)
,
η
(
τ
j
)
), P
+
=
P
(
ξ
(
τ
j
+
)
,
η
(
τ
j
+
)
), Q
+
=
Q
(
ξ
(
τ
j
+
)
,
η
(
τ
j
+
)
), and τ
j
(
j
∈
N
) is the time of the jth jump.


## 3. The Stability of System (6) without Impulsive Effect

In the system (6), if p
=
q
=
λ
=
τ
=
0, that is, the system without impulsive effect, the following system is obtained:
(11)
d
x
d
t
=
x
(
1
-
x
)
-
x
a
+
x
y
,
d
y
d
t
=
y
(
δ
-
β
y
x
)
,
h
h
h
x
(
0
)
>
0
,
y
(
0
)
>
0
.
If set (
x
(
t
)
,
y
(
t
)
) is an arbitrary solution of the system (11) satisfying the initial conditions, then the following lemma is obtained.Lemma 8.
The solutions of the system (11) is bounded, which means ∃
T
>
0 satisfies 0
≤
x
(
t
)
≤
1 and 0
≤
y
(
t
)
≤
δ
/
β for t
≥
T.Obviously, the system (11) exhibits prey isocline L
1
:
y
=
(
1
-
x
)
(
x
+
a
), predator isocline L
2
:
y
=
(
δ
/
β
)
x, nontrivial equilibrium points E
0
(
1,0
), and E
(
x
*
,
y
*
), and here;
(12)
x
*
=
(
β
-
δ
-
a
β
)
+
(
β
-
δ
-
a
β
)
2
+
4
a
β
2
2
β
,
y
*
=
δ
β
x
*
.
Calculating the variational matrix of the equilibrium point in the system (11), we get
(13)
J
(
E
0
)
=
(
-
1
-
1
1
+
a
0
δ
)
.
Obviously, E
0 is saddle point. At E,
(14)
J
(
E
)
=
(
x
*
x
*
+
a
(
1
-
a
-
2
x
*
)
-
x
*
x
*
+
a
δ
2
β
-
δ
)
.
The characteristic equation of J
(
E
) is
(15)
f
(
λ
)
=
λ
2
+
p
λ
+
q
=
0
,
in which p
=
δ
-
(
x
*
/
(
x
*
+
a
)
)
(
1
-
a
-
2
x
*
), q
=
(
δ
x
*
/
(
x
*
+
a
)
)
(
δ
/
β
-
(
1
-
a
-
2
x
*
)
).Thus,(16)
Δ
=
p
2
-
4
q
=
(
δ
-
x
*
x
*
+
a
(
1
-
a
-
2
x
*
)
)
2
-
4
δ
x
*
x
*
+
a
(
δ
β
-
(
1
-
a
-
2
x
*
)
)
,
λ
1
λ
2
=
δ
x
*
x
*
+
a
(
δ
β
-
(
1
-
a
-
2
x
*
)
)
=
δ
x
*
+
a
(
x
*
2
+
a
)
>
0
,
λ
1
+
λ
2
=
-
(
δ
-
x
*
x
*
+
a
(
1
-
a
-
2
x
*
)
)
=
-
2
x
*
2
+
(
a
-
1
+
δ
)
x
*
+
a
δ
x
*
+
a
.
Let P
(
x
)
=
2
x
2
+
(
a
-
1
+
δ
)
x
+
a
δ, and then
(17)
Δ
=
P
2
(
x
*
)
(
x
*
+
a
)
2
-
4
δ
(
x
*
2
+
a
)
x
*
+
a
=
1
(
x
*
+
a
)
2
[
P
2
(
x
*
)
-
4
δ
(
x
*
2
+
a
)
(
x
*
+
a
)
]
.
Based on the above analysis, we can get the following conclusion.Theorem 9.
IfP
(
x
*
)
>
0, the positive equilibrium point E
(
x
*
,
y
*
) of the system (11) is locally asymptotically stable. Specially,(1)
if(
H
1
)
:
0
<
P
(
x
*
)
<
(
4
δ
(
x
*
2
+
a
)
(
x
*
+
a
)
)
1
/
2, E
(
x
*
,
y
*
) is a locally asymptotically stable focus,
(2)
if(
H
2
)
:
P
(
x
*
)
⩾
(
4
δ
(
x
*
2
+
a
)
(
x
*
+
a
)
)
1
/
2, E
(
x
*
,
y
*
) is a locally asymptotically stable node.Next, we discuss the global stability ofE
(
x
*
,
y
*
) about the system (11).Theorem 10.
If(
H
3
)
:
a
+
δ
≥
1 or (
H
4
)
:
1
-
8
a
δ
<
a
+
δ
<
1 is true, then the positive equilibrium E
(
x
*
,
y
*
) of the system (11) is globally asymptotically stable.Proof.
From(
H
3
) or (
H
4
), we have P
(
x
)
=
2
x
2
+
(
a
-
1
+
δ
)
x
+
a
δ
>
0 for x
>
0, and thus P
(
x
*
)
>
0. Structure a Dulac function as follows:
(18)
B
(
x
,
y
)
=
x
+
a
x
y
2
,
x
>
0
,
y
>
0
.
Let f
(
x
,
y
)
=
x
(
1
-
x
)
-
(
x
/
(
a
+
x
)
)
y, g
(
x
,
y
)
=
y
(
δ
-
β
(
y
/
x
)
), and thus
(19)
∂
(
f
B
)
∂
x
+
∂
(
g
B
)
∂
y
=
-
1
x
y
2
[
2
x
2
+
(
a
-
1
+
δ
)
x
+
a
δ
]
∂
(
f
B
)
∂
x
+
∂
(
g
B
)
∂
y
=
-
1
x
y
2
P
(
x
)
⩽
0
.
By the Bendixson-Dulac theorem, there does not exist closed orbit of the system (11) around E. Based on Lemma 8 and Theorem 9, the positive equilibrium E
(
x
*
,
y
*
) is globally asymptotically stable.Remark 11.
If(
H
1
), (
H
3
) or (
H
1
), (
H
4
) are true, then E
(
x
*
,
y
*
) is a globally asymptotically stable focus.Remark 12.
If(
H
2
), (
H
3
) or (
H
2
), (
H
4
) are true, then E
(
x
*
,
y
*
) is a globally asymptotically stable node.
Assume thatE
(
x
*
,
y
*
) is globally asymptotically stable focal point of the system (11), and then the illustration of vector graph of the system is as follows (see Figure 1).Figure 1
Illustration of vector graph of system (11), where a
=
0.05, δ
=
0.5, β
=
0.7.


## 4. The Geometric Analysis of System (6) with Two State Impulses

In this section, we will discuss the existence and stability of periodic solution of system (6) only at focal point situation. So, we assume that the conditions (
H
1
), (
H
3
) or (
H
1
), (
H
4
) are true. According to the practical significance of the integrated pest management model, the condition (
H
5
)
:
h
1
<
(
1
-
p
)
h
2
<
h
2
<
x
* is always given as such. By the analysis of system (6), the curve L
1
:
y
=
(
1
-
x
)
(
x
+
a
) is X-isocline, and the line L
2
:
y
=
(
δ
/
β
)
x is Y-isocline. Let points P, Q, R be the intersection of the curve L
1 and lines x
=
h
1, x
=
(
1
-
p
)
h
2, x
=
h
2, respectively. Obviously E
(
x
*
,
y
*
) is the intersection point of L
1 and L
2. From the previous discussion, we know that the first impulsive set is M
1
=
{
(
x
,
y
)
∣
x
=
h
1
,
0
<
y
≤
y
~
}, and the image set corresponding to M
1 is N
1
=
{
(
x
,
y
)
∣
x
=
h
1
,
0
<
y
≤
y
~
+
λ
}; the second impulsive set is M
2
=
{
(
x
,
y
)
∣
x
=
h
2
,
y
>
0
}, and the image set corresponding to M
2 is N
2
=
{
(
x
,
y
)
∣
x
(
t
+
)
=
(
1
-
p
)
h
2
,
y
(
t
+
)
=
(
1
-
q
)
y
(
t
)
+
τ
}. The structure of the system can be shown as in Figure 2.Figure 2
The structure graph of system (6).Using successor function and geometric theory of impulsive differential equations and according to different positions of orbit initial points, the existence and stability of periodic solution of system (6) are discussed as follows.


### 4.1. The Initial Point onN
1

LetC
0
(
h
1
,
y
C
0
) be an initial point of the orbit of the system (6); if y
C
0
<
y
~, point C
0 is below point P
(
h
1
,
y
~
), then C
0
∈
M
1 (M
1 is impulsive set), and the image point C
0
+ of C
0 must be above point P with n times impulses; therefore, we only need to discuss the cases of y
C
0
>
y
~.The orbitΠ
(
C
0
,
t
) starting from C
0
(
h
1
,
y
C
0
) hits the impulsive set M
1 at point C
1
(
h
1
,
y
C
1
), and then C
1 jumps to point C
11
(
h
1
,
y
C
11
). If y
C
11
<
y
~, C
11 continued to jump to point C
12
(
h
1
,
y
C
12
), and after n times it reaches the point C
1
n
(
h
1
,
y
C
1
n
), where y
C
11
=
y
C
1
+
λ, y
C
12
=
y
C
1
+
2
λ
,
…, y
C
1
n
=
y
C
1
+
n
λ, and y
~
<
y
C
1
n
<
y
~
+
λ. The situation of point C
1
n has three cases as follows.(a) Ify
C
1
n
=
y
C
0(see Figure 3(a)), C
1
n is coincident with C
0, then the curve C
0
C
1
C
11
⋯
C
1
n is closed orbit, and the system (6) exhibits a 1-periodic solution.The orbit starting from the pointC
0 on N
1.
(a)
(b)
(c)(b) Ify
C
1
n
>
y
C
0 (see Figure 3(b)), C
1
n is above C
0; in this time, the successor function of C
0 satisfies f
(
C
0
)
=
y
C
1
n
-
y
C
0
>
0. In the meantime, choose a point D
0
(
h
1
,
y
D
0
) on N
1 satisfying y
D
0
>
y
~
+
λ. The orbit Π
(
D
0
,
t
) starting from D
0 hits the impulsive set at point D
1
(
h
1
,
y
D
1
), and D
1 jumps sometimes to point D
1
m
(
h
1
,
y
D
1
m
), where y
~
<
y
D
1
m
<
y
~
+
λ. Thus, the successor function of D
0 satisfies f
(
D
0
)
=
y
D
1
m
-
y
D
0
<
0.According to Lemma6, the system (6) exhibits a periodic solution, and the initial point of the periodic solution is between C
0 and D
0.(c) Ify
C
1
n
<
y
C
0 (see Figure 3(c)), C
1
n is below C
0; in this time, the successor function of C
0 is f
(
C
0
)
=
y
C
1
n
-
y
C
0
<
0. On the other hand, choose a point D
0
(
h
1
,
y
D
0
) on N
1 satisfying y
~
<
y
D
0
<
y
~
+
ε (ε is a sufficiently small positive number). The orbit Π
(
D
0
,
t
) starting from D
0 hits the impulsive set M
1 at point D
1
(
h
1
,
y
D
1
), and D
1 jumps to point D
11
(
h
1
,
y
D
11
), where y
D
11
=
y
D
1
+
λ. As D
0 sufficiently closed to P, D
1 is sufficiently close to P, then y
D
11
=
y
D
1
+
λ
>
y
D
0. Thus, the successor function of D
0 satisfiesf
(
D
0
)
=
y
D
1
-
y
D
0
>
0.From Lemma6, system (6) exhibits a 1-periodic solution.To sum up the above discussed, we get the following.Theorem 13.
If the initial pointC
0
(
h
1
,
y
c
0
) of the orbit of the system (6) is on N
1 with y
C
0
>
y
~, then the system exhibits a 1-periodic solution.Next, we will discuss the stability of the above periodic solutions.Theorem 14.
Let(
ξ
(
t
)
,
η
(
t
)
) be the T-periodic solution of the system (6) with the initial point C
0
(
h
1
,
η
0
); the closed orbit corresponding to the periodic solution is the curve C
0
C
1
C
11
⋯
C
1
n, if
(20)
|
μ
|
=
|
κ
exp
⁡
{
-
∫
0
T
‍
(
ξ
(
t
)
-
ξ
(
t
)
η
(
t
)
(
ξ
(
t
)
+
a
)
2
+
β
η
(
t
)
ξ
(
t
)
)
d
t
}
|
<
1
,
where
(21)
κ
=
η
0
-
n
λ
η
0
∏
j
=
1
n
‍
1
-
h
1
-
(
1
/
(
h
1
+
a
)
)
(
η
0
-
(
n
-
j
)
λ
)
1
-
h
1
-
(
1
/
(
h
1
+
a
)
)
(
η
0
-
(
n
-
j
+
1
)
λ
)
,
then the periodic solution (
ξ
(
t
)
,
η
(
t
)
) is orbitally asymptotically stable.Proof.
Let the orbitΠ
(
C
0
,
t
) with the initial point C
0
(
h
1
,
η
0
) hit the impulsive set M
1 at C
1
(
ξ
(
T
)
,
η
(
T
)
), and then C
1 jumps to the point C
11
(
ξ
(
τ
1
)
,
η
(
τ
1
)
). The C
11 continued to jump to point C
12
(
ξ
(
τ
2
)
,
η
(
τ
2
)
), and at last, the image point C
12 reaches point C
1
n
(
ξ
(
τ
n
)
,
η
(
τ
n
)
) with n times pulses. Here, η
(
τ
1
)
=
η
(
T
)
+
λ
,
η
(
τ
2
)
=
η
(
T
)
+
2
λ
,
…
,
η
(
τ
n
)
=
η
(
T
)
+
n
λ. For the jth time impulse, obviously ξ
(
τ
j
+
)
=
ξ
(
τ
j
+
1
)
,
η
(
τ
j
+
)
=
η
(
τ
j
+
1
). For the system (6), let P
(
x
,
y
)
=
x
(
1
-
x
)
-
(
x
/
(
a
+
x
)
)
y, Q
(
x
,
y
)
=
y
(
δ
-
β
(
y
/
x
)
), α
(
x
,
y
)
=
0, β
(
x
,
y
)
=
λ, Φ
(
x
,
y
)
=
x
-
h
1, and ξ
(
T
)
=
h
1, η
(
T
)
=
η
0
-
n
λ; therefore we have
(22)
∂
P
∂
x
=
1
-
2
x
-
a
(
a
+
x
)
2
y
,
∂
Q
∂
y
=
δ
-
2
β
y
x
,
∂
α
∂
x
=
∂
α
∂
y
=
0
,
∂
β
∂
x
=
∂
β
∂
y
=
0
,
∂
Φ
∂
x
=
1
,
∂
Φ
∂
y
=
0
.
According to Lemma 7, we get
(23)
κ
j
=
(
(
∂
β
∂
y
·
∂
Φ
∂
x
-
∂
β
∂
x
·
∂
Φ
∂
y
+
∂
Φ
∂
x
)
P
+
h
h
h
+
(
∂
α
∂
x
·
∂
Φ
∂
y
-
∂
α
∂
y
·
∂
Φ
∂
x
+
∂
Φ
∂
y
)
Q
+
)
×
(
∂
Φ
∂
x
P
+
∂
Φ
∂
y
Q
)
-
1
=
P
(
ξ
(
τ
j
+
)
,
η
(
τ
j
+
)
)
P
(
ξ
(
τ
j
)
,
η
(
τ
j
)
)
=
P
(
h
1
,
η
(
T
)
+
j
λ
)
P
(
h
1
,
η
(
T
)
+
(
j
-
1
)
λ
)
=
1
-
h
1
-
(
1
/
(
h
1
+
a
)
)
(
η
0
-
(
n
-
j
)
λ
)
1
-
h
1
-
(
1
/
(
h
1
+
a
)
)
(
η
0
-
(
n
-
j
+
1
)
λ
)
=
1
-
λ
/
(
h
1
+
a
)
1
-
h
1
-
(
1
/
(
h
1
+
a
)
)
(
η
0
-
(
n
-
j
+
1
)
λ
)
,
μ
=
∏
j
=
1
n
‍
κ
j
exp
⁡
[
∫
0
T
‍
(
∂
P
∂
x
+
∂
Q
∂
y
)
d
t
]
=
∏
j
=
1
n
‍
κ
j
exp
⁡
[
∫
0
T
‍
(
1
-
2
x
-
a
(
a
+
x
)
2
y
+
δ
-
2
β
y
x
)
d
t
]
=
∏
j
=
1
n
‍
κ
j
exp
⁡
{
∫
0
T
‍
[
-
x
+
x
y
(
a
+
x
)
2
-
β
y
x
]
(
1
-
x
-
y
x
+
a
)
+
(
δ
-
β
y
x
)
h
h
h
h
h
h
h
h
h
h
h
h
h
h
-
x
+
x
y
(
a
+
x
)
2
-
β
y
x
]
d
t
}
=
∏
j
=
1
n
‍
κ
j
exp
⁡
{
∫
0
T
‍
d
x
x
+
∫
0
T
‍
d
y
y
h
h
h
h
h
h
h
h
h
h
-
∫
0
T
‍
(
x
-
x
y
(
a
+
x
)
2
+
β
y
x
)
d
t
}
=
∏
j
=
1
n
‍
κ
j
η
0
-
n
λ
η
0
exp
⁡
[
-
∫
0
T
‍
(
x
-
x
y
(
a
+
x
)
2
+
β
y
x
)
d
t
]
=
η
0
-
n
λ
η
0
×
∏
j
=
1
n
‍
1
-
h
1
-
(
1
/
(
h
1
+
a
)
)
(
η
0
-
(
n
-
j
)
λ
)
1
-
h
1
-
(
1
/
(
h
1
+
a
)
)
(
η
0
-
(
n
-
j
+
1
)
λ
)
×
exp
⁡
[
-
∫
0
T
‍
(
x
-
x
y
(
a
+
x
)
2
+
β
y
x
)
d
t
]
=
κ
exp
⁡
{
-
∫
0
T
‍
(
ξ
(
t
)
-
ξ
(
t
)
η
(
t
)
(
ξ
(
t
)
+
a
)
2
+
β
η
(
t
)
ξ
(
t
)
)
d
t
}
.
From Lemma 7, if |
μ
|
=
|
κ
exp
⁡
{
-
∫
0
T
‍
(
ξ
(
t
)
-
(
ξ
(
t
)
η
(
t
)
/
(
ξ
(
t
)
+
a
)
2
)
+
(
β
η
(
t
)
/
ξ
(
t
)
)
)
d
t
}
|
<
1, then the periodic solution of the system (6) is orbitally asymptotically stable. This completes the proof.Remark 15.
If1
-
h
1
-
(
η
0
/
(
h
1
+
a
)
)
>
0 and β
≥
1, then the periodic solution with initial point C
0
(
h
1
,
η
0
) (where η
0
>
y
~) is orbitally asymptotically stable.


### 4.2. The Initial Point onN
2

Based on existence and uniqueness theorem of differential equations, there exists a unique pointQ
0
(
(
1
-
p
)
h
2
,
y
Q
0
) on N
2 such that the orbit Π
(
Q
0
,
t
) starting from Q
0 is tangent to N
1 at point P
(
h
1
,
y
~
). Assume that C
0
(
(
1
-
p
)
h
2
,
y
C
0
) is the initial point of the orbit Π
(
C
0
,
t
) of system (6). Next, we will investigate the existence of periodic solution of the system with different positions of C
0 and Q
0. Three cases should be discussed.Case I (
y
C
0
=
y
Q
0
; see Figure4). The initial point C
0 is exactly Q
0.The initial pointC
0 on N
2 (Case  I: y
C
0
=
y
Q
0).
(a)
(b)
(c)The orbitΠ
(
Q
0
,
t
) starting from Q
0 is tangent to N
1 at the point P, and through N
2 hit M
2 at the point Q
1
(
h
2
,
y
Q
1
), and then Q
1 jumps to Q
1
+
(
x
Q
1
+
,
y
Q
1
+
) on N
2. According to (6), the following is obtained:
(24)
x
Q
1
+
=
(
1
-
p
)
h
2
,
y
Q
1
+
=
(
1
-
q
)
y
Q
1
+
τ
.
About the points Q
1
+ and Q
0, there are the following three positional relations.(a)
IfQ
1
+ coincides with Q
0
:
y
Q
1
+
=
y
Q
0 (see Figure 4(a)), then the curve Q
0
P
Q
1
Q
1
+ is closed orbit.
(b)
IfQ
1
+ is below Q
0
:
y
Q
1
+
<
y
Q
0 (see Figure 4(b)), then the successor function of Q
0 satisfies f
(
Q
0
)
=
y
Q
1
+
-
y
Q
0
<
0. In the meantime, take a point S
0
(
(
1
-
p
)
h
2
,
y
S
0
) on N
2 satisfying 0
<
y
S
0
<
ε (ε
>
0 small enough). The orbit Π
(
S
0
,
t
) starting from S
0 hits the impulsive M
2 at the point S
1
(
h
2
,
y
S
1
), and then S
1 jumps to the point S
1
+
(
x
S
1
+
,
y
S
1
+
), where x
S
1
+
=
(
1
-
p
)
h
2
,
y
S
1
+
=
(
1
-
q
)
y
S
1
+
τ. Obviously, the successor function of S
0 is f
(
S
0
)
=
y
S
1
+
-
y
S
0
>
0. From Lemma 6, the system (6) has an order one periodic solution, where the initial point of the periodic solution is between Q
0 and S
0.
(c)
IfQ
1
+ is above Q
0
:
y
Q
1
+
>
y
Q
0 (see Figure 4(c)), the system (6) does not have closed orbit in the area Ω
1
=
{
(
x
,
y
)
∣
h
1
<
x
<
h
2
} in this time.Based on the discussion above, we get the following.Theorem 16.
Assume that the orbitΠ
(
Q
0
,
t
) starting from Q
0
(
(
1
-
p
)
h
2
,
y
Q
0
) is tangent to N
1 at the point P
(
h
1
,
y
~
) and hits the impulsive set M
2 at the point Q
1
(
h
2
,
y
Q
1
), the image point of Q
1 is Q
1
+
(
x
Q
1
+
,
y
Q
1
+
) on N
2. If y
Q
1
+
≤
y
Q
0, the system (6) has 1-periodic solution in the area Ω
1
=
{
(
x
,
y
)
∣
h
1
<
x
<
h
2
}.Case II
(
y
C
0
<
y
Q
0
). The initial point C
0 is below Q
0.The isoclineL
1
:
y
=
(
1
-
x
)
(
x
+
a
) intersects with the phase set N
2 at the point Q
(
(
1
-
p
)
h
2
,
y
Q
), and Q is below Q
0. In this case, we discuss the existence of periodic solution of the system (6) with the example y
C
0
=
y
Q (see Figure 5).The initial pointC
0 on N
2 (Case  II: y
C
0
<
y
Q
0).
(a)
(b)
(c)
(d)
(e)The orbitΠ
(
Q
,
t
) starting from Q moves to the point C
1
(
h
2
,
y
C
1
) on the impulsive set M
2, and C
1 jumps onto C
1
+
(
x
C
1
+
,
y
C
1
+
) on the image set N
2, and then
(25)
y
C
0
=
y
Q
=
(
1
-
(
1
-
p
)
h
2
)
(
a
+
(
1
-
p
)
h
2
)
,
x
C
1
+
=
(
1
-
p
)
h
2
,
y
C
1
+
=
(
1
-
q
)
y
C
1
+
τ
.
(a) Ify
C
1
+
=
y
Q
:
C
1
+ coincide with Q (see Figure 5(a)), the curve Q
C
1
C
1
+ is the closed orbit of the system (6).
(b) Ify
Q
<
y
C
1
+
≤
y
Q
0
:
C
1
+ is between Q and Q
0 (see Figure 5(b)), in this case the successor function of Q is f
(
Q
)
=
y
C
1
+
-
y
Q
>
0. On the other hand, consider the orbit Π
(
Q
0
,
t
) starting from Q
0
,
Π
(
Q
0
,
t
) hits the impulsive set M
2 at Q
1
(
h
2
,
y
Q
1
), and Q
1 jumps onto Q
1
+
(
x
Q
1
+
,
y
Q
1
+
) on N
2. Based on the existence and uniqueness theorem of differential equations, Q
1 must be below C
1, and Q
1
+ must be below C
1
+. Thus y
Q
1
+
<
y
C
1
+
≤
y
Q
0, and the successor function of Q
0 is f
(
Q
0
)
=
y
Q
1
+
-
y
Q
0
<
0. From Lemma 6, the system (6) has 1-periodic solution, and the initial point of the periodic solution is between Q
0 and Q.
Ify
C
1
+
>
y
Q
0
:
C
1
+ is above Q
0, in this case the successor function of Q is f
(
Q
)
=
y
C
1
+
-
y
Q
>
0. There are two different cases (case (c) and case (d)).
(c) Ify
C
1
+
>
y
Q
0 and y
Q
1
+
≤
y
Q
0
:
Q
1
+ is below Q
0 (see Figure 5(c)), f
(
Q
0
)
=
y
Q
1
+
-
y
Q
0
≤
0, the system (6) has closed orbit.
(d) Ify
C
1
+
>
y
Q
0 and y
Q
1
+
>
y
Q
0
:
Q
1
+ is above Q
0 (see Figure 5(d)), the system (6) does not have closed orbit in the area Ω
1
=
{
(
x
,
y
)
∣
h
1
<
x
<
h
2
}.
(e) Ify
C
1
+
<
y
Q
:
C
1
+ is below Q (see Figure 5(e)), in this case the successor function of Q is f
(
Q
)
=
y
C
1
+
-
y
Q
<
0. On the other hand, take a point D
0
(
(
1
-
p
)
h
2
,
y
D
0
) from N
2 satisfying that y
D
0 is sufficiently small number, which is 0
<
y
D
0
<
ε (ε
>
0 small enough). The orbit Π
(
D
0
,
t
) starting from D
0 moves to D
1
(
h
2
,
y
D
1
) on the impulsive set M
2, and D
1 jumps onto D
1
+
(
x
D
1
+
,
y
D
1
+
) on N
2, where x
D
1
+
=
(
1
-
p
)
h
2, y
D
1
+
=
(
1
-
q
)
y
D
1
+
τ. Obviously, the successor function of D
0 is f
(
D
0
)
=
y
D
1
+
-
y
D
0
>
0. From Lemma 6, the system (6) has closed orbit.Based on the discussion above, we get the following.Theorem 17.
Assume that the orbitΠ
(
Q
,
t
) starting from Q
(
(
1
-
p
)
h
2
,
y
Q
) hits the impulsive set M
2 at C
1
(
h
2
,
y
C
1
), and C
1 jumps onto C
1
+
(
x
C
1
+
,
y
C
1
+
) on N
2. The orbit Π
(
Q
0
,
t
) starting from the point Q
0
(
(
1
-
p
)
h
2
,
y
Q
0
) is tangent to N
1 at P
(
h
1
,
y
~
), and hitting the impulsive set M
2 at Q
1
(
h
2
,
y
Q
1
), the image point of Q
1 is Q
1
+
(
x
Q
1
+
,
y
Q
1
+
). Then one has the following. (1)
Ify
C
1
+
≤
y
Q
0, the system (6) has 1-periodic solution.
(2)
Ify
C
1
+
>
y
Q
0 and y
Q
1
+
≤
y
Q
0, the system (6) has 1-periodic solution.Case III (
y
C
0
>
y
Q
0
; see Figure6). The initial point C
0 is above Q
0.Figure 6
The initial pointC
0 on N
2 (Case  III: y
C
0
>
y
Q
0).In this case, the orbitΠ
(
C
0
,
t
) starting from C
0 goes through the isocline L
1 from the left of the line x
=
h
1, hitting the impulsive set M
1 at C
1
(
h
1
,
y
C
1
). The same conclusion can be made as in Section 4.1.Next, we discuss the stability of the periodic solution with the initial point onN
2.Theorem 18.
Assume that(
ξ
(
t
)
,
η
(
t
)
) is the T-periodic solution of the system (6) with initial point C
0
(
(
1
-
p
)
h
2
,
η
0
); and if
(26)
|
μ
|
=
|
κ
exp
⁡
{
-
∫
0
T
‍
(
ξ
(
t
)
-
ξ
(
t
)
η
(
t
)
(
ξ
(
t
)
+
a
)
2
+
β
η
(
t
)
ξ
(
t
)
)
d
t
}
|
<
1
,
where
(27)
κ
=
η
0
-
τ
η
0
·
1
-
(
1
-
p
)
h
2
-
(
η
0
/
(
(
1
-
p
)
h
2
+
a
)
)
1
-
h
2
-
(
(
η
0
-
τ
)
/
(
(
1
-
q
)
(
h
2
+
a
)
)
)
,
the periodic solution (
ξ
(
t
)
,
η
(
t
)
) is orbitally asymptotically stable.Proof.
Assume that the periodic orbitΠ
(
C
0
,
t
) starting from the point C
0
(
(
1
-
p
)
h
2
,
η
0
) moves to the point C
1
(
ξ
(
T
)
,
η
(
T
)
) on impulsive set M
2, and C
1 jumps onto the point C
1
+
(
ξ
(
T
+
)
,
η
(
T
+
)
) on N
2. Therefore, Π
(
C
0
,
T
)
=
C
1, C
1
+
=
φ
2
(
C
1
)
=
C
0, ξ
(
T
+
)
=
(
1
-
p
)
ξ
(
T
), η
(
T
+
)
=
(
1
-
q
)
η
(
T
)
+
τ.
Compared with the system (6), we get
(28)
P
(
x
,
y
)
=
x
(
1
-
x
)
-
x
a
+
x
y
,
Q
(
x
,
y
)
=
y
(
δ
-
β
y
x
)
,
α
(
x
,
y
)
=
-
p
x
,
β
(
x
,
y
)
=
-
q
y
+
τ
,
Φ
(
x
,
y
)
=
x
-
h
2
,
ξ
(
T
)
=
h
2
,
η
(
T
)
=
η
0
-
τ
1
-
q
,
∂
P
∂
x
=
1
-
2
x
-
a
(
a
+
x
)
2
y
,
∂
Q
∂
y
=
δ
-
2
β
y
x
,
∂
α
∂
x
=
-
p
,
∂
α
∂
y
=
0
,
∂
β
∂
x
=
0
,
∂
β
∂
y
=
-
q
,
∂
Φ
∂
x
=
1
,
∂
Φ
∂
y
=
0
.
Thus,
(29)
κ
1
=
(
(
∂
β
∂
y
·
∂
Φ
∂
x
-
∂
β
∂
x
·
∂
Φ
∂
y
+
∂
Φ
∂
x
)
P
+
hhhh
+
(
∂
α
∂
x
·
∂
Φ
∂
y
-
∂
α
∂
y
·
∂
Φ
∂
x
+
∂
Φ
∂
y
)
Q
+
)
×
(
∂
Φ
∂
x
P
+
∂
Φ
∂
y
Q
)
-
1
=
(
1
-
q
)
P
(
ξ
(
T
+
)
,
η
(
T
+
)
)
P
(
ξ
(
T
)
,
η
(
T
)
)
=
(
η
0
(
1
-
p
)
h
2
+
a
(
1
-
p
)
(
1
-
q
)
h
h
h
×
(
1
-
(
1
-
p
)
h
2
-
η
0
(
1
-
p
)
h
2
+
a
)
)
×
(
1
-
h
2
-
η
0
-
τ
(
1
-
q
)
(
h
2
+
a
)
)
-
1
,
μ
=
κ
1
exp
⁡
[
∫
0
T
‍
(
∂
P
∂
x
+
∂
Q
∂
y
)
d
t
]
=
κ
1
exp
⁡
[
∫
0
T
‍
(
1
-
2
x
-
a
(
a
+
x
)
2
y
+
δ
-
2
β
y
x
)
d
t
]
=
κ
1
exp
⁡
{
∫
0
T
‍
[
-
x
+
x
y
(
a
+
x
)
2
-
β
y
x
]
(
1
-
x
-
y
x
+
a
)
+
(
δ
-
β
y
x
)
H
H
H
H
H
H
H
h
-
x
+
x
y
(
a
+
x
)
2
-
β
y
x
]
d
t
}
=

κ
1
exp
⁡
{
∫
0
T
‍
d
x
x
+
∫
0
T
‍
d
y
y
-
∫
0
T
‍
(
x
-
x
y
(
a
+
x
)
2
+
β
y
x
)
d
t
}
=

(
1
-
p
)
(
1
-
q
)
(
1
-
(
1
-
p
)
h
2
-
(
η
0
/
(
(
1
-
p
)
h
2
+
a
)
)
)
1
-
h
2
-
(
(
η
0
-
τ
)
/
(
(
1
-
q
)
(
h
2
+
a
)
)
)
×

exp
⁡
{
∫
0
T
‍
d
x
x
+
∫
0
T
‍
d
y
y
-
∫
0
T
‍
(
x
-
x
y
(
a
+
x
)
2
+
β
y
x
)
d
t
}
=
η
0
-
τ
η
0
·
1
-
(
1
-
p
)
h
2
-
(
η
0
/
(
(
1
-
p
)
h
2
+
a
)
)
1
-
h
2
-
(
(
η
0
-
τ
)
/
(
(
1
-
q
)
(
h
2
+
a
)
)
)
×
exp
⁡
[
-
∫
0
T
‍
(
x
-
x
y
(
a
+
x
)
2
+
β
y
x
)
d
t
]
=
κ
exp
⁡
[
-
∫
0
T
‍
(
x
-
x
y
(
a
+
x
)
2
+
β
y
x
)
d
t
]
.
From Lemma 7, if |
μ
|
=
|
κ
exp
⁡
{
-
∫
0
T
‍
(
ξ
(
t
)
-
(
ξ
(
t
)
η
(
t
)
/
(
ξ
(
t
)
+
a
)
2
)
+
(
β
η
(
t
)
/
ξ
(
t
)
)
)
d
t
}
|
<
1, then the periodic solution of system (6) is orbitally asymptotically stable. This completes the proof.Remark 19.
If|
(
(
1
-
(
1
-
p
)
h
2
)
-
(
η
0
/
(
(
1
-
p
)
h
2
+
a
)
)
)
/
(
1
-
h
2
-
(
(
η
0
-
τ
)
/
(
(
1
-
q
)
(
h
2
+
a
)
)
)
)
|
≤
1 and β
≥
1, the periodic solution of system (6) with initial point C
0
(
(
1
-
p
)
h
2
,
η
0
) is orbitally asymptotically stable.


## 4.1. The Initial Point onN
1

LetC
0
(
h
1
,
y
C
0
) be an initial point of the orbit of the system (6); if y
C
0
<
y
~, point C
0 is below point P
(
h
1
,
y
~
), then C
0
∈
M
1 (M
1 is impulsive set), and the image point C
0
+ of C
0 must be above point P with n times impulses; therefore, we only need to discuss the cases of y
C
0
>
y
~.The orbitΠ
(
C
0
,
t
) starting from C
0
(
h
1
,
y
C
0
) hits the impulsive set M
1 at point C
1
(
h
1
,
y
C
1
), and then C
1 jumps to point C
11
(
h
1
,
y
C
11
). If y
C
11
<
y
~, C
11 continued to jump to point C
12
(
h
1
,
y
C
12
), and after n times it reaches the point C
1
n
(
h
1
,
y
C
1
n
), where y
C
11
=
y
C
1
+
λ, y
C
12
=
y
C
1
+
2
λ
,
…, y
C
1
n
=
y
C
1
+
n
λ, and y
~
<
y
C
1
n
<
y
~
+
λ. The situation of point C
1
n has three cases as follows.(a) Ify
C
1
n
=
y
C
0(see Figure 3(a)), C
1
n is coincident with C
0, then the curve C
0
C
1
C
11
⋯
C
1
n is closed orbit, and the system (6) exhibits a 1-periodic solution.The orbit starting from the pointC
0 on N
1.
(a)
(b)
(c)(b) Ify
C
1
n
>
y
C
0 (see Figure 3(b)), C
1
n is above C
0; in this time, the successor function of C
0 satisfies f
(
C
0
)
=
y
C
1
n
-
y
C
0
>
0. In the meantime, choose a point D
0
(
h
1
,
y
D
0
) on N
1 satisfying y
D
0
>
y
~
+
λ. The orbit Π
(
D
0
,
t
) starting from D
0 hits the impulsive set at point D
1
(
h
1
,
y
D
1
), and D
1 jumps sometimes to point D
1
m
(
h
1
,
y
D
1
m
), where y
~
<
y
D
1
m
<
y
~
+
λ. Thus, the successor function of D
0 satisfies f
(
D
0
)
=
y
D
1
m
-
y
D
0
<
0.According to Lemma6, the system (6) exhibits a periodic solution, and the initial point of the periodic solution is between C
0 and D
0.(c) Ify
C
1
n
<
y
C
0 (see Figure 3(c)), C
1
n is below C
0; in this time, the successor function of C
0 is f
(
C
0
)
=
y
C
1
n
-
y
C
0
<
0. On the other hand, choose a point D
0
(
h
1
,
y
D
0
) on N
1 satisfying y
~
<
y
D
0
<
y
~
+
ε (ε is a sufficiently small positive number). The orbit Π
(
D
0
,
t
) starting from D
0 hits the impulsive set M
1 at point D
1
(
h
1
,
y
D
1
), and D
1 jumps to point D
11
(
h
1
,
y
D
11
), where y
D
11
=
y
D
1
+
λ. As D
0 sufficiently closed to P, D
1 is sufficiently close to P, then y
D
11
=
y
D
1
+
λ
>
y
D
0. Thus, the successor function of D
0 satisfiesf
(
D
0
)
=
y
D
1
-
y
D
0
>
0.From Lemma6, system (6) exhibits a 1-periodic solution.To sum up the above discussed, we get the following.Theorem 13.
If the initial pointC
0
(
h
1
,
y
c
0
) of the orbit of the system (6) is on N
1 with y
C
0
>
y
~, then the system exhibits a 1-periodic solution.Next, we will discuss the stability of the above periodic solutions.Theorem 14.
Let(
ξ
(
t
)
,
η
(
t
)
) be the T-periodic solution of the system (6) with the initial point C
0
(
h
1
,
η
0
); the closed orbit corresponding to the periodic solution is the curve C
0
C
1
C
11
⋯
C
1
n, if
(20)
|
μ
|
=
|
κ
exp
⁡
{
-
∫
0
T
‍
(
ξ
(
t
)
-
ξ
(
t
)
η
(
t
)
(
ξ
(
t
)
+
a
)
2
+
β
η
(
t
)
ξ
(
t
)
)
d
t
}
|
<
1
,
where
(21)
κ
=
η
0
-
n
λ
η
0
∏
j
=
1
n
‍
1
-
h
1
-
(
1
/
(
h
1
+
a
)
)
(
η
0
-
(
n
-
j
)
λ
)
1
-
h
1
-
(
1
/
(
h
1
+
a
)
)
(
η
0
-
(
n
-
j
+
1
)
λ
)
,
then the periodic solution (
ξ
(
t
)
,
η
(
t
)
) is orbitally asymptotically stable.Proof.
Let the orbitΠ
(
C
0
,
t
) with the initial point C
0
(
h
1
,
η
0
) hit the impulsive set M
1 at C
1
(
ξ
(
T
)
,
η
(
T
)
), and then C
1 jumps to the point C
11
(
ξ
(
τ
1
)
,
η
(
τ
1
)
). The C
11 continued to jump to point C
12
(
ξ
(
τ
2
)
,
η
(
τ
2
)
), and at last, the image point C
12 reaches point C
1
n
(
ξ
(
τ
n
)
,
η
(
τ
n
)
) with n times pulses. Here, η
(
τ
1
)
=
η
(
T
)
+
λ
,
η
(
τ
2
)
=
η
(
T
)
+
2
λ
,
…
,
η
(
τ
n
)
=
η
(
T
)
+
n
λ. For the jth time impulse, obviously ξ
(
τ
j
+
)
=
ξ
(
τ
j
+
1
)
,
η
(
τ
j
+
)
=
η
(
τ
j
+
1
). For the system (6), let P
(
x
,
y
)
=
x
(
1
-
x
)
-
(
x
/
(
a
+
x
)
)
y, Q
(
x
,
y
)
=
y
(
δ
-
β
(
y
/
x
)
), α
(
x
,
y
)
=
0, β
(
x
,
y
)
=
λ, Φ
(
x
,
y
)
=
x
-
h
1, and ξ
(
T
)
=
h
1, η
(
T
)
=
η
0
-
n
λ; therefore we have
(22)
∂
P
∂
x
=
1
-
2
x
-
a
(
a
+
x
)
2
y
,
∂
Q
∂
y
=
δ
-
2
β
y
x
,
∂
α
∂
x
=
∂
α
∂
y
=
0
,
∂
β
∂
x
=
∂
β
∂
y
=
0
,
∂
Φ
∂
x
=
1
,
∂
Φ
∂
y
=
0
.
According to Lemma 7, we get
(23)
κ
j
=
(
(
∂
β
∂
y
·
∂
Φ
∂
x
-
∂
β
∂
x
·
∂
Φ
∂
y
+
∂
Φ
∂
x
)
P
+
h
h
h
+
(
∂
α
∂
x
·
∂
Φ
∂
y
-
∂
α
∂
y
·
∂
Φ
∂
x
+
∂
Φ
∂
y
)
Q
+
)
×
(
∂
Φ
∂
x
P
+
∂
Φ
∂
y
Q
)
-
1
=
P
(
ξ
(
τ
j
+
)
,
η
(
τ
j
+
)
)
P
(
ξ
(
τ
j
)
,
η
(
τ
j
)
)
=
P
(
h
1
,
η
(
T
)
+
j
λ
)
P
(
h
1
,
η
(
T
)
+
(
j
-
1
)
λ
)
=
1
-
h
1
-
(
1
/
(
h
1
+
a
)
)
(
η
0
-
(
n
-
j
)
λ
)
1
-
h
1
-
(
1
/
(
h
1
+
a
)
)
(
η
0
-
(
n
-
j
+
1
)
λ
)
=
1
-
λ
/
(
h
1
+
a
)
1
-
h
1
-
(
1
/
(
h
1
+
a
)
)
(
η
0
-
(
n
-
j
+
1
)
λ
)
,
μ
=
∏
j
=
1
n
‍
κ
j
exp
⁡
[
∫
0
T
‍
(
∂
P
∂
x
+
∂
Q
∂
y
)
d
t
]
=
∏
j
=
1
n
‍
κ
j
exp
⁡
[
∫
0
T
‍
(
1
-
2
x
-
a
(
a
+
x
)
2
y
+
δ
-
2
β
y
x
)
d
t
]
=
∏
j
=
1
n
‍
κ
j
exp
⁡
{
∫
0
T
‍
[
-
x
+
x
y
(
a
+
x
)
2
-
β
y
x
]
(
1
-
x
-
y
x
+
a
)
+
(
δ
-
β
y
x
)
h
h
h
h
h
h
h
h
h
h
h
h
h
h
-
x
+
x
y
(
a
+
x
)
2
-
β
y
x
]
d
t
}
=
∏
j
=
1
n
‍
κ
j
exp
⁡
{
∫
0
T
‍
d
x
x
+
∫
0
T
‍
d
y
y
h
h
h
h
h
h
h
h
h
h
-
∫
0
T
‍
(
x
-
x
y
(
a
+
x
)
2
+
β
y
x
)
d
t
}
=
∏
j
=
1
n
‍
κ
j
η
0
-
n
λ
η
0
exp
⁡
[
-
∫
0
T
‍
(
x
-
x
y
(
a
+
x
)
2
+
β
y
x
)
d
t
]
=
η
0
-
n
λ
η
0
×
∏
j
=
1
n
‍
1
-
h
1
-
(
1
/
(
h
1
+
a
)
)
(
η
0
-
(
n
-
j
)
λ
)
1
-
h
1
-
(
1
/
(
h
1
+
a
)
)
(
η
0
-
(
n
-
j
+
1
)
λ
)
×
exp
⁡
[
-
∫
0
T
‍
(
x
-
x
y
(
a
+
x
)
2
+
β
y
x
)
d
t
]
=
κ
exp
⁡
{
-
∫
0
T
‍
(
ξ
(
t
)
-
ξ
(
t
)
η
(
t
)
(
ξ
(
t
)
+
a
)
2
+
β
η
(
t
)
ξ
(
t
)
)
d
t
}
.
From Lemma 7, if |
μ
|
=
|
κ
exp
⁡
{
-
∫
0
T
‍
(
ξ
(
t
)
-
(
ξ
(
t
)
η
(
t
)
/
(
ξ
(
t
)
+
a
)
2
)
+
(
β
η
(
t
)
/
ξ
(
t
)
)
)
d
t
}
|
<
1, then the periodic solution of the system (6) is orbitally asymptotically stable. This completes the proof.Remark 15.
If1
-
h
1
-
(
η
0
/
(
h
1
+
a
)
)
>
0 and β
≥
1, then the periodic solution with initial point C
0
(
h
1
,
η
0
) (where η
0
>
y
~) is orbitally asymptotically stable.


## 4.2. The Initial Point onN
2

Based on existence and uniqueness theorem of differential equations, there exists a unique pointQ
0
(
(
1
-
p
)
h
2
,
y
Q
0
) on N
2 such that the orbit Π
(
Q
0
,
t
) starting from Q
0 is tangent to N
1 at point P
(
h
1
,
y
~
). Assume that C
0
(
(
1
-
p
)
h
2
,
y
C
0
) is the initial point of the orbit Π
(
C
0
,
t
) of system (6). Next, we will investigate the existence of periodic solution of the system with different positions of C
0 and Q
0. Three cases should be discussed.Case I (
y
C
0
=
y
Q
0
; see Figure4). The initial point C
0 is exactly Q
0.The initial pointC
0 on N
2 (Case  I: y
C
0
=
y
Q
0).
(a)
(b)
(c)The orbitΠ
(
Q
0
,
t
) starting from Q
0 is tangent to N
1 at the point P, and through N
2 hit M
2 at the point Q
1
(
h
2
,
y
Q
1
), and then Q
1 jumps to Q
1
+
(
x
Q
1
+
,
y
Q
1
+
) on N
2. According to (6), the following is obtained:
(24)
x
Q
1
+
=
(
1
-
p
)
h
2
,
y
Q
1
+
=
(
1
-
q
)
y
Q
1
+
τ
.
About the points Q
1
+ and Q
0, there are the following three positional relations.(a)
IfQ
1
+ coincides with Q
0
:
y
Q
1
+
=
y
Q
0 (see Figure 4(a)), then the curve Q
0
P
Q
1
Q
1
+ is closed orbit.
(b)
IfQ
1
+ is below Q
0
:
y
Q
1
+
<
y
Q
0 (see Figure 4(b)), then the successor function of Q
0 satisfies f
(
Q
0
)
=
y
Q
1
+
-
y
Q
0
<
0. In the meantime, take a point S
0
(
(
1
-
p
)
h
2
,
y
S
0
) on N
2 satisfying 0
<
y
S
0
<
ε (ε
>
0 small enough). The orbit Π
(
S
0
,
t
) starting from S
0 hits the impulsive M
2 at the point S
1
(
h
2
,
y
S
1
), and then S
1 jumps to the point S
1
+
(
x
S
1
+
,
y
S
1
+
), where x
S
1
+
=
(
1
-
p
)
h
2
,
y
S
1
+
=
(
1
-
q
)
y
S
1
+
τ. Obviously, the successor function of S
0 is f
(
S
0
)
=
y
S
1
+
-
y
S
0
>
0. From Lemma 6, the system (6) has an order one periodic solution, where the initial point of the periodic solution is between Q
0 and S
0.
(c)
IfQ
1
+ is above Q
0
:
y
Q
1
+
>
y
Q
0 (see Figure 4(c)), the system (6) does not have closed orbit in the area Ω
1
=
{
(
x
,
y
)
∣
h
1
<
x
<
h
2
} in this time.Based on the discussion above, we get the following.Theorem 16.
Assume that the orbitΠ
(
Q
0
,
t
) starting from Q
0
(
(
1
-
p
)
h
2
,
y
Q
0
) is tangent to N
1 at the point P
(
h
1
,
y
~
) and hits the impulsive set M
2 at the point Q
1
(
h
2
,
y
Q
1
), the image point of Q
1 is Q
1
+
(
x
Q
1
+
,
y
Q
1
+
) on N
2. If y
Q
1
+
≤
y
Q
0, the system (6) has 1-periodic solution in the area Ω
1
=
{
(
x
,
y
)
∣
h
1
<
x
<
h
2
}.Case II
(
y
C
0
<
y
Q
0
). The initial point C
0 is below Q
0.The isoclineL
1
:
y
=
(
1
-
x
)
(
x
+
a
) intersects with the phase set N
2 at the point Q
(
(
1
-
p
)
h
2
,
y
Q
), and Q is below Q
0. In this case, we discuss the existence of periodic solution of the system (6) with the example y
C
0
=
y
Q (see Figure 5).The initial pointC
0 on N
2 (Case  II: y
C
0
<
y
Q
0).
(a)
(b)
(c)
(d)
(e)The orbitΠ
(
Q
,
t
) starting from Q moves to the point C
1
(
h
2
,
y
C
1
) on the impulsive set M
2, and C
1 jumps onto C
1
+
(
x
C
1
+
,
y
C
1
+
) on the image set N
2, and then
(25)
y
C
0
=
y
Q
=
(
1
-
(
1
-
p
)
h
2
)
(
a
+
(
1
-
p
)
h
2
)
,
x
C
1
+
=
(
1
-
p
)
h
2
,
y
C
1
+
=
(
1
-
q
)
y
C
1
+
τ
.
(a) Ify
C
1
+
=
y
Q
:
C
1
+ coincide with Q (see Figure 5(a)), the curve Q
C
1
C
1
+ is the closed orbit of the system (6).
(b) Ify
Q
<
y
C
1
+
≤
y
Q
0
:
C
1
+ is between Q and Q
0 (see Figure 5(b)), in this case the successor function of Q is f
(
Q
)
=
y
C
1
+
-
y
Q
>
0. On the other hand, consider the orbit Π
(
Q
0
,
t
) starting from Q
0
,
Π
(
Q
0
,
t
) hits the impulsive set M
2 at Q
1
(
h
2
,
y
Q
1
), and Q
1 jumps onto Q
1
+
(
x
Q
1
+
,
y
Q
1
+
) on N
2. Based on the existence and uniqueness theorem of differential equations, Q
1 must be below C
1, and Q
1
+ must be below C
1
+. Thus y
Q
1
+
<
y
C
1
+
≤
y
Q
0, and the successor function of Q
0 is f
(
Q
0
)
=
y
Q
1
+
-
y
Q
0
<
0. From Lemma 6, the system (6) has 1-periodic solution, and the initial point of the periodic solution is between Q
0 and Q.
Ify
C
1
+
>
y
Q
0
:
C
1
+ is above Q
0, in this case the successor function of Q is f
(
Q
)
=
y
C
1
+
-
y
Q
>
0. There are two different cases (case (c) and case (d)).
(c) Ify
C
1
+
>
y
Q
0 and y
Q
1
+
≤
y
Q
0
:
Q
1
+ is below Q
0 (see Figure 5(c)), f
(
Q
0
)
=
y
Q
1
+
-
y
Q
0
≤
0, the system (6) has closed orbit.
(d) Ify
C
1
+
>
y
Q
0 and y
Q
1
+
>
y
Q
0
:
Q
1
+ is above Q
0 (see Figure 5(d)), the system (6) does not have closed orbit in the area Ω
1
=
{
(
x
,
y
)
∣
h
1
<
x
<
h
2
}.
(e) Ify
C
1
+
<
y
Q
:
C
1
+ is below Q (see Figure 5(e)), in this case the successor function of Q is f
(
Q
)
=
y
C
1
+
-
y
Q
<
0. On the other hand, take a point D
0
(
(
1
-
p
)
h
2
,
y
D
0
) from N
2 satisfying that y
D
0 is sufficiently small number, which is 0
<
y
D
0
<
ε (ε
>
0 small enough). The orbit Π
(
D
0
,
t
) starting from D
0 moves to D
1
(
h
2
,
y
D
1
) on the impulsive set M
2, and D
1 jumps onto D
1
+
(
x
D
1
+
,
y
D
1
+
) on N
2, where x
D
1
+
=
(
1
-
p
)
h
2, y
D
1
+
=
(
1
-
q
)
y
D
1
+
τ. Obviously, the successor function of D
0 is f
(
D
0
)
=
y
D
1
+
-
y
D
0
>
0. From Lemma 6, the system (6) has closed orbit.Based on the discussion above, we get the following.Theorem 17.
Assume that the orbitΠ
(
Q
,
t
) starting from Q
(
(
1
-
p
)
h
2
,
y
Q
) hits the impulsive set M
2 at C
1
(
h
2
,
y
C
1
), and C
1 jumps onto C
1
+
(
x
C
1
+
,
y
C
1
+
) on N
2. The orbit Π
(
Q
0
,
t
) starting from the point Q
0
(
(
1
-
p
)
h
2
,
y
Q
0
) is tangent to N
1 at P
(
h
1
,
y
~
), and hitting the impulsive set M
2 at Q
1
(
h
2
,
y
Q
1
), the image point of Q
1 is Q
1
+
(
x
Q
1
+
,
y
Q
1
+
). Then one has the following. (1)
Ify
C
1
+
≤
y
Q
0, the system (6) has 1-periodic solution.
(2)
Ify
C
1
+
>
y
Q
0 and y
Q
1
+
≤
y
Q
0, the system (6) has 1-periodic solution.Case III (
y
C
0
>
y
Q
0
; see Figure6). The initial point C
0 is above Q
0.Figure 6
The initial pointC
0 on N
2 (Case  III: y
C
0
>
y
Q
0).In this case, the orbitΠ
(
C
0
,
t
) starting from C
0 goes through the isocline L
1 from the left of the line x
=
h
1, hitting the impulsive set M
1 at C
1
(
h
1
,
y
C
1
). The same conclusion can be made as in Section 4.1.Next, we discuss the stability of the periodic solution with the initial point onN
2.Theorem 18.
Assume that(
ξ
(
t
)
,
η
(
t
)
) is the T-periodic solution of the system (6) with initial point C
0
(
(
1
-
p
)
h
2
,
η
0
); and if
(26)
|
μ
|
=
|
κ
exp
⁡
{
-
∫
0
T
‍
(
ξ
(
t
)
-
ξ
(
t
)
η
(
t
)
(
ξ
(
t
)
+
a
)
2
+
β
η
(
t
)
ξ
(
t
)
)
d
t
}
|
<
1
,
where
(27)
κ
=
η
0
-
τ
η
0
·
1
-
(
1
-
p
)
h
2
-
(
η
0
/
(
(
1
-
p
)
h
2
+
a
)
)
1
-
h
2
-
(
(
η
0
-
τ
)
/
(
(
1
-
q
)
(
h
2
+
a
)
)
)
,
the periodic solution (
ξ
(
t
)
,
η
(
t
)
) is orbitally asymptotically stable.Proof.
Assume that the periodic orbitΠ
(
C
0
,
t
) starting from the point C
0
(
(
1
-
p
)
h
2
,
η
0
) moves to the point C
1
(
ξ
(
T
)
,
η
(
T
)
) on impulsive set M
2, and C
1 jumps onto the point C
1
+
(
ξ
(
T
+
)
,
η
(
T
+
)
) on N
2. Therefore, Π
(
C
0
,
T
)
=
C
1, C
1
+
=
φ
2
(
C
1
)
=
C
0, ξ
(
T
+
)
=
(
1
-
p
)
ξ
(
T
), η
(
T
+
)
=
(
1
-
q
)
η
(
T
)
+
τ.
Compared with the system (6), we get
(28)
P
(
x
,
y
)
=
x
(
1
-
x
)
-
x
a
+
x
y
,
Q
(
x
,
y
)
=
y
(
δ
-
β
y
x
)
,
α
(
x
,
y
)
=
-
p
x
,
β
(
x
,
y
)
=
-
q
y
+
τ
,
Φ
(
x
,
y
)
=
x
-
h
2
,
ξ
(
T
)
=
h
2
,
η
(
T
)
=
η
0
-
τ
1
-
q
,
∂
P
∂
x
=
1
-
2
x
-
a
(
a
+
x
)
2
y
,
∂
Q
∂
y
=
δ
-
2
β
y
x
,
∂
α
∂
x
=
-
p
,
∂
α
∂
y
=
0
,
∂
β
∂
x
=
0
,
∂
β
∂
y
=
-
q
,
∂
Φ
∂
x
=
1
,
∂
Φ
∂
y
=
0
.
Thus,
(29)
κ
1
=
(
(
∂
β
∂
y
·
∂
Φ
∂
x
-
∂
β
∂
x
·
∂
Φ
∂
y
+
∂
Φ
∂
x
)
P
+
hhhh
+
(
∂
α
∂
x
·
∂
Φ
∂
y
-
∂
α
∂
y
·
∂
Φ
∂
x
+
∂
Φ
∂
y
)
Q
+
)
×
(
∂
Φ
∂
x
P
+
∂
Φ
∂
y
Q
)
-
1
=
(
1
-
q
)
P
(
ξ
(
T
+
)
,
η
(
T
+
)
)
P
(
ξ
(
T
)
,
η
(
T
)
)
=
(
η
0
(
1
-
p
)
h
2
+
a
(
1
-
p
)
(
1
-
q
)
h
h
h
×
(
1
-
(
1
-
p
)
h
2
-
η
0
(
1
-
p
)
h
2
+
a
)
)
×
(
1
-
h
2
-
η
0
-
τ
(
1
-
q
)
(
h
2
+
a
)
)
-
1
,
μ
=
κ
1
exp
⁡
[
∫
0
T
‍
(
∂
P
∂
x
+
∂
Q
∂
y
)
d
t
]
=
κ
1
exp
⁡
[
∫
0
T
‍
(
1
-
2
x
-
a
(
a
+
x
)
2
y
+
δ
-
2
β
y
x
)
d
t
]
=
κ
1
exp
⁡
{
∫
0
T
‍
[
-
x
+
x
y
(
a
+
x
)
2
-
β
y
x
]
(
1
-
x
-
y
x
+
a
)
+
(
δ
-
β
y
x
)
H
H
H
H
H
H
H
h
-
x
+
x
y
(
a
+
x
)
2
-
β
y
x
]
d
t
}
=

κ
1
exp
⁡
{
∫
0
T
‍
d
x
x
+
∫
0
T
‍
d
y
y
-
∫
0
T
‍
(
x
-
x
y
(
a
+
x
)
2
+
β
y
x
)
d
t
}
=

(
1
-
p
)
(
1
-
q
)
(
1
-
(
1
-
p
)
h
2
-
(
η
0
/
(
(
1
-
p
)
h
2
+
a
)
)
)
1
-
h
2
-
(
(
η
0
-
τ
)
/
(
(
1
-
q
)
(
h
2
+
a
)
)
)
×

exp
⁡
{
∫
0
T
‍
d
x
x
+
∫
0
T
‍
d
y
y
-
∫
0
T
‍
(
x
-
x
y
(
a
+
x
)
2
+
β
y
x
)
d
t
}
=
η
0
-
τ
η
0
·
1
-
(
1
-
p
)
h
2
-
(
η
0
/
(
(
1
-
p
)
h
2
+
a
)
)
1
-
h
2
-
(
(
η
0
-
τ
)
/
(
(
1
-
q
)
(
h
2
+
a
)
)
)
×
exp
⁡
[
-
∫
0
T
‍
(
x
-
x
y
(
a
+
x
)
2
+
β
y
x
)
d
t
]
=
κ
exp
⁡
[
-
∫
0
T
‍
(
x
-
x
y
(
a
+
x
)
2
+
β
y
x
)
d
t
]
.
From Lemma 7, if |
μ
|
=
|
κ
exp
⁡
{
-
∫
0
T
‍
(
ξ
(
t
)
-
(
ξ
(
t
)
η
(
t
)
/
(
ξ
(
t
)
+
a
)
2
)
+
(
β
η
(
t
)
/
ξ
(
t
)
)
)
d
t
}
|
<
1, then the periodic solution of system (6) is orbitally asymptotically stable. This completes the proof.Remark 19.
If|
(
(
1
-
(
1
-
p
)
h
2
)
-
(
η
0
/
(
(
1
-
p
)
h
2
+
a
)
)
)
/
(
1
-
h
2
-
(
(
η
0
-
τ
)
/
(
(
1
-
q
)
(
h
2
+
a
)
)
)
)
|
≤
1 and β
≥
1, the periodic solution of system (6) with initial point C
0
(
(
1
-
p
)
h
2
,
η
0
) is orbitally asymptotically stable.


## 5. Example and Numerical Simulation

In this part, we use numerical simulation to confirm the conclusion obtained above. Leta
=
0.05, δ
=
0.5, β
=
0.7, h
1
=
0.24, h
2
=
0.35, λ
=
0.18, p
=
0.2, q
=
0.2. By calculation, we obtain y
~
=
0.22, P
(
0.24,0.22
), Q
(
0.28,0.2376
) and R
(
0.35,0.26
). Then, we have an example as follows:
(30)
d
x
(
t
)
d
t
=
x
(
t
)
(
1
-
x
(
t
)
)
-
x
(
t
)
0.05
+
x
(
t
)
y
(
t
)
,
d
y
(
t
)
d
t
=
y
(
t
)
(
0.5
-
0.7
y
(
t
)
x
(
t
)
)
,
x
≠
0.24,0.35
or
x
=
0.24
,
y
>
0.22
,
Δ
x
(
t
)
=
0
,
Δ
y
(
t
)
=
0.18
,
x
=
0.24
,
y
≤
0.22
,
Δ
x
(
t
)
=
-
0.2
x
(
t
)
,
Δ
y
(
t
)
=
-
0.2
y
(
t
)
+
τ
,
x
=
0.35
,
x
(
0
)
>
0
,
y
(
0
)
>
0
.Case 1.
Letτ
=
0.15, and the initial point is (
0.24,0.3
). From Figure 7 corresponding to Figure 3, the system exhibits a 1-periodic solution.The periodic solution corresponding to Figure3.
(a)
Phase diagram of system  (30)
(b)
Time series of system  (30)
(c)
Time series of system  (30)Case 2.
Letτ
=
0.1, then (
0.28,0.3056
) is the initial point. From Figure 8, which corresponds to Figure 4(b), the system exhibits a 1-periodic solution in the area Ω
1
=
{
(
x
,
y
)
∣
h
1
<
x
<
h
2
}.There exists a 1-periodic solution in the areaΩ
1
=
{
(
x
,
y
)
∣
h
1
<
x
<
h
2
} corresponding to Figure 4(b).
(a)
Phase diagram of system  (30)
(b)
Time series of system  (30)
(c)
Time series of system  (30)Case 3.
Letτ
=
0.15, we get the initial point (
0.28,0.3056
). It is easy to find that the system has no 1-periodic solution in the area Ω
1
=
{
(
x
,
y
)
∣
h
1
<
x
<
h
2
} from Figure 9 which corresponds to Figure 4(c).1-periodic solution in the areaΩ
1
=
{
(
x
,
y
)
∣
h
1
<
x
<
h
2
} corresponding to Figure 4(c) does not exist.
(a)
Phase diagram of system  (30)
(b)
Time series of system  (30)
(c)
Time series of system  (30)Case 4.
Letτ
=
0.1, the initial point (
0.28,0.2376
) is obtained. From Figure 10 corresponding to Figure 5(b), in the area Ω
1
=
{
(
x
,
y
)
∣
h
1
<
x
<
h
2
}, the system exhibits a 1-periodic solution.There exists a 1-periodic solution in the areaΩ
1
=
{
(
x
,
y
)
∣
h
1
<
x
<
h
2
} corresponding to Figure 5(b).
(a)
Phase diagram of system  (30)
(b)
Time series of system  (30)
(c)
Time series of system  (30)Case 5.
Letτ
=
0.15, we can easily get the initial point (
0.28,0.2376
). From Figure 11 corresponding to Figure 5(d), clearly, there is no 1-periodic solution in the area Ω
1
=
{
(
x
,
y
)
∣
h
1
<
x
<
h
2
}.1-periodic solution in the areaΩ
1
=
{
(
x
,
y
)
∣
h
1
<
x
<
h
2
} corresponding to Figure 5(d) does not exist.
(a)
Phase diagram of system  (30)
(b)
Time series of system  (30)
(c)
Time series of system  (30)Case 6.
Letτ
=
0.15, and the initial point is (
0.28,0.45
). Obviously, we can find a 1-periodic solution from Figure 12 which corresponds to Figure 6.The periodic solution corresponding to Figure6.
(a)
Phase diagram of system  (30)
(b)
Time series of system  (30)
(c)
Time series of system  (30)All the simulations above show agreement with the results in Section4.


## 6. Conclusion

This paper establishes a class of integrated pest management model based on state impulse control. In the initial stage of the occurrence of crop pests, that is, the pest density satisfiesx
(
t
)
≤
h
1, we use environment protection measures to control pests, such as releasing natural enemies. Once the pest density reaches a higher level x
(
t
)
=
h
2, we will adapt a combination of spraying insecticide and releasing natural enemies to control pests. With a short time to finish spraying insecticide and releasing natural enemies which bring out a sharp change in the number of pests and natural enemies, the state impulsive differential system (6) is obtained. Firstly, let the control parameters p, q, λ, τ be zero, we get Holling-Tanner ecosystem without impulsive effects. By constructing Dulac function, we discussed the stability of the positive equilibrium point E
(
x
*
,
y
*
), and the globally asymptotically stable conditions are given for focal points and nodal point, respectively. If the control parameters p, q, λ, τ are larger than zero, the system (6) is semicontinuous pulse dynamic system. The existence, uniqueness, and stability of the periodic solutions are the research difficulties, and we need to consider all the pulse conditions (the value of h
1, h
2) and pulse function and its corresponding qualitative properties of the continuous dynamic system. By introducing the successor function, using impulsive differential geometry theory, we have discussed the existence of periodic solutions of the system (6) with a focus. According to the theory of impulsive differential multiplier Analogue of Poincare Criterion, the conditions of periodic solution with orbit asymptotically stable are given. Since (6) is a two-dimensional dynamical system, geometric method is intuitive and effective. How to study high-dimensional ecological dynamic systems by the geometric theory needs to be resolved in the future.



---

*Source: 963072-2014-03-19.xml*