#### Abstract

Selective harvesting plays an important role on the dynamics of predator-prey interaction. On the other hand, the effect of predator self-limitation contributes remarkably to the stabilization of exploitative interactions. Keeping in view the selective harvesting and predator self-limitation, a discrete-time predator-prey model is discussed. Existence of fixed points and their local dynamics is explored for the proposed discrete-time model. Explicit principles of Neimark–Sacker bifurcation and period-doubling bifurcation are used for discussion related to bifurcation analysis in the discrete-time predator-prey system. The control of chaotic behavior is discussed with the help of methods related to state feedback control and parameter perturbation. At the end, some numerical examples are presented for verification and illustration of theoretical findings.

#### 1. Introduction

In recent years, the prey-predator system as the fundamental structure in dynamics has got much attention from researchers [1, 2]. In these studies, there are numerous components that affect dynamical properties of the prey-predator systems. Through mathematical modeling, study of a dynamical behavior among prey and predator has a lot of importance to get general perspicacity during the last few years [3]. The prey-predator interaction with selective harvesting has been investigated by many authors in differential and difference equations.

In addition, the harvesting has a significant impact on the dynamics of the system [4, 5]. In the harvesting system, the main goal is to establish how much we can harvest without altering dangerously the harvested population. Peng et al. [6] examined local stability and Hopf bifurcation in a delayed predator-prey system with a prey refuge and selective harvesting model. Chakraborty et al. [7] investigated global stability, bifurcation, and optimal control for a stage structured prey-predator fishery model with harvesting. Xiao et al. [8] discussed coexistence, bifurcation, and limit cycle in a ratio-dependent prey-predator interaction with predator harvesting model. Similarly, for some other discussions associated with qualitative behavior of selective harvesting models, we refer the interested reader to [9–16].

Taking into account the selective harvesting and predator self-limitation [17], the following predator-prey model was proposed:where , , and represent population densities for prey population, predator population, and harvesting effort at time , respectively. Furthermore, it is supposed that prey population grows logistically with intrinsic growth rate and is used for carrying capacity of environment. The availability of prey for predation is taken into account by implementing Holling type II functional response in which denotes maximum harvesting rate and is the half-saturation constant. Moreover, represents the catchability coefficient, is the intrinsic growth rate of predator, and is quantity of prey required to support one consumer at equilibrium when equals . On the other hand, is a parameter for stiffness which is used to measure the intensity of reaction between the perceived rent and the effort, and are the price for unit catch and the cost of unit effort, respectively, and represents tax on per unit biomass of the landed fish.

In order to investigate more richer and chaotic behavior of (1), we implement piecewise constant argument method [18] to system (1) as follows:

In this paper, qualitative behavior of model (2) is discussed. Remaining discourse of this paper is summarized as follows. Existence of equilibria for system (2) and local stability analysis for these equilibria are presented in Section 2. An explicit criterion related to emergence of Neimark–Sacker bifurcation in system (2) is discussed in Section 3. In Section 4, an explicit criterion related to occurrence of period-doubling bifurcation is implemented. We discuss chaos control techniques for the discrete-time model under the influence of bifurcations in Section 5. Finally, in Section 6, numerical simulations are provided to illustrate all theoretical discussions.

#### 2. Existence of Equilibria and Stability

The fixed points of system (2) satisfy the following algebraic system:

Simple computation yields the following fixed points of system (2):and an interior fixed point is given bywhere

Furthermore, existence conditions for these fixed points of system (2) are summarized in Table 1.

For , , , , , , , , , , and , the nullclines for positive equilibrium are depicted in Figure 1. In Figure 1, -nullclines, -nullclines, and -nullclines are represented by red, blue, and black portions, respectively, and intersection of these nullclines is .

Now, we state the following lemma related to local asymptotic stability of 3-dimensional discrete-time systems.

Lemma 1 (see [19]). * Assume that , , and are real constants. Then, the necessary and sufficient conditions for all roots of the equationto lie inside the unit disk are*

First, we see the dynamics of fixed point and one can easily observe that this equilibrium point is unstable. Indeed, the variational matrix of system (2) at is given as follows:

Here, , , and are eigenvalues of and always lies outside the unit disk. Furthermore, is a saddle point if or . On the other hand, is a source if and . Assume that and ; then, exists and variational matrix at this fixed point is given as follows:

Now, characteristic polynomial for is given by

Then, obviously is an eigenvalue of , so is also unstable. On the other hand, is a saddle point if or , and it is a source if and only if the following conditions are satisfied:

For fixed point , the variational matrix of system (2) evaluated at is given bywhere

Moreover, characteristic polynomial of is computed as follows:

The following lemma gives dynamics of system (2) with respect to its fixed point .

Lemma 2. *The following statements hold true for equilibrium of system (2):*(i)*Fixed point of system (2) is a sink if and only if**(ii)**Fixed point of system (2) is a saddle point if and only if*(iii)*Fixed point of system (2) is a source if and only if*

In order to illustrate Lemma 2 for particular parametric values, the topological classification of is depicted in Figure 2 in -plane. Moreover, in Figure 2, the blue region represents sink, the red region represents source, and three remaining white regions represent saddle. On the other hand, other parameters are taken as , , , , , , , and .

Assume that and ; then, variational matrix of (2) evaluated at unique positive equilibrium point is computed as follows:

Moreover, the characteristic equation of is given bywhere

Lemma 3. *Suppose that and ; then, unique positive equilibrium point of (2) is a sink if and only ifwhere , , and are listed in (21).*

Stability region of positive equilibrium point is depicted in Figure 3.

#### 3. Neimark–Sacker Bifurcation

This section is related to emergence of Hopf bifurcation about positive fixed of system (2). For this, we apply a direct criterion for occurrence of Hopf bifurcation without finding the eigenvalues of the system under consideration (cf. [20–24]). Generally, the explicit criterion for occurrence of Hopf bifurcation in -dimensional discrete system is given as follows.

Lemma 4 (see [25]). *Consider the following -dimensional discrete system:where represents some bifurcation parameter. Furthermore, take into account the following characteristic polynomial for the Jacobian matrix evaluated about fixed point of map :where , , in which denotes control parameter or some other parameter to be determined. Take into account a sequence of determinants of the form such that , andwhere*

Moreover, it is supposed that the following conditions hold true: Eigenvalue criterion: , , , , , for , 2 (or 1), when is odd (or even, respectively). Transversality criterion: . Nonresonance or resonance criterion: , or resonance condition , where , and ; then, Neimark–Sacker bifurcation occurs at critical value .

Choosing , the following result gives conditions under which (2) undergoes Neimark–Sacker bifurcation when is taken as bifurcation parameter.

Lemma 5. *The positive fixed point of system (2) undergoes Neimark–Sacker bifurcation at if the following conditions hold true:where , , and are given in (21), and can be any possible real root of .*

*Proof. *Choosing and as bifurcation parameter, from Lemma 4 and characteristic equation (20), one has

#### 4. Period-Doubling Bifurcation

This section is related to emergence of flip bifurcation about positive fixed of system (2). For this, we apply a direct criterion for occurrence of flip bifurcation without finding the eigenvalues of the system under consideration. Generally, the explicit criterion for occurrence of flip bifurcation in -dimensional discrete system is given as follows.

Lemma 6 (see [26]). *Consider the following -dimensional discrete system:where and represents some bifurcation parameter. Moreover, assume that conditions (24)–(26) of Lemma 4 are also satisfied. Moreover, it is supposed that the following conditions hold true:* * Eigenvalue criterion: , , , , , 1 (or 2), when is even (or odd, respectively).* * Transversality criterion: , where denotes derivative of at ; then, period-doubling bifurcation occurs at critical value .*

Choosing and as bifurcation parameter, the following result provides us parametric conditions under which system (2) undergoes period-doubling bifurcation whenever is taken as bifurcation parameter.

Lemma 7. *The unique positive fixed point of system (2) undergoes flip bifurcation at if the following conditions hold:where , , and are given in (21), is derivative of at , and can be any possible real root of algebraic equation .*

#### 5. Chaos Control

This section is dedicated to implementation of chaos control methods to system (2). First we apply a chaos control method based on state feedback control methodology and this method is known as OGY method (that is, Ott–Grebogi–Yorke method) [27]. The OGY control method is one of the pioneer methods for controlling chaos in discrete-time systems. On the other hand, it has some deficiencies due to some limitations [28]. In particular, the OGY method may not be applicable to a discrete-time model which is obtained through implementation of Euler approximation [29, 30]. For some other methods related to chaos control in discrete-time models, we refer to [29–40].

