Abstract and Applied Analysis

Volume 2014 (2014), Article ID 104156, 11 pages

http://dx.doi.org/10.1155/2014/104156

Research Article

Control of Hopf Bifurcation and Chaos in a Delayed Lotka-Volterra Predator-Prey System with Time-Delayed Feedbacks

Huitao Zhao,1,2 Yaowei Sun,¹ and Zhen Wang³

¹Department of Mathematics and Information Science, Zhoukou Normal University, Zhoukou, Henan 466001, China

²Department of Applied Mathematics, Kunming University of Science and Technology, Kunming, Yunnan 650093, China

³Institute of Information and System Computation Science, Beifang University of Nationalities, Yinchuan, Ningxia 750021, China

Received 10 March 2014; Accepted 19 May 2014; Published 4 June 2014

Academic Editor: Imran Naeem

Copyright © 2014 Huitao Zhao 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.

Abstract.

A delayed Lotka-Volterra predator-prey system with time delayed feedback is studied by using the theory of functional differential equation and Hassard’s method. By choosing appropriate control parameter, we investigate the existence of Hopf bifurcation. An explicit algorithm is given to determine the directions and stabilities of the bifurcating periodic solutions. We find that these control laws can be applied to control Hopf bifurcation and chaotic attractor. Finally, some numerical simulations are given to illustrate the effectiveness of the results found.

1. Introduction

Lotka-Volterra system is one of the most classical and important systems in the field of mathematical biology. Since the word of Volterra, there have been extensively detailed investigations on Lotka-Volterra system including stability, attractivity, persistence, periodic oscillation, bifurcation and chaos (see [1–6] and the references therein). In particular, the properties of periodic solutions arising from the Hopf bifurcation are of great interest [7–10]. But the study on chaos control of Lotka-Volterra system is scarce.

Reference [3] and the references therein proposed that, for a two-species competition system with delays

when is big enough, the chaotic behavior may occur. For example, Yan and Zhang [9] investigated the following delayed prey-predator system with a single delay:

where , , , , , and are all positive constants. The delay denotes the gestation period of the predator. Their results show that, taking as the bifurcation parameter, when passes through a certain critical value, the positive equilibrium loses its stability and Hopf bifurcation takes place. Furthermore, when takes a sequence of critical values containing the above critical value, the positive equilibrium of system (2) will undergo a Hopf bifurcation. With the further increase of the delay, the system will show the chaotic phenomenon (see Figure 1).

Figure 1: Waveform plot and phase plot of system (2) with .

In the sense of biology, chaotic behavior sometimes is to the disadvantage of virtuous cycle and develop of the ecosystem, so we want to control this chaos phenomenon and create periodic orbits. So far, many researchers have proposed chaos control schemes in recent years [11–16]. For example, Song and Wei in [17] investigated the chaos phenomena of Chen’s system using the method of delayed feedback control. Their results show that, when the controlling parameter to be some value, taking the delay as the bifurcation parameter, when passes through a certain critical value, the stability of the equilibrium will be changed from unstable to stable, chaos vanishes, and a periodic solution emerges.

To the end of controlling chaos in system (2), stimulating by the works of above, we add some delayed feedback terms to system (2), that is, the following delayed feedback control system:

where denote the capture coefficient when (or release coefficient when ). By choosing and as bifurcation parameter, we get the conditions under which Hopf bifurcation occurs. And then, we derive the explicit algorithm for determining the direction of the Hopf bifurcation and stability of the bifurcating periodic solutions. At last we will give some example showing that when is fixed, with increasing, the stability of the positive equilibrium will be changed, chaos vanishes, a periodic solution occurs.

This paper is organized as follows. In Section 2, we first focus on the stability and Hopf bifurcation of the positive equilibrium. In Section 3, we derive the direction and stability of Hopf bifurcation by using normal form and central manifold theory. Finally in Section 4, numerical simulations are performed to support the stability results.

2. Stability and Hopf Bifurcation Analysis with Delayed Feedback Control

