#### Abstract

In this paper, an ecological model described by a couple of state-dependent impulsive equations is studied analytically and numerically. The theoretical analysis suggests that there exists a semitrivial periodic solution under some conditions and it is globally orbitally asymptotically stable. Furthermore, using the successor function, we study the existence, uniqueness, and stability of order-1 periodic solution, and the boundedness of solution is also presented. The relationship between order- successor function and order- periodic solution is discussed as well, thereby giving the existence condition of an order-3 periodic solution. In addition, a series of numerical simulations are carried out, which not only support the theoretical results but also show the complex dynamics in the model further, for example, the coexistence of multiple periodic solutions, chaos, and period-doubling bifurcation.

#### 1. Introduction

Since mathematical models represented by differential equations were introduced into biological and ecological systems, it has been proved that they are very useful tools to deal with practical problems such as emerging disease [1] and population dynamics between plankton and nutrient [2–4], which also promote the development of the theory of differential equation further. Especially, in the process, there is a significant development in impulsive differential equations [5–9] because many reality systems exhibit the abrupt jumps phenomena in population sizes during the evolution processes. Additionally, a theoretical method is urgently needed to provide a guide for the management of some practical problems such as pest outbreaks, phytoplankton blooms, vaccination, and so on. Accordingly, impulsive differential equation theory may prove to be one of the most potential theories to strategize about the approach of management for these problems, as impulsive differential equation can model abrupt jump behaviors induced by management, for example, some management occurs at fixed time, which can be modeled using the models with impulses at fixed times [10–13].

State-dependent impulsive differential equation, as a branch of impulsive differential equation, has the advantage to model the control behaviors depending on the population state. In recent years, the studies on state-dependent impulsive differential equation have been paid increasing attention [13–16]. Many studies devote to the properties of periodic solutions [17–20], including the existence, stability and periodicity, etc. It has been confirmed that the Poincaré map is a very useful tool to prove the existence of periodic solution in state-dependent impulsive differential equations [21–23]. In addition, Simeonov and Bainov [24] gave analogue of Poincaré criterion in ref. [25], which contributes to the stability of the periodic solution. Based on the Poincaré map, it can be proved theoretically that the positive periodic solution can bifurcate from the trivial periodic solution through a transcritical bifurcation [25]. Particularly, Chen [26] introduced the successor function a in semicontinuous dynamical system described by state-dependent impulsive differential equation to prove the existence of periodic solution. Using the successor function, Dai et al. [27] studied homoclinic bifurcation in a semicontinuous dynamical system [28]. The successor function has been extensively employed to investigate the existence of periodic solution and even stability and bifurcation in state-dependent impulsive differential equation [28–38]. These studies enrich the theory of state-dependent impulsive differential equation further.

In this paper, we will present a prey-predator two-species model with state-dependent impulsive control strategy to study the dynamics induced by impulsive control. Firstly, a model describing population dynamics of prey and predator is given below.where and denote the densities of prey and predator, respectively; is the intrinsic growth rate of the prey; is the carrying capacity of the prey; is the maximum growth rate; represents the half-saturation concentration; is the assimilation efficiency of predator; and is the specific mortality rate of predator. Here, we employ a Holling II functional response to present the interaction between prey and predator because some results suggest that the population dynamic following Holling II may exist in reality [2, 3, 39, 40].

For some species, when their densities are beyond a certain value, some seriously negative effects will emerge, for instance, algal bloom. For this reason, their densities should be kept below the critical value by harvest or/and released natural enemy. Obviously, the occurrence of control will rely on the density of the species. Here, under the assumption that the interaction between prey and predator follows model (1), we investigate the dynamics of impulsive control and control strategies for prey using state-dependent impulsive differential equation, and the model can be described aswhere denotes the critical value below which the density of prey should be kept. is the harvest rate of prey ; is the releasing amounts of predator ; and represents the harvest rate of predator or the level of negative effect caused by the harvest of prey on predator .

