Stability and Hopf bifurcation of a ratio-dependent predator prey model with time delay and stage structure

Transcription

1 Electronic Journal of Qualitative Theory of Differential Equations 6, No 99, 3; doi: 43/ejqtde699 Stability and Hopf bifurcation of a ratio-dependent predator prey model with time delay and stage structure Yan Song, Ziwei Li and Yue Du School of Mathematics and Physics, Bohai University, Jinzhou 3, P R China Received 9 March 6, appeared 6 October 6 Communicated by Hans-Otto Walther Abstract In this paper, a ratio-dependent predator prey model described by Holling type II functional response with time delay and stage structure for the prey is investigated By analyzing the corresponding characteristic equations, the local stability of the coexistence equilibrium of the model is discussed and the existence of Hopf bifurcations at the coexistence equilibrium is established By using the persistence theory on infinite dimensional systems, it is proven that the system is permanent if the coexistence equilibrium exists By introducing some new lemmas and the comparison theorem, sufficient conditions are obtained for the global stability of the coexistence equilibrium Numerical simulations are carried out to illustrate the main results Keywords: ratio-dependence, time delay, stage structure, Hopf bifurcation, stability Mathematics Subject Classification: 34C3, 34D3, 9D5 Introduction Predator prey models are important in the models of multi-species population interactions One of the important objectives in population dynamics is to comprehend the dynamical relationship between predator and prey, which had long been and will continue to be one of the dominant themes in both ecology and mathematical ecology It is well know that the functional response is a key factor in all predator-prey interactions, which describes the number of prey consumed by per predator per unit time Based on experiments, Holling [8] suggested three different kinds of functional responses, ie Holling type I, Holling type II and Holling type III, for different kinds of species to model the phenomena of predation, which made the standard Lotka Volterra system more realistic These functional responses are generally modeled as being a function of prey density only, ie the number of prey that an individual predator kills is a function of prey density only, and ignore the potential effects of predator density So they are usually called prey-dependent functional responses Obviously, this assumption can not explain the dynamics of the system completely when the variations Corresponding author jzsongyan@63com

2 Y Song, Z W Li and Y Du in predator size have an influence on the system Therefore a new theory, so-called predatordependent functional response, has been developed to consider the influence of both prey and predator populations There have been several famous predator-dependent functional response types: Hassel Varley type [3]; Beddington DeAngelis type [6, ]; Crowley Martin type []; and the well-known ratio-dependence type [] In the ratio-dependence theory, it roughly states that the per capita predator growth rate should be a function of the ratio of prey to predator abundance Moreover, as the number of predators often changes slowly (relative to prey number), there is often competition among the predators, and the per capita rate of predation should therefore depend on the numbers of both prey and predator, most probably and simply on their ratio These hypotheses are strongly supported by numerous field and laboratory experiment and observations (see, for example, [3 5, 7]) Let x(t) and y(t) be the densities of the prey and the predator at time t, respectively, a standard predator-prey model with Holling type functional response is of the form (see [8]) x (t) = xg(x) Φ(x)y, y (t) = eφ(x)y dy () In (), the function g(x) represents the growth rate of the prey in the absence of predation and d is the mortality rate of the predator in the absence of prey; the function Φ(x) is called functional response representing the prey consumption per unit time; e is the rate of conversion of nutrients from the prey into the reproduction of the predator But in ratio-dependent predator prey model, model () is described as ẋ(t) = xg(x) Φ ( x ẏ(t) = eφ y ( ) x y, y ) y dy In (), Φ( x y ) is the ratio-dependent predator functional response Many authors have studied predator-prey models with functional response, especially with ratio-dependent functional response Hsu et al [] investigated a predator-prey model with Hassell Varley type functional response It was shown that the predator free equilibrium is a global attractor only when the predator death rate is greater than its growth ability and the positive equilibrium exists if the above relation reverses In cases of practical interest, it was shown that the local stability of the positive steady state implies it global stability with respect to positive solutions For terrestrial predators that from a fixed number of tight groups, it was shown that the existence of an unstable positive equilibrium in the predator-prey model implies the existence of an unique nontrivial positive limit cycle Cantrell and Cosner [9] investigated predator-prey models with Beddington DeAngelis functional response (with or without diffusion) Criteria for permanence and for predator extinction were derived For systems without diffusion or with no-flux boundary conditions, criteria were derived for the existence of a globally stable coexistence equilibrium or, alternatively, for the existence of periodic orbits Kuang and Beretta [], Berezovskaya et al [8] investigated a ratio-dependent predator prey model with Michaelis Menten or Holling II type functional response, respectively In [], the authors proved that if the positive steady state of the system is locally asymptotically stable then the system has no nontrivial positive periodic solutions They also gave sufficient conditions for each of the possible three steady states to be globally asymptotically stable In [8], the authors gave a complete parametric analysis of stability properties and dynamic regimes of the model ()

