1 Introduction

In recent decades,it has been of great interest to investigate the dynamic interaction between predator species and prey species in both ecology and mathematical ecology(see[1]–[3]).Especially,two species predator-prey systems have been investigated by many scholars at home and abroad(see[4]–[10]),since the pioneering works by Lotka[11]and Volterra[12].However,two species predator-prey systems can describe only a small number of the phenomena that are commonly observed in nature(see[13]).Therefore,it is more realistic to consider a predator-prey system with three or more species to understand complex dynamical behaviors of multiple species predator-prey systems in the real world.Inspired by this idea,Yang and Jia[14]proposed the following three species predator-prey system with reserve area for prey and in the presence of toxicity:

where x(t),y(t)and z(t)denote the biomass densities of the prey species in the unreserve areas,the prey species in the reserve areas and the predator species at time t,respectively.K is the carrying capacity of the prey species in the unreserve areas;r1is the intrinsic growth rate of the prey species in the unreserve areas;r2is the birth rate of the prey species in the reserve areas;v1,v2and ϖ are the infection rates of the prey species in the unreserve areas,the prey species in the reserve areas and the predator species by an external toxic substance,respectively;q1and q2are the catchability coefficients of the prey species in the unreserved areas and the predator species,respectively;E is the effort applied to harvest the prey species and the predator species in the unreserve areas;σ1(σ2)is the rate at which the prey species in the unreserve(reserve)areas migrate into the reserve(unreserve)areas;d is the death rate of the predator species;a is the capturing rate of the predator species;β is the rate of conversing the prey species in the unreserve areas into the predator species;b is the half saturation rate of the predator species.

Time delays of one type or another have been incorporated into predator-prey systems by many scholars(see[15]–[21]),since delay differential equations exhibit much more complicated dynamics than ordinary ones.Motivated by the work above and considering that the consumption of the prey species by the predator species throughout its past history governs the present birth rate of the predator species,we incorporate the time delay due to the gestation of the predator species into system(1.1)and get the following predator-prey system with time delay:

where τ is the time delay due to the gestation of the predator species.

肝内胆管细胞癌起病隐匿，早期临床表现缺乏特异性，多数患者确诊时已是晚期。临床应用氟尿嘧啶、吉西他滨、奥沙利铂等药物治疗胆管细胞癌，但总体预后仍不满意。临床亟需新的治疗思路、治疗药物给以突破。

The rest of this work is organized in this pattern.In the next section,the existence of the Hopf bifurcation is investigated.In Section 3,based on the normal form method and center manifold theory,properties of the Hopf bifurcation are investigated.In Section 4,a numerical example is carried out in order to support the obtained theoretical predictions.The final section gives our conclusion.

2 Local Stability and Hopf Bifurcation

If

then(1.2)has a unique positive equilibrium P∗(x∗,y∗,z∗),where

are satis fied.

with

在这个课例中，老师在完成新课的教学后，没有以大功告成的心态放松对学生的观察，而是积极地寻找学生存在的动作问题，把握住最佳的教学契机，以设计复习课的形式为学生进行再次的示范和讲解，在学生的纠错练习中加深了他们对动作的理解和认识。

where

最后是比喻修辞格的使用。这是一处明喻。将妻子比作棋子，二者有“被困、被限制”的共同点，表现了妻子被生活、被思念围困的状态。此处以比喻之手法，写尽妻子想挣脱生活现状而不得的无力之感，这是千百年来所有囿于世俗将全身心的希望都寄予了丈夫借以改变自身处境的妻子的境遇。

Hence,by Routh-Hurwitz criterion,P∗(x∗,y∗,z∗)is asymptotically stable when the condition(H1)holds,that is,

The linearized system of system(1.2)at P∗(x∗,y∗,z∗)is given by