In this section, by analyzing the characteristic equation of the linearized system of system (3) at the positive equilibrium, we investigate the stability of the positive equilibrium and the existence of the local Hopf bifurcations occurring at the positive equilibrium. To guarantee that system (3) has always a positive equilibrium, throughout this section, we assume that the coefficients of system (3) satisfy the following condition:

(H1).

Clearly, under the hypothesis (H1), system (3) has a unique positive equilibrium , where

Let , ; then system (3) can be rewritten as the following equivalent system:

Thus, the positive equilibrium of system (3) is transformed into the equilibrium of system (5). Linearizing system (5) at the equilibrium yields the following linear system:

The characteristic equation of system (6) is

where , , , , . Multiplying by both sides of (8), we have

Thus, is a root of (8) if and only if satisfies the following equation:

Separating the real and imaginary parts, we have

It follows that

where , , , , , . From , we have

where , , , .

Denote ; (12) becomes

Let

Since , we conclude that if , then (13) has at least one positive root. If all the parameters of system (5) are given, it is easy to get the roots of (13) by using a computer. Suppose the following.

(H2)Equation (13) has at least one positive real root. Without loss of generality, we assume that (14) has four positive roots, defined by , , , , respectively. Then (12) has four positive roots

From (11), we have

If we denote

where , then is a pair of purely imaginary roots of (8) with .

Define

Let be the root of (8) near satisfying , . Substituting into (8) and taking the derivative with respect to , we have

It follows that

Then

Through tedious calculating, we can get

where

Denote ; then , and we have

and note that

In order to get the main results, it is necessary to make the following assumption.

(H3). Then, if (H3) holds, the transversality conditions hold.

Note that when , (8) becomes

From [9], we know that all the roots of (26) have negative real parts; hence, the positive equilibrium is locally asymptotically stable for . Then, we can employ a result of Ruan and Wei [18] to analyze (8). For the convenience of the reader, we state it as follows.

Lemma 1. Consider the exponential polynomial

where and are constants. As vary, the sum of the order of the zeroes of on the open right half plane can change only if a zero appears on or crosses the imaginary axis.

From Lemma 1 and the above assumption, we can obtain the following theorem.

Theorem 2. Suppose that and hold; then the following results hold true.

(i)The positive equilibrium of system (3) is asymptotically stable for ;(ii)system (3) exhibits Hopf bifurcation at the positive equilibrium for .

3. Direction and Stability of the Hopf Bifurcation

In this section, we obtain the conditions under which a family periodic solutions bifurcate form the steady state at the critical value of . As pointed out by Hale and Verduyn Lunel [19] and Hassard et al. [20], it is interesting to determine the direction, stability, and period of these periodic solutions bifurcating from the steady state. Following the ideal of [20], we derive the explicit formulae for determining the properties of the Hopf bifurcation at the critical value of using the normal form and the center manifold theory.

For the sake of simplicity of notation, we denote the critical values as , and when , we denote the pair of purely imaginary roots of (8) as . Let ; then is the Hopf bifurcation value of system (5). In the following, we consider the equivalent system (6). Let ; then the system (5) can be rewritten as a functional differential equation in :

where , For ,

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 denotes Dirac-delta function. For , define

Then when , the system

is equivalent to the system (29), where , . For , define

and a bilinear inner product

where , and 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 . It follows from the definition of , , and that

Thus, we can easily obtain , .

Similarly, let be the eigenvector of corresponding to . By the definition of , we can compute .

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

Thus, we can choose

such that , .

In the following, we first compute the coordinates to describe the center manifold at . 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 , since , we have

where

From (40) and (41), we have

In addition, ; then

By the definition of , we have

Substituting , , , and into the above equation and comparing the coefficients with (43), we get

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

where

Notice that near the origin on the center manifold , we have

thus, we have

From (48), for , we have

Comparing the coefficients with (48) gives that

From (51), (53), and the definition of , we can get

Notice that ; we have

where is a constant vector. In the same way, we can also obtain

where is also a constant vector.

In what follows, we will compute and . From the definition of and (51), we have

From (48) and (49), we have