The rest of the paper is organized as follows. In Section 2, some notations, basic definitions, and lemmas are given. In Section 3, a semitrivial solution and its stability are studied; we also study the existence, uniqueness, and stability of an order-1 periodic solution; moreover, the boundedness of solution and some propositions of periodic solution are discussed as well. Whereafter, a series of numerical results are carried out to study the dynamics of model (2) further in Section 4. Finally, the paper ends with the conclusion in Section 5.

#### 2. Preliminaries

In model (1), there exist three equilibria, , , and , where and . When the condition holds, is a stable positive focus or node. We assume that the following condition holds throughout the paper: (*H*) .

The vertical isocline and the horizontal isocline can be defined as follows:

Furthermore, let be the phase set, and denotes the impulsive set. Under condition (*H*), we assume that intersects with the phase set at point and the impulsive set at point in the first quadrant when . We define and .

Let be any solution of model (2); then, we define the positive trajectory initializing at point for as .

In this paper, some notions are always mentioned, including phased point, impulsive point, and successor point. From works in [20, 23, 26], we state the following definition.

*Definition 1. *Take the point ; set as the trajectory initializing at point in model (2). If there exists such that the trajectory intersects the impulsive set at when , then a point exists in the phase set , where point . The point is called the phased point of point , and the point is called the impulsive point of point . The point is called the order- successor point of point , and the vector is called an impulsive line.

*Definition 2. *Set point ; then, we define a number axis in the phase , and set point as the origin of the number axis . , set the direction of vector as the positive direction of number axis , and set as the coordinate of point in the number axis . Then, a map can be defined as , where and is the order- successor point of point . The map is called the order- successor function, that is, .

*Definition 3. *Take a point in the phase set ; then, is a solution of model (2). If there exists a positive integer such that is the smallest integer for , then the solution is called order- periodic solution.

Lemma 1 (see [26]). *The successor function is continuous, where the point .**Obviously, we can obtain the following lemma by Lemma 1.*

Lemma 2. *If there exist two points and in model (2) and such that , then model (2) has an order-1 periodic solution.**For the point , we always assume the point is the order- successor point of the point , and the point is the impulsive point of point , where .*

#### 3. Main Results

When the values of parameters in model (1) are given, the dynamics of model (2) only depend on the parameters , , , and . In this paper, we are interested in the effect of impulsive control on dynamics of model (2). Hence, we will investigate the dynamical behaviors of model (2) under the condition that the parameters of model (1) are fixed.

For points and , it is obvious that(1) if .(2) if and .

In this section, we just show the proof with respect to case in the figure of illustration. The proof corresponding to other cases is similar, which has been omitted.

By model (1), it is obvious that implies . Additionally, always holds when .

In model (2) with , it is obvious that the following semitrivial periodic solution exists for when the predator ,where the period . Set . Then, we can get the following theorem.

Theorem 1. *Under condition ( H), if and , then semitrivial periodic solution is globally orbitally asymptotically stable in model (2) for any and .*

*Proof. *For , solutions of model (2) will enter the space eventually. Hence, we study the stability of the semitrivial periodic solution in this space.

For any , set and . Then, we assume that the trajectory intersects with impulsive set at point , and under condition (*H*). Set ; then, for any point where , we have under condition (*H*) because and (see Figure 1(a)), which suggests that the distanceholds for all , where . Hence, according to the definition of orbital stability in ref. [40], the semitrivial periodic solution is orbitally stable.

For any , take a sequence of number , where and , such that . Then, the sequence of number is a strictly decreasing positive sequence because and . Hence, there exists a , such that , but . Set , and let point be the order-1 successor point of point if , where . Then, for any point (i.e., ), its order-1 successor point is below the point (see Figure 1(b)), that is, .

Because , there exists a positive integer such that for . According to above analysis, there exists a positive integer such that when , which contradicts with . Hence, , which means that the semitrivial periodic solution is orbitally attractive. So, the semitrivial periodic solution is orbitally asymptotically stable. Since the attraction domain is , the semitrivial periodic solution is globally orbitally asymptotically stable.

This completes the proof.

Suppose that is a positive constant, and then we have the following theorem for any .

**(a)**

**(b)**

