Journal of Function Spaces
Volume 2019 (2019), Article ID 3948621, 13 pages
https://doi.org/10.1155/2019/3948621

Research Article

Integrability and Multiple Limit Cycles in a Predator-Prey System with Fear Effect

Xinli Li¹ and Ming Zhang²

¹Logistic School, Linyi University, Linyi 276000, China
²College of Information Science and Engineering, Linyi University, Linyi 276000, China

Correspondence should be addressed to Xinli Li; lixinli0539@sina.com and Ming Zhang; zhangming@lyu.edu.cn

Received 22 February 2019; Accepted 11 June 2019; Published 27 June 2019

Academic Editor: Gen Q. Xu

Copyright © 2019 Xinli Li and Ming Zhang. 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.

Abstract.

In this paper, a class of predator-prey systems with fear effect is investigated, the integrability conditions of the origin and the positive equilibrium are obtained, and the fact that three limit cycles can be bifurcated from the positive equilibrium is proved, so bistable phenomenon can occur for this system.

1. Introduction

Predator-prey systems which were considered to be among the most important models in ecological systems have been investigated intensively and extensively since they were proposed; there have been many good results for these systems. The first model for predator-prey problem was Lotka-Volterra model; this system was improved gradually year by year. Many new models were constructed with the development of biological system.

Functional response which was introduced by Holling has a great effect on predator-prey system. In 1965, the famous Holling type II functional response of predators [1] was proposed by Holling. As application, the population dynamics of predator-prey systems with the Holling type II functional response have been studied by many authors; for such a model, the result about the existence of a unique stable limit cycle has also been proved. Considering the reality, many more complicated functional responses were suggested one by one, and predator-prey systems have also been modelled gradually. For example, some monotone response functions dependent on prey have been investigated by May [2], Seo and DeAngelis [3], and Huang et al. [4]. Furthermore, Zhu et al. [5], Ruan and Xiao [6], Freedman and Wolkowicz [7], and Wolkowicz [8] studied some nonmonotone response functions dependent on prey. Some interesting bifurcation phenomenon was found for those systems. Functional responses dependent on prey only have been considered intensively; some functional responses dependent on both prey and predators were also considered in recent years, and the famous functional responses dependent on both prey and predators are the Beddington-DeAngelis functional responses, see [9–13]), and ratio dependent functional response [14, 15].

Although the functional responses become more and more complicated for predator-prey systems, they can only reflect the direct killing. There is profound lack of the functional responses and many realities which cannot be concluded in those systems. For example, in 2016, Zou and Wang considered a predator-prey model incorporating the cost of fear to explore the impact that fear can have on population dynamics in predator-prey systems, dynamical behaviors of this system was discussed, and a limit cycle which can be bifurcated from the positive equilibrium by Hopf bifurcation was proved [16]. But the number of limit cycles is not considered in their paper. In fact there can be many limit cycles in predator-prey system which can lead to multiple stable phenomenons. For example, the system of Leslie-type predator-prey schemes with a nonmonotonic functional response and Allee effect allows the existence of three limit cycles under certain conditions over the parameters; by Hopf bifurcation [17], the first two cycles are infinitesimal ones generated, and then the third one can be bifurcated from a homoclinic loop. The singular perturbation theory which was developed by Dumortier and Roussarie can be used to study the canard phenomenon. The authors studied the canard phenomenon for predator-prey systems with response functions of Holling types [18] by using this method. For a Gause type predator-prey system with Holling type III functional response and Allee effect on prey, multiple limit cycles were considered in [19] by multiple Hopf bifurcations. In [20], a criterion to calculate the multiplicity of a multiple focus was developed by the authors for general predator-prey systems. Recently, the number of limit cycles in a Leslie-Gower-type predator-prey system with weak Allee effect on prey was considered in [21].

