罗超，张晓丹

北京科技大学数理学院，北京

收稿日期：2018年5月1日；录用日期：2018年5月17日；发布日期：2018年5月25日

本文研究了具有饱和发生率和logistic增长的随机病毒感染模型，当病毒感染细胞基本再生数 R 0 < 1 时，确定未感染细胞的有界性和病毒感染细胞的灭绝性；当病毒感染细胞基本在生数 R 0 > 1 时，构建适当的Lyapunov函数，确定平稳分布的充分条件。它显示了病毒感染在机体内持久性存在。

关键词 :饱和发生率，Logistic增长，平稳分布

Copyright © 2018 by authors and beplay安卓登录

This work is licensed under the Creative Commons Attribution International License (CC BY).

http://creativecommons.org/licenses/by/4.0/

自Nowak和Bangham提出病毒感染模型以来，随着科学的进步，病毒感染模型的机理逐渐被熟悉，模型不断被丰富和完善，为医学发展和疾病预防和控制提供重要的参考价值。许多文献中所提出病毒感染模型均在理想状态和确定参数来研究系统的稳定性分析 [ 1 ] - [ 8 ] 。在文献 [ 3 ] 中，Bonhoeffer等提出的病毒感染模型：

{ d x d t = s − d x − β x y d y d t = β x y − μ y

然而在一些生态数学模型中，无论环境波动如何，确定性系统都有一定的局限性，对于数值拟合并非理想 [ 9 ] 。自然界任何事物不可能处在绝对稳定的状态，必然受到环境的波动影响，因此一些学者将随机波动因素考虑在传染病及病毒模型中 [ 9 ] - [ 16 ] ，本文根据 [ 3 ] 中的模型，考虑饱和感染发生率，未感染细胞的logistic增长，及未感染细胞的死亡率和免疫发生率的波动影响，提出如下病毒感染模型

{ d x d t = s + r x ( 1 − x K ) − d x − k x y x + y + σ 1 x d B 1 ( t ) d t d y d t = k x y x + y − ( a + μ ) y + σ 2 y d B 2 ( t ) d t (1.1)

其中，x是未感染细胞，y是感染细胞，s是宿主细胞以常数生产，r为宿主细胞最大生产率，K为体内细胞最大承载量，d为自然死亡率，k为饱和发生率，a为免疫细胞作用下的感染细胞死亡率，μ为感染细

胞自然死亡率。 σ 1 2 , σ 2 2 为高斯白噪音强度， B 1 ( t ) , B 2 ( t ) 为标准的维纳过程在概率空间 ( Ω , F , ( F t ) t > 0 , P ) 上。

本文令病毒感染细胞基本在生数 R 0 = k a + μ 。接下来将讨论模型(1)的灭绝性和平稳分布性质。

考虑如下n维随机微分方程

d X ( t ) = f ( X ( t ) , t ) d t + g ( X ( t ) , t ) d B (t)

对于 t > t 0 ，和任意初值 X ( t 0 ) = X 0 ∈ R n 。 B ( t ) 是在全概率空间 ( Ω , F , { F t } t > 0 , P ) 上的标准维纳过程。在 R n × [ t 0 , + ∞ ) 上，定义在 C 2 , 1 ( R n × [ t 0 , + ∞ ) ) 全体的非负函数 V ( X ) 在X处均有连续的二阶微分和在t处

有连续的一阶微分。上述方程的微分算子L被定义为：

L = ∂ ∂ t + ∑ i = 1 n f i ( X , t ) ∂ ∂ x i + 1 2 ∑ i , j = 1 n [ g T ( x i , t ) g ( x j , t ) ] i j ∂ 2 ∂ x i ∂ x j

若L是 V ∈ C 2 , 1 ( R n × [ t 0 , ∞ ) ) 的微分算子，则有：

L V ( X , t ) = V t ( X , t ) + V X ( X , t ) f ( X , t ) + 1 2 t r a c e [ g T ( X , t ) V X X ( X , t ) g ( X , t ) ] ,

其中 V t = ∂ V ∂ t ,   V X = ( ∂ V ∂ x i , ⋯ , ∂ V ∂ x n ) ,   V X X = ( ∂ 2 V ∂ x i ∂ x j ) n × n 。若 X ( t ) ∈ R n ，应用伊藤公式，有：

d V ( x ( t ) , t ) = L V ( x ( t ) , t ) d t + V X ( x ( t ) , t ) g ( x ( t ) , t ) d B ( t ) .

定理1 对于任意的初值 ( x ( 0 ) , y ( 0 ) ) ∈ R + 2 ，系统(2)有唯一正解 ( x ( t ) , y ( t ) ) 且以概率1存在。

证明类似于毛学荣等人 [ 18 ] ，故省略。

定理2 若 R 0 < 1 时，若 ( x ( t ) , y ( t ) ) 是模型(1)带有任意初值 ( x ( 0 ) , y ( 0 ) ) ∈ R + 2 的解，则有：

lim t → ∞ 1 t ∫ 0 t ( x ( s ) − K ( r − d ) 2 r ) 2 d s ≤ s K r + K 2 ( r − d ) 2 4 r 2

lim t → ∞ ln y ( t ) t ≤ ( a + μ ) ( R 0 − 1 ) − 1 2 σ 2 2

即感染细胞灭绝。

证明 (i)由于模型(2)的第一个方程为：

d x d t = s + r x ( 1 − x K ) − d x − k x y x + y + σ 1 x d B 1 ( t ) d t (1)

由(1)可知 d x d t ≤ s + r x ( 1 − x K ) − d x + σ 1 x d B 1 ( t ) d t 。

将上式不等式两边同时在区间 [ 0 , t ] 上积分，因此有：

x ( t ) − x ( 0 ) ≤ s t + ∫ 0 t [ r x ( s ) ( 1 − ( s ) K ) − d ( s ) ] d s + ∫ 0 t σ 1 x ( s ) d B 1 ( s ) (2)

(2)式两边同时除以t，有：

x ( t ) t − x ( 0 ) t ≤ s + 1 t ∫ 0 t [ r x ( s ) ( 1 − ( s ) K ) − d ( s ) ] d s + 1 t ∫ 0 t σ 1 x ( s ) d B 1 (s)

1 t ∫ 0 t ( x ( s ) − K ( r − d ) 2 r ) 2 d s ≤ s K r + K 2 ( r − d ) 2 4 r 2 + K r [ x ( t ) t − x ( 0 ) t ] + 1 t K r ∫ 0 t σ 1 x ( s ) d B 1 ( s ) (3)

由定理1可知： lim t → + ∞ x ( 0 ) t = 0 ， lim t → + ∞ x ( t ) t = 0 ；令 G 1 ( t ) = ∫ 0 t σ 1 x ( s ) d B 1 ( s ) ，则 G 1 ( t ) 是鞅。有 〈 G 1 , G 1 〉 t = ∫ 0 t σ 1 2 x ( s ) d s ≤ sup t > 0 ( σ 1 2 x ( t ) ) t ，根据鞅的强大数定律 [ 19 ] ，则有 lim t → + ∞ G 1 ( t ) t = 0 。

将(3)式两边同时取极限 t → ∞ ，则有：

lim t → + ∞ 1 t ∫ 0 t ( x ( s ) − K ( r − d ) 2 r ) 2 d s ≤ s K r + K 2 ( r − d ) 2 4 r 2

(ii) 定义函数 V 2 = ln y (t)

d V 2 = 1 y ( k x ( t ) y ( t ) x ( t ) + y ( t ) − ( a + μ ) y ( t ) ) − 1 2 σ 2 2 + σ 2 d B 2 ( t ) = k y ( t ) x ( t ) + y ( t ) − ( a + μ ) − 1 2 σ 2 2 + σ 2 d B 2 ( t ) ≤ k − ( a + μ ) − 1 2 σ 2 2 + σ 2 d B 2 ( t ) = ( a + μ ) ( k a + μ − 1 ) − 1 2 σ 2 2 + σ 2 d B 2 ( t ) : = ( a + μ ) ( R 0 − 1 ) − 1 2 σ 2 2 + σ 2 d B 2 ( t ) (4)

将上式两边同时积分在区间 [ 0 , t ] ，有：

ln y ( t ) − ln y ( 0 ) = [ ( a + μ ) ( R 0 − 1 ) − 1 2 σ 2 2 ] t + ∫ 0 t σ 2 d B 2 (s)

其中令 G 2 ( t ) = ∫ 0 t σ 2 d B 2 ( s ) ，则 G 2 ( t ) 是鞅。有 〈 G 2 , G 2 〉 t = ∫ 0 t σ 2 2 d s = σ 2 2 t ，根据鞅的强大数定律 [ 19 ] ，则有 lim t → + ∞ G 2 ( t ) t = 0 ，故有如下：

ln y ( t ) t = ( a + μ ) ( R 0 − 1 ) − 1 2 σ 2 2 + ln y ( 0 ) t + 1 t ∫ 0 t σ 2 d B 2 ( s ) (5)

当(5)式两边取极限 t → + ∞ ，又由定理1可知，初值 y ( 0 ) 有界，故 lim t → + ∞ y ( 0 ) t = 0 ，因此 lim t → + ∞ ln y ( t ) t ≤ ( a + μ ) ( R 0 − 1 ) − 1 2 σ 2 2 。

即当 R 0 < 1 ，病毒感染细胞灭绝。

在这部分，我们引入一个有用的引理。在h维欧氏空间中， X ( t ) 是齐次马尔科夫程且满足如下方程

d X ( t ) = f ( X ) d t + ∑ i j g i ( X ) d B i (t)

则扩散矩阵定义为：

Q ( X ) = ( H p q ( x ) ) ,   H p q ( x ) = g i p ( x ) g i q (x)

因此，有如下引理：

引理1 [ 17 ] [ 20 ] 如果存在一个具有光滑边界 ∂ U 的有界开域 U ∈ R + 2 ，具有如下性质：

1) 在域U及其邻域中，存在常数 J > 0 满足 ∑ p , q = 1 2 H p q ( x ) ϕ p ϕ q ≥ J ϕ 2 ， x ∈ U ， ϕ ∈ R + 2 ， ϕ = ϕ 1 2 + ϕ 2 2 ；

2) 若 X 0 ∈ R + 2 \ U ，对于任一个紧集 K ∈ R + 2 ，存在一个路径从 X 0 出发到达集合U，其平均时间 E τ 是有限的，即 E τ X 0 ∈ K < + ∞ 。

