Table of Contents Author Guidelines Submit a Manuscript
Mathematical Problems in Engineering
Volume 2009 (2009), Article ID 543187, 16 pages
Research Article

Analysis of an Impulsive Predator-Prey System with Monod-Haldane Functional Response and Seasonal Effects

1Department of Mathematics, Kyungpook National University, Daegu 702-701, South Korea
2Department of Mathematics, Pusan National University, Pusan 609-735, South Korea

Received 5 December 2008; Revised 23 March 2009; Accepted 27 April 2009

Academic Editor: J. Rodellar

Copyright © 2009 Hunki Baek 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.


For a class of impulsive predator-prey systems with Monod-Haldane functional response and seasonal effects, we investigate conditions for the local and global stabilities of prey-free solutions and for the permanence of the system by using the Flquet theory of impulsive differential equations and comparison techniques. In addition, we numerically analyze the phenomena caused by seasonal effects and impulsive perturbation. It will be applicable to the controllability for the population of prey and predator.

1. Introduction

One of main goals in population dynamics is to understand the dynamical relationship between predator and prey. Such relationship can be represented as the functional response which refers to the change in the density of prey attached per unit time per predator as the prey density changes. Holling [1] gave three different kinds of functional response for different kinds of species to model the phenomena of predation. These functional responses are monotonic in the first quadrant. But, some experiments and observations indicate that nonmonotonic response occurs at a level [2]. To model such an inhibitory effect, Andrews [2] had suggested a function called the Monod-Haldane function. Sokol and Howell [3] had proposed a simplified Monod-Haldane function having the form Thus, the Monod-Haldane type predator-prey model can be described by the following differential equation [1, 47]: where and represent population densities of prey and predator at time . All parameters are positive constants. Usually, is the intrinsic growth rate of the prey, is the carrying capacity of the prey, the constant is the death rate of the predator, is the rate of conversion of a consumed prey to a predator and measures the level of prey interference with predation.

As well known, there are a number of factors in the environment, which vary with changing seasons and affect various parameters in the ecological models. In fact, countless organisms live in seasonally or diurnally forced environments. For the reason, it is necessary and important to consider models with periodic ecological parameters. There are several ways to apply periodic perturbation in an ecological model [810]. In this paper we consider the intrinsic growth rate in the model (1.3) as periodically varying function of time due to seasonal variation. The seasonality is superimposed as follows: where the parameter represents the degree of seasonality, is the magnitude of the perturbation in , and is the angular frequency of the fluctuation caused by seasonality.

Moreover, there are still some other perturbations such as fire, flood, harvesting seasons, and so forth, that are not suitable to be considered continually. These impulsive perturbations bring sudden change to the model. Thus, it is natural to assume that these perturbations act instantaneously, that is, in the form of impulse [1116].

In this paper, with the ideas discussed above, we consider the following predator-prey system with periodic constant impulsive immigration of the predator and periodic variation in the intrinsic growth rate of the prey: where , is the period of the impulsive immigration or stock of the predator, is the size of immigration or stock of the predator. The parameters and represent the amplitude and frequency of the forcing term, respectively. When and , system (1.5) coincides with (1.3).

Impulsive control methods can be found in almost every field of applied sciences. The theoretical investigation and its application can be found in Bainov and Simeonov [17] and Lakshmikantham et al. [18]. Furthermore, the impulsive differential equations dealing with population dynamics are literate in [9, 13, 14, 16, 1923]. Especially, Zhang et al. [22] studied system (1.5). They investigated the abundance of complex dynamics of system (1.5) and suggested a more executable way for observing chaos and coexistence of attractors by using numerical simulations. Thus, the purpose of this paper is to investigate dynamical behaviors of system (1.5) theoretically. In this context, we find conditions for the local and global stabilities of pest-free periodic solutions and for the permanence of system (1.5) by using the Floquet theory and comparison techniques. Moreover, numerical simulations for the effects of periodic forcing and impulsive perturbations are illustrated.

2. Notations and Lemmas

In this section, we give some notations, definitions and lemmas which will be useful for our main results.