3 Stability and Hopf bifurcation of a predator prey model 3 Beretta and Kuang [7], Xiao and Li [8] investigated a ratio-dependent predator-prey model with Michaelis Menten functional response and time delay, respectively In [7], the authors made use of a rather novel and non-trivial way of constructing proper Lyapunov functions to obtain some new and significant global stability or convergence results In [8], the authors studied the effect of time delay on local stability of the interior equilibrium and investigated conditions on the delay and parameters so that the interior equilibrium of the model is conditionally stable or unstable It was also shown that the interior equilibrium cannot be absolutely stable for all parameters Hsu et al [9] investigated a ratio-dependent one-prey two-predators model It was shown that the dynamites outcome of the interactions are very sensitive to parameter values and initial dates, which reveal far richer dynamics compared to similar prey dependent models However, it is assumed in these works that each individual prey admits the same risk to be attacked by predator This assumption is obviously unrealistic for many animals In natural world, the growth of species often has its development process, immature and mature, while in each stage of its development, it always shown different characteristic For instance, the mature species have preying capacity, while the immature species are raised by their parents and not able to prey Hence, stage-structured models may be more realistic Aiello and Freedman [] proposed and studied stage structured single-species population model with time delay Chen et al [] proposed and discussed a stage structured singlespecies population model without time delay Based on the ideas above, many authors have studied different kinds of biological models with stage structure Among these models, there are many factors that affect dynamical properties of predator-prey system such as the ratiodependent functional response, stage structure, and time delay, etc, especially the joint effect of these factors (see, for example, [3, 4, 4, 5, 7, 9, 3]) In order to analyze the effect of stage structure for prey on the dynamics of ratio-dependent predator-prey system, in [9], the authors proposed and studied the following differential system x (t) = ax (t) r x (t) bx (t), x (t) = bx (t) b x(t) a x (t)y(t) my(t) + x (t), ( ) a y x (t) (t) = y(t) r + my(t) + x (t) Sufficient conditions were derived for the uniform persistence and the global asymptotic stability of nonnegative equilibria of the model However, time delay is an important factor in biological models, since time delay could cause a stable equilibrium to become unstable and cause the species to fluctuate The main purpose of this paper is to study the effect of stage structure for the prey and time delay on the dynamics of a ratio-dependent predator-prey system described by Holling type II functional response To do so, we study the following differential system x (t) = rx (t) re d τ x (t τ) d x (t), x (t) = re dτ x (t τ) d x(t) ax (t)y(t) x (t) + my(t), y (t) = bx (t)y(t) x (t) + my(t) d 3y(t) (3)

4 4 Y Song, Z W Li and Y Du In (3), x (t) and x (t) represent the densities of the immature and the mature prey at time t, respectively; y(t) represents the density of the predator at time t; τ is the maturity of prey; r is the birth rate of the immature prey; d and d are the death rates of the immature and mature prey, respectively; re dτ x (t τ) represents the quantity which the immature born at time t τ can survive at time t; d 3 is the death rate of the predator; a is the capturing rate of the predator; b a is the conversion rate of nutrients into the reproduction of the predator; all the parameters are positive The initial conditions for system (3) take the form x (θ) = ϕ (θ), x (θ) = ϕ (θ), y(θ) = ϕ 3 (θ), θ [ τ, ], ϕ () >, ϕ () >, ϕ 3 () >, (4) where (ϕ (θ), ϕ (θ), ϕ 3 (θ)) C([ τ, ], R 3 + ), R3 + = {(x, x, x 3 ) x i, i =,, 3} In order to ensure the initial continuous, we suppose further that x () = τ re d s ϕ (s)ds The organization of this paper is as follows In the next section, we introduce some lemmas which will be essential to our proofs and discussions In Section 3, by analyzing the corresponding characteristic equations, the local stability of the coexistence equilibrium of system (3) is discussed Furthermore, the conditions for the existence of Hopf bifurcations at the coexistence equilibrium are obtained In Section 4, by using persistence theory on infinite dimensional systems, we prove that system (3) is permanent when the coexistence equilibrium exists In Section 5, by using comparison argument, the global stability of the coexistence equilibrium of system (3) is discussed In Section 6, numerical simulations are carried out to illustrate the main results A brief conclusion is given in Section 7 to conclude this work Preliminaries In this section, we introduce some lemmas which will be useful in next section By the fundamental theory of functional differential equations [5], it is well known that system (3) has a unique solution (x (t), x (t), y(t)) satisfying initial conditions (4) Further, it is easy to show that all solutions of system (3) with initial conditions (4) are defined on [, + ) and remain positive for all t Lemma All positive solutions of system (3) satisfying initial conditions (4) are ultimately bounded Proof We know that all solutions of system (3) are positive Hence we study only in the domain R 3 + = {(x, x, x 3 ) x i >, i =,, 3} is Let V(t) = bx (t) + bx (t) + ay(t), then he derivative of V(t) along solution of system (3) V(t) brx (t) bd x (t) bd x (t) ad 3 y(t) µv(t) + b(r + d )x (t) bd x (t) µv(t) + b(r + d ) 4d,

