Research Article | Open Access

# Stability and Bifurcation of Two Kinds of Three-Dimensional Fractional Lotka-Volterra Systems

**Academic Editor:**Yuncai Wang

#### Abstract

Two kinds of three-dimensional fractional Lotka-Volterra systems are discussed. For one system, the asymptotic stability of the equilibria is analyzed by providing some sufficient conditions. And bifurcation property is investigated by choosing the fractional order as the bifurcation parameter for the other system. In particular, the critical value of the fractional order is identified at which the Hopf bifurcation may occur. Furthermore, the numerical results are presented to verify the theoretical analysis.

#### 1. Introduction

In recent years, fractional calculus has attracted much attention of researchers. It has been pointed out that fractional calculus plays an outstanding role in modelling and simulation of systems, such as viscoelastic systems, dielectric polarization, electromagnetic waves, heat conduction, robotics, and biological systems. In fact, fractional derivatives provide an excellent instrument for the description of memory and hereditary properties of various materials and processes in comparison with the classical integer-order counterparts. Therefore, it may be more important and useful to investigate the fractional systems.

Traditionally, the fractional differential equation defined by mathematicians is a Riemann-Liouville fractional derivative [1]. But this definition is less popular because of the fact that it requires initial conditions to be expressed in terms of fractional integrals and their derivatives. Meanwhile, there is no known physical interpretation for such types of initial conditions. In contrast, the alternative definition of the fractional derivative given by Caputo [2] has the advantage of only requiring the initial conditions given in terms of integer-order derivatives. These initial conditions of integer-order derivatives can be measured accurately and represent well-understood features of a physical situation. In [2], it has been pointed out that Caputo’s derivative is equivalent to the Riemann-Liouville derivative under homogeneous initial conditions and some smoothness conditions. Therefore, Caputo’s definition of fractional derivative is used throughout in this paper.

As is well known, in the field of mathematical biology, the traditional Lotka-Volterra systems are very important mathematical models which describe multispecies population dynamics in a nonautonomous environment. Many important and interesting results on the dynamical behaviors for the Lotka-Volterra systems have been found in [3–9], such as the existence and uniqueness of solutions, the permanence, extinction, global asymptotic behavior, and bifurcation. Because of the good memory and hereditary properties of fractional derivatives, it is often necessary to study the corresponding fractional systems. Therefore, the dynamical analysis of the fractional Lotka-Volterra systems has attracted a great deal of attention due to its theoretical and practical significance.

Many important results regarding stability of fractional systems have been obtained. For instance, the stability, existence, uniqueness, and numerical solution of the fractional logistic equation are investigated in [10]. The stability and solutions of fractional predator-prey and rabies models are discussed in [11]. In addition, bifurcation properties of fractional systems have been studied in some papers. For example, conditions for the occurrence of Hopf’s bifurcation are explored based on numerical simulations in [12]. The critical values of the fractional order are identified for which Hopf’s bifurcation may occur based on the stability analysis in [13]. Thus, it is significant to study the dynamical behaviors in the fractional population systems.

To the best of our knowledge, some papers have concentrated on the dynamic investigation of the fractional population systems [10, 11]. However, there are few results on bifurcation phenomena of the fractional population systems. Therefore, in the paper, we mainly consider stability and bifurcation in the three-dimensional fractional Lotka-Volterra systems.

Motivated by the above discussions, some dynamical properties of two kinds of three-dimensional fractional Lotka-Volterra systems are investigated in this paper. Existence and uniqueness of solutions are considered. Some sufficient conditions are provided for the asymptotic stability of equilibria. Specifically, bifurcation behaviors are analyzed by formulating the critical values of the fractional order at which Hopf’s bifurcations may take place.

The rest of this paper is organized as follows. In Section 2, a three-dimensional fractional Lotka-Volterra predator-prey system with interspecific competition is introduced. And the asymptotic stability of the system is studied. In Section 3, a three-dimensional fractional Lotka-Volterra predator-prey system is provided, and bifurcation properties are investigated. The numerical results in Section 4 are given to verify the theoretical findings. Finally, the paper is concluded in Section 5.

#### 2. Stability Analysis of a Three-Dimensional Fractional Lotka-Volterra Predator-Prey System with Interspecific Competition

