Dynamical analysis of a delayed predator-prey system with modified Leslie-Gower and Beddington-DeAngelis functional response

Transcription

1 Liu Advances in Difference Equations 014, 014:314 R E S E A R C H Open Access Dynamical analysis of a delayed predator-prey system with modified Leslie-Gower and Beddington-DeAngelis functional response Juan Liu * * Correspondence: liujuan716@163.com Department of Mathematics and Physics, Bengbu College, Bengbu, 33030, P.R. China Abstract The dynamics of a delayed predator-prey system with modified Leslie-Gower and Beddington-DeAngelis functional response is investigated. The main results are given in terms of local stability and local Hopf bifurcation. By regarding the possible combination of the feedback delays of the prey and the predator as a bifurcation parameter, sufficient conditions for the local stability and existence of the local Hopf bifurcation of the system are obtained. In particular, the properties of the local Hopf bifurcation such as direction and stability are determined by using the normal form method and center manifold theorem. Finally, numerical simulations are carried out to illustrate the main theoretical results. Keywords: delays; Hopf bifurcation; predator-prey system; stability; periodic solution 1 Introduction As we all know, one of the dominant themes in mathematical ecology is the dynamic relationship between predators and their prey. One of the important factors which affect the dynamics of biological and mathematical models is the functional response. Especially the Beddington-DeAngelis functional response which is first proposed by Beddington and DeAngelis et al. [1, ]. Beddington-DeAngelis functional response has desirable qualitative features of ratio dependent form but takes care of their controversial behaviors at low densities and it has an extra term in the denominator modeling mutual interference among predators. Predator-prey systems with a Beddington-DeAngelis functional response have beenstudiedbymanyauthors[3 1]. Ko and Ryu studied a diffusive two-competing prey and one-predator system with Beddington-DeAngelis functional response and discussed the stability and uniqueness of coexistence states [3]. Liu and Wang studied the global asymptotic stability of two stage-structured predator-prey systems with Beddington- DeAngelis functional response [6]. Li and Takeuchi studied the permanence, local stability, and global stability of two models of a density dependent predator-prey system with Beddington-DeAngelis functional response by using stability theory and Lyapunov functions [10]. Yu investigated the permanence, local stability, and global stability of the following modified Leslie-Gower predator-prey system with Beddington-DeAngelis functional 014 Liu; licensee Springer. This is an Open Access article distributed under the terms of the Creative Commons Attribution License which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

2 LiuAdvances in Difference Equations 014, 014:314 Page of 19 response [1]: { dxt) dt dyt) dt =r 1 pxt) αyt) 1+bxt)+cyt) )xt), =r βyt) xt)+k )yt), 1) xt) andyt) represent the densities of the prey and the predator at time t,respectively. r 1 and r are the intrinsic growth rates of the prey and the predator, respectively. p is the intra-specific competition coefficient of the prey. α is the maximum value of the per capita reduction rate of the prey due to the predator. β is a measure of the food quantity that the prey provides converted to the predator birth. It is well known that dynamical systems with one delay or multiple delays have been studied by many authors [7, 10, 11, 13 0]. Li and Wang studied the stability and Hopf bifurcation of a delayed three-level food chain model with Beddington-DeAngelis functional response [7]. Bianca et al. studied the stability and Hopf bifurcation of a mathematical framework that consists of a system of a logistic equation with a delay and an equation for the carrying capacity by regarding the delay as a bifurcation parameter [13]. Cui and Yan studied a three-species food chain system with two delays and they analyzed the Hopf bifurcation of the system by taking the sum of the two delays as the bifurcation parameter [17]. In [0], Bianca et al. further studied an economic growth model with two delays and they investigated the existence and properties of Hopf bifurcation of the model by regarding the possible combination of the two delays as the bifurcation parameter. To the best of our knowledge, seldom did authors consider the Hopf bifurcation of the delayed predatorprey system with modified Leslie-Gower and Beddington-DeAngelis functional response. Stimulated by this, in this paper we investigate the Hopf bifurcation of the following delayed predator-prey system with modified Leslie-Gower and Beddington-DeAngelis functional response: { dxt) αyt) 1+bxt)+cyt) )xt), =r dt 1 pxt τ 1 ) dyt) =r dt βyt τ ) xt τ )+k )yt), ) τ 1 0isthefeedbackdelayofthepreyandτ 0isthefeedbackdelayofthe predator. The main purpose of this paper is to consider the effect of the two delays on the dynamics of system ). This paper is organized as follows. In Section, we discuss local stability of the positive equilibrium and existence of the local Hopf bifurcation of system ). In Section 3, direction of the Hopf bifurcation and stability of the bifurcated periodic solutions are determined by using the normal form method and center manifold theorem. In Section 4, some numerical simulations are performed to illustrate our theoretical analysis and a final conclusion is given in Section 5. Local stability and the existence of Hopf bifurcation It is easy to verify that if the condition H 1 ): αr k < r 1 β + r 1 r ck holds, then system )has a unique positive equilibrium E x, y ), y = r x +k) β,andx = B+ B 4AC, A A = bpβ + cpr, B = αr + pβ + cpr k r 1 r c br 1 β, C = αr k r 1 β r 1 r ck. Let xt)=xt) x, ȳt)=yt)+y. Dropping the bars, system ) gets the following form: { dxt) dt = a 11 xt)+a 1 yt)+b 11 xt τ 1 )+f 1, dyt) dt = c 1 xt τ )+c yt τ )+f, 3)

