移动环境下三种群Lotka-Volterra竞争合作系统行波解的存在性
The Existence of Traveling Wave Solutions for Three Species Lotka-Volterra Competitive-Cooperative System under a Shifting Habitat
摘要: 考虑在移动环境下局部扩散三种群Lotka-Volterra竞争合作系统行波解的存在性,并假设此系统的内禀增长率函数恒大于某正常数。通过构造一对有序的上下解并利用单调迭代技巧和波动引理,证明了系统的非负受迫行波的存在性。
Abstract: We consider the existence of traveling wave solutions for Lotka-Volterra competitive-cooperative system with three-species under a shifting habitat, and assume that the intrinsic growth rate functions of this system are greater than the normal numbers. We prove the existence of non-negative forced traveling waves of the system by constructing a pair of upper and lower solutions and using monotonic iterative techniques and the fluctuation lemma.
文章引用:陈碧霞. 移动环境下三种群Lotka-Volterra竞争合作系统行波解的存在性[J]. 应用数学进展, 2024, 13(8): 3989-4000. https://doi.org/10.12677/aam.2024.138380

1. 引言

在生态种群动力学的研究中,如持久生存等问题都可以通过建立适当的反应扩散方程模型来进行研究[1]-[5]。自然界中普遍存在合作、竞争和捕食[6]-[8]等种间关系,并且多物种之间的种间关系作用产生的动力学受到了广大学者的关注。特别地,竞争与合作关系被学者广泛研究,是多种群生态种群动力学的研究基础。因此,研究多种群反应扩散竞争合作系统行波的存在性具有重要的生物学意义。与此同时,种间的相互作用会产生相互耦合的非线性项,这为我们的研究工作带来了困难和挑战。

为了分析三种群Lotka-Volterra竞争合作系统的动力学,众多学者进行了广泛的研究[9]-[14]。Chen等人[9]考虑了三种群Lotka-Volterra竞争扩散模型,证明了非平凡行波解的存在性与稳定性并进行了数值模拟。后来,Mimura等人[10]考虑了一个弱的外来竞争物种入侵一个两个物种强烈竞争的系统的情况,讨论了三种群竞争–扩散系统竞争者介导共存的可能性。除此之外,Yang等人[6]讨论了气候变化下局部扩散Lotka-Volterra合作系统受迫行波的存在性和渐近行为。特别地,Hsu等人[11]建立了离散扩散的三种群Lotka-Volterra竞争合作系统来研究行波的存在性和稳定性,并且假设增长函数都是正常数。受上述研究工作的启发,我们考虑了增长函数为恒正的连续非减函数而非恒为正常数的情形,即研究了如下移动环境下局部扩散三种群Lotka-Volterra竞争合作系统

