1. 引言
蚊子是疟疾和登革热等蚊媒疾病的主要传播者之一,在全球范围内的种类有3600多种。其中,以吸食血液作为食物的雌蚊是登革热、疟疾、黄热病、丝虫病、日本脑炎和寨卡病毒等其他病原体的中间寄主,且伊蚊是虫媒病毒疾病的主要载体之一 [1]。在热带和亚热带地区,蚊媒疾病猖獗,给当地公共卫生的防疫工作带来了极大的影响。目前,已经很多学者对蚊虫传疾病的传播方式和有效控制方法作大量的研究 [2] [3] [4] [5]。其中,有些学者考虑使用杀虫剂来消灭蚊子,但可能受到对环境影响或抗药性演变的关切的限制。因此,也有部分学者考虑不使用杀虫剂消灭野生蚊子的方法。例如,Alphey [6] [7] 等人研究发现不育昆虫技术是减少或消灭野生蚊子的有效方法。在此研究基础上,McLean [8] 等人进一步考虑了一种转基因策略,即释放携带显性致死因子的昆虫,携带显性致死因子的雄蚊同野外雌蚊交配产生的杂合后代在预定条件下死亡,达到减少或消灭野生蚊子的效果。此后,Cai [9] 等人提出蚊子种群动力模型,该模型假设野生蚊子和不育蚊子在互不干涉的情况下都遵循Logistic增长,不育蚊子的出生率就是它的释放率。当不育蚊子被释放到环境中并且与野生蚊子发生相互作用的情况下,此模型系统可由以下微分方程描述:
(1.1)
其中w表示野生蚊子种群数量,g表示不育蚊子种群数量,
表示野生蚊子和不育蚊子的种群数量,
表示每个个体的单位时间交配次数,
表示野生蚊子和不育蚊子与种群密度无关的死亡率,
表示野生蚊子和不育蚊子与种群密度相关的死亡率,
表示每次交配产生的野生后代数量,
表示不育蚊子的释放率。
假设蚊子种群模型存在Allee效应 [10],为简化符号,记
为a,交配率
,释放率
。其中
表示最大的交配率,
是一个线性饱和函数且当野生蚊子数量充分大时趋向于常数b,
表示不育蚊子释放率系数。然而,现实世界中不可避免遭受环境噪声的影响,比如湿度、温度、光照等因素。因此,本文考虑环境噪声因素对蚊子种群的影响,即确定性蚊子种群模型相应的随机模型如下:
(1.2)
其中
,
是两个独立的实值布朗运动,
表示白噪声。由于布朗运动是对环境中连续扰动随机事件的描述,我们进一步考虑系统中存在无法用连续扰动描述的突发事件,并使用有限状态空间的马尔科夫链来模拟环境中离散事件的扰动。假设系数
的强度取决于
,方程(1.2)转换为如下随机模型:
(1.3)
其中,
,
是右连续有限状态的马尔科夫链,是独立的布朗运动,其状态空间为
。其生成元设为
,
如下确定:
(1.4)
2. 预备知识
在本文中,使用小写字母
、
、
分别表示
、
、
的初始值,为方便起见,我们假设
其中s = 1,2。记
,
,
和
,系统(1.3)与(1.4)的随机过程
的相关算子:
(2.1)
定义
表示A的转置,特别地
,
和
是
相对于
的梯度和宽度,
和
分别是式(1.3)的漂移系数和扩散系数,定义如下:
其中
表示包含项a和b的对角矩阵。特别地,
的特殊结构意味着
,
和
可以相互表示。
假设2.1 假设
1)
和
是独立于马尔科夫链
的实值标准布朗运动。
2) 对任意
,
是非负的。
3)
满足局部利普希茨条件;
意味着对于
,
;
,
。此外,对每个
,
在
时是一致连续的,即
。
4) 马尔科夫链或它的生成元
是不可约的,即对任意的
,存在
使
。
假设2.2 用
表示
的不变测度,假设下列条件之一成立:
1) 对于每个
,任意的
,
在w上是非减的,在g上是非增的;
2)
(2.2)
3)
(2.3)
3. 存在性与不变测度
本节主要证明随机蚊子种群模型的正解存在性与不变测度。
定理3.1 对任意
,具有马尔可夫链的随机蚊子种群模型(1.3)存在初始值为
的唯一全局正解,且随机过程
是一个马尔科夫费勒过程。此外,如果
,则
和
,如果
,则
。
证明:容易知道,在
中,具有马尔可夫链的随机蚊子种群模型(1.3)的系数满足局部利普西兹条件,故系统(1.3)在时间
内具有唯一连续解,其中
,
,且解是一个强马尔科夫过程 [11]。定义
故
。
考虑
,利用伊藤公式,获得
因此,
,这意味着
故有
或
,从而,有
。因此,具有马尔可夫链的随机蚊子种群模型(1.3)存在唯一的全局连续解。
下面,关注正解的情况。首先,假设
,对任意
,我们定义以下截断函数:
;
;
。特别地,让
替代
。
成为系统(1.3)与(1.4)的解。注意
并且有
。考虑
其中
,
。令
。
由于
满足局部利普西兹条件,
,所以
。类似地,
。根据伊藤公式,有
其中算子
定义作
,
替换为
。再次运用伊藤公式,获得
因此,由
的定义,如果
,则
从而,获得
且对任意
,
,即
对任意的
和
,存在
,使得
成立,从而,有
结合
,故有
(3.1)
如果
,则选择
,证明
。显然
。
考虑初始值
和
的情况,令
满足充分小,即
(3.2)
对任意
,满足
。令
由于随机过程
具有连续性,故
。根据常数变易法公式 [12],将
定义成如下形式:
(3.3)
其中
,
。
根据(3.2),当
时,有
由上式与(3.3)相结合,获得
由于过程
具有强马尔科夫性,结合(3.1)得
定理得证。
定理3.2 如果满足假设2.1和2.2条件,则下列结论成立
如果
,则
几乎处处成立。且
分布弱收敛于唯一不变概率测度
。
如果
,则
(3.4)
此外,如果对于
,
和
存在且正连续,则
的下界是一个正常数。
如果
,则在空间
中存在一个不变概率测度
与
。此外,对于
,有
。并且有:
i) 对于
的情况,对任意
(3.5)
ii) 对于
的情况,此时存在一个
,有
其中
是
的转移概率。
证明:令
来检查边界,让
成为(1.3)在
时的解,即
根据比较定理,当
时,
,有
(3.6)
对于充分小的
,
是方程(3.6)的相关算子。由于它的非退化性,方程(3.6)的扩散系数是正常返的,且系统(3.6)存在唯一不变测度
。此外,
(3.7)
其中
是测度的总变差范数,
是
的转移概率。由于
,当
时,式(3.7)仍然成立。
当
很小时,
可以用
近似表示。由伊藤公式和随机过程
的遍历性,当
很小时,它的长期增长率
可以由临界值近似得到
(3.8)
所以,符号
决定了
是否收敛于0。
1)
的情况。
运用反证法进行论证。假设
在
上有一个不变概率测度
。根据遍历性,我们进一步得到对于每个初始值
,有
(3.9)
对任意的可测函数h来说是
可积的。由于转移概率密度是连续且为正的,不变测度
是唯一的,等价于
,其中
是在
上的勒贝格测度。因此,如果
,那么当
时,对每个
有
,这意味着对于任意初始值
,式(3.9)成立。根据比较定理,当
时,
。因此,由式(1.3),(3.6)和
,
是非负的,以及在
上w是非减的,有
(3.10)
此外,由于
,
的下界是一个正常数,
,因此存在一个
,有
(3.11)
令
和
通过伊藤公式,得到
其中
。
从而,得到
(3.12)
此外,还有
是在T中一致有界的. (3.13)
根据伊藤公式,得到
结合式(3.12),可以知道,对于
,有
(3.14)
根据式(3.6)和(3.12),获得
因此,有
(3.15)
根据(1.3),可知道
此外,根据式(3.13)的一致有界性,
的线性增长有界性,式(3.10)的几乎肯定收敛以及控制收敛定理,得到
(3.16)
根据(1.3)和(3.13),容易知道,极限
存在,且
(3.17)
根据法图引理,有
根据(3.15)和(3.16),有
(3.18)
令
为一个充分大的常数并且满足
,得到
(3.19)
根据(3.18)和(3.19),得出
(3.20)
因为
是连续和非负的,所以存在一个
,有
因此,式(3.20)可以表示为
(3.21)
类似地,存在一个
,对任意的
有
因此,结合
在g上是非增的,我们知道,存在
,
,有
(3.22)
由式(30.21)和(3.22),得到
因此,获得
(3.23)
根据式(1.3),伊藤公式,遍历性以及(3.23),得到
综合上述,得
这个结果与随机过程在空间
上具有不变概率测度的假设相矛盾。此外,
是随机过程
在
中的唯一不变测度,
是度量为0的狄拉克测度。考虑经验测度
回顾式(3.13),对每一个
,函数族
是紧集。由 [13],
在
时的任何弱极限都是
的不变概率测度。由于
是唯一不变概率测度,且
在T中有界,故
,
。
2)
的情况。
令
成为以下方程的解
其中
是(3.6)的解。根据比较定理,给定
,则
几乎处处成立。再次利用比较定理,获得
根据伊藤公式,
的遍历性,
再由假设2.2的条件(2),类似
情况的证明,得到
其中
是常数,
。根据( [14],定理5.1和5.2]的证明,获得
分布弱收敛于唯一不变概率测度
的结果。
3)
的情况。
这是转移概率的总变差收敛到一个不变测度的证明。由 [15]
,得到
的持久性和
在
上的不变概率测度的存在性。由于扩散项的非退化性,对任意的
,随机过程
具有不可约性和强费勒性质,则随机过程
在
上的转移概率到其不变概率测度的总变差具有收敛性 [11]。
现在,加以考虑收敛速度。假设2.2的条件(3)成立,回顾(3.6),存在一个
,有
(3.24)
因此,由于
在
的一致连续性,存在一个
,有
(3.25)
其中
由于
,根据佛雷德霍姆替代,存在
,有
令
满足充分小,有
定义
,获得
(3.26)
其中
,
,
。由于
则,有
(3.27)
另一方面,由可知式(3.8),存在一个
,满足
(3.28)
因此,
(3.29)
令
和
,根据( [13],命题4.1),存在一些
和
,满足
(3.30)
利用定义的
,由式(3.28)和(3.30),得
(3.31)
其中
,
,
,
。
由于
,应用伊藤公式,由(3.28),得
(3.32)
当g足够小且
时,应用(3.31)和(3.32),通过
估计
。
同理,使用式(3.31),(3.32)和(3.27),根据( [13],定理4.1),存在
和
,有
(3.33)
此外,如果
,对
,
,有
,
(3.34)
利用式(3.33)和(3.34)以及( [14],命题22),获得
(3.35)
其中
,
是
上的紧集,
是正常数,和
由于
时,
的漂移项是非负的,完备集
在
上是紧集 [9]。
根据(3.35),参考文献 [16],有
对任意的
,获得
如果
,由式(3.33)和(3.34),得
其中
和
。
对
,有
即随机过程的转移概率
收敛到不变测度
,并获得指数收敛速度,证明完成。
基金项目
国家自然科学基金(41665006,11561009),广西自然科学基金(2018GXNSFAA380240)。
参考文献