3 LiuAdvances in Difference Equations 014, 014:314 Page 3 of 19 a 11 = αbx y 1 + bx + cy ), a 1 = αx 1 + bx ) 1 + bx + cy ), 4) b 11 = px, c 1 = δy x + k), c = βy x + k, 5) f 1 = g 1 x t)+g xt)yt)+g 3 y t)+g 4 xt)xt τ 1 )+g 5 x t)yt) 6) + g 6 xt)y t)+g 7 x 3 t)+g 8 y 3 t)+, 7) f = h 1 x t τ )+h xt τ )yt)+h 3 xt τ )yt τ ) 8) + h 4 x t τ )yt)+h 5 x t τ )yt τ )+h 6 x 3 t τ )+, 9) with g 1 = bαy 1 + cy ) 1 + bx + cy ) 3, g = α1 + bx + cy )+bcαx y 1 + bx + cy ) 3, 10) g 3 = cαx 1 + bx ) 1 + bx + cy ) 3, g 4 = p, g 5 = bα1 + bx )+bcαy bx cy ) 1 + bx + cy ) 4, 11) g 6 = cα1 + cy )+bcαx cy bx ) 1 + bx + cy ) 4, g 7 = b αy 1 + cy ) 1 + bx + cy ) 4, 1) g 8 = c αx 1 + bx ) 1 + bx + cy ), h 4 1 = βy x + k), h βy 3 = x + k), 13) h 3 = βy x + k), h 4 = βy x + k), h 3 5 = βy x + k), h βy 3 6 = x + k). 14) 4 The linearized system of system 3)is { dx dt = a 11xt)+a 1 yt)+b 11 xt τ 1 ), dy dt = c 1xt τ )+c yt τ ). 15) The characteristic equation of system 15)is λ + A 1 λ + B 1 λe λτ 1 +C 1 λ + C 0 )e λτ + D 0 e λτ 1+τ ) =0, 16) A 1 = a 11, B 1 = b 11, C 0 = a 11 c a 1 c 1, C 1 = c, D 0 = b 11 c. Case 1. τ 1 = τ =0. When τ 1 = τ =0,Eq.16)becomes λ + A 11 λ + A 10 =0, 17) A 10 = C 0 + D 0, A 11 = A 1 + B 1 + C 1.

4 LiuAdvances in Difference Equations 014, 014:314 Page 4 of 19 It is easy to obtain A 11 = r 1 + r + αy 1+cy ) 1+bx +cy ) >0,A 10 =a 11 + b 11 )c a 1 c 1.Fromthe expressions of a 11 and b 11,wegeta 11 + b 11 = r 1 αy 1+cy ) 1+bx +cy ) <0.ThenwecangetA 10 >0. Therefore, if the condition H 1 ) holds, then the positive equilibrium of system )without delay is locally asymptotically stable. Case. τ 1 >0,τ =0. When τ 1 >0,τ =0,Eq.16)becomes λ + A 1 λ + A 0 +B 1 λ + B 0 )e λτ 1 =0, 18) A 1 = A 1 + C 1, A 0 = C 0, B 1 = B 1, B 0 = D 0. Let λ = i 1 1 > 0) be the root of Eq. 18). Then we have { B1 1 cos τ 1 1 B 0 1 sin τ 1 1 = A 1 1, B 1 1 sin τ B 0 1 cos τ 1 1 = 1 A 0, from which one can obtain A 1 B 1 A 0) 1 + A 0 B 0 =0. 19) Since C 0 + D 0 > 0, we can conclude that if we have the condition H ): C 0 D 0 <0,then A 0 < B 0, and further Eq. 19) has a unique positive root A 1 10 = B 1 A 0)+ A 1 B 1 A 0) 4A 0 B 0 ), and the corresponding critical value of the delay is Define τ 1k = 1 arccos B 0 A 1 B 1 )10 A 0B 0 10 B kπ, k =0,1,,... B 0 10 τ 10 = 1 arccos B 0 A 1 B 1 )10 A 0B 0 10 B B 0 When τ 1 = τ 10,thenEq.18) has a pair of purely imaginary roots ±i 10. Differentiating the two sides of Eq. 18), one can obtain Thus, [ ] dλ 1 λ + A 1 = dτ 1 λλ + A 1 λ + A 0 ) + B 1 λb 1 λ + B 0 ) τ 1 λ. [ dλ Re dτ 1 ] 1 τ 1 =τ 10 = A 1 B 1 A 0) 4A 0 B 0 ) B >0. B 0

5 LiuAdvances in Difference Equations 014, 014:314 Page 5 of 19 Based on the analysis above and according to the Hopf bifurcation theorem in [1], we have the following results. Theorem 1 Suppose that the conditions H 1 ) and H ) hold. The positive equilibrium E x, y ) of system ) is asymptotically stable for τ 1 [0, τ 10 ). System ) undergoes a Hopf bifurcation at the positive equilibrium E x, y ) of system ) when τ 1 = τ 10 and a family of periodic solutions bifurcate from the positive equilibrium E x, y ) of system ). Case 3. τ 1 =0,τ >0. Substitute τ 1 =0intoEq.16), then Eq. 16)becomes λ + A 31 λ +B 31 λ + B 30 )e λτ =0, 0) A 31 = A 1 + B 1, B 31 = C 1, B 30 = C 0 + D 0, { B31 sin τ + B 30 cos τ =, B 31 cos τ B 30 sin τ = A 31, from which it follows that 4 + A 31 31) B 4 B 30 =0. 1) Obviously, B 30 >0.Thus,wecanconcludethatEq.1) has a unique positive root A 31 0 = B 31 A )+ 31 B 31 ) +4B 30, and the corresponding critical value of the delay is Define τ k = 1 arccos B 30 A 31 B 31 )0 0 B kπ, k =0,1,,... B 30 0 τ 0 = 1 arccos B 30 A 31 B 31 )0 0 B B 30 When τ = τ 0,thenEq.0) has a pair of purely imaginary roots ±i 0 and similar to Case, we can obtain [ ] dλ 1 A 31 B 31 ) +4B 30 Re = dτ τ =τ 0 B >0. B 30 Based on the analysis above and according to the Hopf bifurcation theorem in [1], we have the following results.