For the application of the OGY method to system (2), the corresponding control system is written as follows:where , , and . Moreover, is some nominal value of which is located in bifurcating or chaotic region and , , and are control parameters. Variational matrix of (31) evaluated at is given by

This variational matrix is similar to Jacobian matrix of original system (2) at its positive fixed point . Since this variational matrix is independent of control parameters , , and , system (31) is not controllable via the OGY method.

Next, we introduce another chaos control method based on random state feedback control strategy as follows:where , and . Moreover, , , and are control parameters. We compute variational matrix of (33) about its positive fixed point as follows:

Keeping in view this variational matrix of control system (33), its characteristic equation is calculated as follows:where

Then, it is easy to see that system (33) is controllable as long as its positive fixed point is a sink. Consequently, one has the following result.

Lemma 8. *The unique fixed point of system (33) is a sink if and only if*

At the end of this section, we present a simple and effective method for larger classes of discrete-time systems. This chaos control method is based on state feedback control and parameter perturbation. We call this method as the hybrid chaos control method [41]. For application of the hybrid control method to system (2), the corresponding control system can be written as follows:where is control parameter. Meanwhile, the Jacobian matrix for system (38) at its positive fixed point is given by

So, the characteristic polynomial of is given bywhere

Then, it is easy to see that system (38) is controllable as long as its positive fixed point is a sink. Consequently, one has the following result.

Lemma 9. *The unique fixed point of system (38) is a sink if and only if*

#### 6. Numerical Simulation

First of all, we choose , , , , , , , , , , and bifurcation parameter . With the variation in bifurcation parameter in interval , system (2) undergoes period-doubling bifurcation at . Moreover, at , , , , , , , , , , and , the unique positive fixed point of (2) is , and characteristic polynomial about this steady state for variational matrix of system (2) is given by

The roots of this characteristic equation are , and . Therefore, existence of period-doubling bifurcation with the help of criterion of eigenvalues is justified. Next, all axioms of Lemma 7 at are satisfied as follows:

The bifurcation diagrams and maximum Lyapunov exponents (MLE) are depicted in Figure 4.

**(a)**

**(b)**

**(c)**

**(d)**

Furthermore, in order to observe the controllability of system (33), we choose , , , , , , , , , , and . For these parametric values, the controllable region (green) for system (33) is depicted in Figure 5. Secondly, for similar parametric values, one can see the effectiveness of hybrid control strategy (38). Then, some simple calculation yields that system (38) is controllable for .

Finally, for verification of Neimark–Sacker bifurcation, we choose , , , , , , , , , , and . Then, system (2) encounters Neimark–Sacker bifurcation and . Consequently, at , , , , , , , , , , and , system (2) has nonhyperbolic equilibrium since for these parametric values, characteristic polynomial of variation matrix of system (2) is given bywith roots and such that . Moreover, at , the conditions of Lemma 5 are fulfilled as follows:

We consider the equation , and one has . Thus, nonresonance condition is also satisfied. Bifurcation diagrams and MLE are shown in Figure 6. Moreover, in order to check the controllability of system (33), we choose , , , , , , , , , , and . For these parametric values, the controllable region (green) for system (33) is depicted in Figure 7. Secondly, for similar parametric values, one can see the effectiveness of hybrid control strategy (38). Then, some simple calculation yields that system (38) is controllable for .

**(a)**

**(b)**

**(c)**

**(d)**

#### 7. Concluding Remarks

There are abundant studies dedicated to the effects of harvesting on population growth. In particular, various classes of prey-predator interaction are studied under the influence of harvesting. Most of these investigations are dedicated to predator-prey interactions with overlapping generations. In order to see effect of selective harvesting on the dynamics of prey-predator interaction with nonoverlapping generation, a discrete-time predator-prey model is studied. Our investigation reals that proposed model undergoes period-doubling and Neimark–Sacker bifurcations around its interior fixed point whenever intrinsic growth rate of predator is taken as bifurcation parameter. Consequently, the discrete-time model has rich dynamical behavior compared to its continuous counterpart. Furthermore, chaotic and bifurcating behavior of the model is controlled with implementation of various chaos control strategies [42–45].

#### Data Availability

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

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Authors’ Contributions

All authors contributed equally to this study.

#### Acknowledgments

This work was supported by the Chinese National Natural Science Foundation under Grant 61876047.