Complexity
Volume 2017 (2017), Article ID 2578043, 12 pages
https://doi.org/10.1155/2017/2578043

Research Article

Delay-Induced Oscillations in a Competitor-Competitor-Mutualist Lotka-Volterra Model

Changjin Xu

Guizhou Key Laboratory of Economics System Simulation, Guizhou University of Finance and Economics, Guiyang 550004, China

Correspondence should be addressed to Changjin Xu; xcj403@126.com

Received 14 January 2017; Accepted 30 March 2017; Published 26 April 2017

Academic Editor: Alicia Cordero

Copyright © 2017 Changjin Xu. 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.

This paper deals with a competitor-competitor-mutualist Lotka-Volterra model. A series of sufficient criteria guaranteeing the stability and the occurrence of Hopf bifurcation for the model are obtained. Several concrete formulae determine the properties of bifurcating periodic solutions by applying the normal form theory and the center manifold principle. Computer simulations are given to support the theoretical predictions. At last, biological meaning and a conclusion are presented.

1. Introduction

The study of dynamical behaviors for tremendous predator-prey models has been a hot issue in population dynamics in the past few decades. Many results have been reported [1–11]. In the real world, any biological or environmental parameters are naturally subject to fluctuation in time. The effects of a periodically varying environment are important for evolutionary theory as the selective forces on systems in a fluctuating environment differ from those in a stable environment. Meanwhile, time delay due to gestation is common example because, generally, the consumption of prey by the predator throughout its past history governs the present birth rate of the predator. Based on all the above point, Lv et al. [12] had investigated the periodic solution of the following competitor-competitor-mutualist Lotka-Volterra model by using Krasnoselskii’s fixed point theorem

where and denote the densities of competing species at time and denotes the density of cooperating species at time . and are -periodic functions . The parameters are the feedback time delay of different species. In detail, one can see [12].

It is well known that the research on the Hopf bifurcation, especially on the stability of bifurcating periodic solutions and direction of Hopf bifurcation, is one of the most important themes on the predator-prey dynamics. There are a great deal of papers which deal with this topic [11, 13–22]. The purpose of this paper is to discuss the stability and the properties of Hopf bifurcation of model (1). To simplify the analysis for model (1), we make the following assumptions: all biological and environmental parameters are constants in time and only the feedback time delay of competing species to the growth of the species itself and the feedback time delay of cooperating species to the growth of the species itself exist and are the same. Then system (1) can be described as the form

In this paper, we consider the effect of time delay on the dynamics of system (2). We not only give the conditions on the stability of the positive equilibrium of (2) and the existence of periodic solutions but also derive the formulae for determining the properties of a Hopf bifurcation.

The remainder of the paper is organized as follows. In Section 2, we investigate the stability of the positive equilibrium and the occurrence of local Hopf bifurcations. In Section 3, the direction and stability of the local Hopf bifurcation are established. In Section 4, numerical simulations are carried out to illustrate the validity of the main results. Biological explanations and some main conclusions are drawn in Section 5.

2. Stability of the Positive Equilibrium and Local Hopf Bifurcations

Consider the realistic implication and actual application of biological system; in this section, we shall only study the stability of the positive equilibrium and the existence of local Hopf bifurcations. It is easy to see that system (2) has a unique positive equilibrium if the condition

holds, where

Let ,  ,  and and still denote by , and then (2) takes the form

where

The linearization of (6) near is given by

whose characteristic equation takes the form

That is,

where

Multiplying on both sides of (8), it is easy to obtain

We need the following lemma to discuss the stability of the positive equilibrium.

Lemma 1 (see [23]). For the transcendental equation

as vary, the sum of orders of the zeros of in the open right half complex plane can change, only if a zero appears on or crosses the imaginary axis.

For , (10) becomes

Obviously, . By the Routh-Hurwitz criteria, it follows that all eigenvalues of (12) have negative real parts if and only if the condition

is fulfilled.

For ,   is a root of (10) if and only if

Separating the real and imaginary parts, we get

It follows from (14) that

According to , then (15) takes the form

It is easy to see that (16) is equivalent to

where

where

Let and denote

