aam Advances in Applied Mathematics 2324-7991 2324-8009 beplay体育官网网页版等您来挑战! 10.12677/aam.2024.139406 aam-96535 Articles 数学与物理 一类格上Lotka-Volterra合作系统受迫行波解的存在性
Existence of Forced Traveling Waves for a Class of the Lattice Lotka-Volterra Cooperative System
吴振宇 长沙理工大学数学与统计学院,湖南 长沙 09 09 2024 13 09 4253 4265 14 8 :2024 7 8 :2024 7 9 :2024 Copyright © 2024 beplay安卓登录 All rights reserved. 2024 This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/ 为了描述移动环境轻度恶化对两个弱合作物种的持久性产生的影响,本文考虑一类带有与时空均相关的恒正内禀增长函数的格上Lotka-Volterra合作系统。通过构造合适的上下解并结合单调迭代的方法证明了系统存在两组受迫行波解。
In order to characterize the effect of mild deterioration of the shifting environment on the two weakly cooperative species persistence, we consider a class of the lattice Lotka-Volterra cooperative systems with a constant positive intrinsic growth function that is spatio-temporally correlated. By constructing suitable upper and slower solutions combined with the method of monotone iteration, we prove that there exist two sets of forced traveling wave solutions for the system.
Lotka-Volterra合作系统,移动环境,受迫行波
Lotka-Volterra Cooperative System
Shifting Environment Forced Travelling Wave
1. 引言

近几十年来,由于全球变暖引起的气候持续恶化,使得地球生态圈受到了严重影响 [1] 。反应扩散方程作为研究空间生态学的一种重要工具,通过研究其行波解,可以用来刻画各种生物、物理现象,有着较强的实际意义 [2] [3] 。特别地,Berestycki通过建立如下模型探究了气候变化对单个物种的扩散影响 [4]

u ( t , x ) t = d 2 u ( t , x ) x 2 + f ( x c t , u ) , t 0 , x , (1.1)

其中 u ( t , x ) 表示物种u在t时刻位置x处的种群密度, d > 0 表示扩散系数, f ( x c t , u ) 表示物种的生长和种群密度u以及环境移动速度c有关。文 [4] 证明了当环境迁移速度很快时,物种将会灭绝。在此之后,许多学者通过考虑方程(1.1)及其推广形式探究了不同情形下单个物种的动力学行为 [5] - [14] 。在自然界生态系统中,物种间普遍存在着合作、竞争、捕食等关系,一个自然的问题是,受种间关系影响时,物种能否在栖息地中持久生存?最近,Yang和Wu等人通过建立如下Lotka-Volterra合作系统研究了当环境严重恶化时两个合作物种的动力学行为 [15]