Substituting (55) and (59) into (57) and noticing that

we obtain

which leads to

where

It follows that

Similarly, substituting (56) and (60) into (58), we can get the formula of , where

Thus, we can determine and . Furthermore, we can determine each . Therefore, each is determined by the parameters and delay in (5). Thus, we can compute the following values:

which determine the quantities of bifurcating periodic solutions in the center manifold at the critical value ; that is, determines the directions of the Hopf bifurcation: if , then the Hopf bifurcation is supercritical (subcritical) and the bifurcation exists for ; determines the stability of the bifurcation periodic solutions: the bifurcating periodic solutions are stable (unstable) if ; and determines the period of the bifurcating periodic solutions: the period increases (decrease) if .

4. Numerical Simulations

In this section, we present numerical results to verify the analytical predictions obtained in the previous sections and use the delayed-feedback controller to control the Hopf bifurcation and chaos of system (3).

For the convenience of the calculation, we take the parameters of system (3) as , , , , . Then system (3) has a positive equilibrium .

From [9], when or , system (3) becomes system (2), which has a bifurcation point , and when , is asymptotically stable and is unstable when , and with the increase of the delay , chaos occurs via periodic-doubling bifurcation (Figure 2(a)), when , , system (3) becomes a chaotic one like illustrated by Figure 1.

(a)

(b)

Figure 2: The bifurcation diagram of : (a) ; (b) , .

Under the delayed-feedback control, if we choose , , from the algorithm of Section 2, we can compute that and ; hence from Theorem 2, when , is asymptotically stable and is unstable when , and when passes through the critical point , a family of periodic orbits will bifurcate from . Furthermore, from the algorithm of Section 3, we can compute that , , , ; hence, based on the formula of the end of Section 3, we have . It is easy to see that , . This means that the Hopf bifurcation is supercritical and the bifurcating periodic solutions are stable (Figure 5(b)).

By comparing the two bifurcation plots of Figure 2(a) and Figure 2(b), we observe that the bifurcation occurrence is postponed as the system is being controlled. Figure 3 shows that when , if , is unstable, and there are periodic orbits bifurcating from (Figure 3(a)); if we take , , becomes stable (Figure 3(b)). Figure 4 shows that when , if , is unstable, and system (3) becomes chaotic (Figure 4(a)); if we take , , then becomes stable and chaos vanishes (Figure 4(b)). Figure 5 shows that when , if , is unstable, and system (3) still exhibits chaotic behavior (Figure 5(a)); if we take , , then becomes unstable and stable periodic orbits bifurcate from (Figure 5(b)).

(a)

(b)

Figure 3: Waveform plot and phase plot of system (3) with : (a) ; (b) , .

(a)

(b)

Figure 4: Waveform plot and phase plot of system (3) with : (a) ; (b) , .

(a)

(b)

Figure 5: Waveform plot and phase plot of system (3) with : (a) ; (b) , .

5. Conclusion

In this paper, we have studied a delayed Lotka-Volterra predator-prey system with time delayed feedback by using the theory of functional differential equation and Hassard’s method. By analyzing the corresponding characteristic equations, the local stability of the positive equilibrium of system (3) was discussed.

We have obtained the estimated length of gestation delay which would not affect the stable coexistence of both prey and predator species at their equilibrium values. The existence of Hopf bifurcation for system (3) at the positive equilibrium was also established. From theoretical analysis it was shown that the larger values of gestation time delay cause fluctuation in individual population density and hence the system becomes unstable.

As the estimated length of delay to preserve stability and the critical length of time delay for Hopf-bifurcation are dependent upon the system parameters, it is possible to impose some control, which will prevent the possible abnormal oscillation in population density. Our results show that if we choose some appropriate parameters, the oscillation can be controlled to a stable equilibrium or a stable periodic orbit; that is to say, we can achieve the ecological equilibrium by adjusting the capture (or release) level.

Conflict of Interests

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

Acknowledgments