The number of limit cycles is closely related not only to the Hilbert 16th problem but also to bistable phenomenon or multiple stable phenomenon. For instance, in [21], the existence conditions for which three limit cycles appear and surround different equilibrium points were proved. Above all, the three limit cycles have different stability which means the phenomenon of multistability. For finitely smooth planar autonomous differential systems, in [22], the author considered the bifurcation theory. Theory of rotated equations is discussed and applied to a population model in [23]. For a family of real planar polynomial ODE systems depending on parameters, the problem of how to determine the systems in the family which become time-reversible after some affine transformation was considered in [24]. Recently, an improvement on the number of limit cycles bifurcating from a nondegenerate center of homogeneous polynomial systems was obtained in [25] by normal form method. Some new perturbation methods were given in [26–28]. The integrability of equivalent systems was also considered in [29, 30]. Therefore, in this paper, we reinvestigated the predator-prey system with fear effect which is written as follows:

By transformation

and time scale

system (1) can be reduced to

And system (4) has three singular points

where

The classification of the three singular points and stability have been discussed in [16], but the integrability and center problem remained. The rest of the paper is organized as follows. In the next section, we prove that the origin of system (4) is integral. In Section 3, the first four Lyapunov constants will be computed; bifurcation of limit cycles and center conditions of (4) are investigated. Section 4 is devoted to discussing the kind of the boundary equilibrium. At last, a conclusion is drawn to illustrate that bistable phenomenon can occur in this system.

2. Integrability Conditions of

For the integrability of origin of system (4), the authors have proved that it is a saddle. In this section, we will compute the saddle values; furthermore, the integrability conditions can be obtained at the same time. We only consider a special case; namely, let ; we study the integrability conditions theoretically. The system (4) becomes

Theorem 1. The first five saddle values at the origin of system (7) are

If

(1)

(2)

(3)

If ,

(4)

(5)

Here are given in the appendix.

Thus, we can get the following integrability conditions of system (7) easily.

Theorem 2. The origin of system (7) is integrability if and only if one of the following conditions holds:

Proof. When condition (I) holds, system (7) becomes

which has an integrating factor

When condition (II) holds, system (7) can be rewritten as

which admits an integrating factor

When condition (III) holds, system (7) can be changed into

which admits an integrating factor

When condition (IV) holds, system (7) can be transformed into

which has an integrating factor

3. Hopf Bifurcation at

According to results in [16], there exists a positive singular point . Let ; system (4) can be transformed into

Here

Suppose there are two characteristic values ; then

namely,

and we have

Without loss of generality, let

and the singular point becomes . From the first equation , we get

Meanwhile, the second equation means that

which yields

Accordingly, system (7) can be changed into

The Jacobin matrix at the of system (34) can be written as follows:

and by linear transformation

system can be shifted to

The Jacobin matrix at the origin of system (37) can be written as follows:

and, with the help of computer algebraic system Maple, we can compute the first five Lyapunov constants for system (37) by using formal series method which was given in the following.

Theorem 3. For system (37), the first five Lyapunov constants at the origin can be given as follows:

where are given in the appendix.

Theorem 3 yields the following.

Proposition 4. For system (37), the first five Lyapunov constants at the origin are zero if and only if one of the following conditions holds:

Furthermore, we have the following.

Theorem 5. The origin of system (37) is a center if and only if one of the conditions in Proposition 4 holds.

Proof. When condition (40) holds, system (37) becomes

which has four invariant algebraic curves

and admits an integrating factor

When condition (41) holds, system is simplified into

which has three invariant algebraic curves and admits an integrating factor

Two examples of phase plane of system (42) when condition (40) or (41) is satisfied are given in Figures 1 and 2.

Figure 1: The phase plane of system (42) when .

Figure 2: The phase plane of system (45) when .

Based on Theorem 3, we can conclude the following.

Theorem 6. The origin of system (37) is fifth order weak foci if and only if one of the following conditions holds:

Proof. Let , and we obtain

submitting into , we have

Solve

and , and we can obtain the results.

Considering the reality, all parameters are positive in a predator-prey system, so we cannot get five limit cycles for this system.

In fact, means that

and, submitting it into , we get Figure 3 which shows that cannot be larger than zero simultaneously.

Figure 3: with .

We can adjust the parameters to satisfy , and all parameters are positive. For example,

Then three limit cycles can be bifurcated from the origin; namely, bistable phenomenon can occur for this system; this is not discussed in previous discussion.

4. Boundary Singular Point