{ u ( t , x ) t = d 1 2 u ( t , x ) x 2 + u ( t , x ) [ r 1 ( x c t ) u ( t , x ) + a 1 v ] , v ( t , x ) t = d 2 2 v ( t , x ) x 2 + v ( t , x ) [ r 2 ( x c t ) v ( t , x ) + a 2 u ] , t 0 , x , (1.2)

其中 u ( t , x ) , v ( t , x ) 分别表示物种 u , v 在t时刻位置x处的种群密度; d 1 > 0 , d 2 > 0 为扩散系数;种群的内禀增长速率 r i ( t , x ) : = r i ( x c t ) , i = 1 , 2 ,反映了栖息地以一个恒正的速度向右移动; a 1 , a 2 代表两物种的合作强度。并且有如下假设:

(A) 0 < a , a , < 1

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

其中条件(A)表示两物种为弱合作关系,条件(R)包含对 r i ( ) 的两个主要要求:单调性和符号变化。具体来说, r i ( ) 非减表示随着时间的推移,栖息地环境呈恶化趋势。而 r i ( ) < 0 < r i ( + ) , i = 1 , 2 表示恶化的程度十分严重。通过构造合适的上下解并结合单调迭代,文 [15] 证明了对于任意给定的栖息地边缘正速度,系统(1.2)都存在一个非减受迫波,并指出无论栖息地移动的有多慢,由于弱合作关系的存在,两个物种仍将灭绝,此受迫波代表一类波速与栖息地移动速度保持相同的行波。关于系统(1.2)的更多结果可以参考 [16] [17] [18]

注意到系统(1.2)是在连续环境下讨论的,而某些物种的生存环境是由无限多个斑块所组成,例如生活在分散树洞中的松鼠、池塘里的青蛙等。为了更好地描述此情形下物种的扩散过程,需要将系统离散化,即研究格上Lotka-Volterra系统。例如,Wang和Pan等人通过建立如下格上Lotka-Volterra竞争系统研究了两个竞争物种在移动环境下的传播问题 [19]

{ u ( t , x ) t = d 1 N 2 [ u ] ( t , x ) + u ( t , x ) [ r 1 ( x c t ) u ( t , x ) b 1 v ] , v ( t , x ) t = d 2 N 2 [ v ] ( t , x ) + v ( t , x ) [ r 2 ( x c t ) v ( t , x ) b 2 u ] , t 0 , x , (1.3)

其中 N 2 [ W ] ( t , x ) = W ( x + 1 , t ) 2 W ( x , t ) + W ( x 1 , t ) , W = u , v b 1 , b 2 表示两物种间的竞争强度,且 r i ( ) 满足如下假设:

(R1) r 1 ( ) 上连续非增的有界函数且 r 1 ( + ) < 0 < r 1 ( ) r 2 ( ) 上连续非减的有界函数且 r 2 ( ) < 0 < r 2 ( + )

这表明两个相互竞争的物种的周围环境的变化是不同的。在此条件下,Wang和Pan通过构造上下解结合单调迭代的方法探究了系统(1.3)的受迫波的存在性以及其间隙形成。随后,Qiao等人进一步地考虑了此受迫波的渐近行为 [20] 。对于系统(1.3)在条件(R)下的受迫波的研究可以参考 [21] - [24]

值得注意的是,自然界中的某些物种不仅能在优质区域内生存(相对快),也能在劣质区域内生存(相当慢),这导致 r i ( ) 还存在下面这种无符号变化的情形:

(R+) r i ( ) 上连续非减的有界函数,且 r i ( + ) r i ( ) > 0 , i = 1 , 2

这意味着虽然栖息地环境呈恶化趋势,但恶化的程度较低。在此情形下,两个弱合作物种是否能在栖息地内持久生存是一个值得研究的问题。受上述内容启发,本文考虑如下格上Lotka-Volterra合作系统:

{ u ( t , x ) t = d 1 N 2 [ u ] ( t , x ) + u ( t , x ) [ r 1 ( x c t ) u ( t , x ) + a 1 v ] , v ( t , x ) t = d 2 N 2 [ v ] ( t , x ) + v ( t , x ) [ r 2 ( x c t ) v ( t , x ) + a 2 u ] , t 0 , x , (1.4)

且假设条件(A)和条件(R+)成立。

记系统(1.4)的受迫波为 ( U ( ξ ) , V ( ξ ) ) , ξ = x c t ,代入到系统(1.4)中可得

{ d 1 N 2 [ U ] ( ξ ) + c U ( ξ ) + U ( ξ ) [ r 1 ( ξ ) U ( ξ ) + a 1 V ( ξ ) ] = 0 , d 2 N 2 [ V ] ( ξ ) + c V ( ξ ) + V ( ξ ) [ r 2 ( ξ ) V ( ξ ) + a 2 U ( ξ ) ] = 0 , ξ , (1.5)

其中 N 2 [ W ] ( t , x ) = W ( ξ + 1 ) 2 W ( ξ ) + W ( ξ 1 ) , W = U , V 。显然,系统(1.5)的解即为系统(1.4)的受迫波。根据(A)和(R+)可知,系统(1.5)的极限系统存在两个正平衡点 ( k 1 , k 2 ) ( l 1 , l 2 ) ,其中

k 1 = r 1 ( + ) + a 1 r 2 ( + ) 1 a 1 a 2 , k 2 = r 2 ( + ) + a 2 r 1 ( + ) 1 a 1 a 2 , l 1 = r 1 ( ) + a 1 r 2 ( ) 1 a 1 a 2 , l 2 = r 2 ( ) + a 2 r 1 ( ) 1 a 1 a 2 .

因此,自然地考虑到系统(1.4)是否存在连接这两个正平衡点的受迫波。具体来说,通过上下解并结合单调迭代的方法可以证明系统(1.4)存在一组非减受迫波,且该受迫波连接这两个正平衡点。进一步地,作如下假设:

(R*) 存在一个正常数 ζ ,使得当 ξ + 时, r i ( + ) r i ( ξ ) = o ( e ς ξ ) , i = 1 , 2

在此附加条件下,通过构造另一组合适的上下解并结合单调迭代可以证明系统(1.4)还存在一组连接 ( l 1 , l 2 ) ( 0 , 0 ) 的受迫波。

本文其余部分安排如下:在第二节中给出相关的预备知识,首先根据有序上下解的定义构造先验集并定义一个恰当的积分算子,然后证明该积分算子在先验集上的单调性以及先验集在此积分算子下的不变性。在第3节构造两组合适的有序上下解以及先验集,结合单调迭代和相关分析技巧得到系统(1.4)的两个受迫波的存在性。在第4节总结本文,结果表明种群栖息地环境轻微恶化时(即 r i ( ) > 0 , i = 1 , 2 ),对于任意给定的栖息地位置,经过足够长时间,该位置上所研究的两弱合作物种都不会灭绝。

2. 预备知识

C ( , ) 上全体连续函数所构成的空间, B C ( , ) 表示 上全体连续有界函数所构成的空间,且被陚以极值范数:

U = sup ξ R | U ( ξ ) | ,

同时,定义 C ( , ) 上的正锥 C + = { U C ( , ) | U ( ξ ) 0 , ξ } 。对于任意的 U 1 , U 2 C ( , ) ,若 U 1 U 2 C + ,我们记为 U 1 U 2 或者 U 2 U 1 。进一步地,对于任意的 U = ( U 1 , U 2 ) , V = ( V 1 , V 2 ) C × C ,若 U 1 V 1 U 2 V 2 ,则记 U V 或者 V U

下面给出系统(1.5)的有序上下解的定义。

定义1 设 ( U ( ξ ) , V ( ξ ) ) B C ( , ) { ξ j } j = 1 N 为一列有限递增点列。若对于任意的 ξ \ { ξ j } ( U ( ξ ) , V ( ξ ) ) 分别满足如下不等式组:

{ c U ( ξ ) ( ) d 1 [ U ( ξ + 1 ) 2 U ( ξ ) + U ( ξ 1 ) ] + U ( ξ ) [ r 1 ( ξ ) U ( ξ ) + a 1 V ( ξ ) ] , c V ( ξ ) ( ) d 2 [ V ( ξ + 1 ) 2 V ( ξ ) + V ( ξ 1 ) ] + V ( ξ ) [ r 2 ( ξ ) V ( ξ ) + a 2 U ( ξ ) ] ,

则称 ( U ( ξ ) , V ( ξ ) ) 为系统(1.5)的上(下)解。进一步地,若上解 ( U ¯ ( ξ ) , V ¯ ( ξ ) ) 和下解 ( U _ ( ξ ) , V _ ( ξ ) ) 满足 ( U ¯ ( ξ ) , V ¯ ( ξ ) ) ( U _ ( ξ ) , V _ ( ξ ) ) ,则称 ( U ¯ ( ξ ) , V ¯ ( ξ ) ) ( U _ ( ξ ) , V _ ( ξ ) ) 为系统(1.5)的一对有序上下解。

根据此定义构造如下非空集合:

Ω : = { ( U , V ) B C ( , ) × B C ( , ) | ( U _ , V _ ) ( U , V ) ( U ¯ , V ¯ ) } .(2.1)

再取参数

μ 1 = max { 2 d 1 + 2 max ξ U ¯ ( ξ ) , c ρ + 1 } μ 2 = max { 2 d 2 + 2 max ξ V ¯ ( ξ ) , c ρ + 1 } ,

以及定义算子

H 1 ( U , V ) ( ξ ) = μ 1 U ( ξ ) + d 1 [ U ( ξ + 1 ) 2 U ( ξ ) + U ( ξ 1 ) ] + U ( ξ ) [ r 1 ( ξ ) U ( ξ ) + a 1 V ( ξ ) ] ,

H 2 ( U , V ) ( ξ ) = μ 2 V ( ξ ) + d 2 [ V ( ξ + 1 ) 2 V ( ξ ) + V ( ξ 1 ) ] + V ( ξ ) [ r 2 ( ξ ) V ( ξ ) + a 2 U ( ξ ) ] .

容易验证如下方程组的解即为系统(1.5)的解。

{ U ( ξ ) = 1 c ξ e μ 1 c ( s ξ ) H 1 ( U , V ) ( s ) d s , V ( ξ ) = 1 c ξ e μ 2 c ( s ξ ) H 2 ( U , V ) ( s ) d s . (2.2)

进一步地,定义算子 F ( U , V ) ( ξ ) = ( F 1 ( U , V ) ( ξ ) , F 2 ( U , V ) ( ξ ) ) ,其中

F i ( U , V ) ( ξ ) = 1 c ξ e μ i c ( s ξ ) H i ( U , V ) ( s ) d s , i = 1 , 2. (2.3)

显然,算子F在Ω中的不动点即为系统(1.5)的解。

为了运用单调迭代,我们首先通过如下引理说明算子F的一些性质。

引理1 下列结论成立:

(i) 当 ( U , V ) Ω 时,F是一个非减算子;

(ii) 若 ( U ( ξ ) , V ( ξ ) ) Ω 关于 ξ 非减,则 F ( U , V ) ( ξ ) 关于 ξ 也非减;

(iii) F ( Ω ) Ω

证明. (i)令 ( U ˜ , V ˜ ) , ( U ^ , V ^ ) Ω 并且满足 ( U ˜ , V ˜ ) ( U ^ , V ^ ) ,则对于任意的 ξ

H 1 ( U ˜ , V ˜ ) ( ξ ) H 1 ( U ^ , V ^ ) ( ξ ) = μ 1 [ U ˜ ( ξ ) U ^ ( ξ ) ] + d 1 { [ U ˜ ( ξ + 1 ) 2 U ˜ ( ξ ) + U ˜ ( ξ 1 ) ] [ U ^ ( ξ + 1 ) 2 U ^ ( ξ ) + U ^ ( ξ 1 ) ] } + U ˜ ( ξ ) [ r 1 ( ξ ) U ˜ ( ξ ) + a 1 V ˜ ( ξ ) ] U ^ ( ξ ) [ r 1 ( ξ ) U ^ ( ξ ) + a 1 V ^ ( ξ ) ] = [ μ 1 2 d 1 U ˜ ( ξ ) U ^ ( ξ ) + r 1 ( ξ ) ] [ U ˜ ( ξ ) U ^ ( ξ ) ] + d 1 { [ U ˜ ( ξ + 1 ) U ^ ( ξ + 1 ) ] + [ U ˜ ( ξ 1 ) U ^ ( ξ 1 ) ] } + a 1 [ U ˜ ( ξ ) V ˜ ( ξ ) U ^ ( ξ ) V ^ ( ξ ) ] 0

以及 H 2 ( U ˜ , V ˜ ) ( ξ ) H 2 ( U ^ , V ^ ) ( ξ ) 0 ,再根据 F 1 F 2 的定义可知

F i ( U ˜ , V ˜ ) ( ξ ) F i ( U ^ , V ^ ) ( ξ ) = 1 c ξ e μ i c ( s ξ ) [ H i ( U ˜ , V ˜ ) ( s ) H i ( U ^ , V ^ ) ( s ) ] d s 0 , i = 1 , 2.

因此F是一个非减算子。

(ii) 假设 ( U ( ξ ) , V ( ξ ) ) Ω 关于 ξ 非减,那么对于任意的 ε > 0 以及任意的 ξ ,有

H 1 ( U , V ) ( ξ + ε ) H 1 ( U , V ) ( ξ ) = μ 1 [ U ( ξ + ε ) U ( ξ ) ] + U ( ξ + ε ) [ r 1 ( ξ ) U ( ξ + ε ) + a 1 V ( ξ + ε ) ] U ( ξ ) [ r 1 ( ξ ) U ( ξ ) + a 1 V ( ξ ) ] + d 1 { [ U ( ξ + ε + 1 ) 2 U ( ξ + ε ) + U ( ξ + ε 1 ) ] [ U ( ξ + 1 ) 2 U ( ξ ) + U ( ξ 1 ) ] } = [ μ 1 2 d 1 U ( ξ + ε ) U ( ξ ) ] [ U ( ξ + ε ) U ( ξ ) ] + d 1 { [ U ( ξ + ε + 1 ) U ( ξ + 1 ) ] + [ U ( ξ + ε 1 ) U ( ξ 1 ) ] } + U ( ξ + ε ) r 1 ( ξ + ε ) U ( ξ ) r 1 ( ξ ) + a 1 [ U ( ξ + ε ) V ( ξ + ε ) U ( ξ ) V ( ξ ) ] 0

以及 H 2 ( U , V ) ( ξ + ε ) H 2 ( U , V ) ( ξ ) 0 ,进而得到

F i ( U , V ) ( ξ + ε ) = 1 c ξ + ε e μ i c ( s ξ ε ) H i ( U , V ) ( s ) d s = 1 c ξ e μ i c ( s ξ ) H i ( U , V ) ( s + ε ) d s 1 c ξ e μ i c ( s ξ ) H i ( U , V ) ( s ) d s = F i ( U , V ) ( ξ ) , i = 1 , 2 ,

F ( U , V ) ( ξ ) 关于 ξ 非减。

(iii) 记 ( U ¯ ( ξ ) , V ¯ ( ξ ) ) ( U _ ( ξ ) , V _ ( ξ ) ) 为系统(1.5)的一对有序上下解,那么由上解及F的定义可知

F 1 ( U ¯ , V ¯ ) ( ξ ) = 1 c ξ e μ 1 c ( s ξ ) H 1 ( U ¯ , V ¯ ) ( s ) d s = 1 c ( ξ ξ 1 + ξ 1 ξ 2 + + ξ N ) e μ 1 c ( s ξ ) H 1 ( U ¯ , V ¯ ) ( s ) d s 1 c ( ξ ξ 1 + ξ 1 ξ 2 + + ξ N ) [ c U ¯ ( s ) + μ 1 U ¯ ( s ) ] e μ 1 c ( s ξ ) d s = U ¯ ( ξ )

以及 F 2 ( U ¯ , V ¯ ) ( ξ ) = 1 c ξ e μ 2 c ( s ξ ) H 2 ( U ¯ , V ¯ ) ( s ) d s V ¯ ( ξ ) ,这表示 F ( U ¯ , V ¯ ) ( U ¯ , V ¯ ) 。类似地,由下解以及F的定义可知

F 1 ( U _ , V _ ) ( ξ ) = 1 c ξ e μ 1 c ( s ξ ) H 1 ( U _ , V _ ) ( s ) d s = 1 c ( ξ ξ 1 + ξ 1 ξ 2 + + ξ N ) e μ 1 c ( s ξ ) H 1 ( U _ , V _ ) ( s ) d s 1 c ( ξ ξ 1 + ξ 1 ξ 2 + + ξ N ) [ c U _ ( s ) + μ 1 U _ ( s ) ] e μ 1 c ( s ξ ) d s = U _ ( ξ )

以及 F 2 ( U _ , V _ ) ( ξ ) V _ ( ξ ) ,这表示 ( U _ , V _ ) F ( U _ , V _ ) 。注意到 F ( U , V ) 是一个非减函数,因此对于任意的 ( U , V ) Ω 都有

( U _ , V _ ) F ( U _ , V _ ) F ( U , V ) F ( U ¯ , V ¯ ) ( U ¯ , V ¯ ) , (2.4)

从而可知 F ( Ω ) Ω 。证毕。

3. 主要结果 3.1. 上下解的构造

由于上述引理的结论是在先验集Ω的前提下得到的,且系统(1.4)的受迫波的边界条件也依赖于此先验集。所以下面构造系统(1.5)的两对有序上下解以便定义两个合适的先验集。

经过直接计算可得如下引理:

引理2 令 ( U ¯ 1 ( ξ ) , V ¯ 1 ( ξ ) ) = ( k 1 , k 2 ) ( U _ 1 ( ξ ) , V _ 1 ( ξ ) ) = ( l 1 , l 2 ) ,则 ( U ¯ 1 ( ξ ) , V ¯ 1 ( ξ ) ) ( U _ 1 ( ξ ) , V _ 1 ( ξ ) ) 是系统(1.5)的一对有序上下解。

接下来构造另外一对有序上下解。

引理3 定义函数 θ i ( c , λ ) = d i ( e λ + e λ 2 ) + r i ( + ) c λ , i = 1 , 2 。令

c i * ( ) : = min λ > 0 d i ( e λ + e λ 2 ) + r i ( + ) λ , i = 1 , 2 ;

c ^ * ( ) : = max { c 1 * ( ) , c 2 * ( ) } ;

c ¯ * ( ) : = min { c 1 * ( ) , c 2 * ( ) } .

根据函数的凸性,容易推出如下情形成立:

(i) 当 c = c i * ( ) , i = 1 , 2 时, θ i ( c , λ ) = 0 存在正根 λ i * ( ) , i = 1 , 2

(ii) 当 c > c ^ * ( ) 时, θ 1 ( c , λ ) θ 2 ( c , λ ) 分别有两个不同正根 λ 1 , λ 3 λ 2 , λ 4 ,且满足

θ 1 ( c , λ ) = { < 0 , λ ( λ 1 , λ 3 ) > 0 , λ [ 0 , λ 1 ) ( λ 3 , + ) , θ 2 ( c , λ ) = { < 0 , λ ( λ 2 , λ 4 ) > 0 , λ [ 0 , λ 2 ) ( λ 4 , + ) ;

(iii) 当 c < c ¯ * ( ) 时, θ i ( c , λ ) = 0 , i = 1 , 2 无实数根。

假设 Λ : = ( λ 1 , λ 3 ) ( λ 2 , λ 4 ) 非空,则存在 γ Λ 并且满足 λ 1 + λ 2 > γ 。那么,对任意的 c > c i * ( ) ,有 θ i ( c , λ ) < 0 , i = 1 , 2 。需要注意的是,集合Λ非空在许多情形下都会成立,我们在此说明当c比较大的时候假设成立。不妨设 λ 2 > λ 1 ,根据 θ 1 ( c , λ 1 ) = 0 可知

θ 1 ( c , λ 2 ) = d 1 ( e λ 2 + e λ 2 2 ) + r 1 ( + ) c λ 2 = d 1 ( e λ 2 + e λ 2 2 ) d 1 ( e λ 1 + e λ 1 2 ) c ( λ 2 λ 1 )

这表示当c充分大时, θ 1 ( c , λ 2 ) < 0 ,根据引理3可知集合Λ非空。此外,由于 λ 1 + λ 2 > γ 以及 θ i ( c , λ ) < 0 , i = 1 , 2 ,所以存在充分大的正数 ξ 1 , ξ 2 分别满足

k 1 θ 1 ( c , γ ) + a 1 k 1 e λ 2 ξ 1 + a 1 e ( λ 2 + λ 1 γ ) ξ 1 = 0 , k 2 θ 2 ( c , γ ) + a 2 k 2 e λ 1 ξ 2 + a 2 e ( λ 2 + λ 1 γ ) ξ 2 = 0.

由此定义两个连续函数:

U ¯ 2 ( ξ ) : = min { k 1 , e λ 1 ξ + q k 1 e γ ξ } , V ¯ 2 ( ξ ) : = min { k 2 , e λ 2 ξ + q k 2 e γ ξ } .

其中 q > 1 充分大使得 k i e λ 1 ξ i + q k i e γ ξ i , i = 1 , 2

根据引理3和条件(R*)可知,存在充分大的常数L以及常数 η ( 0 , λ i * ( ) λ i ) , i = 1 , 2 ,使得当 ξ L 时,

r i ( + ) r i ( ξ ) e ζ ξ e η ξ , i = 1 , 2 ,

并且当 c > c ^ * ( ) 时,

θ i ( c , λ ) < 0 , i = 1 , 2 .

M i = 2 θ i ( c , λ i + η ) , i = 1 , 2 ,定义两个连续有界函数:

U _ 2 ( ξ ) = { κ 1 , ξ ξ 3 e λ 1 ξ ( 1 M 1 e η ξ ) , ξ > ξ 3 , V _ 2 ( ξ ) = { κ 2 , ξ ξ 4 e λ 2 ξ ( 1 M 1 e η ξ ) , ξ > ξ 4 .

其中 ξ i + 2 = ln r i ( + ) λ i κ i : = e λ i ξ i + 2 ( 1 M i e η ξ i + 2 ) ξ i + 2 充分大使得 0 < κ i < r i ( ) , i = 1 , 2

引理4 假设(R+)、(R*)成立且 c > c ^ * ( ) 使得 Λ : = ( λ 1 , λ 3 ) ( λ 2 , λ 4 ) 非空,则 ( U ¯ 2 ( ξ ) , V ¯ 2 ( ξ ) ) ( U _ 2 ( ξ ) , V _ 2 ( ξ ) ) 是系统(1.5)的一对有序上下解。

证明. 首先说明 U ¯ 2 ( ξ ) R \ ξ 1 上满足

c U ¯ 2 ( ξ ) d 1 [ U ¯ 2 ( ξ + 1 ) 2 U ¯ 2 ( ξ ) + U ¯ 2 ( ξ 1 ) ] + U ¯ 2 ( ξ ) [ r 1 ( ξ ) U ¯ 2 ( ξ ) + a 1 V ¯ 2 ( ξ ) ] ,

U ¯ 2 ( ξ ) = k 1 时,不等式显然成立。若 U ¯ 2 ( ξ ) = e λ 1 ξ + q k 1 e γ ξ ,意味着 ξ ξ 1 。根据 θ 1 ( c , λ 1 ) = 0 以及 θ 1 ( c , γ ) < 0 可知

c ( e λ 1 ξ + q k 1 e γ ξ ) + d 1 [ e λ 1 ( ξ + 1 ) + q k 1 e γ ( ξ + 1 ) 2 ( e λ 1 ξ + q k 1 e γ ξ ) + e λ 1 ( ξ 1 ) + q k 1 e γ ( ξ 1 ) ] + ( e λ 1 ξ + q k 1 e γ ξ ) [ r 1 ( ξ ) ( e λ 1 ξ + q k 1 e γ ξ ) + a 1 ( e λ 1 ξ + q k 1 e γ ξ ) ] c ( e λ 1 ξ + q k 1 e γ ξ ) + d 1 [ e λ 1 ( ξ + 1 ) + q k 1 e γ ( ξ + 1 ) 2 ( e λ 1 ξ + q k 1 e γ ξ ) + e λ 1 ( ξ 1 ) + q k 1 e γ ( ξ 1 ) ] + ( e λ 1 ξ + q k 1 e γ ξ ) [ r 1 ( + ) ( e λ 1 ξ + q k 1 e γ ξ ) + a 1 ( e λ 1 ξ + q k 1 e γ ξ ) ] q k 1 e γ ξ θ 1 ( c , γ ) + a 1 e ( λ 2 + λ 1 ) ξ + a 1 q k 1 e ( λ 2 + γ ) ξ e γ ξ [ q k 1 θ 1 ( c , γ ) + q a 1 e ( λ 2 + λ 1 γ ) ξ + a 1 q k 1 e λ 2 ξ ] = q e γ ξ [ k 1 θ 1 ( c , γ ) + a 1 e ( λ 2 + λ 1 γ ) ξ + a 1 k 1 e λ 2 ξ ] = 0 ,

因此 U ¯ 2 ( ξ ) \ ξ 1 上满足上解的定义。类似地可推出

c V ¯ 2 ( ξ ) d 1 [ V ¯ 2 ( ξ + 1 ) 2 V ¯ 2 ( ξ ) + V ¯ 2 ( ξ 1 ) ] + V ¯ 2 ( ξ ) [ r 2 ( ξ ) V ¯ 2 ( ξ ) + a 2 U ¯ 2 ( ξ ) ]

\ ξ 2 上成立。从而 ( U ¯ 2 ( ξ ) , V ¯ 2 ( ξ ) ) 是系统(1.5)的一组上解。

接下来验证 U _ 2 ξ ξ 3 时满足下解的定义。当 ξ > ξ 3 时,由条件(R*)结合 θ 1 ( c , λ 1 ) = 0 可知

d 1 [ e λ 1 ( ξ + 1 ) ( 1 M 1 e η ( ξ + 1 ) ) 2 e λ 1 ξ ( 1 M 1 e η ξ ) + e λ 1 ( ξ 1 ) ( 1 M 1 e η ( ξ 1 ) ) ] + e λ 1 ξ ( 1 M 1 e η ξ ) [ r 1 ( ξ ) e λ 1 ξ ( 1 M 1 e η ξ ) + a 1 e λ 1 ξ ( 1 M 2 e η ξ ) ] c λ 1 e λ 1 ξ + c M 1 ( λ 1 + η ) e ( λ 1 + η ) ξ M 1 e ( λ 1 + η ) ξ [ θ 1 ( c , λ 1 + η ) + r 1 ( ξ ) r 1 ( + ) e λ 1 ξ ( 1 M 1 e η ξ ) ] + e λ 1 ξ [ r 1 ( ξ ) r 1 ( + ) e λ 1 ξ ( 1 M 1 e η ξ ) ] e λ 1 ξ [ M 1 e η ξ θ 1 ( c , λ 1 + η ) e η ξ e λ 1 ξ ] e λ 1 ξ [ M 1 e η ξ θ 1 ( c , λ 1 + η ) 2 e η ξ ] = e ( λ 1 + η ) ξ [ M 1 θ 1 ( c , λ 1 + η ) 2 ] = 0.

ξ < ξ 3 时,容易验证 U _ 2 ( ξ ) = κ 1 满足下解的定义。类似地可以证明

c V _ 2 ( ξ ) d 2 [ V _ 2 ( ξ + 1 ) 2 V _ 2 ( ξ ) + V _ 2 ( ξ 1 ) ] + V _ 2 ( ξ ) [ r 2 ( ξ ) V _ 2 ( ξ ) + a 2 U _ 2 ( ξ ) ]

\ ξ 4 上成立。因此 ( U _ 2 ( ξ ) , V _ 2 ( ξ ) ) 是系统(1.5)的一组下解。

另外通过下解的构造可知 U _ 2 min { r 1 ( + ) , e λ 1 ξ } V _ 2 min { r 2 ( + ) , e λ 2 ξ } ,从而 ( U ¯ 2 ( ξ ) , V ¯ 2 ( ξ ) ) ( U _ 2 ( ξ ) , V _ 2 ( ξ ) ) ,所以 ( U ¯ 2 ( ξ ) , V ¯ 2 ( ξ ) ) ( U _ 2 ( ξ ) , V _ 2 ( ξ ) ) 为系统(1.5)的一对有序上下解。证毕。

最后,根据所构造的上下解定义如下非空集合:

Ω 1 = { ( U , V ) B C ( , ) × B C ( , ) | ( U _ 1 , V _ 1 ) ( U , V ) ( U ¯ 1 , V ¯ 1 ) } , Ω 2 = { ( U , V ) B C ( , ) × B C ( , ) | ( U _ 2 , V _ 2 ) ( U , V ) ( U ¯ 2 , V ¯ 2 ) } .

需要指出的是,将引理1中的Ω换成 Ω i , i = 1 , 2 ,引理的结论仍然成立

3.2. 受迫波的存在性

在本节,我们证明系统(1.4)存在两个受迫波。

定理1 假设(A)和(R+)成立,则对于任意的 c > 0 ,系统(1.4)存在一组非减受迫波 ( U 1 ( ξ ) , V 1 ( ξ ) ) ,且满足如下边界条件:

lim ξ ( U 1 ( ξ ) , V 1 ( ξ ) ) = ( l 1 , l 2 ) , lim ξ + ( U 1 ( ξ ) , V 1 ( ξ ) ) = ( k 1 , k 2 ) .

证明. 首先根据算子 F ( Ω 1 ) Ω 1 构造如下迭代:

( U 1 1 , V 1 1 ) = F ( U ¯ 1 , V ¯ 1 ) , ( U 1 n + 1 , V 1 n + 1 ) = F ( U 1 n , V 1 n ) , n 1 ,

根据引理1,F是将 Ω 1 映射到 Ω 1 的非减算子,那么有

( U _ 1 , V _ 1 ) ( U 1 n + 1 , V 1 n + 1 ) ( U 1 n , V 1 n ) ( U ¯ 1 , V ¯ 1 ) , n 1 ,

又因为 U ¯ 1 V ¯ 1 是非减的,结合引理1可知,对每个给定的n, U 1 n ( ξ ) V 1 n ( ξ ) 也是非减的函数。令

lim n ( U 1 n , V 1 n ) = ( U 1 , V 1 ) ,

显然 ( U 1 , V 1 ) 是一个非减函数且

( U _ 1 , V _ 1 ) ( U 1 , V 1 ) ( U ¯ 1 , V ¯ 1 ) .

下面我们证明 ( U 1 , V 1 ) 是算子F在 Ω 1 中的不动点。由于 ( U 1 n , V 1 n ) 逐点收敛于 ( U 1 , V 1 ) ,根据 H 1 H 2 的定义可知, H 1 ( U 1 n , V 1 n ) H 2 ( U 1 n , V 1 n ) 分别逐点收敛到 H 1 ( U 1 , V 1 ) H 2 ( U 1 , V 1 ) 。此外,对任意的 n 1 都有

| H 1 ( U 1 n , V 1 n ) ( ξ ) | k 1 ( μ + 4 d 1 + r 1 ( + ) + k 1 + a 1 k 2 )

以及

| H 2 ( U 1 n , V 1 n ) ( ξ ) | k 2 ( μ 2 + 4 d 2 + r 2 ( + ) + k 2 + a 2 k 1 ) ,

根据勒贝格控制收敛定理可知

U 1 ( ξ ) = lim n + U 1 n ( ξ ) = lim n + F 1 ( U 1 n ( ξ ) , V 1 n ( ξ ) ) = lim n + 1 c ξ e μ 1 c ( s ξ ) H 1 ( U 1 n 1 , V 1 n 1 ) ( s ) d s = 1 c ξ e μ 1 c ( s ξ ) H ( U 1 , V 1 ) ( s ) d s = F 1 ( U 1 , V 1 ) ( ξ ) ,

类似地, V 1 ( ξ ) = F 2 ( U 1 , V 1 ) ( ξ ) 。因此 ( U 1 ( ξ ) , V 1 ( ξ ) ) 是算子F在 Ω 1 中的不动点,即为系统(1.5)的解。

下面证明 ( U 1 ( ξ ) , V 1 ( ξ ) ) 满足相应的边界条件。因为 U 1 ( ξ ) V 1 ( ξ ) 上的单调有界函数,所以存在常数 A 1 , A 2 , B 1 , B 2 满足

lim ξ ( U 1 ( ξ ) , V 1 ( ξ ) ) = ( A 1 , A 2 ) , lim ξ + ( U 1 ( ξ ) , V 1 ( ξ ) ) = ( B 1 , B 2 ) ,

运用L’Hôpital法则进行简单计算可知

A 1 = lim ξ U 1 ( ξ ) = lim ξ 1 c ξ e μ 1 c ( s ξ ) H 1 ( U 1 , V 1 ) ( s ) d s = 1 c lim ξ c H 1 ( U 1 , V 1 ) ( ξ ) μ 1 = A 1 + A 1 [ r 1 ( ) A 1 + a 1 A 2 ] μ 1 ,

以及

A 2 = A 2 + A 2 [ r 2 ( ) A 2 + a 2 A 1 ] μ 2 ,

( U 1 ( ξ ) , V 1 ( ξ ) ) Ω 1 ,所以 0 < l 1 A 1 k 1 0 < l 2 A 2 k 2 ,因此 { A 1 = l 1 A 2 = l 2

类似地可得

B 1 = lim ξ + V 1 ( ξ ) = lim ξ + 1 c ξ e μ 1 c ( s ξ ) H 1 ( U 1 , V 1 ) ( s ) d s = 1 c lim ξ + c H 1 ( U 1 , V 1 ) ( ξ ) μ 1 = B 1 + B 1 [ r 1 ( + ) B 1 + a 1 B 2 ] μ 1 ,

以及 B 2 = B 2 + B 2 [ r 2 ( + ) B 2 + a 2 B 1 ] μ 2 。又 ( U 1 ( ξ ) , V 1 ( ξ ) ) Ω 1 ,所以 0 < B 1 k 1 0 < B 2 k 2 。故 { B 1 = k 1 B 2 = k 2 。证毕。

定理2 假设(A)、(R+)、(R*)成立,并且 c > c ^ * ( ) 使得 Λ : = ( λ 1 , λ 3 ) ( λ 2 , λ 4 ) 非空,则系统(1.4)存在一个非负受迫波,且满足如下边界条件,

lim ξ ( U 2 ( ξ ) , V 2 ( ξ ) ) = ( l 1 , l 2 ) , lim ξ + ( U 2 ( ξ ) , V 2 ( ξ ) ) = ( 0 , 0 ) .

证明. 类似于定理1的证明过程并结合引理4,容易得到存在 ( U 2 ( ξ ) , V 2 ( ξ ) ) Ω 2 使得

F ( U 2 , V 2 ) ( ξ ) = ( U 2 ( ξ ) , V 2 ( ξ ) ) ,

( U 2 ( ξ ) , V 2 ( ξ ) ) 为系统(1.4)的解。此外,由 ( U _ 2 ( ξ ) , V _ 2 ( ξ ) ) ( 0 , 0 ) 可知 ( U 2 ( ξ ) , V 2 ( ξ ) ) 是非负的。

接着我们证明 ( U 2 ( ξ ) , V 2 ( ξ ) ) 的相应边界条件。首先根据 ( U ¯ 2 ( ξ ) , V ¯ 2 ( ξ ) ) ( U _ 2 ( ξ ) , V _ 2 ( ξ ) ) 的定义可知:

lim ξ ( U _ 2 ( ξ ) , V _ 2 ( ξ ) ) = ( 0 , 0 ) , lim ξ + ( U ¯ 2 ( ξ ) , V ¯ 2 ( ξ ) ) = ( 0 , 0 ) ,

( U 2 ( ξ ) , V 2 ( ξ ) ) Ω 2 ,因此

lim ξ ( U 2 ( ξ ) , V 2 ( ξ ) ) = ( 0 , 0 ) .

最后我们验证 lim ξ + ( U 2 ( ξ ) , V 2 ( ξ ) ) = ( l 1 , l 2 ) 。令

C 1 = lim ξ sup U 2 ( ξ ) , D 1 = lim ξ sup V 2 ( ξ ) , C 2 = lim ξ inf U 2 ( ξ ) , D 2 = lim ξ inf V 2 ( ξ ) .

根据 ( U 2 ( ξ ) , V 2 ( ξ ) ) Ω 2 可知, min { r 1 ( ) , r 2 ( ) } < C i , D i max { k 1 , k 2 } , i = 1 , 2 。从而根据波动引理 [25] ,存在一列满足 lim n + σ n = 的序列 { σ n } n = 1 以及满足 lim n + τ n = 的序列 { τ n } n = 1 ,使得

lim n + U 2 ( σ n ) = C 1 , lim n + V 2 ( σ n ) = D 1 , lim n + U 2 ( τ n ) = C 2 , lim n + V 2 ( τ n ) = D 2 ,

以及

lim n + U 2 ( σ n ) = lim n + U 2 ( τ n ) = lim n + V 2 ( σ n ) = lim n + V 2 ( τ n ) = 0.

根据系统(1.5)可得

{ c U 2 ( σ n ) = d 1 [ U 2 ( σ n + 1 ) 2 U 2 ( σ n ) + U 2 ( σ n 1 ) ] + U 2 ( σ n ) [ r 1 ( ξ ) U 2 ( σ n ) + a 1 V 2 ( σ n ) ] , c V 2 ( σ n ) = d 2 [ V 2 ( σ n + 1 ) 2 V 2 ( σ n ) + V 2 ( σ n 1 ) ] + V 2 ( σ n ) [ r 2 ( ξ ) V 2 ( σ n ) + a 2 U 2 ( σ n ) ] ,

{ c U 2 ( τ n ) = d 1 [ U 2 ( τ n + 1 ) 2 U 2 ( τ n ) + U 2 ( τ n 1 ) ] + U 2 ( τ n ) [ r 1 ( ξ ) U 2 ( τ n ) + a 1 V 2 ( τ n ) ] , c V 2 ( τ n ) = d 2 [ V 2 ( τ n + 1 ) 2 V 2 ( τ n ) + V 2 ( τ n 1 ) ] + V 2 ( τ n ) [ r 2 ( ξ ) V 2 ( τ n ) + a 2 U 2 ( τ n ) ] ,

令上述的每个等式中的 n + ,并且对第一个等式组中的两个等式的两边同时取上极限以及对第二个等式组中的两个等式的两边同时取下极限可得

{ 0 C 1 [ r 1 ( ) C 1 + a 1 D 1 ] 0 D 1 [ r 2 ( ) D 1 + a 2 C 1 ] { 0 C 2 [ r 1 ( ) C 2 + a 1 D 2 ] 0 D 2 [ r 2 ( ) D 2 + a 2 C 2 ] .

又因为 C i , D i > min { r 1 ( ) , r 2 ( ) } > 0 , i = 1 , 2 ,且 C 1 C 2 , D 1 D 2

所以

{ r 1 ( ) C 2 + a 1 D 2 r 1 ( ) C 1 + a 1 D 1 , r 2 ( ) D 2 + a 2 C 2 r 2 ( ) D 1 + a 2 C 1 ,

{ C 1 C 2 a 1 ( D 1 D 2 ) , D 1 D 2 a 2 ( C 1 C 2 ) ,

不难发现,该不等式组仅存在以下两种情况:

(i) C 1 = C 2 D 1 = D 2 。此时,不等式组显然成立,并且有

{ r 1 ( ) C 1 + a 1 D 1 = 0 r 2 ( ) D 1 + a 2 C 1 = 0 { r 1 ( ) C 2 + a 1 D 2 = 0 r 2 ( ) D 2 + a 2 C 2 = 0 .

计算可知

{ C 1 = C 2 = l 1 , D 1 = D 2 = l 2 .

(ii) C 1 C 2 D 1 D 2 ,即 C 1 > C 2 , D 1 > D 2 。通过计算可得

{ C 1 C 2 D 1 D 2 a 1 , D 1 D 2 C 1 C 2 a 2 ,

从而可知 a 1 a 2 1 ,与条件(A)矛盾。

因此 lim ξ + ( U 2 ( ξ ) , V 2 ( ξ ) ) = ( l 1 , l 2 ) 。证毕。

4. 总结

本文探究了一类格上Lotka-Volterra合作系统受迫行波解的存在性。通过构造合适的上下解并结合单调迭代的方法证明了系统(1.4)存在两组非负的受迫波,并且根据 r i ( ) > 0 , i = 1 , 2 以及这两组受迫波在 处都趋于 ( l 1 , l 2 ) 可知,当栖息地环境轻度恶化时,对于任意给定的栖息地位置,即使经过足够长时间,在该位置上所考虑的两弱合作物种都不会灭绝。此外,本文仅研究了系统(1.4)受迫波的存在性,还可以进一步考虑其唯一性、稳定性、渐近行为等,以及如何去掉定理2中的技术性条件(R*),将是后续研究工作的主题。

References Parmesan, C. (2006) Ecological and Evolutionary Responses to Recent Climate Change. Annual Review of Ecology, Evolution, and Systematics, 37, 637-669. >https://doi.org/10.1146/annurev.ecolsys.37.091305.110100 楼元. 空间生态学中的一些反应扩散方程模型[J]. 中国科学: 数学, 2015(10): 1619-1634. Fang, J., Lou, Y. and Wu, J. (2016) Can Pathogen Spread Keep Pace with Its Host Invasion? SIAM Journal on Applied Mathematics, 76, 1633-1657. >https://doi.org/10.1137/15m1029564 Li, B., Bewick, S., Shang, J. and Fagan, W.F. (2014) Persistence and Spread of a Species with a Shifting Habitat Edge. SIAM Journal on Applied Mathematics, 74, 1397-1417. >https://doi.org/10.1137/130938463 Hu, C. and Li, B. (2015) Spatial Dynamics for Lattice Differential Equations with a Shifting Habitat. Journal of Differential Equations, 259, 1967-1989. >https://doi.org/10.1016/j.jde.2015.03.025 Hu, H. and Zou, X. (2017) Existence of an Extinction Wave in the Fisher Equation with a Shifting Habitat. Proceedings of the American Mathematical Society, 145, 4763-4771. >https://doi.org/10.1090/proc/13687 Berestycki, H. and Fang, J. (2018) Forced Waves of the Fisher-KPP Equation in a Shifting Environment. Journal of Differential Equations, 264, 2157-2183. >https://doi.org/10.1016/j.jde.2017.10.016 Hu, H., Yi, T. and Zou, X. (2019) On Spatial-Temporal Dynamics of a Fisher-KPP Equation with a Shifting Environment. Proceedings of the American Mathematical Society, 148, 213-221. >https://doi.org/10.1090/proc/14659 Qiao, S., Li, W. and Wang, J. (2022) Multi-type Forced Waves in Nonlocal Dispersal KPP Equations with Shifting Habitats. Journal of Mathematical Analysis and Applications, 505, Article ID: 125504. >https://doi.org/10.1016/j.jmaa.2021.125504 Guo, J., Poh, A.A.L. and Wu, C. (2023) Forced Waves of Saturation Type for Fisher-KPP Equation in a Shifting Environment. Applied Mathematics Letters, 140, Article ID: 108573. >https://doi.org/10.1016/j.aml.2023.108573 Guo, J., Guo, K. and Shimojo, M. (2024) Uniqueness and Stability of Forced Waves for the Fisher-KPP Equation in a Shifting Environment. Nonlinear Analysis, 247, Article ID: 113607. >https://doi.org/10.1016/j.na.2024.113607 Meng, Y., Yu, Z. and Zhang, L. (2024) Existence, Uniqueness and Stability of Forced Waves for Asymptotical KPP Equations with the Nonlocal Dispersal in a Shifting Habitat. Discrete and Continuous Dynamical Systems—B, 29, 2382-2398. >https://doi.org/10.3934/dcdsb.2023182 Hu, H., Yi, T. and Zou, X. (2019) On Spatial-Temporal Dynamics of a Fisher-KPP Equation with a Shifting Environment. Proceedings of the American Mathematical Society, 148, 213-221. >https://doi.org/10.1090/proc/14659 Yi, T., Chen, Y. and Wu, J. (2020) Asymptotic Propagations of Asymptotical Monostable Type Equations with Shifting Habitats. Journal of Differential Equations, 269, 5900-5930. >https://doi.org/10.1016/j.jde.2020.04.025 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 Yuan, Y., Wang, Y. and Zou, X. (2019) Spatial Dynamics of a Lotka-Volterra Model with a Shifting Habitat. Discrete and Continuous Dynamical Systems—B, 24, 5633-5671. Wu, C. and Xu, Z. (2021) Propagation Dynamics in a Heterogeneous Reaction-Diffusion System under a Shifting Environment. Journal of Dynamics and Differential Equations, 35, 493-521. >https://doi.org/10.1007/s10884-021-10018-0 Wang, H., Pan, C. and Ou, C. (2021) Existence, Uniqueness and Stability of Forced Waves to the Lotka‐Volterra Competition System in a Shifting Environment. Studies in Applied Mathematics, 148, 186-218. >https://doi.org/10.1111/sapm.12438 Wang, H., Pan, C. and Ou, C. (2020) Existence of Forced Waves and Gap Formations for the Lattice Lotka-Volterra Competition System in a Shifting Environment. Applied Mathematics Letters, 106, Article ID: 106349. >https://doi.org/10.1016/j.aml.2020.106349 Qiao, S., Zhu, J. and Wang, J. (2021) Asymptotic Behaviors of Forced Waves for the Lattice Lotka-Volterra Competition System with Shifting Habitats. Applied Mathematics Letters, 118, Article ID: 107168. >https://doi.org/10.1016/j.aml.2021.107168 Meng, Y., Yu, Z. and Zhang, S. (2021) Spatial Dynamics of the Lattice Lotka-Volterra Competition System in a Shifting Habitat. Nonlinear Analysis: Real World Applications, 60, Article ID: 103287. >https://doi.org/10.1016/j.nonrwa.2020.103287 Zhu, J., Wang, J. and Dong, F. (2022) Spatial Propagation for the Lattice Competition System in Moving Habitats. Zeitschrift für angewandte Mathematik und Physik, 73, Article No. 92. >https://doi.org/10.1007/s00033-022-01735-7 Guo, J., Guo, K. and Shimojo, M. (2023) Forced Waves for Diffusive Competition Systems in Shifting Environments. Nonlinear Analysis: Real World Applications, 73, Article ID: 103880. >https://doi.org/10.1016/j.nonrwa.2023.103880 Huang, B. and Dai, B. (2024) Spatial Dynamics of a Lattice Lotka-Volterra Competition Model with a Shifting Habitat. Journal of Nonlinear Modeling and Analysis, 6, 161-183. Wu, J. and Zou, X. (2001) Traveling Wave Fronts of Reaction-Diffusion Systems with Delay. Journal of Dynamics and Differential Equations, 13, 651-687. >https://doi.org/10.1023/a:1016690424892
Baidu
map