DDNS Discrete Dynamics in Nature and Society 1607-887X 1026-0226 Hindawi Publishing Corporation 218785 10.1155/2012/218785 218785 Research Article A Lyapunov Function and Global Stability for a Class of Predator-Prey Models Wang Xiaoqin 1 Ma Huihai 2 Liu Junli 1 Faculty of Science Shaanxi University of Science and Technology Xi'an 710021 China sust.edu.cn 2 College of Electrical and Information Engineering Shaanxi University of Science and Technology Xi'an 710021 China sust.edu.cn 2012 30 12 2012 2012 29 11 2012 14 12 2012 2012 Copyright © 2012 Xiaoqin Wang and Huihai Ma. 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 construct a new Lyapunov function for a class of predation models. Global stability of the positive equilibrium states of these systems can be established when the Lyapunov function is used.

1. Introduction

The dynamics of predator-prey systems are often described by differential equations, which represent time continuously. A common framework for such a model is  (1.1)dNdt=Nf(N)-Pg(N,P),dPdt=h[g(N,P),P]P, where N and P are prey and predator densities, respectively, f(N) is the prey growth rate, g(N) is the functional response, for example, the prey consumption rate by an average single predator, and h[g(N),P] the per capita growth rate of predators (also known as the “predator numerical response"), which obviously increases with the prey consumption rate. The most widely accepted assumption for the numerical response with predator density restricting is as follows: (1.2)h[g(N,P),P]=εg(N,P)-δP-β, where β is a per capita predator death rate, ε the conversion efficiency of food into offspring, δ the density dependent rate . And prototype of the prey growth rate f(N) is the logistic growth (1.3)f(N)=r(1-NK), where K>0 is the carrying capacity of the prey.

When (1.4)g(N,P)=1-e-aN, where a is the efficiency of predator capture of prey, model (1.1) is called Ivlev-type predation model, due originally to Ivlev . And Ivlev-type functional response is classified to the prey-dependent; that is, g is independent of predator P .

Both ecologists and mathematicians are interested in the Ivlev-type predator-prey model and much progress has been seen in the study of the model . Of them, Xiao  gave global analysis of the following model: (1.5)dNdt=rN(1-N)-(1-e-aN)P,dPdt=P(1-e-aN-d-δP).

But, in paper , the author gave complex process to prove the global asymptotical stability of the positive equilibrium.

In this paper, we will establish a new Lyapunov function to prove the global stability of the positive equilibrium of model (1.5).

Our paper is organized as follows. In the next section, we discuss the existence, uniqueness of the positive equilibrium, and establish a new Lyapunov function to model (1.5). In Section 3, we will give some examples to show the robustness of our Lyapunov function.

2. Main Results

First of all, it is easy to verify that model (1.5) has two trivial equilibria (belonging to the boundary of R+2, that is, at which one or more of populations has zero density or is extinct), namely, E0=(0,0) and E1=(1,0). For the positive equilibrium, set (2.1)rN(1-N)-(1-e-aN)P=0,P(1-e-aN-δP-d)=0, which yields (2.2)N(1-N)=1rδ((1-e-aN)(1-e-aN-d)).

We have the following Lemma regarding the existence of the positive equilibrium.

Lemma 2.1 (see [<xref ref-type="bibr" rid="B8">8</xref>]).

Suppose 1-e-a>d. Model (1.5) has a unique positive equilibrium E*=(N*,P*) if either of the following inequalities holds:

d2(1-2e-a/2);

d<2(1-2e-a/2),    (a/rδ)(2(1-2e-a/2)-d))e-a/21+(2/a)ln(1/2-d/4) and (1/rδ)((1/2+d/4)2-d(1/2+d/4))<-(1/a)ln(1/2+d/4)-1/a2  .

Lemma 2.2.

Let 1-e-a>d, then Ω={(N,P)  0N1,  0P(1-e-a-d)/δ} is a region of attraction for all solutions of model (1.5) initiating in the interior of the positive quadrant R+2.

Proof.

Let (N(t),P(t)) be any solution of model (1.5) with positive initial conditions. Note that dN/dtN(1-N), by a standard comparison argument, we have (2.3)limtsupN(t)1. Then, (2.4)dPdt=P(1-e-aN-d-δP)P(1-e-a-d-δP). Similarly, since 1-e-a>d, we have (2.5)limtsupP(t)1-e-a-dδ. On the other hand, for all (N,P)Ω, we have dN/dt|N=0=0 and dP/dt|P=0=0. Hence, Ω is a region of attraction. As a consequence, we will focus on the stability of the positive equilibrium E* only in the region Ω.