In this section, we consider the boundary singular point of system (4). The boundary singular point is

and the Jacobin matrix at of system (4) is

It is easy to testify that when , system (4) has two eigenvalues with different signs which means is a saddle. When , system (4) has two eigenvalues with negative signs which means is a stable node. When , namely, , system (4) becomes and one eigenvalue is zero. System (4) becomes

Shifting the of system (4) to origin, system (4) becomes

the Jacobin matrix at origin system (56) is

By transformation

and time scale

system (56) becomes

Based on the theorem in [31], we have the following.

Theorem 7. The origin of system (56) is a saddle-node singular point.

Proof. From , we get a solution by implicit function theorem; submitting it into , we have

Based on the theorem in [31], the origin of system (56) is a saddle-node singular point.

An numerical simulation of phase portrait of system (56) when is given in Figure 4 which shows that the origin of system (56) is a saddle-node singular point.

Figure 4: When , the is a saddle node.

5. Conclusion

In this paper, a class of predator-prey systems with fear effect is investigated; the integrability conditions of and positive equilibrium are obtained and proved by finding their integrating factors. Furthermore, we compute the first five focal values at the positive equilibrium; considering the reality, when we choose all parameters to be positive, the fact that at least three limit cycles can be bifurcated from the positive equilibrium is proved. Consequently, bistable phenomenon can occur for this system.

Appendix

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Both authors contributed equally to the writing of this paper and read and approved the final manuscript.

Acknowledgments

We are thankful for the support of Shandong Provincial Natural Science Foundation, China (No. ZR2014FL017).

References

C. S. Holling, “The functional response of predators to prey density and itsrole inmimicry and population regulation,” Memoirs of the Entomological Society of Canada, vol. 97, no. S45, pp. 5–60, 1965.

R. M. May, “Limit cycles in predator-prey communities,” Science, vol. 177, no. 4052, pp. 900–902, 1972.

G. Seo and D. L. DeAngelis, “A predator–prey model with a holling type I functional response including a predator mutual interference,” Journal of Nonlinear Science, vol. 21, no. 6, pp. 811–833, 2011.

J. Huang, S. Ruan, and J. Song, “Bifurcations in a predator-prey system of Leslie type with generalized Holling type III functional response,” Journal of Differential Equations, vol. 257, no. 6, pp. 1721–1752, 2014.

H. Zhu, S. A. Campbell, and G. S. Wolkowicz, “Bifurcation analysis of a predator-prey system with nonmonotonic functional response,” SIAM Journal on Applied Mathematics, vol. 63, no. 2, pp. 636–682, 2003.

S. Ruan and D. Xiao, “Global analysis in a predator-prey system with nonmonotonic functional response,” SIAM Journal on Applied Mathematics, vol. 61, no. 4, pp. 1445–1472, 2001.

H. I. Freedman and G. S. K . Wolkowicz, “Predator-prey systems with group defence: the paradox of enrichment revisited,” Bulletin of Mathematical Biology, vol. 48, no. 5-6, pp. 493–508, 1986.

G. S. K. Wolkowicz, “Bifurcation analysis of a predator-prey system involving group defence,” SIAM Journal on Applied Mathematics, vol. 48, no. 3, pp. 592–606, 1988.

R. S. Cantrell and C. Cosner, “On the dynamics of predator-prey models with the Beddington-DeAngelis functional response,” Journal of Mathematical Analysis and Applications, vol. 257, no. 1, pp. 206–222, 2001.

J. R. Beddington, “Mutual interference between parasites or predat ors and its effect on searching efficiency,” Journal of Animal Ecology, vol. 44, no. 1, pp. 331–340, 1975.

D. L. DeAngelis, A. H. Goldstein, and R. V. O'Neill, “A model for trophic interaction,” Ecology, vol. 56, no. 4, pp. 881–892, 1975.

T. W. Hwang, “Global analysis of the predator-prey system with Beddington-DeAngelis functional response,” Journal of Mathematical Analysis and Applications, vol. 281, no. 1, pp. 395–401, 2003.

T. W. Hwang, “Uniqueness of limit cycles of the predator-prey system with Beddington-DeAngelis functional response,” Journal of Mathematical Analysis and Applications, vol. 290, no. 1, pp. 113–122, 2004.