Theorem 2. *Under condition ( H), if and , then there exists a unique order-1 periodic solution in model (2) for any and .*

*Proof. *For any , Figure 2(a) displays the relation between point and its order-1 successor point under condition (*H*), where . Obviously, is an order-1 periodic solution if . If , we have and , where is the successor point of point . According to Lemmas 1 and 2, there exists an order-1 periodic solution in model (2).

If , we take a point in phase set . Then, let the point be the order-1 successor point of point and point be the impulsive point of point , and (see Figure 2(a)), so . By , there exists an order-1 periodic solution in model (2) using Lemmas 1 and 2.

Hence, there exists an order-1 periodic solution for any and when and .

Suppose there are two order-1 periodic solutions in model (2), which are and , respectively. Let point and be the phase point and impulsive point of , respectively, and point and be the phase point and impulsive point of , respectively. Without loss of generality, set . Then, the relation among , , and is (i) (); (ii) (); and (iii) ().

Figure 2(b) shows case (i); it is obvious that and cannot coexist by disjoint property of trajectory of model (1). Hence, the order-1 periodic solution is unique for this case.

For case (ii), because , we only need to focus on the region . We rewrite the model (1) in the region as the following initial-value problem:where .

Obviously, the function is continuously differentiable and monotone decreasing in on the space . Then, the solution of model (6) is equivalent to the following integral equation:By equation (7), setwhere ; then, and . Because , we have on the space by the disjoint property of trajectory of model (1). Hence,On the other hand, because both and are order-1 periodic solutions, we have by model (2), which is a contradiction with inequality (9). Hence, the order-1 periodic solution is unique for this case.

For case (iii), the proof is similar with case (ii), so it is omitted.

Hence, there exists a unique order-1 periodic solution in model (2) for any and when and .

This completes the proof.

From Theorem 2, we can find that the existence and stability of order-1 periodic solution depend on parameter . Suppose point is the phase point of the order-1 periodic solution, then we can get the following proposition.

**(a)**

**(b)**

Proposition 1. *Under condition ( H), when and , is a monotonously increasing function with respect to .*

*Proof. *Let be the angle between the impulsive line and the -axis. Then, , so . Hence, is a monotonously decreasing function in . Given an order-1 periodic solution , let points and be the phased point and impulsive point of , respectively; then, guarantees the existence of the order-1 periodic solution .

For any , when , the order-1 periodic solution will disappear. Let point be the order-1 successor point of point ; then, because is a monotonously decreasing function. In addition, and . Therefore, there exists an order-1 periodic solution whose phase point belongs to . By the uniqueness of order-1 periodic solution, is a monotonously increasing function in .

This completes the proof.

By Theorem 2, we know that the trajectory is an order-1 periodic solution of model (2) when . Then, we can get the following theorem.

Theorem 3. *Under condition ( H), if , then the order-1 periodic solution is globally orbitally asymptotically stable in model (2) for any and when .*

*Proof. *When , the order-1 periodic solution of model (2) must belong to space under condition (*H*). denotes an order-1 periodic solution of model (2), whose phased point and impulsive point are point and point , respectively. In addition, the solutions of model (2) will enter the space eventually when . Hence, we only need to study the stability of order-1 periodic solution on the space .

For any , we assume that the trajectory meets impulsive set at point , where (see Figure 3(a)); then, . Take ; for any point , where and , we assume that the trajectory intersects with impulsive set at point ; then, ; otherwise, another order-1 periodic solution exists above , which contradicts with the uniqueness of order-1 periodic solution. Hence, the distanceholds for all according to equation (9).

On the other hand, for any , given , when , we still haveusing the same method. Hence, according to definition of orbital stability in ref. [41], the order-1 periodic solution is orbitally stable.

For any , take a numerical sequence , where and , such that . Due to disjoint of impulsive line, there are two cases on the numerical sequence when and : (1) is a strictly decreasing positive sequence; (2) is a strictly increasing positive sequence.

For case (1) (see Figure 3(b)), we have and . Let ; there exists a point when , where , such that point is the order-1 successor point of point . For any point (i.e., ), its order-1 successor point is below the point . Let point be the order-1 successor point of , then .

