HOPF BIFURCATION ANALYSIS OF A PREDATOR-PREY SYSTEM WITH NON-SELECTIVE HARVESTING AND TIME DELAY

LI Zhen-wei,LI Bi-wen,LIU Wei,WANG Gan

(School of Mathematics and Statistics,Hubei Normal University,Huangshi 435002,China)

HOPF BIFURCATION ANALYSIS OF A PREDATOR-PREY SYSTEM WITH NON-SELECTIVE HARVESTING AND TIME DELAY

LI Zhen-wei,LI Bi-wen,LIU Wei,WANG Gan

(School of Mathematics and Statistics,Hubei Normal University,Huangshi 435002,China)

In this paper,we mainly study the Hopf bifurcation and the stability of modifi ed predator-prey biological economic system with nonselective harvesting and time delay.By using the stability and bifurcation theory of diff erential-algebraic system,the conditions for stability of the positive equilibrium point are obtained,let time delay as bifurcation parameter,the existence of Hopf bifurcation and direction of Hopf bifurcation are obtained.We have improved the Leslie-Gower predator-prey system,make the system which we established more practical,so the conclusions are made more scientifi c.

stability;Hopf bifurcation;time delay;non-selective;predator-prey system; periodic solution

1 Introduction

In recent years,the increasingly serious problem of environmental degradation and resource shortage,made the analysis and modeling of biological systems more interested. The predator-prey system played a crucial role among the relationships between the biological population,and it naturally attracted much attention both for mathematicians and biologists,especially on predator-prey systems with or without time delay.As we know, delay differential equation models exhibit much more complicated dynamics than differential equation models without delay,see[1–12].A lot of researchers studied the dynamics of predator-prey models with harvesting and obtained many dynamic behaviors,such as stability of equilibrium,Hopf bifurcation,periodic solution,Bogdanov-Takens bifurcation, Neimark-Sacker bifurcation,and so on,see[10–15].

In[16],Lucas studied the dynamic properties of the following Leslie-Gower predatorprey system

where x and y denote prey and predator population densities at time t,respectively,a,d, and k are positive constants that represent the prey intrinsic growth rate,predator mortality rate,and the maximum value of the per capita reduction rate of x due to y,respectively.

At present,economic profit is a very important factor for merchants,governments and even every citizen,so it is necessary to research biologicaleconomic systems,which are often described by differential-algebraic equations or differentialdifference-algebraic equations.

In 1954,Gordon[13]studied the effect of the harvest effort on ecosystem from an economic perspective and proposed the following economic theory:

Net Economic Revenue(NER)=Total Revenue(TR)-Total Cost(TC).

This provides theoretical evidence for the establishment of diff erential-algebraic equation.

Based on the economic theory as mentioned above and system(1.1),Liu and Fu[12] considered the following Leslie-Gower predator-prey system

They investigated the Hopf bifurcation of the above system without considering the effect of time delay and the harvesting of predator.

As is known to all,delay differential equation models exhibit much more complicated dynamics than ordinary differential equation models,see[1–12],as was pointed by Kuang [17]that any modelof species dynamics without delays is an approximation at best.When we considered the model with non-selective harvesting,namely at the same time there are also the harvesting of predator and harvesting of the prey in the model,it will be more in line with the actualsituation of the predator-prey systems.

Motivated by the above discussion,in this paper,by choosing the time delay as a bifurcation parameter and consider the predator-prey systems with non-selective harvesting, we investigate a modifi ed Leslie-Gower predator-prey systems with non-selective harvesting and time delay described by the following system

where p1＞ 0 and p2＞ 0 are harvesting reward per unit harvesting effort for unit prey and predator,respectively;c1and c2are harvesting cost per unit harvesting effort for prey and predator,respectively;m is the economic profit per unit harvesting effort.

In this paper,we mainly discuss the effects ofeconomic profit on the dynamics ofsystem (1.3)in the region R3+={(x,y,E)|x ＞ 0,y ＞ 0,E ＞ 0)}.

For convenience,we let

where Xt=(x,y)T.

The rest ofthe paperisarranged asfollows:in Section 2,the localstability ofthe positive equilibrium points are investigated by the corresponding characteristic equation of system (1.3).In Section 3,by using the normal form and Hopf bifurcation theorem,we study the Hopfbifurcation ofthe nonnegative equilibrium depending on the parameter where we show that the positive equilibrium loses its stability and system(1.3)exhibits Hopfbifurcation in the second section.In Section 4,the theoretical result is supplied by a numerical example. Finally,this paper ends with a brief discussion.

2 Local Stability Analysis of System

