1. 引言
捕食是指某物种消耗另一物种活体的全部或部分身体,从而从中直接获得营养以维持自己生命的现象。前者我们称为捕食者,后者称为食饵。在种群动力学中,Lotka-Volterra模型是最为经典和重要的两个种群之间相互作用的动力学模型,该模型分别由意大利数学家Volterra于1923年用来解释鱼群变化规律和美国种群学家Lotka于1921年在研究化学反应时提出。而捕食者–食饵模型是Lotka-Volterra模型的一类特殊情况,模型如下:
该捕食系统是指一个种群以另外一个种群为食饵,其中
分别表示种群x和y的内禀增长率,
反映了以下事实,即在环境容纳量一定的条件下,x的增大将会使每一个体平均的生活条件降低,从而影响种群的相对增长率,因此该项被称为密度制约项,
同理。而
则表示在t时刻,y个捕食者在单位时间内吃掉的食饵总量。其中
,表示种群x为食饵,种群y为捕食者,举一个典型的例子就是草原上的猎豹和斑马之间的关系。
传染病在捕食者-食饵系统中的传播,是种群生态学与传染病动力学的一种结合,是目前生物数学研究的热点问题之一,在害虫防治以及保护生态物种方面具有极其重要的实际意义。本文以疾病在受到密度制约的健全食饵中进行传播且捕食者以感染的食饵的捕食–食饵模型 [1] 为基础进行了进一步的讨论。
(1.1)
其中,
表示健全食饵种群的数量,
表示感染疾病的食饵种群的数量,
表示捕食者的数量;
,表示食饵的总数;
,表示内禀增长率,其中b表示出生率,d表示死亡率;
表示修正后食饵的出生函数,
;
表示修正后食饵死亡函数,其中
;K表示食饵的环境容量;b表示食饵感染率;m表示捕食者捕食食饵的捕食率;e表示捕食者的死亡率。
2. 模型改进
通过对相关文章的阅读与分析,发现系统(1.1)中存在一定的局限性,其中的捕食者的捕食目标只是感染的食饵。考虑到捕食者在捕食过程中,并不会区分食饵是否感染疾病,所以改进模型如下:
(2.1)
其中,
为功能反应函数,表示被捕食者捕获的食饵的数量。与模型(1.1)类似,当
时,说明食饵种群的出生率与食饵种群密度无关,但食饵种群的死亡率与食饵种群密度有关;当
时,说明食饵种群的出生率与食饵种群密度有关,但食饵种群的死亡率与食饵种群密度无关;当
时,说明食饵种群的出生率与死亡率都与食饵种群密度有关。
由系统(2.1)的第一个方程可以看出:当
且
时,有
,为了保证解的非负性,让初始值满足
,且
,
,
。
由
,可得
,即
;
,有
;同理,
。即可以得到
是方程组(2.1)的正不变集。
3. 平衡点及稳定性
3.1. 平衡点
系统(2.1)具有平衡点
,
,
,其中
,
。
当
时,系统(2.1)存在边界平衡点
。当
时,系统将退化为一般的捕食–食饵系统,本文将不做讨论。
3.2. 局部稳定性
系统(1.1)满足任意非零初始条件的解,
都是不稳定的 [1] 。由此可以推广到系统(2.1)上也有同样的结论:系统(2.1)满足任意非零初始条件的解,
都是不稳定的。
系统(2.1)在平衡点
的Jacobian矩阵为:
(3.1)
定义基本再生数
,当
且
时,平衡点
局部稳定;否则,平衡点
不稳定。
系统(2.1)在平衡点
的Jacobian矩阵为:
(3.2)
当
且
时,由Routh-Hurwitz判据可知
,考察以下矩阵:
(3.3)
可知,矩阵
的特征值均具有负实部,从而平衡点
局部稳定;当
且
时,平衡点
不稳定。
3.3. 全局稳定性
定理3.1:当
且
时,平衡点
在G内全局渐近稳定。
证明:由
,可得
,又
,则
(3.4)
当
时,
;由
和式(3.4)可知
,即
(3.5)
显然,系统(2.1)在M中的最大不变子集就是M本身。由Lasalle不变集原理 [2] 可知,当
时,(2.1)在G中的解都有
。由此可知系统(2.1)的极限系统为:
(3.6)
由文献 [3] ,极限系统理论,直接可以证明该极限系统:当
且
时,平衡点
在G内全局渐近稳定。
定理3.2:当
且
时,平衡点
在G内全局渐近稳定。
证明:先考虑下面的方程组
(3.7)
作变换
,得到
(3.8)
方程组(3.8)满足初始条件的解
。
由方程组(3.8)的第一个方程可知:满足上述非零初始条件的解都有当
时,
。
当
且t充分大时,极限系统(3.7)在N中的解以
为极限。由此可知系统(3.8)的平衡点
在N中全局渐近稳定,进而有等价系统(3.7)的平衡点
在G内全局渐近稳定。
对于系统(2.1),将前两个方程相加得到一个新的方程组
(3.9)
根据比较原理可得,方程组(3.9)的解,对
,存在T,使得
时有
,从而当
且t充分大时,方程组(3.9)的第三个方程有:
,由此可知当t充分大时
,也就是说,系统(3.9)有极限系统如下:
(3.10)
与定理3.1的证明类似,系统(3.10)的平衡点
全局渐近稳定,由极限系统理论可知(3.9)的平衡点
也全局渐近稳定。
综上所述,当
且
时,平衡点
在G内全局渐近稳定。
3.4. 正平衡点及极限环
考虑系统(2.1)的正平衡点
的稳定性,则可将系统(2.1)转化为
(3.11)
系统(3.11)可以存在正平衡点
,其中
,
当,且。
由系统(2.1)的前两个方程相减得:
(3.12)
所以
(3.13)
代入,得
(3.14)
其中。
所以
(3.15)
再将方程(3.14)化为:
(3.16)
当,即时,方程(3.16)存在唯一正根;此时,系统(2.1)有唯一的正平衡点,其中。由定理3.1的证明过程可知,
(3.17)
那么可以与(3.11)的证明方法类似,根据文献 [3] ,文献 [4] 直接得出时,系统(3.11)在G内至少存在一个包含正平衡点的稳定极限环。
4. 结论
本文首先提出了一类疾病在食饵中进行传播的无差别捕食的捕食–食饵系统。其次,讨论了边界平衡点的稳定性:当且时,平衡点在G内全局渐近稳定;当且时,平衡点在G内全局渐近稳定。最后,讨论正平衡点时,将系统转化为二维系统后,证明了极限系统在G内至少存在一个包含正平衡点的稳定极限环。由此,我们可以知道,在捕食者不区分食饵是否感染而进行捕食的时候,当且时,最后整个系统会只剩下健康未染病的食饵;当且时,整个系统会只剩下食饵,而捕食者则灭绝。在实际生活中可以通过向系统中投放感染的食饵,从而对捕食者和食饵的数量进行控制,比如在蛇–鼠这一系统中,用来灭蛇,或者控制蛇类与鼠的数量。