留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

具有时滞的离散Lotka-Volterra合作系统波前解的非线性稳定性

闫瑞 刘桂荣 李晓翠

闫瑞, 刘桂荣, 李晓翠. 具有时滞的离散Lotka-Volterra合作系统波前解的非线性稳定性[J]. 应用数学和力学, 2023, 44(4): 461-470. doi: 10.21656/1000-0887.430172
引用本文: 闫瑞, 刘桂荣, 李晓翠. 具有时滞的离散Lotka-Volterra合作系统波前解的非线性稳定性[J]. 应用数学和力学, 2023, 44(4): 461-470. doi: 10.21656/1000-0887.430172
YAN Rui, LIU Guirong, LI Xiaocui. Nonlinear Stability of Traveling Wavefronts for a Discrete Cooperative Lotka-Volterra System With Delays[J]. Applied Mathematics and Mechanics, 2023, 44(4): 461-470. doi: 10.21656/1000-0887.430172
Citation: YAN Rui, LIU Guirong, LI Xiaocui. Nonlinear Stability of Traveling Wavefronts for a Discrete Cooperative Lotka-Volterra System With Delays[J]. Applied Mathematics and Mechanics, 2023, 44(4): 461-470. doi: 10.21656/1000-0887.430172

具有时滞的离散Lotka-Volterra合作系统波前解的非线性稳定性

doi: 10.21656/1000-0887.430172
基金项目: 

国家自然科学基金项目 11971279

国家自然科学基金项目 12101034

详细信息
    作者简介:

    闫瑞(1981—),女,副教授,硕士(E-mail: yanrui@sxufe.edu.cn)

    通讯作者:

    李晓翠(1982—),女,博士(通讯作者. E-mail: xiaocuili@mail.buct.edu.cn)

  • 中图分类号: O175.7