In this section,we discuss the local stability of a positive equilibrium for system(1.3). Now,we try to find all possible positive equilibrium points of system(1.3).A point Y0= (x0,y0,E0)is an equilibrium point of system(1.3)if and only if Y0satisfy the following equations

From(2.1),we can easy get E0satisfy

where Based on the root and coeffi cient relationship of equation and γ3＜ 0,we can find at least one positive root E0,so system(1.3)has at least one positive equilibrium point,where r1＞ E0,r2＞ E0.

Now,we derive the formula for determining the properties of the positive equilibrium point ofsystem(1.3).As in[13],first we consider the localparametric ψ ofthe third equation of system(1.3),which is defined as follows

where

is a smoothing mapping,that is

Then we can obtain the parametric system of system(1.3)as follows (

Noticing that g(ψ(Z(t)))=0,so we can get the linearized system of parametric system(2.3) at(0,0)as follows

From(2.4),we can obtain the characteristic equation of the linearized system of parametric system(2.2)at(0,0)as follows

By eq.(2.5),when τ=0,it is obvious that,then,two roots of eq. (2.5)has always negative teal parts,i.e.,the positive equilibrium point of system(1.3)is locally asymptotically stable.

Now,based on the above discussion,we study the local stability around the positive equilibrium point for system(1.3)and the existence of Hopf bifurcation occurring at the positive equilibrium point when τ＞ 0.

If iω is a root of eq.(2.5),and substituting iω (ω is a positive real number)into eq. (2.5),and separating the real and imaginary parts,two transcendental equations can be obtained as follows

Since sin(ω τ)2+cos(ω τ)2=1 and adding(2.6)and(2.7),we obtain

Substituting ω0into(2.6)and solving for τ,we get

Thus when τ= τn,the characteristic equation(2.5)has a pair of purely imaginary roots iω0.

Lemma 2.1Denote by λn(τ)= ηn(τ)+iωn(τ)the root of(2.5)such that ηn(τn)=0, ωn(τn)= ω0,n=0,1,2,···.Then the following transversality condition η′n(τn)is satisfied.

ProofDifferentiating eq.(2.5)with respect to τ,we obtain

Noting that

The proof is completed.

From the above analysis and[17,18],we have the following results.

Theorem 2.1(i)For system(1.3),its positive equilibrium point Y0is locally asymptotically stable for τ∈ [0,τ0)and unstable for τ＞ τ0.

(ii)System(1.3)undergoes Hopf bifurcation at the positive equilibrium point Y0for τ= τn,n=0,1,2,···.

3 Direction and the Stability of Hopf Bifurcation

In this section,we investigate the direction ofbifurcation and the stability ofbifurcation periodic orbits from the positive equilibrium point Y0of system(1.3)at τ= τ0by using the normalform approach theory and center manifold theory introduced by Hassard[15].

Now,we re-scare the time by

for simplicity,we continue to use Z said ˉZ,then the parametric system(2.3)of system (1.3)is equivalent to the following functional differential equation(FDE)system in C= C([?1,0],R2),

where Z(T)=(y1(t),y2(t))T,and Lμ:C → R,f:R × C → R are given,respectively,by

where

and φ =(φ1,φ2) ∈ C.By the Riesz representation theorem,there exists a matrix whose components are bounded variation functions θ ∈ [?1,0]such that

where

Then system(3.1)can be rewritten as

For ψ ∈ C1([0,1],(R2)?),the adjoint operator A?of A as

where ηTis the transpose of the matrix η.

For φ ∈ C1([?1,0],R2)and for ψ ∈ C1([0,1],(R2)?),in order to normalize the eigenvectors of operator A and adjoint operator A?,we define a bilinear inner product

where η(θ)= η(θ,0).It is easy to verify that A(0)and A?are a pair of adjoint operators.

By the discussion in Section 2,we know that ±iωτ0are eigenvalues of A(0).Thus they are also eigenvalues of A?.Next we calculate the eigenvector q(θ)of A associated to the eigenvalue iω τ0and the eigenvector q?(s)of A?associated to the eigenvalue ?iω τ0.Then it is not diffi cult to show that

where

Moreover,〈q?(s),q(θ)〉=1 and 〈q?(s),ˉq(θ)〉=0.

In the reminder of this section,we use the same notations as those in[15].We fi rst compute the coordinates to describe the center manifold C0at μ =0.Define

On the center manifold C0,we have

In fact,z and ˉz are local coordinates for center manifold C0in the direction of q and ˉq?. Note that W is real if ztis real.We consider only real solutions.For the solution zt∈ C0, since μ =0 and(3.3),we have

rewrite it as

where

From(3.3)and(3.8),we have

Rewrite(3.11)as

where

Substituting the corresponding series into(3.12)and comparing the coeffi cient,we obtain

Notice that

and(3.6)we obtain

According to(3.8)and(3.9),we know that

where

By(3.7),it follows that

That is,

Comparing the coeffi cients with(3.10),it follows that

Now we compute W20(θ)and W11(θ).From(3.11)and(3.15),we have that for θ ∈ [?1,0),

Comparing the coeffi cients with(3.13),we can obtain that

Substituting the above equalities into(3.14),it follows that

Solving(3.18),we have

In what follows,we seek appropriate E and F in(3.19).From(3.11)and(3.15),we have

where

Substituting(3.19)–(3.21)into(3.14)and noting that

We obtain

It is easy to obtain E and F from(3.22)and(3.23),that is

Therefore we can compute the following quantitieswhich determine the direction of Hopf bifurcation and stability of bifurcated periodic solutions of system(1.3)at the critical value τ0.

Theorem 3.1(i)The direction of Hopfbifurcation is determined by the sign ofμ2:the Hopfbifurcation is supercritical if μ2＞ 0 and the Hopf bifurcation is subcriticalif μ2＜ 0.

(ii)The stability of bifurcated periodic solution is determined by β2:the periodic solution are stable if β2＞ 0 and unstable if β2＜ 0.

(iii)The period ofbifurcation periodic solution is determined by t2:the period increase if t2＞ 0,decrease if t2＜ 0.

4 Numerical Simulations

As an example we consider the differential-algebraic predator-prey system(1.3)with the parameters r1=1.6,r2=1.3,b=k=m=0.5,p1=7,p2=6,c1=5,c2=3,that is,

And by the discussions in Section 2 and Section 3,we determine the stability of the positive equilibrium point and Hopf bifurcation.Here,for convenience,we only discuss one of the positive equilibrium point Y0of system(4.1),and others positive equilibrium points of system(4.1)can be similar studied.we can easily get Y0=(2.0053,3.1480,0.0256),and by computing,we get ω0=0.9942, τ0=0.6473.So by Theorem 2.1,the equilibrium point Y0is asymptotically stable when τ∈ [0,τ0)=[0,0.6473)and unstable when τ＞ 0.6473.

When τ=0,we can easily show that the positive equilibrium point

is asymptotically stable.

By the theory of Hassard[15],as it is discussed in former section,we also determine the direction of Hopf bifurcation and the other properties of bifurcating periodic solution.By computing,we can obtain the following values C1(0)=0.5303 ? 0.4428i, λ′(τ0)=1.6352+ 1.1431i,it follows that μ2= ?0.3243 ＜ 0,β2=1.0607 ＞ 0,t2=1.2643 ＞ 0,from which and Theorem 3.1 we conclude that the Hopf bifurcation of system(4.1)occurring at τ0=0.6473 is subcriticaland the bifurcating periodic solution exists when τcross τ0to the left and the bifurcating periodic solution is unstable.

By Theorem 3.1,the positive equilibrium point Y0of system(4.1)is locally asymptotically stable when τ=0.62 ＜ τ0as is illustrated by computer simulation in Fig.1.And periodic solutions occur from Y0when τ=0.682 ＞ τ0as is illustrated by computer simulation in Fig.2.

5 Discussion

Nowadays,economic profit is a very important factor for governments,merchants,and even citizen,and the harvested biologicalresources in the predator-prey systems are usually sold as commodities in the market in order to achieve the economic interest.So modelling and qualitative analysis for bio-economic system are necessary.

Compared with most other researches on dynamics ofpredator-prey population,see[1, 5,12,18],the main contribution ofthis paper lies in the following aspect.The predator-prey system we consider incorporate delay and non-selective harvesting,which could make our model more realistic and the analysis result in this paper is more scientific.So our paper provide a new ideal and a effi cacious method for the qualitative analysis of the Hopf bifurcation of the differential-algebraic biologicaleconomic system.In addition,stage structure, diffusion effects,disease effects may be incorporated into our bio-economic system,which would make the bio-economic system exhibit much more complicated dynamics.

[1]Li K,Wei J.Stability and Hopf bifurcation analysis of a prey-predator system with two delays[J]. Chaos Solitions Fract.,2009,42:2606–2613.

[2]Ma Y F.Global Hopf bifurcation in the Leslie-Gower predator-prey model with two delays[J]. Nonlinear Anal.Real World Appl.,2012,13:370–375.

[3]Sadhukhan D,Mondal B,Maiti M.Discrete age-structured population model with age dependent harvesting and its stability analysis[J].Appl.Math.Comput.,2008,201:631–639.

[4]Huo H F,Li W T.Existence and global stability of periodic solutions of a discrete predator-prey system with delays[J].Appl.Math.Comput.,2004,153:337–351.

[5]Zhang G D，Shen Y，Chen B S.Bifurcation analysis in a discrete diff erential-algebraic predator-prey system[J].Appl.Math.Model.,2014,38:4835–4848.

[6]Li X,Ruan S,Wei J.Stability and bifurcation in delay-diff erential equations with two delays[J].J. Math.Anal.Appl.,1999,236:254–280.

[7]Pan K,Li B W.Existence of positive periodic solution for two-patches predator-prey impulsive diff usion delay systems with functional response[J].J.Math.,.2010,30(1):183–190.

[8]Zhao H,Wang L,Ma C.Hopf bifurcation in a delayed Lotka-Volterra predator-prey system[J].Nonl. Anal.RWA.,2008,9(1):114–127.

[9]Wei J,Ruan S.Stability and bifurcation in a neural model with two delays[J].Phys.D.,1999,130: 255–272.

[10]Celik C.The stability and Hopf bifurcation for a predator-prey system with time delay[J].Chaos Solitons Fract.,2008,37:87–99.

[11]Pan S X.Asymptotic spreading in a Lotka-Volterra predator-prey system[J].J.Math.Anal.Appl., 2013,407:230–236.

[12]Liu W,Fu C J,Chen B S.Hopf bifurcation and center stability for a predator-prey biological economic modelwith prey harvesting[J].Commun.Nonlinear.Sci.,doi:10.1016/j.cnsns.2012.02.025.

[13]Gordon H S.Economic theory of a common property resource:the fi shery[J].J.Polit.,1954,62(2): 124–142.

[14]Li P L,Yu CC,Zeng XW.The qualitative analysis of a class ofpredator-prey system with functional responses[J].J.Math.,2006,26(2):217–222.

[15]Hassard N D,Kazarinoff Y H,Wan Y H.Theory and application of Hopf bifurcation[M].Cambridge: Lodon Mathematical Society Lecture Notes,vol.41,Cambridge Univ.Press,1981.

[16]Lucas W F.Modules in applied mathematics:diff erential equation models[M].New York:Springer, 1983.

[17]Kuang Y.Delay diff erential equations with applications in population dynamics[M].Boston:Academic Press,1993.

[18]Zhang G D,Shen Y,Yin Q,Sun J W.Passivity analysis for memristor-based recurrent neural networks with discrete and distributed delays[J].Neural Networks,2015,61:49–58.

[19]Zhang X,Zhang Q L.Bifurcation analysis and control of a class of hybrid biological economic models[J].Nonlinear Anal HS.,2009,3:578–87.

[20]Chen B S,Liao X X,Liu Y Q.Normal forms and bifurcations for the diff erential-algebraic systems (in Chinese)[J].Acta Math.Appl.Sinica.,2000,23(3):429–443.

[21]Lv X,Lu S P,Yan P.Existence and global attractivity of positive periodic solutions of Lotka-Volterra predator-prey systems with deviating arguments[J].Nonl.Anal.Real World Appl.,2010, 11:574–583.

一類帶無選擇性捕獲和時滯的捕食食餌系統(tǒng)的Hopf分支分析

李震威,李必文,劉 煒,汪 淦
(湖北師范學院數(shù)學與統(tǒng)計學院,湖北 黃石 435002)

本文主要研究了一個改進的帶時滯和無選擇捕獲函數(shù)的捕食-食餌生態(tài)經(jīng)濟系統(tǒng)的穩(wěn)定性和Hopf分支. 利用微分代數(shù)系統(tǒng)的穩(wěn)定性理論和分支理論, 得到了系統(tǒng)正平衡點穩(wěn)定性的條件, 以及當時滯τ 作為分支參數(shù)時系統(tǒng)產(chǎn)生Hopf分支的條件. 對Leslie-Gower捕食-食餌模型進行了一定程度的完善, 使得建立的模型更符合實際情況,因此得到的結(jié)論也更加科學.

穩(wěn)定性;Hopf分支; 時滯; 無選擇性; 捕食食餌系統(tǒng); 周期解

34D20;34K18;34C23

O29;O193

tion:34D20;34K18;34C23

A < class="emphasis_bold">Article ID:0255-7797(2017)02-0257-14

0255-7797(2017)02-0257-14

?Received date:2014-11-16 Accepted date:2015-02-26

Foundation item:Supported by the Funding Program of Higher School Outstanding Youth Scientifi c and Technological Innovation Team in Hubei of China(T201412).

Biography:Li Zhenwei(1991–),male,born at Qianjiang,Hubei,major in ordinary diff erential equations and control theory.