Let , , and . Denote the set of all of nonnegative integers and the right hand side of the model (1.5). Let , then is said to be in a class if

(1) is continuous on , and exists.(2) is locally Lipschitzian in .

Definition 2.1. Let . The upper right derivatives of with respect to the impulsive differential system (1.5) is defined as

Remark 2.2. () The solution of system (1.5) is a piecewise continuous function . That is, is continuous on and exists.
() The smoothness properties of guarantees the global existence and uniqueness of solution of system (1.5); see [18] for the details.

We will use the following important comparison theorem on an impulsive differential equation [18].

Lemma 2.3 (Comparison theorem). Suppose and is continuous on and for , exists, is nondecreasing. Let be the maximal solution of the scalar impulsive differential equation existing on . Then implies that , where is any solution of (2.2).

We now indicate a special case of Lemma 2.3 which provides estimations for the solution of a system of differential inequalities. For this, we let denote the class of real piecewise continuous (real piecewise continuously differentiable) functions defined on .

Lemma 2.4 (see [18]). Let the function satisfy the inequalities where and and are constants, and is a strictly increasing sequence of positive real numbers. Then, for ,

Similar result can be obtained when all conditions of the inequalities in Lemmas 2.3 and 2.4 are reversed. Using Lemma 2.4, it is possible to prove that the solutions of the Cauchy problem (2.3) with strictly positive initial value remain strictly positive.

Lemma 2.5. The positive octant is an invariant region for system (1.5).

Proof. Let be a solution of system (1.5) with a strictly positive initial value . By Lemma 2.4, we can obtain that, for , where and . Thus, and remain strictly positive on .

Now, we give the basic properties of the following impulsive differential equation Then we can easily obtain the following results.

Lemma 2.6. () , , and is a positive periodic solution of (2.7).
() is the solution of (2.7) with , and .
() All nonnegative solutions of (2.7) tend to . That is, as .

It is from Lemma 2.6() that the general solution of (2.7) can be synchronized with the positive periodic solution of (2.7) for sufficiently large and we can obtain the complete expression for the prey-free periodic solution of system (1.5):

The boundness of the solutions of system (1.5) was proven by [22] as follows.

Theorem 2.7 (see [22]). There is an such that for all large enough, where is a solution of the model (1.5).

3. Main Theorems

3.1. Stability of Pest-Free Periodic Solutions

First, we present a condition which guarantees the stability for the prey-free periodic solution .

Theorem 3.1. If , then the pest-free periodic solution is locally asymptotically stable. Moreover, if , then the solution is globally asymptotically stable.

Proof. The local stability of the periodic solution can be determined by considering the behavior of small amplitude perturbations of the solution. Define . Then they may be written as where satisfies and , where is the identity matrix. The linearization of the third and fourth equation of system (1.5) becomes Note that all eigenvalues of are and . Since , the condition is equivalent to the equation According to Floquet theory [17, 18], is locally stable.
Next, to prove the global stability of , let be any solution of system (1.5). Simple calculations yield that . So, under the condition , we know that maintains the local stability and can take a sufficiently small number satisfying It follows from the first equation in (1.5) that for . Then, from Lemma 2.3, we have , where is a solution of the following impulsive differential equation: Since as , for any with large enough. For simplicity we may assume that for all . Since for , we can obtain from Lemmas 2.3 and 2.6 that for sufficiently large. Without loss of generality, we may assume that (3.7) holds for all . From (1.5) and (3.7) we obtain Integrating (3.8) on , we get and thus which implies that as . Further, for , we obtain which implies that as . Now, take a sufficiently small numbers satisfying to prove that as . Since , we may assume that for all . It follows from the second equation in (1.5) that, for , Thus, by Lemma 2.3, we induce that , where is the solution of (2.7) with changed into . Therefore, by taking sufficiently small , we obtain from Lemma 2.6 and (3.7) that tends to as .