This research is supported by Science and Technology Department of Henan Province (no.122300410417), Education Department of Henan Province (no.13A110108), 2013 Narure Science Foundation of Ningxia (no. NZ13096), 2013 Higher educational scientific research project of Ningxia (no. NGY2013086).

References

J. Zhen and Z. Ma, “Stability for a competitive Lotka-Volterra system with delays,” Nonlinear Analysis. Theory, Methods & Applications, vol. 51, no. 7, pp. 1131–1142, 2002.

Y. Saito, “The necessary and sufficient condition for global stability of a Lotka-Volterra cooperative or competition system with delays,” Journal of Mathematical Analysis and Applications, vol. 268, no. 1, pp. 109–124, 2002.

K. Gopalsamy and X.-Z. He, “Oscillations and convergence in an almost periodic competition system,” Acta Applicandae Mathematicae, vol. 46, no. 3, pp. 247–266, 1997.

X. H. Tang and X. Zou, “Global attractivity of non-autonomous Lotka-Volterra competition system without instantaneous negative feedback,” Journal of Differential Equations, vol. 192, no. 2, pp. 502–535, 2003.

Y. Song and Y. Peng, “Stability and bifurcation analysis on a logistic model with discrete and distributed delays,” Applied Mathematics and Computation, vol. 181, no. 2, pp. 1745–1757, 2006.

H. Zhao and Y. Lin, “Hopf bifurcation in a partial dependent predator-prey system with delay,” Chaos, Solitons and Fractals, vol. 42, no. 2, pp. 896–900, 2009.

Y. Song, Y. Peng, and J. Wei, “Bifurcations for a predator-prey system with two delays,” Journal of Mathematical Analysis and Applications, vol. 337, no. 1, pp. 466–479, 2008.

X.-P. Yan, “Stability and Hopf bifurcation for a delayed prey-predator system with diffusion effects,” Applied Mathematics and Computation, vol. 192, no. 2, pp. 552–566, 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.

J. Zhang, Z. Jin, J. Yan, and G. Sun, “Stability and Hopf bifurcation in a delayed competition system,” Nonlinear Analysis. Theory, Methods & Applications, vol. 70, no. 2, pp. 658–670, 2009.

H. Chen, “Chaos control and global synchronization of Liu chaotic systems using linear balanced feedback control,” Chaos, Solitons and Fractals, vol. 40, no. 1, pp. 466–473, 2009.

W. Chen, “Dynamics and control of a financial system with time-delayed feedbacks,” Chaos, Solitons and Fractals, vol. 37, pp. 1198–1207, 2008.

D. Ghosh, A. R. Chowdhury, and P. Saha, “Bifurcation continuation, chaos and chaos control in nonlinear Bloch system,” Communications in Nonlinear Science and Numerical Simulation, vol. 13, no. 8, pp. 1461–1471, 2008.

A. El-Gohary and A. S. Al-Ruzaiza, “Chaos and adaptive control in two prey, one predator system with nonlinear feedback,” Chaos, Solitons and Fractals, vol. 34, no. 2, pp. 443–453, 2007.

J. Lü and S. Zhang, “Controlling Chen’s chaotic attractor using backstepping design based on parameters identification,” Physics Letters A, vol. 286, pp. 148–152, 2001.

S. Boccaletti, C. Grebogi, Y.-C. Lai, H. Mancini, and D. Maza, “The control of chaos: theory and applications,” Physics Reports, vol. 329, no. 3, pp. 103–197, 2000.

Y. Song and J. Wei, “Bifurcation analysis for Chen's system with delayed feedback and its application to control of chaos,” Chaos, Solitons & Fractals, vol. 22, no. 1, pp. 75–91, 2004.

S. Ruan and 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.

J. K. Hale and S. M. Verduyn Lunel, Introduction to Functional-Differential Equations, vol. 99 of Applied Mathematical Sciences, Springer, New York, NY, USA, 1993.

B. D. Hassard, N. D. Kazarinoff, and Y. H. Wan, Theory and Applications of Hopf Bifurcation, vol. 41 of London Mathematical Society Lecture Note Series, Cambridge University Press, New York, NY, USA, 1981.