Nonlinear Stability of Traveling Wavefronts for a Discrete Cooperative Lotka-Volterra System With Delays

  • 摘要: 反应扩散模型的行波解的稳定性是一个很重要的研究课题.该文主要研究了一类具有时滞的离散Lotka-Volterra合作系统波前解的全局非线性稳定性.具体来讲, 当初值在无穷远处指数衰减到有较大波速的波前解而在其他位置可以任意大时, 运用L2-加权能量方法、比较原理和挤压技术可以得到该系统的此类波前解是指数渐近稳定的, 并解决了离散扩散算子及时滞共同作用下建立能量估计的问题.总之, 将加权能量方法推广到带有时滞的离散系统中, 丰富了相关的研究内容.
  • 由专著[1]可知在实际应用中, 许多物理学、化学、生物学、传染病学等问题都可以归结为反应扩散方程.所以, 研究反应扩散方程的相关内容对生产实践和科学理论具有重要的意义.近年来, 关于具有时滞的反应扩散系统的研究是一个重要的课题.因此, 本文主要研究如下具有时滞的离散Lotka-Volterra合作系统:

    $$ \left\{\begin{array}{l} u_t(x, t)=d_1 \mathcal{D}[u(x, t)]+r_1 u(x, t)\left[1-a_1 u(x, t)+b_1 v(x, t-\tau)\right], \\ v_t(x, t)=d_2 \mathcal{D}[v(x, t)]+r_2 v(x, t)\left[1+b_2 u(x, t-\tau)-a_2 v(x, t)\right], \end{array}\right. $$ (1)

    其中x∈$ \mathbb{R} $, t∈$ \mathbb{R}^{+}$, $ \mathcal{D} $[w(x, t)]=w(x+1, t)+w(x-1, t)-2w(x, t), w=u, v, 所有系数都是正的且时滞τ>0.这里u(x, t)和v(x, t)分别代表x位置和t时刻的人口密度.其初值条件为

    $$ \begin{cases}u(x, s)=u_0(x, s), & (x, s) \in \mathbb{R} \times[-\tau, 0], \\ v(x, s)=v_0(x, s), & (x, s) \in \mathbb{R} \times[-\tau, 0] .\end{cases} $$ (2)

    事实上, 系统(1)可以认为是如下系统的空间离散化形式:

    $$ \left\{\begin{array}{l} u_t(x, t)=d_1 u_{x x}(x, t)+r_1 u(x, t)\left[1-a_1 u(x, t)+b_1 v(x, t-\tau)\right], \\ v_t(x, t)=d_2 v_{x x}(x, t)+r_2 v(x, t)\left[1+b_2 u(x, t-\tau)-a_2 v(x, t)\right] . \end{array}\right. $$ (3)

    显然, 系统(1)有4个可能的平衡点

    $$ (0, 0), \left(0, \frac{1}{a_2}\right), \left(\frac{1}{a_1}, 0\right), \left(\frac{b_1+a_2}{a_1 a_2-b_1 b_2}, \frac{a_1+b_2}{a_1 a_2-b_1 b_2}\right):=\left(k_1, k_2\right) . $$

    易知, 当a1a2-b1b2>0时, ki>0(i=1, 2), 在系统(3)对应的无扩散系统中, 平衡点(0, 0)是不稳定的, 而平衡点(k1, k2)是稳定的.

    在离散方程和系统的研究中, 行波解的相关研究也是热点问题之一.对于系统(1) 的连接(0, 0)与(k1, k2)的行波解指的是具有如下形式:

    $$ (u(x, t), v(x, t))=(\phi(\xi), \psi(\xi)), \xi:=x+c t, $$

    且满足

    $$ \left\{\begin{array}{l} c \phi^{\prime}(\xi)=d_1 \mathcal{D}[\phi(\xi)]+r_1 \phi(\xi)\left[1-a_1 \phi(\xi)+b_1 \psi(\xi-c \tau)\right], \\ c \psi^{\prime}(\xi)=d_2 \mathcal{D}[\psi(\xi)]+r_2 \psi(\xi)\left[1+b_2 \phi(\xi-c \tau)-a_2 \psi(\xi)\right], \\ (\phi, \psi)(-\infty)=(0, 0), (\phi, \psi)(+\infty)=\left(k_1, k_2\right) \end{array}\right. $$

    的解, 其中D[φ(ξ)]:=φ(ξ+1)+φ(ξ-1)-2φ(ξ), φ=ϕ, ψ. 另外, 如果ϕψ是单调的, 称(ϕ(ξ), ψ(ξ))为一个波前解.

    关于不同类型的扩散系统的行波解存在性研究已有了大量的结论, 可参考文献[2-8].在关于行波解的研究中, 除存在性外, 行波解稳定性也是一个比较有趣和困难的问题.关于这方面的研究可参考文献[9-15].对于带有时滞的反应扩散方程, Schaaf首先通过谱分析方法研究了该类方程的行波解的线性稳定性[16].Mei等[17]则利用加权L2能量估计的方法证明了具有时滞的Nicholson方程波前解的非线性稳定性.接着, Mei等在文献[18-20]中应用加权能量法和比较原理,研究了一般的带有时滞的单稳反应扩散方程行波解的全局稳定性.之后, Yu等[21]和Zhang等[22]将加权能量法证明波前解稳定性的方法推广到不同的非局部扩散系统.此外, 对于离散扩散方程的行波解的稳定性的研究可以参考文献[23-26].对于离散扩散系统, 最近, Chen等[27]及Su等[28]分别运用加权能量结合比较原则证明了一个三种群离散扩散竞争系统的单稳波前解的非线性稳定性.而在文献[29]中, Hsu等运用不同的比较定理, 研究了一些离散反应扩散系统行波解的稳定性, 这些结论可以广泛应用到诸如多种群Lotka-Volterra合作模型、流行病模型、三种群Lotka-Volterra竞争模型等.但是, 目前关于系统(1)的波前解稳定性并没有任何的结论, 因此, 受文献[17-20]的启发, 本文将通过加权能量的方法研究该系统波前解的全局稳定性.

    本文的主要内容安排如下: 第1节给出了一些预备知识和符号, 并给出了波前解稳定性的主要结果;第2节给出了稳定性的证明;最后给出了本文的结论.

    在本节中, 我们将引入一些记号以方便叙述并给出主要结论.首先, 在本文中, C>0代表一个一般的常数, Ci>0代表具体的常数, I是一个区间, 通常I=$ \mathbb{R} $.其次, L2(I)表示I上的平方可积函数构成的空间, 若函数f(x)且其导数dif/dxi(i=1, 2, …, k)都属于L2(I), 则这些函数构成Sobolev空间Hk(I)(k≥0).进一步, Lω2(I)表示加权的L2空间(ω(x)>0), 其范数为

    $$ \|f\|_{L_\omega^2}=\left(\int_I \omega(x)|f(x)|^2 \mathrm{~d} x\right)^{1 / 2}, $$

    Hωk(I)为加权Sobolev空间, 其范数为

    $$ \|f\|_{H_\omega^k(I)}=\left(\sum\limits_{i=0}^k \int_I \omega(x)\left|\frac{\mathrm{d}^i}{\mathrm{~d} x^i} f(x)\right|^2 \mathrm{~d} x\right)^{1 / 2} \cdot $$

    T为一个正常数, $ \mathcal{B} $是一个Banach空间.C([0, T]; $ \mathcal{B} $)表示由[0, T]上的$ \mathcal{B} $-值连续函数构成的空间.L2([0, T]; $ \mathcal{B} $)表示由[0, T]上的$ \mathcal{B} $-值L2-函数构成的空间.[0, ∞)上的空间可以类似地给出.

    接下来, 对于λ>0, c>0, 定义

    $$ \Delta_i(c, \lambda)=d_i\left(\mathrm{e}^\lambda+\mathrm{e}^{-\lambda}-2\right)-c \lambda+r_i, \quad i=1, 2. $$

    引理1[30]   存在唯一常数c*>0使得

    1) 对于任意c>c*, 1(c, λ)=0有两个实根0 < λ1(c) < λ2(c), 且当λ∈(λ1(c), λ2(c)) 时, 1(c, λ) < 0, 当λ∈$ \mathbb{R} $+\[λ1(c), λ2(c)] 时, 1(c, λ)>0;

    2) 对于任意c>c*, 2(c, λ)=0有两个实根0 < λ3(c) < λ4(c), 且当λ∈(λ3(c), λ4(c))时, 2(c, λ) < 0, 当λ∈$ \mathbb{R} $+\[λ3(c), λ4(c)]时, 2(c, λ)>0;

    3) 对于任意0≤c < c*及所有λ>0, δi(c, λ)>0, i=1, 2.

    另外, 当a1a2-b1b2>0, d1d2, r1r2时, 对于c>c*, 系统(1)存在连接(0, 0)和(k1, k2)的波前解.对于0 < c < c*, 连接(0, 0)和(k1, k2)的波前解不存在.

    在文献[30]的基础之上, 本文将进一步研究系统(1)波前解的稳定性.首先给出如下假设:

    (A1) a1a2-b1b2>0, d1d2, r1r2, 2r1-d1>0, 2r2-d2>0;

    (A2) 4a1r1k1>2r1(1+2b1k2)+b1r1k1+b2r2k2+d1, 4a2r2k2>2r2(1+2b2k1)+b1r1k1+b2r2k2+d2.

    注1   令a1=2, a2=2.5, b1=1, b2=1, r1=1.5, r2=2, d1=0.1, d2=0.2, 通过计算可知假设(A2)是成立的.

    为了给出加权能量函数及适当的估计, 在此我们定义如下两个关于β的函数:

    $$ \begin{aligned} & f_1(\beta)=4 a_1 r_1 k_1-2 r_1\left(1+2 b_1 k_2\right)-b_1 r_1 k_1-b_2 r_2 k_2-d_1\left(\mathrm{e}^\beta+\mathrm{e}^{-\beta}-1\right), \\ & f_2(\beta)=4 a_2 r_2 k_2-2 r_2\left(1+2 b_2 k_1\right)-b_1 r_1 k_1-b_2 r_2 k_2-d_2\left(\mathrm{e}^\beta+\mathrm{e}^{-\beta}-1\right) . \end{aligned} $$

    由假设(A2), 易得

    $$ \begin{aligned} & f_1(0)=4 a_1 r_1 k_1-2 r_1\left(1+2 b_1 k_2\right)-b_1 r_1 k_1-b_2 r_2 k_2-d_1>0, \\ & f_2(0)=4 a_2 r_2 k_2-2 r_2\left(1+2 b_2 k_1\right)-b_1 r_1 k_1-b_2 r_2 k_2-d_2>0 . \end{aligned} $$

    再由fi(β)(i=1, 2)的连续性, 知存在β0>0,使得f1(β0)>0和f2(β0)>0.之后再定义另外两个关于ξ的函数:

    $$ \begin{aligned} & g_1(\xi)=-2 r_1\left(1+2 b_1 k_2\right)-b_1 r_1 k_1-b_2 r_2 k_2+4 a_1 r_1 \phi(\xi)-d_1\left(\mathrm{e}^{\beta_0}+\mathrm{e}^{-\beta_0}-1\right), \\ & g_2(\xi)=-2 r_2\left(1+2 b_2 k_1\right)-b_1 r_1 k_1-b_2 r_2 k_2+4 a_2 r_2 \psi(\xi)-d_2\left(\mathrm{e}^{\beta_0}+\mathrm{e}^{-\beta_0}-1\right), \end{aligned} $$

    其中(ϕ(ξ), ψ(ξ))是系统(1)的波前解.易得

    $$ \lim\limits _{\xi \rightarrow+\infty} g_1(\xi)=f_1\left(\beta_0\right)>0, \lim\limits _{\xi \rightarrow+\infty} g_2(\xi)=f_2\left(\beta_0\right)>0 . $$

    于是存在足够大的ξ0>0,使得

    $$ g_i\left(\xi_0\right)>0, \quad i=1, 2 . $$ (4)

    由上文, 我们进一步定义如下的关于β0ξ0的加权函数:

    $$ \omega(\xi)= \begin{cases}\mathrm{e}^{-\beta_0\left(\xi-\xi_0\right)}, & \xi \leqslant \xi_0, \\ 1, & \xi>\xi_0 . \end{cases} $$

    最后, 我们给出本文的主要结论, 即系统(1)的波前解的指数稳定性.

    定理1   假设(A1)和(A2)成立.令

    $$ \tilde{c}=\max \left\{c^*, \beta_0^{-1} c_1, \beta_0^{-1} c_2\right\}, $$

    其中

    $$ \begin{aligned} & c_1=2 r_1\left(1+2 b_1 k_2\right)+b_1 r_1 k_1+b_2 r_2 k_2+d_1\left(\mathrm{e}^{\beta_0}+\mathrm{e}^{-\beta_0}-1\right), \\ & c_2=2 r_2\left(1+2 b_2 k_1\right)+b_1 r_1 k_1+b_2 r_2 k_2+d_2\left(\mathrm{e}^{\beta_0}+\mathrm{e}^{-\beta_0}-1\right) . \end{aligned} $$

    c>$\tilde{c} $时, 对于系统(1)给定的波前解(ϕ(x+ct), ψ(x+ct)), 如果初始值及其扰动满足

    $$ 0 \leqslant u_0(x, s) \leqslant k_1, 0 \leqslant v_0(x, s) \leqslant k_2, \quad(x, s) \in \mathbb{R} \times[-\tau, 0], $$

    u0(x, s)-ϕ(x+cs)∈C([-τ, 0]; Hω1($ \mathbb{R} $))和v0(x, s)-ψ(x+cs)∈C([-τ, 0]; Hω1($ \mathbb{R} $)), 那么Cauchy问题(1)、(2)的解满足

    $$ \begin{aligned} & 0 \leqslant u(x, t) \leqslant k_1, 0 \leqslant v(x, t) \leqslant k_2, \quad(x, t) \in \mathbb{R} \times \mathbb{R}_{+}, \\ & u(x, t)-\phi(x+c t) \in C\left([0, +\infty) ; H_\omega^1(\mathbb{R})\right), v(x, t)-\psi(x+c t) \in C\left([0, +\infty) ; H_\omega^1(\mathbb{R})\right) . \end{aligned} $$

    另外, 解(u(x, t), v(x, t))指数收敛到(ϕ(x+ct), ψ(x+ct)), 即

    $$ \begin{array}{ll} \sup\limits _{x \in \mathbb{R}}|u(x, t)-\phi(x+c t)| \leqslant C \mathrm{e}^{-\mu t}, & t \geqslant 0, \\ \sup\limits _{x \in \mathbb{R}}|v(x, t)-\psi(x+c t)| \leqslant C \mathrm{e}^{-\mu t}, & t \geqslant 0, \end{array} $$

    其中常数μ>0.

    本节将通过加权能量的方法证明本文的主要结论.首先, 类似文献[22, 27]中的证明, 我们可以给出系统(1)的存在性定理和比较原理.

    引理2   假设a1a2-b1b2>0, 2r1-d1>0, 2r2-d2>0成立, 若初始条件满足

    $$ (0, 0) \leqslant\left(u_0(x, s), v_0(x, s)\right) \leqslant\left(k_1, k_2\right), \quad(x, s) \in \mathbb{R} \times[-\tau, 0], $$

    那么Cauchy问题(1)、(2)的解存在唯一且满足

    $$ (0, 0) \leqslant(u(x, t), v(x, t)) \leqslant\left(k_1, k_2\right), \quad(x, t) \in \mathbb{R} \times \mathbb{R}_{+}. $$

    令(u(x, t), v(x, t))和(u(x, t), v(x, t))分别是系统(1)带有初值(u0(x, s), v0(x, s))和(u0(x, s), v0(x, s))的解.如果

    $$ (0, 0) \leqslant\left(\underline{u}_0(x, s), \underline{v}_0(x, s)\right) \leqslant\left(\bar{u}_0(x, s), \bar{v}_0(x, s)\right) \leqslant\left(k_1, k_2\right), \quad(x, s) \in \mathbb{R} \times[-\tau, 0], $$

    那么

    $$ (0, 0) \leqslant(\underline{u}(x, t), \underline{v}(x, t)) \leqslant(\bar{u}(x, t), \bar{v}(x, t)) \leqslant\left(k_1, k_2\right), \quad(x, t) \in \mathbb{R} \times \mathbb{R}_{+} . $$

    若(u0(x, s), v0(x, s))满足

    $$ (0, 0) \leqslant\left(u_0(x, s), v_0(x, s)\right) \leqslant\left(k_1, k_2\right), \quad(x, s) \in \mathbb{R} \times[-\tau, 0], $$

    $$ \left\{ \begin{array}{*{35}{l}} u_{0}^{+}(x, s)=\max \left\{ {{u}_{0}}(x, s), \phi (x+cs) \right\}, \\ u_{0}^{-}(x, s)=\min \left\{ {{u}_{0}}(x, s), \phi (x+cs) \right\}, \\ v_{0}^{+}(x, s)=\max \left\{ {{v}_{0}}(x, s), \psi (x+cs) \right\}, \\ v_{0}^{-}(x, s)=\min \left\{ {{v}_{0}}(x, s), \psi (x+cs) \right\}, \\ \end{array} \right.\ \ \ \ \ \ \ \ (x, s)\in \mathbb{R}\times [-\tau , 0], $$ (5)

    那么

    $$ \left\{\begin{array}{l} 0 \leqslant u_0^{-}(x, s) \leqslant u_0(x, s) \leqslant u_0^{+}(x, s) \leqslant k_1, \\ 0 \leqslant u_0^{-}(x, s) \leqslant \phi(x+c s) \leqslant u_0^{+}(x, s) \leqslant k_1, \\ 0 \leqslant v_0^{-}(x, s) \leqslant v_0(x, s) \leqslant v_0^{+}(x, s) \leqslant k_2, \\ 0 \leqslant v_0^{-}(x, s) \leqslant \psi(x+c s) \leqslant v_0^{+}(x, s) \leqslant k_2, \end{array} \quad(x, s) \in \mathbb{R} \times[-\tau, 0] .\right. $$ (6)

    设(u-(x, t), v-(x, t))和(u+(x, t), v+(x, t))分别是系统(1)的初值为(u0-(x, s), v0-(x, s))和(u0+(x, s), v0+(x, s))相对应的解, 由引理2可得

    $$ \left\{ \begin{array}{*{35}{l}} 0\le {{u}^{-}}(x, t)\le u(x, t)\le {{u}^{+}}(x, t)\le {{k}_{1}}, \\ 0\le {{u}^{-}}(x, t)\le \phi (x+ct)\le {{u}^{+}}(x, t)\le {{k}_{1}}, \\ 0\le {{v}^{-}}(x, t)\le v(x, t)\le {{v}^{+}}(x, t)\le {{k}_{2}}, \\ 0\le {{v}^{-}}(x, t)\le \psi (x+ct)\le {{v}^{+}}(x, t)\le {{k}_{2}}, \\ \end{array} \right.\ \ \ \ \ \ \ \ (x, t)\in \mathbb{R}\times {{\mathbb{R}}_{+}}. $$ (7)

    接下来将分三步来证明主要的结果.

    首先, 证明u+(x, t)收敛到ϕ(x+ct). 令

    $$ \begin{aligned} & U(\xi, t)=u^{+}(x, t)-\phi(x+c t), U_0(\xi, s)=u_0^{+}(x, s)-\phi(x+c s), \\ & V(\xi, t)=v^{+}(x, t)-\psi(x+c t), V_0(\xi, s)=v_0^{+}(x, s)-\psi(x+c s), \end{aligned} $$

    其中ξ=x+ct.由式(5)和(6), 有

    $$ (0, 0) \leqslant(U(\xi, t), V(\xi, t)) \leqslant\left(k_1, k_2\right), (0, 0) \leqslant\left(U_0(\xi, s), V_0(\xi, s)\right) \leqslant\left(k_1, k_2\right) \text {. } $$

    通过计算易得U(ξ, t), V(ξ, t)满足

    $$ \left\{ \begin{align} & {{U}_{t}}+c{{U}_{\xi }}={{d}_{1}}\mathcal{D} [U(\xi , t)]-{{r}_{1}}U\left[ 2{{a}_{1}}\phi -1-{{b}_{1}}\psi (\xi -c\tau )-{{b}_{1}}V(\xi -c\tau , t-\tau ) \right]- \\ & \ \ \ \ \ \ \ {{a}_{1}}{{r}_{1}}{{U}^{2}}+{{b}_{1}}{{r}_{1}}\phi V(\xi -c\tau , t-\tau ), \\ & {{V}_{t}}+c{{V}_{\xi }}={{d}_{2}}\mathcal{D}[V(\xi , t)]-{{r}_{2}}V\left[ 2{{a}_{2}}\psi -1-{{b}_{2}}\phi (\xi -c\tau )-{{b}_{2}}U(\xi -c\tau , t-\tau ) \right]- \\ & \ \ \ \ \ \ \ {{a}_{2}}{{r}_{2}}{{V}^{2}}+{{b}_{2}}{{r}_{2}}\psi U(\xi -c\tau , t-\tau ), \\ \end{align} \right. $$ (8)

    式(8)的方程两端分别乘以e2μtω(ξ)U(ξ, t)和e2μtω(ξ)V(ξ, t), 其中μ>0待定, 可得

    $$ \begin{aligned} & \left(\frac{1}{2} \mathrm{e}^{2 \mu t} \omega U^2\right)_t+\left(\frac{c}{2} \mathrm{e}^{2 \mu t} \omega U^2\right)_{\xi}-\mathrm{e}^{2 \mu t} \omega U d_1[U(\xi+1, t)+U(\xi-1, t)]+ \\ & \quad \mathrm{e}^{2 \mu t} \omega U^2\left\{2 d_1-\mu-\frac{c}{2} \frac{\omega^{\prime}}{\omega}+r_1\left[2 a_1 \phi-1-b_1 \psi(\xi-c \tau)-b_1 V(\xi-c \tau, t-\tau)\right]\right\}= \\ & \quad-a_1 r_1 \mathrm{e}^{2 \mu t} \omega U^3+b_1 r_1 \phi \mathrm{e}^{2 \mu t} \omega U V(\xi-c \tau, t-\tau), \end{aligned} $$ (9)
    $$ \begin{aligned} & \left(\frac{1}{2} \mathrm{e}^{2 \mu t} \omega V^2\right)_t+\left(\frac{c}{2} \mathrm{e}^{2 \mu t} \omega V^2\right)_{\xi}-\mathrm{e}^{2 \mu t} \omega V d_2[V(\xi+1, t)+V(\xi-1, t)]+ \\ & \mathrm{e}^{2 \mu t} \omega V^2\left\{2 d_2-\mu-\frac{c}{2} \frac{\omega^{\prime}}{\omega}+r_2\left[2 a_2 \psi-1-b_2 \phi(\xi-c \tau)-b_2 U(\xi-c \tau, t-\tau)\right]\right\}= \\ & \quad-a_2 r_2 \mathrm{e}^{2 \mu t} \omega V^3+b_2 r_2 \psi \mathrm{e}^{2 \mu t} \omega V U(\xi-c \tau, t-\tau) . \end{aligned} $$ (10)

    注意到UHω1VHω1, 那么由文献[27]可得

    $$ \left.\left\{\mathrm{e}^{2 \mu t} \omega(\xi) U^2(\xi, t)\right\}\right|_{\xi=-\infty} ^{\infty}=0, \left.\left\{\mathrm{e}^{2 \mu t} \omega(\xi) V^2(\xi, t)\right\}\right|_{\xi=-\infty} ^{\infty}=0 . $$

    因此, 对式(9)和式(10)分别关于ξt在$ \mathbb{R} $×[0, t]上积分可得

    $$ \begin{aligned} & \mathrm{e}^{2 \mu t}\|U(t)\|_{L_\omega^2}^2-2 d_1 \int_0^t \int_{\mathrm{R}} \mathrm{e}^{2 \mu s} \omega U(\xi, s)[U(\xi+1, s)+U(\xi-1, s)] \mathrm{d} \xi \mathrm{d} s+ \\ & \quad \int_0^t \int_{\mathbb{R}} \mathrm{e}^{2 \mu s} \omega U^2\left\{4 d_1-2 \mu-c \frac{\omega^{\prime}}{\omega}+2 r_1\left[2 a_1 \phi-1-b_1 \psi(\xi-c \tau)-b_1 V(\xi-c \tau, s-\tau)\right]\right\} \mathrm{d} \xi \mathrm{d} s \leqslant \\ & \quad\left\|U_0(0)\right\|_{L_\omega^2}^2+2 b_1 r_1 \int_0^t \int_{\mathbb{R}} \phi \mathrm{e}^{2 \mu s} \omega U V(\xi-c \tau, s-\tau) \mathrm{d} \xi \mathrm{d} s, \end{aligned} $$ (11)
    $$ \begin{aligned} \mathrm{e}^{2 \mu t} & \|V(t)\|_{L_\omega^2}^2-2 d_2 \int_0^t \int_{\mathbb{R}} \mathrm{e}^{2 \mu s} \omega V(\xi, s)[V(\xi+1, s)+V(\xi-1, s)] \mathrm{d} \xi \mathrm{d} s+ \\ & \int_0^t \int_{\mathbb{R}} \mathrm{e}^{2 \mu s} \omega V^2\left\{4 d_2-2 \mu-c \frac{\omega^{\prime}}{\omega}+2 r_2\left[2 a_2 \psi-1-b_2 \phi(\xi-c \tau)-b_2 U(\xi-c \tau, s-\tau)\right]\right\} \mathrm{d} \xi \mathrm{d} s \leqslant \\ & \left\|V_0(0)\right\|_{L_\omega^2}^2+2 b_2 r_2 \int_0^t \int_{\mathbb{R}} \psi \mathrm{e}^{2 \mu s} \omega V U(\xi-c \tau, s-\tau) \mathrm{d} \xi \mathrm{d} s . \end{aligned} $$ (12)

    由Cauchy-Schwarz不等式可得

    $$ \begin{aligned} & 2 \int_0^t \int_{\mathbb{R}} \mathrm{e}^{2 \mu s} \omega W(\xi, s) W(\xi \pm 1, s) \mathrm{d} \xi \mathrm{d} s \\ & \quad \int_0^t \int_{\mathbb{R}} \mathrm{e}^{2 \mu s} \omega\left[W^2(\xi, s)+W^2(\xi \pm 1, s)\right] \mathrm{d} \xi \mathrm{d} s= \\ & \quad \int_0^t \int_{\mathbb{R}} \mathrm{e}^{2 \mu s} \omega W^2(\xi, s) \mathrm{d} \xi \mathrm{d} s+\int_0^t \int_{\mathbb{R}} \frac{\omega(\xi \mp 1)}{\omega(\xi)} \mathrm{e}^{2 \mu s} \omega W^2(\xi, s) \mathrm{d} \xi \mathrm{d} s, \end{aligned} $$ (13)

    其中W=U, V.类似地, 对于式(11)和式(12)右端的第二项, 有

    $$ \begin{aligned} & 2 b_1 r_1 \int_0^t \int_{\mathbb{R}} \mathrm{e}^{2 \mu s} \omega \phi U V(\xi-c \tau, s-\tau) \mathrm{d} \xi \mathrm{d} s \leqslant \\ & b_1 r_1 \int_0^t \int_{\mathbb{R}} \mathrm{e}^{2 \mu s} \omega \phi U^2(\xi, s) \mathrm{d} \xi \mathrm{d} s+b_1 r_1 \mathrm{e}^{2 \mu \tau} \int_0^t \int_{\mathbb{R}} \mathrm{e}^{2 \mu s} \omega(\xi+c \tau) \phi(\xi+c \tau) V^2(\xi, s) \mathrm{d} \xi \mathrm{d} s+ \\ & \\ & b_1 r_1 \mathrm{e}^{2 \mu \tau} \int_{-\tau}^0 \int_{\mathbb{R}} \mathrm{e}^{2 \mu s} \omega(\xi+c \tau) \phi(\xi+c \tau) V_0^2(\xi, s) \mathrm{d} \xi \mathrm{d} s, \end{aligned} $$ (14)
    $$ \begin{aligned} & 2 b_2 r_2 \int_0^t \int_{\mathbb{R}} \mathrm{e}^{2 \mu s} \omega \psi V U(\xi-c \tau, s-\tau) \mathrm{d} \xi \mathrm{d} s \leqslant \\ & b_2 r_2 \int_0^t \int_{\mathbb{R}} \mathrm{e}^{2 \mu s} \omega \psi V^2(\xi, s) \mathrm{d} \xi \mathrm{d} s+b_2 r_2 \mathrm{e}^{2 \mu \tau} \int_0^t \int_{\mathbb{R}} \mathrm{e}^{2 \mu s} \omega(\xi+c \tau) \psi(\xi+c \tau) U^2(\xi, s) \mathrm{d} \xi \mathrm{d} s+ \\ & b_2 r_2 \mathrm{e}^{2 \mu \tau} \int_{-\tau}^0 \int_{\mathbb{R}} \mathrm{e}^{2 \mu s} \omega(\xi+c \tau) \psi(\xi+c \tau) U_0^2(\xi, s) \mathrm{d} \xi \mathrm{d} s. \end{aligned} $$ (15)

    将式(13)—(15)代入到式(11)和(12)中, 由于对于ξ∈$ \mathbb{R} $有ω(ξ+)/ω(ξ)≤1, 因此

    $$ \begin{aligned} & \mathrm{e}^{2 \mu t}\|U(t)\|_{L_\omega^2}^2+ \\ & \int_0^t \int_{\mathbb{R}} \mathrm{e}^{2 \mu s} \omega U^2\left\{4 d_1-2 \mu-c \frac{\omega^{\prime}}{\omega}-b_1 r_1 \phi-d_1\left[2+\frac{\omega(\xi-1)}{\omega(\xi)}+\frac{\omega(\xi+1)}{\omega(\xi)}\right]\right\} \mathrm{d} \xi \mathrm{d} s+ \\ & \int_0^t \int_{\mathbb{R}} 2 \mathrm{e}^{2 \mu s} \omega U^2 r_1\left[2 a_1 \phi-1-b_1 \psi(\xi-c \tau)-b_1 V(\xi-c \tau, s-\tau)\right] \mathrm{d} \xi \mathrm{d} s \leqslant \\ & \left\|U_0(0)\right\|_{L_\omega^2}^2+b_1 r_1 k_1 \mathrm{e}^{2 \mu \tau} \int_{-\tau}^0\left\|V_0(s)\right\|_{L_\omega^2}^2 \mathrm{~d} s+ \\ & b_1 r_1 \mathrm{e}^{2 \mu \tau} \int_0^t \int_{\mathbb{R}} \mathrm{e}^{2 \mu s} \omega(\xi+c \tau) \phi(\xi+c \tau) V^2(\xi, s) \mathrm{d} \xi \mathrm{d} s, \end{aligned} $$ (16)
    $$ \begin{aligned} & \mathrm{e}^{2 \mu t}\|V(t)\|_{L_\omega^2}^2+ \\ & \int_0^t \int_{\mathbb{R}} \mathrm{e}^{2 \mu s} \omega V^2\left\{4 d_2-2 \mu-c \frac{\omega^{\prime}}{\omega}-b_2 r_2 \psi-d_2\left[2+\frac{\omega(\xi-1)}{\omega(\xi)}+\frac{\omega(\xi+1)}{\omega(\xi)}\right]\right\} \mathrm{d} \xi \mathrm{d} s+ \\ & \int_0^t \int_{\mathbb{R}} 2 \mathrm{e}^{2 \mu s} \omega V^2 r_2\left[2 a_2 \psi-1-b_2 \phi(\xi-c \tau)-b_2 U(\xi-c \tau, s-\tau)\right] \mathrm{d} \xi \mathrm{d} s \leqslant \\ & \left\|V_0(0)\right\|_{L_\omega^2}^2+b_2 r_2 k_2 \mathrm{e}^{2 \mu \tau} \int_{-\tau}^0\left\|U_0(s)\right\|_{L_\omega^2}^2 \mathrm{~d} s+ \\ & b_2 r_2 \mathrm{e}^{2 \mu \tau} \int_0^t \int_{\mathbb{R}} \mathrm{e}^{2 \mu s} \omega(\xi+c \tau) \psi(\xi+c \tau) U^2(\xi, s) \mathrm{d} \xi \mathrm{d} s . \end{aligned} $$ (17)

    结合式(16)与式(17), 可得

    $$ \begin{aligned} & \mathrm{e}^{2 \mu t}\left(\|U(t)\|_{L_\omega^2}^2+\|V(t)\|_{L_\omega^2}^2\right)+ \\ & \quad \int_0^t \int_{\mathbb{R}} \mathrm{e}^{2 \mu s} \omega(\xi)\left(B_{\mu, \omega}^1(\xi, s) U^2(\xi, s)+B_{\mu, \omega}^2(\xi, s) V^2(\xi, s)\right) \mathrm{d} \xi \mathrm{d} s \leqslant \\& \quad C_0\left(\left\|U_0(0)\right\|_{L_\omega^2}^2+\left\|V_0(0)\right\|_{L_\omega^2}^2+\int_{-\tau}^0\left\|U_0(s)\right\|_{L_\omega^2}^2 \mathrm{~d} s+\int_{-\tau}^0\left\|V_0(s)\right\|_{L_\omega^2}^2 \mathrm{~d} s\right), \end{aligned} $$ (18)

    其中

    $$ \begin{aligned} & B_{\mu, \omega}^1(\xi, t)=A_\omega^1(\xi, t)-2 \mu-b_2 r_2\left(\mathrm{e}^{2 \mu \tau}-1\right) \psi(\xi+c \tau) \frac{\omega(\xi+c \tau)}{\omega(\xi)}, \\ & B_{\mu, \omega}^2(\xi, t)=A_\omega^2(\xi, t)-2 \mu-b_1 r_1\left(\mathrm{e}^{2 \mu \tau}-1\right) \phi(\xi+c \tau) \frac{\omega(\xi+c \tau)}{\omega(\xi)}, \\ & A_\omega^1(\xi, t)=4 d_1-c \frac{\omega^{\prime}}{\omega}+2 r_1\left[2 a_1 \phi-1-b_1 \psi(\xi-c \tau)-b_1 V(\xi-c \tau, t-\tau)\right]-b_1 r_1 \phi- \\ & b_2 r_2 \psi(\xi+c \tau) \frac{\omega(\xi+c \tau)}{\omega(\xi)}-d_1\left[2+\frac{\omega(\xi-1)}{\omega(\xi)}+\frac{\omega(\xi+1)}{\omega(\xi)}\right], \\ & A_\omega^2(\xi, t)=4 d_2-c \frac{\omega^{\prime}}{\omega}+2 r_2\left[2 a_2 \psi-1-b_2 \phi(\xi-c \tau)-b_2 U(\xi-c \tau, t-\tau)\right]-b_2 r_2 \psi- \\ & b_1 r_1 \phi(\xi+c \tau) \frac{\omega(\xi+c \tau)}{\omega(\xi)}-d_2\left[2+\frac{\omega(\xi-1)}{\omega(\xi)}+\frac{\omega(\xi+1)}{\omega(\xi)}\right] \cdot \end{aligned} $$

    接下来, 我们通过下面的引理证明Aωi(ξ, t)(i=1, 2)大于某些正整数.

    引理3   假设(A1)与(A2)成立.当c>$ \tilde{c}$及ξ∈$ \mathbb{R} $, t≥0时, 存在正整数Ci, i=1, 2, 使得

    $$ A_\omega^i(\xi, t) \geqslant C_i, \quad i=1, 2 . $$

    证明   首先, 我们证明存在常数C1>0使得Aω1(ξ, t)≥C1.当ξ < ξ0-1时, 此时ω(ξ+)/ω(ξ)≤1.由前面的讨论及直接计算可得

    $$ A_\omega^1(\xi, t) \geqslant c \beta_0-2 r_1\left(1+2 b_1 k_2\right)-b_1 r_1 k_1-b_2 r_2 k_2-d_1\left(\mathrm{e}^{\beta_0}+\mathrm{e}^{-\beta_0}-1\right)+d_1 \cdot $$

    类似地, 当ξ0-1≤ξξ0时, 此时ω(ξ+)/ω(ξ)≤1.然后可得

    $$ A_\omega^1(\xi, t) \geqslant c \beta_0-2 r_1\left(1+2 b_1 k_2\right)-b_1 r_1 k_1-b_2 r_2 k_2-d_1\left(\mathrm{e}^{\beta_0}+\mathrm{e}^{-\beta_0}-1\right)+d_1 \mathrm{e}^{-\beta_0}. $$

    ξ0 < ξξ0+1时, 此时ω(ξ+)/ω(ξ)=1.然后可得

    $$ A_\omega^1(\xi, t) \geqslant 2 r_1\left[2 a_1 \phi\left(\xi_0\right)-1-2 b_1 k_2\right]-b_1 r_1 k_1-b_2 r_2 k_2-d_1\left(\mathrm{e}^{\beta_0}+\mathrm{e}^{-\beta_0}-1\right)+d_1 \mathrm{e}^{-\beta_0}. $$

    ξ>ξ0+1时, 此时ω(ξ+)/ω(ξ)=1.然后可得

    $$ A_\omega^1(\xi, t) \geqslant 2 r_1\left[2 a_1 \phi\left(\xi_0\right)-1-2 b_1 k_2\right]-b_1 r_1 k_1-b_2 r_2 k_2-d_1\left(\mathrm{e}^{\beta_0}+\mathrm{e}^{-\beta_0}-1\right)+d_1 \cdot $$

    综上所述, 令C1=d1e-β0, 可得Aω1(ξ, t)≥C1.

    类似地, 可证存在常数C2>0使得Aω2(ξ, t)≥C2.

    在引理3基础之上, 我们将进一步给出Bμ, ωi(ξ, t)(i=1, 2)的正下界.

    引理4   设μ1>0是如下方程

    $$ \bar{C}-2 \mu-h\left(\mathrm{e}^{2 \mu \tau}-1\right)=0 $$

    的唯一实根, 其中C=min{C1, C2}和h=max{b1r1k1, b2r2k2}. 取0 < μ < μ1, 那么存在C3>0使得

    $$ B_{\mu, \omega}^i(\xi, t) \geqslant C_3, \quad i=1, 2, \xi \in \mathbb{R}, t \geqslant 0. $$ (19)

    证明   由于(0, 0)≤(ϕ(ξ), ψ(ξ))≤(k1, k2)及ω(ξ+)/ω(ξ)≤1(ξ∈$ \mathbb{R} $), 因此, 当0 < μ < μ1时, 可得

    $$ \begin{aligned} & B_{\mu, \omega}^1(\xi, t)=A_\omega^1(\xi, t)-2 \mu-b_2 r_2\left(\mathrm{e}^{2 \mu \tau}-1\right) \psi(\xi+c \tau) \frac{\omega(\xi+c \tau)}{\omega(\xi)} \geqslant \\ & \quad \bar{C}-2 \mu-h\left(\mathrm{e}^{2 \mu \tau}-1\right):=C_3>0 . \end{aligned} $$

    类似地,可证Bμ, ω2(ξ, t)≥C3.

    然后, 将式(19)代入式(18), 我们可得如下估计.

    引理5   假设(A1)和(A2)成立.对于c>$\tilde{c} $, 存在正常数C4使得

    $$ \left\{\begin{array}{l} \mathrm{e}^{2 \mu t}\|U(t)\|_{L_\omega^2}^2 \leqslant C_4\left[\left\|U_0(0)\right\|_{L_\omega^2}^2+\left\|V_0(0)\right\|_{L_\omega^2}^2+\int_{-\tau}^0\left(\left\|U_0(s)\right\|_{L_\omega^2}^2+\left\|V_0(s)\right\|_{L_\omega^2}^2\right) \mathrm{d} s\right], \\ \mathrm{e}^{2 \mu t}\|V(t)\|_{L_\omega^2}^2 \leqslant C_4\left[\left\|U_0(0)\right\|_{L_\omega^2}^2+\left\|V_0(0)\right\|_{L_\omega^2}^2+\int_{-\tau}^0\left(\left\|U_0(s)\right\|_{L_\omega^2}^2+\left\|V_0(s)\right\|_{L_\omega^2}^2\right) \mathrm{d} s\right] . \end{array}\right. $$ (20)

    之后, 我们给出UξVξ在加权空间中的L2估计.类似地, 对式(8)方程两端关于ξ求导并对所得方程的两边分别乘以e2μtω(ξ)Uξ(ξ, t)和e2μtω(ξ)Vξ(ξ, t), 可得

    $$ \begin{aligned} & \left(\frac{1}{2} \mathrm{e}^{2 \mu t} \omega U_{\xi}^2\right)_t+\left(\frac{c}{2} \mathrm{e}^{2 \mu t} \omega U_{\xi}^2\right)_{\xi}-\mathrm{e}^{2 \mu t} \omega U_{\xi} d_1\left[U_{\xi}(\xi+1, t)+U_{\xi}(\xi-1, t)\right]+ \\ & \mathrm{e}^{2 \mu t} \omega U_{\xi}^2\left\{2 d_1-\mu-\frac{c}{2} \frac{\omega^{\prime}}{\omega}+r_1\left[2 a_1 \phi-1-b_1 \psi(\xi-c \tau)-b_1 V(\xi-c \tau, t-\tau)\right]\right\}+ \\ & r_1 \mathrm{e}^{2 \mu t} \omega U U_{\xi}\left[2 a_1 \phi^{\prime}-b_1 \psi^{\prime}(\xi-c \tau)-b_1 V_{\xi}(\xi-c \tau, t-\tau)\right]= \\ & \quad-2 a_1 r_1 \mathrm{e}^{2 \mu t} \omega U U_{\xi}^2+b_1 r_1 \phi^{\prime} \mathrm{e}^{2 \mu t} \omega U_{\xi} V(\xi-c \tau, t-\tau)+b_1 r_1 \phi \mathrm{e}^{2 \mu t} \omega U_{\xi} V_{\xi}(\xi-c \tau, t-\tau), \\ & \left(\frac{1}{2} \mathrm{e}^{2 \mu t} \omega V_{\xi}^2\right)_t+\left(\frac{c}{2} \mathrm{e}^{2 \mu t} \omega V_{\xi}^2\right)_{\xi}-\mathrm{e}^{2 \mu t} \omega V_{\xi} d_2\left[V_{\xi}(\xi+1, t)+V_{\xi}(\xi-1, t)\right]+ \\ & \mathrm{e}^{2 \mu t} \omega V_{\xi}^2\left\{2 d_2-\mu-\frac{c}{2} \frac{\omega^{\prime}}{\omega}+r_2\left[2 a_2 \psi-1-b_2 \phi(\xi-c \tau)-b_2 U(\xi-c \tau, t-\tau)\right]\right\}+ \\ & r_2 \mathrm{e}^{2 \mu t} \omega V V_{\xi}\left[2 a_2 \psi^{\prime}-b_2 \phi^{\prime}(\xi-c \tau)-b_2 U_{\xi}(\xi-c \tau, t-\tau)\right]= \\ & -2 a_2 r_2 \mathrm{e}^{2 \mu t} \omega V V_{\xi}^2+b_2 r_2 \psi^{\prime} \mathrm{e}^{2 \mu t} \omega V_{\xi} U(\xi-c \tau, t-\tau)+b_2 r_2 \psi \mathrm{e}^{2 \mu t} \omega V_{\xi} U_{\xi}(\xi-c \tau, t-\tau) . \end{aligned} $$

    重复上述证明过程可得引理6.

    引理6  假设(A1)与(A2)成立.对于c>$ \tilde{c}$, 存在正常数C5使得

    $$ \left\{\begin{array}{l} \mathrm{e}^{2 \mu t}\left\|U_{\xi}(t)\right\|_{L_\omega^2}^2 \leqslant C_5\left[\left\|U_0(0)\right\|_{H_\omega^1}^2+\left\|V_0(0)\right\|_{H_\omega^1}^2+\int_{-\tau}^0\left(\left\|U_0(s)\right\|_{H_\omega^1}^2+\left\|V_0(s)\right\|_{H_\omega^1}^2\right) \mathrm{d} s\right], \\ \mathrm{e}^{2 \mu t}\left\|V_{\xi}(t)\right\|_{L_\omega^2}^2 \leqslant C_5\left[\left\|U_0(0)\right\|_{H_\omega^1}^2+\left\|V_0(0)\right\|_{H_\omega^1}^2+\int_{-\tau}^0\left(\left\|U_0(s)\right\|_{H_\omega^1}^2+\left\|V_0(s)\right\|_{H_\omega^1}^2\right) \mathrm{d} s\right] . \end{array}\right. $$

    由引理5和6, 可得到如下估计.

    引理7   假设(A1)与(A2)成立.对于c>$\tilde{c} $, 有

    $$ \left\{\begin{array}{l} \|U(t)\|_{H_\omega^1} \leqslant C \mathrm{e}^{-\mu t}\left[\left\|U_0(0)\right\|_{H_\omega^1}^2+\left\|V_0(0)\right\|_{H_\omega^1}^2+\int_{-\tau}^0\left(\left\|U_0(s)\right\|_{H_\omega^1}^2+\left\|V_0(s)\right\|_{H_\omega^1}^2\right) \mathrm{d} s\right]^{1 / 2}, \\ \|V(t)\|_{H_\omega^1} \leqslant C \mathrm{e}^{-\mu t}\left[\left\|U_0(0)\right\|_{H_\omega^1}^2+\left\|V_0(0)\right\|_{H_\omega^1}^2+\int_{-\tau}^0\left(\left\|U_0(s)\right\|_{H_\omega^1}^2+\left\|V_0(s)\right\|_{H_\omega^1}^2\right) \mathrm{d} s\right]^{1 / 2} . \end{array}\right. $$

    由于ω(ξ)≥1, 运用Sobolev嵌入定理H1($ \mathbb{R} $)$ \hookrightarrow$C($ \mathbb{R} $), 我们得到如下衰减结论.

    引理8   假设(A1)与(A2)成立.对于c>$\tilde{c} $, 有

    $$ \begin{cases}\sup\limits _{x \in \mathbb{R}} \mid & u^{+}(x, t)-\phi(x+c t) \mid= \\ & \sup\limits _{\xi \in \mathbb{R}}|U(\xi, t)| \leqslant\|U(t)\|_{H^1} \leqslant\|U(t)\|_{H_\omega^1} \leqslant C \mathrm{e}^{-\mu t}, \quad t \geqslant 0, \\ \sup\limits _{x \in \mathbb{R}} \mid & v^{+}(x, t)-\psi(x+c t) \mid= \\ & \sup\limits _{\xi \in \mathbb{R}}|V(\xi, t)| \leqslant\|V(t)\|_{H^1} \leqslant\|V(t)\|_{H_\omega^1} \leqslant C \mathrm{e}^{-\mu t}, \quad t \geqslant 0 .\end{cases} $$

    其次, 重复上述证明过程可得引理9.

    引理9   假设(A1)与(A2)成立.对于c>$\tilde{c} $, 有

    $$ \begin{cases}\sup\limits _{x \in \mathbb{R}}\left|u^{-}(x, t)-\phi(x+c t)\right| \leqslant C \mathrm{e}^{-\mu t}, & t \geqslant 0, \\ \sup\limits _{x \in \mathbb{R}}\left|v^{-}(x, t)-\psi(x+c t)\right| \leqslant C \mathrm{e}^{-\mu t}, & t \geqslant 0 .\end{cases} $$

    最后, 运用挤压定理并结合引理8与9, 可证明定理1成立.

    本文运用L2-加权能量方法、比较原理和挤压技术得到了具有时滞的离散Lotka-Volterra合作系统波前解的全局指数稳定性,并解决了离散扩散算子及时滞共同作用下建立能量估计的问题.此稳定性结果可以有助于理解解(u, v)渐近地表现为以正速度c传播的行波解,并且随着时间t的推移,这两个物种趋于合作水平.

  • [1] 叶其孝, 李正元, 王明新, 等. 反应扩散方程引论[M]. 2版. 北京: 科学出版社, 2011.

    YE Qixiao, LI Zhengyuan, WANG Mingxin, et al. Introduction to Reaction-Diffusion Equations[M]. 2nd ed. Beijing: Science Press, 2011. (in Chinese))
    [2] HUANG J H, ZOU X F. Travelling wave fronts in diffusive and cooperative Lotka-Volterra system with delays[J]. Journal of Mathematical Analysis and Applications, 2002, 271: 455-466. doi: 10.1016/S0022-247X(02)00135-X
    [3] WU J H, ZOU X F. Traveling wave fronts of reaction diffusion systems with delay[J]. Journal of Dynamics and Differential Equations, 2001, 13: 651-687. doi: 10.1023/A:1016690424892
    [4] LV G Y, WANG M X. Traveling wave front in diffusive and competitive Lotka-Volterra system with delays[J]. Nonlinear Analysis: Real World Applications, 2010, 11(3): 1323-1329. doi: 10.1016/j.nonrwa.2009.02.020
    [5] GUO J S, WU C H. Traveling wave front for a two-component lattice dynamical system arising in competitionmodels[J]. Journal of Differential Equations, 2012, 252(8): 4357-4391. doi: 10.1016/j.jde.2012.01.009
    [6] LI K, HUANG J H, LI X, et al. Traveling wave fronts in a delayed lattice competitive system[J]. Applicable Analysis, 2017, 97(6): 982-999.
    [7] 张秋, 陈广生. 一类具有非线性发生率与时滞的非局部扩散SIR模型的临界波的存在性[J]. 应用数学和力学, 2019, 40(7): 713-727. doi: 10.21656/1000-0887.390208

    ZHANG Qiu, CHEN Guangsheng. Existence of critical traveling waves for nonlocal dispersal SIR models with delay and nonlinear incidence[J]. Applied Mathematics and Mechanics, 2019, 40(7): 713-727. (in Chinese)) doi: 10.21656/1000-0887.390208
    [8] 谷雨萌, 黄明迪. 一类时间周期的时滞竞争系统行波解的存在性[J]. 应用数学和力学, 2020, 41(6): 658-668. doi: 10.21656/1000-0887.400275

    GU Yumeng, HUANG Mingdi. Existence of periodic traveling waves for time-periodic Lotka-Volterra competition systems with delay[J]. Applied Mathematics and Mechanics, 2020, 41(6): 658-668. (in Chinese)) doi: 10.21656/1000-0887.400275
    [9] BRAMSON M. Convergence of Solutions of the Kolmogorov Equation to Traveling Waves[M]. Memoirs of the American Mathematical Society, 1983.
    [10] KIRCHGÄSSNER K. On the nonlinear dynamics of travelling fronts[J]. Journal of Differential Equations, 1992, 96(2): 256-278. doi: 10.1016/0022-0396(92)90153-E
    [11] TSAI J C, SNEYD J. Existence and stability of traveling waves in bufferedsystems[J]. SIAM Journal on Mathematical Analysis, 2005, 66: 237-265.
    [12] MA S W, ZHAO X Q. Global asymptotic stability of minimal fronts in monostable lattice equations[J]. Discrete and Continuous Dynamical Systems, 2008, 21: 259-275. doi: 10.3934/dcds.2008.21.259
    [13] WANG Y, CAO X Y, MA Z H, et al. Global stability of noncritical traveling front solutions of Fisher-type equations with degenerate nonlinearity[J]. Journal of Mathematical Physics, 2021, 62(5): 051506. doi: 10.1063/5.0043893
    [14] ZHOU Y H, YAN Z M, JI S G. Global stability of traveling waves with non-monotonicity for population dynamicsmodel[J]. Tokyo Journal of Mathematics, 2021, 44: 383-396.
    [15] LIU C C, MEI M, YANG J Q. Global stability of traveling waves for nonlocal time-delayed degenerate diffusion equation[J]. Journal of Differential Equations, 2022, 306: 60-100. doi: 10.1016/j.jde.2021.10.027
    [16] SCHAAF K W. Asymptotic behavior and traveling wave solutions for parabolic functional differentialequations[J]. Transactions of the American Mathematical Society, 1987, 302: 587-615.
    [17] MEI M, SO J W H, LI M Y, et al. Asymptotic stability of traveling waves for the Nicholson's blowflies equation with diffusion[J]. Proceedings of the Royal Society of Edinburgh Section A: Mathematics, 2004, 134(3): 579-594. doi: 10.1017/S0308210500003358
    [18] LIN C K, MEI M. On travellingwavefronts of the Nicholson's blowflies equations with diffusion[J]. Proceedings of the Royal Society of Edinburgh Section A: Mathematics, 2010, 140(1): 135-152. doi: 10.1017/S0308210508000784
    [19] MEI M, LIN C K, LIN C T, et al. Traveling wavefronts for time-delayed reaction-diffusion equation, Ⅰ: local nonlinearity[J]. Journal of Differential Equations, 2009, 247(2): 495-510. doi: 10.1016/j.jde.2008.12.026
    [20] MEI M, LIN C K, LIN C T, et al. Traveling wavefronts for time-delayed reaction-diffusion equation, Ⅱ: nonlocal nonlinearity[J]. Journal of Differential Equations, 2009, 247(2): 511-529. doi: 10.1016/j.jde.2008.12.020
    [21] YU Z X, XU F, ZHANG W G. Stability of invasion traveling waves for a competition system with nonlocaldispersals[J]. Applicable Analysis, 2017, 96(7): 1107-1125. doi: 10.1080/00036811.2016.1178242
    [22] ZHANG G B, DONG F D, LI W T. Uniqueness and stability of traveling waves for a three-species competition system with nonlocaldispersal[J]. Discrete and Continuous Dynamical Systems(Series B), 2019, 24(4): 1511-1541. doi: 10.3934/dcdsb.2018218
    [23] MA S, ZOU X F. Existence, uniqueness and stability of traveling waves in a discrete reaction-diffusionmonostable equation with delay[J]. Journal of Differential Equations, 2005, 217(1): 54-87. doi: 10.1016/j.jde.2005.05.004
    [24] GUO S J, ZIMMER J. Stability of travelingwavefronts in discrete reaction-diffusion equations with nonlocal delay effects[J]. Nonlinearity, 2015, 28(2): 463-492. doi: 10.1088/0951-7715/28/2/463
    [25] HSU C H, LIN J J, YANG T S. Stability for monostable wave fronts of delayed lattice differential equations[J]. Journal of Dynamics and Differential Equations, 2017, 29: 323-342. doi: 10.1007/s10884-015-9447-9
    [26] YU Z X, HSU C H. Wave propagation and its stability for a class of discrete diffusion systems[J]. Zeitschrift für Angewandte Mathematik und Physik, 2020, 71: 194. doi: 10.1007/s00033-020-01423-4
    [27] CHEN G S, WU S L, HSU C H. Stability of travelingwavefronts for a discrete diffusive competition system with three species[J]. Journal of Mathematical Analysis and Applications, 2019, 474(2): 909-930. doi: 10.1016/j.jmaa.2019.01.079
    [28] SU T, ZHANG G B. Stability of travelingwavefronts for a three-component Lotka-Volterra competition system on a lattice[J]. Electronic Journal of Differential Equations, 2018, 57(2018): 1-16.
    [29] HSU C H, LIN J J. Stability analysis of traveling wave solutions for lattice reaction-diffusion equations[J]. Discrete and Continuous Dynamical Systems(Series B), 2020, 25(5): 1757-1774. doi: 10.3934/dcdsb.2020001
    [30] 朱福国. 格上时滞Lotka-Volterra合作系统的波前解[J]. 生物数学学报, 2012, 27(1): 150-156. https://www.cnki.com.cn/Article/CJFDTOTAL-SWSX201201018.htm

    ZHU Fuguo. Traveling wavefrons of delayed Lotka-Volterra system on lattice[J]. Journal of Biomathematics, 2012, 27(1): 150-156. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SWSX201201018.htm
  • 加载中
计量
  • 文章访问数:  521
  • HTML全文浏览量:  193
  • PDF下载量:  57
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-05-23
  • 修回日期:  2022-06-28
  • 刊出日期:  2023-04-01

目录

/

返回文章
返回