6 LiuAdvances in Difference Equations 014, 014:314 Page 6 of 19 Theorem Suppose that the condition H 1 ) holds. The positive equilibrium E x, y ) of system ) is asymptotically stable for τ [0, τ 0 ). System ) undergoes a Hopf bifurcation at the positive equilibrium E x, y ) of system ) when τ = τ 0 and a family of periodic solutions bifurcate from the positive equilibrium E x, y ) of system ). Case 4. τ 1 = τ = τ >0. When τ 1 = τ = τ,theneq.16)becomes λ + A 41 λ +B 41 λ + B 40 )e λτ + C 40 e λτ =0, ) A 41 = A 1, B 41 = B 1 + C 1, B 40 = C 0, C 40 = D 0. Multiplying by e λτ,eq.)becomes B 41 λ + B 40 + λ + A 41 λ ) e λτ + C 40 e λτ =0. 3) Let λ = i > 0) be the root of Eq. 3). Then { C 40 ) cos τ+ A 41 sin τ= B 40, C 40 ) sin τ A 41 cos τ= B 41. It follows that cos τ= B 40 A 41 B 41 ) + B 40 C A 41 C 40 Then we have, cos τ= B 40 A 41 B 41 ) + B 40 C A 41 C e e e 41 + e 40 =0, 4) e 40 = C 40 C 40 B 40), e43 =A 41 B 41, 5) e 41 =B 40 C 40 A 41 B 41 B 40 ) A 41 B 40 B 41 C 40 ) A 41 C 40, 6) e 4 = A 4 41 C 40 B 41A 41 B 40 B 41 C 40 ) B 40 A 41 B 41 ). 7) Let = v,theneq.4)becomes v 4 + e 43 v 3 + e 4 v + e 41 v + e 40 =0. 8) The discussion of the roots of Eq. 8) is similar to that in []. Denote f v)=v 4 + e 43 v 3 + e 4 v + e 41 v + e 40. 9)

7 LiuAdvances in Difference Equations 014, 014:314 Page 7 of 19 From Eq. 9), we have Set f v)=4v 3 +3e 43 v +e 4 v + e 41. 4v 3 +3e 43 v +e 4 v + e 41 =0. 30) Let z = 3e v.theneq.9)becomes Define z 3 + α 41 z + α 40 =0, 31) α 41 = 1 e e 43, α 40 = 1 3 e e 4e 43 + e 41. ) ) 3 α40 α41 γ 41 = +, γ 4 = 1 + 3i, 3) 3 z 1 = 3 α 40 + γ α 40 γ 41, 33) z = 3 α 40 + γ 41 γ α 40 γ 41 γ4, 34) z 3 = 3 α 40 + γ 41 γ4 + 3 α 40 γ 41 γ 4, 35) v i = z i 3e 43, i =1,,3. 36) 4 Based on the conditions H 1 )andh ), we can conclude that e 40 >0.Thenwehavethe following results according to Lemma. in []. Lemma 1 For Eq.9), i) if γ 41 0, then Eq.9) has positive root if and only if v 1 >0and f v 1 )<0; ii) if γ 41 <0, then Eq.9) has positive root if and only if there exists at least one v {v 1, v, v 3 }, such that v >0and f v ) 0. In what follows, we assume that H 41 ): the coefficients satisfy one of the following conditions in α -β :α ) γ 41 0, v 1 >0,andf v 1 )<0;β ) γ 41 <0andthereexistsatleastone v {v 1, v, v 3 },suchthatv >0andf v ) 0. If the condition H 41 )holds,eq.8) has at least one positive root v 0.Thus,Eq.4)has at least one positive root 0 = v 0 and the corresponding critical value of the delay is Let τ k = 1 arccos B 40 A 41 B 41 )0 + B 40C A jπ, k =0,1,,... C 40 0 τ 0 = 1 arccos B 40 A 41 B 41 )0 + B 40C A C 40