For τ＞ 0,substituting λ =iω (ω ＞ 0)into(2.2),we obtain{

It follows that

whereLet ω2= ω∗.Then

Based on the discussion of the distribution of the roots of(2.6)in[3],we suppose that

(2)实行橙色预警机制，执行标准到位：①基层队组在生产组织过程中，达到预警范围而没有处理时，专业科室在日常检查时要书面通知整改，不做考核；②因队组管理不善或科室超前管理不到位而产生影响的，对科队两级进行责任追究。

Then there existssuch that(2.2)has a pair of purely imaginary roots±iω0.Thus,

Differentiating(2.2)with respect to τ,we have

and

where

Hence,the transversality condition is satis fied if the condition

holds.Thus,according to the Hopf bifurcation theorem in[22],we have the following.

从分析结果可以看出，电机的温度分布中，最高温度产生在绕组部位为130℃。对于一般电机来说完全不是问题，也就是说该电机完全可以满足长时工作的需求。

Theorem 2.1 For system(1.2),if the conditions(H1)–(H3)hold,then the infected equilibrium P∗(x∗,y∗,z∗)of system(1.2)is asymptotically stable for τ∈ [0,τ0)and system(1.2)undergoes a Hopf bifurcation at the positive equilibrium P∗(x∗,y∗,z∗)when τ= τ0,where τ0 is defined in(2.7).

The characteristic equation of(2.1)is

3 Stability of the Bifurcating Periodic Solutions

Let τ= τ0+µ,µ ∈ R.Then µ =0 is the Hopf bifurcation value for system(1.2).Define the space of continuous real-valued functions as C=C([−1,0],R3).Let

The system(1.2)transforms to functional differential equation in C as

where ut=(u1(t),u2(t),u3(t))T ∈ C=C([−1,0],R3),

and

where

《天津职业技术师范大学学报》是天津职业技术师范大学主办的面向国内外公开发行的综合性学术期刊，刊号CN 12－1423/Z，ISSN 2095-0926。本刊主要刊登机械工程、自动化工程、电子工程、计算机科学与技术、职业技术教育科学及基础科学包含的各学科领域、各专业门类的理论与应用研究学术论文（包括综述报告）、调查报告和问题讨论等。欢迎大家踊跃投稿和订阅。

and

with

By the representation theorem,there exists a 3 × 3 matrix function η(θ,µ),θ ∈ [−1,0]such that

In view of(3.2),we choose

where δ is the Dirac delta function.

When τ=0,(2.2)becomes

For ϕ ∈ C([−1,0],R3),define

and

For φ ∈ C1([0,1],(R3)∗),define

(1)地震技术的优势是毋庸置疑的，但勘探部署应根据探区地质条件的复杂程度合理设置地震测线网度，不宜过大。对勘探程度较低的探区采用“先地震后钻探；先大测网扫描选定有利目标区，然后在有利目标区一定程度加密控制”的煤层气地震勘探思路无疑是正确的和科学的；

where ut(θ)=u(t+ θ)for θ∈ [−1,0].

Then system(3.1)is equivalent to

7）在本项目成功举办了“以标准化信息化为抓手，以创精品工程为目标，全面提升中国十七冶集团工程建设管理水平，打造中冶管廊第一品牌”为主题的管廊工程观摩交流会，得到了一致好评，起到了较好的示范和引领作用；

and a bilinear inner product

Safety and happiness walk hand in hand while civilization and beauty comesideby side.

where η(θ)= η(θ,0).Then A(0)and A∗ are adjoint operators.

Since ±iω0τ0are the eigenvalues of A(0)and A∗,respectively.We need to compute eigenvectors of A(0)and A∗ corresponding to+iω0τ0and −iω0τ0,respectively.Suppose that ρ(θ)=(1,ρ2,ρ3)Teiω0τ0θ is the eigenvector of A(0)belonging to+iω0τ0and ρ∗(s)=is the eigenvector of A∗(0)belonging to −iω0τ0.By a direct computation,we can get

From(3.5)we can get

Then we choose

虽然南通鹏越纺织能紧跟时代发展需求，创新营销的方式方法，使得自身的发展不断壮大。但也在调研中发现了不少问题，主要表现如下：

such that 〈ρ∗,ρ〉=1.

Next,we can obtain the coefficients by using the method introduced in[22]and a computation process similar as that in[17],[23]and[24]:

with

政府的任何一个站(办)所都要求员工具备一定的电脑操作水平，懂得基本的软件操作，经济发展办与城建办的员工更是需要具备一定的专业素养，能够较为系统地统计和分析数据，否则大部分人将存在“搭便车”的可能。既然短时间内无法打破干部职工内部老龄化的局面，那就必须正视它，在每年新招年轻干部职工的基础上必须对老干部职工特别是年龄介于35-45岁之间的干部职工进行学习再造。基层领导应当解放思想、开拓创新，利用各种媒介广泛动员，不断完善人才培养继续教育机制。与此同时，老干部职工自己应当有“不服老、不服输”的精神风貌，要敢于学习新事物、接受新事物。老干部职工因为工作经验丰富，对其进行学习再教育具有重要意义。

where E1and E2are given by the following equations,respectively

Further,we have

Then we can get the following coefficients which determine the properties of the Hopf bifurcation:

我院恶性肿瘤住院患者营养状态与营养干预情况调查…………………………………………………… 鲁良秀等（10）：1416

In conclusion,we have the following results.

Theorem 3.1 For system(1.2),Ifµ2＞ 0(µ2＜ 0),then the Hopf bifurcation is supercritical(subcritical);If β2 ＜ 0(β2 ＞ 0),then the bifurcating periodic solutions are stable(unstable);If T2＞0(T2＜0),then the bifurcating periodic solutions increase(decrease).

4 Numerical Example

In this section,we present a numerical example of system(1.2)to illustrate our obtained analytical findings.Let r1=1.5,K=4,σ1=2,σ2=2,v1=0.4,a=1,b=0.25,q1=0.1,E=1,r2=1,v2=0.4,β=0.8,d=0.01,ϖ=0.4,q2=0.2.Then we get the following special case of system(1.2):

By a direct computation using Matlab software package,we obtain the unique positive equilibrium P∗(0.8026,1.1112,1.6284).Further,we obtain

徐备说，前二十几年，国内印刷行业的发展其实是粗放式的，可以总结为利润非常高、门槛非常低，比较偏重于硬件投入，缺少对流程、软件、人才、管理的投入。那么未来的二十年会怎样？从国外的例子来看，印刷业会随着提高效率、稳定品质、降低成本这几个软环境的改变而改变，这也是爱克发希望跟中国的印刷企业一起去做的事情。比方说现在谈的最多的是智能化、工业互联网等，但印刷行业整体来看可能连工业2.0都还没有达到。这就需要各方一起努力，共同推动标准化、数据化、工业化建设。只有当这些基本要求都达到了，工业2.0才能在印刷行业里实现，而后进阶更高级的3.0、4.0。

Let τ=0.83 ∈ [0,τ0).We can see that the positive equilibrium P∗(0.8026,1.1112,1.6284)is locally asymptotically stable,which can be shown in Fig.4.1.This reveals that the densities of the three species in system(4.1)will tend to stabilization and this situation does not vary with the time delay τ∈ [0,τ0).However,once the time delay τ passes through the critical value τ0,then the positive equilibrium P∗(0.8026,1.1112,1.6284)will lose its stability and a Hopf bifurcation occurs.This property can be illustrated by Fig.4.2.From this we know that the densities of the three species in system(4.1)will oscillate in the vicinity of the three species,respectively.In addition,we also obtain λ′(τ0)=0.4011 − 0.0367i and C1(0)=−3.4254+1.6580i by some complex computations.Thus,based on(3.6),we haveµ2=8.54＞ 0,β2=−6.8508＜0 and T2=−0.5666＜0.It follows that the Hopf bifurcation at τ0is supercritical and stable,and the periodic solutions decrease.

5 Conclusions

In this paper,a delayed predator-prey system with reserve area for prey and in the presence of toxicity is proposed by introducing the time delay due to the gestation of the predator species into the model considered in[14].Compared with the work in[14],we mainly investigate the effect of the time delay on the proposed system.

The main results are given in terms of local stability of the positive equilibrium and the local Hopf bifurcation of the system.It is proved that the positive equilibrium is locally asymptotically stable when the time delay is suitable small and a Hopf bifurcation occurs once the value of the time delay is larger than the critical value τ0.Through the numerical simulation,we know that the species in the system can coexist in an oscillatory mode when the certain conditions are satis fied.This is valuable from the viewpoint of biology.

References

[1]Baek H.A food chain system with Holling type IV functional response and impulsive perturbations.Comput.Math.Appl.,2010,60(5):1152–1163.

[2]Cui G H,Yan X P.Stability and bifurcation analysis on a three-species food chain system with two delays.Commun.Nonlinear Sci.Numer.Simul.,2011,16(9):3704–3720.

[3]Meng X Y,Huo H F,Zhang X B,Xiang H.Stability and Hopf bifurcation in a three-species system with feedback delays.Nonlinear Dynam.,2011,64(4):349–364.

[4]Yu S B.Global stability of a modified Leslie-Gower model with Beddington-DeAngelis functional response.Adv.Difference Equ.,2014,84:1–14.

[5]Yue Q.Permanence for a modified Leslie-Gower predator-prey model with Beddington-DeAngelis functional response and feedback controls.Adv.Difference Equ.,2015,81:1–8.

[6]Zhu Y L,Wang K.Existence and global attractivity of positive periodic solutions for a predator-prey model with modified Leslie-Gower Holling-type II schemes.J.Math.Anal.Appl.,2011,384:400–408.

[7]Etoua R M,Rousseau C.Bifurcation analysis of a generalized Gause model with prey harvesting and a generalized Holling response function of type III.J.Differential Equations,2010,249(9):2316–2356.

[8]Yang W S.Global asymptotical stability and persistent property for a diffusive predator-prey system with modified Leslie-Gower functional response.Nonlinear Anal.Real World Appl.,2013,14(3):1323–1330.

[9]Banerjee M,Banerjee S.Turing instabilities and spatio-temporal chaos in ratio-dependent Holling-Tanner model.Math.Biosci.,2012,236(1):64–76.

[10]Liu J.Dynamical analysis of a delayed predator-prey system with modified Leslie-Gower and Beddington-DeAngelis functional response.Adv.Difference Equ.,2014,314:1–19.

[11]Lotka A J.Elements of Physical Biology.New York:Williams and Wilkins,1925.

[12]Volterra V.Variazionie fluttuazioni del numero dindividui in specie animali conviventi.Mem.Acad Licei,1926,2(2):31–113.

[13]Do Y,Baek H,Lim Y,Lim D.A three-species food chain system with two types of functional response.Abstract and Applied Analysis,2011,Article ID 934569,16pp.

[14]Yang H,Jia J W.Harvesting of a predator-prey model with reserve area for prey and in the presence of toxicity.J.Appl.Math.Comput.,2017,53:693–708.

[15]Zhang J F.Bifurcation analysis of a modified Holling-Tanner predator-prey model with time delay.Appl.Math.Model.,2012,36(3):1219–1231.

[16]Celik C.Stability and Hopf Bifurcation in a delayed ratio dependent Holling-Tanner type model.Appl.Math.Comput.,2015,255:228–237.

[17]Upadhyay R K.Agrawal R.Dynamics and responses of a predator-prey system with competitive interference and time delay.Nonlinear Dynam.,2016,83(1-2):821–837.

[18]Kar T K,Ghorai A.Dynamic behaviour of a delayed predator-prey model with harvesting.Appl.Math.Comput.,2011,217(22):9085–9104.

[19]Wei F Y,Fu Q Y.Hopf bifurcation and stability for predator-prey systems with Beddington-DeAngelis type functional response and stage structure for prey incorporating refuge.Appl.Math.Model.,2016,40(1):126–134.

[20]Xu R.Global stability and Hopf bifurcation of a predator-prey model with stage structure and delayed predator response.Nonlinear Dynam.,2012,67(2):1683–1693.

[21]Biswas S,Sasmal S K,Samanta S,Saifuddin M,Pal N,Chattopadhyay J.Optimal harvesting and complex dynamics in a delayed eco-epidemiological model with weak Allee effects.Nonlinear Dynam.,2017,87(3):1553–1573.

[22]Hassard B D,KazarinoffN D,Wan Y H.Theory and Applications of Hopf Bifurcation.in:London Mathematical Society Lecture Note Series,41.Cambridge-New York:Cambridge University Press,1981.

[23]Xu C J.Delay-induced oscillations in a competitor-competitor-mutualist Lotka-Volterra model.Complexity,2017,Article ID 2578043,12pp.

[24]Bianca C,Ferrara M,Guerrini L.The Cai model with time delay:existence of periodic solutions and asymptotic analysis.App.Math.Inf.Sci.,2013,7(1):21–27.