It is easy to obtain that Set

Let Then (21) becomes

where and  

Define ,   By (22), then we obtain

By the discussion above, we can obtain the expression of , say

where is a function with respect to . Substitute (24) into (15); then we can easily get the expression of , say

where is a function with respect to . Thus we obtain

If all the coefficients of system (2) are given, it is easy to use computer to calculate the roots of (26) (say ). Then from (24), we derive

Let be a root of (10) near , , and Due to functional differential equation theory, for every , there exists such that is continuously differentiable in for . Substituting into the left hand side of (10) and taking derivative with respect to , we have

Then

where

In order to obtain the main results in this paper, it is necessary to make the following assumption:

In view of Lemma 1, it is easy to obtain the following result on stability and bifurcation of system (2).

Theorem 2. Suppose that – hold; then

(i)for system (2), its positive equilibrium is asymptotically stable for ;(ii)system (2) undergoes a Hopf bifurcation at the positive equilibrium when ; that is, system (2) has a branch of periodic solutions bifurcating from the positive equilibrium near .

3. Direction and Stability of the Hopf Bifurcation

In this section, we shall analyze the direction, stability, and period of these periodic solutions bifurcating from the positive equilibrium at these critical value of , by using techniques from normal form and center manifold theory [24]. Throughout this section, we always assume that system (2) undergoes Hopf bifurcation at the positive equilibrium for ,  , and then are corresponding purely imaginary roots of the characteristic equation at the positive equilibrium .

For convenience, let and , where is defined by (27) and , drop the bar for the simplification of notations, and then system (4) can be written as an FDE in as

where and ,  ,  , and and are given by

respectively, where and

From the discussion in Section 2, we know that if , then system (31) undergoes a Hopf bifurcation at the positive equilibrium and the associated characteristic equation of system (31) has a pair of simple imaginary roots .

By the representation theorem, there is a matrix function with bounded variation components ,  , such that

In fact, we can choose

where is the Dirac delta function.

For , define

Then (31) is equivalent to the abstract differential equation

where ,  

For , define

For and , define the bilinear form

where , , and are adjoint operators. By the discussions in Section 2, we know that are eigenvalues of , and they are also eigenvalues of .

Suppose that is the eigenvector of corresponding to ; then . It follows from definition and (32) and (35) that

Then we can obtain

Similarly, let be the eigenvector of corresponding to ; then . It follows from the definition that

Hence

In order to assure and , we need to determine the value of . From (40), we have

Thus we can choose

Next, we use the same notations as those in Hassard et al. [24] and we first compute the coordinates to describe the center manifold at . Let be the solution of (31) when .

Define

on the center manifold , and we have

where

and and are local coordinates for center manifold in the direction of and . Noting that is also real if is real, we consider only real solutions. For solutions of (31)

which we write in abbreviated form as

where

Hence we have

where

Thus we obtain

For unknown ,  ,  ,  ,   in , we still need to compute them.

In view of (38) and (47), we have

where

Comparing the coefficients, we obtain

We know that for

Comparing the coefficients of (60) with (57) gives that

From (58), (61), and the definition of , we get

Noting that , we have

where is a constant vector.

Similarly, from (59), (62), and the definition of , we have

where is a constant vector.

In what follows, we shall seek appropriate and   in (64) and (66), respectively. It follows from the definition of , (61), and (62) that

where . From (58), we have

Noting that

and substituting (64) and (69) into (67), we have

That is,

It follows that

where

Similarly, substituting (65) and (70) into (68), we have

That is,

It follows that

where

From (64), (66), (74), and (78), we can calculate and derive the following values:

These formulae give a description of the Hopf bifurcation periodic solutions of (31) at ,   on the center manifold. From the discussion above, we have the following result.

Theorem 3. The direction of the Hopf bifurcation is forward (backward) if ; the bifurcating periodic solutions are orbitally asymptotically stable with asymptotical phase (unstable) if ; the periods of the bifurcating periodic solutions increase (decrease) if

4. Numerical Examples