5 Stability and Hopf bifurcation of a predator prey model 5 where µ = min{d, d, d 3 } Therefore we derive that t V(t) e (V() µt + b(r + d ) 4d ) e µs ds = e µt V() + b(r + d ) ( e µt ) 4d µ b(r + d ) 4d µ The proof of Lemma is completed Lemma ([6]) Consider the following system (t + ) u (t) = au(t τ) bu(t) cu (t) here a, c, τ >, b, and u(t) > for t [ τ, ], we have (i) if a < b, then lim t + u(t) = ; (ii) if a > b, then lim t + u(t) = a b c Lemma 3 Consider the following system u (t) = ru (t) re d τ u (t τ) d u (t), u (t) = re d τ u (t τ) d u (t), () here r, d, d, τ > and u i (t) > (i =, ) for t [ τ, ], we have lim u (t) = r e dτ ( e dτ ), lim t + d d u (t) = re dτ t + d Proof It is easy to see that system () has two equilibria F (, ) and F (û, û ), where û = r e d τ ( e d τ ) d d, û = re d τ d, and easily show that F is unstable and F is locally asymptotically stable By the second equation of system () and Lemma, we derive that lim u (t) = re dτ = û t + d Therefore the limit equation of the first equation of system () takes the form which implies that u (t) = r e dτ ( e dτ ) d u (t), d lim u (t) = r e dτ ( e dτ ) = û, t + d d that is, the equilibrium F is globally asymptotically stable This proves Lemma 3 Lemma 4 ([8]) Consider the following system ( ) a u (t) = b + mu(t) d We have that lim t + u(t) = a bd md u(t), a, b, m, d > if a > bd and lim t + u(t) = if a < bd

6 6 Y Song, Z W Li and Y Du Lemma 5 Consider the following system u (t) = re dτ u(t τ) d u (t) apu(t) u(t) + mp with r, d, d, τ, m, a, P >, u(t) > for t [ τ, ], we have lim t + u(t) = u if mre d τ > a, where u = U + U + 4V, d U = re dτ d mp, V = d P(mre dτ a) Proof It is easy to know that u(t) > for all t > For any t >, we have re d τ u(t τ) d u (t) a m u(t) < u (t) < re d τ u(t τ) d u (t) By Lemma, we know that there exists a t > such that Then we get that u := mre d τ a md ε < u(t) < re dτ d + ε =: u for all t t re d τ u(t τ) d u (t) apu(t) u + mp < u (t) < re d τ u(t τ) d u (t) apu(t) u + mp By the comparison theorem and Lemma, there exists a t > t such that ap u +mp u := re dτ ε < u(t) < re dτ ap u +mp + ε =: u for all t t d d and < u < u < u(t) < u < u for all t t Continuing this process, we derive the sequence {u n } n= and {u n} n= with where < u < u < < u n < u n < < u < u, t > t n, u n := re d τ ap u n +mp u n +mp d ε, u n = re dτ ap d + ε By the bounded monotonic principle, we know that the limit of the sequence {u n } n= and {u n } n= exists Denote u = lim n u n and u = lim n u n, then we easily know that u = u and lim t + u(t) = u = u =: u, where u = U + U + 4V, d U = re dτ d mp, V = d P(mre dτ a)

7 Stability and Hopf bifurcation of a predator prey model 7 3 Local stability and Hopf bifurcation It is easy to show that system (3) always has a trivial equilibrium E (,, ) and a predatorextinction equilibrium E ( x, x, ), where x = r e d τ ( e d τ ) d d, x = re dτ d Further, if < b d 3 < mbre d τ a E (x, x, y ), where holds, then system (3) has a unique coexistence equilibrium x = r( e d τ ) d x, x = mbre dτ + ad 3 ab mbd, y = b d 3 x In this section, we are only concerned with the local stability of the coexistence equilibrium and the existence of Hopf bifurcation for system (3), since the biological meaning of the coexistence equilibrium implies that immature prey and mature prey and predator all exist For the coexistence equilibrium E (x, x, y ), the characteristic equation of (3) has the form where (λ + d )[λ + A λ + A + (B λ + B )e λτ ] =, (3) A = b(mbre d τ + ad 3 ab) + (b d 3 )(ab ad 3 + mbd 3 ) mb, A = d 3(b d 3 )(mbre d τ + ad 3 ab) + ad 3 (b d 3 ) mb, B = re d τ, B = rd 3(b d 3 )e d τ b Clearly, λ = d is a negative real root of Eq(3) Other two roots of (3) are given by the roots of equation When τ =, Eq(3) becomes By calculation, we know that λ + A λ + A + (B λ + B )e λτ = (3) λ + (A + B )λ + A + B = A + B = mbd 3(b d 3 ) + ad 3 + b (mr a) mb, A + B = d 3(b d 3 )(mbr ab + ad 3 ) mb > Hence, E is locally asymptotically stable if mbd 3 (b d 3 ) + ad 3 + b (mr a) > and unstable if mbd 3 (b d 3 ) + ad 3 + b (mr a) <