Consider a three-dimensional fractional Lotka-Volterra system with the initial values , , where ; especially when , the system (1) is a classical integer-order system. All constant coefficients and () can be arbitrary positive real numbers. represents the density of prey species at time , and , represent the densities of predator species at time . In this case, system (1) can be regarded as a fractional Lotka-Volterra predator-prey system with interspecific competition.

In the following, existence and uniqueness of solutions for system (1) are given. In addition, the important results related to the stability of the fractional systems are presented to provide the theoretical bases for the further study.

Here, the fractional Lotka-Volterra system (1) can be rewritten in the form where , , and

*Definition 1. *For , let be the set of continuous column vectors on the interval . The norm of is given by .

Theorem 2. *System (2) has a unique solution if .*

*Proof. *Let , then implies . In addition, take and ; the following inequality holds:
where , and and are positive and satisfy , as a result of . Based on Theorems 2.1 and 2.2 in [14], system (2) has a unique solution.

Theorem 3 (see [15]). *The linear autonomous system is asymptotically stable if and only if
**
where , , and ; denotes the set of all eigenvalues of the matrix .*

Theorem 4. *Let be an equilibrium of the nonlinear system (1), then the equilibrium is locally asymptotically stable if
**
where ; denotes the set of all eigenvalues of the Jacobian matrix .*

*Proof. *The proof follows from Theorem 3 and [11].

In the following, the stability of system (1) is investigated by giving some appropriate conditions. The asymptotic stability of the equilibria is demonstrated based on Theorem 4. Through simple calculation, the equilibria of system (1) are obtained and denoted as where , and .

Because of the fact that all constant coefficients of system (1) are positive, , , and are in contradiction with the actual situation; hence the asymptotical stability of other five equilibria will be studied in detail.

Theorem 5. *For the three-dimensional fractional Lotka-Volterra system (1), the following results can be obtained.*(a)* is unstable;*(b)* is locally asymptotically stable if ;*(c)* is locally asymptotically stable if ;*(d)* is locally asymptotically stable if ;*(e)* is locally asymptotically stable if .*

*Proof. *For , its Jacobian matrix is
and the eigenvalues of satisfy , and ; hence the equilibrium is unstable.

For , its Jacobian matrix is
and the eigenvalues of satisfy , , and ; hence the equilibrium is locally asymptotically stable.

For , use the notations below:
and its characteristic equation is
Based on the condition from (c), the following formulas can be easily got
where , and . Then, the following results can be obtained:
Hence the equilibrium is locally asymptotically stable.

Similarly, it can be readily derived that the equilibrium is locally asymptotically stable.

For , let ; the Jacobian matrix of can be written as
and its characteristic equation is
where , and . For simplicity, the equivalent characteristic equation is introduced as follows:
On the basis of the above equivalent substitutions, the following inequalities can be gained:
where , and . Using the proof by contradiction, it can be concluded that the eigenvalues of satisfy
Hence the equilibrium is locally asymptotically stable.

For the further dynamic investigation of the fractional population systems, the other fractional Lotka-Volterra systems will be considered in the following section. Particularly, bifurcation properties for the system will be studied in detail.

#### 3. Bifurcation Analysis of a Three-Dimensional Fractional Lotka-Volterra Predator-Prey System

Consider a three-dimensional fractional Lotka-Volterra system: with the initial values , , where , , and the other constant coefficients are positive. represents the density of prey species at time , and , represent the densities of predator species at time . In this case, system (19) can be regarded as a fractional Lotka-Volterra predator-prey system.

On the basis of Theorem 2, it is not difficult to prove that system (19) has a unique solution in a similar way.

It is clear that there are eight equilibria for system (19). Here, we focus on the bifurcation investigation of the equilibrium which can be called a positive equilibrium when some conditions are satisfied. The equilibrium is obtained as where , and . And its Jacobian matrix can be expressed as Furthermore, the eigenvalues of satisfy the characteristic equation where , , and .

In the following, by choosing the fractional order as the bifurcation parameter and analyzing the associated characteristic equation (22) of system (19) at the positive equilibrium, we investigate the bifurcation phenomena of the positive equilibrium of system (19) and obtain the conditions under which system (19) undergoes a Hopf bifurcation.

Proposition 6. *The positive equilibrium of system (19) is locally asymptotically stable if and only if all the following conditions are satisfied:*(i)*,*(ii)*, and*(iii)*.*