8 LiuAdvances in Difference Equations 014, 014:314 Page 8 of 19 When τ = τ 0,Eq.3)hasapairofpurelyimaginaryroots±i 0. Differentiating both sides of Eq. 3)withrespecttoτ,weget Thus, [ ] dλ 1 λ + A 41 )e λτ + B 41 = dτ C 40 λe λτ λ 3 + A 41 λ )e τ λτ λ. [ dλ Re dτ ] 1 = P 4RQ 4R + P 4IQ 4I τ=τ 0 Q 4R +, Q 4I P 4R = A 41 cos τ sin τ B 41, 37) P 4I = A 41 sin τ cos τ 0 0, 38) Q 4R = A 41 0 cos τ C ) 3 sin τ0 0, 39) Q 4I = A 41 0 sin τ C ) 3 cos τ ) Obviously, if the condition H 4 ): P 4R Q 4R +P 4I Q 4I 0 holds, the transversality condition is satisfied. Based on the discussion above and according to the Hopf bifurcation theorem in [1] we have the following results. Theorem 3 Suppose that the conditions H 1 ), H ), H 41 ), and H 4 ) hold. The positive equilibrium E x, y ) of system ) is asymptotically stable for τ [0, τ 0 ). System ) undergoes a Hopf bifurcation at the positive equilibrium E x, y ) of system ) when τ = τ 0 and a family of periodic solutions bifurcate from the positive equilibrium E x, y ) of system ). Case 5. τ 1 >0,τ 0, τ 0 ). We consider Eq. 16)withτ in its stable interval and τ 1 is considered as the bifurcation parameter. Let λ = i 1 1 > 0) be the root of Eq. 16). Then { M51 sin τ M 5 cos τ 1 1 = N 51, M 51 cos τ 1 1 M 5 sin τ 1 1 = N 5, with M 51 = B 1 1 D 0 sin τ 1, M 5 = D 0 cos τ 1, 41) N 51 = 1) C1 1 sin τ 1 C 0 cos τ 1, 4) N 5 = A 1 1 C 1 1 cos τ 1 + C 0 sin τ 1. 43) Then one can obtain ) ) e e51 1 cos τ 1 + e 5 sin τ 1 =0, 44)

9 LiuAdvances in Difference Equations 014, 014:314 Page 9 of 19 e 50 1 ) = 1 ) 4 + A 1 B 1 + C 1 ) 1 ) + C 0 D 0, 45) e 51 1 ) =A1 C 1 C 0 ) 1), 46) e 5 1 ) = C1 1 ) 3 +B1 D 0 A 1 C 0 ) 1. 47) In order to give the main results in this paper, we make the following assumption. H 51 ): Eq. 44) has at least finite positive roots. If the condition H 51 ) holds, we denote the roots of Eq. 44)as 11, 1,..., 1k.Foreveryfixed 1i i =1,,...,k), the corresponding critical value of the delay is τ j) 1i = 1 { } M51 N 5 + M 5 N 51 1i arccos M51 + M 5 1 = 1i + jπ 1i, i =1,,...,k, j =0,1,,... Let τ10 j) = min{τ 1i i =1,,...,k, j =0,1,,...}.Whenτ 1 = τ10,eq.16)hasapairofpurely imaginary roots ±i1 for τ 0, τ 0 ). Differentiating Eq. 16) withrespecttoτ 1,onecan obtain Thus, [ ] dλ 1 = λ + A 1 + B 1 e λτ 1 +C 1 τ C 0 τ C 1 )e λτ τ D 0 e λτ 1+τ ). dτ 1 B 1 λ e λτ 1 + D 0 λe λτ 1+τ ) [ dλ dτ 1 ] 1 τ 1 =τ 10 = P 5RQ 5R + P 5I Q 5I Q 5R +, Q 5I P 5R = A 1 +C 1 τ C 1 τ C 0 ) cos τ 1 48) + B 1 + τ D 0 cos τ 1 ) cos τ 10 1 τ D 0 sin τ 1 sin τ 10 1, 49) P 5I =1 C 1 τ C 1 τ C 0 ) sin τ 1 50) B 1 + τ D 0 cos τ 1 ) τ D 0 sin τ 1 cos τ 10 1, 51) Q 5R = D 0 1 sin τ 1 B ) 1 ) 1 cos τ D 01 cos τ 1 sin τ 10 1, 5) Q 5I = ) B 1 1 D0 1 sin τ 1 ) sin τ D 01 cos τ 1 cos τ ) Obviously, if the condition H 5 ): P 5R Q 5R +P 5I Q 5I 0holds,thenRe[ dλ dτ 1 ] 1 τ 1 =τ10 0.Thus, according to the Hopf bifurcation theorem in [1] we have the followingresult. Theorem 4 Suppose that the conditions H 1 ), H 51 ), and H 5 ) hold and τ 0, τ 0 ). The positive equilibrium E x, y ) of system ) is asymptotically stable for τ 1 [0, τ10 ). System ) undergoes a Hopf bifurcation at the positive equilibrium E x, y ) of system ) when τ 1 = τ 10 and a family of periodic solutions bifurcate from the positive equilibrium E x, y ) of system ).

10 Liu Advances in Difference Equations 014, 014:314 Page 10 of 19 Case 6. τ >0,τ 1 0, τ 10 ). We consider Eq. 16)withτ 1 in its stable interval and τ is considered as the bifurcation parameter. Let λ = i > 0) be the root of Eq. 16). Then with { M61 sin τ + M 6 cos τ = N 61, M 61 cos τ M 6 sin τ = N 6, M 61 = C 1 D 0 sin τ 1, M 6 = C 0 + D 0 cos τ 1, 54) N 61 = ) B1 sin τ 1, N 6 = A 1 B 1 cos τ 1. 55) It follows that ) ) e 60 + e61 cos τ1 + e 6 sin τ 1 =0, 56) ) ) e 60 = 4 + A 1 + B 1 ) ) C 1 C 0 + D 0, 57) ) ) e 51 =A1 B 1 C0 D 0, 58) ) ) e 5 = B1 3 +C1 D 0. 59) SimilartoCase5,wemakethefollowingassumption.H 61 ): Eq. 56) has at least finite positive roots. We denote the roots of Eq. 56)as 1,,..., k.foreveryfixed i i = 1,,...,k), the corresponding critical value of the delay is τ j) i = 1 { } M61 N 6 + M 6 N 61 i arccos M61 + M 6 = i + jπ i, i =1,,...,k, j =0,1,,... Let τ0 j) = min{τ i i =1,,...,k, j =0,1,,...}.Whenτ = τ0,eq.16)hasapairofpurely imaginary roots ±i for τ 1 0, τ 10 ). Differentiating Eq. 16) withrespecttoτ,onecan obtain Thus, [ ] dλ 1 = λ + A 1 +B 1 τ 1 B 1 λ)e λτ 1 + C 1 e λτ τ 1 D 0 e λτ 1+τ ) dτ [ dλ dτ ] 1 τ =τ 0 C 1 λ + C 0 λ)e λτ + D 0 λe λτ 1+τ ) = P 6RQ 6R + P 6I Q 6I Q 6R +, Q 6I τ λ. P 6R = C 1 τ 1 D 0 cos τ 1 ) cos τ 0 + τ 1D 0 sin τ 1 sin τ 0 60) + A 1 + B 1 cos τ 1 + τ 1B 1 sin τ 1, 61) P 6I = τ 1 D 0 cos τ 1 C 1) sin τ 0 + τ 1D 0 sin τ 1 cos τ 0 6)