8 8 Y Song, Z W Li and Y Du If λ = iω(ω > ) is a purely imaginary root of Eq(3), separating real and imaginary parts, we have ω A = B ω sin(ωτ) + B cos(ωτ), A ω = B ω cos(ωτ) + B sin(ωτ) Eliminating sin(ωτ) and cos(ωτ), we obtain the equation with respect to ω ω 4 + (A B A )ω + A B = (33) Since A >, B <, A + B >, then A B > Therefore, if B + A A < A B, Eq (33) has no positive real roots Accordingly, by [, Theorem 34], we see that if mbd 3 (b d 3 ) + ad 3 + b (mr a) > and B + A A < A B hold, then E is locally asymptotically stable for all τ < d ln mbr a(b d 3 ) If B + A A > A B, Eq(33) has two positive real roots denoted by ω + = (B + A A ) +, ω = (B + A A ), respectively, where = (B + A A ) 4(A B ) Denote kπ + arccos (B A B )ω+ A B τ (k) B + = ω + +B, ω + kπ + arccos (B A B )ω A B τ (k) B = ω +B, ω k =,,, In the following we verify transversality condition of Eq (3) respect to τ, it follows that Differentiating (3) with [λ + A + B e λτ τ(b λ + B )e λτ ] dλ dτ λ(b λ + B )e λτ = By direct calculation, we derive that ( ) dλ = λ + A + B e λτ τ(b λ + B )e λτ dτ λ(b λ + B )e λτ λ + A = λ(λ + A λ + A ) + B λ(b λ + B ) τ λ, ( ) dλ { } ωi + A Re = Re dτ ωi( ω + A + A ωi) + B ωi(b ωi + B ) λ=ωi { } d Re λ sign dτ λ=ωi = ω + A B A B ω + B, { ( ) } dλ = sign Re = sign{ω + A dτ B A } λ=ωi

9 Stability and Hopf bifurcation of a predator prey model 9 Therefore { } d Re λ sign dτ { } d Re λ sign dτ λ=ω + i λ=ω i = sign { ω+ + A B A } { } = sign >, = sign { ω + A B A } { = sign } < Summarizing the above discussion, we have the following theorem on the local stability of E and Hopf bifurcations of system (3) Theorem 3 Assume that < b d 3 < mbre d τ a For system (3), we have the following (i) If mbd 3 (b d 3 ) + ad 3 + b (mr a) > and B + A A < A B, then the coexistence equilibrium E is locally asymptotically stable for all τ < d ln mbr a(b d 3 ) (ii) If mbd 3 (b d 3 ) + ad 3 + b (mr a) > and B + A A > A B, then there exists a τ = τ () +, such that E is stable for τ < τ () + (3) undergoes a Hopf bifurcation at E when τ = τ () + and unstable for τ > τ() + Furthermore, system (iii) If mbd 3 (b d 3 ) + ad 3 + b (mr a) < and B + A A < A B, then the coexistence equilibrium E is unstable for all τ < d ln mbr a(b d 3 ) (iv) If mbd 3 (b d 3 ) + ad 3 + b (mr a) < and B + A A > A B, then there exists a τ = τ (), such that E is unstable for τ < τ () and stable for τ > τ() Furthermore, system (3) undergoes a Hopf bifurcation at E when τ = τ () 4 Permanence In this section, we are concerned with the permanence of system (3) Definition 4 System (3) is said to be permanent (uniformly persistent) if there are positive constants m and M such that each positive solution of system (3) (x (t), x (t), y(t)) satisfies m lim t + inf x i(t) lim t + sup x i(t) M, i =,, m lim inf y(t) lim sup y(t) M t + t + In order to study the permanence of system (3), we present the persistence theory on infinite dimensional systems from [6] Let X be a complete metric space with metric d The distance d(x, Y) of a point x X from a subset Y of X is defined by d(x, Y) = inf d(x, y) y Y Assume that X X, X X and X X = φ Also, assume that T(t) is a C semigroup on X satisfying T(t) : X X, T(t) : X X (4) Let T b (t) = T(t) X and A b be the global attractor for T b (t)

