APP下载

非局部扩散三种群竞争系统的双稳行波解的Lyapunov稳定性

2022-06-06张国宝郝玉财

关键词:行波波速平衡点

张国宝,郝玉财

(西北师范大学 数学与统计学院,甘肃 兰州 730070)

研究具有非局部扩散的三种群Lotka-Volterra竞争系统[1-2]

(1)

其中,u(x,t),v(x,t),w(x,t)分别表示三个种群在空间位置x和时刻t的种群密度,bij>0是种群j和种群i之间的竞争系数,ri>0是种群i的自然增长率,di>0是种群i的扩散系数,Ji是种群i的迁徙率.为了简单起见,通过时间和空间尺度变换,将各种群的环境承载力标准化为1.Ji*υ-υ称为非局部扩散算子,表示长距离的扩散机制,其中Ji*υ是R上的标准空间卷积,即

由系统(1)可见,种群v和种群w对食物资源有不同的偏好,即种群v和w之间没有竞争,但是,种群u与种群v,w都存在竞争关系.核函数Ji可以看作概率密度函数.假设核函数Ji和系统(1)的竞争系数满足以下条件:

(H1)Ji∈C(R),Ji(x)=Ji(-x)≥0,且

(H3)01.

条件(H3)表明种群v和w相对于种群u是弱竞争者,但将种群v和w放在一起可能会在与种群u的竞争中胜出.

非局部扩散系统存在于许多应用学科,如种群生态学、流行病学、相变以及神经网络模型等[3-7].系统(1)是如下具有随机扩散的Lotka-Volterra竞争系统的非局部扩散情形:

(2)

显然,系统(1)和系统(2)对应的空间齐次系统为

(3)

常微分方程系统(3)有一个平凡平衡点(0,0,0),4个半平凡平衡点(1,0,0),(0,1,0),(0,0,1),(0,1,1)以及一个正平衡点(u*,v*,w*),其中,

由条件(H3)可知,常数平衡点(0,1,1)和(1,0,0)是线性稳定的,其他平衡点是不稳定的,因此,该系统具有双稳非线性.假定目前的生态系统中已有种群v和w,现在引入新的种群u,一个很自然的问题是哪个种群将在竞争中获胜?为了回答这个问题,我们研究系统(1)的连接平衡点(0,1,1)和(1,0,0)的行波解.

系统(1)的行波解具有形式

(u,v,w)(x,t)=(φ,ψ,φ)(x+ct),

(4)

其中,x∈R,t≥0,c∈R是非零传播速度.将(4)式代入(1)式,可得系统(1)的波廓系统为

其中,ξ:=x+ct,且

基于对入侵者u的入侵过程的兴趣,(φ,ψ,φ)满足如下边界条件:

2017年,在b21>1,b31>1,b12+b13<1的假设下,Dong等[8]建立了系统(1)的连接平衡点(0,1,1)和(1,0,0)的单稳行波解的存在性和渐近行为,Zhang等[2]进一步研究了系统(1)的单稳行波解的唯一性和稳定性.最近,He等[9]考虑了系统(1)的行波解的最小波速,并给出了最小波速线性抉择的条件.特别地,Dong等[1]研究了系统(1)具有非局部各向异性扩散的双稳行波解,他们首先通过截断问题的极限论证建立了行波解的存在性和单调性,然后利用Ikehara引理得到了行波解的渐近行为,并在此基础上研究了波廓和波速的唯一性.到目前为止,系统(1)的双稳行波解的稳定性依然是未知的,因此,文中主要研究系统(1)的双稳行波解的稳定性.

1 预备知识和主要结论

(5)

其中,x∈R,t>0,且f1,f2,f3定义如下:

显然,在区域{(u,v,w):u≥0,v≤1,w≤1}中,f1,f2,f3满足如下条件:

同时,系统(5)对应的行波解(φ,ψ,φ)满足

(6)

且具有渐近边界条件

命题1[1]假定条件(H1)~(H3)成立,则系统(5)存在连接平衡点(0,0,0)和(1,1,1)的波速为c∈R的非降行波解(φ,ψ,φ)(x+ct),且上述双稳行波解在c≠0时是严格单调的.

定理1(Lyapunov稳定性) 假定条件(H1)~(H3)成立,(φ,ψ,φ)(x+ct)是系统(5)的双稳行波解,其中c≠0.设(u,v,w)(x,t)是系统(5)的具有初值(u,v,w)(x,0)=(u0,v0,w0)(x)的解,其中

则(φ,ψ,φ)(x+ct)是Lyapunov稳定的.即对任意的ε>0,存在常数γ>0,若

||(u0,v0,w0)(·)-(φ,ψ,φ)(·)||<γ

(9)

成立,则

2 Lyapunov稳定性

本节利用上下解方法研究系统(5)的具有非零波速的双稳行波解的Lyapunov稳定性.

定义1假定函数(u+,v+,w+)满足

(11)

其中,x∈R,t>0,则(u+,v+,w+)称为系统(5)的上解.类似地,通过改变不等号的方向可以定义系统(5)的下解.

设E-和E+分别是f1,f2,f3在(0,0,0)和(1,1,1)处的Jacobi矩阵,即

λ±=max{Reλ:λ∈σ(M±)}<0,

(12)

M±φ±=λ±φ±.

(13)

设E=∂(f1,f2,f3)/∂(u,v,w)是取值于(u,v,w)处的Jacobi矩阵.若取ν>0足够小,则

由(7)式可知,存在ξ1=ξ1(ν)>0,使得

定义函数

由(6)式可得

Γ(x,t)=σe-βt(βφ′(ξ)+δ(I1-I2)),

其中

由(18)和(19)式可知,存在常数ξ2>0,使得

接下来证明:对任意的(x,t)∈R×[0,+∞)和σ>0,若δ足够小,则Γ(x,t)≥0.令ξ0:=max{ξ1,ξ2+1},并且假设δ0≤ν/2σ.分以下三种情形进行讨论:

情形1η>ξ0.在此情形下,η>ξ2+1,因此,由(22)和(23)式可得

显然,ξ≥η>ξ1,所以由(17)式可知

由(13),(15)和(22)式可得

因此

情形2η<-ξ0.和情形1类似地讨论可知Γ(x,t)≥0.

显然L3仅依赖于v.进一步有I2≤3L3.

定义

因此Γ(x,t)≥0只需取δ0≤βL4/(|c|+L1+3L3).因此,(11)式的第一个不等式成立.

同理可证其他不等式也成立.】

定理1的证明令

再由(9)式可得

其中,x∈R,t≥0.因此,由比较原理可得

这说明|u(x,t)-φ(x+ct)|<2σL5.因此,由σ=ε/(2L5)可得

||u(·,t)-φ(·)||<2σL5=ε.

同理可证

综上,定理1自然成立.】

猜你喜欢

行波波速平衡点
2013-12-16巴东MS5.1地震前后波速比异常特征
带有超二次位势无限格点上的基态行波解
一类非局部扩散的SIR模型的行波解
用Riccati方程的新解求Fitzhugh-Nagumo方程的新行波解
土层剪切波速与埋深间的统计关系研究
基于实测波速探讨地震反射波法超前预报解译标志
一类带治疗项的非局部扩散SIR传染病模型的行波解
灰岩声波波速和力学参数之间的关系研究
探寻中国苹果产业的产销平衡点
电视庭审报道,如何找到媒体监督与司法公正的平衡点