一类集合生态系统的一致持久性
The Uniform Persistence of a Meta-Ecosystem
DOI: 10.12677/PM.2020.1011126, PDF, HTML, XML, 下载: 395  浏览: 1,492 
作者: 刘 菁, 杨喜陶:湖南科技大学数学与计算科学学院,湖南 湘潭
关键词: 集合生态系统一致持久性时滞比较原理Meta-Ecosystem Uniform Persistence Delay Comparison Theorem
摘要: 本文研究了一类集合生态系统内部双向资源交换的动力学行为,利用比较原理研究它的一致持久性,得到了该系统一致持久性的充分条件。
Abstract: In this paper we study the dynamic behavior of bi-directional resource exchange within a meta-ecosystem. By using the comparison principle, some sufficient conditions are determined that guarantee the uniform persistence of the model.
文章引用:刘菁, 杨喜陶. 一类集合生态系统的一致持久性[J]. 理论数学, 2020, 10(11): 1061-1069. https://doi.org/10.12677/PM.2020.1011126

1. 引言

生物种群的持续生存问题一直是很多生物数学家所关注的重要课题之一,很多学者研究系统的一致持久性 [1] - [10]。持久性的动力系统方法为处理这类数学问题提供了非常有用的工具,例如平均Lyapunov函数 [3],正常或外部Lyapunov指数 [3] [4] [5],不变概率测度 [3] [5] [6] [7],或链递归,与莫尔斯分解或无环理论 [6] [8] [9] 结合使用。

最近,Messan和Kopp等人为了研究集合生态系统中两个生态系统之间的双向资源交换的动态产出,在 [11] 中建立了如下生态模型:

{ d P d t = r p P [ 1 P K p + a p Q ] b p Q P d Q d t = r q Q [ 1 Q K q + a q P ] b q Q P

其中, P ( t ) , Q ( t ) 是在时间t进行资源交换的两个相邻生态系统的产出, r i , i = p , q 是生态系统中生物量的内在增长率, K i 是在没有资源交换的情况下生态系统i的承载能力; a i 表示从贡献者生态系统转移到接收系统i的资源数量; b i 表示生态系统层面上供体生态系统i的资源枯竭率。 [11] 指出上述系统当

β i = b i r i K j < 1 ( i j )

时为一致持久的。

由于受气候、环境等因素变化的影响, r i , K i , a i , b i 不一定是常数,而是会随着时间的推移而发生变化,即为时间t的函数。

同时,如 [12] [13] 中指出,时滞现象是生物学中另外一种重要现象,说明了群体在成长或消化吸收时的延迟规律。由于种群的相互影响普遍存在滞后现象,近年来具有时滞阶段结构的生物模型得到了广泛的研究。考虑到进行双向资源交换后生物种群的发育时间,生态系统的产出受到生物种群发育情况的影响,即生态系统的资源产出对生物种群数量存在时间滞后现象。基于以上两点考虑,建立如下带有时滞的模型:

{ d x d t = x ( t ) [ r 1 ( t ) r 1 ( t ) x ( t τ 1 ( t ) ) K 1 ( t ) + a 1 ( t ) y ( t ) b 1 ( t ) y ( t ) ] d y d t = y ( t ) [ r 2 ( t ) r 2 ( t ) y ( t τ 2 ( t ) ) K 2 ( t ) + a 2 ( t ) x ( t ) b 2 ( t ) x ( t ) ] , (1)

其中 r i ( t ) , K i ( t ) , a i ( t ) , b i ( t ) , i = 1 , 2 是定义在R上的非负连续有界函数。据我们所知,系统(1)的一致持久性未有相应结果,受 [14] [15] 启发,本文将研究这一问题。我们首先给出一些记号与定义:对定义于R上的非负连续有界函数 g ( t ) ,令

g u = sup t R g ( t ) , g l = inf t R g ( t ) , g = lim sup x g ( t ) , g = lim inf x g ( t ) , τ = max ( τ 1 u , τ 2 u )

φ = ( φ 1 , φ 2 ) C + = { ( φ 1 , φ 2 ) C [ τ , 0 ] , R 2 | φ i ( θ ) 0 , φ i ( 0 ) > 0 , i = 1 , 2 } ,由 [12] 定理2.1知,方程有唯一满足初始条件 x 0 = φ 的解

x ( t , 0 , φ ) = ( x 1 ( t , 0 , φ 1 , φ 2 ) , x 2 ( t , 0 , φ 1 , φ 2 ) )

x i ( t , 0 , φ 1 , φ 2 ) > 0 ( t > 0 , i = 1 , 2 )

2. 一致持久性

定义1 系统(1)的解是一致持久的,如果存在正常数 m i , M i ( i = 1 , 2 ) ,对任意 φ C +

m 1 lim inf t x ( t , 0 , φ ) lim sup t x ( t , 0 , φ ) M 1 ,

m 2 lim inf t y ( t , 0 , φ ) lim sup t y ( t , 0 , φ ) M 2 .

为了本节主要结果的证明,先给出以下两个引理:

引理1 x ( t ) , y ( t ) 是定义在 [ T , + ) 上的非负函数,如果

x x ( m n x ) K + a y ,(2)

lim sup t x m n ,

其中 m , n , K , a 为正常数。

证明:令 g ( t ) = K + a y ,显然 g ( t ) K

由比较原理知:

x [ e T t m g ( s ) d s ( T t n g ( u ) e T u m g ( s ) d s d u + c ˜ ) ] 1 ,

其中 c ˜ = 1 x ( T ) 1 ,令 t + lim sup t x m n

引理2 x ( t ) , y ( t ) 是定义在 [ T , + ) 上的非负函数,如果

x x ( m n x ) K + a y , (3)

lim inf t x m n ,

其中 m , n , K , a 为正常数。

证明类似于引理1,从略。

引理3 M 1 , M 2 > 0 ,使得 φ C + lim sup t x ( t , 0 , φ ) M 1 lim sup t y ( t , 0 , φ ) M 2 。其中

M 1 = { r 1 K 1 exp ( r 1 τ 1 u ) r 1 , r 1 a 1 < K 1 b 1 ( r 1 a 1 + K 1 b 1 ) 2 exp ( r 1 τ 1 u ) 4 a 1 b 1 r 1 , r 1 a 1 K 1 b 1 ,

M 2 = { r 2 K 2 exp ( r 2 τ 2 u ) r 2 , r 2 a 2 < K 2 b 2 ( r 2 a 2 + K 2 b 2 ) 2 exp ( r 2 τ 2 u ) 4 a 2 b 2 r 2 , r 2 a 2 K 2 b 2 .

证明:记 x ( t ) = x ( t , 0 , φ ) , y ( t ) = y ( t , 0 , φ ) 。下面考虑两种情况:

情形(1): r 1 a 1 < K 1 b 1 * 时, ε 0 > 0 ,当 ε ( 0 , ε 0 ) 时,有

( r 1 + ε ) ( a 1 + ε ) < ( K 1 + ε ) ( b 1 * ε ) . (4)

ε ( 0 , ε 0 ) , T 1 > 0 使得当 t > T 1 时:

r 1 ε < r 1 ( t ) < r 1 + ε , K 1 ε < K 1 ( t ) < K 1 + ε (5)

a 1 ε < a 1 ( t ) < a 1 + ε , b 1 ε < b 1 ( t ) < b 1 + ε (6)

t > T 1 时,由(1)和(5)知

x ( t ) ( r 1 + ε ) x ( t ) .

从而

x ( t ) x ( t τ 1 ( t ) ) exp ( ( r 1 + ε ) τ 1 u ) .

代入(1),由(5)及(6)得当 t > T 1 时,

d x d t x ( t ) [ ( r 1 + ε ) ( r 1 ε ) x ( t ) exp ( ( r 1 + ε ) τ 1 u ) ( K 1 + ε ) + ( a 1 + ε ) y ( t ) ( b 1 * ε ) y ( t ) ] = x ( t ) g 1 ( ε ) + g 2 ( ε ) y ( t ) g 3 ( ε ) y 2 ( t ) g 4 ( ε ) x ( t ) ( K 1 + ε ) + ( a 1 + ε ) y ( t )

其中

g 1 ( ε ) = ( r 1 + ε ) ( K 1 + ε ) , g 2 ( ε ) = ( r 1 + ε ) ( a 1 + ε ) ( K 1 + ε ) ( b 1 * ε ) , g 3 ( ε ) = ( a 1 + ε ) ( b 1 * ε ) , g 4 ( ε ) = ( r 1 * ε ) exp ( ( r 1 + ε ) τ 1 u )

从而

x x ( t ) g 1 ( ε ) g 4 ( ε ) x ( t ) ( K 1 + ε ) + ( a 1 + ε ) y ( t ) .

ε 的任意性及引理1,

lim sup t x ( t ) r 1 K 1 exp ( r 1 τ 1 u ) r 1 * = M 1 .

情形(2): r 1 a 1 K 1 b 1 * 时, ε > 0 , T 2 > 0 ,当 t > T 2 时(4)与(5)成立。此时有:

x x ( t ) g 1 ( ε ) + g 2 ( ε ) y ( t ) g 3 ( ε ) y 2 ( t ) g 4 ( ε ) x ( t ) ( K 1 + ε ) + ( a 1 + ε ) y ( t ) = x ( t ) g 1 ( ε ) g 3 ( ε ) ( y ( t ) g 2 ( ε ) 2 g 3 ( ε ) ) 2 + g 2 2 ( ε ) 4 g 3 ( ε ) g 4 ( ε ) x ( t ) ( K 1 + ε ) + ( a 1 + ε ) y ( t ) x ( t ) g 1 ( ε ) + g 2 2 ( ε ) 4 g 3 ( ε ) g 4 ( ε ) x ( t ) ( K 1 + ε ) + ( a 1 + ε ) y ( t )

ε 的任意性及引理1,

lim sup t x ( t ) ( r 1 a 1 + K 1 b 1 * ) 2 exp ( r 1 τ 1 u ) 4 a 1 b 1 * r 1 * = M 1 .

综上所述:

lim sup t x ( t ) M 1 = { r 1 K 1 exp ( r 1 τ 1 u ) r 1 * , r 1 a 1 < K 1 b 1 * ( r 1 a 1 + K 1 b 1 * ) 2 exp ( r 1 τ 1 u ) 4 a 1 b 1 * r 1 * , r 1 a 1 K 1 b 1 * .

同理可得

lim sup t y ( t ) M 2 = { r 2 K 2 exp ( r 2 τ 2 u ) r 2 * , r 2 a 2 < K 2 b 2 * ( r 2 a 2 + K 2 b 2 * ) 2 exp ( r 2 τ 2 u ) 4 a 2 b 2 * r 2 * , r 2 a 2 K 2 b 2 * .

引理4 记 s 1 = r 2 K 2 exp ( r 2 τ 2 u ) r 2 * , s 2 = r 1 K 1 exp ( r 1 τ 1 u ) r 1 * 。当

r i * b i s i > 0 ( i = 1 , 2 ) (7)

时, η i > 0 ( i = 1 , 2 ) ,使得 φ C + , lim sup t x ( t , 0 , φ ) η 1 , lim sup t y ( t , 0 , φ ) η 2

证明:由(7)知存在 η 1 > 0 ,使得

( r 1 * b 1 s 1 ( η 1 ) ) K 1 * exp ( r 1 * r 1 M 1 K 1 * b 1 s 1 ( η 1 ) ) τ 1 l r 1 > η 1 . (8)

其中 s 1 ( η 1 ) = r 2 ( K 2 + a 2 η 1 ) exp ( r 2 τ 2 u ) r 2 *

我们断言:对这个 η 1 ,对任意 φ C + ,有 lim sup t x ( t , 0 , φ ) η 1

φ C + ,使得 lim sup t x ( t , 0 , φ ) < η 1 。记 x ( t ) = x ( t , 0 , φ ) , y ( t ) = y ( t , 0 , φ ) 。则存在 ε 0 > 0 ,使得当 ε ( 0 , ε 0 ) 时, T 3 > 0 ,当 T 3 > 0 时有(5)与(6)成立且 x ( t ) < η 1 + ε

由(1)得:

d y d t y ( t ) ( r 2 + ε ) ,

y ( t ) y ( t τ 2 ( t ) ) exp ( ( r 2 + ε ) τ 2 u ) . (9)

由(1),(9)得:当 t > T 时有

d y d t y ( t ) ( ( r 2 + ε ) ( r 2 * ε ) exp ( ( r 2 + ε ) τ 2 u ) y ( t ) ( K 2 + ε ) + ( a 2 + ε ) ( η 1 + ε ) ) .

由比较原理及 ε 的任意性可知:

lim sup t y ( t ) r 2 ( K 2 + a 2 η 1 ) exp ( r 2 τ 2 u ) r 2 * = s 1 ( η 1 ) . (10)

由(1),(10)得: T 1 > T ,当 t > T 1

d x d t x ( t ) ( ( r 1 * ε ) ( r 1 + ε ) M 1 ( K 1 * ε ) ( b 1 + ε ) ( s 1 ( η 1 ) + ε ) ) ,

x ( t ) x ( t τ 1 ( t ) ) exp ( ( r 1 * ε ) ( r 1 + ε ) ( M 1 + ε ) ( K 1 * ε ) ( b 1 + ε ) ( s 1 ( η 1 ) + ε ) ) τ 1 l . (11)

由(1),(11)得:

d x d t x ( t ) ( ( r 1 * ε ) h ( ε ) x ( t ) K 1 * ε ( b 1 + ε ) ( s 1 ( η 1 ) + ε ) ) ,

其中

h ( ε ) = ( r 1 + ε ) exp ( ( ( r 1 * ε ) ( r 1 + ε ) ( M 1 + ε ) K 1 * ε ( b 1 + ε ) ( s 1 ( η 1 ) + ε ) ) τ 1 l )

由比较原理及 ε 的任意性可知:

lim inf t x ( t ) ( r 1 * b 1 s 1 ( η 1 ) ) K 1 * exp ( r 1 * r 1 M 1 K 1 * b 1 s 1 ( η 1 ) ) τ 1 l r 1 > η 1 .

这与(8)及假设 lim sup t x ( t , 0 , φ ) < η 1 矛盾。从而 lim sup t x ( t , 0 , φ ) η 1 。同理存在 η 2 > 0 ,使得 φ C + , lim sup t y ( t , 0 , φ ) η 2

定理1 当(7)成立时,系统(1)的解是一致持久的。

证明:要证存在 m 1 > 0 使得 φ C +

lim inf t x ( t , 0 , φ ) > m 1 . (12)

由(7)知,存在 ε > 0 ,使得

( r 2 * ε ) ( r 1 + ε ) ε ( K 1 * ε ) ( b 2 + ε ) ( r 1 * ε ( r 1 + ε ) ( K 1 + ε ) exp ( ( r 1 + ε ) ( τ 1 u + ε ) ) ) 1 > 0 . (13)

对这个 ε > 0 T 1 使得当 t > T 1 时(5)与(6)成立。若(12)不成立,则存在 φ n C + 使得 lim inf t x ( t , 0 , φ n ) < 1 n 2 。存在 N 0 > 0 使得 1 N 0 < η 1 ( η 1 为引理4所给)。当 n > N 0 时,由引理4,存在 0 < τ 1 ( n ) < t 1 ( n ) < τ 2 ( n ) < < τ k ( n ) < t k ( n ) < 使得 lim k τ k ( n ) =

x ( τ k ( n ) , 0 , φ n ) = 1 n , x ( t k ( n ) , 0 , φ n ) = 1 n 2 . (14)

1 n 2 < x ( t , 0 , φ n ) < 1 n , τ k ( n ) < t < t k ( n ) 。由引理3,存在 T ( φ n ) > 0 使得当 t > T ( φ n ) x ( t , 0 , φ n ) < 2 M 1 , y ( t , 0 , φ n ) < 2 M 2 。从而当 t > T ( φ n ) + τ 时有

x ( t , 0 , φ n ) > x ( t , 0 , φ n ) ( 2 r 1 u M 1 K 1 l 2 b 1 u M 2 ) β 1 x ( t , 0 , φ n ) . (15)

存在 K 1 ( n ) ,当 k > K 1 ( n ) τ k ( n ) > max ( T ( φ n ) , T 1 ) 。由(14)与(15)知

t k ( n ) τ k ( n ) > ln n β 1 , k > K 1 ( n ) .

由(1)知 y ( t , 0 , φ n ) < ( r 2 * + ε ) y ( t , 0 , φ n ) , t T 1

从而 y ( t , 0 , φ n ) y ( t τ 2 ( t ) , 0 , φ n ) exp ( r 2 * + ε ) ( τ 2 * + ε ) 。故有 t ( τ k ( n ) , t k ( n ) )

y ( t , 0 , φ n ) y ( t , 0 , φ n ) ( ( r 2 * + ε ) ( r 2 * ε ) exp ( r 2 * + ε ) ( τ 2 * + ε ) ( K 2 * + ε ) + ( a 2 * + ε ) 1 n y ( t , 0 , φ n ) ) .

由此得

y ( t , 0 , φ n ) ( ( 1 y ( τ k ( n ) ) ( r 2 * ε ) exp ( r 2 * + ε ) ( τ 2 * + ε ) ( r 2 * + ε ) ( ( K 2 * + ε ) + ( a 2 * + ε ) 1 n ) ) exp ( r 2 * + ε ) ( t τ k ( n ) ) + ( r 2 * ε ) exp ( r 2 * + ε ) ( τ 2 * + ε ) ( r 2 * + ε ) ( ( K 2 * + ε ) + ( a 2 * + ε ) 1 n ) ) 1 (16)

k > K 1 ( n ) t τ k ( n ) > t k ( n ) τ k ( n ) 2 > ln n 2 β 1 时,由于

lim n ( 1 2 M 2 ( r 2 * ε ) exp ( r 2 + ε ) ( τ 2 + ε ) ( r 2 + ε ) ( ( K 2 + ε ) + ( a 2 + ε ) 1 n ) ) 1 n r 2 * + ε 2 β = 0 .

存在 N 1 > N 0 ,当 n > N 1 1 n < ε

| ( 1 2 M 2 ( r 2 * ε ) exp ( r 2 + ε ) ( τ 2 + ε ) ( r 2 + ε ) ( ( K 2 + ε ) + ( a 2 + ε ) 1 n ) ) 1 n r 2 * + ε 2 β | < ε . (17)

从而当 n > N 1 , k > K 1 ( n ) , t τ k ( n ) > t k ( n ) τ k ( n ) 2 时,由(16)与(17)知

y ( t , 0 , φ n ) < ( ( r 2 * ε ) exp ( r 2 + ε ) ( τ 2 + ε ) ( r 2 + ε ) ( ( K 2 + ε ) + ( a 2 + ε ) 1 n ) ε ) 1 .

从而由(13),当 n > N 1 k > K 1 ( n ) t ( τ k ( n ) > t k ( n ) τ k ( n ) 2 , t k ( n ) ) 时,

x ( t , 0 , φ n ) > x ( t , 0 , φ n ) ( r 1 * ε ( r 1 + ε ) ε K 1 * ε ( b 1 + ε ) ( ( r 2 * ε ) exp ( r 2 + ε ) ( τ 2 + ε ) ( r 2 + ε ) ( ( K 2 + ε ) + ( a 2 + ε ) ε ) ε ) 1 ) > 0

故有 x ( t k ( n ) , 0 , φ n ) > x ( τ k ( n ) + t k ( n ) τ k ( n ) 2 , 0 , φ n ) > 1 n 2 。矛盾。故(12)成立。同理,存在 m 2 > 0 使得 φ C + , lim inf t y ( t , 0 , φ ) > m 2

注:当 τ = 0 r i ( t ) , K i ( t ) , a i ( t ) , b i ( t ) , i = 1 , 2 为常数时,此时定理1的条件

r 1 * b 1 r 2 K 2 exp ( r 2 τ 2 u ) r 2 * > 0 , r 2 * b 2 r 1 K 1 exp ( r 1 τ 1 u ) r 1 * > 0

变为

b 1 K 2 r 1 < 0 , b 2 K 1 r 2 < 0 .

这与 [1] 中的一致持久的条件相同,因此,我们的结果是合理的。

下面我们给出一个例子来验证我们的理论结果。

例:考虑系统 x , y 的一致持久性:

{ d x d t = x ( t ) [ ( 0.17 + 0.01 sin t ) ( 0.25 + 0.01 sin t ) x ( t e 5 π + π sin t ) 0.9 + 0.1 cos t + ( 0.4 + 0.1 sin t ) y ( t ) ( 0.11 + 0.01 sin t ) y ( t ) ] d y d t = y ( t ) [ ( 0.22 + 0.01 sin t ) ( 0.22 + 0.01 sin t ) y ( t e 6 π + π cos t ) 1 + 0.1 cos t + ( 0.4 + 0.1 sin t ) x ( t ) ( 0.15 + 0.01 cos t ) x ( t ) ]

证明:由原方程组可得:

r 1 * = 0.18 , r 1 * = 0.16 , r 2 * = 0.23 , r 2 * = 0.21

b 1 * = 0.12 , b 2 * = 0.16 , K 1 * = 1 , K 2 * = 1.1 , τ 1 u = e 4 π , τ 2 u = e 5 π .

带入数值求得 r 1 * b 1 r 2 K 2 exp ( r 2 τ 2 u ) r 2 * > 0 r 2 * b 2 r 1 K 1 exp ( r 1 τ 1 u ) r 1 * > 0

根据定理1,从而原方程组是一致持久的。

参考文献

[1] Freedman, H.I. and Waltman, P. (1977) Mathematical Analysis of Some Three-Species Food-Chain Models. Mathe-matical Biosciences, 33, 257-276.
https://doi.org/10.1016/0025-5564(77)90142-0
[2] Gard, T.C. (1980) Persis-tence in Food-Chains with General Interactions. Mathematical Biosciences, 51, 165-174.
https://doi.org/10.1016/0025-5564(80)90096-6
[3] Garay, B.M. and Hofbauer, J. (2003) Robust Permanence for Ecological Differential Equations, Minimax, and Discretizations. SIAM Journal on Mathematical Analysis, 34, 1007-1039.
https://doi.org/10.1137/S0036141001392815
[4] Ashwin, P., Buescu, J. and Stewart, I. (1996) From Attractor to Chaotic Saddle: A Tale of Transverse Instability. Nonlinearity, 9, 703-737.
https://doi.org/10.1088/0951-7715/9/3/006
[5] Schreiber, S. (2000) Criteria for Cr Robust Permanence. Journal of Differential Equations, 162, 400-426.
https://doi.org/10.1006/jdeq.1999.3719
[6] Hirsch, M.W., Smith, H.L. and Zhao, X.-Q. (2001) Chain Transitivity, Attractivity and Strong Repellors for Semidynamical Systems. Journal of Dynamics and Differential Equations, 13, 107-131.
https://doi.org/10.1023/A:1009044515567
[7] Hofbauer, J. and Schreiber, S.J. (2010) Robust Permanence for Structured Populations. Journal of Differential Equations, 248, 1955-1971.
https://doi.org/10.1016/j.jde.2009.11.010
[8] Smith, H.L. and Zhao, X.-Q. (2001) Robust Persistence for Semidynamical Systems. Nonlinear Analysis, 47, 6169-6179.
https://doi.org/10.1016/S0362-546X(01)00678-2
[9] Thieme, H.R. (1993) Persistence under Relaxed Point-Dissipativity (with Application to an Endemic Model). SIAM Journal on Mathematical Analysis, 24, 407-435.
https://doi.org/10.1137/0524026
[10] Thieme, H.R. (2000) Uniform Persistence and Permanence for Non-Autonomous Semiflows in Population Biology. Mathematician Bios, 166, 173-201.
https://doi.org/10.1016/S0025-5564(00)00018-3
[11] Messan, M.R., Kopp, D., Allen, D.C. and Kang, Y. (2018) Dynamics Implications of Bi-Directional Exchange within a Meta-Ecosystem. Mathematical Biosciences, 301, 167-184.
https://doi.org/10.1016/j.mbs.2018.05.006
[12] Hale, J.K. (1977) Theory of Functional Differential Equations. Springer-Verlag, Berlin.
https://doi.org/10.1007/978-1-4612-9892-2
[13] Kuang, Y. (1993) Delay Differential Equations with Applications in Population Dynamics. Academic, New York.
[14] Yang, X.T. (2006) Uniform Persistence for a Predator-Prey System with Delays. Applied Mathematics and Computation, 173, 523-534.
https://doi.org/10.1016/j.amc.2005.04.057
[15] 蒋燕, 杨喜陶. 一类非自治密度依赖的时滞捕食者–食饵系统的一致持久性[J]. 湖南文理学院学报(自然科学版), 2017(29): 9-12.