10 Y Song, Z W Li and Y Du Lemma 4 Suppose that T(t) satisfies (4) and the following conditions: (i) there is a t such that T(t) is compact for t > t ; (ii) T(t) is point dissipative in X; (iii) Ã b = x Ab ω(x) is isolated and has an acyclic covering M, where M = {M, M, M n }; (iv) W s (M i ) X = φ for i =,,, n Then X is a uniform repeller with respect to X, that is, there is an ε > such that for any x X, lim t + inf d(t(t)x, X ) ε Theorem 43 If < b d 3 < mbre d τ a holds, then system (3) is uniformly persistent Proof We need only to prove that the boundaries of R 3 + repel positive solutions of system (3) uniformly Let C + ([ τ, ], R 3 + ) denote the space of continuous functions mapping [ τ, ] into R 3 + Define C = { (ϕ, ϕ, ϕ 3 ) C + ([ τ, ], R 3 +) ϕ (θ), ϕ (θ), θ [ τ, ] }, C = {(ϕ, ϕ, ϕ 3 ) C + ([ τ, ], R 3 +) ϕ (θ) >, ϕ (θ) >, ϕ 3 (θ), θ [ τ, ]}, C = C C, X = C + ([ τ, ], R 3 +), C = int C + ([ τ, ], R 3 +) In the following, we verify that the conditions in Lemma 4 are satisfied By the definition of C and C, it is easy to see that C and C are positively invariant Moreover, the conditions (i) and (ii) in Lemma 4 are clearly satisfied (see for instance [, Theorem 8]) Thus we need only to show that the conditions (iii) and (iv) hold Clearly, corresponding to x (t) = x (t) = y(t) = and x (t) = x, x (t) = x, y(t) =, respectively, there are two constant solutions in C : Ẽ C, Ẽ C satisfying Ẽ = {(ϕ, ϕ, ϕ 3 ) ([ τ, ], R 3 +) ϕ (θ), ϕ (θ), ϕ 3 (θ), θ [ τ, ]}, Ẽ = {(ϕ, ϕ, ϕ 3 ) ([ τ, ], R 3 +) ϕ (θ) x, ϕ (θ) x, ϕ 3 (θ), θ [ τ, ]} We now verify the condition (iii) of Lemma 4 If (x (t), x (t), y(t)) is a solution of system (3) initiating from C, then ẏ (t) = d 3 y(t), which yields y(t) as t + If (x (t), x (t), y(t)) is a solution of system (3) initiating from C with x () >, x () >, then it follows from the first and the second equations of system (3) that By Lemma 3, we get that x (t) = rx (t) re d τ x (t τ) d x (t), x (t) = re d τ x (t τ) d x (t) lim x (t) = r e dτ ( e dτ ) = x, lim t + d d x (t) = re dτ = x t + d that is, (x (t), x (t), y(t)) ( x, x, ) as t + Noting that C C = φ, it follows that the invariant sets Ẽ and Ẽ are isolated Hence, {Ẽ, Ẽ } is isolated and is an acyclic covering satisfying the condition (iii) in Lemma 4 We now verify the condition (iv) of Lemma 4 Here we only show that W s (Ẽ ) C = φ holds since the proof of W s (Ẽ ) C = φ is simple Assume that W s (Ẽ ) C = φ

11 Stability and Hopf bifurcation of a predator prey model Then there is a positive solution of system (3) (x (t), x (t), y (t)) initiating from C with lim t + (x (t), x (t), y (t)) = E ( x, x, ) Therefore we have lim t + x (t) = x, that is, for ε > small enough, there exists a t > such that x ε < x (t) < x + ε for all t > t + τ It follows from the third equation of system (3) that for t > t + τ y (t) [ ] b( x ε) x ε + my (t) d 3 y (t) By Lemma 4 and comparison theorem, we get that lim t + y (t) (b d 3)( x ε) Since ε > is arbitrary small, then we conclude lim t + y (t) (b d 3) x, which contradicts y (t) (t + ) Hence, we have W s (Ẽ ) C = φ By Lemma 4, we are now able to conclude that C repel positive solutions of system (3) uniformly Therefore system (3) is permanent The proof is complete 5 Global stability In this section, we are concerned with the global stability of the coexistence equilibrium of system (3) Theorem 5 The coexistence equilibrium E of system (3) is globally asymptotically stable provided that (i) < b d 3 < mrd 3e d τ a ; (ii) mre d τ > a Proof Let (x (t), x (t), y(t)) be any positive solution of system (3) with initial conditions (4) We derive from the second equation of system (3) that x (t) re d τ x (t τ) d x (t) By comparison theorem and Lemma, we have lim t + x (t) re d τ d Therefore, for any ε >, there exists a T > such that x (t) < re d τ + ε =: N for t > T d We derive from the third equation of system (3) that [ y (t) bn N + my(t) d 3 ] y(t) Since bn N d 3 = N (b d 3 ) >, we get from Lemma 4 and comparison theorem that lim t + y(t) (b d 3)N

12 Y Song, Z W Li and Y Du Then there exists a T > T such that y(t) < (b d 3)N + ε =: P for t > T We derive from the second equation of system (3) that x (t) re dτ x (t τ) d x(t) ap x (t) x (t) + mp Since mre dτ > a, we get from Lemma 5 and comparison theorem that there exists a T 3 > T such that x (t) > z ε =: N for t > T 3, where and z z = U + U + 4V, d U = re dτ d mp, V = d P (mre dτ a) is the positive root for the equation re d τ d x ap x + mp = We derive from the third equation of system (3) that [ y (t) > bn N + my(t) d 3 ] y(t) From Lemma 4 and comparison theorem we get that there exists a T 4 > T 3 such that y(t) > (b d 3)N ε =: P for t > T 4 (5) We derive from the first equation of system (3) that rn re d τ N d x (t) < x (t) < rn re d τ N d x (t), t T 4 Then there exists a T 5 > T 4 such that M := r(n e d τ N ) d ε < x (t) < r(n e dτ N ) d + ε =: M for t > T 5 Hence we have that M < x (t) < M, N < x (t) < N, P < y(t) < P, t > T 5 Replacing (5) into the second equation of (3), we have x (t) < re dτ x (t τ) d x(t) ap x (t) x (t) + mp By Lemma 5 and comparison theorem we get that there exists a T 6 > T 5 such that