Remark 3.2. Since , we can obtain a condition without that is locally stable. In other words, if the condition holds, then the prey-free periodic solution is locally stable. Now, in order to substantiate our theoretical results, let and . It follows from Theorem 3.1 that system (1.5) is locally stable if the condition holds and is globally stable if the condition holds (see Figure 1). Further, Figure 2 indicates that system (1.5) may be globally stable even if .

Figure 1: (a)-(b) Time series of and when and and (c)-(d) time series of and when and .
Figure 2: (a)-(b) Time series of and when and and (c)-(d) time series of and when and .
3.2. Permanence

Before discussing permanence we set up its definition.

Definition 3.3. System (1.5) is permanent if there exist such that any solution of system (1.5) with satisfies

Theorem 3.4. System (1.5) is permanent if .

Proof. Let be a solution of system (1.5) with and . From Theorem 2.7, we may assume that , , , and . Let , . From the second and the forth equations of system (1.5), we obtain that when and . So, it follows from Lemma 2.6 that for all that are large enough. Now we will find an such that for all that are large enough. To do this we will take the following two steps.
Step 1. Since , we can choose , small enough such that and . Suppose that for all . Then it follows from (1.5) that . By Lemma 2.3, we have and as where is a solution of and , . Then there exists such that, for , and Let and . Integrating (3.13) on , , we obtain . Then as which is a contradiction to Theorem 2.7. Hence there exists a such that .
Step 2. If for all , then the proof is complete. If not, let . Then for and, by the continuity of , we have . Suppose that for some . Since and is small enough, we may assume that . Select such that and . Let . Then there are two possible cases for .
Case 1. For all , .
In this case we will show that there exists such that . If not, we have for all . Then for all . It follows from (1.5) that for . By (3.12) with , we have for So we get . Thus, from Lemma 2.3, we obtain for . In the same manner as (3.13), it follows that for . As in Step 1, we have Since , we obtain for all . Integrating this on we obtain Thus which is a contradiction. Now, let . Then for and . So, we have, for , . By the integration of it on , we can get that . For , the same arguments can be continued since .
Case 2. There exists a such that .
Let . Then for and . For , . Integrating this on , we can get that . For , the same arguments can be continued since .

So we know that for since . Moreover, the same arguments as Step 2 can be continued since . Thus we obtain for . Therefore, we complete the proof.

Remark 3.5. We illustrate Figures 5 and 9 as numerical examples of Theorem 3.4. Now assume that to improve slightly the result of Theorem 3.4. Then, by (1.5), we obtain that Note that since . Thus, using these facts and the method used to prove Theorem 3.4, we get the following Corollary.

Corollary 3.6. Suppose that . Then system (1.5) is permanent if .

4. Numerical Analysis on Seasonal Effects and Impulsive Perturbations

In this section, we will investigate the influence of periodic forcing and impulsive perturbations. We fix the parameters and an initial condition as follows: From [5], we know that system (1.3) has an unstable positive equilibrium and a unique stable limit cycle. Zhang et al. [22] have already studied the abundant dynamical behaviors of system (1.5) under the same setting of parameters. Here, we focus on the permanence condition and the control strategy for the population densities of the prey and the predator. Using Theorems 3.1 and 3.4, we display the permanence and stability regions for fixed parameter shown in Figure 3. First, we assume . It follows from Theorems 3.1 and 3.4 that the equation is the border line (see Figure 3(a)), which completely separates the permanence and stability regions. If we consider the prey as a pest and the predator as an enemy of the prey, we should release the predator frequently and abundantly in the aim of exterminating the pest. Or if we regard the prey as a rare variety to be protected, the impulsive perturbations of the predator should occur at rare intervals and their quantity should be very small. Figure 3(a) explains well such situation. For instance, let and ; it follows from Theorems 3.1 and 3.4 that if the prey and the predator coexist, but the prey is extinct after (see Figure 4 ).

Figure 3: Regions induced from Theorems 3.1 and 3.4. (a) , (b) , (c) , (d) .
Figure 4: Bifurcation diagrams of system (1.5) with and . (a) is plotted for ; (b) is plotted for .
Figure 5: (a) Phase portrait of system (1.5) with and and , (b) and (c) Time series of and .

