Research Article  Open Access
Zhihua Chen, Qamar Din, Muhammad Rafaqat, Umer Saeed, Muhammad Bilal Ajaz, "DiscreteTime PredatorPrey Interaction with Selective Harvesting and Predator SelfLimitation", Journal of Mathematics, vol. 2020, Article ID 6737098, 13 pages, 2020. https://doi.org/10.1155/2020/6737098
DiscreteTime PredatorPrey Interaction with Selective Harvesting and Predator SelfLimitation
Abstract
Selective harvesting plays an important role on the dynamics of predatorprey interaction. On the other hand, the effect of predator selflimitation contributes remarkably to the stabilization of exploitative interactions. Keeping in view the selective harvesting and predator selflimitation, a discretetime predatorprey model is discussed. Existence of fixed points and their local dynamics is explored for the proposed discretetime model. Explicit principles of Neimark–Sacker bifurcation and perioddoubling bifurcation are used for discussion related to bifurcation analysis in the discretetime predatorprey 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 preypredator 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 preypredator 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 preypredator 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 predatorprey system with a prey refuge and selective harvesting model. Chakraborty et al. [7] investigated global stability, bifurcation, and optimal control for a stage structured preypredator fishery model with harvesting. Xiao et al. [8] discussed coexistence, bifurcation, and limit cycle in a ratiodependent preypredator 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 selflimitation [17], the following predatorprey 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 halfsaturation 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 perioddoubling bifurcation is implemented. We discuss chaos control techniques for the discretetime 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 3dimensional discretetime 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. PeriodDoubling 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, perioddoubling bifurcation occurs at critical value .
Choosing and as bifurcation parameter, the following result provides us parametric conditions under which system (2) undergoes perioddoubling 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 discretetime 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 discretetime model which is obtained through implementation of Euler approximation [29, 30]. For some other methods related to chaos control in discretetime 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 discretetime 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 perioddoubling 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 perioddoubling 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 preypredator interaction are studied under the influence of harvesting. Most of these investigations are dedicated to predatorprey interactions with overlapping generations. In order to see effect of selective harvesting on the dynamics of preypredator interaction with nonoverlapping generation, a discretetime predatorprey model is studied. Our investigation reals that proposed model undergoes perioddoubling and Neimark–Sacker bifurcations around its interior fixed point whenever intrinsic growth rate of predator is taken as bifurcation parameter. Consequently, the discretetime 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.
References
 C. Çelik, “Hopf bifurcation of a ratiodependent predatorprey system with time delay,” Chaos, Solitons & Fractals, vol. 42, no. 3, pp. 1474–1484, 2009. View at: Publisher Site  Google Scholar
 L. Deng, X. Wang, and M. Peng, “Hopf bifurcation analysis for a ratiodependent predatorprey system with two delays and stage structure for the predator,” Applied Mathematics and Computation, vol. 231, pp. 214–230, 2014. View at: Publisher Site  Google Scholar
 A. A. Berryman, “The orgins and evolution of predatorprey theory,” Ecology, vol. 73, no. 5, pp. 1530–1535, 1992. View at: Publisher Site  Google Scholar
 R. P. Gupta and P. Chandra, “Bifurcation analysis of modified LeslieGower predatorprey model with MichaelisMenten type prey harvesting,” Journal of Mathematical Analysis and Applications, vol. 398, no. 1, pp. 278–295, 2013. View at: Publisher Site  Google Scholar
 T. K. Kar and A. Ghorai, “Dynamic behaviour of a delayed predatorprey model with harvesting,” Applied Mathematics and Computation, vol. 217, no. 22, pp. 9085–9104, 2011. View at: Publisher Site  Google Scholar
 M. Peng, Z. Zhang, X. Wang, and X. Liu, “Hopf bifurcation analysis for a delayed predatorprey system with a prey refuge and selective harvesting,” Journal of Applied Analysis and Computation, vol. 8, no. 3, pp. 982–997, 2018. View at: Google Scholar
 K. Chakraborty, S. Jana, and T. K. Kar, “Global dynamics and bifurcation in a stage structured preypredator fishery model with harvesting,” Applied Mathematics and Computation, vol. 218, no. 18, pp. 9271–9290, 2012. View at: Publisher Site  Google Scholar
 D. Xiao, W. Li, and M. Han, “Dynamics in a ratiodependent predatorprey model with predator harvesting,” Journal of Mathematical Analysis and Applications, vol. 324, no. 1, pp. 14–29, 2006. View at: Publisher Site  Google Scholar
 Z. Gui and W. Ge, “The effect of harvesting on a predator–prey system with stage structure,” Ecological Modelling, vol. 187, no. 23, pp. 329–340, 2005. View at: Publisher Site  Google Scholar
 J. Huang, Y. Gong, Y. Gong, and S. Ruan, “Bifurcation analysis in a predatorprey model with constantyield predator harvesting,” Discrete & Continuous Dynamical Systems—B, vol. 18, no. 8, pp. 2101–2121, 2013. View at: Publisher Site  Google Scholar
 S. O. Sowole, A. Ibrahim, D. Sangare, and A. O. Lukman, “Mathematical model for measles disease with control on the susceptible and exposed compartments,” Open Journal of Mathematical Scineces, vol. 4, no. 1, pp. 60–75, 2020. View at: Publisher Site  Google Scholar
 C. M. D. Simarmata, N. Susyanto, N. Susyanto, I. J. Hammadi, and C. Rahmaditya, “A mathematical model of smoking behaviour in Indonesia with densitydependent death rate,” Open Journal of Mathematical Sciences, vol. 4, no. 1, pp. 118–125, 2020. View at: Publisher Site  Google Scholar
 M. Tahir, G. Zaman, G. Zaman et al., “The stability analysis and control transmission of mathematical model for Ebola virus,” Open Journal of Mathematical Analysis, vol. 3, no. 2, pp. 91–102, 2019. View at: Publisher Site  Google Scholar
 A. Bukkuri, “A mathematical model of the effects of ascorbic acid on the onset of neurodegenerative diseases,” Open Journal of Mathematical Sciences, vol. 3, no. 1, pp. 300–309, 2019. View at: Publisher Site  Google Scholar
 M. N. Hasan, M. H. A. Biswas, M. H. A. Biswas, and M. S. Uddin, “A mathematical model for fish management in the Sundarbans ecosystem,” Open Journal of Mathematical Analysis, vol. 3, no. 2, pp. 42–49, 2019. View at: Publisher Site  Google Scholar
 A. Bukkuri, “A mathematical model showing the potential of vitamin c to boost the innate immune response,” Open Journal of Mathematical Scineces, vol. 3, no. 3, pp. 245–255, 2019. View at: Publisher Site  Google Scholar
 T. K. Kar, U. K. Pahari, and K. S. Chaudhuri, “Conservation of a preypredator fishery with predator self limitation based on continuous fishing effort,” Journal of Applied Mathematics and Computing, vol. 19, no. 12, pp. 311–326, 2005. View at: Publisher Site  Google Scholar
 L. Dai, Nonlinear Dynamics of Piecewise Constant Systems and Implementation of Piecewise Constant Arguments, World Scientific Publishing Co. Pte. Ltd., Singapore, 2008.
 E. Camouzis and G. Ladas, Dynamics of Third–Order Rational Difference Equations with Open Problems and Conjectures, Chapman and Hall/CRC, New York, NY, USA, 2008.
 Q. Din, A. A. Elsadany, and H. Khalil, “NeimarkSacker bifurcation and chaos control in a fractionalorder plantherbivore model,” Discrete Dynamics in Nature and Society, vol. 2017, Article ID 6312964, 15 pages, 2017. View at: Publisher Site  Google Scholar
 A. E. Matouk, A. A. Elsadany, E. Ahmed, and H. N. Agiza, “Dynamical behavior of fractionalorder HastingsPowell food chain model and its discretization,” Communications in Nonlinear Science and Numerical Simulation, vol. 27, no. 1–3, pp. 153–167, 2015. View at: Publisher Site  Google Scholar
 Q. Din and W. Ishaque, “Bifurcation analysis and chaos control in discretetime ecoepidemiological models of pelicans at risk in the Salton Sea,” International Journal of Dynamics and Control, vol. 8, no. 1, pp. 132–148, 2019. View at: Publisher Site  Google Scholar
 W. Ishaque, Q. Din, M. Taj, and M. A. Iqbal, “Bifurcation and chaos control in a discretetime predatorprey model with nonlinear saturated incidence rate and parasite interaction,” Advances in Difference Equations, vol. 2019, no. 1, 2019. View at: Publisher Site  Google Scholar
 S. Qin, J. Zhang, W. Du, and J. Yu, “Neimark–Sacker bifurcation in a new three–dimensional discrete chaotic system,” ICIC–EL, vol. 10, no. 4, pp. 1–7, 2016. View at: Google Scholar
 G. Wen, “Criterion to identify Hopf bifurcations in maps of arbitrary dimension,” Physical Review E, vol. 72, no. 2, Article ID 026201, 2005. View at: Publisher Site  Google Scholar
 G. Wen, S. Chen, and Q. Jin, “A new criterion of perioddoubling bifurcation in maps and its application to an inertial impact shaker,” Journal of Sound and Vibration, vol. 311, no. 12, pp. 212–223, 2008. View at: Publisher Site  Google Scholar
 E. Ott, C. Grebogi, and J. A. Yorke, “Controlling chaos,” Physical Review Letters, vol. 64, no. 11, pp. 1196–1199, 1990. View at: Publisher Site  Google Scholar
 S. Lynch, Dynamical Systems with Applications Using Mathematica, Birkhauser, Boston, MA, USA, 2007.
 Q. Din, “Bifurcation analysis and chaos control in discretetime glycolysis models,” Journal of Mathematical Chemistry, vol. 56, no. 3, pp. 904–931, 2018. View at: Publisher Site  Google Scholar
 Q. Din, T. Donchev, and D. Kolev, “Stability, bifurcation analysis and chaos control in chlorine dioxideiodinemalonic acid reaction,” MATCH Communications in Mathematical and in Computer Chemistry, vol. 79, no. 3, pp. 577–606, 2018. View at: Google Scholar
 Q. Din, “Complexity and chaos control in a discretetime preypredator model,” Communications in Nonlinear Science and Numerical Simulation, vol. 49, pp. 113–134, 2017. View at: Publisher Site  Google Scholar
 Q. Din, “NeimarkSacker bifurcation and chaos control in HassellVarley model,” Journal of Difference Equations and Applications, vol. 23, no. 4, pp. 741–762, 2017. View at: Publisher Site  Google Scholar
 Q. Din, Ö. A. Gümüş, and H. Khalil, “Neimarksacker bifurcation and chaotic behaviour of a modified hostparasitoid model,” Zeitschrift für Naturforschung A, vol. 72, no. 1, pp. 25–37, 2017. View at: Publisher Site  Google Scholar
 Q. Din, “Controlling chaos in a discretetime preypredator model with Allee effects,” International Journal of Dynamics and Control, vol. 6, no. 2, pp. 858–872, 2018. View at: Publisher Site  Google Scholar
 Q. Din, “Qualitative analysis and chaos control in a densitydependent hostparasitoid system,” International Journal of Dynamics and Control, vol. 6, no. 2, pp. 778–798, 2018. View at: Publisher Site  Google Scholar
 Q. Din, “A novel chaos control strategy for discretetime Brusselator models,” Journal of Mathematical Chemistry, vol. 56, no. 10, pp. 3045–3075, 2018. View at: Publisher Site  Google Scholar
 Q. Din, “Stability, bifurcation analysis and chaos control for a predatorprey system,” Journal of Vibration and Control, vol. 25, no. 3, pp. 612–626, 2019. View at: Publisher Site  Google Scholar
 Q. Din and M. A. Iqbal, “Bifurcation analysis and chaos control for a discretetime enzyme model,” Zeitschrift für Naturforschung A, vol. 74, no. 1, pp. 1–14, 2019. View at: Publisher Site  Google Scholar
 Q. Din and M. Hussain, “Controlling chaos and NeimarkSacker bifurcation in a hostparasitoid model,” Asian Journal of Control, vol. 21, no. 3, pp. 1202–1215, 2019. View at: Publisher Site  Google Scholar
 M. A. Abbasi and Q. Din, “Under the influence of crowding effects: stability, bifurcation and chaos control for a discretetime predatorprey model,” International Journal of Biomathematics, vol. 12, no. 4, Article ID 1950044, 2019. View at: Publisher Site  Google Scholar
 X. S. Luo, G. Chen, B. Hong Wang, and J. Qing Fang, “Hybrid control of perioddoubling bifurcation and chaos in discrete nonlinear dynamical systems,” Chaos, Solitons & Fractals, vol. 18, no. 4, pp. 775–783, 2003. View at: Publisher Site  Google Scholar
 D. Jana and G. P. Samanta, “Significance of prey harvesting in preypredator system in discrete time scale using interval parameters,” International Journal of Ecological Economics & Statistics, vol. 39, no. 2, pp. 46–60, 2018. View at: Google Scholar
 D. Jana, S. Dutta, and G. P. Samanta, “Interplay between reproduction and age selective harvesting: a case study of Hilsa (Tenualosa ilisha) fish at Sundarban estuary of northern Bay of Bengal, India,” International Journal of Biomathematics, vol. 12, no. 2, Article ID 1950023, 2019. View at: Publisher Site  Google Scholar
 S. Sharma and G. P. Samanta, “Dynamical behaviour of ageselective harvesting of a preypredator system,” International Journal of Dynamics and Control, vol. 6, no. 2, pp. 550–560, 2018. View at: Publisher Site  Google Scholar
 D. Pal, G. P. Samanta, and G. S. Mahapatra, “Selective harvesting of two competing fish species in the presence of toxicity with time delay,” Applied Mathematics and Computation, vol. 313, pp. 74–93, 2017. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Zhihua Chen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.