aam Advances in Applied Mathematics 2324-7991 2324-8009 beplay体育官网网页版等您来挑战! 10.12677/aam.2024.139417 aam-96918 Articles 数学与物理 一类移动环境下Fisher-KPP方程受迫行波的存在性及其渐近行为
The Existence and Asymptotic Behavior of Forced Traveling Waves for a Class of Fisher-KPP Equation under a Shifting Habitat
陈碧霞 长沙理工大学数学与统计学院,湖南 长沙 09 09 2024 13 09 4370 4380 23 8 :2024 17 8 :2024 17 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/ 考虑在移动环境下Fisher-KPP方程受迫行波的存在性及其渐近行为,并假设此方程的内禀增长率函数恒大于某正常数。利用单调迭代的技巧证明了方程的非减受迫行波和非负受迫行波的存在性,进一步研究了两种受迫波的渐近行为。
In this paper, we consider the existence and the asymptotic behavior of the Fisher-KPP equation in the shifting habitat, and assume that the intrinsic growth rate function of this equation is always greater than a normal number. Using the technique of monotone iteration to prove the existence of non-decreasing and non-negative forced waves of the equation, we further study the asymptotic behavior of two forced waves.
移动环境,Fisher-KPP方程,单调迭代,渐近行为
Shifting Habitat
Fisher-KPP Equation Monotonic Iteration Asymptotic Behavior
1. 引言

在生态种群动力学的研究中,诸如种群的数量变化、渐近传播以及物种入侵等问题都可以通过建立合适的反应扩散方程模型来进行研究 [1] - [5] 。许多应用数学学者 [6] - [13] 利用反应扩散方程来描述环境变化对物种持续性的影响。特别地,Li等人 [6] 建立如下反应扩散方程

u ( t , x ) t = d 2 u ( t , x ) x 2 + u ( t , x ) [ r ( x c t ) u ( t , x ) ] , t > 0 , x , (1)

来研究在气候变化的影响下物种在移动栖息地灭绝和持久生存的条件以及渐近传播速度。其中 u ( t , x ) 表示种群在时间t位置x的密度大小,常数 d > 0 表示种群扩散系数, r ( ) 表示内禀增长率函数,以速度 c > 0 变化。Li等人假设 r ( ) 满足非减连续且 r ( ) < 0 < r ( + ) 的条件。Hu和Zou [8] 证明了对任意给定的 c > 0 ,方程(1)均存在一个连接0和 r ( + ) 的非减受迫行波 u ( t , x ) = φ ( x c t ) 。其中,受迫行波是指其波速与环境移动速度c相同的波。Berestycki和Fang [9] 考虑了方程(1)中的反应项为 f ( x c t , u ) 的情形,运用PDE的理论方法证明了强迫行波的完全存在性和多重性以及它们的吸引性,并考虑了行波在 + 处的衰减速度。Yang等人 [12] 还讨论了气候变化下局部扩散Lotka-Volterra合作系统受迫行波的存在性及其渐近行为。此外,Hu等人 [13] 把内禀增长率函数的条件减弱为: r ( ) 非减连续且满足 r ( ) 0 < r ( + ) ,并研究了移动环境下非局部扩散Lotka-Volterra型合作系统的受迫行波的存在性。

受上述研究工作的启发,我们将探讨种群受气候变化影响较弱或种群栖息地环境轻微恶化的情形下局部扩散Fisher-KPP方程受迫行波的存在性及其渐近行为,且假设如下条件成立:

(A1) 函数 r ( ) 上非减有界且连续,并满足 r ( + ) > r ( ) > 0

(A2) 函数 r ( ) 上连续可微,且 r ( ± ) 存在。

首先将方程(1)的受迫行波解记作 φ ( ξ ) , ξ : = x c t 。那么,方程(1)可转化为如下非自治常微分方程

c φ ( ξ ) = d φ ( ξ ) + φ ( ξ ) [ r ( ξ ) φ ( ξ ) ] . (2)

此外,方程(2)对应的极限方程分别为

c φ ( ξ ) = d φ ( ξ ) + φ ( ξ ) [ r ( + ) φ ( ξ ) ] (3)

c φ ( ξ ) = d φ ( ξ ) + φ ( ξ ) [ r ( ) φ ( ξ ) ] (4)