To study the seasonal effects on the prey, we suppose . In this case, the area of unknown region increases according to the value of (see Figure 3). Let and be in the unknown region. Then we know from Figures 6 and 7 that two species coexist if , otherwise the prey is extinct. One of interesting things in Figures 4, 6, and 7 is that the fluctuations of the density of the predator are smaller than those of the prey for small , even though there is the influx of the predator periodically. However, the fluctuations of the prey diminishe according to . These facts suggest that releasing the predator on purpose to increase its density may have no effect. On the other hand, abundant releasing and frequent loosing of the predator can cause destroying the prey and eventually the ecosystem. As shown in Figure 8, we figure out that the longer the period is, the larger the permanence region is and the smaller the stability region is if is not a positive integer.

Figure 6: Bifurcation diagrams of system (1.5) with and . (a) is plotted for , (b) is plotted for .
Figure 7: Bifurcation diagrams of system (1.5) with and . (a) is plotted for , (b) is plotted for .
Figure 8: Regions induced from Theorems 3.1 and 3.4. (a) ; (b) ; (c) ; (d) .
Figure 9: (a) Phase portrait of system (1.5) with , , and ; (b), (c) time series of and .

It means that we should release the predator within a short period to eradicate the prey. On the contrary, impulsive perturbations of the predator should be occurred at long-time intervals for the coexistence of the prey and the predator. If we choose , we can see the coexistence of the prey and the predator as shown in Figure 9.

5. Conclusion

In this paper, we have investigated the effects of periodic forcing in the intrinsic growth rate of the prey and impulsive perturbations of the predator, on the Monod-Haldane functional response predator-prey model. We have shown that there exists an asymptotically stable prey-free periodic solution if the impulsive period is less than a critical value and have also found parameter regions where system (1.5) is permanent. In addition, we have discussed biological roles of parameters , , and in Section 4. Although Theorems 3.1 and 3.4 cannot completely classify the dynamical properties of system (1.5), numerical results supply the border lines which separate two areas having different dynamical properties.