11 Liu Advances in Difference Equations 014, 014:314 Page 11 of 19 + B 1 sin τ 1 τ 1B 1 cos τ 1, 63) Q 6R = C 0 + D 0 cos τ 1 ) sin τ 0 + D 0 sin τ 1 C 1 ) ) cos τ 0, 64) Q 6I = C 0 + D 0 cos τ 1 ) cos τ 0 D 0 sin τ 1 C 1 ) ) sin τ 0. 65) Similar as in Case 5, we can conclude that if the condition H 6 ): P 6R Q 6R + P 6I Q 6I 0 holds, then Re[ dλ dτ ] 1 τ =τ0 0. Thus, according to the Hopf bifurcation theorem in [1]we have the following results. Theorem 5 Suppose that the conditions H 1 ), H 61 ), and H 6 ) hold and τ 1 0, τ 10 ). The positive equilibrium E x, y ) of system ) is asymptotically stable for τ [0, τ0 ). System ) undergoes a Hopf bifurcation at the positive equilibrium E x, y ) of system ) when τ = τ0 and a family of periodic solutions bifurcate from the positive equilibrium E x, y ) of system ). 3 Direction and stability of bifurcated periodic solutions In this section, we investigate the direction of the Hopf bifurcation and the stability of bifurcated periodic solutions of system )withrespecttoτ for τ 1 0, τ 10 ). Let τ = τ 0 + μ, μ R, so that the Hopf bifurcation occurs at μ = 0. Without loss of generality, we assume that τ 1 < τ 0,τ 1 0, τ 10 ). Let u 1 t)=xt) x, u t)=yt) y, and rescale the time delay t t/τ ), then system )becomes ut)=l μ u t + Fμ, u t ), 66) ut)=u 1 t), u t)) T R and L μ : C R, F : R C R are given, respectively, by L μ φ = τ0 + μ) A φ0) + B φ τ ) ) 1 τ 0 + C φ 1) and Fμ, φ)= τ 0 + μ) F 1, F ) T, with φθ)=φ 1 θ), φ θ)) T C[ 1, 0], R ), ) ) ) A a 11 a 1 =, B b 11 0 =, C 0 0 = c 1 c and F 1 = g 1 φ1 0) + g φ 1 0)φ 0) + g 3 φ 0) + g 4φ 1 0)φ 1 τ ) 1 τ0 + g 5 φ1 0)φ 0) 67) + g 6 φ 1 0)φ 0) + g 7φ1 3 0) + g 8φ 3 0) +, 68) F = h 1 φ 1 1) + h φ 1 1)φ 0) + h 3 φ 1 1)φ 1) + h 4 φ 1 1)φ 0) 69) + h 5 φ1 1)φ 1) + h 6 φ1 3 1) +. 70)

12 Liu Advances in Difference Equations 014, 014:314 Page 1 of 19 Using the Riesz representation theorem, there exists a matrixfunctionηθ, μ), θ [ 1, 0] R whose elements are of bounded variation, such that L μ φ = 0 1 In fact, we can choose dηθ, μ)φθ), φ C [ 1, 0], R ). 71) τ0 + μ)a + B + C ), θ =0, τ0 ηθ, μ)= + μ)b + C ), θ [ τ 1 τ,0), 0 τ0 + μ)c, θ 1, τ 1 τ0 ), 0, θ = 1. For φ C[ 1, 0]), we define Aμ)φ = { dφθ), 1 θ <0, dθ 0 1 dηθ, μ)φθ), θ =0 and Rμ)φ = { 0, 1 θ <0, Fμ, φ), θ =0. Then system 66) canbetransformedintothefollowingoperatorequation: ut)=aμ)u t + Rμ)u t, 7) u t = ut + θ)=u 1 t + θ), u t + θ)) for θ [0, 1]. For ϕ C 1 [0, 1], R ) ), we define the adjoint operator A of A: A ϕ)= { dϕs) ds, 0<s 1, 0 1 dηt s, μ)ϕ s), s =0, and a bilinear inner product: ϕs), φθ) = ϕ T 0)φ0) 0 θ= 1 θ ξ=0 ϕ T ξ θ) dηθ)φξ) dξ, 73) ηθ)=ηθ,0). Let ρθ) =1,ρ ) T e i τ 0 θ be the eigenvectors of A0) corresponding to i τ 0 and ρ s) =D1, ρ )ei τ 0 s be the eigenvectors of A corresponding to i τ 0.Byasimple computation, we can get c 1 ρ = i c e iτ 0, ρ = i + c e iτ 0 a 11 + b 11 e iτ 1. From Eq. 73), we choose D = [ 1+ρ + ρ + τ 1b 11 e iτ 1 + ρ τ 0 e iτ 0 c 1 + c 3 ρ ) ] 1, such that ρ, ρ =1, ρ, ρ =0.