Because and , there exists a positive integer such that for . According to above analysis, there exists a positive integer such that when , which contradicts with . So, . Likewise, also holds for case (2), which says the order-1 periodic solution is orbitally attractive.

Hence, the order-1 periodic solution is orbitally asymptotically stable. Since the attraction domain is , the order-1 periodic solution is globally orbitally asymptotically stable.

This completes the proof.

**(a)**

**(b)**

Proposition 2. *Under condition ( H), when and , if there exists an order- periodic solution in model (2), setting points as the phase points of the order- periodic solution, then .*

*Proof. *If , it is obvious that there is no order- according to the proofs in the previous theorems. Hence, we only need to consider . Let points be the impulsive point of the order- periodic solution and ; then, . Because , for any , if , then by the disjoint of trajectory of model (1).

For , the relationship among , , and is one of the following cases because of the disjoint of the trajectory of model (1) and the disjoint of the impulsive line of model (2): For , we get , so using . Hence, we claimWhen , it is obvious that inequality (12) holds. Now, suppose inequality (12) holds when . In the following, we will prove that inequality (12) also holds when . Because , we have . Consequently, , that is, inequality (12) holds when . Hence, inequality (12) holds when .

In this case, obviously, for an order- periodic solution, we have when , which contradicts to . Hence, the order- periodic solution does not exist. Using similar procedure, it is proved that the order- periodic solution does not exist for cases , , and . So, there is no order- periodic solution in model (2) when . Therefore, when an order- periodic solution exists in model (2).

This completes the proof.

By Proposition 2, the following theorem will exist.

Theorem 4. *Under condition ( H), if and , then there exists a globally orbitally asymptotically stable and unique order-1 periodic solution or a unique order-1 periodic solution and an order-2 solution coexists in model (2) for any , and there is no order- periodic solution.*

*Proof. *By Proposition 2, it is obvious that there is no order- periodic solution in model (2) because for any . If order-2 periodic solution does not exist, then a globally orbitally asymptotically stable and unique order-1 periodic solution exists using Theorems 2 and 3. If an order-2 periodic solution exists, then the order-2 periodic solution coexists with a unique order-1 periodic solution by Theorem 2.

This completes the proof.

Theorem 5. *Under condition ( H), if , for given , then all the non-negative solutions of model (2) which start in are ultimately bounded, and the bound depends on parameter , and .*

*Proof. *By model (2) and the space , obviously, and hold. In addition, set and . Then, , so we can obtain by . Hence, , which obviously depends on , , , and .

This completes the proof.

From Theorem 5, it is not difficult to find that there is a rough estimate of the bound, and we have the following remark.

*Remark 1. *Under condition (*H*), if , then and .

Proposition 3. *Under condition ( H), if there exists a point such that in model (2) when , then there exists an order- periodic solution, where the positive integer is a factor of the positive integer , including 1 and .*

*Proof. *Suppose positive integer is not a factor of the positive integer , then let be the remainder of and be the quotient of . Because there exists an order- periodic solution in model (2), without loss of generality, assume . Obviously, we have , but . Hence, the equation , which is a contradiction with because . So, the positive integer is a factor of the positive integer . When the positive integer is a factor of the positive integer , is obvious if there exists an order- periodic solution in model (2).

This completes the proof.

Proposition 4. *Under condition ( H), if there exists a point below point such that holds when , , , and , then model (2) has an order-3 periodic solution when , and the order-3 periodic solution coexists with an order-1 periodic solution.*

*Proof. *Setting , we have and . Hence, there exists a point such that because , so an order- periodic solution exists by Proposition 3. However, owing to , it is obvious that there is no order-1 periodic solution in space by the proof of Theorem 2. Thus, the periodic solution is the order-3 one. Furthermore, because , there exists an order-1 periodic solution according to the proof of Theorem 2, and its phased point is above point . Therefore, the order-3 periodic solution coexists with an order-1 periodic solution.

This completes the proof.