注意到 r ( + ) > r ( ) > 0 ,那么极限方程(3)和(4)分别有正平衡点 r ( + ) r ( ) 。因此,我们将考虑连接两个正平衡点之间行波的存在性,即探讨对任意 c > 0 方程(2)是否存在满足边界条件

lim ξ φ ( ξ ) = r ( ) , lim ξ + φ ( ξ ) = r ( + ) (5)

的非减的解。又由于 r ( ) > 0 ,我们还将在一定条件下探寻方程(2)满足边界条件

lim ξ φ ( ξ ) = r ( ) , lim ξ + φ ( ξ ) = 0 (6)

的非负解。为达到研究目的,我们还需要如下假设条件:

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

我们的研究结果将表明在种群受气候变化影响较弱或环境恶化程度较轻时(即物种内禀增长率恒正时),种群在任何一个区域内都能持久生存( φ ( ) = r ( ) > 0 )。

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

2. 预备知识

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

B C ( , ) = { u C ( , ) | sup ξ | u ( ξ ) | < } ,

对任意的 u , v C ( , ) ,如果 u v C + ,我们记 u v v u

α = 2 r ( + ) ,则方程 d λ 2 c λ + α = 0 有如下两个实根:

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

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

Δ h : = d h c h + α h ,

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

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

定义2.1 若 φ ¯ ( ξ ) , φ _ ( ξ ) B C ( , ) 满足 φ ¯ , φ _ , φ ¯ , φ _ L ( , ) φ ¯ φ _ φ ¯ φ _ \ { ξ j } 上连续( ξ j 为一有限递增点列),且有 φ ¯ ( ξ j + ) φ ¯ ( ξ j ) φ _ ( ξ j + ) φ _ ( ξ j ) ,此时,若有不等式

c φ ¯ ( ξ ) d φ ¯ ( ξ ) + φ ¯ ( ξ ) [ r ( ξ ) φ ¯ ( ξ ) ] (7)

c φ _ ( ξ ) d φ _ ( ξ ) + φ _ ( ξ ) [ r ( ξ ) φ _ ( ξ ) ] (8)

\ { ξ j } 上成立,则称 φ ¯ ( ξ ) φ _ ( ξ ) 为方程(2)的一对有序上下解。

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

Γ = { φ | φ B C ( , ) , φ _ φ φ ¯ } .

对任意 φ Γ ,定义如下算子:

H ( φ ) ( ξ ) : = α φ ( ξ ) + φ ( ξ ) [ r 1 ( ξ ) φ ( ξ ) ] .

接下来,我们定义算子F,其中 F = Δ 1 H 。下面我们先讨论算子F的一些性质。

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

证明 一方面,取 φ ˜ , φ ^ Γ 且满足 φ ˜ φ ^ ,则对任何 ξ

H ( φ ˜ ) ( ξ ) H ( φ ^ ) ( ξ ) = ( α + r ( ξ ) φ ˜ ( ξ ) φ ^ ( ξ ) ) ( φ ˜ ( ξ ) φ ^ ( ξ ) ) 0.

从而,对任意的 ξ ,有 F ( φ ˜ ) ( ξ ) F ( φ ^ ) ( ξ ) 。即F是一个非减算子。

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

φ _ F ( φ ) φ ¯ .

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

F ( φ _ ) = Δ 1 H ( φ _ ) Δ 1 Δ ( φ _ ) = φ _ , F ( φ ¯ ) = Δ 1 H ( φ ¯ ) Δ 1 Δ φ ¯ = φ ¯ . (9)

又由F是一个非减算子,所以对任意的 φ Γ

F ( φ _ ) F ( φ ) F ( φ ¯ ) . (10)

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

引理2.2 若 φ ( ξ ) Γ 非减,则 F ( φ ) ( ξ ) 关于 ξ 也非减。

证明 若 φ ( ξ ) Γ 是一个关于 ξ 的非减函数,则对任意 s ζ > 0 ,有

H ( φ ) ( s + ζ ) H ( φ ) ( s ) = ( φ ( s + ζ ) φ ( s ) ) [ α + r ( s + ζ ) φ ( s + ζ ) φ ( s ) ] + ( r ( s + ζ ) r ( s ) ) φ ( s ) 0.

对于任意的 ξ ,注意到当h为连续有界函数时,有 ( Δ 1 h ( s + ζ ) ) ( ξ ) = ( Δ 1 h ( s ) ) ( ξ + ζ ) ,那么

