1. 引言
近年来,全世界持续感染丙型肝炎病毒(HCV),并有罹患慢性肝病,肝硬化和肝细胞癌风险的人口数量逐步增加。由此,丙型肝炎病毒感染已成为一个严重的全球公共卫生问题,备受关注 [1] - [8]。
最早的丙型肝炎病毒感染模型(HCV模型)是由诺伊曼等人提出 [1],具体模型如下
(1.1)
其中
分别为t时刻的未感染肝细胞数量,感染肝细胞量和病毒载量,且
表示细胞总数。S表示未感染的肝细胞可能通过移行或分化成肝细胞前体细胞,以S的构成率增加数量;a表示感染细胞的最大增值速率;b表示未感染细胞的死亡率;
表示使用抗病毒药物治疗所降低的感染率;q表示通过非细胞溶解过程自发治愈感染的肝细胞的治愈率;
表示未感染细胞的最大增值速率;
表示肝细胞游离病毒感染率;
表示感染病毒细胞死亡率;
表示病毒生产率;p表示细胞产生游离病毒的速率;c表示清除游离病毒的速率。
本文考虑以下三个具有混合切换的随机丙型肝炎病毒感染模型的动力学性质。
模型一:与状态无关的马尔科夫状态切换HCV模型
(1.2)
其初始值
,对于实数
,有
。因此,
是一个具有状态空间
的连续时间马尔科夫链,且其具有如下转移概率
(1.3)
其中
,且Q-矩阵
;
。
模型二:仍然是模型(1.2),但此时
是一个具有状态空间
的跳跃过程,对任意
和
,其过渡内核如下所示
(1.4)
无论何时,都有
,其中
。
模型三:考虑到依赖于状态的随机环境和随机扰动的影响,随机HCV模型如下
(1.5)
2. 具有独立于状态的随机HCV模型的灭绝性与持久性
在随机HCV模型(1.2)和(1.3)中,连续时间马尔科夫链
的转移概率与状态无关。令
,
,其他变量
的定义与它类似。假设
·
对所有的
,
是局部利普希茨连续的,且存在常数
,使得
;
·
连续时间马尔科夫链
是不可约和正递归的,且其具有不变概率测度
。
下面的引理表明,随机HCV模型(1.2)和(1.3)的唯一解在正象限中,并揭示了肝细胞的总数量有一个上界。
引理2.1 假设
成立,令
。则随机丙型肝炎病毒感染模型(1.2)和(1.3) 有唯一强解
,其解的初始值
。此外,仍有
(2.1)
Proof. 首先由随机模型(1.2)可以得到
(2.2)
又因为
,故有
(2.3)
令
,则有
特别地
对
应用伊藤链式法,有
从而,获得
(2.4)
引理证明完毕。
推论2.2 假设
成立,则随机丙型肝炎病毒感染系统(1.2)的解
具有不变概率测度。
证明:首先注意到
是一个费勒过程。由(2.1)式,可得
(2.5)
对任意
,令
,
是起始点为
的
的过渡内核。对任意
和
,定义概率测度
则对任意的
,由切比雪夫不等式和(2.5)可知,存在充分大的
,使得
因为
是
的紧子集,所以,
是紧的。由Krylov-Bogoliubov’s定理 [9] [10] 可知,
具有不变概率测度。推论证明完毕。
定理2.3 假设
和
成立,并进一步假设存在具有
的
:
和
使得
(2.6)
和
(2.7)
则有
(2.8)
和
(2.9)
证明:设
,由(1.2)和(2.6),可得
(2.10)
从而,获得
(2.11)
因此,利用马氏链的强遍历性定理,除了有
之外,还获得
因此,有
,a.s.
下证
,a.s.,注意到
对
应用伊藤链式法,有
(2.12)
又因为
,a.s.,所以对任意
,存在具有
的
和
,使得
这就意味着
结合(2.12)式可得:只要
就有
(2.13)
由此证得
,a.s.
根据(2.8)有
(2.14)
由(2.2)可以得到下式
之间,通过BDG不等式和(4.1)式可知,存在常数
,使得
(4.11)
对任意常数
和每个实数
,令
由切比雪夫不等式及(4.11)式,可得
又因为级数
收敛,于是利用波莱尔—康特立引理,可以得到
由C的任意性可知(4.9)式成立。
下证(4.10)式成立 令
显然,
,另一方面,(4.11)式意味着
(4.12)
因此,由( [15], Theorem 5)知,(4.10)式成立。
类似地,可以得到以下结果
(4.13)
定理4.1 假设
和
及引理3.1中的条件成立,并进一步假设存在具有
的
和
,使得
。
(a) 若
为非减的,且有
(4.14)
则有
(4.15)
和
(4.16)
(b) 若
为非增的,且有
(4.17)
则有
和
(4.18)
证明:为了使证明更加简化,记
由伊藤公式及
,可得
(4.19)
由引理4.2知
(4.20)
令
为
的二次变差,则有
由此,根据连续鞅的强大数定理,可得
(4.21)
结合(4.20)和(4.21)及
非降(或非增)的性质,并利用引理3.1和连续时间马尔科夫链的强遍历性定理,即得
这意味着对任意
充分小,存在
和
,使得
(4.22)
令
,且对任意
,
由(4.21)得
这就意味着必定会存在
和
,使得
(4.23)
而后,利用随机积分的迭代对数定律 [16] [17],可以得到
因此,必定存在
和
,使得
(4.24)
由常数变易公式,由(4.23)和(4.24)式可知,对任意和
,
(4.25)
显然,经观察可知,必定会存在实数
和
,使得
(4.26)
将(4.26)式代入(4.25)式,并由
的任意性,可得(4.15)式成立。
结合(4.15),根据(4.2)和(4.13)式可得
(4.27)
结合引理3.1即得(4.16)和(4.18)式。
基金项目
国家自然科学基金(41665006、11561009),广西自然科学基金(2018GXNSFAA380240)。