![MathJax Logo](/templates/jsp/_style2/_tandf/pb2/images/math-jax.gif)
Abstract
In this paper, a general predator–prey model with state-dependent impulse is considered. Based on the geometric analysis and Poincaré map or successor function, we construct three typical types of Bendixson domains to obtain some sufficient conditions for the existence of order-1 periodic solutions. At the same time, the existing domains are discussed with respect to the system parameters. Moreover, the Analogue of Poincaré Criterion is used to obtain the asymptotic stability of the periodic solutions. Finally, to illustrate the results, an example is presented and some numerical simulations are carried out.
1. Introduction
Impulsive control is widely utilized in biology [Citation6,Citation7,Citation25]. For instance, in the treatment of HIV-1 infection, the combination of vaccination and stimulation delays viral rebound following ART (Antiretroviral therapy) discontinuation [Citation2,Citation11,Citation12]. Compared with the process of disease, the impact of taking drugs, vaccination and stimulation is short enough to be assumed that the therapy leads to an impulsive effect. Usually, from the economy point of view, the corresponding treatment is proposed to receive when the infected system reaches a critical state, on which the amount decreases to 350 or 500 mm
. Similarly, in the controlling of insects and other animals, to minimize damage to non-target organisms, and to preserve the environment, the integrated pest management (IPM) was introduced [Citation3,Citation17], which is a strategy that reduce pests to tolerable levels by releasing natural enemies, spraying pesticides, isolating or harvesting the pests when they reach a threshold. Likewise, in fishery industry, for the sake of ecological balance and maximum economic benefit, people usually take a sustainable strategy by harvesting and releasing according to the densities of certain species. For example, sometimes we need suppress the predator-fish by harvesting, and to stimulate the prey-fish by releasing, and other time it's just on the opposite purpose [Citation9]. The discontinuity of human intervention leads to population changing very rapidly, which can be taken into account impulsive differential equations.
The impulsive differential equations have even been intensively studied in recent two decades, especially for state-dependent impulsive equations, which means the impulsive control takes place when a certain state is met [Citation4,Citation5,Citation8,Citation10–24,Citation26]. However, because of the discontinuity of the system, the issue becomes more complicated but not less important than a continuous one. So it is still challenging to probe the dynamical properties of such systems.
In [Citation21], Wang et al. consider an impulsive Kolmogorov predator–prey model with non-selective harvesting given by
(1)
(1) where
and
denote the densities of the prey and predator species at time t. Suppose one expects to perform an economical harvesting on the prey and predator (e.g. two kinds of food fish). When the density of prey increases to h, one harvests the prey and the predator at a rate p and q respectively. Because the reduction of prey may prevent predators from getting enough growth, and the harvesting on predator may also make it be on the verge of extinction, therefore to maintain the ecological balance among two populations and to obtain a sustainable economic benefit, one can release the predator on a scare of τ.
The authors incorporated two discriminants and
to discuss the existence, non-existence and multiplicity as well as the stability of order-1 periodic solutions. Motivated by the work of [Citation20,Citation21], we consider the following state-dependent impulsive predator–prey population ecological model:
(2)
(2) where
, h and τ are in the same meaning as those in (Equation1
(1)
(1) ); n and m represent the rates at which the prey and the predator are harvested, respectively.
Assume that h is positive constant, and ,
. The initial value conditions satisfy
(3)
(3) Moreover, the following assumptions hold:
,
,
for x>0.
There exist two positive numbers w<k such that
and
, and
.
In this paper, based on the theory of Bendixson domains established in [Citation20,Citation21], some results on the existence and stability of order-1 periodic solutions for system (Equation2(2)
(2) ) are established. Since the isolines of system (Equation2
(2)
(2) ) are different from that of system (Equation1
(1)
(1) ), the discriminants
and
are not always valid for system (Equation2
(2)
(2) ). We classify the structure of Bendixson domains into three categories, and combine them with Poincaré maps and successor functions, which help us to obtain the existence and location of impulsive order-1 periodic solution(s). Further, by using the Analogue of Poincaré Criterion, the stability of the order-1 periodic solution for system (Equation2
(2)
(2) ) is considered. In the last, we focus on the numerical simulations for the theory results.
2. Preliminaries
Let the impulse function be I, then . For the sake of simplicity, it is also denoted by
. Further, we denote the system without impulse that corresponding to (Equation2
(2)
(2) ) by
(4)
(4) Let impulse set
and image set
, where
. Obviously
. Let
be the positive trajectory of the system (Equation4
(4)
(4) ) starting from the point Q and let
be a negative trajectory reaching the point Q. For any
, if
intersects with M for the first time at a point
, and is consecutively mapped to
, then we define Poincaré map
and successor function F as follows:
(5)
(5) where Q also represents its coordinate
. For convenience, if a point A is above point B, that is
, then we write A>B. Obviously,
,
The geometric interpretation of Poincaré map and successor function F is shown in the Figure .
Definition 2.1
If , that is
, then the ring formed by trajectory
and line segment
is called the order-1 periodic orbit.
It is known from the properties of U and V, for system (Equation4(4)
(4) ), that the existence and uniqueness of the solutions hold true, and the solutions are continuously dependent on and differentiable with respect to the initial value. It is also easy to verify that the solution of (Equation4
(4)
(4) ) with a positive initial value must be positive.
Lemma 2.1
Under the Assumption –
, the system (Equation4
(4)
(4) ) has two unstable saddle points
,
and an asymptotically stable node or focus
, where
is the solution of the equation
. There is no positive periodic solution in the first quadrant for system (Equation4
(4)
(4) ).
Proof.
Obviously, the point is an equilibrium. Since
, the point
is a boundary equilibrium of the system (Equation4
(4)
(4) ). In addition, according to assumption
, the equation
determine an implicit function, named
. From the differential rule for implicit function, we have
which means that the function
is decreasing with x, and
because of
. Since
and w<k, so the curve of
must intersect with the line x = w at the point
, where
is the solution of the equation
. Hence
is a positive equilibrium of the system (Equation4
(4)
(4) ).
The Jacobian matrix along the system (Equation4(4)
(4) ) is
Calculating the eigenvalues of the Jacobian matrix at
,
and
, respectively, it is known that the equilibria
,
are unstable saddle points and
is an asymptotically stable node or focus.
Taking a Dulac function , we have
Therefore, the unique positive equilibrium point
is globally asymptotically stable and there is no positive periodic solution in the first quadrant.
When , we called the function P,
, F are well defined. Following the discussion in [Citation21] and [Citation20], we have
Lemma 2.2
If P,, F are well defined on a non-tangent interval
(except the end point), then P,
, F are continuous on
.
Further similar analysis to Lemma 2.2 in [Citation21] shows that the vector field of the system (Equation4(4)
(4) ) is counterclockwise.
Let the isoline intersects with M and N at points
and
respectively. Then
(6)
(6) Because of
and
, we have
. The possible location of impulse set, image set, isolines and related characterized points is shown in Figure .
Figure 2. The possible location of impulse section, image section, isoclines and some characterized points: (a); (b)
.
![Figure 2. The possible location of impulse section, image section, isoclines and some characterized points: (a)h_<h<x∗=w; (b) h_<x∗<h.](/cms/asset/d564c07d-9e70-4d7e-a291-ebe5c590c363/tjbd_a_1937721_f0002_oc.jpg)
Definition 2.2
For system (Equation2(2)
(2) ), suppose a domain D is composed of impulse set M, image set N, trajectory
and
, and satisfies
there is no singularity in D;
the trajectory
intersects with N and M at A and
successively; the trajectory
intersects with N, M at B and
successively;
line segments AB and
have no tangent points with the trajectory of the system (Equation4
(4)
(4) ) except for the end points.
If , then D is called a parallel rectangular domain (see Figure a).
Figure 3. The illustration for three categories of Bendexion domain: (a) parallel domain; (b) sub-parallel domain; (c) semi-ring domain.
![Figure 3. The illustration for three categories of Bendexion domain: (a) parallel domain; (b) sub-parallel domain; (c) semi-ring domain.](/cms/asset/38d1b128-29e5-42ea-982f-0499d2ccd7d8/tjbd_a_1937721_f0003_oc.jpg)
If , and
is tangent to N at B, then D is called a sub-parallel rectangular domain (see Figure b).
If ,
intersects with N at A and
successively, and
intersects with N at B and
successively, and
is tangent to M at
, then D is called a semi-circle domain (see Figure c).
Lemma 2.3
[Citation20,Citation21]
For parallel and sub-parallel domains, if F is defined on the non-tangent interval and such that
, then there is an order-1 periodic solution starting from
. For semi-ring domains, if
or
, then there is an order-1 periodic solution starting from AB. If
or
, then there is an order-1 periodic solution starting from
. If
or
, then there is no order-1 periodic solution in D (see Figure ).
3. Main results
Let ,
.
Lemma 3.1
For any , if
, then
; if
, then
.
Proof.
Obviously, . Assume that
. On the one hand, we have
. On the other hand,
. Thus
, which means
. Similarly, when
, we have
, that is
.
Further, if is well defined, then
(7)
(7) Next, we discuss the existence and location of order-1 periodic solutions of the system (Equation2
(2)
(2) ) under the case
and
, respectively.
3.1. ![](//:0)
![](//:0)
In this case, for any ,
is well defined since
is globally asymptotically stable.
Theorem 3.1
Assume –
hold and
. Then there must exist a positive order-1 periodic solution for system (Equation2
(2)
(2) ).
Proof.
To specify the structure and location of the order-1 periodic solution and probe the effect of parameter m and τ, we discuss case by case.
Case 1. .
. The domain formed by trajectory
and
, segments
and
is sub-parallel. It follows from the structure of the trajectories of system (Equation2
(2)
(2) ) and the monotonic increasing of impulse function I that
. Thus
. According to Lemma 2.3, there is an order-1 periodic solution in this sub-parallel domain.
. Because
is well defined, the domain composed of trajectories
and
, segments
and
is a parallel one. It follows from (Equation7
(7)
(7) ) that
, which implies that there is a positive order-1 periodic solution starting from the segment
.
Case 2. .
Take the trajectory as a graph of a function, then is decreasing with respect to x for
. Therefore
. It follows from Lemma 3.1 that
no matter
or
, and hence
. Moreover
, which leads to
. So in the parallel domain composed of trajectories
,
and segments
and
, there is a positive order-1 periodic solution initiating from the segment
.
Theorem 3.2
Suppose that . Then there is a semi-trivial periodic solution
, where
.
Proof.
When ,
is decreasing for any Q<T, and it follows that
. If
, then
for Q<T. Moreover, similar to the discussion of Lemma 2.2 in [Citation21], we have
(8)
(8) From the assumption
–
, it gives
for
. Thus for the successor function F, both
and
are satisfied for Q<T. Since 0<m<1 and the solutions of system (Equation4
(4)
(4) ) are positive with positive initial value. Henceforth,
holds true because of the continuity of F.
Remark 3.1
In fact, even since , the inequality
is always true under the case that
is well defined and
. In view of the uniqueness of the solution, the trajectory
intersects with the isoline on the right side of M, therefore
holds, which means
.
When m = 0, the result is similar to the case .
Corollary 3.1
Suppose m = 0. Then the following statements are true. (i) If , then there is an order-1 periodic solution initiating above T. (ii) If
, then the order-1 periodic solution is starting from
, which is below T.
3.2. ![](//:0)
![](//:0)
In this case, may be not well defined.
If is well defined,
is defined for any
. Similar to the proof of Theorem 3.1, there must exist an order-1 periodic solution for system (Equation2
(2)
(2) ).
If is not well defined (it is only possible when
), we only consider the case that
and there are at least two intersection points, because the impulse control is hardly valid when there is unique or no intersection[Citation20,Citation21]. Let
From Lemma 2.3, we have the following theorems.
Theorem 3.3
Suppose that ,
and
hold. Then there is an order-1 periodic solution starting from segment
, which lies in the semi-ring domain.
Proof.
Obviously, under the conditions given by the theorem, holds true. It follows from Lemma 2.3 that there is an order-1 periodic solution in the semi-ring domain, which is composed of trajectories
,
and segments
and
.
Theorem 3.4
Suppose that and
hold. Then the system (Equation2
(2)
(2) ) has no positive periodic solution.
Proof.
Since , we have
. From Lemma 3.1, no matter
or
, the inequality
can be obtained from
and
. The trajectories starting from the points below
on N are all mapped onto segment
under the effect of impulse, so the trajectory initiating on
will no longer reach M, which are determined by system (Equation4
(4)
(4) ). From Lemma 2.3, there is no positive periodic solution for system (Equation4
(4)
(4) ), so does for system (Equation2
(2)
(2) ).
Theorem 3.5
Suppose holds. If
, then there is an order-1 periodic solution which departures from segment
. Specifically, if
, then there is an order-1 periodic solution initiating from segment
.
Proof.
Equation (Equation7(7)
(7) ) gives
. Since
, that is
, and hence
holds. It is known from Lemma 2.3 that there is an order-1 periodic solution initiating from the segment
. Specifically, if
, then
. Therefore there is an positive order-1 periodic solution initiating from segment
.
Remark 3.2
Theorems 3.2–3.4 present the sufficient conditions for the existence or non-existence of order-1 periodic solutions to system (Equation2(2)
(2) ). We have not established the existence conditions of order-1 periodic solution when
and
, the simulations show that there maybe exist or not exist a positive order-1 periodic solution.
Remark 3.3
When is well defined, the results are the same as that under the case
.
Theorem 3.6
Suppose and
be not well defined. Then there is a semi-trivial order-1 periodic solution under the case
, and when
, there maybe exist an positive order-1 periodic solution.
Proof.
Since , we have
is well defined for
. Further,
holds true, which means there exists a semi-trivial order-1 periodic solution. When
, from Equation (Equation8
(8)
(8) ), the inequality
does not necessarily hold for
. Therefore, when
, there maybe exist an positive order-1 periodic solution besides the semi-trivial periodic solution.
Corollary 3.2
Suppose m = 0 and be well defined. Then the following statements are true. (i) If
, then there is an order-1 periodic solution initiating above T. (ii) If
, the order-1 periodic solution is initiating from
, which is below T.
Corollary 3.3
Suppose that m = 0 and be not well defined. (i) If
and
, there is an order -1 periodic solution starting above
. (ii) If
, there is an order -1 periodic solution starting from
.
Remark 3.4
When ,
maybe or not well defined. If
is well defined, for any
,
is well defined. If
is not well defined, for any
,
is well defined. When
, the trajectory
is decreasing with x on
, however, it is not necessarily decreasing on
under the case
. Hence, it is still possible that there is a point Q such that
, that is,
. If
is well defined and
, then it is still possible that there is a positive order-1 periodic solution starting from the segment QT. If
is not well defined and
, that is,
, then there may exist a positive order-1 periodic solution which initiating from the segment
.
Next, we summarize the results in two tables as follows, which are listed according to whether is well defined, where the symbol ‘par.’ and ‘semi-tri.’ are the abbreviations of ‘parallel’ and ‘semi-trivial’. The check mark ‘
’ indicates the claim is true(see Table and Table ).
Table 1. The existence and location of order-1 periodic solution when is well defined.
Table 2. The existence and location of order-1 periodic solution when is not well defined.
4. Stability of the order-1 positive periodic solution
Consider the general plane state-dependent impulsive differential system
(9)
(9) where
,
is continuously differentiable,
is a smooth function and
.
Lemma 4.1
Analogue of Poincaré Criterion [Citation1,Citation8,Citation15].
If , then the order-1 periodic solution of system (9) is orbitally asymptotically stable, where
(10)
(10)
(11)
(11) and
are calculated at the point
,
and
are calculated at the point
, and
is the time of the jth jump.
Theorem 4.1
Let be an order-1 ω-periodic solution of the system (Equation2
(2)
(2) ). If
, then
is orbitally asymptotically stable, where
(12)
(12)
Proof.
Suppose intersects with M, N at the point
and
, respectively, where
. Corresponding to system (Equation9
(9)
(9) ), we have
,
,
, and
and
. Therefore
It follows that
(13)
(13) Moreover, we have
(14)
(14) and
(15)
(15) Further, the last items of (Equation14
(14)
(14) ) give
(16)
(16) Hence, it follows from (Equation14
(14)
(14) )–(Equation16
(16)
(16) ) that
(17)
(17) Therefore, according to (Equation12
(12)
(12) ), (Equation13
(13)
(13) ), (Equation14
(14)
(14) ) and (Equation17
(17)
(17) ), we have
According to Lemma 4.1, the order-1 periodic solution
is orbitally asymptotically stable.
Remark 4.1
If and
, then there exists an order-1 periodic solution
, which must be orbitally asymptotically stable, in view of
and V <0, Y >0 for
, and hence
.
5. Example and simulations
Consider the state-dependent impulsive differential system
(18)
(18) Obviously, corresponding to system (Equation2
(2)
(2) ), we have
From
, it follows that
. It is not difficult to figure out
If k>h and
, then the assumptions
and
hold.
In the following, some simulations are carried out by Matlab software, which are somewhat corresponding to Table and Table .
5.1. ![](//:0)
is well defined
Taking the parameters r = 1, k = 10, , a = 0.188, n = 0.3 and h = 5.5, the computation shows that
,
,
and
, which means
.
Let and m = 0. Then
and
are satisfied. As can be seen from Figure (a) that there is an order-1 periodic solution in a sub-parallel domain. Let
and m = 0.1. Then we have
and
. There is an order-1 periodic solution initiating below T, which is in a parallel domain (Figure b). If
and m = 0.1, then the condition
holds true, it can be seen from Figure (c) that there exists an order-1 periodic solution in a parallel domain. Let
and m = 0.1. It is known from Figure (d) that there is a semi-trivial periodic solution in the system.
Figure 4. The order-1 periodic solution of system (Equation18(18)
(18) ) under the case
, where r = 1, k = 10,
, a = 0.188, n = 0.3,h = 5.5. (a) m = 0,
such that
and
; (b)m = 0.1,
such that
and
; (c) m = 0.1,
such that
and
; (d) when
and m = 0.1, there is a semi-trivial periodic solution.
![Figure 4. The order-1 periodic solution of system (Equation18(18) dxdt=rx1−xk−βxydydt=y(βx−α),x<h,Δx=−nx(t)Δy=−my(t)+τ,x=h.(18) ) under the case h_<h<x∗, where r = 1, k = 10, β=0.025, a = 0.188, n = 0.3,h = 5.5. (a) m = 0, τ=15 such that τm>yT and F(T)>0; (b)m = 0.1, τ=7 such that τm>yT and F(T)<0; (c) m = 0.1, τ=2 such that 0<τm<yT and F(T)<0; (d) when τ=0 and m = 0.1, there is a semi-trivial periodic solution.](/cms/asset/f395416c-5bdb-462e-bbcf-c4f01705a030/tjbd_a_1937721_f0004_oc.jpg)
When ,
may be well defined or not well defined.
Let r = 0.5, k = 15, , a = 0.165, m = 0.03, n = 0.428, h = 7 be fixed, the calculation shows that
,
and
, which means
is well defined. Compared with Figure , Figure illustrates the existence of the order-1 periodic solutions which are the same as
.
Figure and Figure show that there must exist an positive order-1 periodic solution for system (Equation2(2)
(2) ) under the case
is well defined and
.
5.2. ![](//:0)
is not well defined
If is not well defined, then
.
Taking r = 0.5, k = 15, , a = 0.27, n = 0.3, m = 0.5 and h = 11, then
,
, which means
. The calculation gives
,
and
is not well defined. Further,
and
are obtained.
Let . Then
,
and
are satisfied. It can be seen from Figure (a) that there is an order-1 periodic solution which locates in a semi-ring domain. Take
, then
and
hold. As can be seen from Figure (b) that there does not exist periodic solution because the trajectories of the system ultimately approach the equilibrium
. While
, we can verify that
,
and
. We can see from Figure (c) that an order-1 periodic solution exists.
Figure 5. The existence of order-1 periodic solution of system (5.1) when and
is well defined, where r = 0.5, k = 15,
, a = 0.165, m = 0.03, n = 0.428 and h = 7. (a)
such that
and
; (b)
such that
and
; (c)
such that
; (d)
.
![Figure 5. The existence of order-1 periodic solution of system (5.1) when h_<x∗<h and F(T) is well defined, where r = 0.5, k = 15, β=0.025, a = 0.165, m = 0.03, n = 0.428 and h = 7. (a) τ=7 such that τm>yT and F(T)>0; (b) τ=1 such that τm>yT and F(T)<0; (c) τ=0.05 such that 0<τm≤yT; (d) τ=0.](/cms/asset/cac8a95b-b833-4a99-9c3f-bf9c6f33621a/tjbd_a_1937721_f0005_oc.jpg)
Figure 6. The existence of order-1 periodic solution of system (Equation18(18)
(18) ) when
and
is not well defined, where r = 0.5,k = 15,
,a = 0.27, n = 0.3, m = 0.5 and h = 11. (a)
such that
,
and
; (b)
such that
and
; (c)
such that
,
.
![Figure 6. The existence of order-1 periodic solution of system (Equation18(18) dxdt=rx1−xk−βxydydt=y(βx−α),x<h,Δx=−nx(t)Δy=−my(t)+τ,x=h.(18) ) when h_<x∗<h and F(T) is not well defined, where r = 0.5,k = 15,β=0.025,a = 0.27, n = 0.3, m = 0.5 and h = 11. (a) τ=30 such that τm≥yR−,I(R)>R− and PN(R+)>R−; (b) τ=7 such that τ≥YR− and τm<yR−; (c) τ=2 such that 0<τ<YR−,I(R)<R−.](/cms/asset/4ac312ab-e8d6-4533-a6db-bbd871404776/tjbd_a_1937721_f0006_oc.jpg)
Let r = 0.5, k = 15, , a = 0.165, n = 0.3 and h = 9 be fixed. Then we have
,
and
. Take
and m = 0.04, then
. When
, it can be seen from Figure (a) that the system has a limit cycle. Changing the initial value
, the trajectory approaches the equilibrium
after finite times of impulse impact. Let
and m = 0.5. Then
, and the system has a semi-trivial order-1 periodic solution (see Figure b). Given a small
, even since the initial value of y is zero, there is still a positive order-1 periodic solution (see Figure c).
Figure 7. The existence of periodic solution of system (Equation18(18)
(18) ) when
and
is not well defined, where r = 0.5, k = 15,
, a = 0.165, n = 0.3 and h = 9. (a)
, m = 0.04, such that
,
is taken by
and
,respectively; (b)
, m = 0.5, such that
,
, and
; (c)
, m = 0.5, such that
, where
.
![Figure 7. The existence of periodic solution of system (Equation18(18) dxdt=rx1−xk−βxydydt=y(βx−α),x<h,Δx=−nx(t)Δy=−my(t)+τ,x=h.(18) ) when h_<x∗<h and F(T) is not well defined, where r = 0.5, k = 15, β=0.025, a = 0.165, n = 0.3 and h = 9. (a) τ=0, m = 0.04, such that I(R)>R−, (x0,y0) is taken by (6.3,5.5) and (6.3,0.5),respectively; (b) τ=0, m = 0.5, such that I(R)<R−, (x0,y0)=(6.3,1.5), and limϵ→0+F(ϵ)=0; (c) τ=0.05, m = 0.5, such that I(R)<R−, where (x0,y0)=(6.3,0).](/cms/asset/811c459b-3a68-4714-8477-66beaea726ed/tjbd_a_1937721_f0007_oc.jpg)
From Figure , it can be seen that the existence of periodic solutions depends on the initial values, parameters and the structure of the system without impulse. Moreover, Figure (b) and Figure (c) show that a small τ can stimulate a positive order-1 periodic solution.
5.3. Stability of the order-1 periodic solution
Let be an order-1 ω-periodic solution of the system (Equation18
(18)
(18) ). Corresponding to system (Equation2
(2)
(2) ), we have
(19)
(19)
Proposition 5.1
Assume k>h and hold. If
, then the order-1 periodic solution
is orbitally asymptotically stable, where
is defined in
Discretizing the above integral, (Equation19(19)
(19) ) becomes
(20)
(20) where
, which is the step size, dividing the interval
into N subintervals.
Let r = 1, k = 10, , a = 0.188, n = 0.3, m = 0.03, h = 5.5, which means that
and
is well defined. According to Theorem 3.1, there exists an order-1 periodic solution
for system (Equation18
(18)
(18) ). Take the uniform step size
. When
, the computation gives
,
,
and
. Following Equation (Equation18
(18)
(18) ), we have
. Therefore, the order-1 periodic solution is orbitally asymptotically stable(see Figure a). Let
. Then
,
,
and N = 97, and hence
, which can be seen from Figure (b) that the order-1 periodic solution is orbitally asymptotically stable.
Figure 8. The stability of order-1 periodic solution of system under the case that and
is well defined, where r = 1, k = 10,
, a = 0.188, n = 0.3, m = 0.03, h = 5.5 and
such that
: (a)
such that
; (b)
such that
.
![Figure 8. The stability of order-1 periodic solution of system under the case that h_<h<x∗ and F(T) is well defined, where r = 1, k = 10, β=0.025, a = 0.188, n = 0.3, m = 0.03, h = 5.5 and Δti=0.05 such that |μ¯1|<1: (a) τ=30 such that F(T)>0; (b) τ=7 such that F(T)<0.](/cms/asset/bb35adb2-aa2e-4879-b9d3-6627d846e88e/tjbd_a_1937721_f0008_oc.jpg)
Fix r = 0.5, k = 15, and n = 0.3. Then
. Let a = 0.165, m = 0.03,
, h = 9. Then
is not well defined, and
,
,
. take
. The calculation shows
. Further, let a = 0.27, m = 0,
, h = 11 and
, the computation shows that
,
,
and N = 212, and hence
. Obviously, there is an orbitally asymptotically stable periodic solution of system, respectively, which can be seen from Figure (a) and Figure (b).
6. Conclusion and discussion
In the paper, the existence and existing domain, along with the stability of order-1 periodic solutions are considered. Further, we focus on the effects of parameters τ and m on the existence and stability of periodic solutions, which are performed theoretically and numerically. We construct three kinds of Bendixson domains with different structures, which overcome the difficulty, to some extent, that caused by the discontinuity of impulsive impact.
Our results show that the system (Equation4(4)
(4) ) without a positive periodic solution may have periodic oscillations under the effects of pulses. In this sense, the pulse action can ‘activate’ the state of periodic oscillation. When
,
must be well defined, and when
,
is not necessarily well defined. If
is well defined, then system (Equation2
(2)
(2) ) has an order-1 periodic solution, and when
, there must exist a positive order-1 periodic solution. If
is not well defined and the parameters τ and m are taken moderate values, the impulse effect will be invalid ultimately (see Theorem 3.4). If τ and
are small enough, the population of the predator species, corresponding to y of the periodic solution, is also small. When
is large enough, the periodic solutions locate in the sub-parallel domain (when
is well defined) or in the semi-ring domain (when
is not well defined), corresponding to a larger scale of predator species. When
, no matter
or
holds, there exist a semi-trivial periodic solution. However, when
, there may still exist other positive order-1 periodic solutions (see Figure (a)).
From biological control point of view, when the density of the prey reaches a threshold h, the two species are harvested. When , the system has a semi-trivial periodic solution, which means the predator becomes extinct, and the ecosystem is destroyed. Comparing Figure (c), Figure (c), Figure (b) with Figure (d), Figure (d), Figure (c), respectively, it can be seen that a small τ can ‘activate’ an positive order-1 periodic solution, which means the releasing of a small scale of predator can ‘stimulate’ an dynamic equilibrium of the ecosystem. When
, the dynamic equilibrium is stable, and we can implement the control of harvesting or releasing at every ω time, so as to maintain the ecological balance between the two species and to ensure the sustainable economic benefits.
Disclosure statement
No potential conflict of interest was reported by the author(s).
Additional information
Funding
References
- D. Bainov and P. Simeonov, Impulsive Differential Equations: Periodic Solutions and Applications, Longman Scientific & Technical, England, 1993.
- E. Borducchi, C. Cabral, K. Stephenson, J. Liu, and H. Dan, Ad26/mva therapeutic vaccination with tlr7 stimulation in siv-infected rhesus monkeys, Nature 540 (2016), pp. 284–287.
- J. Chavez, D. Jungmann, and S. Siegmund, A comparative study of integrated pest management strategies based on impulsive control, J. Biological Dyn. 12 (2018), pp. 318–341.
- T. Cheng, S. Tang, and R.A. Cheke, Threshold dynamics and bifurcation of a state-dependent feedback nonlinear control SIR model, J. Comput. Nonlinear Dyn. 14 (2019), pp. 209.
- R. Hakl, M. Pinto, and V. Tkachenko, Almost periodic evolution systems with impulse action at state-dependent moments, J. Math. Anal. Appl. 446 (2017), pp. 1030–1045.
- Z. He, C. Li, L. Chen, and Zh. Cao, Dynamic behaviors of the FitzHugh–Nagumo neuron model with state-dependent impulsive effects, Neural Netw. 121 (2019), pp. 497–511.
- X. Hu, J. Li, and X. Feng, Threshold dynamics of a HCV model with virus to cell transmission in both liver with CTL immune response and the extrahepatic tissue, J. Biological Dyn. 15 (2021), pp. 19–34.
- G. Jiang, Q. Lu, and L. Qian, Complex dynamics of a Holling type II prey–predator system with state feedback control, Chaos Solit. Fract. 31 (2007), pp. 448–461.
- T. Kar and S. Chattopadhyay, A focus on long-run sustainability of a harvested prey predator system in the presence of alternative prey, Comp. Rend. Biol. 333 (2010), pp. 841–849.
- L. Li, Ch. Li, and W. Zhang, Existence of solution, pulse phenomena and stability criteria for state-dependent impulsive differential equations with saturation, Commun. Nonlinear Sci. Numer. Simul. 77 (2019), pp. 312–323.
- M. Liao, Y. Liu, Sh. Liu, and A. Meyad, Stability and Hopf bifurcation of HIV-1 model with Holling II infection rate and immune delay, J. Biological Dyn. (2021). Available at DOI:https://doi.org/10.1080/17513758.2021.1895334.
- L. Nie, J. Shen, and C. Yang, Dynamic behavior analysis of SIVS epidemic models with state-dependent pulse vaccination, Nonlinear Anal. Hybrid Syst. 27 (2018), pp. 258–270.
- L. Nie, Zh. Teng, H. Lin, and J. Peng, The dynamics of a Lotka-Volterra predator-prey model with state dependent impulsive harvest for predator, BioSystems 98 (2009), pp. 67–72.
- M. Nowak and R. May, Virus Dynamics: Mathematical Principles of Immunology and Virology, Oxford University Press, Oxford, 2000.
- P. Simeonov and D. Bainov, Orbital stability of periodic solutions of autonomous systems with impulse effect, Int. J. Syst. Sci.19 (1989), pp. 2561–2585.
- S. Tang and R. Cheke, State-dependent impulsive models of integrated pest management (IPM) strategies and their dynamic consequences, Math. Biol. 50 (2005), pp. 257–292.
- X. Tang and X. Fu, On period-k solution for a population system with state-dependent impulsive effect, J. Appl. Anal. Comput. 7 (2017), pp. 39–56.
- S. Tang, X. Tan, and J. Yang, Periodic solution bifurcation and spiking dynamics of impacting predator–prey dynamical model, Int. J. Bifurcation Chaos 28 (2018), pp. 1850147. Article ID 1850147.
- J. Wang, H. Chen, Y. Li, and X. Zhang, The geometrical analysis of a predator-prey model with multi-state dependent impulses, J. Appl. Anal. Comput. 8 (2018), pp. 427–442.
- H. Wang, B. Dai, and Q. Xiao, Existence of order-1 periodic solutions for a viral infection model with state-dependent impulsive control, Adv. Differ. Eq. 2019 (1) (2019). Available at https://doi.org/10.1186/s13662-019-1967-x.
- H. Wang, Ch. Ou, and B. Dai, The existence and stability of order-1 periodic solutions for an impulsive Kolmogorov predator–prey model with non-selective harvesting, J. Appl. Anal. Comput. 11(2021), pp. 1348–1370.
- A. Wang, Y. Xiao, and R. Smith, Using non-smooth models to determine thresholds for microbial pest management, J. Math. Biol. 78 (2019), pp. 1389–1424.
- Ch. Xu, M. Liao, P. Li, Y. Guo, Q. Xiao, and Sh. Yuan, Influence of multiple time delays on bifurcation of fractional-order neural networks, Appl. Math. Comput. 361 (2019), pp. 565–582.
- G. Zeng, L. Chen, and L. Sun, Existence of periodic solution of order one of planar impulsive autonomous system, J. Comput. Appl. Math. 186 (2006), pp. 466–481.
- Q. Zhang, B. Tang, and S. Tang, Vaccination threshold size and backward bifurcation of SIR model with state-dependent pulse control, J. Theoretical Biol. 455 (2018), pp. 75–85.
- Y. Zhou, Ch. Li, and H. Wang, Stability analysis on state-dependent impulsive hopfield neural networks via fixed-time impulsive comparison system method, Neurocomputing 316 (2018), pp. 20–29.