时滞比率依赖种群模型Hopf-Fold分支现象分析
2014-09-17夏晶,郭爽
夏 晶 ,郭 爽
(大庆师范学院)
0 引言
1977年,Freedman和 Waltman在文献[1]中提出了一类三维捕食-食饵模型:

这里x(t),y(t)和z(t)分别表示t时刻食饵,捕食者和顶层捕食者的数量;g(x)是食饵的内部增长函数;p(x)和q(x)分别是捕食者和顶层捕食者的功能反应增长函数;h,s>0分别是捕食者和顶层捕食者的死亡率;e,m>0分别是食饵和捕食者的转换率.许多作者对该模型进行了研究,具体参见文献[2-7].该文主要针对顶层捕食者方程中含有时滞项的情形进行分析,对该模型能否发生Hopf-Fold分支现象进行讨论.
1 Hopf-Fold分支
该文研究的模型如下:

对方程(2)进行非量纲化处理,得到:

其中x0(θ)=φ1(θ)≥0,y0(θ)=φ2(θ)≥0,z0(θ)=φ3(θ)≥ 0,θ∈ [-τ,0],x(0)>0,y(0)>0,z(0)>0,‖φ‖ =max{|φ(θ)|:θ∈[τ,0]},φ(θ)=(φ1,φ2,φ3)∈C([-1,0],R3)及这里a,k,c,m,d,f,r,s,e都是正参数.假设方程(3)的正平衡点存在,并记为E*(x*,y*,z*),其中

在E*处考虑相应的线性化系统,特征值λ满足如下的特征方程:

其中:a2=-m11-m22,a1=m11m22-m12m21,b1=m23m32andb0=m11m23n32.以及

当m11=0时,有b0=0,易知λ=0是方程(4)的一个零根,同时:

这就意味着当a1+b1<0时,λ=0是方程(4)的一个单零特征根.
接下来把λ=iω代入到方程(4)中,分离实虚部得:

将其化简得:

令ω2=l及h(l)=l2+(a22-2a1)l+(a21-b21),所以方程(4)就至少有一个正根ω0,同时由(5)式有:

记

所以(±iω0,τ0)是方程(4)的解.
记λ(τ)=α(τ)+iω(τ)是方程(4)满足α(τ0)=0,ω(τ0)=ω0的根,定义,于是当p=p*,τ=τ0时,有下述定理成立:
定理:若a1+b1<0,当p=p*,τ=τ0,时,方程(4)的特征根除了λ=0和λ=±iω0外,其余特征根都具有负实部.
证明 根据前面的分析,当p=p*,τ=τ0时,方程(4)有一个单零特征根和一对纯虚根.假设方程(4)还具有正实部的特征根,并记为λ=α0+iβ0,那么λ=α(τ)+iβ(τ)就是方程(4)在p=p*和τ=τ0时满足α(τ0)=α0>0和β(τ0)=β0的根.于是存在一个正数0<ζ<τ0,使得当τ∈(τ0-ζ,τ0)时,有α(τ)>0 成立.
当τ=0时,D(λ,0)=λ3+a2λ2+(a1+b1)λ=0有一个单零特征根和一对纯虚根.进一步:

矛盾.证毕.
2 数值模拟
选择以下参数值进行数值模拟:a=0.312,b=0.171,c=0.360,r=0.312,s=0.123,d=0.396,l=0.571,通过简单的……