In the following, we devote to the global stability of the positive equilibrium E*=(N*,P*) for model (1.5) by constructing a new Lyapunov function which is motivated by the work of Hsu .

Theorem 2.3.

If a2, the positive equilibrium E*=(N*,P*) of model (1.5) is globally asymptotically stable in the region Ω.

Proof.

For model (1.5), we construct a Lyapunov function of the form (2.6)V(N,P)=N*N1-e-aξ-d-δP*1-e-aξdξ+P*Pη-P*ηdη. Note that V(N,P) is non-negative, V(N,P)=0 if and only if (N,P)=(N*,P*). Furthermore, the time derivative of V along the solutions of (1.5) is (2.7)dVdt=1-e-aN-d-δP*1-e-aNdNdt+P-P*PdPdt. Substituting the expressions of dN/dt and dP/dt defined in (1.5) into (2.7), we can obtain (2.8)dVdt=(1-e-aN-d-δP*)(rN(1-N)1-e-aN-P*)-δ(P-P*)2  =(1-e-aN-d-δP*)(rN(1-N)1-e-aN-rN*(1-N*)1-e-aN*)-δ(P-P*)2. Define (2.9)ϕ(N)=1-e-aN-d-δP*, then (2.10)ϕ(N*)=0,ϕ(N)=1+ae-aN>0. So, (2.11)(N-N*)(ϕ(N)-ϕ(N*))>0. Then we can get (2.12)dVdt0. If (2.13)(1-e-aN-d-δP*)(rN(1-N)1-e-aN-rN*(1-N*)1-e-aN*)0 holds, which is equivalent to (2.14)(N-N*)(rN(1-N)1-e-aN-rN*(1-N*)1-e-aN*)0. Set φ(N)=rN(1-N)/(1-e-aN), we obtain φ(0)=0 and (2.15)φ(N)=r((1-2N)(1-e-aN)-aN(1-N)e-aN)(1-e-aN)2. And set ψ(N)=(1-2N)(1-e-aN)-aN(1-N)e-aN, we can get ψ(0)=0 and (2.16)ψ(N)=e-aN(a2N(1-N)+2)-2,ψ′′(N)=ae-aN(a2N2-a(a+2)N+a-2).

In view of a2, it follows that ψ′′(N)0 and ψ(N)ψ(0)=0 in the region Ω. Then ψ(N)0 is always true. It follows that φ(N)0, that is, V0. Consequently, the function V(N,P) satisfies the asymptotic stability theorem . Hence, E*=(N*,P*) is globally asymptotically stable. This completes the proof.

3. Applications

In this paper, we construct a new Lyapunov function for proving the global asymptotical stability of model (1.5). The new Lyapunov function is useful not only to model (1.5), but also to other models.

In this section, we will give some examples to show the robustness of the Lyapunov function (2.6). The parameters of the following models are positive and have the same ecological meanings with those of in model (1.5).

Example 3.1.

Considering the following Ivlev predator-prey model incorporating prey refuges (see ): (3.1)dNdt=rN(1-NK)-(1-e-a(1-m)N)P,dPdt=(1-e-a(1-m)N-δP-d)P, where m[0,1) is a refuge protecting of the prey. We can choose a Lyapunov functional as follows: (3.2)V(N,P)=N*N1-e-a(1-m)ξ-δP*-d1-e-a(1-m)ξdξ+P*Pη-P*ηdη.

The proof is similar to that of the Section 2.

Example 3.2.

Considering the following predator-prey model with Rosenzweig functional response (see ): (3.3)dNdt=rN(1-NK)-bNμP,dPdt=(cNμ-δN-d)P, where μ(0,1] is the victim’s competition constant. We can choose a Lyapunov functional as follows: (3.4)V(N,P)=N*NcNξ-δP*-dbNξdξ+P*Pη-P*ηdη.

We omit the proof here.

Example 3.3.

Considering the following model (1.1) with Holling-type functional response (see ): (3.5)dNdt=rN(1-NK)-bg(N)P,dPdt=(g(N)-δP-d)P, where g(N)=αN/(β+N) is known as a Holling type-II function, g(N)=αN2/(β+N2) as a Holling type-III function and g(N)=αN2/(β+ωN+N2) as a Holling type-IV function. We choose a Lyapunov function: (3.6)V(N,P)=N*Ng(ξ)-δP*-dbg(ξ)dξ+P*Pη-P*ηdη. For more details, we refer to .

