- About this Journal ·
- Abstracting and Indexing ·
- Advance Access ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
Abstract and Applied Analysis
Volume 2014 (2014), Article ID 835310, 16 pages
Hopf Bifurcation and Global Periodic Solutions in a Predator-Prey System with Michaelis-Menten Type Functional Response and Two Delays
1Department of Applied Mathematics, Kunming University of Science and Technology, Kunming, Yunnan 650093, China
2Department of Mathematics and Information Science, Zhoukou Normal University, Zhoukou, Henan 466001, China
Received 27 November 2013; Accepted 31 March 2014; Published 11 May 2014
Academic Editor: Yuming Chen
Copyright © 2014 Yunxian Dai 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.
We consider a predator-prey system with Michaelis-Menten type functional response and two delays. We focus on the case with two unequal and non-zero delays present in the model, study the local stability of the equilibria and the existence of Hopf bifurcation, and then obtain explicit formulas to determine the properties of Hopf bifurcation by using the normal form method and center manifold theorem. Special attention is paid to the global continuation of local Hopf bifurcation when the delays .
In , Xu and Chaplain studied the following delayed predator-prey model with Michaelis-Menten type functional response: with initial conditions where , and denote the densities of the prey, predator, and top predator population, respectively. are positive constants. , and are nonnegative constants. denote the delay in the negative feedback of the prey, predator, and top predator crowding, respectively. , are constant delays due to gestation; that is, mature adult predators can only contribute to the production of predator biomass. . are continuous bounded functions in the interval . The authors proved that the system is uniformly persistent under some appropriate conditions. By means of constructing suitable Lyapunov functional, sufficient conditions are derived for the global asymptotic stability of the positive equilibrium of the system.
Time delays of one type or another have been incorporated into systems by many researchers since a time delay could cause a stable equilibrium to become unstable and fluctuation. In [2–12], authors showed effects of two delays on dynamical behaviors of system.
It is well known that periodic solutions can arise through the Hopf bifurcation in delay differential equations. However, these periodic solutions bifurcating from Hopf bifurcations are generally local. Under some circumstances, periodic solutions exist when the parameter is far away from the critical value. Therefore, global existence of Hopf bifurcation is a more interesting and difficult topic. A great deal of research has been devoted to the topics [13–21]. In this paper, let in (1); we consider Hopf bifurcation and global periodic solutions of the following system with two unequal and nonzero delays: with initial conditions Our goal is to investigate the possible stability switches of the positive equilibrium and stability of periodic orbits arising due to a Hopf bifurcation when one of the delays is treated as a bifurcation parameter. Special attention is paid to the global continuation of local Hopf bifurcation when the delays .
This paper is organized as follows. In Section 2, by analyzing the characteristic equation of the linearized system of system (3) at positive equilibrium, the sufficient conditions ensuring the local stability of the positive equilibrium and the existence of Hopf bifurcation are obtained . Some explicit formulas determining the direction and stability of periodic solutions bifurcating from Hopf bifurcations are demonstrated by applying the normal form method and center manifold theory due to Hassard et al.  in Section 3. In Section 4, we consider the global existence of these bifurcating periodic solutions  with two different delays. Some numerical simulation results are included in Section 5.
2. Stability of the Positive Equilibrium and Local Hopf Bifurcations
In this section, we first study the existence and local stability of the positive equilibrium and then investigate the effect of delay and the conditions for existence of Hopf bifurcations.
There are at most four nonnegative equilibria for system (3): where and satisfy where is a nonnegative equilibrium point if there is a positive solution of (6), and is a nonnegative equilibrium point if there is a positive solution of (7).
Let be the arbitrary equilibrium point, and let ; still denote by , respectively; then the linearized system of the corresponding equations at is as follows: where all the others of , and are .
The characteristic equation for system (8) is where
We consider the following cases.
(1) . The characteristic equation reduces to There are always a positive root and two negative roots of (12); hence is a saddle point.
(2) . Equation (10) takes the form There is a positive root if ; hence, is a saddle point. If , is locally asymptotically stable.
Lemma 1. Consider the exponential polynomial where and are constants. As vary, the sum of the order of the zeros of on the open right half plane can change only if a zero appears on or crosses the imaginary axis.
By using Lemma 1, we can easily obtain the following results.
Lemma 2. If is a nonnegative equilibrium point, then(1) is unstable if ;(2) is locally asymptotically stable if , and .
Proof. (1) is a root of (14); if , then is unstable.
(2) Clearly, is not a root of (14); we should discuss the following equation instead of (14): Assume that with is a solution of (16). Substituting into (16) and separating the real and imaginary parts yield which implies If , that is , there is no real root of (16). Hence there is no purely imaginary root of (18). When , (16) reduces to If and , both roots of (19) have negative real parts. Thus, by using Lemma 1, when , and , is locally asymptotically stable.
The associated characteristic equation of system (3) is Let(H5). By Routh-Hurwitz criterion, we have the following.
Theorem 3. For , assume that hold. Then when , the positive equilibrium of system (3) is locally asymptotically stable.
Case b. Consider
The associated characteristic equation of system (3) is
We want to determine if the real part of some root increases to reach zero and eventually becomes positive as varies. Let be a root of (21); then we have Separating the real and imaginary parts, we have It follows that where .
Denoting , (24) becomes
Let we have
If , then . We can know that (25) has at least one positive root.
If , we obtain that when , (25) has no positive roots for . On the other hand, when , the following equation has two real roots: . Because of and are the local minimum and the local maximum of , respectively. By the above analysis, we immediately obtain the following.
Without loss of generality, we assume that (25) has three positive roots, defined by , respectively. Then (24) has three positive roots: From (23) we have Thus, if we denote where ; then is a pair of purely imaginary roots of (21) corresponding to . Define
According to (35), we have where . Notice that , then we have the following lemma.
Lemma 5. Suppose that and , where is defined by (26); then has the same sign with .
Theorem 6. For , suppose that hold. (i)If and , then all roots of (10) have negative real parts for all , and the positive equilibrium is locally asymptotically stable for all .(ii)If either or , and , then has at least one positive roots, and all roots of (23) have negative real parts for , and the positive equilibrium is locally asymptotically stable for .(iii)If (ii) holds and , then system (3)undergoes Hopf bifurcations at the positive equilibrium for .
Case c. Consider
The associated characteristic equation of system (3) is
Similar to the analysis of Case , we get the following theorem.
Theorem 7. For , suppose that hold. (i)If and , then all roots of (38) have negative real parts for all , and the positive equilibrium is locally asymptotically stable for all .(ii)If either or , , and , then has at least one positive root , and all roots of (38) have negative real parts for , and the positive equilibrium is locally asymptotically stable for .(iii)If (ii) holds and , then system (3) undergoes Hopf bifurcations at the positive equilibrium for ,
where where ; ; then is a pair of purely imaginary roots of (38) corresponding to . Define
Case d. Consider
Let be a root of (41); then we have Separating the real and imaginary parts, we have It follows that where
Let When , (41) has a pair of purely imaginary roots for .
In the following, we assume that(H6).
Theorem 8. For , , , suppose that is satisfied. If and , then the positive equilibrium is locally asymptotically stable for . System (3) undergoes Hopf bifurcations at the positive equilibrium for .
3. Direction and Stability of the Hopf Bifurcation
In Section 2, we obtain the conditions under which system (3) undergoes the Hopf bifurcation at the positive equilibrium . In this section, we consider with and regard as a parameter. We will derive the explicit formulas determining the direction, stability, and period of these periodic solutions bifurcating from equilibrium at the critical values by using the normal form and the center manifold theory developed by Hassard et al. . Without loss of generality, denote any one of these critical values by , at which (43) has a pair of purely imaginary roots and system (3) undergoes Hopf bifurcation from .
Throughout this section, we always assume that . Let and . Then is the Hopf bifurcation value of system (3). System (3) may be written as a functional differential equation in where , and where , and Obviously, is a continuous linear function mapping into . By the Riesz representation theorem, there exists a matrix function , whose elements are of bounded variation such that In fact, we can choose where is Dirac-delta function. For , define Then when , the system is equivalent to where , . For , define and a bilinear inner product where . Let ; then and are adjoint operators. By the discussion in Section 2, we know that are eigenvalues of . Thus, they are also eigenvalues of . We first need to compute the eigenvector of and corresponding to and , respectively. Suppose that is the eigenvector of corresponding to . Then . From the definition of , and we can easily obtain , where and . Similarly, let be the eigenvector of corresponding to . By the definition of , we can compute From (57), we have Thus, we can choose such that .
In the remainder of this section, we follow the ideas in Hassard et al.  and use the same notations as there to compute the coordinates describing the center manifold at . Let be the solution of (48) when . Define On the center manifold , we have where and are local coordinates for center manifold in the direction of and . Note that is real if is real. We consider only real solutions. For the solution of (48), since , we have where By (62), we have , and then It follows together with (50) that