1. 引言
2011年,Ajraldi等人 [1] 和Braza [2] 提出了具有群体效应的捕食者–食饵模型,即用带有食饵平方根的功能反应函数去刻画具有这种群体效应的捕食者–食饵模型。许多学者在这类具有群体效应的捕食–食饵系统的研究方面取得了很大的进展 [3] [4] [5]。
在自然界中,食饵种群常常建立避难所为了达到躲避被捕食者捕食的目的。为了更好的符合现实生态环境,在系统中增加对食饵避难所的考虑是必要的。在文献 [6] 中作者建立了一类具有食饵保护区的捕食者–食饵模型,模型如下:
(1)
其中,
和
分别代表猎物和捕食者的密度函数,r,k,a,m分别代表食饵的内禀增长率,食饵的环境容纳量,食饵的消耗率以及食饵的庇护率,b和d分别是食饵的转化率和捕食者的自然死亡率。
本文主要研究模型(1)的Hopf分支问题。主要内容如下:首先分析了模型(1)平衡点的稳定性,接下来通过中心流行定理,规范型理论和文献 [7] 的方法分析Hopf分支的方向以及分支周期解的稳定性,最后用数值模拟来验证所得到的结论。
2. Hopf分支的存在性
通过计算可知
是系统(1)的一个正平衡点,
要使
存在,其中
。
定理1若
(H1)
则
是不稳定的,若
,则
是渐近稳定的,其中
证明:系统(1)在
处的Jacobi矩阵为
,
,
对应的特征方程为
(1.1)
其中
如果
,
成立,则(1.1)的根具有负实部,如果
,我们可以得到平衡点
是局部渐近稳定的。
下面,我们选择避难所m作为分支参数分析系统(1)在
处的Hopf分支。
引理1 [6] 若条件
成立,则模型(1)在
处产生Hopf分支。
证明:由计算可知,当
时,
,由于
,
有一对纯虚特征值
。
现在我们验证横截性条件。事实上,如果
是特征方程(1.1)的根,则
,
,因此,
,
。
这表明横截性条件满足。因此,由文献 [8] 的Poincaré-Andronov-Hopf分支定理可得,当m经过
时,系统(1)在
处产生Hopf分支。
3. Hopf分支的方向与稳定性
这一节我们将研究在条件
下,模型(1)在
处发生Hopf分支的方向及由分支产生的周期解的稳定性。
作变换
,
,则平衡点
平移到原点,为了方便起见,变换后仍用x和y表示
和
,则模型变为
(2)
对于模型(2)用Taylor展式可化简为
,
其中
以及
,
,
,
,
。
定义矩阵
,
其中,
,
,容易得到
,
那么
,
当
时,有
,
,
通过变换
,模型可重写为
,
其中
再用极坐标
,
变换,系统等价于
对于模型在
处作Taylor展式有
为了确定Hopf分支方向和分支周期解的稳定性,需要计算下面系数的符号,
其中所有的偏导数都取值于分支点
,并且
故
定义一阶Lyapunov系数为
。
由Poincaré-Andronov-Hopf分支定理 [8],
以及上面的计算可得如下定理。
定理2假设
成立,当
时,系统(1)在
处发生Hopf分支,
(i) 若
,则在
处产生Hopf分支的分支周期解是稳定的,分支方向为超临界。
(ii) 若
,则在
处产生Hopf分支的分支周期解是不稳定的,分支方向为亚临界。
Figure 1. When
, E* is asymptotically stable
图1. 当
时,E*是渐近稳定的
Figure 2. When
, E* is unstable
图2. 当
时,E*是不稳定的
4. 数值模拟
在这一节我们借助Matlab软件通过数值模拟来验证前面所得到的结果。
在系统(1)中选取参数
,
,
,
,
。经计算,
。当
时,正平衡点
是渐近稳定的,如图1。所示。当m减小到临界值
时,正平衡点
失稳并且经历Hopf分支,进一步计算可得
,由定理3可以得到Hopf分支的分支周期解是轨道渐近稳定的,方向为超临界。当
时,正平衡点
是不稳定的,且系统(1)存在一个稳定的极限环,如图2所示。