In this section, numerical simulations are carried out to investigate the effect of the time delay on dynamical behaviors of system (2) as well as illustrate our theoretical results. For convenience, we assume that some parametric values of system (2) are kept as ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  and . Then system (2) takes the form

which has a positive equilibrium and all the conditions given in Theorem 2 are satisfied. When , the positive equilibrium is asymptotically stable. Fix and by means of Matlab 7.0, we get ,  , and  – Thus we can derive the following values: –,  ,  ,  and Furthermore, it follows that and . Thus the positive equilibrium is stable when . When , we can easily plot the time series of every population and phase portrait of the system and find that the solution of (81) with initial value , , and would tend to a positive equilibrium (see Figures 1(a)–1(j), in detail); that is, the positive equilibrium is asymptotically stable. When passes through the critical value , the positive equilibrium loses its stability and a Hopf bifurcation occurs; that is, a family of periodic solutions bifurcate from the positive equilibrium . Figures 2(a)–2(j) are plotted by fixing the time delay and initial value ,  ,  and . It is shown that a Hopf bifurcation occurs from the positive equilibrium . In view of and , we find that the direction of the Hopf bifurcation is and these bifurcating periodic solutions from at are stable.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

Figure 1: Dynamical behavior of system (81) with . The positive equilibrium is asymptotically stable. The initial value .

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

Figure 2: Dynamical behavior of system (81) with . Hopf bifurcation occurs from the positive equilibrium . The initial value .

5. Biological Meanings and Conclusions

In this paper, the dynamics including the local stability of the positive equilibrium and the local Hopf bifurcation of a competitor-competitor-mutualist Lotka-Volterra model are investigated. It is showed that the positive equilibrium of system (2) is asymptotically stable for all when the conditions – are satisfied. This reveals that the densities of competing species and cooperating species will tend to stabilization; in other words, the densities of competing species and cooperating species will tend to , respectively, and this situation does not vary with the delay . When the delay increases, the positive equilibrium loses its stability and a sequence of Hopf bifurcations occur from the positive equilibrium ; that is, a family of periodic orbits bifurcate from the positive equilibrium . From this we can find that competing species and cooperating species can coexist and keep in periodic mode; that is, the densities of competing species and cooperating species will oscillate in the vicinity of , and , respectively. Finally, the direction of Hopf bifurcation and the stability of the bifurcating periodic orbits are discussed by applying the normal form theory and the center manifold theorem. A numerical example is also included to verify our theoretical findings.

Conflicts of Interest

The author declares that there are no conflicts of interest.

Acknowledgments

This work is supported by National Natural Science Foundation of China (no. 61673008 and no. 11261010) and Project of High-Level Innovative Talents of Guizhou Province ([2016]5651).

References

G. Aggelis, D. V. Vayenas, V. Tsagou, and S. Pavlou, “Prey-predator dynamics with predator switching regulated by a catabolic repression control mode,” Ecological Modelling, vol. 183, no. 4, pp. 451–462, 2005.

N. Bairagi and D. Jana, “On the stability and Hopf bifurcation of a delay-induced predator-prey system with habitat complexity,” Applied Mathematical Modelling. Simulation and Computation for Engineering and Environmental Systems, vol. 35, no. 7, pp. 3255–3267, 2011.

K. Chakraborty, M. Chakraborty, and T. K. Kar, “Bifurcation and control of a bioeconomic model of a prey-predator system with a time delay,” Nonlinear Analysis: Hybrid Systems, vol. 5, no. 4, pp. 613–625, 2011.

K. pada Das, K. Kundu, and J. Chattopadhyay, “A predator-prey mathematical model with both the populations affected by diseases,” Ecological Complexity, vol. 8, no. 1, pp. 68–80, 2011.

C.-H. Hsu, C.-R. Yang, T.-H. Yang, and T.-S. Yang, “Existence of traveling wave solutions for diffusive predator-prey type systems,” Journal of Differential Equations, vol. 252, no. 4, pp. 3040–3075, 2012.

J. Jiao, S. Cai, and L. Chen, “Analysis of a stage-structured predatory-prey system with birth pulse and impulsive harvesting at different moments,” Nonlinear Analysis: Real World Applications, vol. 12, no. 4, pp. 2232–2244, 2011.