Example 3.4.

Considering the following diffusive Ivlev-type predator-prey model (see ): (3.7)Ndt=rN(1-N)-(1-e-aN)P+d12N,Pdt=(ε(1-e-aN)-d)P+d22P, where the nonnegative constants d1 and d2 are the diffusion coefficients of N and P, respectively. 2=2/x2+2/y2, the usual Laplacian operator in two-dimensional space, is used to describe the Brownian random motion.

Model (3.7) is to be analyzed under the following non-zero initial conditions (3.8)N(x,y,0)0,P(x,y,0)0,(x,y)Ω=[0,Lx]×[0,Ly], and zero-flux boundary conditions: (3.9)Nn=Pn=0. In the above, n is the outward unit normal vector of the boundary Ω.

In order to give the proof of the global stability, we construct a Lyapunov function: (3.10)E(t)=ΩV(N(t),P(t))dxdy, where (3.11)V(N,P)=N*Nε(1-e-aξ)-d1-e-aξdξ+P*Pη-P*ηdη.

Then, differentiating E(t) with respect to time t along the solutions of model (3.7), we can obtain (3.12)dE(t)dt=ΩdVdtdxdy+Ω(VNd12N+VPd22P)dxdy. Using Green's first identity in the plane, and considering the zero-flux boundary conditions (3.9), one can show that (3.13)dE(t)dt=ΩdVdtdxdy-d12VN2Ω[(Nx)2+(Py)2]dxdy-d22VP2Ω[(Px)2+(Py)2]dxdy  ΩdVdtdxdy. The remaining arguments are rather similar as Theorem 2.3.

Alonso D. Bartumeus F. Catalan J. Mutual interference between predators can give rise to Turing spatial patterns Ecology 2002 83 1 28 34 10.1890/0012-9658(2002)083[0028:MIBPCG]2.0.CO;2 Abrams P. Ginzburg L. The nature of predation: prey dependent, ratio dependent or neither? Trends in Ecology and Evolution 2000 15 8 337 341 10.1016/S0169-5347(00)01908-X Hsu S. B. On global stability of a predator-prey system Mathematical Biosciences 1978 39 1-2 1 10 MR0472126 10.1016/0025-5564(78)90025-1 ZBL0383.92014 Ivlev V. S. Experimental Ecology of the Feeding of Fishes 1961 Yale University Press Kooij R. E. Zegeling A. A predator-prey model with Ivlev's functional response Journal of Mathematical Analysis and Applications 1996 198 2 473 489 10.1006/jmaa.1996.0093 MR1376275 ZBL0851.34030 Sugie J. Two-parameter bifurcation in a predator-prey system of Ivlev type Journal of Mathematical Analysis and Applications 1998 217 2 349 371 10.1006/jmaa.1997.5700 MR1492094 ZBL0894.34025 Feng J. Chen S. Global asymptotic behavior for the competing predators of the Ivlev types Mathematica Applicata 2000 13 4 85 88 MR1787751 ZBL1037.92029 Xiao H. B. Global analysis of Ivlev's type predator-prey dynamic systems Applied Mathematics and Mechanics 2007 28 4 461 470 10.1007/s10483-007-0406-1 MR2325161 ZBL1231.37054 Collings J. B. Bifurcation and stability analysis of a temperature-dependent mitepredator-prey interaction model incorporating a prey refuge Bulletin of Mathematical Biology 1995 57 1 63 76 Rosenzweig M. L. Paradox of enrichment: destabilization of exploitation ecosystems in ecological time Science 1971 171 3969 385 387 10.1126/science.171.3969.385 Holling C. S. The functional response of predators to prey density and its role in mimicry and population regulation Memoirs of the Entomological Society of Canada 1965 97 S45 5 60 10.4039/entm9745fv Ruan S. Xiao D. Global analysis in a predator-prey system with nonmonotonic functional response SIAM Journal on Applied Mathematics 2000/01 61 4 1445 1472 10.1137/S0036139999361896 MR1813689 Wang W. Zhang L. Wang H. Li Z. Pattern formation of a predator-prey system with Ivlev-type functional response Ecological Modelling 2010 221 2 131 140 10.1016/j.ecolmodel.2009.09.011 Merkin D. R. Afagh F. F. Smirnov A. L. Introduction to the Theory of Stability 1997 24 New York, NY, USA Springer Texts in Applied Mathematics MR1418401