Feedback control effect on the Lotka-Volterra prey-predator system with discrete delays

Advances in Difference Equations, Dec 2017

In this paper, we study a Lotka-Volterra prey-predator system with feedback control. We establish sufficient conditions under which a unique positive equilibrium is globally stable. Further, we show that a suitable feedback control on predator species can make prey species that is on the brink of extinction become globally stable, but under the conditions of small feedback control on predator, the prey species still extinct, whereas the predator species is stable at certain values. Several examples are presented to show the feasibility of the main results.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

https://advancesindifferenceequations.springeropen.com/track/pdf/10.1186/s13662-017-1410-0?site=advancesindifferenceequations.springeropen.com

Feedback control effect on the Lotka-Volterra prey-predator system with discrete delays

Shi et al. Advances in Difference Equations Feedback control effect on the Lotka-Volterra prey-predator system with discrete delays Chunling Shi 1 Xiaoying Chen 1 Yiqin Wang 0 0 Department of Science Research, Fuzhou Fujian Institute of Education , Fuzhou, Fujian 350001 , P.R. China 1 College of Zhicheng, Fuzhou University , Fuzhou, Fujian 350002 , P.R. China In this paper, we study a Lotka-Volterra prey-predator system with feedback control. We establish sufficient conditions under which a unique positive equilibrium is globally stable. Further, we show that a suitable feedback control on predator species can make prey species that is on the brink of extinction become globally stable, but under the conditions of small feedback control on predator, the prey species still extinct, whereas the predator species is stable at certain values. Several examples are presented to show the feasibility of the main results. Lotka-Volterra system; discrete delays; global stability; extinction; feedback controls - where x(t) is the prey population density at time t, x(t) is the predator population density at time t, r is the intrinsic growth rate of prey, r is the intrinsic growth rate for predator species, aii for i = ,  are intraspecific competition rates, and aij for i = j, i, j = ,  are interspecific competition rates for prey and predator species, τ is the time of catching prey, and τ is maturation delay of predator. The dynamical behavior of prey-predator system such as (.) has been investigated by many authors, and many excellent results concerned with permanence, extinction and persistence or uniform persistence, global stability, and almost periodic solutions are obtained (see, for example, [–]). In all prey-predator systems, there is a common competition between prey and predator species, and in general this competition is one sided, that is, the loss occurs only in the prey populations. In some cases, poor natural environment leads to a decrease in the birth rate of the prey population. Meanwhile, with the large predatory, the prey population is less and less and then tends to zero. Recently, many scholars have done research on the ecosystem with feedback controls (see [–] and the references therein). In particular, Gopalsamy and Weng [] introduced a feedback control variable into a two-species competitive system and discussed the existence of the globally attractive positive equilibrium of the system with feedback controls. Hu et al. [] considered the extinction of a nonautonomous Lotka-Volterra competitive system with pure delays and feedback controls, and by simulation they found that suitable feedback control variables can transform extinct species into permanent. Therefore, a natural and important question is that whether a proper control only on the predator population can make the extinct prey species become permanent, and the prey-predator can coexist in a certain pattern. Motivated by the above work, in this paper, we consider the following predator-prey Lotka-Volterra delay system with feedback control: where u(t) is the indirect control variable, ri > , τi ≥ , aij > , i, j = , , τ = max{τ, τ}, c > , e > , d > , x˙(t) = x(t) r – ax(t) – ax(t – τ) , x˙(t) = x(t) r + ax(t – τ) – ax(t) – cu(t) , u˙ (t) = –eu(t) + dx(t), xi(θ ) = φi(θ ), u(θ ) = ψ (θ ), θ ∈ [–τ , ], (.) (.) where φi(θ ) (i = , ) and ψ (θ ) are nonnegative and bounded continuous functions on [–τ , ]. The aim of this paper is, by using the method of multiple Lyapunov functionals [, ] and by developing a new analysis technique [], to obtain sufficient conditions under which a unique positive equilibrium is globally stable. This paper is organized as follows. In the next section, as preliminaries, some assumptions and lemmas are introduced. In Section , the main results of this paper are stated and proved. Finally, several examples together with their numerical simulations show the feasibility of the main results and the considerable effects of feedback controls to extinction of prey species. 2 Preliminaries Throughout this paper, we introduce the following hypotheses: (H) (H) (H) (H) (H) r > a , r a r < a , r a rr > aa+ ced , r < r a a + ced , a > a , a a a > a , a a aa+ ced < rr < aa < aa . tl→im∞ x(t) = r . a By the method of Lyapunov functions (see [, –]), one can show that if the coefficients of system (.) satisfy (H), then system (.) possesses a unique positive equilibrium (x¯, x¯) = ( araa–+araa , araa++raaa ), which is globally attractive, that is, all positive solutions of system (.) satisfy If the coefficients of system (.) satisfy (H), then system (.) is extinct, that is, all positive solutions of system (.) satisfy Now, we state the following lemmas, which are useful in the proof of the main results. Lemma . Suppose that assumption (H) holds. It is not difficult to verity that system (.) has a unique positive equilibrium x∗ = x∗ = e(ra – ra) + rcd , e(aa + aa) + cda e(ra + ra) e(aa + aa) + cda , u∗ = de x∗. Lemma . Let (x(t), x(t), u(t))T be a solution of system (.) with initial condition (.). Then (x(t), x(t), u(t))T is positive and bounded for all t ≥ . Proof Obviously, the solution (x(t), x(t), u(t))T of system (.) with initial condition (.) is positive for all t ≥ . By the first equation of system (.) we have x˙(t) ≤ x(t) r – ax(t) . By a standard comparison principle and basic ODE theory it follows that lim supt→∞ x(t) ≤ ar . Hence, for any ε >  sufficiently small, there exists T >  such that x(t) ≤ ar + ε, t > T. By the second equation of system (.) we have r + ε – ax(t) x˙(t) ≤ x(t) r + a a for all t > T + τ, which implies that lim supt→∞ x(t) ≤ r+aa(ar +ε) . Setting ε → , it follows that lim supt→∞ x(t) ≤ ra+ra . Therefore, for any ε >  small enough, there exists T > aa T + τ >  such that x(t) ≤ ra+ra + ε, t > T . By the third equation of system (.), aa u˙ (t) ≤ –eu(t) + d ra + ra + ε aa for all t > T , which means that lim supt→∞ u(t) ≤ de ( ra+ra + ε). Letting ε → , it follows that aa d ra+ra . This completes the proof of Lemma .. lim supt→∞ u(t) ≤ e aa 3 Main results Now, we give our main results. Theorem . Let (x(t), x(t), u(t))T be a solution of system (.). If condition (H) or (H) holds, then the unique positive equilibrium of system (.) is globally asymptotically stable, that is, tl→im∞ xi(t) = xi∗, i = , , ηi xi – xi∗ – xi∗ ln xxi∗i + cηd u – u∗ , where η = , and η is a positive constant to be determined. Note that system (.) can be rewritten as x˙(t) = x(t) –a x(t) – x∗ – a x(t – τ) – x∗ , x˙(t) = x(t) a x(t – τ) – x∗ – a x(t) – x∗ – c u(t) – u∗ , u˙ (t) = –e u(t) – u∗ + d x(t) – x∗ . Calculating the derivative of V(t) along the solution (x(t), x(t), u(t))T of system (.), we have V˙ (t) = x(t) – x∗ –a x(t) – x∗ – a x(t – τ) – x∗ + η x(t) – x∗ a x(t – τ) – x∗ – a x(t) – x∗ – c u(t) – u∗ + cdη u(t) – u∗ –e u(t) – u∗ + d x(t) – x∗ ≤ –a x(t) – x∗  – aη x(t) – x∗  – a x(t) – x∗ x(t – τ) – x∗ + ηa x(t) – x∗ x(t – τ) – x∗ . (.) θ a + θ b, θ > , from (.) it follows that By the inequality ab ≤  V˙ (t) ≤ –a x(t) – x∗  – aη x(t) – x∗  + a  θ θ x(t) – x∗  +  x(t – τ) – x∗  + ηa θ x(t – τ) – x∗  + θ x(t) – x∗  = –a x(t) – x∗  – aη x(t) – x∗  + aθ x(t) – x∗  + aθ x(t – τ) – x∗  + ηθa x(t – τ) – x∗  + ηaθ x(t) – x∗ . (.) (.) Let V(t) = ηθa t t–τ x(s) – x∗  ds + aθ  t t–τ x(s) – x∗  ds. Calculating the derivative of V(t), we obtain V˙ (t) = ηθa x(t) – x∗  – ηθa x(t – τ) – x∗  + aθ x(t) – x∗  – aθ x(t – τ) – x∗ .   Define V (t) = V(t) + V(t). It follows from (.) and (.) that V˙ (t) ≤ – a – aθ – ηa θ x(t) – x∗  – aη – aθ – ηaθ x(t) – x∗ . Denote δ = a – aθ – ηθa and δ = aη – aθ – ηaθ . Then taking η = aa and θ = θ = aa+aa , we have aa δ = a(aa – aa) , aa + aa δ = a(aa – aa) . aa Then, (H) or (H) shows that δi > , i = , . It is easy to see that V˙ (t) ≤ –δ x(t) – x∗  – δ x(t) – x∗ . Therefore, V (t) is nonincreasing. By Lemma ., x˙i(t), i = , , are bounded. So |xi(t) – xi∗|, i = , , are uniformly continuous on [, +∞). Integrating both sides of (.) on the interval [T , t), we have V(t) = ηθa t t–τ x(s) – x∗  ds + aθ  t t–τ x(s) – x∗  ds ≤ ηθa τ M + x∗  + aθ τ M + x∗  < +∞. Therefore xi(s) – xi∗  ds < V δ(Ti ) < +∞, i = , . x(s) – x∗  ds + δ x(s) – x∗  ds < V (T ). It follows from Lemma . and the initial condition φi that xi(t), i = , , are bounded for t ∈ R, that is, there exists M >  such that  < xi(t) < M, i = , , t ∈ R. Obviously, V(T ) is bounded, and (.) (.) (.) (.) (.) From this inequality it follows that (xi(s) – xi∗) ∈ L[, +∞), i = , . By Barbalat’s lemma (see []) we conclude that lim xi(t) – xi∗  = , i = , , t→∞ and therefore tl→im∞ xi(t) = xi∗, i = , . Corollary . Let (x(t), x(t), u(t))T be a solution of system (.). If condition (H) holds, then the unique positive equilibrium of system (.) is globally asymptotically stable, that is, tl→im∞ xi(t) = xi∗, i = , , Remark . If (H) holds, then systems (.) and (.) are globally stable. Theorem . implies that the feedback control keeps the property of stability of system (.) but only changes the position of the unique positive equilibrium. That is, feedback control of system (.) leads to the number of the prey population increased (x∗ = e(ae(raa+–araa))++rcdcad > x¯ = araa–+araa ) and the number of the predator population decreased (x∗ = e(aa+aa)+cda < x¯ = araa++raaa ). e(ra+ra) Remark . If (H) holds, system (.) is extinct. (H) implies (H), and system (.) has a unique positive equilibrium, which is globally asymptotically stable. Theorem . implies that the proper feedback control on predator species can change extinct prey species to be permanent. Theorem . Assume that (H) holds. Let (x(t), x(t), u(t))T be a solution of system (.). Then By (.) we can choose positive constants α > , β > , γ >  such that rr < βα < aa (.) and ec < γβ < αad–ββa < ard–rar . Thus there exists δ >  such that –βr + αr < –δ < , βa + γ d – αa < , βc – γ e < . (.) Consider a Lyapunov functional of the form V (t) = xα(t)x–β (t) exp γ u(t) – βa x(s) ds – αa x(s) ds . t t–τ t t–τ Calculating the derivative of V along the solution of system (.), we have V˙ (t) = V (t) +α r – ax(t) – ax(t – τ) – β r + ax(t – τ) – ax(t) – cu(t) + γ –eu(t) + dx(t) – βax(t) + βax(t – τ) – αax(t) + αax(t – τ) = V (t) –βr + αr + (–aα – βa)x(t) + (βa + γ d – αa)x(t) + (βc – γ e)u(t) . V˙ (t) ≤ –δV (t). V (t) ≤ V () exp(–δt). From inequalities (.) and (.) we obtain Integrating this inequality from  to t, we have By similar arguments as in the proof of Theorem ., there exists M >  such that  < xi(t) < M, i = , , t ∈ R. So V () = xα()x–β () exp γ u() – βa x(s) ds – αa x(s) ds  –τ  –τ ≤ xα()x–β () exp γ u() + βaMτ + αaMτ ≤ +∞. On the other hand, V (t) ≥ xα(t)x–β (t) exp –βa x(s) ds – αa ≥ xα(t)x–β (t) exp[–βaMτ – αaMτ ]. t t–τ t t–τ x(s) ds (.) (.) (.) (.) Combining inequalities (.), (.), and (.), we have xα(t) ≤ xβ (t) exp(βaMτ + αaMτ )V () exp(–δt). (.) Inequality (.) means that δ x(t) ≤ λ exp – α t ,  where λ = [Mβ exp(βaMτ + αaMτ )V ()] α . Hence, we obtain that Theorem . If condition (H) holds, then the equilibrium (x∗∗, x∗∗, u∗∗) = (, eaer+cd , eadr+cd ) of system (.) is globally asymptotically stable, that is, Proof Let (x(t), x(t), u(t))T be a solution of system (.). It follows from Theorem . that limt→∞ x(t) = , and it is easy to obtain that system (.) has an equilibrium (x∗∗, x∗∗, u∗∗) = (, ea+cd , eadr+cd ). Therefore, we only need to verify er Calculating the derivative of V(t) along the solution of system (.), we obtain that V(t) = x(t) – x∗∗ ax(t – τ) – a x(t) – x∗∗ – c u(t) – u∗∗ + dc u(t) – u∗∗ –e u(t) – u∗∗ + d x(t) – x∗∗ ≤ –a x(t) – x∗∗  + a M + x∗∗ x(t – τ) – Calculating the derivative of V(t) along the solution of system (.), we have V(t) ≤ a M + x∗∗ x(t) – a M + x∗∗ x(t – τ). (.) (.) (.) Define V (t) = V(t) + V(t). It follows from (.) and (.) that V (t) ≤ –a x(t) – x∗∗  + a M + x∗∗ x(t) – Integrating both sides of (.) on the interval [T , t), we obtain V (t) + a T t x(s) – x∗∗  ds + ce d T t u(s) – u∗∗  ds ≤ V (T ) + a M + x∗∗ (.) T t tl→im∞ x(t) – x∗∗  = , tl→im∞ u(t) – u∗∗  = , that is, Remark . By comparative analysis of (H) and (H) note that when the feedback control repression ced remains small, feedback control on predator species has no influence on the extinction of system (.). 4 Examples In this section, we give examples to illustrate the results obtained. Example . We consider the following equations: x˙(t) = x(t)  – x(t) – x(t – ) , x˙(t) = x(t)  + x(t – ) – x(t) . (.) Consider the initial conditions (x(θ ), x(θ ))T = (., .)T and (, .)T for all θ ∈ [–, ] and t ∈ [, ]. Obviously, rr =  > aa =  , aa =  > aa =  . Then condition (H) holds, and there exists a unique positive equilibrium (x∗, x∗) = (  ,  ), which is globally asymptotically stable (see Figure ). Example . We introduce feedback control to the predator species of system (.): x˙(t) = x(t)  – x(t) – x(t – ) , x˙(t) = x(t)  + x(t – ) – x(t) – u(t) , Consider the initial conditions (x(θ ), x(θ ), u(θ ))T = (., ., .)T and (, ., )T for all θ ∈ [–, ] and t ∈ [, ]. For system (.), condition (H) holds; then system (.) has a unique positive equilibrium (x∗, x∗, u∗) = (  ,  ,  ), which is globally asymptotically  stable. Therefore, the feedback control variable can change the position of equilibrium and retain the stable property (Figure ). Example . We consider the equation x˙(t) = x(t)  – x(t) – x(t – ) , x˙(t) = x(t)  + x(t – ) – x(t) . Consider the initial conditions (x(θ ), x(θ ))T = (., .)T and (, .)T for all θ ∈ [–, ] and t ∈ [, ]. It is easy to see that rr =  < aa =  , so that condition (H) holds. Then system (.) has an equilibrium (x∗∗, x∗∗) = (,  ), which means that the prey species will be extinct whereas the predator species will be globally asymptotically stable (see Figure ). (.) (.) x˙(t) = x(t)  – x(t) – x(t – ) , x˙(t) = x(t)  + x(t – ) – x(t) – u(t) , Consider the initial conditions (x(θ ), x(θ ), u(θ ))T = (., ., .)T and (, ., )T for all θ ∈ [–, ] and t ∈ [, ]. By calculation, rr =  , a+ ced =  , aa =  , aa =  . a+ ced < a a r < aa < aa , so that condition (H) or (H) holds. Therefore, system (.) has a unique r positive equilibrium (x∗, x∗, u∗) = (  ,  ,  ), which is globally asymptotically stable. The result shows that the proper feedback control on predator species can change extinct prey species to be permanent, and the label shows that the number of the prey species has increased, however, is larger than that of the predator species (see Figure ). Example . Finally, we consider the system x˙(t) = x(t)  – x(t) – x(t – ) , x˙(t) = x(t)  + x(t – ) – x(t) – u(t) , Consider the initial conditions (x(θ ), x(θ ), u(θ ))T = (., ., .)T and (, ., )T for all θ ∈ [–, ] and t ∈ [, ], so that condition (H) holds. System (.) has an equilibrium (x∗∗, x∗∗, u∗∗) = (,  ,  ), which shows that the prey species is extinct whereas the predator species still has a positive equilibrium, which is globally asymptotically stable (see Figure ). (.) (.) Figure 5 Dynamic behavior of model (4.5) with the initial value (x1, x2, u) = (0.3, 0.2, 0.3) and (1, 0.8, 1). Acknowledgements This work is supported by the National Natural Science Foundation of Fujian Province (2015J01012) and the Foundation of Fujian Education Bureau (JA13361). Competing interests The authors declare that they have no competing interests. Authors’ contributions All authors read and approved the final manuscript. Publisher’s Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. 1. Saito , Y , Hara, T , Ma, W: Necessary and sufficient conditions for permanence and global stability of a Lotka-Volterra system with two delays . J. Math. Anal. Appl . 236 , 534 - 556 ( 1999 ) 2. Xu , R , Chen, LS : Persistence and global stability for a delayed nonautonomous predator-prey system without dominating instantaneous negative feedback . J. Math. Anal. Appl . 262 , 50 - 61 ( 2001 ) 3. Muroya , Y: Permanence and global stability in a Lotka-Volterra predator-prey system with delays . Appl . Math. Lett. 16 , 1245 - 1250 ( 2003 ) 4. Song , X , Chen , L : Persistence and global stability for nonautonomous predator-prey system with diffusion and time delay . Comput. Math. Appl . 35 , 33 - 40 ( 1998 ) 5. Cui , J: The effect of dispersal on permanence in a predator-prey population growth model . Comput. Math. Appl . 44 , 1085 - 1097 ( 2002 ) 6. Chen , FD , Chen, LJ , Xie, XD: On a Leslie-Gower predator-prey model incorporating a prey refuge . Nonlinear Anal., Real World Appl . 10 , 2905 - 2908 ( 2009 ) 7. Chen , FD , Ma, ZZ , Zhang, HY: Global asymptotical stability of the positive equilibrium of the Lotka-Volterra prey-predator model incorporating a constant number of prey refuges . Nonlinear Anal., Real World Appl . 13 , 2790 - 2793 ( 2012 ) 8. Chen , FD , Xie, X , Miao , Z , Pu, L: Extinction in two species nonautonomous nonlinear competitive system . Appl. Math. Comput . 274 , 119 - 124 ( 2016 ) 9. Ahmad , S , Stamova, IM: Partial persistence and extinction in N-dimensional competitive systems . Nonlinear Anal . 60 , 821 - 836 ( 2005 ) 10. Zu , L , Jiang, DQ, O 'Regan , D , Ge , B: Periodic solution for a non-autonomous Lotka-Volterra predator-prey modal with random perturbation . J. Math. Anal. Appl . 430 , 428 - 437 ( 2015 ) 11. Li , Z , Chen, FD , He, MX: Permanence and global attractivity of a periodic predator-prey system with mutual interference and impulses . Commun. Nonlinear Sci. Numer . Simul. 17 , 444 - 453 ( 2012 ) 12. Teng , ZD , Yu, YH: The extinction in nonautonomous prey-predator Lotka-Volterra systems . Acta Math. Appl. Sin . 15 ( 4 ), 401 - 408 ( 1999 ) 13. Gomez , L , Ortega, R: The periodic predator-prey Lotka-Volterra models . Adv. Differ. Equ . 1 , 403 - 423 ( 1996 ) 14. Kuang , Y: Delay Differential Equations with Applications in Population Dynamics . Academic Press, New York ( 1993 ) 15. Hou , Z: Vanishing components in autonomous competitive Lotka-Volterra systems . J. Math. Anal. Appl . 359 , 302 - 310 ( 2009 ) 16. Kar , TK , Ghorai, A : Dynamic behaviour of a delayed predator-prey model with harvesting . Appl. Math. Comput . 217 , 9085 - 9104 ( 2011 ) 17. Jana , S , Kar, TK: A mathematical study of a prey-predator model in relevance to pest control . Nonlinear Dyn . 74 , 667 - 683 ( 2013 ) 18. Liu , M , Mandal, PS: Dynamical behavior of a one-prey two-predator model with random perturbations . Commun. Nonlinear Sci. Numer . Simul. 28 , 123 - 137 ( 2015 ) 19. Huo , HF , Li, WT : Positive periodic solutions of a class of delay differential system with feedback control . Appl. Math. Comput . 148 , 35 - 46 ( 2004 ) 20. Muroya , Y: Global stability of a delayed nonlinear Lotka-Volterra system with feedback controls and patch structure . Appl. Math. Comput . 239 , 60 - 73 ( 2014 ) 21. Weng , PX: Existence and global stability of positive periodic solution of periodic integrodifferential systems with feedback controls . Comput. Math. Appl . 40 , 747 - 759 ( 2000 ) 22. Gopalsamy , K , Weng, PX: Global attractivity in a competition system with feedback controls . Comput. Math. Appl . 45 , 665 - 676 ( 2003 ) 23. Shi , CL , Li , Z , Chen, FD: Extinction in nonautonomous Lotka-Volterra competitive system with infinite delay and feedback controls . Nonlinear Anal., Real World Appl . 13 , 2214 - 2226 ( 2012 ) 24. Li , Z , Han, M , Chen, FD : Influence of feedback controls on an autonomous Lotka-Volterra competitive system with infinite delays . Nonlinear Anal., Real World Appl . 14 , 402 - 413 ( 2013 ) 25. Nie , LF , Teng, ZD , Hu, L , Peng, JG : Permanence and stability in non-autonomous predator-prey Lotka-Volterra system with feedback controls . Comput. Math. Appl . 58 , 436 - 448 ( 2009 ) 26. Hu , HX , Teng, ZD , Gao, SJ: Extinction in a nonautonomous Lotka-Volterra competitive system with pure-delays and feedback controls . Nonlinear Anal., Real World Appl . 10 , 2508 - 2520 ( 2009 ) 27. Zeeman , ML : Extinction in competitive Lotka-Volterra systems . Proc. Am. Math. Soc . 123 , 87 - 96 ( 1995 ) 28. Liu, SQ , Chen, LS , Luo, GL , Jiang, YL : Asymptotic behaviors of competitive Lotka-Volterra system with stage structure . J. Math. Anal. Appl . 271 , 124 - 138 ( 2002 ) 29. He , M , Chen , F , Li , Z: Permanence and global attractivity of an impulsive delay logistic model . Appl . Math. Lett. 62 , 92 - 100 ( 2016 ) 30. Gopalsamy , K: Stability and Oscillations in Delay Different Equations of Population Dynamics . Kluwer Academic, Dordrecht ( 1992 )


This is a preview of a remote PDF: https://advancesindifferenceequations.springeropen.com/track/pdf/10.1186/s13662-017-1410-0?site=advancesindifferenceequations.springeropen.com

Chunling Shi, Xiaoying Chen, Yiqin Wang. Feedback control effect on the Lotka-Volterra prey-predator system with discrete delays, Advances in Difference Equations, 373,