13 Stability and Hopf bifurcation of a predator prey model 3 where x (t) < z + ε =: N for t > T 6, (5) z = U + U + 4V, d U = re dτ d mp, V = d P (mre dτ a) Replacing (5) into the third equation of (3), we have [ y (t) bn N + my(t) d 3 ] y(t) By Lemma 4 and comparison theorem we get that there exists a T 7 > T 6 such that y(t) < (b d 3)N + ε =: P for t > T 7 (53) Replacing (53) into the second equation of (3), we have x (t) re dτ x (t τ) d x(t) ap x (t) x (t) + mp By Lemma 5 and comparison theorem we get that there exists a T 8 > T 7 such that where x (t) > z 3 ε =: N for t > T 8, (54) z3 = U 3 + U3 + 4V 3, d U 3 = re dτ d mp, V 3 = d P (mre dτ a) Replacing (54) into the third equation of (3), we have [ y (t) > bn N + my(t) d 3 ] y(t) By Lemma 4 and comparison theorem we get that there exists a T 9 > T 8 such that y(t) > (b d 3)N ε =: P for t > T 9 Replacing (5) and (54) into the first equation of (3), we have rn re d τ N d x (t) < x (t) < rn re d τ N d x (t), t T 9 Then there exists a T > T 9 such that M := r(n e d τ N ) d ε < x (t) < r(n e dτ N ) d + ε =: M for t > T Hence we have that < M < M < x (t) < M < M, < N < N < x (t) < N < N, < P < P < y(t) < P < P, t > T

14 4 Y Song, Z W Li and Y Du Continuing this process, we derive the six sequences {M n } n=, {M n} n=,{n n} n=, {N n} n=, {P n } n=, {P n} n= with < M < M < < M n < x (t) < M n < < M < M, < N < N < < N n < x (t) < N n < < N < N, < P < P < < P n < y(t) < P n < < P < P, t > T 4n+ Since M, N, P, M, N, P >, then the sequences {M n } n=, {N n} n=, {P n} n= are bounded and decrease, {M n } n=, {N n} n=, {P n} n= are bounded and increase, so there exist constants M, M, N, N, P, P such that lim t + M n = M, lim t + M n = M, lim t + N n = N, lim t + N n = N, lim t + P n = P, lim t + P n = P Easily know that M M, N N, P P Next we prove that M = M, N = N, P = P From above discussion, we have then P n = (b d 3)N n + ε, P m = (b d 3)N m ε, P n P m = (b d 3)(N n N m ) + ε But N n N n = re dτ d mp n + (re d τ d mp n ) + 4d P n (mre d τ a) d re dτ d mp n + (re d τ d mp n ) + 4d P n (mre d τ a) + ε d = re dτ d mp n + (re d τ + d mp n ) 4ad P n d re dτ d mp n + (re d τ + d mp n ) 4ad P n + ε d = { d m(p n P d n ) + (re d τ + d mp n ) 4ad P n } (re d τ + d mp n ) 4ad P n + ε = d = P n P n < P n P n = { d m(p n P n ) + d m(p n P n )[re d τ +d m(p n P n )+4ad (P n P n )] (re d τ +d mp n ) 4ad P n + (re d τ +d mp n ) 4ad P n } + ε { } mre m d τ +d m (P n +P n ) 4a + ε (re d τ +d mp n ) 4ad P n + (re d τ +d mp n ) 4ad P n { } m (mre dτ a) + d m (P n + P n ) + ε re d τ + d mp n + re d τ + d mp n a(p n P n ) re d τ + d m(p n + P n ) + ε < a re d τ (P n P n ) + ε = a re d τ (b d 3 )(N n N n ) ( + ε + a ) re d τ