The authors would like to thank the referee for carefully reading of the manuscript and suggesting improvements. This work was supported by the Korea Science and Engineering Foundation (KOSEF) Grant funded by the Korea government (MEST) (no. R01-2008-000-20088-0).


  1. C. S. Holling, “The functional response of predator to prey density and its role in mimicry and population regulations,” Mem. Ent. Sec. Can, vol. 45, pp. 1–60, 1965. View at Google Scholar
  2. J. F. Andrews, “A mathematical model for the continuous culture of macroorganisms untilizing inhibitory substrates,” Biotechnology and Bioengineering, vol. 10, pp. 707–723, 1968. View at Google Scholar
  3. W. Sokol and J. A. Howell, “Kineties of phenol oxidation by ashed cell,” Biotechnology and Bioengineering, vol. 23, pp. 2039–2049, 1980. View at Google Scholar
  4. S.-B. Hsu and T.-W. Huang, “Global stability for a class of predator-prey systems,” SIAM Journal on Applied Mathematics, vol. 55, no. 3, pp. 763–783, 1995. View at Google Scholar · View at MathSciNet
  5. 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. View at Google Scholar · View at MathSciNet
  6. E. Sáez and E. González-Olivares, “Dynamics of a predator-prey model,” SIAM Journal on Applied Mathematics, vol. 59, no. 5, pp. 1867–1878, 1999. View at Google Scholar · View at MathSciNet
  7. J. Sugie, R. Kohno, and R. Miyazaki, “On a predator-prey system of Holling type,” Proceedings of the American Mathematical Society, vol. 125, no. 7, pp. 2041–2050, 1997. View at Google Scholar · View at MathSciNet
  8. J. M. Cushing, “Periodic time-dependent predator-prey systems,” SIAM Journal on Applied Mathematics, vol. 32, no. 1, pp. 82–95, 1977. View at Google Scholar · View at MathSciNet
  9. S. Gakkhar and R. K. Naji, “Chaos in seasonally perturbed ratio-dependent prey-predator system,” Chaos, Solitons & Fractals, vol. 15, no. 1, pp. 107–118, 2003. View at Google Scholar · View at MathSciNet
  10. G. C. W. Sabin and D. Summers, “Chaos in a periodically forced predator-prey ecosystem model,” Mathematical Biosciences, vol. 113, no. 1, pp. 91–113, 1993. View at Publisher · View at Google Scholar
  11. G. J. Ackland and I. D. Gallagher, “Stabilization of large generalized Lotka-Volterra foodwebs by evolutionary feedback,” Physical Review Letters, vol. 93, no. 15, Article ID 158701, 4 pages, 2004. View at Publisher · View at Google Scholar
  12. G. Jiang and Q. Lu, “The dynamics of a prey-predator model with impulsive state feedback control,” Discrete and Continuous Dynamical Systems. Series B, vol. 6, no. 6, pp. 1301–1320, 2006. View at Google Scholar · View at MathSciNet
  13. X. Liu and L. Chen, “Complex dynamics of Holling type II Lotka-Volterra predator-prey system with impulsive perturbations on the predator,” Chaos, Solitons & Fractals, vol. 16, no. 2, pp. 311–320, 2003. View at Google Scholar · View at MathSciNet
  14. B. Liu, Y. Zhang, and L. Chen, “Dynamic complexities in a Lotka-Volterra predator-prey model concerning impulsive control strategy,” International Journal of Bifurcation and Chaos in Applied Sciences and Engineering, vol. 15, no. 2, pp. 517–531, 2005. View at Google Scholar · View at MathSciNet
  15. K. Negi and S. Gakkhar, “Dynamics in a Beddington-DeAngelis prey-predator system with impulsive harvesting,” Ecological Modelling, vol. 206, no. 3-4, pp. 421–430, 2007. View at Publisher · View at Google Scholar
  16. X. Song and Y. Li, “Dynamic behaviors of the periodic predator-prey model with modified Leslie-Gower Holling-type II schemes and impulsive effect,” Nonlinear Analysis: Real World Applications, vol. 9, no. 1, pp. 64–79, 2008. View at Google Scholar · View at MathSciNet
  17. D. D. Bainov and P. S. Simeonov, Impulsive Differential Equations: Periodic Solutions and Applications, vol. 66 of Pitman Monographs and Surveys in Pure and Applied Mathematics, Longman Science & Technical, Harlo, UK, 1993.
  18. V. Lakshmikantham, D. D. Baĭnov, and P. S. Simeonov, Theory of Impulsive Differential Equations, vol. 6 of Series in Modern Applied Mathematics, World Scientific, Teaneck, NJ, USA, 1989. View at MathSciNet
  19. H. Baek, “Dynamic complexites of a three-species Beddington-DeAngelis system with impulsive control strategy,” Acta Applicandae Mathematicae, pp. 1–16, 2008. View at Publisher · View at Google Scholar
  20. W. Wang, H. Wang, and Z. Li, “The dynamic complexity of a three-species Beddington-type food chain with impulsive control strategy,” Chaos, Solitons & Fractals, vol. 32, no. 5, pp. 1772–1785, 2007. View at Google Scholar · View at MathSciNet
  21. Z. Xiang and X. Song, “The dynamical behaviors of a food chain model with impulsive effect and Ivlev functional response,” Chaos, Solitons & Fractals, vol. 39, no. 5, pp. 2282–2293, 2009. View at Publisher · View at Google Scholar
  22. S. Zhang, D. Tan, and L. Chen, “Chaos in periodically forced Holling type IV predator-prey system with impulsive perturbations,” Chaos, Solitons & Fractals, vol. 27, no. 4, pp. 980–990, 2006. View at Google Scholar · View at MathSciNet
  23. S. Zhang and L. Chen, “A study of predator-prey models with the Beddington-DeAnglis functional response and impulsive effect,” Chaos, Solitons & Fractals, vol. 27, no. 1, pp. 237–248, 2006. View at Google Scholar · View at MathSciNet