13 Liu Advances in Difference Equations 014, 014:314 Page 13 of 19 In the remainder of this section, following the algorithms given in [1]andusingasimilar computation process to that in [3], we obtain the following coefficients that determine the properties of the Hopf bifurcation: g 0 =τ0 [g D 1 + g ρ ) 0) + g 3 ρ ) 0) ) + g4 ρ 1) τ ) 1 τ0 + ρ h1 ρ 1) 1) ) 74) + h ρ 1) 1)ρ ) 0) + h 3 ρ 1) 1)ρ ) 1) )], 75) g 11 = τ0 [g D 1 + g ρ ) 0) + ρ ) 0) ) +g 3 ρ ) 0) ρ ) 0) + g 4 ρ 1) τ 1 τ 0 ) 76) + ρ 1) τ )) 1 τ0 + ρ h1 ρ 1) 1)ρ 1) 1) + h ρ 1) 1) ρ ) 0) 77) + ρ 1) 1)ρ ) 0) ) + h 3 ρ 1) 1) ρ ) 1) + ρ 1) 1)ρ ) 1) ))], 78) g 0 =τ0 [g D 1 + g ρ ) 0) + g 3 ρ ) 0) ) + g4 ρ 1) τ ) 1 τ0 + ρ h1 ρ 1) 1) ) 79) + h ρ 1) 1) ρ ) 0) + h 3 ρ 1) 1) ρ ) 1) )], 80) [ g 1 =τ0 D g 1 W 1) 1) 11 0) + W 0 0)) + g W 1) 11 0)ρ) 0) + 1 W 1) 0 0) ρ) 0) 81) + W ) 11 0) + 1 W ) 0 0) ) + g 3 W 1) 11 0)ρ) 0) + W ) 0 0) ρ) 0) ) 8) + g 4 W 1) 11 0)ρ1) τ 1 τ 0 + W 1) 11 τ ) 1 τ0 + 1 W 1) 0 ) + 1 W 1) 0 0) ρ1) τ 1 τ )) 1 τ0 τ 0 ) 83) + g 5 ρ ) 0) + ρ ) 0) ) + g 6 ρ ) 0) ) +ρ ) 0) ρ ) 0) ) +3g 7 84) +3g 8 ρ ) 0) ) ρ ) 0) + ρ h 1 W 1) 11 1)ρ1) 1) + W 1) 0 1)) 85) + h W 1) 11 1)ρ) 0) + 1 W 1) 0 1) ρ) 0) + W ) 11 0)ρ1) 1) 86) + 1 ) W ) 0 0) ρ1) 1) + h 3 W 1) 11 1)ρ) 1) + 1 W 1) 0 1) ρ) 1) 87) + W ) 11 1)ρ1) 1) + 1 ) W ) 0 1) ρ1) 1) + h 4 ρ 1) 1) ) ρ ) 0) 88) +ρ 1) 1)ρ ) 0) ρ 1) 1) ) + h 5 ρ 1) 1) ) ρ ) 1) 89) +ρ 1) 1)ρ ) 1) ρ 1) 1) ) +3h 6 ρ 1) 1) ) )] ρ 1) 1), 90) with W 0 θ)= ig 0ρ0) τ 0 W 11 θ)= ig 11ρ0) τ 0 e iτ 0 θ + iḡ 0 ρ0) 3τ0 e iτ 0 θ + E 0 e iτ 0 θ, 91) e iτ 0 θ + iḡ 11 ρ0) τ0 e iτ 0 θ + E 11, 9)

14 Liu Advances in Difference Equations 014, 014:314 Page 14 of 19 E 0 and E 11 can be computed by the following equations, respectively: with ) i E 0 = 11 b 11 e iτ 1 1 ) a 1 E 1) 0 c 1 e iτ 0 i c e iτ 0 E ), 0 93) ) 1 ) a 11 + b 11 a 1 E 1) 11 E 11 = c 1 c E ), 11 94) E 1) 0 = g 1 + g ρ ) 0) + g 3 ρ ) 0) ) + g4 ρ 1) τ ) 1 τ0, 95) E ) 0 = h 1 ρ 1) 1) ) + h ρ 1) 1)ρ ) 0) + h 3 ρ 1) 1)ρ ) 1), 96) E 1) 11 =g 1 + g ρ ) 0) + ρ ) 0) ) +g 3 ρ ) 0) ρ ) 0) 97) + g 4 ρ 1) τ ) 1 τ0 + ρ 1) τ )) 1 τ0, 98) E ) 11 =h 1ρ 1) 1)ρ 1) 1) + h ρ 1) 1) ρ ) 0) + ρ 1) 1)ρ ) 0) ) 99) + h 3 ρ 1) 1) ρ ) 1) + ρ 1) 1)ρ ) 1) ). 100) Therefore, we can calculate the following values: C 1 0) = i τ 0 g 11 g 0 g 11 g 0 3 ) + g 1, 101) μ = Re{C 10)} Re{λ τ0 10) )}, β =Re { C 1 0) }, 103) T = Im{C 10)} + μ Im{λ τ0 )} τ0. 104) Based on the discussion above, we can obtain the following results for system ). Theorem 6 The direction of the Hopf bifurcation is determined by the sign of μ : if μ >0 μ <0),the Hopf bifurcation is supercritical subcritical). The stability of bifurcating periodic solutions is determined by the sign of β : if β <0β >0),the bifurcating periodic solutions are stable unstable). The period of the bifurcating periodic solutions is determined by the sign of T : if T >0T <0),the period of the bifurcating periodic solutions increases decreases). 4 Numerical example In order to verify the analytic results obtained above and depict the Hopf bifurcation phenomenonofsystem ), we give some numerical simulations in this section. The study of system ) in this paper is restricted only to a theoretical analysis, therefore, for the choice of the value of the parameters in system ), we only consider the conditions mentioned in Section and the simulation effect. We hope that it may be helpful for experimental studies of the real situation. To this end, we choose a set of parameters randomly which can