*Proof. *For the characteristic equation (22), the root , and satisfy the equation . It is clear that if and only if the conditions (ii) and (iii) hold. Based on Theorem 4, Proposition 6 proves to be true.

In addition, by analyzing the condition (iii) of Proposition 6 in detail, the following results can be gained.

Proposition 7. *With respect to system (19), if and , the following statements can be obtained.*(a)*If **, the equilibrium ** is locally asymptotically stable, for any **.*(b)*If **, the equilibrium ** is locally asymptotically stable if and only if **, where **.*(c)*If **, the equilibrium ** is unstable for any **.*

* Proof. *The conclusions (a) and (c) are obvious. For the statement (b), due to , the equation has two complex roots , and their real part is . Then . Besides, according to the condition if and only if . Based on Theorem 4, it is concluded that Proposition 7 is true.

*Remark 8. *It is apparent that the critical value satisfies . When , is locally asymptotically stable; when , and specially , is unstable. That is to say, it has verified that fractional differential equations are, at least, as stable as their integer-order counterparts [4].

*Remark 9. *Under the situation of statement (b), a bifurcation phenomenon must happen at the critical value . However, it is difficult to confirm precise bifurcation type. As an interesting bifurcation behavior, Hopf’s bifurcation is expected to take place.

According to Proposition 7, if some appropriate conditions about the constant coefficients of system (19) can be found so that statement (b) is satisfied, system (19) will undergo a bifurcation phenomenon. And the critical value of the bifurcation parameter can be expressed by the constant coefficients of system (19). From this, the following theorem is specifically proposed.

Theorem 10. *With respect to system (19), if the following conditions are satisfied:*(i)*, , ;*(ii)*,**then the positive equilibrium is locally asymptotically stable if and only if , where
*

*Proof. *According to the condition (i), the equilibrium can be expressed as
For (22), the following results can be obtained
where and . Obviously, the above conclusions satisfy statement (b) of Proposition 7, then it can be derived that
Hence, the positive equilibrium of system (19) is locally asymptotically stable if and only if .

According to the statement of Theorem 10, it can be concluded that the positive equilibrium is locally asymptotically stable if and only if . At , the Hopf bifurcation is expected to take place. As increases above the critical value , the positive equilibrium is unstable and a limit cycle is expected to appear in the proximity of due to the Hopf bifurcation phenomenon.

The analysis of periodic solutions in fractional dynamical systems is a very recent and promising research topic. As a consequence, the nonexistence of exact periodic solutions in time invariant fractional systems is obtained [16]. As an application, it is emphasized that the limit cycle observed in numerical simulations of a simple fractional neural network cannot be an exact periodic solution of the system [17]. In addition, there are some other papers providing the numerical evidences of limit cycles.

*Remark 11. *Even though exact periodic solutions do not exist in autonomous fractional systems [16, 17], limit cycles have been observed by numerical simulations in many systems such as a fractional neural system [13], a fractional Van der Pol system [18], fractional Chua and Chen’s systems [19, 20], and a fractional financial system [21].

#### 4. Numerical Simulation

In this paper, an Adams-type predictor-corrector method is used for the numerical solutions of fractional differential equations. This method has been introduced in [22, 23] and further investigated in [24–27]. In order to verify the theoretical analysis, the following numerical results are given.

For system (1), the approximate solutions are displayed in Figure 1 for the step size 0.005 and different values of , , respectively. Taking , and , and choosing the initial values , and , the equilibrium is . Then, Figure 1 shows that the equilibrium is locally asymptotically stable. Namely, the fifth conclusion of Theorem 5 is verified. Similarly, the other conclusions of Theorem 5 can be confirmed.

For system (19), the approximate solutions are displayed in Figures 2, 3, and 4 for the step size 0.001 and different values of , , and . Taking , and , and choosing the initial values , the positive equilibrium is , and the critical value is . Indeed, Figures 2–4 present the fact that the positive equilibrium is locally asymptotically stable when , and when increases across , an asymptotically stable limit cycle appears in a neighborhood of the positive equilibrium .

**(a)**

**(b)**

**(c)**

#### 5. Conclusion

