1. 介绍
在生态学的研究过程中,由于捕食者和食饵是普遍存在的,因此捕食者与被捕食者之间的相互作用及其关系是一个研究的热点问题之一。第一个捕食食饵模型是由Lotka [1] 和Volterra [2] 提出的,自从著名的Lotka-Volterra模型被提出以来,捕食者与食饵之间关系研究一直是生物学家和数学家等学者们研究的重点。为了更真实地反应捕食者和食饵之间的关系,更贴合真实的生态系统,学者们在捕食–食饵模型中引入了多种功能反应,并对其做了大量的研究,得到了丰富的结论 [3] [4] [5] [6]。
然而,在我们建立模型时有一个很重要的影响因素Allee效应,也是值得考虑的。Allee效应是由Allee [7] 提出来的,并以他的名字命名,他对种群的影响很显著。Allee效应表示的是个体适应性与种群密度之间的正相关关系 [8] [9]。近些年来,将Allee效应引入到捕食–食饵模型中引起了学者们的关注,并对其做了大量研究 [3] [10] - [16],得到了一些新的结论,丰富了生物数学的研究。
但是在我们的生态系统中,生态系统的种群存在是很复杂的,有未成熟种群逐渐转化为成熟种群,这种阶段变化可以更好地反应种群的生长和相互关系,更能真实地体现生态系统的瞬息万变。此外,很多学者做了很多关于具有阶段结构的捕食–食饵系统的动力学研究 [17] - [23]。Xiao在 [17] 研究了带有时滞和Allee效应具有阶段结构的捕食–食饵模型,证明了平衡点的稳定性以及产生分支的情况。蒋达清在参考文献 [24] 中提出了具有Holling Type II的阶段结构的捕食–食饵模型,得到了模型在存在遍历平稳分布以及捕食者灭绝的结论。针对此,我们提出了具有阶段结构和Allee效应的捕食–食饵的确定性模型:
(1.1)
其中
分别表示食饵种群,未成熟捕食者种群,成熟捕食者种群在t时刻的密度,m表示捕获率,n是捕食者的食物转化率,
和e分别表示未成熟和成熟捕食者的死亡率,
表示未成熟捕食者变成成熟捕食者速度,其中
均为正常数。
然而,在现实生活中,大多数现象无法用确定性的规律来解释,而总是受到环境噪声的影响。当把环境噪声引入到模型中时,我们发现它可以更好地预测未来生态系统模型的发展。许多学者研究了带有环境噪声的随机捕食–食饵模型 [5] [6] [24] [25] [26],Zhao在文献 [26] 中考虑了具有阶段结构的随机捕食–食饵模型,给出了系统存在平稳分布的充分条件,同时给出了捕食者灭绝的条件。因此,为了更好地讨论具有阶段结构和Allee效应的捕食食饵模型,我们引入了环境噪声,模型如下:
(1.2)
其中
,
,
为标准布朗运动,
为噪声强度。本文的研究主要分为以下几个部分:第二部分,介绍与本篇文章相关的预备知识;第三部分,证明系统(1.2)存在全局唯一正解;第四部分,建立了系统(1.2)存在平稳分布的条件;第五部分,证明了系统(1.2)捕食者
和食饵x灭绝的条件;最后对全文做了总结和期望。
2. 预备知识
在本部分,我们介绍了一些关于随机微分方程的基础知识和基本定理,在我们接下来的证明过程都会用到。
表示一个带有过滤的完全概率空间,其中
满足条件(即
包含所有
-空集,它是递增和右连续的)。我们引入定义
和
考虑一个d维的随机微分方程的马尔可夫过程,
(2.1)
初始值
,
表示定义在完全概率空间
上的d维标准布朗运动。用
表示定义在
上的非负函数
,使得它们在X上连续可微两次。公式(2.1)的微分算子L由 [27] 定义
(2.2)
如果将L作用在
,则可以得到
,
其中
,
。利用Itô’s公式 [27] 可得到,如果
,则有
(2.3)
定义2.1. [28] 过渡概率函数
被称为时间齐次函数(相应的Markov过程称为时间齐次函数),如果
是独立于s,
,
,
,
表示在
中的Borel集的
代数。
令
为随机微分方程中描述的
中的正则时齐马尔可夫过程
.
的扩散矩阵为
,
。
引理2.1. [28] 马尔可夫过程
有一个唯一的遍历平稳分布
,如果存在一个有界开域
,常规边界
,具有下列性质:
(B.1):存在一个正数M使
,
,
。
(B.2):存在一个非负C2-V函数,使LV在任何
上是负的。
3. 正解的存在性和唯一性
在本章节中,我们主要研究系统(1.2)存在正解并证明其正解是唯一的,对于系统(1.2)我们可以清楚的知道
为食饵,未成熟捕食者和成熟捕食者的密度,则它们均为非负常数,且其系数满足局部Lipschitz条件,所以系统(1.2)存在局部解
,
,其中
表示爆炸时间。因此,我们只需要证明系统(1.2)存在唯一正解以及正解是全局的,即证
。为此,我们给出如下定理:
定理3.1. 对于任意初始值
,系统(1.2)在
时存在唯一的全局正解,并且解始终以概率1在
。
证明:对于一个足够大的非负数
,使得
,
以及
始终在
这个区间内,对于每个整数
,定义停时 [27]
在这篇文章中
,当
时,显然可以知道
是递增的。令
,
(a.s.)。如果我们可以证明
,a.s.则可以得到
,
a.s.对于
。因此,为了完成证明,我们只需证明
a.s.如果
,则存在常数
,
,使
因此存在一个常数
,使
,
(3.1)
定义一个C2函数
(3.2)
对上述V函数利用Itô’s公式 [27] 可得
(3.3)
令
,公式(3.3)则有
(3.4)
其中K为正常数,因此可以得到:
(3.5)
对公式(3.5)从0到
两边积分,并取期望可得
(3.6)
对于
,设
,由公式(3.1)可以得到
。对于每个
,
或
或
等于z或
。因此,可以得到
(3.7)
(3.8)
其中
表示
的指标函数,当
时,可以得到
这就导致了矛盾,故我们可以得到
a.s.证毕。
4. 遍历平稳分布的存在
因为当随机扰动加入确定性模型时,系统的稳定性会被破坏,从而模型会出现随机弱稳定的特性,从而引起了很多研究者们的关注 [29] [30] [31] [32]。在本节中,我们通过构造李雅普诺夫函数,证明了系统(1.2)存在平稳分布。
定理4.1. 如果
,
,则对于系统(1.2)的任意给定的初始值
,都存在唯一的遍历平稳分布。
证明:为了证明系统(1.2)存在遍历平稳分布,只需证明引理2.1两个条件成立即可。通过题意,显然可以得到系统(1.2)的扩散矩阵为
(4.1)
对于任意的
,
,其中
,引理2.1中的条件(B.1)得证。现在我们证明条件(B.2)。为了证明条件(B.2),我们构造李雅普诺夫函数
,
,
,
,
其中
是一个足够小的正常数,
为正常数,其中
(4.2)
(4.3)
分别对上述V函数利用Itô’s公式可得,
(4.4)
则有
则对(4.4)利用不等式可得
(4.5)
(4.6)
(4.7)
(4.8)
(4.9)
下面我们定义一个C2的函数
,则有
则根据(4.5)~(4.9)可以得到
(4.10)
为了证明引理2.1中的条件(B.2),我们只需要构造一个有界开集
,定义有界开集
为
,证明
在
成立,其中
是足够小的常数,且满足如下条件:
(4.11)
(4.12)
(4.13)
(4.14)
(4.15)
(4.16)
(4.17)
(4.18)
(4.19)
下面,为了证明
,我们把
分为6个区域,
,
,
,
,
,
显然,
,接下来证明对于任意的
有
,这等价于分别在上述六个域上证明它。
情况1:对于任意的
,我们可以得到
(4.20)
因此可以得到对于任意的
有
。
情况2:对于任意的
,我们得到
(4.21)
因此可以得到对于任意的
有
。
情况3:对于任意的
,我们得到
(4.22)
因此可以得到对于任意的
有
。
情况4:对于任意的
,我们得到
(4.23)
因此可以得到对于任意的
有
。
情况5:对于任意的
,我们得到
(4.24)
因此可以得到对于任意的
有
。
情况6:对于任意的
,我们得到
(4.25)
因此可以得到对于任意的
有
。
综上,通过公式(4.17)~(4.25),我们可以得到对于足够小的
,对所有的
,
,因此引理(3.1)中的(B.2)满足,故通过引理2.1我们可知系统(1.2)具有唯一的遍历平稳分布。
5. 灭绝
这一部分,我们讨论了捕食者和食饵均灭绝的情况。
定理5.1. 若
是系统(1.2)的解,且其初始值为
,如果满足条件
则食饵x,未成熟捕食者
和成熟捕食者
均灭绝。
证明:对于系统(1.2)的第一个式子利用Itô’s公式,可以得到
(5.1)
对公式(5.1)从0到t积分,并两边同时除t
(5.2)
对公式(5.2)取上确界,利用局部鞅的强大数定理 [27],可以知道
几乎是处处成立的,
(5.3)
即
,食饵灭绝。
因此,对于任意小的
,存在
和一个集合
使得
,对于
,
,
成立。接下来,证明捕食者
和
灭绝,我们对
利用Itô’s积分
(5.4)
对公式(5.4)两边从0到t积分,并两边同时除t,
(5.5)
对公式(5.5)取上确界,利用局部鞅的强大数定理 [27],可以知道
,
几乎是处处成立的,则有
(5.6)
即
,
捕食者
和
灭绝。
6. 结论
本文主要研究了具有Allee效应和阶段结构的捕食–食饵模型随机动力学行为。首先,我们通过构造恰当的李雅普诺夫函数,建立了系统(1.2)存在全局唯一正解,并得到了系统存在遍历平稳分布的充分条件。此外,我们得到了食饵、未成熟捕食者和成熟捕食者灭绝的条件。得到了一系列的结果:
当
,
时,系统存在唯一的遍历平稳分布;当
时,食饵和捕食者均会灭绝。
当然,在整个生态系统中,会受到很多自然因素和人类的干预,因此捕食–食饵模型还是有很多值得研究的问题。例如,我们可以考虑脉冲效应等,这些都可以在以后的研究中加以思考和解决。