本文研究了按比率p移除易感植株为控制措施的一类Lotka-Volterra Filippov型植物传染病模型。讨论了无病平衡点和地方病平衡点的存在性和稳定性以及伪焦点的存在性。通过构造合适的非光滑的Lyapunov函数利用广义Lyapunov方法研究该模型在不同经济阈值下的全局动力学,得到三个结果。 In this paper, a Lotka-Volterra Filippov plant infectious disease model was studied with the removal of susceptible plants at a ratio of p as control measure. The existence and stability of disease-free equilibrium and endemic equilibrium and the existence of false focus were discussed. By constructing appropriate non-smooth Lyapunov functions and using the generalized Lyapunov method, the global dynamics of the model under different economic thresholds were studied, and three results were obtained.
本文研究了按比率p移除易感植株为控制措施的一类Lotka-Volterra Filippov型植物传染病模型。讨论了无病平衡点和地方病平衡点的存在性和稳定性以及伪焦点的存在性。通过构造合适的非光滑的Lyapunov函数利用广义Lyapunov方法研究该模型在不同经济阈值下的全局动力学,得到三个结果。
植物传染病模型,伪焦点,广义Lyapunov方法,全局动力学
Xuan Leng, Huimin Chen
School of Mathematics and Statistics, Changsha University of Science and Technology, Changsha Hunan
Received: Mar. 22nd, 2021; accepted: Apr. 11th, 2021; published: Apr. 27th, 2021
In this paper, a Lotka-Volterra Filippov plant infectious disease model was studied with the removal of susceptible plants at a ratio of p as control measure. The existence and stability of disease-free equilibrium and endemic equilibrium and the existence of false focus were discussed. By constructing appropriate non-smooth Lyapunov functions and using the generalized Lyapunov method, the global dynamics of the model under different economic thresholds were studied, and three results were obtained.
Keywords:Plant Disease Model, Pseudo-Focus, Generalized Lyapunov Approach, Global Dynamics
Copyright © 2021 by author(s) and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
由植物病原体所引起的植物传染病会造成巨大的作物损失,如何有效地控制植物传染病引起了科学家和各国政府的高度重视。为了有效地控制植物传染病,减少对环境的危害,了解植物传染病的传播动力学是至关重要的。大量的理论和实证研究表明:建立微分方程描述的数学模型是研究传染病传播动力学的一个非常有效的方法,通过对所建模型的分析研究能帮助我们理解影响疾病传播的潜在机理,为控制策略的指定提供重要建议和依据。例如Jeger等人 [
由于化学残留物对环境有直接的负面影响,所以不提倡使用化学方法来控制植物传染病。近年来采用多种控制措施相结合的疾病控制策略开始进入人们视线。直到1972年,疾病综合治理(IDM)的概念 [
Zhao等人 [
{ d S d t = P ( S ) − β S I − η 1 S + ψ ( I ) p S , d I d t = β S I − η 2 I − ψ ( I ) ν S ,
其中 P ( S ) = α 表示植株的常数种植比率, ψ ( I ) 是一个与经济阈值ET有关的控制函数,具体表达式为:
ψ ( I ) = { 0 , I < E T , 1 , I > E T .
模型具有五种类型的平衡,其中一种类型的平衡,用植物初始数的迭代方程来讨论全局结构。对于其他四种类型的平衡,通过构造适当的Lyapunov函数得到每种类型的有界全局焦点。Wang等人 [
基于传染病没爆发时追求经济效益最大化以及当传染病爆发时经济损失最小化的原则,本文在Wang和Zhao研究的Lotka-Volterra Filippov型植物传染病模型上,考虑一类以移除易感植株为控制策略的植物传染病模型,为了追求经济效益最大化,我们在传染病还没有发现时,可以采取按常数 A 与按比率 p 种植易感植株的方式获得最大经济效益,考虑到环境容纳量的原因,我们假设种植比率 p 小于易感植株的自然死亡率 η 1 。传染病在经济阈值为 σ 1 爆发时,为了追求经济损失最小化,将不再按比率 p 种植易感植株达到控制传染病传播和及时止损的目的。用右端不连续微分方程描述一个具有经济阈值的Lotka-Volterra Filippov型植物传染病模型,具体模型如下:
{ d S d t = A − β S I − η 1 S + φ ( I ) p S , d I d t = β S I − η 2 I . (1.1)
其中 A 表示易感植株的常数种植数量; β 表示传染病的感染率; η 1 表示易感植株的自然死亡率; η 2 表示感染植株的死亡率,包括自然死亡率和因病死亡率。 p 表示易感植株的种植比率,这里我们假设 p < η 1 。 φ ( I ) 是一个控制函数,函数值与易感植株的数量有关,具体表达为:
φ ( I ) = { 0 , I > σ 1 , 1 , I < σ 1 .
定义
Σ = { ( S , I ) ∈ ℝ + 2 | I = σ 1 }
和
G 1 = { ( S , I ) ∈ ℝ + 2 | I < σ 1 } , G 2 = { ( S , I ) ∈ ℝ + 2 | I > σ 1 }
其中 ℝ + 2 = { S > 0 , I < 0 } 。模型(1.1)在 G 1 上的向量场方程 F 1 为:
{ d S d t = A − β S I − η 1 S + p S , d I d t = β S I − η 2 I .
模型(1.1)在 G 2 上的向量场方程 F 2 为:
{ d S d t = A − β S I − η 1 S , d I d t = β S I − η 2 I .
在本小节中,将介绍本文中要用到的一些基本的理论知识,其中关于Filippov系统的定性理论参考文献 [
引理2.1 [
d d t V ( x ( t ) ) = 〈 ζ ( t ) , V ( x ( t ) ) 〉 , ∀ ζ ∈ ∂ V ( x ( t ) ) .
下面我们给出一个关于非光滑Lyapunov函数的LaSalle不变原理的简短描述,具体细节可以参考文献 [
d d t V ( x ( t ) ) = { a ∈ ℝ | ∃ ξ ∈ F ( x ) 使 得 ζ · ξ = a , ∀ ζ ∈ ∂ V ( x ( t ) ) }
引理2.2 [
Z V = { x ∈ ℝ | ∃ ξ ∈ F ( x ) 使 得 〈 ζ · ξ 〉 = 0 , ∀ ζ ∈ ∂ V ( x ) }
M 是 Z V ∩ L l 的最大弱不变集, x ( t ) 是系统满足初值条件 x ( 0 ) = x 0 ∈ L l 的解。有当 t → + ∞ 时 d i s t ( x ( t ) , M ) → 0 。特别地,当 M = x * 并且 l → + ∞ 时, x * 是全局渐近稳定的。
在本小节中,我们主要研究模型(1.1)在区域 G i , i = 1 , 2 上的动力学行为。在开始讨论本节之前,首先给出两个引理保证模型(1.1)解的正性和有界性,说明模型(1.1)的建立符合生物学意义。
引理3.1 [
由模型(1.1)的第一个方程,得到 d S d t | S = 0 = A > 0 所以当 S 0 > 0 有 S ( t ) > 0 。接下来用反证法证明 I ( t ) > 0 。假设存在 t 0 使得 I ( t 0 ) ≤ 0 ,那么存在 t 1 > 0 使得 I ( t 1 ) = 0 并且对于所有的 t ∈ [ 0 , t 1 ) ,有 I ( t ) > 0 。由模型(1.1)的第二个方程可以得到对于所有的 t ∈ [ 0 , t 1 ) 有
d I d t = I ( β S − η 2 ) ≥ − η 2 I ,
因此当 t ∈ [ 0 , t 1 ) 时有 I ( t ) ≥ I 0 e − η 2 t 。特别地,有 I ( t 1 ) ≥ I 0 e − η 2 t 1 > 0 这和 I ( t 1 ) = 0 相矛盾。证毕。
引理3.2 [
设 ( S ( t ) , I ( t ) ) 为模型(1.1)满足初始条件 S ( 0 ) = S 0 > 0 和 I ( 0 ) = I 0 > 0 的解,且解的最大存在区间为 [ 0,T ) ,其中 T ∈ [ 0 , + ∞ ) 。令 N = S + I 。有
N ′ = A − η 1 S + φ p S − η 2 I ≤ A − η 1 S + p S − η 2 I ≤ A − ξ N ,
其中 ξ = min { η 1 − p , η 2 } 。因此
d d t ( N ( t ) e ξ t ) = e ξ t ( N ′ ( t ) + ξ N ( t ) ) ≤ A e ξ t (2.1)
同时对不等式(2.1)两边进行0到T积分
N ( T ) e ξ T − N ( 0 ) = ∫ 0 T d d t ( N ( t ) e ξ t ) d t ≤ ∫ 0 T A e ξ t d t = A ξ ( e ξ T − 1 )
其中 N ( 0 ) = S ( 0 ) + I ( 0 ) = S 0 + I 0 ,当 N ( 0 ) ≤ A ξ ,有
N ( T ) e ξ T ≤ N ( 0 ) + A ξ ( e ξ T − 1 ) ≤ A ξ e ξ T (2.2)
根据不等式(2.2)得到 N ( T ) ≤ A ξ ,这说明解 ( S ( t ) , I ( t ) ) 在 [ 0,T ) 上有界。利用延拓定理可得解 ( S ( t ) , I ( t ) ) 在 [ 0 , + ∞ ) 上有界。
这部分我们讨论 F i 在区域 G i 上的全局动力行为,其中 i = 1 , 2 。首先,通过计算得到子系统 F i 的无病平衡点 E 0 i = ( S 0 i , 0 ) , i = 1 , 2 为
E 01 = ( A η 1 − p , 0 ) , E 02 = ( A η 1 , 0 ) .
应用下一代矩阵方法,计算子系统 F i 的基本再生数 R i , i = 1 , 2 。
R 1 = A β η 2 ( η 1 − p ) , R 2 = A β η 1 η 2 .
对于子系统 F i 而言,当 R i < 1 时,不存在地方病平衡点。当 R i > 1 时,存在唯一的地方病平衡点 E i = ( S i , I i ) , i = 1 , 2 。通过计算,子系统的无病平衡点 E i 为:
E 1 = ( S 1 , I 1 ) = ( η 2 β , A β − η 2 ( η 1 − p ) β η 2 ) , E 2 = ( S 2 , I 2 ) = ( η 2 β , A β − η 1 η 2 β η 2 ) .
我们得到如下结果。
定理3.3 对于子系统 F i 而言,当 R i < 1 时,无病平衡点 E 0 i 是全局渐近稳定的;当 R i > 1 时,地方病平衡点 E i 是全局渐近稳定的,其中 i = 1 , 2 。
证明 当 R i < 1 时,考虑Lyapunov函数
V ( S , I ) = ( S − S 0 i − S 0 i ln I S 0 i ) + I , i = 1 , 2.
函数 V ( S , I ) 沿着子系统 F i 求导
d V d t = − ( η 1 − ϱ i p ) ( S − S 0 i ) 2 S + β S 0 i I − η 2 I = − ( η 1 − ϱ i p ) ( S − S 0 i ) 2 S − ( 1 − R i ) η 2 I (3.1)
其中
ϱ i = { 1 , i = 1 , 0 , i = 2.
当 R i < 1 时, d V d t ≤ 0 。 d V d t = 0 当且仅当 S = S 0 i , I = 0 。根据LaSalle不变原理,此时子系统 F i 的无病平衡点 E 0 i 是全局渐近稳定的。
当 R i > 1 时,可以将系统(1.1)的平衡点 E i 表示成
E i = ( S i , I i ) = ( η 2 β , A β − η 2 ( η 1 − ϱ i p ) β η 2 ) , i = 1 , 2.
我们将子系统 F i 用如下形式表示:
{ d S d t = − ( η 1 − ϱ i p + β I ) ( S − S i ) − S i ( I − I i ) , d I d t = β ( S − S i ) I . , i = 1 , 2.
考虑Lyapunov函数
V ( S , I ) = 1 2 ( S − S i ) 2 + S i ( I − I i − I i ln I I i ) , i = 1 , 2.
函数 V ( S , I ) 沿着子系统 F i 求导
d V d t = − ( η 1 − ϱ i p + β I ) ( S − S i ) 2 − β S i ( S − S i ) ( I − I i ) + β S i ( S − S i ) ( I − I i ) = − ( η 1 − ϱ i p + β I ) ( S − S i ) 2 ≤ 0 ,
d V d t = 0 当且仅当 S = S i , I = I i 。根据LaSalle不变原理,此时子系统 F i 的地方病平衡点 E i 是全局渐近稳定的。
这一部分主要我们分析模型(1.1)在不连续边界 σ 1 上的动力学行为,其中不连续边界 σ 1 可以表示为 H ( S , I ) = I − σ 1 。包括滑模区域、在滑模区域上所定义的滑模方程以及是否存在伪平衡点。
根据滑模区域的定义,当 ( ∇ H , F 1 ) ≥ 0 并且 ( ∇ H , F 2 ) ≤ 0 时 σ 1 上存在滑模区域。解不等式 β S I − η 2 I ≥ 0 并且 β S I − η 2 I ≤ 0 ,得 ∑ = { ( S , I ) ∈ σ 1 | S = η 2 β } 是一个单点集。 E p = ( S p , σ 1 ) = ( η 2 β , σ 1 ) 是系统的一个伪焦点。
我们利用微分包含定义不连续边界 σ 1 上的向量场方程,向量场方程如下:
F S ( S , I ) = λ F 1 ( S , I ) + ( 1 − λ ) F 2 ( S , I ) = ( A − β σ 1 S − η 1 S + λ p S β σ 1 S − η 2 σ 1 ) , λ ∈ ( 0 , 1 ) .
因为 ( 0 , 0 ) ∈ F S ( η 2 β , σ 1 ) ,得 λ = − A + β σ 1 S + η 1 S p S 。伪焦点 E p 存在的条件为 λ ∈ ( 0 , 1 ) ,解得 A β − η 1 η 2 β η 2 < σ 1 < A β − η 1 ( η 2 − p ) β η 2 ,即 I 2 < σ 1 < I 1 当时,伪焦点 E p 存在。
下面根据 I 1 , I 2 , σ 1 的大小关系,分三种情况分析系统(1.1)平衡点 E 1 , E 2 的实虚情况以及伪焦点 E p 的存在性。
1) 当 σ 1 > I 1 时,此时平衡点 E 1 是实平衡点,平衡点 E 2 是虚平衡点,用 E 1 R , E 1 V 表示,伪焦点 E p 不存在。
2) 当 I 1 > σ 1 > I 2 时,此时平衡点 E 1 , E 2 都是虚平衡点,用 E 1 V , E 2 V 表示,伪焦点 E p 存在。
3) 当 I 2 > σ 1 时,此时 E 1 是一个虚平衡点,平衡点 E 2 是实平衡点,用 E 1 V , E 2 R 表示,伪焦点 E p 不存在。
定理4.1 当 σ 1 > I 1 时,系统(1.1)存在唯一的实平衡点 E 1 R 并且 E 1 R 是全局渐近稳定的。
证明 设 ( S ( t ) , I ( t ) ) 是系统(1.1)满足初值条件 S ( 0 ) = S 0 > 0 , I ( 0 ) = I 0 > 0 的解,存在可测函数 γ ( t ) ∈ [ 0 , 1 ] 使得对几乎处处的 t ∈ [ 0 , + ∞ ) 有
{ d S d t = A − β S I − η 1 S + γ p S , d I d t = β S I − η 2 I . (3.1)
其中
{ γ = 1 , I < σ 1 , γ ∈ ( 0 , 1 ) , I = σ 1 , γ = 0 , I > σ 1 .
将系统(3.1)配凑成如下形式:
{ d S d t = − ( η 1 − γ p + β I ) ( S − S 1 ) − S 1 ( I − I γ ) , d I d t = β ( S − S 1 ) I . (3.2)
其中 I γ = A η 2 − η 1 − γ p β 。
构建如下的Lyapunov函数:
V ( S , I ) = { V 1 ( S , I ) , I ≤ σ 1 , V 2 ( S , I ) + V 1 ( S , E T 1 ) − V 2 ( S , E T 1 ) , I > σ 1 .
其中
V 1 ( S , I ) = 1 2 ( S − S 1 ) 2 + S 1 ( I − I 1 − I 1 ln I I 1 ) , (3.3)
和
V 2 ( S , I ) = 1 2 ( S − S 2 ) 2 + S 2 ( I − I 2 − I 2 ln I I 2 ) . (3.4)
对于 I > σ 1 区域,根据已知条件 σ 1 > I 1 得到此时 I > σ 1 > I 2 ,所以有 ∂ V 2 ∂ I = S 2 I − I 2 I > 0 ,因此 V 2 ( S , I ) > V 2 ( S , σ 1 ) 。构建的Lyapunov函数 V ( S , I ) 是正则的,所以
∂ V ( S , I ) = { ζ ∈ R 2 | ζ = ( S − S 1 , S 1 I − I λ I ) , λ ∈ [ 0 , 1 ] } ,
其中
I λ = λ I 1 + ( 1 − λ ) I 2 = A η 2 − η 1 − λ p β .
让 λ = γ ,有 ζ = ( S − S 1 , S 1 I − I λ I ) = ( S − S 1 , S 1 I − I γ I ) ∈ ∂ V ( S , I ) 。函数 V ( S , I ) 沿着系统(3.2)求导有
d V d t = − ( η 1 − γ p + β I ) ( S − S 1 ) 2 − β S 1 ( S − S 1 ) ( I − I γ ) + β S 1 S ( S − S 1 ) ( I − I γ ) = − ( η 1 − γ p + β I ) ( S − S 1 ) 2 ≤ 0. (3.5)
并且,最大弱不变集 Z V = E 1 ,根据引理2.2,系统(1.1)的实平衡点 E 1 R 是全局渐近稳定的。证毕。
定理4.2 当 I 2 > σ 1 > I 1 时,系统(1.1)的伪焦点 E p 存在并且是全局渐近稳定的。
证明 设 ( S ( t ) , I ( t ) ) 是系统(1.1)满足初值条件 S ( 0 ) = S 0 > 0 , I ( 0 ) = I 0 > 0 的解,将系统(3.1)配凑成如下形式:
{ d S d t = − ( η 1 − γ p + β I ) ( S − S p ) − β S p ( I − I γ ) , d I d t = β ( S − S p ) I . (3.6)
考虑如下的Lyapunov函数
V ( S , I ) = { V 1 ( S , I ) + V 3 ( S , σ 1 ) − V 1 ( S , σ 1 ) , I < σ 1 , V 3 ( S , I ) , I = σ 1 , V 2 ( S , I ) + V 3 ( S , σ 1 ) − V 2 ( S , σ 1 ) , I > σ 1 .
其中
V 3 ( S , I ) = 1 2 ( S − S p ) 2 + S p ( I − σ 1 − σ 1 ln I σ 1 ) . (3.7)
在已知条件 I 1 > σ 1 > I 2 下,当 I > σ 1 时有 I > I 2 ,所以有 ∂ V 2 ∂ I = S 2 I − I 2 I > 0 ,因此 V 2 ( S , I ) > V 2 ( S , σ 1 ) 。当 I < σ 1 时有 I < I 1 ,所以有 ∂ V 1 ∂ I = S 1 I − I 1 I < 0 ,因此 ( S , I ) > V 1 ( S , σ 1 ) 。构造的Lyapunov函数 V ( S , I ) 是正则的。计算
∂ V ( S , I ) = { ζ ∈ R 2 | ζ = ( S − S p , S p I − I λ I ) , λ ∈ [ 0 , 1 ] } .
让 λ = γ 有 ζ = ( S − S p , S p I − I λ I ) = ( S − S p , S p I − I γ I ) ∈ ∂ V ( S , I ) 。函数 V ( S , I ) 沿着系统(3.6)求导
d V d t = − ( η 1 − γ p + β I ) ( S − S p ) 2 − β S p ( S − S p ) ( I − I γ ) + β S p ( S − S p ) ( I − I γ ) = − ( η 1 − γ p + β I ) ( S − S p ) 2 ≤ 0. (3.8)
并且,最大弱不变集 Z V = E p ,根据引理2.2,系统(1.1)的伪焦点 E p 是全局渐近稳定的。证毕。
定理4.3 当 I 2 > σ 1 时,系统(1.1)存在唯一的实平衡点 E 2 R 并且 E 2 R 是全局渐近稳定的。
证明 设 ( S ( t ) , I ( t ) ) 是系统(1.1)满足初值条件 S ( 0 ) = S 0 > 0 , I ( 0 ) = I 0 > 0 的解,将系统(3.1)配凑成如下形式:
{ d S d t = − ( η 1 − γ p + β I ) ( S − S 2 ) − β S 2 ( I − I γ ) , d I d t = β ( S − S 2 ) I . (3.9)
构建如下的Lyapunov函数:
V ( S , I ) = { V 1 ( S , I ) + V 2 ( S , σ 1 ) − V 1 ( S , σ 1 ) , I < σ 1 , V 2 ( S , I ) , I ≥ σ 1 .
在已知条件 I 2 > σ 1 下,当 I < σ 1 时有 I < I 1 ,所以有 ∂ V 1 ∂ I = S 1 I − I 1 I < 0 ,因此 V 1 ( S , I ) > V 1 ( S , σ 1 ) ,构造的Lyapunov函数 V ( S , I ) 是正则的。所以
∂ V ( S , I ) = { ζ ∈ R 2 | ζ = ( S − S 2 , S 2 I − I λ I ) , λ ∈ [ 0 , 1 ] } .
让 λ = γ ,有 ζ = ( S − S 2 , S 2 I − I λ I ) = ( S − S 2 , S 2 I − I γ I ) ∈ ∂ V ( S , I ) 。函数 V ( S , I ) 沿着系统(3.9)求导
V ( S , I ) d V d t = − ( η 1 − γ p + β I ) ( S − S 2 ) 2 − β S 2 ( S − S 2 ) ( I − I γ ) + β S 2 ( S − S 2 ) ( I − I γ ) = − ( η 1 − γ p + β I ) ( S − S 2 ) 2 ≤ 0. (3.10)
并且,最大弱不变集 Z V = E 2 ,根据引理2.2,系统(1.1)的实平衡点 E 2 R 是全局渐近稳定的。证毕。
本小节中,通过数值模拟检验所得结论的正确性。
令 A = 8 , β = 0.05 , η 1 = 0.3 , η 2 = 0.8 , p = 0.1 此时 I 1 = A β − η 2 ( η 1 − p ) β η 2 = 6 , I 2 = A β − η 2 η 1 β η 2 = 4 。当经济阈值 σ 1 = 8 > I 1 满足定理4.1的条件,根据定理4.1实平衡点 I 2 > σ 1 是全局渐近稳定的(数值模拟图见图1(a))。当经济阈值 I 1 > σ 1 = 5 > I 2 满足定理4.2的条件,根据定理4.2伪焦点 E p 是全局渐近稳定的(数值模拟图见图1(b))。当经济阈值 σ 1 = 3 < I 2 满足定理4.3的条件,根据定理4.3实平衡点 E 2 是全局渐近稳定的(数值模拟图见图1(c))。
图1. 系统解的数值模拟图
为追求经济效益最大化和经济损失最小化的目的,本文研究了按比率 p 移除易感植株为控制措施的一类Lotka-Volterra Filippov型植物传染病模型。讨论了无病平衡点和地方病平衡点的存在性和稳定性以及焦点的存在性。通过构造合适的非光滑的Lyapunov函数利用广义Lyapunov方法研究该模型的全局动力学,得到三个结论。当 σ 1 > I 1 时,系统存在唯一的全局渐近稳定的实平衡点 E 1 ,形成地方病。当 I 1 > σ 1 > I 2 时,系统的伪焦点 E p 存在并且是全局渐近稳的,形成地方病。当 I 2 > σ 1 时,系统存在唯一的全局渐近稳定的实平衡点 E 2 ,形成地方病。
本文的理论和数值结果表明,传染病最终不会消失,而是变成地方性疾病,这样的结果是合适现实的。按比率 p 移除易感植株最终会减少染病植株的数量,我们可以采取这一措施使染病植株的数量达到人们可接受的范围内,最大程度减少经济损失。解的最终趋势依赖与经济阈值 σ 1 ,这说明合理的经济阈值对控制传染病的传播有重要的意义。
冷 轩,陈会敏. 基于广义Lyapunov方法的植物传染病模型的全局动力学Global Dynamics for a Plant Disease Model Based on Generalized Lyapunov Approach[J]. 应用数学进展, 2021, 10(04): 1086-1095. https://doi.org/10.12677/AAM.2021.104117