若上述条件成立，则马尔科夫过程 X ( t ) 有唯一的平稳分布 π ( ⋅ ) 。且若 f ( ⋅ ) 关于测度π为可积函数，则：

P { 1 T ∫ 0 T f ( X ( t ) ) d t = f R + 2 ( u ) π ( d u ) } = 1

引理2 当 σ 1 = σ 2 = 0 时，随机系统(2)即为确定系统，其平衡点为 ( x * , y * )

x * = ( a + μ ) ( K r − K d − r ) + K 2 ( ( a + μ ) ( r − d ) − k ) 2 + 4 ( a + μ ) 2 K r s 2 r ( a + μ ) y * = ( k − a − μ ) ( a + μ ) ( K r − K d − r ) + K 2 ( ( a + μ ) ( r − d ) − k ) 2 + 4 ( a + μ ) 2 K r s 2 r ( a + μ ) 2

定理3 若 R 0 > 1 时，以下条件成立：

1) r K > k ( y * + 5 ) 2 ( x * + y * ) ， a + μ > k ( y * + 5 ) 2 ( x * + y * ) + σ 2 2

2) 0 < M < min ( m 1 ( x * ) 2 , m 2 ( y * ) 2 )

m 1 = r K − k ( y * + 5 ) 2 ( x * + y * ) ,   m 2 = a + μ − k ( y * + 5 ) 2 ( x * + y * ) − σ 2 2