F ( φ ) ( ξ + ζ ) = [ Δ 1 H ( φ ) ( s ) ] ( ξ + ζ ) = [ Δ 1 H ( φ ) ( s + ζ ) ] ( ξ ) [ Δ 1 H ( φ ) ( s ) ] ( ξ ) = F ( φ ) ( ξ ) .

证毕。

方程(2)可写为

Δ φ = H ( φ ) (11)

因此,若映射F在Γ中存在一个不动点,即存在 φ Γ 使得

φ = F ( φ )

成立,那么该不动点必是(11)的解。若该不动点还满足边界条件(5)或(6),则必是方程(2)的受迫行波。因此,我们通过选取两对恰当的有序上下解构造出两个先验集。

引理2.3 若 φ ( ξ ) Γ 非减,则 F ( φ ) ( ξ ) 关于 ξ 也非减。设 φ ¯ 1 ( ξ ) = r ( + ) φ _ 1 ( ξ ) = r ( ) ,则 φ ¯ 1 ( ξ ) φ _ 1 ( ξ ) 是(2)的一对有序上下解。

证明 由 r ( ξ ) 的单调性和 φ ¯ 1 ( ξ ) 的定义可知

c ( r ( + ) ) d 1 ( r ( + ) ) + r ( + ) [ r 1 ( ξ ) r ( + ) ] ,

φ ¯ 1 ( ξ ) 满足定义2.1的所有条件,因此 φ ¯ 1 ( ξ ) = r ( + ) 是(2)的上解。

同理可得 φ _ 1 ( ξ ) = r ( ) 是(2)的下解。证毕。

c * ( ) : = 2 d r ( + ) .

接下来,我们在 c > c * ( ) 的条件下构造另一对有序上下解。方程(2)在原点处的线性化方程对应的特征方程为

Φ ( c , λ ) : = d λ 2 c λ + r ( + ) = 0.

容易看出,当 c = c * ( ) Φ ( c , λ ) = 0 有唯一正根 λ 0 ;当 c > c * ( ) Φ ( c , λ ) = 0 有两个正根。方便起见,我们记最小正根为 λ 1 。由假设(A*)可知,存在 K > 0 使得当 ξ K r ( + ) r ( ξ ) e ρ ξ 成立。因此,对于正数 ω min { ρ , λ 0 λ 1 } ,当 ξ K 时有

r ( + ) r ( ξ ) e ρ ξ e ω ξ .

此外,当 c > c * ( ) Φ ( c , λ 1 + ω ) < 0 ,令 M : = 2 Φ ( c , λ 1 + ω ) ,可取 ω 足够小使得 1 ω ln 1 M > max { K , ln r ( + ) λ 1 } 成立。定义如下两个有界连续函数