Although it is not proved that point must exist in Proposition 4, it provides a method to search an order-3 periodic solution for numerical analysis. Furthermore, the existence of order-3 periodic solution means that order- periodic solution and chaos exist in model (2) by Sarkovskii’s theorem [42] and Li and Yorke’s theorem [43].

#### 4. Numerical Results

In the previous section, we have analyzed the properties of periodic solution in model (2), including the existence, uniqueness, stability, and boundedness. In this section, we will further show the complex dynamics induced by impulsive control using the numerical simulations. Since the focus of this paper is the role of impulsive control in the population dynamics, the parameters of model (1) are fixed: , and we can obtain the equilibrium with index +1. Obviously, condition (*H*) always holds under this parameter set.

By Theorem 1, there only exists a semitrivial periodic solution in model (2) when and . Given , and , we get a semitrivial periodic solution as shown in Figure 4(a), and its period is about 1.309. In addition, we also obtain the theoretical value of the period by , which is the same as the numerical one. When and , the point is fixed at . The numerical simulation shows that all the successor points of any trajectory except the semitrivial solution are below initial point in model (2) and these trajectories converge towards the semitrivial periodic solution (see Figure 4(b)).

**(a)**

**(b)**

When , however, the semitrivial periodic solution disappears. Then, an order-1 periodic solution emerges and it should be unique by Theorem 2. In phase set , because the order-1 successor point of any point above point is below point , we compute the order-1 successor function where for different values of (see Figure 5(a)), where . In the rest of this section, we will employ this control parameter set unless otherwise specified.

**(a)**

**(b)**

From Figure 5(a), it is not difficult to find that the equation has a unique root, which means that there exists a unique order-1 periodic solution. Hence, we carry out some numerical simulations with respect to these roots of , and the results verify that these solutions are order-1 periodic solutions as shown in Figure 5(b). But we find that some of these order-1 periodic solutions are unstable, which indicates that there may exist order- periodic solution for some values of .

In order to illustrate the existence of order- periodic solution, we further compute order-1 successor function, order-2 successor function, and order-3 successor function, and the results are very interesting and are shown in Figure 6(a). From Figure 6(a), it is easy to find that the equation has a unique root, the equation has three roots, and the equation has seven roots. Furthermore, we can find that the root of is the common root of three equations: , , and , which means that there exists an order-1 periodic solution corresponding to this root. Except this root, the equation still has two roots which indicate an order-2 periodic solution exists in virtue of the uniqueness of order-1 periodic solution. For order-3 successor function , the equation has six roots except for the root of .

**(a)**

**(b)**

**(c)**

Except for the root of , it is not difficult to find that and are positive or negative definite in a small neighborhood of each root of , which means that there may exist two order-3 periodic solutions in model (2). Thus, we perform some numerical simulations to show the solutions corresponding to these roots. The results show that a stable order-3 periodic solution coexists with an unstable order-3 periodic solution, an unstable order-2 periodic solution, and an unstable order-1 periodic solution (see Figure 6(b)). Additionally, the numerical solutions for prey are given to further exhibit these periodic solutions, as shown in Figure 6(c).

The above results indicate that the order-1 periodic solutions are not always stable. Using order-1 and order-2 successor functions, we compute with respect to the phased point of order-1 periodic solution to simulate its variation as parameter increases (see Figure 7(a)). From Figure 7(a), the order-1 periodic solution is stable before parameter reaches line , and the order-1 periodic solution is still stable near when parameter is beyond . This result is in accord with Theorem 3. Keeping parameter increasing, the order-1 periodic solution loses its stability until it is beyond . It is well worth noting that model (2) only has order-1 and order-2 periodic solution when parameter is between and , which agrees with Theorem 4. Furthermore, Figure 7(a) demonstrates that the order-1 periodic solution is unique, which supports the result in Theorem 2. Additionally, Figure 7(a) also shows that is a monotonously increasing function with respect to , which is consistent with Proposition 1.

**(a)**

**(b)**