带有初值 X ( 0 ) ∈ R + 2 随机微分方程(2)，存在一个平稳分布 π ( ⋅ ) ，其遍历性为：

P { 1 T ∫ 0 T f ( X i ( t ) ) d t = ∫ R + 2 f ( x , y ) π ( d x , d y ) } = 1 ,   i = 1 , 2 .

证明 (i) 由模型(1)可知，其扩散矩阵为：

Q ( X ) = ( σ 1 2 x 2 0 0 σ 2 2 y 2 )

令 G = min ( x , y ) ∈ U ( σ 1 2 x 2 , σ 2 2 y 2 ) ，有如下

∑ p , q = 1 2 H p q ( x ) ϕ p ϕ q = σ 1 2 x 2 ϕ 1 2 + σ 2 2 y 2 ϕ 2 2 ≥ G ( ϕ 1 2 + ϕ 2 2 ) ,   ϕ 1 , ϕ 2 ∈ R + 2

因此引理1中的(1)成立。

(ii) 存在一个邻域U和一个非负函数 V ∈ C 2 ，使得 X 0 ∈ R + 2 \ U 时， L V < 0 。由于 R 0 > 1 ，模型(1.1)的确定性系统的平衡点为 ( x * , y * ) ，即为引理2中的平衡点。定义V函数如下：

V = x − x * − x * ln x x * + y − y * − y * ln y y * + ( y − y * ) 2

令 V 11 = x − x * − x * ln x x * ， V 12 = y − y * − y * ln y y * ， V 13 = ( y − y * ) 2

运用伊藤公式，则有