Y. Song and X. Zou, “Bifurcation analysis of a diffusive ratio-dependent predator–prey model,” Nonlinear Dynamics, vol. 78, no. 1, pp. 49–70, 2014.

Y. Song and X. Zou, “Spatiotemporal dynamics in a diffusive ratio-dependent predator-prey model near a Hopf-Turing bifurcation point,” Computers & Mathematics with Applications, vol. 67, no. 10, pp. 1978–1997, 2014.

X. Wang, L. Zanette, and X. Zou, “Modelling the fear effect in predator–prey interactions,” Journal of Mathematical Biology, vol. 73, no. 5, pp. 1179–1204, 2016.

P. Aguirre, E. González-Olivares, and E. E. Sáez, “Three limit cycles in a Leslie-Gower predator-prey model with additive Allee effect,” SIAM Journal on Applied Mathematics, vol. 69, no. 5, pp. 1244–1262, 2009.

C. Li and H. Zhu, “Canard cycles for predator-prey systems with Holling types of functional response,” Journal of Differential Equations, vol. 254, no. 2, pp. 879–910, 2013.

E. González-Olivares and A. Rojas-Palma, “Multiple limit cycles in a Gause type predator-prey model with Holling type III functional response and Allee effect on prey,” Bulletin of Mathematical Biology, vol. 73, no. 6, pp. 1378–1397, 2011.

D. Xiao and H. Zhu, “Multiple focus and hopf bifurcations in a predator-prey system with nonmonotonic functional response,” SIAM Journal on Applied Mathematics, vol. 66, no. 3, pp. 802–819, 2006.

E. González-Olivares, A. Rojas-Palma, and B. González-Yañez, “Multiple limit cycles in a Leslie–Gower-type predator–prey model considering weak Allee effect on prey,” Nonlinear Analysis: Modelling and Control, vol. 22, no. 3, pp. 347–365, 2017.

M. Han, L. Sheng, and X. Zhang, “Bifurcation theory for finitely smooth planar autonomous differential systems,” Journal of Differential Equations, vol. 264, no. 5, pp. 3596–3618, 2018.

M. Han, X. Hou, L. Sheng, and C. Wang, “Theory of rotated equations and applications to a population model,” Discrete and Continuous Dynamical Systems - Series A, vol. 38, no. 4, pp. 2171–2185, 2018.

M. Han, T. Petek, and V. G. Romanovski, “Reversibility in polynomial systems of ODE's,” Applied Mathematics and Computation, vol. 338, pp. 55–71, 2018.

P. Yu, M. Han, and J. Li, “An improvement on the number of limit cycles bifurcating from a nondegenerate center of homogeneous polynomial systems,” International Journal of Bifurcation and Chaos, vol. 28, no. 6, Article ID 1850078, 2018.

Y. Liu and F. Li, “Double bifurcation of nilpotent focus,” International Journal of Bifurcation and Chaos, vol. 25, no. 3, Article ID 1550036, p. 10, 2015.

F. Li and M. Wang, “Bifurcation of limit cycles in a quintic system with ten parameters,” Nonlinear Dynamics, vol. 71, no. 1-2, pp. 213–222, 2013.

T. Liu, Y. Liu, and F. Li, “A kind of bifurcation of limit cycles from a nilpotent critical point,” Journal of Applied Analysis and Computation, vol. 8, no. 1, pp. 10–18, 2018.

Y. Yan, Y. Pan, F. Lu, and Z. Zhou, “On the integrability and equivalence of the abel equation and some polynomial equations,” Journal of Nonlinear Modeling and Analysis, vol. 1, no. 2, pp. 207–220, 2019.

F. Li, Y. Liu, Y. Liu, and P. Yu, “Bi-center problem and bifurcation of limit cycles from nilpotent singular points in Z2-equivariant cubic vector fields,” Journal of Differential Equations, vol. 265, no. 10, pp. 4965–4992, 2018.

Z. Zhang, T. Ding, W. Huang, and Z. Dong, The Qualitative Theory of Differential Equations, Science Press, Bejing, China, 1997.