一类n-斑块Lotka-Volterra型竞争系统的定性分析
Qualitative Analysis for a Lotka-Volterra n-Patch Competition System
摘要: 本文研究了扩散和竞争能力对n个斑块两个物种的Lotka-Volterra竞争系统的动力学行为的联合影响。主要利用单调动力系统理论等方法建立了半平凡稳态解和共存稳态解的渐近稳定性。此外,在论文中还引入了一种新的矩阵主特征值的定义。
Abstract: In this paper, we investigate the joint effects of diffusion and competition ability on the dynamics of a n-patch two-species Lotka-Volterra competition system. The asymptotic stability of semi-trivial and coexistence steady states is established mainly by using the theory of monotone dynamical systems. Moreover, we point out that a new definition of principal eigenvalue of matrix is introduced in this paper.
文章引用:李娟. 一类n-斑块Lotka-Volterra型竞争系统的定性分析[J]. 理论数学, 2022, 12(2): 240-251. https://doi.org/10.12677/PM.2022.122028

1. 引言

在生态学中生物种群与环境间的相互作用以及种群与种群间的相互作用可以利用动力学方法研究,并由此产生种群动力学。而种群动力学的研究,一直是生态学研究的核心问题之一。生物体的扩散和带有不均匀资源分配的空间异质环境之间的联系在近些年得到广泛关注。对扩散稳定策略的研究已经成为种群动力学的基本研究目标。离散栖息地上两个竞争物种的扩散率对稳定性的影响是种群动力学研究中重要的一方面,且结合数学的微分方程可以合理化的解释种群动力学。近年来,离散栖息地上两个物种的Lotka-Volterra竞争系统的全局动力学行为已经得到广泛研究(见 [1] [2] [3] [4] [5] )。

2019年,J. Xiang和Y. Fang在 [6] 中考虑了如下两个斑块的竞争模型:

{ d u 1 d t = d ( u 2 u 1 ) q u 1 + u 1 ( 1 u 1 v 1 ) d u 2 d t = d ( u 1 2 u 2 ) + q u 1 q u 2 + u 2 ( 1 u 2 v 2 ) d v 1 d t = D ( v 2 v 1 ) q v 1 + v 1 ( 1 u 1 v 1 ) d v 2 d t = D ( v 1 2 v 2 ) + q v 1 q v 2 + v 2 ( 1 u 2 v 2 )

对如上模型,作者们得出存在一个中间扩散率是进化稳定的。

2020年,H. Jiang,K.-Y. Lam和Y. Lou在 [7] 中研究了三个斑块两个竞争物种模型的动力学行为,并说明了定向河流网络模块的拓扑结构是如何影响扩散演变。对带有不同拓扑结构的河流网络模块,物种在斑块间的扩散是双向,沿着河流方向的漂流是单向,这是对流和无对流的结合。模型一假设有两个上游斑块和一个下游斑块,并且两个上游斑块间没有直接联系。模型二假设有上游,中游和下游三个斑块。模型三假设有一个上游斑块和两个下游斑块,并且两个下游斑块间没有直接联系。这三种模型都假定上中下游三个斑块的承载力由大到小,这意味上游斑块承载力更大,更有利于物种生存。

并且证明当对流较小时,对三个模型来说都是两个竞争物种中拥有较低扩散率的物种获胜。当对流较大时,对模型一和二来说没有奇异策略,并且两个竞争物种中拥有较高扩散率的物种获胜。对模型三来说扩散率为零和无穷是收敛稳定策略,并且存在奇异策略。对中间范围的对流来说,模型一和二存在奇异策略。而对模型三这个策略可能不存在。

2021年,S. Chen等人在 [1] 中考虑如下n个斑块两个物种的Lotka-Volterra竞争模型:

{ u i = d 1 j = 1 n ( a i j u j a j i u i ) + u i ( p i u i c v i ) , i = 1 , 2 , , n , t > 0 , v i = d 2 j = 1 n ( a i j v j a j i v i ) + v i ( q i b u i v i ) , i = 1 , 2 , , n , t > 0.

在弱竞争以及矩阵 L = ( L i j ) n × n 不可约且其加权有向图是循环平衡的假设下,给出Lotka-Volterra竞争

斑块模型关于全局动力学行为的分类。

本文在原有结论的基础上,对两个竞争物种在n个斑块上分别具有不同资源时,研究扩散率的确切位置,对系统的动力学行为给出更进一步的分类。本文的结果是对以有结论 [1] 的推广。

2. 预备知识和基本引理

本文主要研究如下n个斑块两个物种的Lotka-Volterra竞争模型:

{ u i = d 1 j = 1 n ( a i j u j a j i u i ) + u i ( p i u i c v i ) , i = 1 , 2 , , n , t > 0 , v i = d 2 j = 1 n ( a i j v j a j i v i ) + v i ( q i b u i v i ) , i = 1 , 2 , , n , t > 0 , (2.1)

其中 u = ( u 1 , u 2 , , u n ) v = ( v 1 , v 1 , , v n ) 分别代表两个竞争物种在n个斑块的种群密度; n 2 是一个整数; d 1 , d 2 > 0 分别是物种u和v的扩散率; b , c > 0 代表两个物种的种间竞争系数; p i , q i > 0 分别代表物

u i , v i 在斑块i的承载能力;以及矩阵 A = ( a i j ) n × n 描述斑块间的移动模式,其中 a i j 0 ( i j ) 表示斑块j到斑块i的移动程度。

在本文中,向量 u = ( u 1 , u 2 , , u n ) > ( ) 0 表示u每个分量都是正的(非负的)。记 L = ( L i j ) n × n 为模型的

连接矩阵,其中

L i j = { a i j , i j , k i n a k i , i = j . (2.2)

因此,本文考虑如下改写的模型:

{ u i = d 1 j = 1 n L i j u j + u i ( p i u i c v i ) , i = 1 , 2 , , n , t > 0 , v i = d 2 j = 1 n L i j v j + v i ( q i b u i v i ) , i = 1 , 2 , , n , t > 0 , u ( 0 ) = u 0 ( ) 0 , v ( 0 ) = v 0 ( ) 0 . (2.3)

为了描述本文关于系统(2.3)的结果,我们作出如下基本假设:

(A1) A = ( a i j ) n × n (或L)对称,且对所有的 i j a i j > 0

(A2) p = ( p 1 , p 2 , , p n ) > 0 q = ( q 1 , q 2 , , q n ) > 0

我们记 w ( d , r ) = ( w 1 , w 2 , , w n ) 是如下方程组的唯一正解:

d j = 1 n L i j w j + w i ( r i w i ) = 0 , i = 1 , 2 , , n . (2.4)

(2.4)的全局动力学行为由 [8] [9] [10] [11] 给出。为了简化符号,以下记

u * ( d 1 ) = w ( d 1 , p ) v * ( d 2 ) = w ( d 2 , q ) ,

其中 u * ( d 1 ) = ( u 1 * , u 2 * , , u n * ) , v * ( d 2 ) = ( v 1 * , v 2 * , , v n * )

显然(2.3)总是有两个半平凡稳态解,以下我们用 ( u * ( d 1 ) , 0 ) ( 0 , v * ( d 2 ) ) 表示。为了精确刻画其线性稳定性,对系统(2.3)在 Q = + × + 中定义如下三个子集:

Σ u : = { ( d 1 , d 2 ) Q | ( u * ( d 1 ) , 0 ) 线 } , Σ v : = { ( d 1 , d 2 ) Q | ( 0 , v * ( d 2 ) ) 线 } , Σ : = { ( d 1 , d 2 ) Q | ( u * ( d 1 ) , 0 ) ( 0 , v * ( d 2 ) ) 线 } . (2.5)

定义2.1 给定一个常数 d > 0 和一个向量 h = ( h 1 , h 2 , , h n ) n ,我们定义 μ 1 ( d , h ) 为如下特征值问题的主特征值:

d j = 1 n L i j ψ j + h i ψ i + μ ψ i = 0 , i = 1 , 2 , , n . (2.6)

现在我们定义Q的另外三个子集:

Σ u , 0 : = { ( d 1 , d 2 ) Q | μ 1 ( d 2 , q b u * ( d 1 ) ) = 0 } , Σ v , 0 : = { ( d 1 , d 2 ) Q | μ 1 ( d 1 , p c v * ( d 2 ) ) = 0 } , Σ 0 , 0 : = { ( d 1 , d 2 ) Q | μ 1 ( d 2 , q b u * ( d 1 ) ) = μ 1 ( d 1 , p c v * ( d 2 ) ) = 0 } . (2.7)

受 [12] 的启发,定义

L u : = inf d 1 > 0 i = 1 n q i i = 1 n u i * ( d 1 ) ( 0 , ) , S u : = sup d 1 > 0 max 1 i n q i u i * ( d 1 ) ( 0 , ) , (2.8)

L v : = inf d 2 > 0 i = 1 n p i i = 1 n v i * ( d 2 ) ( 0 , ) , S v : = sup d 2 > 0 max 1 i n p i v i * ( d 2 ) ( 0 , ) . (2.9)

下面我们描述当改变b (或c)时,集合 Σ u Σ u ,0 (或 Σ v Σ v ,0 )在 d 1 d 2 -平面上如何变化。为了刻画对 b > 0 集合 Σ u 的特征,对每个 b > 0 定义,

I b : = { d 1 > 0 | i = 1 n ( q i b u i * ( d 1 ) ) < 0 } = I b 0 I b 1 , (2.10)

其中

I b 0 : = { d 1 > 0 | q b u * ( d 1 ) ( ) 0 } , I b 1 : = { d 1 I b | max 1 i n ( q i b u i * ( d 1 ) ) > 0 } . (2.11)

同样地为了刻画对 c > 0 集合 Σ v 的特征,对每个 c > 0 定义,

I c : = { d 2 > 0 | i = 1 n ( p i c v i * ( d 2 ) ) < 0 } = I c 0 I c 1 , (2.12)

其中

I c 0 : = { d 2 > 0 | p c v * ( d 2 ) ( ) 0 } , I c 1 : = { d 2 I c | max 1 i n ( p i c v i * ( d 2 ) ) > 0 } . (2.13)

Σ u ¯ Σ v ¯ 分别为 Σ u Σ v 在Q中的闭包,则 Σ u = Σ u ¯ \ Σ u Σ u , 0 Σ v = Σ v ¯ \ Σ v Σ v , 0

如下给出本文需要的基本引理。

引理 2.1设(A1)成立,且 r = ( r 1 , r 2 , , r n ) > 0 以及r的坐标不全相等。则如下结论成立:

i) 对任意 d > 0 i = 1 n r i i = 1 n w i ( d , r ) i = 1 n r i = i = 1 n w i ( d , r ) 当且仅当 r 1 = r 2 = = r n

ii) w ( d , r ) 连续依赖于 d > 0 ,此外,

w i ( d , r ) { r i d 0 + , i = 1 , 2 , , n , i = 1 n r i / n d , i = 1 , 2 , , n ;

iii) 对任意 d > 0 max 1 i n w i ( d , r ) < max 1 i n r i

证明:首先证明引理2.1-(i)。对(2.4)两边乘以 1 / w i ( d , r ) 并把所有方程相加,可推出

d ( i = 1 n L i i + i = 1 n j = 1 , j i n L i j w j ( d , r ) w i ( d , r ) ) + i = 1 n r i i = 1 n w i ( d , r ) = 0.

由(2.2)和(A1),有

d ( i = 1 n j = 1 , j i n a i j + i = 1 n j = 1 , j i n a i j w j ( d , r ) w i ( d , r ) ) + i = 1 n r i i = 1 n w i ( d , r ) = 0. (2.14)

b i j = a i j w j ( d , r ) w i ( d , r ) ,易证对于任意 1 i , j n b i j + b j i 2 a i j 。因此,

i = 1 n j = 1 , j i n a i j w j ( d , r ) w i ( d , r ) i = 1 n j = 1 , j i n a i j . (2.15)

由(2.14)和(2.15),我们得出

i = 1 n r i i = 1 n w i ( d , r ) . (2.16)

注意到对任意 i j b i j + b j i = 2 a i j 当且仅当 w 1 ( d , r ) = w 2 ( d , r ) = = w n ( d , r ) 。因此,

i = 1 n r i = i = 1 n w i ( d , r ) r 1 = r 2 = = r n .

接下来证明引理2.1-(ii)。 w ( d , r ) 关于d的连续依赖性可由隐含函数定理(见 [13],命题3.6)得到。然后对(2.4)令 d 0 + ,由(2.16)可得 w i ( d , r ) r i i = 1 , 2 , , n

对(2.4)两边乘以1/d并令 d ,我们得出

0 = lim d j = 1 n L i j w j ( d , r ) = j = 1 n L i j ( lim d w j ( d , r ) ) , i = 1 , 2 , , n . (2.17)

由假设(A1),(2.16)和(2.17)可知,有 lim d w 1 ( d , r ) = = lim d w n ( d , r ) 0 。注意到对所有 1 j n ,有 i = 1 n L i j = 0 。然后把(2.4)所有方程相加并令 d ,可得出

w i ( d , r ) i = 1 n r i / n , i = 1 , 2 , , n .

最后,我们通过反证法验证引理2.1-(iii)。若存在某个 i 0 使得

w i 0 ( d , r ) = max 1 i n w i ( d , r ) max 1 i n r i r i 0 , (2.18)

则由(2.2),(2.4)和(2.18),可得出

j = 1 n L i 0 j w j ( d , r ) = j = 1 , j i 0 n a i 0 j w i 0 ( d , r ) + j = 1 , j i 0 n a i 0 j w j ( d , r ) 0 ,

这表明

w i 0 ( d , r ) j = 1 , j i 0 n a i 0 j w j ( d , r ) j = 1 , j i 0 n a i 0 j max 1 j n j i 0 w j ( d , r ) w i 0 ( d , r ) .

显然 w 1 ( d , r ) = w 2 ( d , r ) = = w n ( d , r ) 。根据假设r的坐标不全相等可知,这不可能成立。从而完成了引理2.1证明。

引理2.2 假设(A1)成立,则(2.6)的主特征值 μ 1 ( d , h ) 由如下公式给出

μ 1 ( d , h ) = inf ψ n ψ 0 d i = 1 n j = 1 n L i j ψ i ψ j i = 1 n h i ψ i 2 i = 1 n ψ i 2 (2.19)

且主特征值 μ 1 ( d , h ) 的特征向量为正。

证明:为了方便起见,记

R ( ψ ) = d i = 1 n j = 1 n L i j ψ i ψ j i = 1 n h i ψ i 2 i = 1 n ψ i 2 .

注意到

i = 1 n j = 1 n L i j ψ i ψ j = ( ψ 1 , ψ 2 , , ψ n ) ( k 1 a 1 k a 12 a 1 n a 12 k 2 a 2 k a 2 n a 1 n a 2 n k n a n k ) ( ψ 1 ψ 2 ψ n ) 0. (2.20)

R ( ψ ) i = 1 n h i ψ i 2 i = 1 n ψ i 2 max 1 i n | h i | ,

即对任意 ψ n \ { 0 } R ( ψ ) 下有界。因此, inf ψ n ψ 0 R ( ψ ) >

步骤1:定义 μ 1 ( d , h ) = inf ψ n ψ 0 R ( ψ ) ,证明 μ 1 ( d , h ) 是特征值。由(2.20)和

i = 1 n j = 1 n L i j ψ i ψ j = 0 ψ 1 = ψ 2 = = ψ n , (2.21)

我们在 n 中引入一个内积

( ϕ , ψ ) * = d i = 1 n j = 1 n L i j ϕ i ψ j + i = 1 n ( h i + 2 max 1 i n | h i | ) ϕ i ψ i ,

从而在 n 中产生一个范数 ψ * = ( ψ , ψ ) * 。以及记 ψ = ( i = 1 n ψ i 2 ) 1 2

现在取极小化序列 { ψ m } m = 1 n ,我们有当 m 时, R ( ψ m ) μ 1 ( d , h ) 。由于 R ( ψ m ) = R ( ψ m ψ m ) ,所以不失一般性,我们假设对任意 m 1 ψ m = 1 。因此,

R ( ψ m ) = d i = 1 n j = 1 n L i j ψ i m ψ j m i = 1 n h i ( ψ i m ) 2 = ( ψ m * ) 2 2 max 1 i n | h i | .

因为 { R ( ψ m ) } m = 1 有界,所以 { ψ m * } m = 1 有界,即 { ψ m } m = 1 有界。从而存在子列(仍记为 { ψ m } ),当 m 时,有

i = 1 n ( ψ i m ) 2 i = 1 n ( ψ i ) 2 ψ m * ψ * ,

其中 ψ n ψ = 1 。因此,

R ( ψ ) = ( ψ * ) 2 2 max 1 i n | h i | ψ m 2 = lim m ( ( ψ m * ) 2 2 max 1 i n | h i | ψ m 2 ) = lim m R ( ψ m )

这表明 R ( ψ ) = inf ψ n ψ 0 R ( ψ )

下面我们断言 ( μ 1 ( d , h ) , ψ ) 是一个特征对,会完成步骤1的证明。为了证明断言成立,定义对任意 ϕ n f ( s ) = R ( ψ + s ϕ ) 。因此, f ( 0 ) = inf s f ( s ) f ( 0 ) = 0 。由

0 = f ( 0 ) = ( d i = 1 n j = 1 n L i j ( ϕ i ψ j + ϕ j ψ i ) 2 i = 1 n h i ϕ i ψ i ) i = 1 n ( ψ i ) 2 ( i = 1 n ( ψ i ) 2 ) 2 2 ( d i = 1 n j = 1 n L i j ψ i ψ j i = 1 n h i ( ψ i ) 2 ) i = 1 n ϕ i ψ i ( i = 1 n ( ψ i ) 2 ) 2 ,

可以推出

d i = 1 n j = 1 n L i j ϕ i ψ j i = 1 n h i ϕ i ψ i = ( d i = 1 n j = 1 n L i j ψ i ψ j i = 1 n h i ( ψ i ) 2 ) i = 1 n ϕ i ψ i = μ 1 ( d , h ) i = 1 n ϕ i ψ i .

从而由 ϕ 的任意性,我们有 d j = 1 n L i j ψ j + h i ψ i + μ 1 ( d , h ) ψ i = 0 ,( i = 1 , 2 , , n )。因此断言得证。

步骤2:我们证明 μ 1 ( d , h ) 的特征向量 ψ > 0 。首先证明 ψ ( ) 0 ,由(2.2)和 a i j 0 ( i j ),我们得到 i = 1 n j = 1 n L i j ψ i ψ j i = 1 n j = 1 n L i j | ψ i | | ψ j |

因此,

μ 1 ( d , h ) R ( | ψ | ) d i = 1 n j = 1 n L i j ψ i ψ j i = 1 n h i ( ψ i ) 2 i = 1 n ( ψ i ) 2 = R ( ψ ) = μ 1 ( d , h ) . (2.22)

从而由(2.22)和 ψ = 1 ,可推得 ψ ( ) 0 。接下来我们断言 ψ > 0 。反证法,假设存在某个 i 0 ,使得 ψ i 0 = 0 。注意到 ψ i 0 满足 d j = 1 n L i 0 j ψ j + h i 0 ψ i 0 + μ 1 ( d , h ) ψ i 0 = 0 。因此由假设(A1),我们得到对所有 j i 0 ψ j = 0 。由 ψ 0 可知,这不可能成立。从而完成了该引理的证明。

为了进一步刻画主特征值 μ 1 ( d , h ) 的性质,需引入如下特征值问题:

j = 1 n L i j ϕ j + λ h i ϕ i = 0 , i = 1 , 2 , , n , (2.23)

其中 h n ,h的坐标不全相等以及 h 1 , h 2 , , h n 符号不全相同。若(2.23)存在一个正解,我们称 λ 是主特征值。

引理2.3 [13] 假设(A1)和 i = 1 n h i < 0 成立。设 { ϕ m } n 中的一个序列且满足 i = 1 n ( ϕ i m ) 2 = 1 i = 1 n h i ( ϕ i m ) 2 > 0 ,则存在一个常数 c 0 > 0 ,使得对所有m, i = 1 n j = 1 n L i j ϕ i m ϕ j m c 0

引理2.4 [13] 假设(A1)成立以及 h 1 , h 2 , , h n 的符号不全相同,则问题(2.23)存在正主特征值 λ 1 = λ 1 ( h ) 当且仅当 i = 1 n h i < 0 ,并且正主特征值 λ 1 ( h ) 由如下式子决定

1 λ 1 ( h ) = sup ϕ n i , j = 1 n L i j ϕ i ϕ j 0 i = 1 n h i ϕ i 2 i = 1 n j = 1 n L i j ϕ i ϕ j . (2.24)

该引理的证明与 [13] 中的证明类似,这里省略。

为了证明命题2.1,我们先给出如下引理。

引理2.5 [13] 假设(A1)成立, λ 是一个正参数并且 i = 1 n h i < 0 h 1 , h 2 , , h n 的符号不全相同。则

j = 1 n L i j ψ j + λ h i ψ i + σ ψ i = 0 , i = 1 , 2 , , n , (2.25)的主特征值 σ 1 > 0 当且仅当 0 < λ < λ 1 ,其中 λ 1 = λ 1 ( h ) 是(2.23)的正主特征值。

从而若 i = 1 n h i < 0 h 1 , h 2 , , h n 的符号不全相同,则如下结论成立:

{ σ 1 < 0 λ > λ 1 , σ 1 = 0 λ = λ 1 , σ 1 > 0 λ < λ 1 .

命题 2.1假设(A1)成立,则(2.6)的主特征值 μ 1 ( d , h ) 连续依赖于 d > 0 。此外, μ 1 ( d , h ) 有如下性质:

i) 若 i = 1 n h i 0 h 0 ,则对所有 d > 0 μ 1 ( d , h ) < 0

ii) 若 i = 1 n h i < 0 h 1 , h 2 , , h n 的符号不全相同,则

{ μ 1 ( d , h ) < 0 d < 1 / λ 1 ( h ) , μ 1 ( d , h ) = 0 d = 1 / λ 1 ( h ) , μ 1 ( d , h ) > 0 d > 1 / λ 1 ( h ) ;

iii) 若 h ( ) 0 ,则对所有 d > 0 μ 1 ( d , h ) > 0

命题2.1-(i)和(iii)的证明是基本的,这里略去。命题2.1-(ii)的证明可由引理2.5得出。

引理2.6 [14] 系统(2.3)中 ( u * ( d 1 ) , 0 ) , ( 0 , v * ( d 2 ) ) ( 0 , 0 ) 的线性稳定性,分别由 μ 1 ( d 2 , q b u * ( d 1 ) ) , μ 1 ( d 1 , p c v * ( d 2 ) ) min { μ 1 ( d 1 , p ) , μ 1 ( d 2 , q ) } 的符号决定。

该引理的证明类似于 [14] 中的推论2.10,这里省略。

3. 主要结果

本文的主要结果如下:

定理3.1假设(A1)和(A2)成立且p和q的坐标至少有一个不全相等, L u , S u , L v S v 的定义如(2.8)和(2.9),以下结论对系统(2.3)成立:

i) 对 Σ u ,我们有如下特征:

Σ u = { 0 < b L u , { ( d 1 , d 2 ) | d 1 I b , d 2 > d 2 * ˜ ( d 1 ) } Q L u < b < S u , Q b S u , (3.1)

其中 I b , I b 0 I b 1 定义如(2.10)和(2.11),且在 I b d 2 * ˜ ( d 1 ) 定义如下:

d 2 * ˜ ( d 1 ) : = { 0 d 1 I b 0 , 1 λ 1 ( q b u * ( d 1 ) ) > 0 d 1 I b 1 , (3.2)

其中 λ 1 定义如(2.24)。因此 Σ u 当且仅当 b > L u

ii) 对 Σ v ,我们有如下特征:

Σ v = { 0 < c L v , { ( d 1 , d 2 ) | d 2 I c , d 1 > d 1 * ˜ ( d 2 ) } Q L v < c < S v , Q c S v , (3.3)

其中 I c , I c 0 I c 1 定义如(2.12)和(2.13),且在 I c d 1 * ˜ ( d 2 ) 定义如下:

d 1 * ˜ ( d 2 ) : = { 0 d 2 I c 0 , 1 λ 1 ( p c v * ( d 2 ) ) > 0 d 2 I c 1 , (3.4)

其中 λ 1 定义如(2.24)。因此, Σ v 当且仅当 c > L v

iii) 对 Σ u ,0 ,我们有如下特征:

Σ u , 0 = { 0 < b < L u b S u , { ( d 1 , d 2 ) | q L u u * ( d 1 ) } b = L u , Σ u { ( d 1 , d 2 ) | q b u * ( d 1 ) } L u < b < S u ; (3.5)

iv) 对 Σ v ,0 ,我们有如下特征:

Σ v , 0 = { 0 < c < L v c S v , { ( d 1 , d 2 ) | p L v v * ( d 2 ) } c = L v , Σ v { ( d 1 , d 2 ) | p c v * ( d 2 ) } L v < c < S v . (3.6)

证明:首先证明定理3.1-(i)成立。由引理2.6可知 Σ u = { ( d 1 , d 2 ) | μ 1 ( d 2 , q b u * ( d 1 ) ) > 0 } 。假设 d 1 I b ,其中 I b 定义如(2.10),则 i = 1 n ( q i b u i * ( d 1 ) ) 0 。由命题2.1-(i)可知对所有 d 2 ,有 μ 1 ( d 2 , q b u * ( d 1 ) ) 0 ,即对所有 d 2 > 0 ( d 1 , d 2 ) Σ u 。所以 ( d 1 , d 2 ) Σ u 表明 d 1 I b 。下面刻画对所有 b > 0 ,集合 I b 的特征。若 b L u ,则由(2.8)中 L u 定义可知对所有 d 1 > 0 ,我们有 i = 1 n ( q i b u i * ( d 1 ) ) 0 。因此 I b = Σ u = 。此外,

I b 1 L u < b < S u . (3.7)

d 1 I b 1 ,则存在某个 i 0 使得 L u < b q i 0 b u i 0 * ( d 1 ) > 0 成立。因此

L u < b < q i 0 u i 0 * ( d 1 ) S u .

另一方面,若 L u < b < S u ,则存在某个 d 1 > 0 j 0 使得

q j 0 b u j 0 * ( d 1 ) > 0 i = 1 n ( q i b u i * ( d 1 ) ) < 0 ,

这表明 d 1 I b 1 。这就完成了(3.7)的证明。

下面我们断言 I b 有如下分解:

{ I b = b L u , I b = I b 0 I b 1 + L u < b < S u , I b = I b 0 = + b S u . (3.8)

为了证明(3.8)成立,我们只需证明若 b S u ,则 I b = I b 0 = + 。根据(2.8)中 S u 的定义和 b S u ,我们得出对所有 d 1 > 0 ,有 q b u * ( d 1 ) 0 。因此为了证明 I b = I b 0 = R + 成立,只需证明若 b S u ,则对所有 d 1 > 0 q b u * ( d 1 ) 0 成立。容易得出 b > S u 结论成立。

接下来,我们只需通过假设 b = S u 证明对所有 d 1 > 0 q b u * ( d 1 ) 0 成立即可。若 p 1 = p 2 = = p n = p * ,则由假设可知q的坐标不全相等以及对所有 d 1 > 0 ,我们有 u i * ( d 1 ) = p * ( i = 1 , 2 , , n )。因此对所有 d 1 > 0 q b u * ( d 1 ) 的坐标不全相等。若 q 1 = q 2 = = q n = q * ,则由假设可知p的坐标不全相等且对所有 d 1 > 0 u * ( d 1 ) 的坐标不全相等。因此对所有 d 1 > 0 q b u * ( d 1 ) 的坐标不全相等。如果p和q的坐标都不全相等,我们用反证法证明 q b u * ( d 1 ) 0 成立。若存在某个 d 1 > 0 使得 q b u * ( d 1 ) 0 ,则由(2.8)中 S u 的定义,可得出对所有 i = 1 , 2 , , n d 1 > 0 ,有

q i u i * ( d 1 ) = S u q i u i * ( d 1 ) max 1 i n q i u i * ( d 1 ) .

因此对所有 d 1 d 1 u * ( d 1 ) u * ( d 1 ) ( ) 0 成立。但由引理2.1可知,这与 i = 1 n u i * ( d 1 ) 不能在某个有限

d 1 ( 0, ) 处达到全局最小相矛盾。从而(3.8)得证。所以由(3.8)和命题2.1-(iii),我们得出对所有 ( d 1 , d 2 ) Q ,有

b S u μ 1 ( d 2 , q b u * ( d 1 ) ) > 0. (3.9)

下面假设 L u < b < S u d 1 I b = I b 0 I b 1 成立。若 d 1 I b 0 ,则由命题2.1-(iii)可知对所有 d 2 > 0 = d 2 * ˜ ( d 1 ) ,有 μ 1 ( d 2 , q b u * ( d 1 ) ) > 0 。若 d 1 I b 1 ,则由命题2.1-(ii)可知对所有 d 2 > 1 / λ 1 ( q b u * ( d 1 ) ) = d 2 * ˜ ( d 1 ) ,有 μ 1 ( d 2 , q b u * ( d 1 ) ) > 0 。因此, ( d 1 , d 2 ) Σ u 当且仅当 d 1 I b d 2 > d 2 * ˜ ( d 1 ) ,从而(3.2)得证。由(3.2)和(3.7),可推出当 b ( L u , S u ) 时,有 Σ u Q 成立。这就完成了定理3.1-(i)的证明。定理3.1-(ii)的证明与定理3.1-(i)相似,这里省略。

最后我们证明定理3.1-(iii)。若 b L u ,则对所有 d 1 > 0 ,有 i = 1 n ( q i b u i * ( d 1 ) ) 0 。因此由命题2.1-(i),

我们得出对所有 ( d 1 , d 2 ) Q ,有 μ 1 ( d 2 , q b u * ( d 1 ) ) 0 并且 μ 1 ( d 2 , q b u * ( d 1 ) ) = 0 (即 ( d 1 , d 2 ) Σ u ,0 )当且仅当 b = L u q L u u * ( d 1 ) 。这表明如果 b < L u ,则 Σ u , 0 = 成立。因此只需证明 b = L u 即可。若存在某个 d 1 * > 0 ,使得 q L u u * ( d 1 * ) ,则 d 1 * 唯一。事实上,如果p的坐标不全相等且这样的 d 1 * > 0 存在,则由 u i * ( d 1 * ) 所满足的方程可得到对所有 i = 1 , 2 , , n

d 1 * u i * ( d 1 * ) ( p i u i * ( d 1 * ) ) j = 1 n L i j u j * ( d 1 * ) q i ( L u p i q i ) L u j = 1 n L i j q j ,

显然 d 1 * 唯一。若 p 1 = p 2 = = p n = p * ,则由假设可知q的坐标不全相等并且对所有 d 1 > 0 ,有 u i * ( d 1 ) = p * ( i = 1 , 2 , , n )。这表明这样的 d 1 * 不存在且 Σ u , 0 = 。若 b S u ,则由(3.9)可推出 Σ u , 0 =

因此,为了完成定理3.1-(iii)的证明,我们只需证明当 L u < b < S u 时结论成立即可。由定理3.1-(i)可知当 L u < b < S u 时,有 Σ u Q 。假设 ( d 1 , d 2 ) Σ u ,0 \ Σ u ,我们断言 q b u * ( d 1 ) 0 。事实上,若 q b u * ( d 1 ) 0 ,则由 μ 1 ( d 2 , q b u * ( d 1 ) ) = 0 和命题2.1-(i)和(iii),可得出 i = 1 n ( q i b u i * ( d 1 ) ) < 0 以及 q b u * ( d 1 ) 坐标的符号不全相同。从而根据引理2.4和命题2.1-(ii),我们推出 1 / d 2 = λ 1 ( q b u * ( d 1 ) ) > 0 且对所有 d 2 > d 2 μ 1 ( d 2 , q b u * ( d 1 ) ) > 0 (即 ( d 1 , d 2 ) Σ u )。这表明 ( d 1 , d 2 ) Σ u ,但这与 ( d 1 , d 2 ) Σ u ,0 \ Σ u 矛盾,因此 q b u * ( d 1 ) 0 。通过与上述 b = L u 证明类似的方法,我们可以得到若存在某个 d 1 ,使得 q b u * ( d 1 ) ,则 d 1 唯一。这样就完成了定理3.1-(iii)的证明。定理3.1-(iv)的证明与定理3.1-(iii)相似,这里省略。

参考文献

[1] Chen, S., Shi, J., Shuai, Z. and Wu, Y. (2022) Global Dynamics of a Lotka-Volterra Competition Patch Model. Non-linearity, 35, 817-842.
[2] Cheng, C.-Y., Lin, K.-H. and Shih, C.-W. (2019) Coexistence and Extinction for Two Competing Species in Patchy Environments. Mathematical Biosciences and Engineering, 16, 909-946.
https://doi.org/10.3934/mbe.2019043
[3] Gourley, S.A. and Kuang, Y. (2005) Two-Species Competition with High Dispersal: The Winning Strategy. Mathematical Biosciences and Engineering, 2, 345-362.
https://doi.org/10.3934/mbe.2005.2.345
[4] Lin, K.-H., Lou, Y., Shih, C.-W. and Tsai, T.-H. (2014) Global Dy-namics for Two-Species Competition in Patchy Environment. Mathematical Biosciences and Engineering, 11, 947-970.
https://doi.org/10.3934/mbe.2014.11.947
[5] Slavík, A. (2020) Lotka-Volterra Competition Model on Graphs. SIAM Journal on Applied Dynamical Systems, 19, 725-762.
https://doi.org/10.1137/19M1276285
[6] Xiang, J. and Fang, Y. (2019) Evolutionarily Stable Dispersal Strategies in a Two-Patch Advective Environment. Discrete and Continuous Dynamical Systems-B, 24, 1875-1887.
https://doi.org/10.3934/dcdsb.2018245
[7] Jiang, H., Lam, K.-Y. and Lou, Y. (2020) Are Two-Patch Models Sufficient? The Evolution of Dispersal and Topology of River Network Modules. Bulletin of Mathematical Biology, 131, 1-42.
https://doi.org/10.1007/s11538-020-00803-1
[8] Chen, S., Shi, J., Shuai, Z. and Wu, Y. (2019) Spectral Monotonicity of Perturbed Quasi-Positive Matrices with Applications in Population Dynamics.
[9] Cosner, C. (1996) Variability, Vagueness and Comparison Methods for Ecological Models. Bulletin of Mathematical Biology, 58, 207-246.
https://doi.org/10.1007/BF02458307
[10] Li, M.Y. and Shuai, Z. (2010) Global-Stability Problem for Coupled Systems of Differential Equations on Networks. Journal of Differential Equations, 248, 1-20.
https://doi.org/10.1016/j.jde.2009.09.003
[11] Lu, Z.Y. and Takeuchi, Y. (1993) Global Asymptotic Behavior in Single-Species Discrete Diffusion Systems. Journal of Mathematical Biology, 32, 67-77.
https://doi.org/10.1007/BF00160375
[12] He, X. and Ni, W.-M. (2016) Global Dynamics of the Lotka-Volterra Competition-Diffusion System: Diffusion and Spatial Heterogeneity I. Communications on Pure and Applied Mathe-matics, 69, 981-1014.
https://doi.org/10.1002/cpa.21596
[13] Cantrell, R.S. and Cosner, C. (2003) Spatial Ecology via Reaction-Diffusion Equations. Wiley Series in Mathematical and Computational Biology. John Wiley & Sons Ltd., Chichester, 428 p.
[14] Lam, K.-Y. and Ni, W.-M. (2012) Uniqueness and Complete Dynamics in Heterogeneous Competi-tion-Diffusion Systems. SIAM Journal on Applied Dynamical Systems, 6, 1695-1712.
https://doi.org/10.1137/120869481