1. Introduction
Pest control is a serious challenge facing agricultural production today, and it is also a key factor affecting agricultural production efficiency, increasing income, ensuring the quality of agricultural products, and the sustainable development of agriculture. A more effective control strategy is to use the biological control method of releasing natural enemies. Consequently, studying predator-prey systems [
1] is essential to mitigate or manage the types and quantities of pests that pose a threat to human development, thereby reducing economic losses. However, in reality, pest outbreaks result in countless crop losses each year. These outbreaks not only cause significant harm to agricultural production but also facilitate the spread of diseases. Consequently, pest control has become a pressing concern for agricultural and disease control authorities.
Integrated pest management (IPM) is currently recognized as one of the most effective strategies for pest control [
2,
3,
4,
5]. It capitalizes on the predatory behavior of natural enemies to regulate pest populations. In specific agricultural ecosystems, the natural enemy populations often fall short in comparison to pest numbers. Consequently, pest populations on farmlands tend to exhibit a natural growth trend due to the abundance of habitat and food resources. Considering the economic significance of farmland, it becomes necessary to artificially introduce natural enemies. Recognizing that complete eradication of pests is typically unattainable. Instead, efforts are focused on limiting pest populations below the economic threshold, the point at which they cause significant economic harm [
6,
7,
8,
9].
Artificial release of natural enemies is one of the most effective pest control methods, but it also introduces complexities to the ecosystem. As pests are preyed upon by natural enemies, they exhibit anti-predator behavior in response to these threats [
10,
11,
12,
13]. However, natural enemies are not always capable of completely suppressing pest populations. The outcome depends on the density ratio of the two and the effectiveness of anti-predator behavior.
One prevalent form of anti-predator behavior involves role reversals between predator and prey, where juvenile predators may be attacked by adult prey, thus reducing the survival rate of the predators. In this process, prey also incur losses as attacking juvenile predators increases their own risk of predation. For instance, when left alone, Crossaster papposus, which feeds on sea urchins, is attacked and consumed by multiple sea urchins. The specific process of anti-predator behavior is as follows: a leading sea urchin attempts to approach Crossaster papposus and attaches itself to one of its tentacles. However, this process carries risks, as the leading sea urchin may become food for Crossaster papposus if it fails. Once the first sea urchin starts chewing on the tentacle, other sea urchins also attach themselves to the remaining tentacles, gradually consuming Crossaster papposus bit by bit. When prey invests excessively in anti-predator behavior, the costs for both predator and prey become significant.
During the long-term evolutionary process, the emergence of anti-predator behavior can exacerbate the complexity of predation relationships, which poses challenges to the study of pest-natural enemy ecosystems in farmland [
14]. To capture the complexity of these systems, a two-dimensional differential equation approach, known as the ODE model, is employed to depict the trend of changes in the number of two populations. This mathematical modeling technique for pest management decisions [
15,
16] builds upon the classic predator-prey model, which utilizes differential equations to represent the average growth rates of natural enemies and pests.
Furthermore, it is found that the relative growth rate described by linear function is not adaptive; considerable attention and research have been devoted to predator-prey models that incorporate various functional response functions [
17,
18,
19,
20,
21,
22]. Scholars have recognized the significance of considering the influence of anti-predator behavior. In this regard, Ives and Dobson [
13] proposed the following predator-prey model, aiming to account for the impact of anti-predator behavior
Here, the variables
and
represent the density of the prey and predator populations at time
t, respectively. The parameter
represents the intrinsic growth rate of the pest, while
m represents the death rate of the predator. The carrying capacity of the pest population is denoted by
K. Additionally, the nonnegative constant
v signifies the extent of the pest’s investment in anti-predation behavior, with higher values of
v resulting in lower predation rates. The efficiency of the anti-predator behavior is represented by
;
c represents the conversion rate of prey to predators;
corresponds to the Holling II functional response.
In pursuit of optimal economic development and ecological preservation, impulsive control techniques are employed to reduce pest populations below a specific economic threshold. State-dependent feedback control strategies have recently emerged as a valuable approach in integrated pest management. In many practical scenarios, the decision to apply impulsive control is often contingent upon the state of the pest population, indicating that control tactics are implemented only when the model reaches a certain state. When pest populations reach the economic threshold, measures are taken to impulsively control their numbers to a certain extent. Subsequently, pest populations begin to grow again until they reach the economic threshold once more. This cyclical process of applying impulsive control is repeated to consistently keep pest populations below the economic threshold. Consequently, both economic threshold and state feedback control are considered to offer a rational framework for describing the dynamics of predation systems.
To gain a better understanding of the effects of anti-predator behavior on the dynamics of two populations [
23,
24,
25], we propose an impulsive semi-dynamical system [
26,
27,
28]
Here,
represents the economic threshold. The parameter
p denotes the killing rate of pests when the pesticide is applied, with
. Additionally,
represents the number of released predators. In the event of immediate application of control tactics,
and
instantly becomes
and
, respectively. In addition, to maintain the inhibitory effect of natural enemies on pests, it is necessary to retain the number of natural enemies as much as possible. Therefore, insecticides with no or minimal impact on natural enemies were selected.
Based on model (
2), the IPM strategy or intermittent control strategy is as follows: When pest density increases to the threshold
(i.e.,
), chemical and biological control measures are implemented to target both populations. On the other hand, when the pest density is updated to a value less than the prescribed threshold
(i.e.,
), no control policies are implemented, and model (
2) reverts to the free model (
1). The purpose of this article is to analyze the factors affecting pest outbreaks by studying the population dynamics of model (
2) and to use reasonable pulse control to suppress the number of pests below the economic threshold ET. Its innovation is reflected in the impact of anti-predator behavior on both pests and natural enemies, as well as the combination of anti-predator behavior, functional response, and pulse control.
The remaining sections of the paper are organized as follows: In
Section 2, the main properties of an ODE model (
1) will be introduced.
Section 3 focuses on the exact domain of the impulsive set and phase set of model (
2). Additionally, we provide the definition of the Poincaré map.
Section 4 and
Section 5 are dedicated to the investigation of the boundary and internal order-1 periodic solutions, as well as the order-k periodic solution (
) of the proposed model. Furthermore, we explore the complex dynamics of the model, including multiple bifurcation phenomena and chaos, through computational simulations.
Section 6 conclude the paper by discussing our findings and their implications for pest control strategies.
2. Basic Properties of an ODE Model (1)
To investigate the dynamics of the state-dependent feedback control model (
2), it is essential to examine the basic properties of model (
1). The following are the two isoclines of model (
1)
here we define
. When
, the natural enemies become extinct, and
x satisfies the equation
. Considering the condition
, it indicates that the prey investment
v satisfies
; we obtain two boundary equilibriums:
and
. The values of
and
are
it is evident that
, and
.
When
, model (
1) admits an internal equilibrium
with
in order to ensure that the internal equilibrium
is positive, we need the inequality
to hold. The necessary and sufficient condition for this inequality to have a solution is
and
. Thus, when
or when the investment
, the system (
1) does not possess an internal equilibrium. Additionally,
is equivalent to the condition
. To simplify calculations, we define
. Therefore, the internal equilibrium
can also be expressed as follows:
In particular, when
, model (
1) exhibits a zero equilibrium
and a boundary equilibrium
. This demonstrates that without human interference, as time approaches infinity
, the predator population gradually tends towards extinction. Unless stated otherwise, in order to guarantee the existence of a positive internal equilibrium, we assume that
.
To perform a qualitative analysis of the dynamics of model (
2) under varying values of
v, we can define it as a hybrid dynamic system that combines discrete events and continuous system interactions. Specifically, when
, the model corresponds to a classical Lotka-Volterra predator-prey model [
29] with a Holling II functional response. This variant exhibits rich dynamics depending on the relative positions of the equilibria and the choice of economic threshold. These dynamics may include multiple globally stable equilibria, limit cycles, and transcritical bifurcations, among others.
On the other hand, when
, the predation rate decreases due to the prey’s investment in anti-predator behavior. This investment indirectly affects the prey’s per capita growth rate, either positively or negatively. Based on these considerations, this paper aims to theoretically and numerically examine the dynamics of model (
2).
If
, we have reformulated model (
1) as
The properties of model (
3) are presented in [
30]. In this paper, we primarily focus on analyzing the dynamics of model (
1) when
. By examining the relationship between stability and eigenvalues at equilibriums in a planar system, we obtain the following theorem.
Theorem 1. (i) If , model (1) admits no positive internal equilibrium. In this case, is unstable, and is stable; (ii) If and and , then is a stable focus; (iii) If and , then is a stable node. Here we denote Proof. Consider equations
and
in the model (
1) as P and Q, and take the derivative of x and y to obtain the Jacobian matrix, which is defined as
Substituting
into the Jacobian matrix obtains
it is easy to get that
. If
, then
and
, which indicates that
is an unstable focus or node. If
, then
, which indicates that
is a saddle. Then substituting
into the Jacobian matrix gives
It is evident that
. If
, then
, indicating that
is a saddle. On the other hand, if
, then
and
, implying that
is stable, as shown in
Figure 1A. Analogously, the internal equilibrium
is substituted into the Jacobian matrix as follows
Here, we have
. Therefore, when
and
and
,
is a stable focus, as illustrated in
Figure 2. On the other hand, when
and
and
,
is a stable node, as shown in
Figure 3. □
In addition, for the case of
, as shown in
Figure 1B, due to the small value of
, the trajectory of the system is not completely in the first quadrant, which is clearly not in line with the actual biological significance, so it will not be discussed.
3. Establishment of Poincaré Map
Define
. Considering its biological significance, the subsequent research will focus on the values of
x and
y within this set. Moreover, given the condition
, we define an open set
. Now, based on the third and fourth equations of model (
2), we can proceed with the following definitions
By combining
and
, we can find their intersection point, denoted as
, here
Analogously, we can obtain the intersection point of
and
, denoted as
, here
The impulsive set
M is a closed subset of
; we can define
based on model (
2). Consequently, the phase set can be defined as follows
where
; an impulsive point of
is denoted as
. We assume that the initial value
belongs to
N in this paper.
Next, we discuss the precise domain of the pulse and phase sets; this can be achieved by excluding the points in the impulsive set M that do not experience impulsive effects. Since the impulsive set M represents the largest interval in the ordinate, we can obtain the exact domain of the impulsive set by removing those points.
In Case (i) of Theorem 1, model (
2) possesses a unique stable boundary equilibrium, denoted as
, as illustrated in
Figure 4A. If
, a trajectory
within the system touches the line
, and the point of tangency is denoted as
, which represents the intersection of
and
. It is evident that trajectory
will also intersect the line
at a certain point, denoted as
. It is worth noting that
is satisfied. Based on these considerations, the impulsive set is
and the phase set is
where
. After such a definition, any solutions of model (
2) that originate from the interior of trajectory
will either not reach the impulsive set or intersect with it, leading to a mapping to the phase set
N after a single impulsive effect. If
, the solution starting from the phase set may directly approach the stable equilibrium
without reaching the impulsive set. In this case, the exact domain of the impulsive set and phase set cannot be determined.
In Case (ii) of Theorem 1, model (
2) possesses an unique stable focus, we consider two cases
Case
, as shown in
Figure 4B, the trajectory
is tangent to
, as observed in the phase portrait. It can be seen that the point of tangency
represents the intersection of
and
. The intersection of trajectory
and
on the left side is denoted as
. If
, employing the same analysis method as in case (i), we can define
and
. On the other hand, if
, the trajectory
must intersect the line
at two points, marked as
and
respectively. When
, any solutions of the system originating from
will not reach the impulsive set. Therefore, in this scenario, the impulsive set can be defined as
M and
where
. For Case
, similar to Case (i), we can derive the impulsive set
and the phase set
.
In Case (iii) of Theorem 1, it is also discussed in two cases, and . In case , if , it is impossible to determine whether the solution of the system will reach the impulsive set. Therefore, the impulsive set and phase set cannot be defined in this case. In case , if , then we derive the impulsive set and the phase set .
Based on the above analysis, we have determined the exact domains of the impulsive set and phase set for model (
2) in all cases. The results are summarized in
Table 1. In the table, the cases of the same pulse set and phase set are defined as
. Note that since the pulse set of case
cannot be determined, this case will not be discussed in the following sections.
Consider an arbitrary solution with an initial value . This solution undergoes a single impulsive effect at time , a second impulsive effect at time , and so on, until the kth impulsive effect at time . We denote the points associated with these impulsive effects as and . Thus, we can express as a continuous function: . Consequently, we have . Since and lie on the same trajectory, can be uniquely determined by , and can be uniquely determined by . Specifically, we can express as . According to the Cauchy-Lipschitz theorem, the function f is continuously differentiable in cases and .
In order to determine the domain of the Poincaré map, we need to establish the condition under which a solution initiating from does not experience impulsive effects. In case , it is straightforward to determine the domain and range of the Poincaré map as and , respectively. In Case , we observe that when , the initial point lies within the trajectory . In this scenario, the trajectory originating from the initial point will not reach the impulsive set. Consequently, the domain is , and the range is . To ensure the continuous occurrence of impulsive effects, it suffices to select the initial point from .
Therefore, except for the condition
in case
, any solutions of model (
2) initiating from
will reach
after a single impulsive effect. We can therefore obtain the domain and range of the Poincaré map as follows:
The lemma regarding the Poincaré map within the phase set domain for cases and is presented.
Lemma 1. The Poincaré map generated by the sequence of impulsive points in model (2) can be defined as follows:where is defined by Formula (4). For Cases
and
, since the function
f is continuously differentiable with
,
is also continuously differentiable with
. Therefore, the fixed point of
corresponds to an order-1 periodic solution within model (
2).
5. Existence and Stability of Order-k Positive Periodic Solution
The previous section has systematically described the boundary solution of the model (
2), and our primary focus of this section is on examining the properties of the internal solution when
, especially the existence and stability of the order-1 periodic solution as well as the existence of the order-k periodic solution. This aim can be achieved by demonstrating the presence and stability of the fixed point of the Poincaré map
.
Firstly, there is a threshold condition for the existence of order-1 periodic solution in case .
Theorem 3. In case , possesses a fixed point. As a result, the model (2) exhibits an order-1 periodic solution. Proof. For this particular case, the trajectory starting from crosses , and the point where they intersect is labeled as . Where is the tangent point between curve and line . Since both point and point T lie on the same curve, it can be concluded that .
When a single impulsive effect occurs,
T is mapped to
, and obtain
from the pulse function. Accordingly, since
is determined by
, we assume the existence of
such that
. Then, when
, the point
and the point
coincide, establishing a fixed point for
. Consequently, the trajectory
is treated as an order-1 periodic solution for model (
2).
If
, as illustrated in
Figure 7A, the position of the impulsive point
is higher than that of point
. Consequently, we have
. This leads to the following inequality:
The trajectory originating from
crosses
, and the point where they intersect is labeled as
. Since any two trajectories of the plane system cannot intersect, it follows that the position of point
is lower than that of point
T. After undergoing a pulse effect,
is mapped to
. As
, we have
. Consequently, the position of point
is lower than that of point
. Hence,
. This leads to the following inequality:
By combining (
8) and (
9), we conclude that there possesses a fixed point within the interval
for
.
If
, the position of the impulsive point
is lower than that of point
. Consequently, we have
. This leads to the following inequality:
Select a point
on the line
that satisfies the condition
. The trajectory originating from
crosses
, and the point where they intersect is labeled as
, where
. After undergoing a pulse effect,
is mapped to
. Now, considering
, we obtain the following inequality:
Combining (
10) and (
11), we conclude that there possesses a fixed point within the interval
for
.
The above process illustrates that no matter what value
takes, the system (
2) always has an order-1 periodic solution. □
Next, explain the conditions for the existence of an order-1 periodic solution in case . Before that, define a key point , which is the impulsive point that appears on the trajectory and initiates from point after a single impulsive effect.
Theorem 4. In case , if or , it can be observed that the Poincaré map possesses a fixed point. As a consequence, model (2) exhibits an order-1 periodic solution. Proof. For this situation, let’s mark
as the second intersection of
and
. If
, as depicted in
Figure 7B, since
is increasing on the right area of
, the position of
is higher than that of
. In addition,
is mapped to
after undergoing a pulse effect. Considering
, it follows that the position of point
is also higher than that of point
.
If
, then the trajectory
forms an order-1 periodic solution of the system (
2). If
, we have
. This leads to the following inequality:
Furthermore, the trajectory originating from
crosses
, and the point where they intersect is labeled as
. By virtue of the uniqueness of the solution, we deduce that the position of
is lower than that of
. After undergoing a pulse effect,
is mapped to
. Clearly, the position of
is also lower than that of
, hence
. This yields
Combining (
12) and (
13), we conclude that there possesses a fixed point within the interval
for
.
If , two cases arise based on the relative positions of points and . In the first case, when , the point is positioned above . As demonstrated in the previous proof, possesses a fixed point in this case.
In the second case, when , two sub-cases can be considered. If , the conclusion is evidently valid. On the other hand, if , the proof process from case in Theorem 3, indicates possesses a fixed point for this case.
In summary, the order-1 periodic solution of model (
2) always exists. □
Theorems 3 and 4 have confirmed the existence of the order-1 periodic solution in the system (
2), and then we will prove the stability of this periodic solution.
Define
to represent an order-1 periodic solution with a period of
T. The points where this periodic solution intersects with the lines
and line
are denoted as follows:
then it can be calculated.
and
Hence, the Poincaré criterion serves as a valuable tool to analyze the stability of the order-1 periodic solution
in model (
2).
Theorem 5. The solution of model (2) is asymptotically stable when Proof. For case
, considering the proof methodology of Theorem 2,
can be expressed as
Looking back at Theorem 3, it is evident when
,
holds, which implies that the points
and
coincide. Consequently, we obtain
, indicating that
is stable in this case.
When , point lies above , and it is evident that point is also positioned above . Moreover, . Thus, for , we have and . Thus, it holds true that . Similarly, when , we have .
In summary, when Formula (
14) is satisfied, i.e.,
. Regardless of the initial value, the solution of system (
2) will approach solution
, so the solution
is asymptotically stable. Similarly, under the conditions of Theorem 4,
is also asymptotically stable for case
. □
Next, we will investigate the existence of order
k (
) periodic solutions for model (
2).
Theorem 6. In Case , if , then model (2) only possesses an order-1 periodic solution. However, if and , then the system (2) can have either an order-1 periodic solution or an order-2 periodic solution. Proof. If
, according to Theorem 3, a fixed point must exist in
. The trajectory originating from
crosses
, and the point where they intersect is labeled as
. Since any two trajectories cannot intersect, the position of
is lower than that of
T, and
will be mapped to
after undergoing a pulse effect. Since
, the impulsive point
of point
satisfies
, and thus, the position of point
is lower than that of point
. After undergoing another pulse effect, the next impulsive point,
, is obtained, and the position of
is lower than that of
. By analogy, the position of point
is lower than that of point
(
). Therefore, the sequence of impulsive points monotonically decreases and gradually approaches the point
, yielding
Furthermore, since each trajectory does not intersect with others, the system (
2) only possesses an order-1 periodic solution and does not have an order-2 periodic solution.
If
and
, it is evident that the conclusion is valid. If
, The trajectory originating from
crosses
, and the point where they intersect is labeled as
, satisfying the inequality
. The point
is mapped to the point
after undergoing a pulse effect and satisfies the inequality
. Repeat the above steps as the control strategy is implemented multiple times. It can be observed that the position of
is between
and
, and the position of
is between
and
(
). If the series of impulsive points is divided into odd and even series, two possibilities arise.
Based on
and
, the odd series
and even series
are increasing and decreasing, respectively, on the interval
. Both series are also convergent. Furthermore, since all impulsive segments do not intersect each other, either has a fixed point or exhibits a stable period two-point cycle for
. This corresponds to the order-1 or order-2 periodic solutions of the model (
2), respectively. □
According to the aforementioned analysis, we can also conclude that for case
, if
, then model (
2) only possesses an order-1 periodic solution. If
and
, then model (
2) has either an order-1 periodic solution or an order-2 periodic solution.
Now, we will proceed to discuss the conditions under which the order-3 periodic solution exists.
Theorem 7. Let . In case , if and , then the order-3 periodic solution of model (2) exists. Proof. According to Theorem 3, when
, the
possesses a fixed point on
. To prove the model (
2) admits an order-3 periodic solution, it should be demonstrated that there is a
y on interval
makes
and
. According to the requirements stipulated, yield
and
Considering the property of
being continuous, we can say that there is a number
in the interval
that makes
. Hence, model (
2) has an order-3 periodic solution. □
Theorem 7 indicates that a discrete mapping
has a periodic point with a minimum period of 3. According to the Period Three theorem and the Sarkovskii theorem, it is shown that mapping
also has periodic points with a minimum period of any positive integer. The reason is as follows: From the reasoning process of the Sarkovsskii theorem, it can be inferred that there is such a sequence rearranged by all positive integers:
This sequence is called the Sarkovsskii sequence; periodic points will operate in the order of the sequence. The Sarkovsskii theorem states that when
m is sorted before
n in the Sarkovsskii sequence, if there are periodic points with period
m, then there must be periodic points with period
n, and 3 is the first digit of the Sarkovsskii sequence, before all positive integers. According to the Sarkovsskii theorem, if there is a period of 3 in the mapping
during the iteration process, there must be a period point with a period of any positive integer.
This means the existence of order-
k periodic solutions (
) in model (
2). This observation suggests the presence of chaos in model (
2). Compared with reference [
13], our study not only emphasizes the phase diagrams and sets of each equilibrium point but also highlights the existence and stability of periodic solutions under pulse control.
For cases
and
, we have established the existence and stability of the order-1 periodic solution, as well as the conditions for the existence of the order-
k periodic solution (
), in model (
2). However, to enhance our understanding of model dynamics, numerical simulations are performed to validate the accuracy of our theoretical findings.
As shown in
Figure 8, we observe a bifurcation behavior of model (
2) by selecting parameter
p as bifurcation parameter. It is evident that as
p increases, the system undergoes abrupt transitions from an order-
periodic solution to an order-
k periodic solution (
) through period-decreasing bifurcations accompanied by chaotic regions. Furthermore, when
, period-halving bifurcations occur, leading to the emergence of an order-1 periodic solution. This demonstrates the significant influence of the key parameter
p on the dynamics of model (
2). Moreover, it indicates that the growth patterns and quantity trends of pests are dependent on the initial quantity
.
The bifurcation diagram of the parameter
reveals the intricate dynamics of model (
2), as illustrated in
Figure 9. When
, it does not show the solution of model (
2). However, at
, a transcritical bifurcation occurs, giving rise to an order-1 periodic solution. In addition, this solution will become unstable in pace as the parameter
increases further, leading to a period-doubling bifurcation and the subsequent generation of an order-2 periodic solution in model (
2). As there is a continued increase in
, the period-doubling bifurcations drive the system (
2) into a chaotic regime. Notably, the system (
2) demonstrates distinct transitions from an order-
k periodic solution to an order-
periodic solution (where
) through period-adding bifurcations involving chaotic regions. These solutions validate that both populations can coexist in different states, thus further corroborating the conclusions of the above theorems.
Figure 9 clearly illustrates the transition from regular dynamics to chaos, where a period-doubling bifurcation occurs as the bifurcation parameter values increase. The originally stable periodic orbit constantly becomes unstable, and new stable orbits appear, but the emergence of new stable orbits does not mean that the original orbit disappears. It can only be said that the original orbit exists in an unstable form in the new stable state. As the number of iterations increases, a blurry area eventually appears, indicating that the system has entered a completely chaotic state. And
Figure 8 can be seen as the inverse process of the period-doubling bifurcation.
Next, the correctness of the above theorems is verified, and schematic diagrams of the order-1, order-2, order-3, and even order-4 periodic solutions of the system (
2) are presented in conjunction with the bifurcation diagram. Here, we mainly discuss case
. As shown in
Figure 10,
Figure 11,
Figure 12 and
Figure 13. In addition, in order to further understand the population dynamics of prey and predators, we have also provided the phase portrait and time series. This also once again demonstrates the complex, dynamic behavior of the system. Here, we choose the order-2 and order-4 periodic solutions as examples to illustrate.
For the given parameter values, it can be observed from
Figure 11 that the order-2 periodic solution exhibits distinct characteristics. It consists of an active period with short trajectories and a quiescent period with long trajectories. During the long stretches of trajectory where both populations reside, the pest population undergoes a period of rapid growth until an outbreak occurs, as the number of predators is relatively low. However, with the introduction of predators, a quiescent period comes to an end, and the active period begins following a single impulsive effect. The population density of the pest immediately adjusts to
and experiences rapid growth until reaching
on short notice. At this point, control tactics are implemented, causing the quantities of both populations to return to their long stretches of trajectory, thus initiating a new cycle.
Similarly, by varying the values of the key parameters
p and
, we can obtain the order-4 periodic solution of model (
2), as illustrated in
Figure 13. The phase portrait and time series exhibit a more intricate pattern, featuring two quiescent periods and two active periods with long trajectories and short trajectories, respectively. These active and quiescent periods alternate, giving rise to highly complex dynamics. It is evident that random perturbations play a vital role in the occurrence of pest outbreaks. Small changes in the key parameters can directly impact the population dynamics and outbreak patterns of the pests. This further emphasizes the sensitivity of the system to variations in parameter values, highlighting the need for effective control measures to mitigate and manage pest outbreaks.
The above descriptions of bifurcation and chaos all emphasize that the complex dynamics of the model pose challenges to the comprehensive control strategy of pests. This also indicates that analyzing biological problems and explaining biological phenomena and related conclusions from a mathematical perspective is still extremely difficult. In addition to using predatory models to describe and analyze the common biological activity of pest control, other biological activities such as virus outbreaks and vaccination can also be explained using SIR mathematical models. For example, in an improved SIR model in [
31], the author and his team linked the classical SIR model with pulse vaccination, analyzed the dynamic properties of disease-free periodic solution and stable positive T-periodic solution, and described the bifurcation phenomenon and chaos in the pulse model. In addition, there are also models in [
32,
33] that combine mathematics with biology, which reveal that biological systems can be simplified into computable models, thereby helping us understand and predict biological processes.
6. Conclusions
To enhance the efficacy of pest control, impulsive tactics such as releasing natural enemies and using insecticides are employed when the pest density increases to . In light of this, these control strategies dependent on pest status prove more suitable for practical implementation. However, in previous models, the detrimental impacts of pests’ anti-predator behavior on themselves were often disregarded. Recognizing the significance of this aspect, the Lotka-Volterra model incorporates the effect of anti-predator behavior on both parties, and its primary objective is to offer more rational pest control strategies and illustrate the dynamics of pest-natural enemy ecosystems more accurately.
The phase portrait of different equilibria was examined as the parameters
m and
v underwent changes. Subsequently, the relationship between the equilibrium
and the line
was utilized to discuss the precise domain of the phase set and impulsive set in cases
to
. So as to facilitate the analysis of the periodic solution in subsequent studies, these cases where the domain is identical are defined as
, as illustrated in
Table 1. Building upon the threshold control process outlined in the impulsive semi-dynamical system, it was observed that the ordinates of two adjacent impulsive points satisfy
, leading to the establishment of the corresponding Poincaré map.
The stability condition about the boundary order-1 periodic solution
is determined by Theorem 2. For the given parameters, an illustration depicting the global stability of solution
when
is provided. In this scenario, the number of pests will exhibit periodic oscillations, while predators will gradually approach 0. This behavior is depicted in
Figure 5. As the parameters
a and
p are altered,
is no longer stable, leading to the coexistence of pests and natural enemies in a stable focus; see
Figure 6 for details.
Where , we provide the threshold conditions sequentially for the existence of order-k () periodic solutions; this suggests the presence of chaos within the system. Notably, these conditions also indicate that pests and natural enemies can coexist under specific parameter values. Numerical results demonstrate that the system exhibits diverse and intricate dynamic properties, including transcritical bifurcations, period-decreasing bifurcations, period-adding bifurcations, chaotic regions, and more.
One of the highlights of this article is the bifurcation and chaos of the model (
2). When related to other biological models [
14], it can be found that we all express the stability and threshold conditions of the periodic solution and also elaborate on the switching process of the periodic solution under different bifurcations. However, the innovation of this article is to discuss chaos under different sensitivity parameters and provide periodic solutions to various orders of the system in combination with chaos.
Significant development has been achieved in the theoretical exploration of predator-prey systems [
34,
35,
36,
37] over recent years. However, this paper introduces several innovative contributions that set it apart from previous research. The key distinctions are emphasized in the following aspects: (1) Consideration of the reciprocal influence of anti-predator behavior on both populations. By acknowledging the strategic actions taken by pests to minimize predation rates, the dynamics of both populations are affected, establishing a symbiotic relationship between the species; (2) Precise delineation of the domains of the phase set and impulsive set based on distinct threshold conditions. This enables a more accurate understanding of when impulsive control strategies should be employed and how the system behaves accordingly; (3) Provision of multiple bifurcation diagrams for essential parameters. These diagrams provide a comprehensive depiction of the intricate and diverse dynamics exhibited by system (
2), facilitating a deeper comprehension of its behavior across varying parameter values; (4) Illustration of order-
k periodic solutions, specifically focusing on order-2 and order-4 solutions. These solutions shed light on the periodic behavior exhibited by the system, contributing to a comprehensive understanding of its temporal patterns.
Collectively, these distinct characteristics presented in this article contribute to the progress in our comprehension of predator-prey systems and provide valuable insights into the intricate dynamics of model (
2). However, there are still some areas for improvement in the article, such as assuming that harvesting has a meaning. In such a situation, if the harvesting can be either a constant number, a proportional number, or a periodic number of the density of the prey over a given time interval, what can change in the system’s behavior? Additionally, this study did not account for the impact of limited resources, and the utilization of linear pulse control fails to optimize the economic conditions of farmland. Exploring the impact of resource constraints and investigating alternative control strategies to enhance the economic viability of agricultural settings will be key research objectives in the future.