15 Liu Advances in Difference Equations 014, 014:314 Page 15 of 19 Figure 1 E is stable for τ 1 =.35<τ 10 = Figure E is unstable for τ 1 =3.5>τ 10 = Figure 3 E is stable for τ = 1.08 < τ 0 = describe the Hopf bifurcation phenomenon of system ) and get the following system: { dxt) 0.yt) xt)+0.5yt) )xt), =1 0.5xt τ dt 1 ) dyt) =1. 0.3yt τ ) dt xt τ )+0.6 )yt), 105) which has a positive equilibrium E 1.058, 6.638). Then we obtain A 10 =0.7513>0, A 11 =1.705>0andC 0 D 0 = < 0. Thus, the conditions H 11 )andh 1 )holds. For τ 1 >0,τ = 0. By a simple computation, we get 10 =0.3848,τ 10 =.7488.From Theorem 1, thepositiveequilibriume 1.058, 6.638) of system 105) is asymptotically stable when τ 1 < τ 10.ThisisillustratedbyFigure1. As can be seen from Figure 1, when τ 1 =.35<τ 10 =.7488, the positive equilibrium E 1.058, 6.638) of system 105) is asymptotically stable. Once τ 1 passes through the critical value τ 10,thepositiveequilibrium E 1.058, 6.638) of system 105) loses stability and a family of periodic solutions bifurcate from the positive equilibrium E 1.058, 6.638), which can be shown as in Fig-

16 Liu Advances in Difference Equations 014, 014:314 Page 16 of 19 Figure 4 E is unstable for τ = 1.51 > τ 0 = Figure 5 E is stable for τ =1.04<τ 0 = Figure 6 E is unstable for τ = > τ 0 = ure. AsshowninFigure, wechooseτ 1 =3.5>τ 10 =.7488,thenE 1.058, 6.638) is unstable and a Hopf bifurcation occurs and a family of periodic solutions bifurcate from E 1.058, 6.638). Similarly, we obtain 0 =1.455andτ 0 =1.1705whenτ 1 =0, τ > 0. According to Theorem, thepositiveequilibriume 1.058, 6.638) of system 105) is asymptotically stable when the value of τ is below the critical value τ 0 and E 1.058, 6.638) becomes unstable when the value of τ is above the critical value τ 0. This property is illustrated by Figures 3-4. Consider the case τ 1 = τ = τ > 0. By some complex computations, we obtain 0 = and τ 0 =1.106.FromTheorem3, we can conclude that the positive equilibrium E 1.058, 6.638) of system 105) is asymptotically stable when τ [0, 1.106) and a Hopf bifurcation occurs when τ >1.106=τ 0.Figure5 shows that the positive equilibrium E 1.058, 6.638) of system 105) is asymptotically stable when τ =1.04 [0, 1.106). Then the positive equilibrium E 1.058, 6.638) of system 105)becomesunstablewhen τ = > = τ 0, which can be shown in Figure 6.

17 Liu Advances in Difference Equations 014, 014:314 Page 17 of 19 Figure 7 E is stable for τ 1 =.56<τ 10 =.7138 and τ =0.75. Figure 8 E is unstable for τ 1 = > τ 10 =.7138 and τ =0.75. Figure 9 E is stable for τ =1.05<τ 0 = and τ 1 =0.85. Now we consider τ 1 >0,τ =0.75 0, τ 0 ). We can obtain 1 =1.133andthenwe obtain τ 10 =.7138.Figure7 shows that the positive equilibrium E 1.058, 6.638) of system 4) is asymptotically stable when τ =.56<.7138=τ10 and Figure 8 shows that there is a Hopf bifurcation occurs at the positive equilibrium E 1.058, 6.638) of system 105) and a family of periodic solutions bifurcate from E 1.058, 6.638) when τ = > τ 10 =.7138.Similarly,wehave =0.4950andτ 0 =1.0988whenτ >0, τ 1 =0.85 0, τ 10 ). The corresponding waveform and plots are shown in Figures In addition, we obtain λ τ 0 )= i and C 10) = i by some complex computations. Further we have μ = > 0, β = <0,T = 7.569<0. According to Theorem 6, we can conclude that the Hopf bifurcation with respect to τ with τ 1 =0.85 0, τ 10 ) is supercritical, the bifurcating periodic solutions are stable and decrease. Since the bifurcating periodic solutions of system 105) are stable, we know that the two species in system 105) can coexist in an oscillatory mode.