15 Stability and Hopf bifurcation of a predator prey model 5 Taking n, we get that That is, N N a(b d ( 3) mrd 3 e d τ (N N) + ε + a ) re d τ ( (N N) a(b d ) ( 3) mrd 3 e d ε + a ) τ re d τ Since a(b d 3) mrd 3 e = mrd 3e dτ a(b d 3 ) > d τ mrd 3 e d, then N = N Therefore P = P τ With M n = r(n n e dτ N n ) + ε, M d n = r(n n e dτ N n ) ε, d we have Taking n, we get that M = M Hence M n M n = r( + e d τ )(N n N n ) d + ε lim x (t) = M = M, t + lim x (t) = N = N, t + lim y(t) = P = P t + It is easy to know that M = M = x, N = N = x, P = P = y Therefore lim t + x (t) = x, lim t + x (t) = x, lim t + y(t) = y, that is, E (x, x, y ) is globally attractive If mre dτ > a, then mbd 3 (b d 3 ) + ad 3 + b (mr a) > and E is locally asymptotically stable Therefore we conclude that E is globally asymptotically stable The proof is complete 6 Numerical simulations Now we give some numerical simulations to illustrate the main results Example 6 In (3), we let r = 7, a = 69, b = 35, d = 7, d =, d 3 =, m = 7 It is easy to verify that the conditions < b d 3 < mbre d τ a, mbd 3 (b d 3 ) + ad 3 + b (mr a) > and B + A A > A B hold According to Theorem 3 (ii), there exists a τ 586 such that the coexistence equilibrium E of system (3) is locally asymptotically stable if τ < τ and unstable if τ > τ Furthermore, system (3) undergoes a Hopf bifurcation at E if τ = τ Numerical simulations illustrate these results (see Figure 6 and Figure 6) Example 6 In (3), we let r = 7, a = 69, b = 35, d = 7, d =, d 3 =, m = 4 It is easy to verify that the conditions < b d 3 < mbre d τ a, mbd 3 (b d 3 ) + ad 3 + b (mr a) < and B + A A > A B hold According to Theorem 3 (iv), there exists a τ 373 such that the coexistence equilibrium E of system (3) is locally asymptotically stable if τ > τ and unstable if τ < τ Furthermore, system (3) undergoes a Hopf bifurcation at E if τ = τ Numerical simulations illustrate these results (see Figure 63 and Figure 64) Example 63 In (3), we let r = 8, a =, b =, d = 4, d = 8, d 3 =, m = 5, τ =, then < b d 3 < mbre d τ a and mre dτ > a hold By Theorem 43 we see that system (3) is uniformly persistent By Theorem 5 we see that the coexistence equilibrium E (434, 6584, 658) of system (3) is globally asymptotically stable Numerical simulations illustrate these results (see Figure 65)

16 6 Y Song, Z W Li and Y Du 3 x 5 x y values of x (t)x (t)andy(t) time t y(t) x (t) x (t) Figure 6: Numerical simulation shows that E is locally asymptotically stable for τ = < τ

17 Stability and Hopf bifurcation of a predator prey model x x y values of x (t)x (t)andy(t) time t Figure 6: Numerical simulation shows that E is unstable for τ = 598 > τ,which yields a Hopf bifurcation

18 8 Y Song, Z W Li and Y Du 3 x 5 x y values of x (t)x (t)andy(t) time t y(t) 5 x (t) x (t) Figure 63: Numerical simulation shows that E is locally asymptotically stable for τ = > τ

19 Stability and Hopf bifurcation of a predator prey model 9 5 x x y values of x (t)x (t)andy(t) time t Figure 64: Numerical simulation shows that E is unstable for τ = < τ, which yields a Hopf bifurcation

20 Y Song, Z W Li and Y Du 7 x 6 x y values of x (t)x (t)andy(t) time t 3 5 y(t) x (t) x (t) 3 4 Figure 65: Numerical simulation shows that E is globally asymptotically stable

21 Stability and Hopf bifurcation of a predator prey model 7 Conclusions It is well-known that many species go through two or more life stages as they proceed from birth to death Delay is common in population dynamics Any biological or environmental parameters are naturally subject to fluctuation in time Researches show that a system with time delays exhibits more complicated dynamics than that without time delay since time delay could bring a switch in the stability of equilibria and induce various oscillations and periodic solution Gourley and Kuang [4] investigated a general predator-prey model with stage structure for the predator and constant maturation time delay It was shown that if the juvenile death rate is nonzero, then for small and large values of maturation time delay, the population dynamics takes the simple form of a globally attractive steady state It was also shown that if the functional response function takes the Holling I type and the resource is dynamics, as in nature, there is a window in maturation time delay parameter that generates sustainable oscillatory dynamics In this paper, we have investigated a ratio-dependent predator-prey model described by Holling type II functional response with time delay and stage structure for the prey By analyzing the corresponding characteristic equations, the sufficient conditions for the local stability of the coexistence equilibrium and the existence of Hopf bifurcations are obtained By means of the persistence theory on infinite dimensional systems, it is proven that the system is permanent if the coexistence equilibrium is feasible By introducing some new lemmas and the comparison theorem, sufficient conditions are obtained for the global stability of the coexistence equilibrium We have shown the effect of stage structure and time delay on the dynamics of a ratio-dependent predator-prey system In system (3), the delay τ is the time taken from birth to maturity, d is the death rate of the immature prey, thus e dτ is the surviving rate of each immature prey before reaching maturity By Theorem 43, we see that system (3) is uniformly persistent if the birth rate into the immature prey population, the rate of immature prey becoming mature prey, and the conversion rate and the half saturation rate of the predator are high and the capturing rate of the predator and the death rates of both the immature prey and the predator are low enough satisfying the condition < b d 3 < mbre d τ a By Theorem 5, we see that the coexistence equilibrium is globally asymptotically stable under somewhat stronger assumptions than those in Theorem 43 on the uniformly persistent of system (3) By the discussion of Theorem 3(ii) and (iv), we can see that under some conditions the equilibrium E changes its stability and a periodic solution through Hopf bifurcation occurs when the delay τ passes through a critical value This implies that the time delay is able to cause a periodic evolution of the prey and predator populations and alter the dynamics of system (3) significantly Acknowledgements This work was supported by the Natural Science Foundation of China (667373) and the Natural Science Foundation of Liaoning Province of China (LN46) References [] W G Aiello, H I Freedman, A time delay model of single species growth with stage structure, Math Biosci (99), url