T. K. Kar and A. Ghorai, “Dynamic behaviour of a delayed predator-prey model with harvesting,” Applied Mathematics and Computation, vol. 217, no. 22, pp. 9085–9104, 2011.

W. Ko and K. Ryu, “Coexistence states of a nonlinear Lotka-Volterra type predator-prey model with cross-diffusion,” Nonlinear Analysis: Theory, Methods & Applications, vol. 71, no. 12, pp. 1109–1115, 2009.

Q. M. Liu and R. Xu, “Periodic solution for a delayed three-species food-chain system with Holling type-II functional response,” International Journal of Mathematics and Mathematical Sciences, vol. 64, pp. 4057–4070, 2003.

Y. Saito, J. Sugie, and Y.-H. Lee, “Global asymptotic stability for predator-prey models with environmental time-variations,” Applied Mathematics Letters, vol. 24, no. 12, pp. 1973–1980, 2011.

Y. Song and J. Wei, “Local Hopf bifurcation and global periodic solutions in a delayed predator-prey system,” Journal of Mathematical Analysis and Applications, vol. 301, no. 1, pp. 1–21, 2005.

X. Lv, P. Yan, and S. Lu, “Existence and global attractivity of positive periodic solutions of competitor-competitor-mutualist Lotka-Volterra systems with deviating arguments,” Mathematical and Computer Modelling, vol. 51, no. 5-6, pp. 823–832, 2010.

M. Xiao and J.-D. Cao, “Stability and Hopf bifurcation in a delayed competitive web sites model,” Physics Letters A, vol. 353, no. 2-3, pp. 138–150, 2006.

C. Xu, X. Tang, and M. Liao, “Stability and bifurcation analysis of a delayed predator-prey model of prey dispersal in two-patch environments,” Applied Mathematics and Computation, vol. 216, no. 10, pp. 2920–2936, 2010.

C. Xu, X. Tang, M. Liao, and X. He, “Bifurcation analysis in a delayed Lokta-Volterra predator-prey model with two delays,” Nonlinear Dynamics, vol. 66, no. 1-2, pp. 169–183, 2011.

R. Xu and Z. E. Ma, “Stability and Hopf bifurcation in a ratio-dependent predator-prey system with stage structure,” Chaos, Solitons & Fractals, vol. 38, no. 3, pp. 669–684, 2008.

X.-P. Yan and W.-T. Li, “Bifurcation and global periodic solutions in a delayed facultative mutualism system,” Physica D. Nonlinear Phenomena, vol. 227, no. 1, pp. 51–69, 2007.

X.-P. Yan and C.-H. Zhang, “Hopf bifurcation in a delayed Lokta-Volterra predator-prey system,” Nonlinear Analysis: Real World Applications, vol. 9, no. 1, pp. 114–127, 2008.

C. B. Yu and J. J. Wei, “Stability and bifurcation analysis in a basic model of the immune response with delays,” Chaos, Solitons & Fractals, vol. 41, no. 3, pp. 1223–1234, 2009.

S. Yuan and F. Zhang, “Stability and global Hopf bifurcation in a delayed predator-prey system,” Nonlinear Analysis: Real World Applications, vol. 11, no. 2, pp. 959–977, 2010.

M. L. Zeeman, “Hopf bifurcations in competitive three-dimensional Lotka-Volterra systems,” Dynamics and Stability of Systems, vol. 8, no. 3, pp. 189–217, 1993.

X. Zhou, X. Shi, and X. Song, “Analysis of nonautonomous predator-prey model with nonlinear diffusion and time delay,” Applied Mathematics and Computation, vol. 196, no. 1, pp. 129–136, 2008.

S. G. Ruan and J. J. Wei, “On the zeros of transcendental functions with applications to stability of delay differential equations with two delays,” Dynamics of Continuous, Discrete & Impulsive Systems A: Mathematical Analysis, vol. 10, no. 6, pp. 863–874, 2003.

B. Hassard, D. Kazarino, and Y. Wan, Theory and Applications of Hopf Bifurcation, Cambridge University Press, Cambridge, UK, 1981.