φ ¯ 2 ( ξ ) = { r ( + ) , ξ ξ 1 , e λ 1 ξ , ξ > ξ 1 ,

φ _ 2 ( ξ ) = { e λ 1 ξ 2 ( 1 M e ω ξ 2 ) , ξ ξ 2 , e λ 1 ξ ( 1 M e ω ξ ) , ξ > ξ 2 ,

其中 ξ 1 = 1 λ 1 ln r ( + ) ξ 2 ( 1 ω ln 1 M , 1 ω ln λ 1 ( λ 1 + ω ) M ) 使得 e λ 1 ξ 2 ( 1 M e ω ξ 2 ) ( 0 , r ( ) )

引理2.4 假设 c > c * ( ) 且(A*)成立,则 φ ¯ 2 ( ξ ) φ _ 2 ( ξ ) 是(2)的一对有序上下解。

证明 若 ξ > ξ 1 ,则 φ ¯ 2 ( ξ ) = e λ 1 ξ ,计算可得

d λ 1 2 e λ 1 ξ c λ 1 e λ 1 ξ + e λ 1 ξ [ r ( ξ ) e λ 1 ξ ] = e λ 1 ξ [ d λ 1 2 c λ 1 + r ( ξ ) e λ c ξ ] e λ 1 ξ [ d λ 1 2 c λ 1 + r ( + ) ] = 0.

这表明 φ ¯ 2 ( ξ ) ξ > ξ 1 时使得不等式(7)成立。当 ξ < ξ 1 时,易证 φ ¯ 2 ( ξ ) = r ( + ) 也满足不等式(7)。此外,易知 φ ¯ 2 ( ξ 1 + ) φ ¯ 2 ( ξ 1 ) 。因此 φ ¯ 2 ( ξ ) 是方程(2)的一个上解。

下面证明 φ _ 2 ( ξ ) 是方程(2)的一个下解。当 ξ > ξ 2 时,有 φ _ 2 ( ξ ) = e λ 1 ξ ( 1 M e ω ξ ) φ _ 2 ( ξ ) = λ 1 e λ 1 ξ + ( λ 1 + ω ) M e ( λ 1 + ω ) ξ φ _ 2 ( ξ ) = λ 1 2 e λ 1 ξ ( λ 1 + ω ) 2 M e ( λ 1 + ω ) ξ 。于是

d φ _ 2 ( ξ ) + c φ _ 2 ( ξ ) + φ _ 2 ( ξ ) [ r ( ξ ) φ _ 2 ( ξ ) ] = [ d λ 1 2 c λ 1 + r ( ξ ) e λ 1 ξ ( 1 M e ω ξ ) ] e λ 1 ξ [ d ( λ 1 + ω ) 2 c ( λ 1 + ω ) + r ( ξ ) e λ 1 ξ ( 1 M e ω ξ ) ] M e ( λ 1 + ω ) ξ = [ r ( + ) + r ( ξ ) e λ 1 ξ ( 1 M e ω ξ ) ] e λ 1 ξ [ Φ ( c , λ 1 + ω ) r ( + ) + r ( ξ ) e λ 1 ξ ( 1 M e ω ξ ) ] M e ( λ 1 + ω ) ξ e λ 1 ξ [ M e ω ξ Φ ( c , λ 1 + ω ) e ω ξ e λ 1 ξ ] e λ 1 ξ [ M e ω ξ Φ ( c , λ 1 + ω ) 2 e ω ξ ] = 0.

这表明 φ _ 2 ( ξ ) ξ > ξ 2 时使得不等式(8)成立。当 ξ < ξ 2 时,显然 φ _ 2 ( ξ ) = e λ 1 ξ 2 ( 1 M e ω ξ 2 ) 也满足不等式(8)。此外,易知 φ _ 2 ( ξ 2 + ) φ _ 2 ( ξ 2 ) ,而且根据 φ ¯ 2 ( ξ ) φ _ 2 ( ξ ) 的具体表达式可知 φ ¯ 2 ( ξ ) φ _ 2 ( ξ ) 。因此 φ ¯ 2 ( ξ ) φ _ 2 ( ξ ) 为方程(2)的一对有序上下解。证毕。

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

Γ 1 : = { φ | φ B C ( , ) , φ _ 1 φ φ ¯ 1 } ,

Γ 2 : = { φ | φ B C ( , ) , φ _ 2 φ φ ¯ 2 } .

3. 行波解的存在性

定理3.1 若(A1)成立,则对任意 c > 0 ,方程(2)总存在一个满足边界条件(5)的非减受迫行波。

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

φ 1 ( 1 ) = F ( φ ¯ 1 ) , φ 1 ( n + 1 ) = F ( φ 1 ( n ) ) , n 1.

由于 φ ¯ 1 ( ξ ) Γ 1 上的非减函数,结合引理2.1和引理2.2得到对所有的 n 1 φ ¯ 1 ( ξ ) 也是 上的非减函数且满足不等式

φ ¯ 1 ( ξ ) φ 1 ( 1 ) ( ξ ) φ 1 ( 2 ) ( ξ ) φ 1 ( n ) ( ξ ) φ 1 ( n + 1 ) ( ξ ) φ _ 1 ( ξ ) .

从而,存在一个有界非减函数 φ 1 ( ξ ) 使得 lim n φ 1 ( n ) ( ξ ) = φ 1 ( ξ ) 。很容易看出,对所有的 n 1 ξ

| H ( U 1 ( n ) ) ( ξ ) | r ( + ) ( 2 r ( + ) + α ) .

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

φ 1 ( ξ ) = lim n φ 1 ( n + 1 ) ( ξ ) = lim n F ( φ 1 ( n ) ) ( ξ ) = lim n ( Δ 1 H ( φ 1 ( n ) ) ( ξ ) ) = 1 d ( λ + λ ) [ ξ e λ ( ξ η ) H ( φ 1 ) ( η ) d η + ξ + e λ + ( ξ η ) H ( φ 1 ) ( η ) d η ] = F ( φ 1 ) ( ξ ) .

φ 1 ( ξ ) Γ 1 是算子F的不动点,也就是说 φ 1 ( ξ ) 是方程(2)的解。

接下来证明 φ 1 ( ξ ) 满足边界条件(5)。因 φ 1 ( ξ ) 上的有界非减函数,记 A 1 : = lim ξ φ 1 ( ξ ) B 1 : = lim ξ + φ 1 ( ξ ) 。显然有

0 < r ( ) A 1 r ( + ) , 0 < r ( ) B 1 r ( + ) .

lim ξ H ( φ 1 ) ( ξ ) = A 1 ( α + r ( ) A 1 ) , lim ξ + H ( φ 1 ) ( ξ ) = B 1 ( α + r ( + ) B 1 ) ,

利用L’Hôpital法则可得

A 1 = lim ξ φ 1 ( ξ ) = lim ξ ( Δ 1 H ( φ 1 ) ( ξ ) ) = lim ξ 1 d ( λ + λ ) [ ξ e λ ( ξ η ) H ( φ 1 ) ( η ) d η + ξ + e λ + ( ξ η ) H ( φ 1 ) ( η ) d η ]

= lim ξ 1 d ( λ + λ ) ( H ( φ 1 ) ( ξ ) λ + H ( φ 1 ) ( ξ ) λ + ) = A 1 + A 1 [ r ( ) A 1 ] α .

于是 A 1 = r ( ) 。同理可得

B 1 = lim ξ + φ 1 ( ξ ) = B 1 + B 1 [ r ( + ) B 1 ] α .

所以 B 1 = r ( + ) 。证毕。

定理3.2 若(A1)和(A*)成立,则当 c > c * ( ) 时方程(1)存在满足边界条件(6)的非负受迫行波。

证明 类似于定理3.1的证明过程可得行波的存在性。由于此时 φ 2 ( ξ ) Γ 2 上的连续有界函数但不具备单调性,因此方程(2)的解 φ 2 ( ξ ) 也不具备单调性。

接下来证明 φ 2 ( ξ ) 满足边界条件(6)。根据 φ _ 2 ( ξ ) φ ¯ 2 ( ξ ) 的定义可知 φ 2 ( ξ ) 为正解,且有

lim ξ + φ _ 2 ( ξ ) = 0 , lim ξ + φ ¯ 2 ( ξ ) = 0.

于是 lim ξ + φ 2 ( ξ ) = 0

下面我们证明 lim ξ φ 2 ( ξ ) = r ( ) 。记

A 2 : = lim sup ξ φ 2 ( ξ ) , B 2 : = lim inf ξ φ 2 ( ξ ) .

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

lim n φ 2 ( s n ) = A 2 , lim n φ 2 ( t n ) = B 2 ,

lim n φ 2 ( s n ) = lim n φ 2 ( t n ) = 0.

根据 φ 2 ( ξ ) = F ( φ 2 ) ( ξ ) 可得

φ 2 ( ξ ) = λ + φ 2 ( ξ ) 1 d ξ e λ ( ξ η ) H ( φ 2 ) ( η ) d η .

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

0 < φ 2 ( η ) < A 2 + ϵ , r ( ) < r ( η ) < r ( ) + ϵ .

因此,当 n > N 1 时有

φ 2 ( s n ) = λ + φ 2 ( s n ) 1 d s n e λ ( s n η ) H ( φ 2 ) ( η ) d η λ + φ 2 ( s n ) 1 d s n e λ ( s n η ) ( A 2 + ϵ ) ( α + r ( ) A 2 ) d η = λ + φ 2 ( s n ) + 1 d λ ( A 2 + ϵ ) ( α + r ( ) A 2 ) .

n ,我们有

λ + A 2 + 1 d λ ( A 2 + ϵ ) ( α + r ( ) A 2 ) 0.

ϵ 的任意性得

λ + A 2 + 1 d λ A 2 ( α + r ( ) A 2 ) 0.

注意到 A 2 > 0 ,上不等式可推得

A 2 r ( ) . (12)

类似地,对任意 ϵ ( 0 , B 2 ) N 2 > 0 使得当 η ( , t N 2 ) 时有

B 2 ϵ < φ 2 ( η ) r ( + ) , r ( ) < r ( η ) < r ( ) + ϵ .

因此,当 n > N 2 时有

φ 2 ( t n ) = λ + φ 2 ( t n ) 1 d t n e λ ( t n η ) H ( φ 2 ) ( η ) d η λ + φ 2 ( t n ) 1 d t n e λ ( t n η ) ( B 2 ϵ ) ( α + r ( ) B 2 + ϵ ) d η = λ + φ 2 ( t n ) + 1 d λ ( B 2 ϵ ) ( α + r ( ) B 2 + ϵ ) .

n ,我们有

λ + B 2 + 1 d λ ( B 2 ϵ ) ( α + r ( ) B 2 + ϵ ) 0.

ϵ 的任意性得

λ + B 2 + 1 d λ B 2 ( α + r ( ) B 2 ) 0.

注意到 B 2 > 0 ,上不等式可推得

B 2 r ( ) . (13)

注意到 B 2 A 2 ,由(12)和(13)可知 A 2 = B 2 = r ( ) 。因此, lim ξ φ 2 ( ξ ) = r ( ) 。证毕。

4. 行波解的渐近行为

为了获得更多的在 ± 处的行波的渐近信息,我们令

φ ( ξ ) = u 0 ( ξ ) for < ξ < + ,

它是满足边界条件(5)或(6)的方程(2)的解,对(2)关于 ξ 求导,那么我们知道 φ ( ξ ) = : θ 满足

d θ ( ξ ) + c θ ( ξ ) + [ r ( ξ ) 2 u 0 ] θ + r ( ξ ) u 0 = 0. (14)

定理4.1 假设(A1)和(A2)成立,那么存在正常数 P 1 , P 2 使得方程(1)的行波具有以下渐近性质当 ξ 时有

φ 1 ( ξ ) = r ( ) + ( P 1 + o ( 1 ) ) e c + c 2 + 4 d r ( ) 2 d ξ .

且当 ξ + 时有

φ 1 ( ξ ) = r ( + ) ( P 2 + o ( 1 ) ) e c c 2 + 4 d r ( + ) 2 d ξ .

证明 当 ξ 时方程(14)的极限方程可以写为

d ϕ ( ξ ) + c ϕ ( ξ ) r ( ) ϕ ( ξ ) = 0. (15)

注意到在(15)中我们用到了 r ( ) = 0 。事实上,应用L’Hôpital’s法则,我们可以得到

lim ξ r ( ξ ) = lim ξ e ξ r ( ξ ) e ξ = lim ξ e ξ [ r ( ξ ) + r ( ξ ) ] e ξ = lim ξ [ r ( ξ ) + r ( ξ ) ] ,

这表明 r ( ) = 0 。同理可得 r ( + ) = 0

方程(15)有两个形式如下的线性无关的解

ϕ 1 ( ξ ) = e c c 2 + 4 d r ( ) 2 d ξ , ϕ 2 ( ξ ) = e c + c 2 + 4 d r ( ) 2 d ξ .

联立(14)和(15),我们发现当 ξ θ 具有如下性质

θ ( ξ ) = p 1 [ 1 + o ( 1 ) ] ϕ 1 ( ξ ) + q 1 [ 1 + o ( 1 ) ] ϕ 2 ( ξ ) .

由于 lim ξ θ ( ξ ) = 0 ,所以必须有 p 1 = 0 。因此,对于 ξ

φ 1 ( ξ ) = θ ( ξ ) = q 1 [ 1 + o ( 1 ) ] e c + c 2 + 4 d r ( ) 2 d ξ .

对上式从 ξ 积分,可以得到当 ξ 时,对某个常数 P 1 > 0

φ 1 ( ξ ) = r ( ) + ( P 1 + o ( 1 ) ) e c + c 2 + 4 d r ( ) 2 d ξ .

类似地,当 ξ + 时,可以得到对某个常数 P 2 > 0

φ 1 ( ξ ) = r ( + ) ( P 2 + o ( 1 ) ) e c c 2 + 4 d r ( + ) 2 d ξ .

证毕。

定理4.2 假设(A1)、(A2)和(A*)成立,那么存在正常数 Q 1 , Q 2 使得方程(1)的行波具有以下渐近性质。当 ξ 时有

φ 2 ( ξ ) = r ( ) ( Q 1 + o ( 1 ) ) e c + c 2 + 4 d r ( ) 2 d ξ .

且当 ξ + 时有

φ 2 ( ξ ) = ( Q 2 + o ( 1 ) ) e c + c 2 4 d r ( + ) 2 d ξ .

证明 运用定理4.1的证明方法可得当 ξ 时,对某个常数 C 1

φ 2 ( ξ ) = r ( ) + ( C 1 + o ( 1 ) ) e c + c 2 + 4 d r ( ) 2 d ξ .

由于对任意 ξ lim ξ φ 2 ( ξ ) = r ( ) φ 2 ( ξ ) r ( ) 。我们可以得到常数 C < 0 ,记 Q 1 : = C > 0 。由此可得

φ 2 ( ξ ) = r ( ) ( Q 1 + o ( 1 ) ) e c + c 2 + 4 d r ( ) 2 d ξ .

ξ + 时,(14)的极限方程如下

d ψ ( ξ ) + c ψ ( ξ ) + r ( + ) ψ ( ξ ) = 0. (16)

方程(16)有两个形式如下的独立的解

ψ 1 ( ξ ) = e c c 2 4 d r ( + ) 2 d ξ , ψ 2 ( ξ ) = e c + c 2 4 d r ( + ) 2 d ξ .

类似地,我们可以得到当 ξ + 时,对某个常数 D 1 , D 2

φ 2 ( ξ ) = ( D 1 + o ( 1 ) ) e c c 2 4 d r ( ) 2 d ξ ( D 2 + o ( 1 ) ) e c + c 2 4 d r ( ) 2 d ξ .

进一步可得对某个常数D有

φ 2 ( ξ ) = ( D + o ( 1 ) ) e c + c 2 4 d r ( ) 2 d ξ .

由于对任意 ξ lim ξ + φ 2 ( ξ ) = 0 φ 2 ( ξ ) 0 。我们可以得到常数 D < 0 ,记 Q 2 : = D > 0 ,由此可得

φ 2 ( ξ ) = ( Q 2 + o ( 1 ) ) e c + c 2 4 d r ( ) 2 d ξ .

证毕。

5. 总结

本文主要研究了一类移动环境下随机扩散Fisher-KPP方程受迫波的存在性及其渐近行为,主要通过构造上下解和先验集,运用单调迭代技巧证明了该方程非减受迫行波和非负受迫行波的存在性。当环境移动速度 c > 0 时,证明了方程(1)存在一个连接 r ( ) r ( + ) 两个正平衡点的非减受迫行波。当 c > c * ( ) 时,证明了方程(1)存在一个连接 r ( ) 和0的非负受迫行波。并研究了这两种行波在 ± 处的渐近行为。同时注意到这两个受迫波的边界条件中都有 lim ξ φ ( ξ ) = r ( ) 。由于 r ( ) 恒正, r ( ) > 0 ,这表明种群栖息地环境是轻微恶化的,因此可知该种群在任一固定的区域上都可以持久生存。

References 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 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. Cantrell, R.S. and Cosner, C. (2004) Spatial Ecology via Reaction‐Diffusion Equations. Wiley. >https://doi.org/10.1002/0470871296 倪维明. 浅谈反应扩散方程[J]. 数学传播, 2016, 34(4): 17-26. 楼元. 空间生态学中的一些反应扩散方程模型[J]. 中国科学: 数学, 2015, 45(10): 1619-1634. 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 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 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 Li, W., Wang, J. and Zhao, X. (2018) Spatial Dynamics of a Nonlocal Dispersal Population Model in a Shifting Environment. Journal of Nonlinear Science, 28, 1189-1219. >https://doi.org/10.1007/s00332-018-9445-2 Wu, C., Wang, Y. and Zou, X. (2019) Spatial-Temporal Dynamics of a Lotka-Volterra Competition Model with Nonlocal Dispersal under Shifting Environment. Journal of Differential Equations, 267, 4890-4921. >https://doi.org/10.1016/j.jde.2019.05.019 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 Hu, H., Deng, L. and Huang, J. (2021) Traveling Wave of a Nonlocal Dispersal Lotka-Volterra Cooperation Model under Shifting Habitat. Journal of Mathematical Analysis and Applications, 500, Article 125100. >https://doi.org/10.1016/j.jmaa.2021.125100 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 Smith, H.L. (2011) An Introduction to Delay Differential Equations with Applications to the Life Sciences. Springer.
Baidu
map