22 Y Song, Z W Li and Y Du [] R Arditi, L R Ginzburg, Coupling in predator prey dynamics: ratio-dependence, J Theoret Biol 39(989), 3 36 url [3] R Arditi, L R Ginzburg, H R Akcakaya, Variation in plankton densities among lakes: a case for ratio-dependent models, American Nat 38(99), url [4] R Arditi, N Perrin, H Saiah, Functional response and heterogeneities: an experiment test with cladocerans, Oikos 6(99), url [5] R Arditi, H Saiah, Empirical evidence of the role of heterogeneity in ratio-dependent consumption, Ecology 73(99), url [6] J R Beddington, Mutual interference between parasites or predators and its effect on searching efficiency, J Animal Ecol 44(975), url [7] E Beretta, Y Kuang, Global analysis in some delayed ratio-dependent predator prey systems, Nonlinear Anal 3(998), url [8] F Berezovskaya, G Karev, R Arditi, Parametric analysis of the ratio-dependent predator prey model, J Math Biol 43(), 46 url [9] R S Cantrell, C Cosner, On the dynamics of predator prey models with the Beddington DeAngelis functional response, J Math Anal Appl 57(), 6 url [] L S Chen, X Y Song, Z Lu, Mathematical models and methods in ecology, Sichuan Science and Technology Press, Sichuan, 3 [] P H Crowley, E K Martin, Functional response and interference within and between year classes of a dragonfly population, J N Am Benthol Soc 8(989), url [] D L DeAngelis, R A Goldstein, R V ONeill, A model for trophic interaction, Ecology 56(975), url [3] L W Deng, X D Wang, M Peng, Hopf bifurcation analysis for a ratio-dependent predator prey system with two delays and stage structure for the predator, Appl Math Comput 3(4), 4 3 url [4] S A Gourley, Y Kuang, A stage structured predator prey model and its dependence on maturation delay and death rate, J Math Biol 49(4), 88 url [5] J K Hale, Theory of functional differential equations, Springer, New York, 976 url [6] J Hale, P Waltman, Persistence in infinite-dimensional systems, SIAM J Math Anal (989), url [7] I Hanski, The functional response of predator: worries about scale, Trends Ecol Evol 6(99), 4 4 url [8] C S Holling, The functional response of predator to prey density and its role in mimicry and population regulation, Mem Ent Soc Can 97(965), 6 url [9] S B Hsu, T W Hwang, Y Kuang, Rich dynamics of a ratio-dependent one prey two predator model, J Math Biol 43(), url

23 Stability and Hopf bifurcation of a predator prey model 3 [] S B Hsu, T W Hwang, Y Kuang, Global dynamics of a predator prey model with Hassell Varley type functional response, Discrete Contin Dyn Syst Ser B (8), url [] Y Kuang, Delay differential equations with applications in population dynamics, Academic Press, New York, 993 url [] Y Kuang, E Beretta, Global qualitative analysis of a ratio-dependent predator prey system, J Math Biol 36(998), url [3] M P Massell, C C Varley, New inductive population model for insect parasites and its bearing on biological control, Nature 3(969), url [4] A K Misra, B Dubey, A ratio-dependent predator prey model with delay and harvesting, J Biol Syst 8(), url [5] P J Pal, P K Mandal, K K Lahiri, A delayed ratio-dependent predator-prey model of interacting populations with Holling type III functional response, Nonlinear Dyn 76(4), url [6] X Y Song, L S Chen, Optional harvesting and stability for a two species competitive system with stage structure, Math Biosci 7(), url [7] W Y Wang, L J Pei, Stability and Hopf bifurcation of a delayed ratio-dependent predator prey system, Acta Mech Sin 7(), url [8] D Xiao, W Li, Stability and bifurcation in a delayed ratio-dependent predator prey system, Proc Edin Math Soc 45(), 5 url [9] R Xu, M A J Chaplain, F A Davidson, Persistence and global stability of a ratiodependent predator prey model with stage structure, Appl Math Comput 55(4), url [3] R Xu, Q Gan, Z E Ma, Stability and bifurcation analysis on a ratio-dependent predatorprey model with time delay, J Comput Appl Math 3(9), 87 3 url