{ u( t,x ) t = d 1 2 u( t,x ) x 2 +u( t,x )[ r 1 ( xct )u( t,x )+hv( t,x ) a 1 w( t,x ) ] v( t,x ) t = d 2 2 v( t,x ) x 2 +v( t,x )[ r 2 ( xct )v( t,x )+ku( t,x ) a 2 w( t,x ) ] w( t,x ) t = d 3 2 w( t,x ) x 2 +w( t,x )[ r 3 ( xct )w( t,x ) b 1 u( t,x ) b 2 v( t,x ) ] (1)

行波解的存在性,且假设如下条件成立:

(H1) d 1 d 2 d 3 a 1 < a 2 h,k, a i , b i ( 0,1 ),i=1,2

(H2) r i ( )> a i ( i=1,2 ) r 1 ( )( b 1 + b 2 k )+ r 2 ( )( b 2 + b 1 h )> r 3 ( )( 1hk )>0

(H3) 函数 r i ( ) 上非减有界且连续,并满足 0< r i ( )< r i ( + )=1( i=1,2,3 )

其中 c>0 x u( t,x ) v( t,x ) w( t,x ) 分别表示3个物种在时刻t位置x处的种群密度, d i >0( i=1,2,3 ) 表示种群扩散速率, a 1 , a 2 , b 1 , b 2 表示物种间的竞争速率, h,k 表示种间合作率。在系统(1)中,种群uv是相互合作的,而种群w与种群uv是竞争的。

首先将系统(1)的受迫行波解记作 ( ϕ 1 ( ξ ), ϕ 2 ( ξ ), ϕ 3 ( ξ ) ),ξ=xct ,代入系统(1)可得

{ c ϕ 1 ( ξ )= d 1 ϕ 1 ( ξ )+ ϕ 1 ( ξ )[ r 1 ( ξ ) ϕ 1 ( ξ )+h ϕ 2 ( ξ ) a 1 ϕ 3 ( ξ ) ], c ϕ 2 ( ξ )= d 2 ϕ 2 ( ξ )+ ϕ 2 ( ξ )[ r 2 ( ξ ) ϕ 2 ( ξ )+k ϕ 1 ( ξ ) a 2 ϕ 3 ( ξ ) ], c ϕ 3 ( ξ )= d 3 ϕ 3 ( ξ )+ ϕ 3 ( ξ )[ r 3 ( ξ ) ϕ 3 ( ξ ) b 1 ϕ 1 ( ξ ) b 2 ϕ 2 ( ξ ) ], (2)

另外,系统(2)对应的极限方程组分别为

{ c ϕ 1 ( ξ )= d 1 ϕ 1 ( ξ )+ ϕ 1 ( ξ )[ r 1 ( + ) ϕ 1 ( ξ )+h ϕ 2 ( ξ ) a 1 ϕ 3 ( ξ ) ] c ϕ 2 ( ξ )= d 2 ϕ 2 ( ξ )+ ϕ 2 ( ξ )[ r 2 ( + ) ϕ 2 ( ξ )+k ϕ 1 ( ξ ) a 2 ϕ 3 ( ξ ) ] c ϕ 3 ( ξ )= d 3 ϕ 3 ( ξ )+ ϕ 3 ( ξ )[ r 3 ( + ) ϕ 3 ( ξ ) b 1 ϕ 1 ( ξ ) b 2 ϕ 2 ( ξ ) ] (3)

{ c ϕ 1 ( ξ )= d 1 ϕ 1 ( ξ )+ ϕ 1 ( ξ )[ r 1 ( ) ϕ 1 ( ξ )+h ϕ 2 ( ξ ) a 1 ϕ 3 ( ξ ) ], c ϕ 2 ( ξ )= d 2 ϕ 2 ( ξ )+ ϕ 2 ( ξ )[ r 2 ( ) ϕ 2 ( ξ )+k ϕ 1 ( ξ ) a 2 ϕ 3 ( ξ ) ], c ϕ 3 ( ξ )= d 3 ϕ 3 ( ξ )+ ϕ 3 ( ξ )[ r 3 ( ) ϕ 3 ( ξ ) b 1 ϕ 1 ( ξ ) b 2 ϕ 2 ( ξ ) ], (4)

通过简单计算,我们发现系统(2)的极限系统存在12个平衡点:

E 1 =( 0,0,0 ), E 2 =( 1,0,0 ), E 3 =( 0,1,0 ), E 4 =( 0,0,1 ), E 5 =( 0, 1 a 2 1 a 2 b 2 , 1 b 2 1 a 2 b 2 ), E 6 =( 1 a 1 1 a 1 b 1 ,0, 1 b 1 1 a 1 b 1 ), E 7 =( 1+h 1hk , 1+k 1hk ,0 ):=( k 1 + , k 2 + ,0 ), E 8 =( u + D , v + D , w + D ), E 9 =( 0, r 2 ( ) a 2 r 3 ( ) 1 a 2 b 2 , r 3 ( ) b 2 r 2 ( ) 1 a 2 b 2 ), E 10 =( r 1 ( ) a 1 r 3 ( ) 1 a 1 b 1 ,0, r 3 ( ) b 1 r 1 ( ) 1 a 1 b 1 ), E 11 =( r 1 ( )+h r 2 ( ) 1hk , r 2 ( )+k r 1 ( ) 1hk ,0 ), E 12 =( u D , v D , w D ).

其中

{ u + =[ a 1 a 1 b 2 + a 2 b 2 + a 2 hh1 ], v + =[ a 2 a 2 b 1 + a 1 b 1 + a 1 kk1 ], w + =[ b 1 + b 2 + b 1 h+ b 2 k+hk1 ],

{ u =[ r 1 ( )( 1 a 2 b 2 ) r 2 ( )( h+ a 1 b 2 )+ r 3 ( )( a 1 + a 2 h ) ], v =[ r 2 ( )( 1 a 1 b 1 ) r 1 ( )( k+ a 2 b 1 )+ r 3 ( )( a 2 + a 1 k ) ], w =[ r 3 ( )( 1hk )+ r 1 ( )( b 1 + b 2 k )+ r 2 ( )( b 2 + b 1 h ) ],

D= a 1 b 1 + a 2 b 2 + a 2 b 1 h+ a 1 b 2 k+hk1.

由于种群uv是相互合作的,而种群w与种群uv是竞争的,所以我们只考虑连接平衡点

E 11 ( r 1 ( )+h r 2 ( ) 1hk , r 2 ( )+k r 1 ( ) 1hk ,0 ) E 4 ( 0,0,1 ) 的行波解。

我们将 ( ϕ 1 , ϕ 2 , ϕ 3 ) 替换为 ( ϕ 1 , ϕ 2 ,1 ϕ 3 ) ,系统(2)转变为如下合作系统

{ c ϕ 1 ( ξ )= d 1 ϕ 1 ( ξ )+ ϕ 1 ( ξ )[ r 1 ( ξ ) a 1 ϕ 1 ( ξ )+h ϕ 2 ( ξ )+ a 1 ϕ 3 ( ξ ) ], c ϕ 2 ( ξ )= d 2 ϕ 2 ( ξ )+ ϕ 2 ( ξ )[ r 2 ( ξ ) a 2 ϕ 2 ( ξ )+k ϕ 1 ( ξ )+ a 2 ϕ 3 ( ξ ) ], c ϕ 3 ( ξ )= d 3 ϕ 3 ( ξ )+( 1 ϕ 3 ( ξ ) )[ 1 r 3 ( ξ ) ϕ 3 ( ξ )+ b 1 ϕ 1 ( ξ )+ b 2 ϕ 2 ( ξ ) ], (5)

此时 E 11 ( r 1 ( )+h r 2 ( ) 1hk , r 2 ( )+k r 1 ( ) 1hk ,0 ) 变为 E ( r 1 ( )+h r 2 ( ) 1hk , r 2 ( )+k r 1 ( ) 1hk ,1 ):=( k 1 , k 2 , k 3 ) E 4 ( 0,0,1 ) 变为 E + =( 0,0,0 ) 。此时边界条件转变为

lim ξ ( ϕ 1 ( ξ ), ϕ 2 ( ξ ), ϕ 3 ( ξ ) )=( k 1 , k 2 , k 3 ), lim ξ+ ( ϕ 1 ( ξ ), ϕ 2 ( ξ ), ϕ 3 ( ξ ) )=( 0,0,0 ), (6)

我们将在一定条件下寻找系统(5)满足边界条件(6)的非负解,为此,我们还需要如下技术性假设条件:

(H*) 存在常数 ρ>0 使得当 ξ+ r( + )r( ξ )=o( e ρξ ) 成立。

本文的剩余部分安排如下:在第2节中我们将给出一些预备知识,定义算子并验证它的一些性质,同时构造一对恰当的上下解。最后,在第3节中我们将利用单调迭代技巧结合波动引理证明非负受迫行波的存在性。

2. 预备知识

首先,我们引入一些函数空间。设空间 C( , ) 上所有连续函数组成, C + 表示由所有非负连续函数组成的空间,记

对任意的 u,vC( , ) ,如果 uv C + ,我们记 uv vu 。进一步,对任意的 u=( u 1 , u 2 , u 3 ) v=( v 1 , v 2 , v 3 )C( , )×C( , )×C( , ) ,如果 u 1 v 1 u 2 v 2 u 3 v 3 ,则记 uv vu

β i = a i +2 k i + r i ( )( i=1,2 ) β 3 =2+ b 1 k 1 + + b 2 k 2 + r 3 ( ) ,则方程 d i λ 2 cλ+ β i =0( i=1,2,3 ) 分别有两个实根:

λ i = c c 2 +4 d i β i 2 d i <0,   λ i+ = c+ c 2 +4 d i β i 2 d i >0.

定义二阶微分算子 Δ i 和它们的逆 Δ i 1 分别为

Δ i h:= d i h c h + β i h,

( Δ i 1 h )( ξ ):= 1 d i ( λ i+ λ i ) [ ξ e λ i ( ξη ) h( η )dη+ ξ + e λ i+ ( ξη ) h( η )dη ].

不难验证对任给连续有界函数h都有 Δ i ( Δ i 1 h )=h 。此外,若 h 也是连续有界函数,那么 Δ i 1 ( Δ i h )=h

定义2.1 ( ϕ ¯ 1 ( ξ ), ϕ ¯ 2 ( ξ ), ϕ ¯ 3 ( ξ ) ),( ϕ _ 1 ( ξ ), ϕ _ 2 ( ξ ), ϕ _ 3 ( ξ ) )BC×BC×BC 满足 ( ϕ ¯ 1 , ϕ ¯ 2 , ϕ ¯ 3 ) ( ϕ _ 1 , ϕ _ 2 , ϕ _ 3 ) ( ϕ ¯ 1 , ϕ ¯ 2 , ϕ ¯ 3 ) ( ϕ _ 1 , ϕ _ 2 , ϕ _ 3 ) L × L × L ϕ ¯ i ( ξ ) ϕ _ i ( ξ )( i=1,2,3 ) ( ϕ ¯ 1 , ϕ ¯ 2 , ϕ ¯ 3 ) ( ϕ _ 1 , ϕ _ 2 , ϕ _ 3 ) \{ ξ j } 上连续( ξ j 为一有限递增点列),且有 ϕ ¯ i ( ξ j + ) ϕ ¯ i ( ξ j ) ϕ _ i ( ξ j + ) ϕ _ i ( ξ j ) ,此时,若有不等式组

{ c ϕ ¯ 1 ( ξ ) d 1 ϕ ¯ 1 ( ξ )+ ϕ ¯ 1 ( ξ )[ r 1 ( ξ ) a 1 ϕ ¯ 1 ( ξ )+h ϕ ¯ 2 ( ξ )+ a 1 ϕ ¯ 3 ( ξ ) ] c ϕ ¯ 2 ( ξ ) d 2 ϕ ¯ 2 ( ξ )+ ϕ ¯ 2 ( ξ )[ r 2 ( ξ ) a 2 ϕ ¯ 2 ( ξ )+k ϕ ¯ 1 ( ξ )+ a 2 ϕ ¯ 3 ( ξ ) ] c ϕ ¯ 3 ( ξ ) d 3 ϕ ¯ 3 ( ξ )+( 1 ϕ ¯ 3 ( ξ ) )[ 1 r 3 ( ξ ) ϕ ¯ 3 ( ξ )+ b 1 ϕ ¯ 1 ( ξ )+ b 2 ϕ ¯ 2 ( ξ ) ] (7)

{ c ϕ _ 1 ( ξ ) d 1 ϕ _ 1 ( ξ )+ ϕ _ 1 ( ξ )[ r 1 ( ξ ) a 1 ϕ _ 1 ( ξ )+h ϕ _ 2 ( ξ )+ a 1 ϕ _ 3 ( ξ ) ] c ϕ _ 2 ( ξ ) d 2 ϕ _ 2 ( ξ )+ ϕ _ 2 ( ξ )[ r 2 ( ξ ) a 2 ϕ _ 2 ( ξ )+k ϕ _ 1 ( ξ )+ a 2 ϕ _ 3 ( ξ ) ] c ϕ _ 3 ( ξ ) d 3 ϕ _ 3 ( ξ )+( 1 ϕ _ 3 ( ξ ) )[ 1 r 3 ( ξ ) ϕ _ 3 ( ξ )+ b 1 ϕ _ 1 ( ξ )+ b 2 ϕ _ 2 ( ξ ) ] (8)

\{ ξ j } 上成立,则称 ( ϕ ¯ 1 ( ξ ), ϕ ¯ 2 ( ξ ), ϕ ¯ 3 ( ξ ) ) ( ϕ _ 1 ( ξ ), ϕ _ 2 ( ξ ), ϕ _ 3 ( ξ ) ) 为系统(5)的一对有序上下解。

给定一对有序上下解,可以构造先验集 Γ

Γ={ ( ϕ 1 , ϕ 2 , ϕ 3 )| ϕ i BC( , ), ϕ _ i ϕ i ϕ ¯ i ,i=1,2,3 }.

对任意 ( ϕ 1 , ϕ 2 , ϕ 3 )Γ ,定义如下算子 H=( H 1 , H 2 , H 3 )

H 1 ( ϕ 1 , ϕ 2 , ϕ 3 )( ξ ):= β 1 ϕ 1 ( ξ )+ ϕ 1 ( ξ )[ r 1 ( ξ ) a 1 ϕ 1 ( ξ )+h ϕ 2 ( ξ )+ a 1 ϕ 3 ( ξ ) ], H 2 ( ϕ 1 , ϕ 2 , ϕ 3 )( ξ ):= β 2 ϕ 2 ( ξ )+ ϕ 2 ( ξ )[ r 2 ( ξ ) a 2 ϕ 2 ( ξ )+k ϕ 1 ( ξ )+ a 2 ϕ 3 ( ξ ) ], H 3 ( ϕ 1 , ϕ 2 , ϕ 3 )( ξ ):= β 3 ϕ 3 ( ξ )+( 1 ϕ 3 ( ξ ) )[ 1 r 3 ( ξ ) ϕ 3 ( ξ )+ b 1 ϕ 1 ( ξ )+ b 2 ϕ 2 ( ξ ) ].

接下来,我们定义映射 F=( F 1 , F 2 , F 3 ) ,其中 F i ( ϕ 1 , ϕ 2 , ϕ 3 )= Δ i 1 H i ( ϕ 1 , ϕ 2 , ϕ 3 )( i=1,2,3 ) Δ 1 =( Δ 1 1 , Δ 2 1 , Δ 3 1 ) ,对所有的 ( ϕ 1 , ϕ 2 , ϕ 3 )Γ 。下面先讨论映射F的一些性质。

引理2.1 F是一个非减算子,且 F( Γ )Γ

证明 一方面,取 ( ϕ ˜ 1 , ϕ ˜ 2 , ϕ ˜ 3 ),( ϕ ^ 1 , ϕ ^ 2 , ϕ ^ 3 )Γ 且满足 ( ϕ ˜ 1 , ϕ ˜ 2 , ϕ ˜ 3 )( ϕ ^ 1 , ϕ ^ 2 , ϕ ^ 3 ) ,则对任何 ξ

H 1 ( ϕ ˜ 1 , ϕ ˜ 2 , ϕ ˜ 3 )( ξ ) H 1 ( ϕ ^ 1 , ϕ ^ 2 , ϕ ^ 3 ) =( ϕ ˜ 1 ( ξ ) ϕ ^ 1 ( ξ ) )[ β 1 + r 1 ( ξ ) a 1 ϕ ˜ 1 ( ξ ) ϕ ^ 1 ( ξ )+h ϕ ˜ 2 ( ξ )+ a 1 ϕ ˜ 3 ( ξ ) ]   +h( ϕ ˜ 2 ( ξ ) ϕ ^ 2 ( ξ ) ) ϕ ^ 1 ( ξ )+ a 1 ( ϕ ˜ 3 ( ξ ) ϕ ^ 3 ( ξ ) ) ϕ ^ 1 ( ξ )0.

同理可得

H 2 ( ϕ ˜ 1 , ϕ ˜ 2 , ϕ ˜ 3 )( ξ ) H 2 ( ϕ ^ 1 , ϕ ^ 2 , ϕ ^ 3 )0

H 3 ( ϕ ˜ 1 , ϕ ˜ 2 , ϕ ˜ 3 )( ξ ) H 3 ( ϕ ^ 1 , ϕ ^ 2 , ϕ ^ 3 )0.

H( ϕ ˜ 1 , ϕ ˜ 2 , ϕ ˜ 3 )H( ϕ ^ 1 , ϕ ^ 2 , ϕ ^ 3 ) 。从而,对任意的 ξ ,有 F( ϕ ˜ 1 , ϕ ˜ 2 , ϕ ˜ 3 )( ξ )F( ϕ ^ 1 , ϕ ^ 2 , ϕ ^ 3 )( ξ ) 。即F是一个非减算子。

另一方面,由 F( Γ ) 的定义可知,我们只需证对所有 ( ϕ 1 , ϕ 2 , ϕ 3 )Γ 都有

( ϕ _ 1 , ϕ _ 2 , ϕ _ 3 )F( ϕ 1 , ϕ 2 , ϕ 3 )( ϕ ¯ 1 , ϕ ¯ 2 , ϕ ¯ 3 ).

由于 ( ϕ ¯ 1 ( ξ ), ϕ ¯ 2 ( ξ ), ϕ ¯ 3 ( ξ ) ) ( ϕ _ 1 ( ξ ), ϕ _ 2 ( ξ ), ϕ _ 3 ( ξ ) ) 是一对有序上下解,结合文献[15]的引理3.2我们有

F i ( ϕ _ 1 , ϕ _ 2 , ϕ _ 3 )= Δ i 1 H i ( ϕ _ 1 , ϕ _ 2 , ϕ _ 3 ) Δ i 1 Δ i ( ϕ _ 1 , ϕ _ 2 , ϕ _ 3 )= ϕ _ i ,i=1,2,3,

F i ( ϕ ¯ 1 , ϕ ¯ 2 , ϕ ¯ 3 )= Δ i 1 H i ( ϕ ¯ 1 , ϕ ¯ 2 , ϕ ¯ 3 ) Δ i 1 Δ i ( ϕ ¯ 1 , ϕ ¯ 2 , ϕ ¯ 3 )= ϕ ¯ i ,i=1,2,3,

F( ϕ _ 1 , ϕ _ 2 , ϕ _ 3 )( ϕ _ 1 , ϕ _ 2 , ϕ _ 3 ),  F( ϕ ¯ 1 , ϕ ¯ 2 , ϕ ¯ 3 )( ϕ ¯ 1 , ϕ ¯ 2 , ϕ ¯ 3 ). (9)

又由F是一个非减算子,所以对任意的 ( ϕ 1 , ϕ 2 , ϕ 3 )Γ

F( ϕ _ 1 , ϕ _ 2 , ϕ _ 3 )F( ϕ 1 , ϕ 2 , ϕ 3 )F( ϕ ¯ 1 , ϕ ¯ 2 , ϕ ¯ 3 ). (10)

结合(9)和(10)得到 F( Γ )Γ 。证毕。

系统(5)可写为

Δ i ( ϕ 1 , ϕ 2 , ϕ 3 )= H i ( ϕ 1 , ϕ 2 , ϕ 3 ), i=1,2,3. (11)

因此,若映射F Γ 中存在一个不动点,即存在 ( ϕ 1 , ϕ 2 , ϕ 3 )Γ 使得

( ϕ 1 , ϕ 2 , ϕ 3 )=F( ϕ 1 , ϕ 2 , ϕ 3 ),

则该不动点必是(11)的解。若该不动点还满足边界条件(6),则必是系统(5)的受迫行波。这就是我们的研究目标。为此,最后我们通过选取一对恰当的有序上下解构造出一个先验集。

对于任意给定的 c>0 ,定义连续函数

Φ i ( c,λ ):= d i λ 2 cλ+1 a i =0,λ,i=1,2.

显然 Φ i ( c,λ ),i=1,2 有以下性质

(i) Φ i ( c,0 )>0

(ii) 对任意的 c>0 ,都有 lim λ Φ i ( c,λ )=+

(iii) 对任意的 λ>0 ,都有 lim c Φ i ( c,λ )=

(iv) Φ i ( c,λ ) c =λ<0 2 Φ i ( c,λ ) λ 2 =2 d i >0

根据这些性质可以得到如下引理

引理2.2

c i * ( ):=2 d i ( 1 a i ) ,i=1,2,

若(H1)成立,则有

(i) 当 c= c i * 时, Φ i ( c,λ )=0 有唯一正根 λ i 0

(ii) 当 c> c i * 时, Φ i ( c,λ )=0 分别有两个不同正根 λ i ± ,并满足:当 λ( λ i , λ i + ) 时,有 Φ i ( c,λ )<0 ;当 λ[ 0, λ i )( λ i + ,+ ) 时,有 Φ i ( c,λ )>0

(iii) 当 c< c i * 时,对于任意的 λ Φ i ( c,λ )>0 。此外,当 c> c 1 * 时,我们有 λ 2 λ 1 < λ 1 + λ 2 +

c> c 1 * 时,存在 η>0 使得 λ 1 <η λ 1 <min{ λ 1 + , λ 2 + , λ 1 + λ 2 } 。那么对任意 c> c 1 * ( ) Φ i ( c,η λ 1 )<0,i=1,2 。因此存在充分大的正数 ξ 1 ξ 2 分别满足

( a 1 +h ) e ( λ 1 + λ 2 η λ 1 ) ξ 1 +( a 1 +h ) k 1 + e λ 2 ξ 1 + k 1 + Φ 1 ( c,η λ 1 )=0

k e ( λ 1 + λ 2 η λ 1 ) ξ 2 +k k 2 + e λ 1 ξ 2 + k 2 + Φ 2 ( c,η λ 1 )=0.

由此,定义如下有界连续函数

ϕ ¯ 1 ( ξ )={ e λ 1 ξ +q k 1 + e η λ 1 ξ ,ξ ξ 1 , k 1 + ,                     ξ< ξ 1 ,     ϕ ¯ i ( ξ )={ e λ 2 ξ +q k i + e η λ 1 ξ ,ξ ξ i , k i + ,                     ξ< ξ i ,   i=2,3.

其中 k 1 + = 1+h 1hk k 2 + = 1+k 1hk k 3 + = k 3 =1

引理2.3 对任意 c> c 1 * ,当 q>1 足够大时, ( ϕ ¯ 1 , ϕ ¯ 2 , ϕ ¯ 3 )( ξ ) 是系统(5)的一个上解。

证明 方便起见,我们记

A 1 ( ξ ):= d 1 ϕ ¯ 1 ( ξ )+c ϕ ¯ 1 + ϕ ¯ 1 ( ξ )[ r 1 ( ξ ) a 1 ϕ ¯ 1 ( ξ )+h ϕ ¯ 2 ( ξ )+ a 1 ϕ ¯ 3 ( ξ ) ], A 2 ( ξ ):= d 2 ϕ ¯ 2 ( ξ )+c ϕ ¯ 2 + ϕ ¯ 2 ( ξ )[ r 2 ( ξ ) a 2 ϕ ¯ 2 ( ξ )+k ϕ ¯ 1 ( ξ )+ a 2 ϕ ¯ 3 ( ξ ) ], A 3 ( ξ ):= d 3 ϕ ¯ 3 ( ξ )+c ϕ ¯ 3 +( 1 ϕ ¯ 3 ( ξ ) )[ 1 r 3 ( ξ ) ϕ ¯ 3 ( ξ )+ b 1 ϕ ¯ 1 ( ξ )+ b 2 ϕ ¯ 2 ( ξ ) ].

要证明 ( ϕ ¯ 1 ( ξ ), ϕ ¯ 2 ( ξ ), ϕ ¯ 3 ( ξ ) ) 是系统(5)的一个上解,只需证明 A i ( ξ )0( i=1,2,3 ) 即可。下面先证明 A 1 ( ξ )0

(i) 当 ξ< ξ 1 时, ϕ ¯ 1 ( ξ )= k 1 + ϕ ¯ i ( ξ ) k i + ,i=2,3 ,则有

A 1 ( ξ ) k 1 + ( 1 a 1 k 1 + +h k 2 + + a 1 k 3 + )=0.

(ii) 当 ξ ξ 1 时, ϕ ¯ 1 ( ξ )= e λ 1 ξ +q k 1 + e η λ 1 ξ ϕ ¯ i ( ξ ) e λ 2 ξ +q k i + e η λ 1 ξ ,i=2,3 。由假设(H1)可知 a 1 <1 ,那么 k 1 + +h k 2 + + a 1 k 3 + <0 ,因此

A 1 ( ξ )= d 1 ( λ 1 ) 2 e λ 1 ξ c λ 1 e λ 1 ξ +( e λ 1 ξ +q k 1 + e η λ 1 ξ )[ r 1 ( ξ ) a 1 ]             + d 1 q k 1 + ( η λ 1 ) 2 e η λ 1 ξ cq k 1 + η λ 1 e η λ 1 ξ + ϕ ¯ 1 ( ξ )[ ϕ ¯ 1 ( ξ )+h ϕ ¯ 2 ( ξ )+ a 1 ϕ ¯ 3 ]( ξ )           e λ 1 ξ Φ 1 ( c, λ 1 )+q k 1 + e η λ 1 ξ Φ 1 ( c,η λ 1 )+ ϕ ¯ 1 ( ξ )[ ϕ ¯ 1 ( ξ )+h ϕ ¯ 2 ( ξ )+ a 1 ϕ ¯ 3 ( ξ ) ]          q k 1 + e η λ 1 ξ Φ 1 ( c,η λ 1 )+ ϕ ¯ 1 ( ξ )[ e λ 1 ξ +( a 1 +h ) e λ 2 ξ +q e η λ 1 ξ ( k 1 + +h k 2 + + a 1 k 3 + ) ]          q k 1 + e η λ 1 ξ Φ 1 ( c,η λ 1 )+( e λ 1 ξ +q k 1 + e η λ 1 ξ )( a 1 +h ) e λ 2 ξ           e η λ 1 ξ [ q k 1 + Φ 1 ( c,η λ 1 )+q( a 1 +h ) e ( λ 1 + λ 2 η λ 1 ) ξ 1 +q( a 1 +h ) k 1 + e λ 2 ξ 1 ]=0.

接下来证明 A 2 ( ξ )0 .

(i) 当 ξ< ξ 2 时, ϕ ¯ 2 ( ξ )= k 2 + ϕ ¯ i ( ξ ) k i + ,i=1,3 ,则有

A 2 ( ξ ) k 2 + ( 1 a 2 k 2 + +k k 1 + + a 2 k 3 + )=0.

(ii) 当 ξ ξ 2 时, ϕ ¯ 2 ( ξ )= e λ 2 ξ +q k 2 + e η λ 1 ξ ϕ ¯ 1 ( ξ ) e λ 1 ξ +q k 1 + e η λ 1 ξ ϕ ¯ 3 ( ξ ) e λ 2 ξ +q k 3 + e η λ 1 ξ 。由假设(H1)可知 a 2 <1 ,那么 k 2 + +k k 1 + + a 2 k 3 + <0 ,因此,类似 ξ ξ 1 A 1 ( ξ )0 的证明过程可得

A 2 ( ξ )0.

最后证明 A 3 ( ξ )0 。由于 ϕ i ϕ ¯ i k i + 0< r 3 ( ξ )1 ,其中 i=1,2,3 ,因此存在正常数 M>0 使得

( 1 ϕ ¯ 3 ( ξ ) )[ 1 r 3 ( ξ ) ϕ ¯ 3 ( ξ )+ b 1 ϕ ¯ 1 ( ξ )+ b 2 ϕ ¯ 2 ( ξ ) ]M.

(i) 当 ξ< ξ 3 时, ϕ ¯ 3 ( ξ )= k 3 + =1 ,则有

A 3 ( ξ )=( 1 k 3 + )[ 1 r 3 ( ξ ) k 3 + + b 1 ϕ ¯ 1 ( ξ )+ b 2 ϕ ¯ 2 ( ξ ) ]=0.

(ii) 当 ξ ξ 3 时, ϕ ¯ 3 ( ξ )= e λ 2 ξ +q k 3 + e η λ 1 ξ 。由假设(H1)可知 d 3 d 2 ,那么

A 3 ( ξ )= d 3 ϕ ¯ 3 ( ξ )+c ϕ ¯ 3 +( 1 ϕ ¯ 3 ( ξ ) )[ 1 r 3 ( ξ ) ϕ ¯ 3 ( ξ )+ b 1 ϕ ¯ 1 ( ξ )+ b 2 ϕ ¯ 2 ( ξ ) ]           d 2 ( λ 2 ) 2 e λ 2 ξ c λ 2 e λ 2 ξ + d 2 q k 3 + ( η λ 1 ) 2 e η λ 1 ξ cq k 3 + η λ 1 e η λ 1 ξ +M          = e λ 2 ξ Φ 2 ( c, λ 2 )+q k 3 + e η λ 1 ξ ( d 2 ( η λ 1 ) 2 cη λ 1 )( 1 a 2 ) e λ 2 ξ +M0.

q>1 足够大时最后一个不等式成立。

综上所述,对任意的 ξ ,都有 A i ( ξ )0,i=1,2,3 。因此, ( ϕ ¯ 1 , ϕ ¯ 2 , ϕ ¯ 3 )( ξ ) 是系统(5)的一个上解,证毕。

接下来我们构造系统(5)的下解。由假设条件(H*)可知,存在 K,C>0 使得当 ξK

1 r i ( ξ )C e ρξ ,i=1,2 成立。定义连续函数

ϕ _ i ( ξ )={ e λ i ξ q k i + e η λ i ξ ,ξ ξ i * , r i ( ) a i ,        ξ< ξ i * ,   i=1,2,      ϕ _ 3 ( ξ )0.

其中q η 都是正常数,并满足

1<η<min{ λ 1 + / λ 1 , λ 2 + / λ 2 ,2,1+ρ/ λ 1 ,1+ρ/ λ 2 }.

引理2.4 对任意 c> c 1 * ,当 q>1 足够大时, ( ϕ _ 1 , ϕ _ 2 , ϕ _ 3 )( ξ ) 是系统(5)的一个下解。

证明 不妨记

B 1 ( ξ ):= d 1 ϕ _ 1 ( ξ )+c ϕ _ 1 + ϕ _ 1 ( ξ )[ r 1 ( ξ ) a 1 ϕ _ 1 ( ξ )+h ϕ _ 2 ( ξ )+ a 1 ϕ _ 3 ( ξ ) ], B 2 ( ξ ):= d 2 ϕ _ 2 ( ξ )+c ϕ _ 2 + ϕ _ 2 ( ξ )[ r 2 ( ξ ) a 2 ϕ _ 2 ( ξ )+k ϕ _ 1 ( ξ )+ a 2 ϕ _ 3 ( ξ ) ], B 3 ( ξ ):= d 3 ϕ _ 3 ( ξ )+c ϕ _ 3 +( 1 ϕ _ 3 ( ξ ) )[ 1 r 3 ( ξ ) ϕ _ 3 ( ξ )+ b 1 ϕ _ 1 ( ξ )+ b 2 ϕ _ 2 ( ξ ) ].

要证 ( ϕ _ 1 , ϕ _ 2 , ϕ _ 3 )( ξ ) 是系统(5)的下解,需证 B i ( ξ )0( i=1,2,3 )

下面先证明 B 1 ( ξ )0

(i) 当 ξ< ξ 1 * 时, ϕ _ 1 ( ξ )= r 1 ( ) a 1 >0 ,则有

B 1 ( ξ )= ϕ _ 1 ( ξ )[ r 1 ( ξ ) r 1 ( )+h ϕ _ 2 ( ξ ) ]>0.

(ii) 当 ξ ξ 1 * 时, ϕ _ 1 ( ξ )= e λ 1 ξ q k 1 + e η λ 1 ξ e λ 1 ξ ϕ _ 3 ( ξ )=0 。那么

B 1 ( ξ )= d 1 ( λ 1 ) 2 e λ 1 ξ c λ 1 e λ 1 ξ +( e λ 1 ξ q k 1 + e η λ 1 ξ )[ 1 a 1 ]              d 1 q k 1 + ( η λ 1 ) 2 e η λ 1 ξ +cq k 1 + η λ 1 e η λ 1 ξ + ϕ _ 1 ( ξ )[ r 1 ( ξ )1 ϕ _ 1 ( ξ )+h ϕ _ 2 ( ξ ) ]          = e λ 1 ξ Φ 1 ( c, λ 1 )q k 1 + e η λ 1 ξ Φ 1 ( c,η λ 1 )+ ϕ _ 1 ( ξ )[ r 1 ( ξ )1 ϕ _ 1 ( ξ )+h ϕ _ 2 ( ξ ) ]          q k 1 + e η λ 1 ξ Φ 1 ( c,η λ 1 )+( e λ 1 ξ q k 1 + e η λ 1 ξ )[ C e ρξ e λ 1 ξ ]           e η λ 1 ξ [ q k 1 + Φ 1 ( c,η λ 1 )( C e ( λ 1 +ρη λ 1 )ξ + e ( 2η ) λ 1 ξ ) ]           e η λ 1 ξ [ q k 1 + Φ 1 ( c,η λ 1 )( C e ( λ 1 +ρη λ 1 ) ξ 1 * + e ( 2η ) λ 1 ξ 1 * ) ]0.

q>1 足够大时最后一个不等号成立。接下来证明 B 2 ( ξ )0

(i) 当 ξ< ξ 2 * 时, ϕ _ 2 ( ξ )= r 2 ( ) a 2 >0 ,则有

B 2 ( ξ )= ϕ _ 2 ( ξ )[ r 2 ( ξ ) r 2 ( )+k ϕ _ 1 ( ξ ) ]>0.

(ii) 当 ξ ξ 2 * 时, ϕ _ 2 ( ξ )= e λ 2 ξ q k 2 + e η λ 2 ξ e λ 2 ξ ϕ _ 3 ( ξ )=0 ,那么,与 ξ ξ 1 * B 1 ( ξ )0 的证明过程类似,当 q>1 足够大时有

B 2 ( ξ )0.

最后对于任意 ξ ,都有 ϕ _ 3 ( ξ )0 ,因此 B 3 ( ξ )=0

综上所述,当q足够大时,对任意的 ξ ,都有 B i ( ξ )0( i=1,2,3 ) ,又由 ( ϕ ¯ 1 , ϕ ¯ 2 , ϕ ¯ 3 )( ξ ) ( ϕ _ 1 , ϕ _ 2 , ϕ _ 3 )( ξ ) 的构造易知 ( ϕ ¯ 1 , ϕ ¯ 2 , ϕ ¯ 3 )( ξ )( ϕ _ 1 , ϕ _ 2 , ϕ _ 3 )( ξ ) 。因此, ( ϕ _ 1 , ϕ _ 2 , ϕ _ 3 )( ξ ) 是系统(5)的一个下解,证毕。

由上述两个引理,我们得到如下先验集:

Γ:={ ( ϕ 1 , ϕ 2 , ϕ 3 )| ϕ 1 , ϕ 2 , ϕ 3 BC( , ), ϕ _ i ϕ i ϕ ¯ i ,i=1,2,3 }.

3. 行波解的存在性

定理3.1 若(H1)~(H3)以及(H*)成立,那么对任意 c> c 1 * ,当 q>1 足够大时,系统(5)总存在一个满足边界条件(6)的非负受迫行波。

证明 首先构造如下迭代序列:

( ϕ 1 ( 1 ) , ϕ 2 ( 1 ) , ϕ 3 ( 1 ) )=F( ϕ ¯ 1 , ϕ ¯ 2 , ϕ ¯ 3 ),  ( ϕ 1 ( n+1 ) , ϕ 2 ( n+1 ) , ϕ 3 ( n+1 ) )=F( ϕ 1 ( n ) , ϕ 2 ( n ) , ϕ 3 ( n ) ),n1.

由于 ϕ ¯ 1 ( ξ ), ϕ ¯ 2 ( ξ ), ϕ ¯ 3 ( ξ ) 上的有界连续函数,结合引理2.1得到对所有的 n1 ϕ 1 ( n ) ( ξ ) ϕ 2 ( n ) ( ξ ) ϕ 3 ( n ) ( ξ ) 也是 上的有界连续函数且满足不等式

( ϕ ¯ 1 , ϕ ¯ 2 , ϕ ¯ 3 )( ϕ 1 ( n ) , ϕ 2 ( n ) , ϕ 3 ( n ) )( ϕ 1 ( n+1 ) , ϕ 2 ( n+1 ) , ϕ 3 ( n+1 ) )( ϕ _ 1 , ϕ _ 2 , ϕ _ 3 ).

从而,存在一个有界连续函数 ( ϕ 1 ( ξ ), ϕ 2 ( ξ ), ϕ 3 ( ξ ) )Γ 使得

lim n ( ϕ 1 ( n ) ( ξ ), ϕ 2 ( n ) ( ξ ), ϕ 3 ( n ) ( ξ ) )=( ϕ 1 ( ξ ), ϕ 2 ( ξ ), ϕ 3 ( ξ ) ).

不难看出,对所有的 n1 ξ

| H 1 ( ϕ 1 ( n ) , ϕ 2 ( n ) , ϕ 3 ( n ) )( ξ ) | k 1 + ( β 1 +1+ k 1 + +h k 2 + ),

| H 2 ( ϕ 1 ( n ) , ϕ 2 ( n ) , ϕ 3 ( n ) )( ξ ) | k 2 + ( β 2 +1+ k 2 + +k k 1 + )

| H 3 ( ϕ 1 ( n ) , ϕ 2 ( n ) , ϕ 3 ( n ) ) | β 3 +3+ b 1 k 1 + + b 2 k 2 + .

从而利用Lebesgue’s控制收敛定理可得

ϕ 1 ( ξ )= lim n ϕ 1 ( n+1 ) ( ξ )= lim n F 1 ( ϕ 1 ( n ) , ϕ 2 ( n ) , ϕ 3 ( n ) )( ξ )         = lim n ( Δ 1 1 H 1 ( ϕ 1 ( n ) , ϕ 2 ( n ) , ϕ 3 ( n ) ) )( ξ )         = 1 d 1 ( λ 1+ λ 1 ) [ ξ e λ 1 ( ξη ) H 1 ( ϕ 1 , ϕ 2 , ϕ 3 )( η )dη+ ξ + e λ 1+ ( ξη ) H 1 ( ϕ 1 , ϕ 2 , ϕ 3 )( η )dη ]         = F 1 ( ϕ 1 , ϕ 2 , ϕ 3 )( ξ ).

同理可得

ϕ 2 ( ξ )= F 2 ( ϕ 1 , ϕ 2 , ϕ 3 )( ξ ),   ϕ 3 ( ξ )= F 3 ( ϕ 1 , ϕ 2 , ϕ 3 )( ξ ).

( ϕ 1 ( ξ ), ϕ 2 ( ξ ), ϕ 3 ( ξ ) )Γ 是算子F的不动点,也就是说 ( ϕ 1 ( ξ ), ϕ 2 ( ξ ), ϕ 3 ( ξ ) )Γ 是系统(5)的解。由 ( ϕ _ 1 ( ξ ), ϕ _ 2 ( ξ ), ϕ _ 3 ( ξ ) ) 的构造特点可知 ( ϕ 1 ( ξ ), ϕ 2 ( ξ ), ϕ 3 ( ξ ) ) 为非负解。

根据 ( ϕ _ 1 ( ξ ), ϕ _ 2 ( ξ ), ϕ _ 3 ( ξ ) ) ( ϕ ¯ 1 ( ξ ), ϕ ¯ 2 ( ξ ), ϕ ¯ 3 ( ξ ) ) 的定义知

lim ξ+ ( ϕ _ 1 ( ξ ), ϕ _ 2 ( ξ ), ϕ _ 3 ( ξ ) )=( 0,0,0 ), lim ξ+ ( ϕ _ 1 ( ξ ), ϕ _ 2 ( ξ ), ϕ _ 3 ( ξ ) )=( 0,0,0 ).

于是

lim ξ+ ( ϕ 1 ( ξ ), ϕ 2 ( ξ ), ϕ 3 ( ξ ) )=( 0,0,0 ).

下面我们证明 lim ξ ( ϕ 1 ( ξ ), ϕ 2 ( ξ ), ϕ 3 ( ξ ) )=( k 1 , k 2 , k 3 ) 。记

( P 1 , P 2 , P 3 ):= limsup ξ ( ϕ 1 ( ξ ), ϕ 2 ( ξ ), ϕ 3 ( ξ ) ),  ( Q 1 , Q 2 , Q 3 ):= liminf ξ ( ϕ 1 ( ξ ), ϕ 2 ( ξ ), ϕ 3 ( ξ ) ).

于是 0< Q i P i k i ( i=1,2 ) 0 Q 3 P 3 k 3 =1 。由波动引理(文献[16]中引理A.1),存在满足 lim n s n = 的单调序列 { s n } n=1 和满足 lim n t n = 的单调序列 { t n } n=1 使得

lim n ( ϕ 1 ( s n ), ϕ 2 ( s n ), ϕ 3 ( s n ) )=( P 1 , P 2 , P 3 ), lim n ( ϕ 1 ( t n ), ϕ 2 ( t n ), ϕ 3 ( t n ) )=( Q 1 , Q 2 , Q 3 ),

lim n ( ϕ 1 ( s n ), ϕ 2 ( s n ), ϕ 3 ( s n ) )= lim n ( ϕ 1 ( t n ), ϕ 2 ( t n ), ϕ 3 ( t n ) )=( 0,0,0 ).

根据 ( ϕ 1 , ϕ 2 , ϕ 3 )( ξ )=F( ϕ 1 , ϕ 2 , ϕ 3 )( ξ ) 可得

ϕ i ( ξ )= λ i+ ϕ i ( ξ ) 1 d i ξ e λ i ( ξη ) H i ( ϕ 1 , ϕ 2 , ϕ 3 )( η )dη.

对任意的 ϵ>0 ,存在 N 1 >0 使得当 η( , s N 1 ) 时有

0< ϕ i ( η )< P i +ϵ, r i ( )< r i ( η )< r i ( )+ϵ,i=1,2,

0 ϕ 3 ( η )< P 3 +ϵ, r 3 ( )< r 3 ( η )< r 3 ( )+ϵ.

因此,当 n> N 1 时有

ϕ 1 ( s n )= λ 1+ ϕ 1 ( s n ) 1 d 1 s n e λ 1 ( s n η ) H 1 ( ϕ 1 , ϕ 2 , ϕ 3 )( η )dη λ 1+ ϕ 1 ( s n ) 1 d 1 s n e λ 1 ( s n η ) ( P 1 +ϵ )[ β 1 + r 1 ( ) a 1 P 1 +h P 2 + a 1 P 3 +( a 1 +h )ϵ ]dη = λ 1+ ϕ 1 ( s n )+ 1 d 1 λ 1 ( P 1 +ϵ )[ β 1 + r 1 ( ) a 1 P 1 +h P 2 + a 1 P 3 +( a 1 +h )ϵ ].

n ,我们有

λ 1+ P 1 + 1 d 1 λ 1 ( P 1 +ϵ )[ β 1 + r 1 ( ) a 1 P 1 +h P 2 + a 1 P 3 +( a 1 +h )ϵ ]0.

ϵ 的任意性得

λ 1+ P 1 + 1 d 1 λ 1 P 1 ( β 1 + r 1 ( ) a 1 P 1 +h P 2 + a 1 P 3 )0.

注意到 P 1 >0 ,上述不等式可推得

r 1 ( ) a 1 P 1 +h P 2 + a 1 P 3 0. (12)

同理可得

r 2 ( ) a 2 P 2 +k P 1 + a 2 P 3 0 (13)

( 1 P 3 )( 1 r 3 ( ) P 3 + b 1 P 1 + b 2 P 2 )0. (14)

由于 0 P 3 k 3 =1 ,可分为 P 3 =1 0 P 3 <1 两种情况进行讨论。

P 3 =1 ,那么联立(12)和(13)可得

P 1 k 1 , P 2 k 2 . (15)

0 P 3 <1 ,由(14)可以推出

1 r 3 ( ) P 3 + b 1 P 1 + b 2 P 2 0. (16)

联立(12)、(13)和(16)可得

P 1 u D , P 2 v D , P 3 1 w D . (17)

类似地,对任意的 ϵ( 0,min{ Q 1 , Q 2 , Q 3 } ) ,存在 N 2 >0 使得当 η( , t N 2 ) 时有

Q i ϵ< ϕ i ( η )< Q i +ϵ, r i ( )< r i ( η )< r i ( )+ϵ,i=1,2,

Q 3 ϵ ϕ 3 ( η )< Q 3 +ϵ, r 3 ( )< r 3 ( η )< r 3 ( )+ϵ.

因此,当 n> N 2 时有

ϕ 1 ( t n ) λ 1+ ϕ 1 ( t n )+ 1 d 1 λ 1 ( Q 1 ϵ )[ β 1 + r 1 ( ) a 1 Q 1 +h Q 2 + a 1 Q 3 ( a 1 +h1 )ϵ ].

注意到 Q 1 >0 ,令 n 且由 ϵ 的任意性可得

r 1 ( ) a 1 Q 1 +h Q 2 + a 1 Q 3 0. (18)

同理可得

r 2 ( ) a 2 Q 2 +k Q 1 + a 2 Q 3 0 (19)

( 1 Q 3 )( 1 r 3 ( ) Q 3 + b 1 Q 1 + b 2 Q 2 )0. (20)

由于 0 Q 3 k 3 =1 ,可分为 Q 3 =1 0 Q 3 <1 两种情况进行讨论。

Q 3 =1 ,那么联立(18)和(19)可得

Q 1 k 1 , Q 2 k 2 . (21)

0 Q 3 <1 ,由(20)可以推出

1 r 3 ( ) Q 3 + b 1 Q 1 + b 2 Q 2 0. (22)

联立(18)、(19)和(22)可得

Q 1 u D , Q 2 v D , Q 3 1 w D . (23)

注意到 Q i P i k i + ( i=1,2,3 ) ,结合(15)和(21)可知

P i = Q i = k i ,i=1,2,3.

又由(17)和(23)可知

P 1 = Q 1 = u D , P 2 = Q 2 = v D , P 3 = Q 3 =1 w D .

D<0 ,由(H2)可知 w D <0 ,那么 1 w D >1 。若 D>0 ,通过简单计算并结合(H2)可知 1 w D <0 ,即 1 w D [ 0,1 ] 。这与 0 Q 3 P 3 k 3 =1 矛盾。

综上可知, lim ξ ( ϕ 1 ( ξ ), ϕ 2 ( ξ ), ϕ 3 ( ξ ) )=( k 1 , k 2 , k 3 ) 。证毕。

4. 小结

本文研究了移动环境下三种群Lotka-Volterra竞争合作系统行波解的存在性,主要通过构造上下解和先验集,并结合单调迭代技巧和波动引理证明了该系统非负受迫行波的存在性。同时注意到受迫波的边界条件为 lim ξ ( ϕ 1 ( ξ ), ϕ 2 ( ξ ), ϕ 3 ( ξ ) )=( k 1 , k 2 , k 3 ) ,由于 r( ) 恒正,我们有 k 1 , k 2 , k 3 >0 ,这表明种群栖息地环境是轻微恶化的,由此可知这三个种群在任一固定的区域上都可以持久生存。

参考文献

[1] Fisher, R.A. (1937) The Wave of Advance of Advantageous Genes. Annals of Eugenics, 7, 355-369.
https://doi.org/10.1111/j.1469-1809.1937.tb02153.x
[2] Kolomgorov, A.N., Petrovskii, I.G. and Piskunov, N.S. (1937) Study of a Diffusion Equation That Is Related to the Growth of a Quality of Matter, and Its Application to a Biological Problem. Moscow University Mathematics Bulletin, 1, 1-26.
[3] Cantrell, R.S. and Cosner, C. (2004) Spatial Ecology via Reaction‐Diffusion Equations. John Wiley & Sons.
https://doi.org/10.1002/0470871296
[4] 倪维明. 浅谈反应扩散方程[J]. 数学传播, 2016, 34(4): 17-26.
[5] 楼元. 空间生态学中的一些反应扩散方程模型[J]. 中国科学: 数学, 2015, 45(10): 1619-1634.
[6] Yang, Y., Wu, C. and Li, Z. (2019) Forced Waves and Their Asymptotics in a Lotka-Volterra Cooperative Model under Climate Change. Applied Mathematics and Computation, 353, 254-264.
https://doi.org/10.1016/j.amc.2019.01.058
[7] Dong, F., Li, B. and Li, W. (2021) Forced Waves in a Lotka-Volterra Competition-Diffusion Model with a Shifting Habitat. Journal of Differential Equations, 276, 433-459.
https://doi.org/10.1016/j.jde.2020.12.022
[8] Sharmila, N.B. and Gunasundari, C. (2022) Travelling Wave Solutions for a Diffusive Prey-Predator Model with One Predator and Two Preys. International Journal of Apllied Mathematics, 35, 661-684.
https://doi.org/10.12732/ijam.v35i5.3
[9] Chen, C., Hung, L., Mimura, M. and Ueyama, D. (2012) Exact Travelling Wave Solutions of Three-Species Competition-Diffusion Systems. Discrete and Continuous Dynamical Systems-Series B, 17, 2653-2669.
https://doi.org/10.3934/dcdsb.2012.17.2653
[10] Mimura, M. and Tohma, M. (2015) Dynamic Coexistence in a Three-Species Competition-Diffusion System. Ecological Complexity, 21, 215-232.
https://doi.org/10.1016/j.ecocom.2014.05.004
[11] Hsu, C.-H., Lin, J.-J. and Wu, S.-L. (2019) Existence and Stability of Traveling Wavefronts for Discrete Three Species Competitive-Cooperative Systems. Mathematical Biosciences and Engineering, 16, 4151-4181.
https://doi.org/10.3934/mbe.2019207
[12] Wu, C., Yang, Y. and Wu, Z. (2021) Existence and Uniqueness of Forced Waves in a Delayed Reaction-Diffusion Equation in a Shifting Environment. Nonlinear Analysis: Real World Applications, 57, Article 103198.
https://doi.org/10.1016/j.nonrwa.2020.103198
[13] Yan, R., Liu, G. and Wang, Y. (2022) Stability of Bistable Traveling Wavefronts for a Three Species Competitive-Cooperative System with Nonlocal Dispersal. Japan Journal of Industrial and Applied Mathematics, 39, 515-541.
https://doi.org/10.1007/s13160-021-00497-5
[14] Yang, Z., Zhang, G. and He, J. (2023) Existence and Stability of Traveling Wavefronts for a Three‐Species Lotka-Volterra Competitive‐Cooperative System with Nonlocal Dispersal. Mathematical Methods in the Applied Sciences, 46, 13051-13073.
https://doi.org/10.1002/mma.9232
[15] Hu, H. and Zou, X. (2021) Traveling Waves of a Diffusive SIR Epidemic Model with General Nonlinear Incidence and Infinitely Distributed Latency but without Demography. Nonlinear Analysis: Real World Applications, 58, Article 103224.
https://doi.org/10.1016/j.nonrwa.2020.103224
[16] Smith, H.L. (2011) An Introduction to Delay Differential Equations with Applications to the Life Sciences. Springer.
https://doi.org/10.1007/978-1-4419-7646-8