L V 11 = ( x − x * ) x [ s + r x ( 1 − x K ) − d x − k x y x + y ] + 1 2 σ 1 2 x * = ( x − x * ) [ s x − s x * − r K ( x − x * ) + k y * x * + y * − k x y x + y ] + 1 2 σ 1 2 x * = s ( x − x * ) 2 x x * − r K ( x − x * ) 2 − k x ( x − x * ) ( y − y * ) ( x + y ) ( x * + y * ) + k y ( x − x * ) 2 ( x + y ) ( x * + y * ) + 1 2 σ 1 2 x * ≤ − r K ( x − x * ) 2 + k x * + y * [ 1 2 ( x − x * ) 2 + ( y − y * ) 2 ] + k ( x − x * ) 2 x * + y * + 1 2 σ 1 2 x * ≤ − ( r K − 3 k 2 ( x * + y * ) ) ( x − x * ) 2 + k ( y − y * ) 2 x * + y * + 1 2 σ 1 2 x * (6)

同样，我们有：

L V 12 = ( y − y * ) [ k x y x + y − k y * x * + y * ] + 1 2 σ 2 2 y * = k x ( y − y * ) 2 ( x + y ) ( x * + y * ) + k y ( x − x * ) ( y − y * ) ( x + y ) ( x * + y * ) + 1 2 σ 2 2 y * ≤ k ( y − y * ) 2 x * + y * + k ( x * + y * ) [ 1 2 ( x − x * ) 2 + 1 2 ( y − y * ) 2 ] + 1 2 σ 2 2 y * = 3 k 2 ( x * + y * ) ( y − y * ) 2 + k ( x − x * ) 2 2 ( x * + y * ) + 1 2 σ 2 2 y * (7)

L V 13 = ( y − y * ) [ k x y x + y − ( a + μ ) y ] + 1 2 σ 2 2 y 2 = ( y − y * ) [ k x y x + y − ( a + μ ) y − k x * y * x * + y * + ( a + μ ) y * ] + 1 2 σ 2 2 y 2 = k x y * ( y − y * ) 2 ( x + y ) ( x * + y * ) + k y y * ( x − x * ) ( y − y * ) ( x + y ) ( x * + y * ) − ( a + μ ) ( y − y * ) 2 + 1 2 σ 2 2 y 2

由于 ( c + d ) 2 ≤ 2 ( c 2 + d 2 ) 成立( c , d 均非负)，则：

L V 13 ≤ k y * ( y − y * ) 2 x * + y * + k y * x * + y * [ 1 2 ( y − y * ) 2 + 1 2 ( x − x * ) 2 ]     − ( a + μ ) ( y − y * ) 2 + σ 2 2 ( y − y * ) 2 + σ 2 2 ( y * ) 2 = − ( a + μ − 3 k y * 2 ( x * + y * ) − σ 2 2 ) ( y − y * ) 2 + k y * 2 ( x * + y * ) ( x − x * ) 2 + σ 2 2 ( y * ) 2 (8)

综上：

L V ≤ − ( r K − k ( y * + 5 ) 2 ( x * + y * ) ) ( x − x * ) 2 − ( a + μ − k ( y * + 5 ) 2 ( x * + y * ) − σ 2 2 ) ( y − y * ) 2     + 1 2 σ 1 2 x * + σ 2 2 ( y * ) 2 + 1 2 σ 2 2 y * : − m 1 ( x − x * ) 2 − m 2 ( y − y * ) 2 + M (9)

若 0 < M < min ( m 1 ( x * ) 2 m 2 ( y * ) 2 ) 成立，则椭圆 m 1 ( x − x * ) 2 + m 2 ( y − y * ) 2 = M 全属于 R + 2 。选取椭圆任意一个邻域U使得U的闭包 U ¯ ∈ R + 2 。则对于任意的 X ∈ R + 2 \ U ，有 L V < 0 。因此模型(1.1)有平稳分布。

本文首先考虑具有饱和发生率和logistic增长的随机病毒感染模型。首先当R₀< 1，我们分析了感染细胞的灭绝及宿主细胞的生存。然后，我们构造有效的Lyapunov函数，确立满足平稳分布的充分条件，理论证明感染能在机体内持续存在。

罗 超,张晓丹. 一类带有饱和发生率和Logistic增长的随机病毒感染模型的灭绝性及平稳分布 Extinction and Stationary Distribution of a Classic Stochastic Viral Infection Model with Saturation Rate and Logistic Growth[J]. 应用数学进展, 2018, 07(05): 609-616. https://doi.org/10.12677/AAM.2018.75072