In this paper, two kinds of three-dimensional fractional Lotka-Volterra systems have been studied. The main results are divided into two parts. On the one hand, for system (1), the asymptotic stability of the equilibria is investigated by providing simple and reasonable sufficient conditions. And simulation results prove to be quite consistent with the theoretical findings. On the other hand, for system (19), the conditions which could lead to bifurcation phenomena are obtained. Specifically, the fractional order is chosen as the bifurcation parameter and the expression of the critical value is precisely derived. Furthermore, the numerical result is presented to illustrate that Hopf’s bifurcation can take place.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgment

This work is supported by the National Nature Science Foundation of China under Grant no. 11371049 and the Science Foundation of Beijing Jiaotong University under Grant 2011JBM130.

#### References

- K. S. Miller and B. Ross,
*An Introduction to the Fractional Calculus and Fractional Differential Equations*, Wiley-Interscience Publication, New York, NY, USA, 1993. - I. Podlubny,
*Fractional Differential Equations*, Academic Press, New York, NY, USA, 1999. - S. Ahmad and A. C. Lazer, “Average conditions for global asymptotic stability in a nonautonomous Lotka-Volterra system,”
*Nonlinear Analysis, Theory, Methods and Applications*, vol. 40, no. 1, pp. 37–49, 2000. View at: Google Scholar | Zentralblatt MATH - S. Ahmad and A. C. Lazer, “Average growth and total permanence in a competitive Lotka-Volterra System,”
*Annali di Matematica Pura ed Applicata*, vol. 185, supplement 5, pp. S47–S67, 2006. View at: Publisher Site | Google Scholar | Zentralblatt MATH - P. van den Driessche and M. L. Zeeman, “Three-dimensional competitive Lotka-Volterra systems with no periodic orbits,”
*SIAM Journal on Applied Mathematics*, vol. 58, no. 1, pp. 227–234, 1998. View at: Publisher Site | Google Scholar | Zentralblatt MATH - Z. Teng and L. Chen, “Global asymptotic stability of periodic Lotka-Volterra systems with delays,”
*Nonlinear Analysis: Theory, Methods and Applications*, vol. 45, no. 8, pp. 1081–1095, 2001. View at: Publisher Site | Google Scholar | Zentralblatt MATH - N. Fang and X. X. Chen, “Permanence of a discrete multispecies Lotka-Volterra competition predator-prey system with delays,”
*Nonlinear Analysis: Real World Applications*, vol. 9, no. 5, pp. 2185–2195, 2008. View at: Publisher Site | Google Scholar | Zentralblatt MATH - G. Lu and Z. Lu, “Permanence for two-species Lotka-Volterra cooperative systems with delays,”
*Mathematical Biosciences and Engineering*, vol. 5, no. 3, pp. 477–484, 2008. View at: Publisher Site | Google Scholar | Zentralblatt MATH - X. P. Yan and W. T. Li, “Stability and Hopf bifurcation for a delayed cooperative system with diffusion effects,”
*International Journal of Bifurcation and Chaos*, vol. 18, no. 2, pp. 441–453, 2008. View at: Publisher Site | Google Scholar | Zentralblatt MATH - A. M. A. El-Sayed, A. E. M. El-Mesiry, and H. A. A. El-Saka, “On the fractional-order logistic equation,”
*Applied Mathematics Letters*, vol. 20, no. 7, pp. 817–823, 2007. View at: Publisher Site | Google Scholar | Zentralblatt MATH - E. Ahmed, A. M. A. El-Sayed, and H. A. A. El-Saka, “Equilibrium points, stability and numerical solutions of fractional-order predator-prey and rabies models,”
*Journal of Mathematical Analysis and Applications*, vol. 325, no. 1, pp. 542–553, 2007. View at: Publisher Site | Google Scholar | Zentralblatt MATH - H. A. El-Saka, E. Ahmed, M. I. Shehata, and A. M. A. El-Sayed, “On stability, persistence, and Hopf bifurcation in fractional order dynamical systems,”
*Nonlinear Dynamics*, vol. 56, no. 1-2, pp. 121–126, 2009. View at: Publisher Site | Google Scholar | Zentralblatt MATH - E. Kaslik and S. Sivasundaram, “Nonlinear dynamics and chaos in fractional-order neural networks,”
*Neural Networks*, vol. 32, pp. 245–256, 2012. View at: Publisher Site | Google Scholar | Zentralblatt MATH - K. Diethelm and N. J. Ford, “Analysis of fractional differential equations,”
*Journal of Mathematical Analysis and Applications*, vol. 265, no. 2, pp. 229–248, 2002. View at: Publisher Site | Google Scholar | Zentralblatt MATH - D. Matignon, “Stability results for fractional differential equations with applications to control processing,” in
*Proceedings of the Computational Engineering in Systems Applications*, vol. 2, pp. 963–968, Lille, France, July 2012. View at: Google Scholar - M. S. Tavazoei and M. Haeri, “A proof for non existence of periodic solutions in time invariant fractional order systems,”
*Automatica*, vol. 45, no. 8, pp. 1886–1890, 2009. View at: Publisher Site | Google Scholar | Zentralblatt MATH - E. Kaslik and S. Sivasundaram, “Non-existence of periodic solutions in fractional-order dynamical systems and a remarkable difference between integer and fractional-order derivatives of periodic functions,”
*Nonlinear Analysis: Real World Applications*, vol. 13, no. 3, pp. 1489–1497, 2012. View at: Publisher Site | Google Scholar | Zentralblatt MATH - R. S. Barbosa, J. A. T. MacHado, B. M. Vinagre, and A. J. Calderón, “Analysis of the van der Pol oscillator containing derivatives of fractional order,”
*Journal of Vibration and Control*, vol. 13, no. 9-10, pp. 1291–1301, 2007. View at: Publisher Site | Google Scholar | Zentralblatt MATH - D. Cafagna and G. Grassi, “Bifurcation and chaos in the fractional-order Chen system via a time-domain approach,”
*International Journal of Bifurcation and Chaos*, vol. 18, no. 7, pp. 1845–1863, 2008. View at: Publisher Site | Google Scholar | Zentralblatt MATH - D. Cafagna and G. Grassi, “Fractional-order Chua's circuit: time-domain analysis, bifurcation, chaotic behavior and test for chaos,”
*International Journal of Bifurcation and Chaos*, vol. 18, no. 3, pp. 615–639, 2008. View at: Publisher Site | Google Scholar | Zentralblatt MATH - M. S. Abd-Elouahab, N. E. Hamri, and J. Wang, “Chaos control of a fractional-order financial system,”
*Mathematical Problems in Engineering*, vol. 2010, Article ID 270646, 18 pages, 2010. View at: Publisher Site | Google Scholar | Zentralblatt MATH - K. Diethelm and A. Freed, “On the solution of nonlinear fractional order differential equations used in the modelling of viscoplasticity,” in
*Scientific Computing in Chemical Engineering II-Computational Fluid Dynamics, Reaction Engineering, and Molecular Properties*, F. Keil, W. Mackens, H. Voß, and J. Werther, Eds., pp. 217–224, Springer, Heidelberg, Germany, 1999. View at: Google Scholar - K. Diethelm and A. Freed, “The FracPECE subroutine for the numerical solution of differential equations of fractional order,” in
*Forschung und Wissenschaftliches Rechnen 1998*, S. Heinzel and T. Plesser, Eds., pp. 57–71, Gesellschaft fr Wisseschaftliche Datenverarbeitung, Göttingen, Germany, 1999. View at: Google Scholar - R. Zhao,
*Dynamical Analysis of Fractional-Order Species Models*, College of Science Beijing Jiaotong University, Beijing, China, 2011. - E. Ahmed, A. M. A. El-Sayed, A. E. M. El-Mesiry, and H. A. A. El-Saka, “Numerical solution for the fractional replicator equation,”
*International Journal of Modern Physics C*, vol. 16, no. 7, pp. 1017–1025, 2005. View at: Publisher Site | Google Scholar | Zentralblatt MATH - E. Ahmed, A. M. A. El-Sayed, and H. A. A. El-Saka, “On some Routh-Hurwitz conditions for fractional order differential equations and their applications in Lorenz, Rössler, Chua and Chen systems,”
*Physics Letters A*, vol. 358, no. 1, pp. 1–4, 2006. View at: Publisher Site | Google Scholar | Zentralblatt MATH - A. E. M. El-Mesiry, A. M. A. El-Sayed, and H. A. A. El-Saka, “Numerical methods for multi-term fractional (arbitrary) orders differential equations,”
*Applied Mathematics and Computation*, vol. 160, no. 3, pp. 683–699, 2005. View at: Publisher Site | Google Scholar | Zentralblatt MATH

#### Copyright

Copyright © 2014 Jinglei Tian 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.