18 Liu Advances in Difference Equations 014, 014:314 Page 18 of 19 Figure 10 E is unstable for τ =1.15>τ 0 = and τ 1 = Conclusion In the present paper, a delayed predator-prey system with modified Leslie-Gower and Beddington-DeAngelis functional response is considered. We incorporate the feedback delays of the prey and the predator into the predator-prey system considered in the literature [1] and get a predator-prey system with Beddington-DeAngelis functional response and two delays. The main purpose of the paper is to investigate the effect of the two delays on the system. By choosing the diverse delay as a bifurcation parameter, we show that the complex Hopf bifurcation phenomenon at the positive equilibrium of the system can occur as the diverse delay crosses some critical values. Furthermore, the properties of the Hopf bifurcation such as direction and stability are determined. From the numerical example, we can know that the two species in system )couldcoexistinanoscillatorymode under some certain conditions. This is valuable from the point of view of ecology. It should be pointed out that Gakkhar and Singh [19] have earlier considered the Hopf bifurcation of another modified Leslie-Gower predator-prey system with Holling- II functional response and two delays. But the predator-prey system with a Beddington- DeAngelis functional response considered in this paper is more general and it can reflect the dynamic relationship between the predator and the prey more effectively and it may be more helpful for experimental studies of the real situation. In addition, the global stability of the positive equilibrium and global existence of the Hopf bifurcation are disregarded in the paper. We leave this for further investigation. Competing interests The author declares that she has no competing interests. Author s contributions The author mainly studied the Hopf bifurcation of a predator-prey system with two delays and wrote the manuscript carefully. Acknowledgements The author would like to thank the editor and the two anonymous referees for their constructive suggestions on improving the presentation of the paper. This work was supported by Natural Science Foundation of the Higher Education Institutions of Anhui Province KJ013A003, KJ013B137). Received: May 014 Accepted: 1 September 014 Published: 15 Dec 014 References 1. Beddington, JR: Mutual interference between parasites or predators and its effect on searching efficiency. J. Anim. Ecol. 44, ). DeAngelis, DL, Goldstein, RA, O Neill, RV: A model for trophic interaction. Ecology 56, ) 3. Ko, W, Ryu, K: Analysis of diffusive two-competing-prey and one-predator systems with Beddington-DeAngelis functional response. Nonlinear Anal. 71, ) 4. Liu, SQ, Zhang, JH: Coexistence and stability of predator-prey model with Beddington-DeAngelis functional response and stage structure. J. Math. Anal. Appl. 34, )

19 Liu Advances in Difference Equations 014, 014:314 Page 19 of Liu, H, Liu, M, Wang, K, Wang, Y: Dynamics of a stochastic predator-prey system with Beddington-DeAngelis functional response. Appl. Math. Comput. 19, ) 6. Liu, M, Wang, K: Global stability of stage-structured predator-prey models with Beddington-DeAngelis functional response. Commun. Nonlinear Sci. Numer. Simul. 16, ) 7. Li, WL, Wang, LS: Stability and bifurcation of a delayed three-level food chain model with Beddington-DeAngelis functional response. Nonlinear Anal., Real World Appl. 10, ) 8. Baek, H: Qualitative analysis of Beddington-DeAngelis type impulsive predator-prey models. Nonlinear Anal., Real World Appl. 11, ) 9. Cui, JA, Takeuchi, Y: Permanence, extinction and periodic solution of predator-prey system with Beddington-DeAngelis functional response. J. Math. Anal. Appl. 317, ) 10. Li, HY, Takeuchi, Y: Dynamics of the density dependent predator-prey system with Beddington-DeAngelis functional response. J. Math. Anal. Appl. 374, ) 11. Liu, S, Beretta, E, Breda, D: Predator-prey model of Beddington-DeAngelis type with maturation and gestation delays. Nonlinear Anal., Real World Appl. 11, ) 1. Yu, SB: Global stability of a modified Leslie-Gower model with Beddington-DeAngelis functional response. Adv. Differ. Equ. 014, ) 13. Bianca, C, Ferrara, M, Guerrini, L: Qualitative analysis of a retarded mathematical framework with applications to living systems. Abstr. Appl. Anal. 013,Article ID ) 14. Jana, S, Kar, TK: Modeling and analysis of a prey-predator system with disease in the prey. Chaos Solitons Fractals 47, ) 15. Bianca, C, Guerrini, L: On the Dalgaard-Strulik model with logistic population growth rate and delayed-carrying capacity. Acta Appl. Math. 18, ) 16. Kar, TK, Ghorai, A: Dynamic behaviour of a delayed predator-prey model with harvesting. Appl. Math. Comput. 17, ) 17. Cui, GH, Yan, XP: Stability and bifurcation analysis on a three-species food chain system with two delays. Commun. Nonlinear Sci. Numer. Simul. 16, ) 18. Bianca, C, Guerrini, L: Existence of limit cycles in the Solow model with delayed-logistic population growth. Sci. World J. 014, Article ID ) 19. Gakkhar, S, Singh, A: Complex dynamics in a prey predator system with multiple delays. Commun. Nonlinear Sci. Numer. Simul. 17, ) 0. Bianca, C, Ferrara, M, Guerrini, L: The time delays effects on the qualitative behavior of an economic growth model. Abstr. Appl. Anal. 013,Article ID ) 1. Hassard, BD, Kazarinoff, ND, Wan, YH: Theory and Applications of Hopf Bifurcation. Cambridge University Press, Cambridge 1981). Li, XL, Wei, JJ: On the zeros of a fourth degree exponential polynomial with applications to a neural network model with delays. Chaos Solitons Fractals 6, ) 3. Bianca, C, Ferrara, M, Guerrini, L: The Cai model with time delay: existence of periodic solutions and asymptotic analysis. Appl. Math. Inf. Sci. 7, ) / Cite this article as: Liu: Dynamical analysis of a delayed predator-prey system with modified Leslie-Gower and Beddington-DeAngelis functional response. Advances in Difference Equations 014, 014:314