In addition, our numerical results also indicate that the phased point of order-1 periodic solution rises along phase set as parameter increases. In order to present the stability of order-1 periodic solution, we compute two groups of order- successor points whose initial point is on the sides of the phased point of order-1 periodic solution, respectively. It is revealed that these points monotonously converge towards the phased point of order-1 periodic solution (see Figure 7(b)). Hence, this order-1 periodic solution is stable due to the nonexistence of other periodic solutions.

In order to display complex dynamics, we chose parameter as the controlled parameter and ran numerous simulations using a wide range of parameter to show the bifurcation diagram (see Figure 8(a)). Since the stability of order-1 periodic solution has been discussed above (Figure 7(a)), we will present the dynamics when the order-1 periodic solution is unstable. From Figure 8(a), increased parameter can generate chaos by a cascade of period-doubling bifurcation. Keeping parameter increasing, an order-3 periodic solution bifurcates from chaos. Subsequently, chaos occurs again with the increase of parameter . Finally, chaos disappears again via a cascade of inverse period-doubling bifurcation at higher values in the simulated range of parameter . To confirm the occurrence of chaos further, we calculate the largest Lyapunov exponents in the same range of parameter as Figure 8(a), which agrees with the bifurcation diagram (see Figure 8(b)). In addition, dozens of time-series of population for differently initial values are given in Figure 8(c) to further show chaos, which suggests that all states described by these time-series are in disorder.

**(a)**

**(b)**

**(c)**

#### 5. Conclusion

In this paper, we proposed an ecological model with state-dependent impulsive control strategy to investigate the effect of impulsive control on prey-predator dynamics. Theoretical results suggest that model (2) has a globally orbitally asymptotically stable semitrivial periodic solution when parameter and , where condition (*H*) holds to guarantee existence of positive equilibrium . When , the semitrivial periodic solution disappears, but the theoretical analysis shows that a unique order-1 periodic solution emerges for any and .

Furthermore, we find that the order-1 periodic solution is globally orbitally asymptotically stable when the value of is below . Additionally, when the value of is beyond , theoretical analysis reveals that a globally orbitally asymptotically stable order-1 periodic solution exists or coexistence of a unique order-1 periodic solution and order-2 periodic solution occurs.

Additionally, the existence of order- is theoretically analyzed based on equation , which demonstrates that an order- periodic solution exists when , and is a factor of the positive integer , including 1 and . According to this result, we discuss the existence of order-3 periodic solution. Although we cannot get the exact range of parameter with respect to the existence of order-3 periodic solution, our results imply that the order-3 periodic solution may exist for , which even means that chaos may occur.

In order to study dynamics induced by impulsive control further, a series of numerical simulations are carried out, which agreed with the theoretical results well. Numerical simulations show that order-1, order-2, and order-3 periodic solutions can coexist. Specially, we find the coexistence of double order-3 periodic solutions. Moreover, numerical bifurcation analysis shows that an unstable order-1 periodic solution and a stable order-2 periodic solution can coexist for , but there exists only a stable order-1 periodic solution when is beyond some critical value. In addition, numerical simulations show that chaos occurs. Obviously, dynamics induced by impulsive control is much richer.

The theoretical and numerical results predict that the population can be controlled using the state-dependent impulsive control strategies. Because , the method by impulsive harvesting population was a unique choice which can keep the density of population below the critical value. Our results shows that the control period is much lower when the released amount is much smaller, which means that the control frequency is much higher. But from the viewpoint of management for some ecological problems (e.g., pest outbreak, phytoplankton blooms, etc.), it is always hoped that the control period is much longer. Thus, much more predators are needed to be released. However, our results also indicate that chaos may occur for some values of , which is unfavourable for prediction of control period. Hence, the released amount of predator should be chosen carefully. Obviously, our results are much more interesting, and we expect that these results are helpful to the management of ecological problems.

#### Data Availability

The data used to support the findings of this study are included within the article.

#### Conflicts of Interest

The author declares that there are no conflicts of interest.

#### Acknowledgments

This study was supported by the Zhejiang Provincial Natural Science Foundation of China (grant no. LQ18C030002) and the National Natural Science Foundation of China (grant no. 61901303).