考虑异质环境对带有扩散项的Holling-Tanner捕食–食饵模型正稳态解全局稳定性的影响。首先利用抛物方程的比较定理等方法证明了该模型正稳态解的存在性,进一步通过迭代过程,得出了该模型正稳态解的下界。最后通过构造一种新的李雅普诺夫函数,研究Holling-Tanner捕食-食饵模型正稳态解的稳定性:该模型在环境异质下的正稳态解是全局稳定的。 Consider the influence of heterogeneous environment on the global stability of the positive steady-state solution of the Holling-Tanner predator-prey model with diffusion term. First, the existence of the positive steady-state solution of the model is proved by the comparison theorem of the parabolic equation, and the lower bound of the positive steady-state solution of the model is obtained through the iterative process. Finally, by constructing a new Lyapunov function, the stability of the positive steady-state solution of the Holling-Tanner predator-prey model is studied: the positive steady-state solution of this model is globally stable under heterogeneous environments.
考虑异质环境对带有扩散项的Holling-Tanner捕食–食饵模型正稳态解全局稳定性的影响。首先利用抛物方程的比较定理等方法证明了该模型正稳态解的存在性,进一步通过迭代过程,得出了该模型正稳态解的下界。最后通过构造一种新的李雅普诺夫函数,研究Holling-Tanner捕食–食饵模型正稳态解的稳定性:该模型在环境异质下的正稳态解是全局稳定的。
捕食–食饵模型,异质环境,正稳态解,全局稳定性,Lyapunov函数
Wanhong Zhang
School of Science, Chang’an University, Xi’an Shaanxi
Received: Aug. 17th, 2021; accepted: Sep. 19th, 2021; published: Sep. 26th, 2021
Consider the influence of heterogeneous environment on the global stability of the positive steady-state solution of the Holling-Tanner predator-prey model with diffusion term. First, the existence of the positive steady-state solution of the model is proved by the comparison theorem of the parabolic equation, and the lower bound of the positive steady-state solution of the model is obtained through the iterative process. Finally, by constructing a new Lyapunov function, the stability of the positive steady-state solution of the Holling-Tanner predator-prey model is studied: the positive steady-state solution of this model is globally stable under heterogeneous environments.
Keywords:Predator-Prey Model, Heterogeneous Environment, Positive Steady-State Solution, Global Stability, Lyapunov Function
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/
生物数学模型是人们研究种群规模发展变化的重要手段之一 [
考虑在异质环境下带有扩散项的Holling-Tanner捕食–食饵模型:
{ u t = d 1 ( x ) Δ u + a ( x ) u − u 2 − b u v m + r u , x ∈ Ω , t > 0 , v t = d 2 ( x ) Δ v + b v − b v 2 γ u , x ∈ Ω , t > 0 , ∂ u ( x , t ) ∂ ν = ∂ v ( x , t ) ∂ ν = 0 , x ∈ ∂ Ω , t > 0 , u ( x , 0 ) = u 0 ( x ) > 0 , v ( x , 0 ) = v 0 ( x ) ≥ ( ≠ ) 0 , x ∈ Ω (1.1)
其中, u ( x , t ) 和 v ( x , t ) 分别代表捕食者和食饵的密度函数, Ω 是 R n 中具有光滑边界 ∂ Ω 的有界区域,施加在 ∂ Ω 上的无通量边界条件表示它们生活在一个封闭的环境里。 d 1 ( x ) 和 d 2 ( x ) 分别是u和v的空间依赖扩散系数函数, a ( x ) 是空间异质资源函数,其中 b , r , m , γ 都是常数。
为了证明异质环境下扩散Holling-Tanner模型正稳态解的存在性和全局稳定性,引入如下引理。
引理1 [
{ ∂ ω ∂ t − d Δ ω ≤ ω 1 + β f ( ω ) ( α − ω ) , x ∈ Ω , t ∈ ( T , ∞ ) ∂ ω ∂ n = 0 , x ∈ ∂ Ω , t ∈ [ T , ∞ )
且常数 α > 0 ,那么
lim t → ∞ sup max Ω ¯ ω ( ⋅ , t ) ≤ α (2.1)
引理2 [
φ ′ ( t ) ≤ − ψ ( t ) + h ( t ) , ψ ′ ( t ) ≤ K
(其中 K > 0 )
那么 lim t → ∞ ψ ( t ) = 0 。
系统(1.1)对应的稳态问题为:
{ − d 1 ( x ) Δ u = u ( a ( x ) − u − b v m + r u ) , x ∈ Ω − d 2 ( x ) Δ v = b v ( 1 − v γ u ) , x ∈ Ω ∂ u ( x , t ) ∂ ν = ∂ v ( x , t ) ∂ ν = 0 , x ∈ ∂ Ω
本节给出扩散Holling-Tanner模型正稳态解的存在性结果。
定理1 系统(1.1)存在唯一解 ( u ( x , t ) , v ( x , t ) ) ,并满足:对任意小的 ε 1 , ε 2 (其中 ε 1 > 0 , ε 2 ∈ ( 0 , ε 1 ] ),存在一个常数 T > 0 ,当 t > T 时,有
1) ε 2 ≤ u ( x , t ) , v ( x , t ) ≤ c ¯ + ε 1 , ∀ x ∈ Ω ¯ . (3.1)
2) 存在一个常数 G > 0 ,使得
max t > T 1 ‖ u ( ⋅ , t ) ‖ C 2 + α ( Ω ¯ ) , max t > T 1 ‖ v ( ⋅ , t ) ‖ C 2 + α ( Ω ¯ ) ≤ G . (3.2)
其中
c ¯ : = max Ω ¯ a ( x ) , c _ : = min Ω ¯ a ( x ) . (3.3)
证明 由系统(1.1)的第一个方程, u ( x , t ) 满足
u t − d 1 ( x ) Δ u ≤ u ( max x ∈ Ω ¯ a ( x ) − u )
利用引理1,得
lim sup t → ∞ max Ω ¯ u ( x , t ) ≤ max x ∈ Ω ¯ a ( x ) = c ¯
那么对任意的 ε > 0 ,存在 T > 0 ,使得 u ( x , t ) < c ¯ + ε 。故模型(1.1)的第二个方程 v t = d 2 ( x ) Δ v + b v − b v 2 γ u 中, v ( x , t ) 满足
v t − d 2 ( x ) Δ v ≤ b v ( 1 − v γ ( c ¯ + ε ) )
结合系统(1.1)的边界条件 ∂ u ∂ ν = 0 ,同理,利用引理1得出
lim sup t → ∞ max Ω ¯ v ( x , t ) ≤ c ¯ + ε
故存在 T > 0 ,使得
u ( x , t ) , v ( x , t ) ≤ c ¯ + ε 1 , x ∈ Ω ¯ , t > T .
又因 u ( x , T ) , v ( x , T ) > 0 ,我们可以选择一个很小的常数 ε 2 ,其中 ε 2 ∈ ( 0 , ε 1 ] ,使得 u ( x , T ) , v ( x , T ) > ε 2
故
ε 2 ≤ u ( x , t ) , v ( x , t ) ≤ c ¯ + ε 1 .
因此,问题(1.1)在 [ ε 2 , c ¯ + ε 1 ] × [ ε 2 , c ¯ + ε 1 ] 上存在一个正稳态解 ( u ( x , t ) , v ( x , t ) ) 。
对任意的 T > 0 ,定义 u ¯ 1 = v ¯ 1 = a ¯ + ε 1 , u _ 1 = v _ 1 = ε 2 。显然有
| a ( x ) u 1 − u 1 2 − b u 1 v 1 m + r u 1 − a ( x ) u 2 + u 2 2 + b u 2 v 2 m + r u 2 | ≤ T ( | u 1 − u 2 | + | v 1 − v 2 | ) , (4.1)
| b v 1 − b v 1 2 r u 1 − b v 2 + b v 2 2 r u 2 | ≤ T ( | u 1 − u 2 | + | v 1 − v 2 | ) .
定义如下的迭代序列
{ u ¯ i + 1 = u ¯ i + 1 T ( a ¯ u ¯ i − u ¯ i 2 − b u ¯ i v _ i m + r u ¯ i ) , u _ i + 1 = u _ i + 1 T ( a _ u _ i − u _ i 2 − b u _ i v ¯ i m + r u _ i ) , v ¯ i + 1 = v ¯ i + 1 T ( b v ¯ i − b v ¯ i 2 r u ¯ i ) , v _ i + 1 = v _ i + 1 T ( b v _ i − b v _ i 2 r u _ i ) . (4.2)
定理2 迭代序列(4.2)满足以下的单调性和收敛性:
1) 单调性:常数序列 { u ¯ i } i = 1 ∞ , { u _ i } i = 1 ∞ , { v ¯ i } i = 1 ∞ , { v _ i } i = 1 ∞ 满足
{ 0 < u _ 1 ≤ ⋯ ≤ u _ i ≤ u _ i + 1 ≤ ⋯ ≤ u ¯ i + 1 ≤ u ¯ i ≤ ⋯ ≤ u ¯ 1 , 0 < v _ 1 ≤ ⋯ ≤ v _ i ≤ v _ i + 1 ≤ ⋯ ≤ v ¯ i + 1 ≤ v ¯ i ≤ ⋯ ≤ v ¯ 1 .
2) 收敛性:定义 u ¯ ∞ = lim i → ∞ u ¯ i , u _ ∞ = lim i → ∞ u _ i , v ¯ ∞ = lim i → ∞ v ¯ i , v _ ∞ = lim i → ∞ v _ i 满足
{ a ¯ u ¯ ∞ − u ¯ ∞ 2 − b u ¯ ∞ v _ ∞ m + r u ¯ ∞ = 0 , a _ u _ ∞ − u _ ∞ 2 − b u _ ∞ v ¯ ∞ m + r u _ ∞ = 0 , b v ¯ ∞ − b v ¯ ∞ 2 r u ¯ ∞ = 0 , b v _ ∞ − b v _ ∞ 2 r u _ ∞ = 0. (4.3)
证明 (1)根据定理1 (3.1)式以及迭代序列(4.2)可得
{ u _ 1 ≤ u _ 2 ≤ u ¯ 2 ≤ u ¯ 1 , v _ 1 ≤ v _ 2 ≤ v ¯ 2 ≤ v ¯ 1 . (4.4)
{ a ¯ u ¯ 1 − u ¯ 1 2 − b u ¯ 1 v _ 1 m + r u ¯ 1 ≤ 0 , a _ u _ 1 − u _ 1 2 − b u _ ∞ v ¯ 1 m + r u _ 1 ≥ 0.
{ b v ¯ 1 − b v ¯ 1 2 r u ¯ 1 ≤ 0 b v _ 1 − b v _ 1 2 r u _ 1 ≥ 0
根据(3.3)式有 c _ ≤ c ¯ ,结合迭代序列(4.1)
利用数学归纳法
当 i = 1 时
u ¯ 2 − u _ 2 ≥ u ¯ 1 − u _ 1 + 1 T ( a ¯ u ¯ 1 − u ¯ 1 2 − b u ¯ 1 v _ 1 m + r u ¯ 1 ) − 1 T ( a _ u _ 1 − u _ 1 2 − b u _ 1 v ¯ 1 m + r u _ 1 ) ≥ 0
同理有 v ¯ 2 − v _ 2 ,因此(4.4)式是成立的。
假设 i = N 时也成立,当 i = N + 1 时有
u ¯ N + 1 − u _ N + 1 ≥ u ¯ N − u _ N + 1 T ( a ¯ u ¯ N − u ¯ N 2 − b u ¯ N v _ N m + r u ¯ N ) − 1 T ( a _ u _ N − u _ N 2 − b u _ N v ¯ N m + r u _ N ) ≥ 0
u ¯ N + 1 − u ¯ N ≤ u ¯ N − u ¯ N − 1 + 1 T ( a ¯ u ¯ N − u ¯ N 2 − b u ¯ N v _ N − 1 m + r u ¯ N ) − 1 T ( a ¯ u ¯ N − 1 − u ¯ N − 1 2 − b u ¯ N − 1 v _ N − 1 m + r u ¯ N − 1 ) ≤ 0
成立。同样可以得到:
u _ i + 1 ≥ u _ i , v _ i + 1 ≥ v _ i , v ¯ i + 1 ≥ v ¯ i + 1 , v ¯ i + 1 ≤ v ¯ i .
单调性证明完毕。
(2)在证明(1)中的公式已经表明序列 { u ¯ i } i = 1 ∞ , { u _ i } i = 1 ∞ , { v ¯ i } i = 1 ∞ , { v _ i } i = 1 ∞ 是分别收敛于某个常数的,因此收敛性也得到了证明。
本节得出了问题(1.1)正稳态解的下界,通过迭代过程,得到了问题(1.1)正稳态解的更精确估计。
对于资源分布均匀的捕食–食饵系统,可以采用上下解方法或者构造李雅普诺夫函数来说明正稳态解的全局稳定性 [
E ( t ) : = ∫ Ω ∫ u * ( x ) u ( x , t ) u * ( x ) d 1 ( x ) s − u * ( x ) s d s d x + ξ ∫ Ω ∫ v * ( x ) v ( x , t ) v * ( x ) d 2 ( x ) s − v * ( x ) s d s d x
定理3 假设在 Ω ¯ 上,对任意 0 < α < 1 , a ( x ) > 0 , d i ( x ) > 0 ,我们定义
d i ¯ = max x ∈ Ω ¯ d i ( x ) , d i _ = min x ∈ Ω ¯ d i ( x ) .
若
ξ = d 2 _ d 2 ¯ ( u _ ∞ − ε ) ( u ¯ ∞ − ε ) u _ ∞ 3 u ¯ ∞ d 1 _ d 1 ¯ v _ ∞ v ¯ ∞ 3 1 m + r ( u _ ∞ − ε )
那么方程(1.1)有一个唯一的正稳态解 ( u * , v * ) ,且正稳态解 ( u * , v * ) 是全局稳定的。其中 lim t → ∞ u ( x , t ) = u * ( x ) , lim t → ∞ v ( x , t ) = v * ( x ) 。
证明 设 ( u ( x , t ) , v ( x , t ) ) 是方程(1.1)的解,定义一个函数
E ( t ) : = ∫ Ω ∫ u * ( x ) u ( x , t ) u * ( x ) d 1 ( x ) s − u * ( x ) s d s d x + ξ ∫ Ω ∫ v * ( x ) v ( x , t ) v * ( x ) d 2 ( x ) s − v * ( x ) s d s d x
其中 ξ > 0 ,且当 t ≥ 0 时 E ( t ) ≥ 0 ,根据格林公式
d E ( t ) d t = ∫ Ω u * ( u − u * ) d 1 u u t d x + ξ ∫ Ω v * ( v − v * ) d 1 v v t d x = ∫ Ω u * ( u − u * ) d 1 u ( d 1 Δ u + f ( a ( x ) , u , v ) − u u * d 1 Δ u * − u u * f ( a ( x ) , u * , v * ) ) d x + ξ ∫ Ω v * ( v − v * ) d 2 v ( d 2 Δ v + g ( u , v ) − v v * d 2 Δ v * − v v * g ( u * , v * ) ) d x
= ∫ Ω [ u * ( u − u * ) u ( Δ u − u u * Δ u * ) + u * ( u − u * ) d 1 ( f ( a ( x ) , u , v ) u − f ( a ( x ) , u * , v * ) u * ) ] d x + ξ ∫ Ω [ v * ( v − v * ) v ( Δ v − v v * Δ v * ) + v * ( v − v * ) d 2 ( g ( u , v ) v − g ( u * , v * ) v * ) ] d x ≤ ∫ Ω u * ( u − u * ) d 1 ( f ( a ( x ) , u , v ) u − f ( a ( x ) , u * , v * ) u * ) d x + ξ ∫ Ω v * ( v − v * ) d 2 ( g ( u , v ) v − g ( u * , v * ) v * ) d x
其中 f ( a ( x ) , u , v ) : = a ( x ) u − u 2 − b u v m + r u , g ( u , v ) : = b v − b v 2 γ u 。
d E ( t ) d t ≤ ∫ Ω u * ( u − u * ) d 1 ( − u − b v m + r u + u * + b v * m + r u * ) d x + ξ ∫ Ω v * ( v − v * ) d 2 b ( − v γ u + v γ u * ) d x = ∫ Ω J 1 ( u − u * ) 2 + J 2 ( u − u * ) ( v − v * ) + J 3 ( v − v * ) 2 d 1 d 2 γ u u * ( m + r u ) ( m + r u * ) d x
其中
J 1 : = − d 2 γ u u * 2 [ ( m + r u ) ( m + r u * ) − b r v * ] ,
J 2 : = − b d 2 γ u u * 2 ( m + r u * ) + d 1 ξ b v * 2 ( m + r u ) ( m + r u * ) ,
J 3 : = − d 1 ξ b u * v * ( m + r u ) ( m + r u * ) ,
当 ξ = d 2 _ d 2 ¯ ( u _ ∞ − ε ) ( u ¯ ∞ − ε ) u _ ∞ 3 u ¯ ∞ d 1 _ d 1 ¯ v _ ∞ v ¯ ∞ 3 1 m + r ( u _ ∞ − ε ) 时,存在 0 < δ ≤ 1 ,使得 B 2 < 4 A C ,有
d E ( t ) d t ≤ − ∫ Ω δ ( u − u * ) 2 + δ ( v − v * ) 2 d 1 d 2 γ u u * ( m + r u ) ( m + r u * ) d x = : ψ ( t ) ≤ 0
根据(3.2)式,存在 G 1 > 0 ,使得 | ψ ( t ) | < G 1 。利用引理2有:
lim t → ∞ ψ ( t ) = − ∫ Ω δ ( u − u * ) 2 + δ ( v − v * ) 2 d 1 d 2 γ u u * ( m + r u ) ( m + r u * ) d x = 0
由(3.1)式,存在 t ≥ T 1 ,有 0 < ε 2 ≤ u ( x , t ) ,得到
lim t → ∞ u ( x , t ) = u * ( x ) , lim t → ∞ v ( x , t ) = v * ( x ) . (5.1)
同时,(3.2)式也表明集合 { u ( ⋅ , t ) : t ≥ 1 } 在 C 2 ( Ω ¯ ) 上是相对紧凑的,可以假设
‖ u ( x , t p ) − u ˜ ( x ) ‖ C 2 ( Ω ¯ ) , ‖ v ( x , t p ) − v ˜ ( x ) ‖ C 2 ( Ω ¯ ) → 0 , t p → ∞
可以得出,在 Ω ¯ 上有
u ˜ ( x ) ≡ u * ( x ) , v ˜ ( x ) ≡ v * ( x ) .
因此, lim t → ∞ u ( x , t ) = u * ( x ) , lim t → ∞ v ( x , t ) = v * ( x ) 在 C 2 ( Ω ¯ ) 上一致成立,那么正稳态解 ( u * , v * ) 是全局稳定的,证毕。
针对Holling-Tanner捕食–食饵模型,本文研究了资源空间分散和资源分布不均对种群的影响。本文首先利用抛物方程的比较定理和上下解方法证明了方程(1.1)正稳态解的存在性,进一步通过迭代过程,得出了问题(1.1)正稳态解的下界。最后利用新构造的李雅普诺夫函数,得到Holling-Tanner捕食–食饵模型在异质环境下的正稳态解是全局稳定的。
国家自然科学基金青年项目(52004215);陕西省自然科学基金青年项目(2019JQ-755);陕西省教育厅基金(19JK0462)。
张万红. 异质环境下扩散Holling-Tanner模型正稳态解的全局稳定性Global Stability of the Positive Steady-State Solution of the Diffusion Holling-Tanner Model in a Heterogeneous Environment[J]. 理论数学, 2021, 11(09): 1665-1672. https://doi.org/10.12677/PM.2021.119185