非局部时滞反应扩散方程波前解的指数稳定性*

李盼晓

(西安电子科技大学 数学与统计学院，西安 710071)

摘要：该文考虑了一类具有非局部时滞反应项的空间非局部扩散模型, 主要研究其波前解的渐近稳定性及收敛率．通过构造加权函数，建立了相关线性方程的比较原理，证明了当初始扰动在加权最大范数意义下一致有界，满足初值问题的解将依时间指数收敛到波前解，而且得到其指数收敛率．

关键词：行波解; 非局部时滞; 稳定性

引 言

在自然界中，许多现象都可以通过相应的反应扩散模型进行描述．近些年来，空间非局部微分方程广泛应用于生物学和传染病学等多个领域，并引起了许多学者的极大关注．目前，对于此类方程的研究已经获得了一些成果[1-6]．在种群生物学中，由于种群的随机移动使得个体能够从各个可能的位置来到当前所处的位置，考虑到这一因素,Britton[7-8]首先提出用带有非局部时滞反应项的模型来描述这种现象．后来，Wang等[6]为了研究种群动力学，建立了如下非局部时滞反应扩散方程:

x∈R,t>0，

(1)

其中u(x,t)表示x位置、t时刻的种群密度,d>0是扩散项系数，时滞τ>0，k(·)是非负核函数．在单稳的假设下，Wang等[6]得到了模型(1)波前解的存在性，唯一性及渐近稳定性．在生物学领域，Laplace算子通常用于描述种群个体在栖息区域内的局部迁移和扩散．实际上，许多种群个体在空间区域上却是广泛迁移和大范围扩散的．为了准确描述这一非局部扩散现象，一些研究者引入了积分形式的扩散项[9-13]．特别地，Yu和Yuan[14]将非局部时滞的反应项引进非局部扩散模型:

x∈R,t>0，

(2)

其中d>0,τ>0是常数,

(J*u)(x,t)

最近,Guo和Zimmer[15]将以上各个因素考虑在内，建立了如下非局部时滞反应扩散模型，并得到其波前解的存在性、渐近行为、唯一性及单调性结果：

ut(x,t)=puxx(x,t)+d(J*u-u)(x,t)+f(u(x,t),(h**u)(x,t)),

x∈R,

(3)

其中p≥0,d≥0，且J(·),h(·)为非负核函数,

(h**u)(x,t)

该模型既含有局部Laplace扩散项又含有非局部扩散项，主要描述种群的扩散不仅依赖于空间位置x处、t时刻的种群密度u(x,t)，而且还依赖于整个空间中其他位置上的u值．反应项f(u,h**u)通常用于描述种群个体的出生率和死亡率对种群密度的影响,h**u表示种群密度在过去时间和空间上的加权平均．

当h(x,t)=h(x)δ(t-τ)，其中δ(·)为Dirac δ函数，方程(3)变为以下非局部时滞反应扩散方程：

ut(x,t)=puxx(x,t)+d(J*u-u)(x,t)+

f(u(x,t),(h*u)(x,t-τ)),x∈R,d≥0,p≥0,

(4)

其中

(h*u)(x,t)

模型(4)波前解的存在性、渐近行为和唯一性结果已经得到很好的研究[15]．然而在许多应用领域中，需要知道反应扩散方程的行波解在时间趋于无穷时的性态，即发展方程解的稳定性理论．因此，本文主要研究非局部时滞反应扩散方程(4)在单稳条件下，非临界波前解的指数稳定性．

首先假设非线性函数f(u,v)和h(u)满足以下条件:

(P1)f∈C([0,K]2,R)，f(0,0)=f(K,K)=0,f(u,u)>0，∀u∈(0,K)，对任意的(u,v)∈[0,K]2,∂2f(u,v)≥0，其中K是正常数．

(P2) 存在M>0和σ∈(0,1]，使得0≤∂1f(0,0)u+∂2f(0,0)v-f(u,v)≤M(u+v)1+σ,∀(u,v)∈[0,K]2,且∂1f(K,K)+∂2(K,K)<0．

(P3)J(·)，h(·)都是非负可积函数，且满足

(P4) 存在某个λ0>0(可能等于无穷)使得对任意的λ∈[0,λ0),有

J(x)eλxdx<+∞,h(x)eλxdx<+∞．

从条件(P1)可以看出, 方程(4)有两个平衡点分别为0和K．因此，结合条件(P1)和(P2)，可得

于是u=0是不稳定的,u=K是稳定的．本文不需要∂2f(0,0)>0．

反应扩散模型的稳定性理论在生物学[16-20]和传染病学[21-23]等多个领域都具有重要的作用．到目前为止,反应扩散方程行波解的稳定性已经得到了广泛研究[6，16-20，23-26]，但是对于非局部时滞反应扩散方程行波解稳定性的结果相对较少．Ma和Zhao[19]用挤压技术证明了单稳格微分方程行波解渐近稳定性，并且可以得到指数衰减率．Chen和Guo[24]通过对单稳拟线性方程建立挤压技术，得到波前解的渐近稳定性，但是却不能得到具体的收敛率．另外，Cheng和Yuan[16]通过构造合适的上下解结合比较原理和挤压技术，证明了具有非局部时滞反应扩散模型(2)的非临界波前解的渐近稳定性，然而该模型临界波前解的稳定性尚未得到．在生物学和传染病学中，临界行波解的稳定性对于研究生物入侵和疾病传播等问题非常重要．Guo和Zimmer[17]通过加权能量法和Green函数技巧，证明方程(1)相应空间离散模型波前解的全局指数稳定性，证明当波前解的初始扰动在负无穷远处衰减到零时，初值问题的解在c>c*条件下关于时间指数衰减收敛到波速为c的行波解，在c=c*条件下关于时间代数衰减收敛到波速为c*的行波解．特别地，时滞对行波解的稳定性没有影响．此外，Mei等[18]考虑具有年龄结构和空间扩散的单一物种的动力学行为，研究了一类具有年龄结构的非局部时滞反应扩散模型行波解的稳定性，主要通过加权能量法结合比较原理，证明大初始扰动下行波解是全局渐近稳定性的．

值得说明的是，加权能量法结合比较原理可以证明大波速和小波速行波解的渐近稳定性结果，而且能得到具体的收敛速度，但是需要先对解进行加权L1,L2及L∞能量估计，这使得计算过程变的复杂．为了进一步将稳定性的结果精确化,Ouyang等[20]和Wang等[25]利用权函数结合比较原理，分别证明了空间周期非局部时滞单种群模型和周期部分退化反应扩散系统行波解的指数渐近稳定性．随后,Wu和Chen[23]也用这种方法证明了非局部扩散传染病模型行波解的指数渐近稳定性．由此可知，权函数结合比较原理方法是证明行波解全局指数渐近稳定性的一个有效途径，相对于加权能量法和挤压技术，它在很大程度上简化了计算过程且能得到初值问题的解与行波解之间精确的衰减估计．那么此方法能否应用于模型(4)? 本文将对此问题给出一个肯定的答案．具体地说，本文通过特征值分析法，建立方程(4)相关线性系统的比较原理(后文中引理3)，得到当初始函数u0在加权范数的意义下趋于某一行波解φ时，满足此初值的解将指数收敛到它的行波解．

文章主体结构大体分为5个部分: 首先，给出一些基本知识和主要结论; 其次，通过主特征值分析的方法得出方程(4)线性初值问题的比较原理; 然后，通过权函数结合比较原理的方法证明方程(4)非临界波前解的指数渐近稳定性，并且得到指数收敛率; 第4节主要给出具体实例去验证所得的稳定性结果; 最后对全文工作进行总结．

1 主 要 结 论

本文中将C>0记为一般常数,Ci(i=1,2，…)为特殊常数，令I是一个区间，特别地，取I=R．设T>0是一个实数且B为Banach空间．记C([0,T],B)为[0,T]上的连续函数空间，L2([0,T],B)表示[0,T]上的L2函数空间．那么[0,+∞)上相应的函数空间可以用类似的方法定义．

假设u(x,t)=φ(x+ct)为方程(4)连接平衡点0和K的波前解，其中c为波速,φ为波廓．将u=φ(ξ),ξ=x+ct代入方程(4)得到行波方程：

(5)

其中

定义

显然，B(ξ)和Gj(ξ),j=1,2是非增的且满足

G1(+∞)=∂1f(K,K),B(+∞)=G2(+∞)=∂2f(K,K)．

对于λ∈C且Reλ<λ0，定义函数

显然,G(λ)关于λ∈[0,λ0)是二阶可微的．因此

G′(λ)=yh(y)(eλy-e-λy)dy>0,G″(λ)=y2h(y)(eλy+e-λy)dy>0．

设

Δ(c,λ)=cλ-pλ2-d[H(λ)-1]-∂1f(0,0)-∂2f(0,0)e-λcτG(λ),

(6)

(7)

对于c∈R,λ∈C，有c≥0，Reλ<λ+，若∂2f(0,0)>0，则λ+=λ; 若∂2f(0,0)=0，则λ+=+∞．

为了方便，记limξ→-∞φ(ξ)=φ(-∞), limξ→∞φ(ξ)=φ(+∞)．根据文献[15]中有关行波解存在性的结果，可以得出以下结论．

定理1假设条件(P1)、(P2)、(P3)和(P4)成立．存在最小波速c*>0，使得对任意c≥c*，方程(4)有满足行波方程(5)的波前解φ(x+ct)．当c∈(0,c*)时，方程(4)不存在满足方程(5)的波前解．因此

1) 方程(5)的解φ关于变换是唯一的．

2) 方程(5)的每个解φ是严格单调的,即φ′(ξ)>0，ξ∈R．

3) 方程(5)的每个解φ满足0<φ(ξ)<K，∀ξ∈R．

4) 方程(5)的任何解满足limξ→-∞(φ′(ξ)/φ(ξ))=λ，其中λ是方程(8)的最小正根，

cλ-pλ2-d[H(λ)-1]-∂1f(0,0)-∂2f(0,0)e-λcτG(λ)=0,

(8)

其中λ∈C, Reλ<λ0．

5) 方程(5)的任何解满足limξ→+∞(φ′(ξ)/(K-φ(ξ)))=γ，其中γ是方程(9)的唯一根，

cγ+pγ2+d[H(-γ)-1]+∂1f(K,K)+∂2f(K,K)eγcτG(-γ)=0,

γ∈C, Reγ<λ0．

(9)

注意到当c>c*时，特征方程(8)恰有两个正实根．下面给出两个引理．

引理1存在c*>0和λ*∈(0,λ+)，使得Δ(c*,λ*)=0且Δλ(c*,λ*)=0．因此，

1) 如果0<c<c*，则Δ(c,λ)<0，∀λ≥0．

2) 如果c>c*，则Δ(c,·)=0有两个正实根分别为λ1(c),λ2(c)，且满足0<λ1(c)<λ*<λ2(c)<λ+，使得当λ∈(λ1(c),λ2(c))时，Δ(c,λ)>0;当λ∈(-∞,λ1(c))∪(λ2(c),λ+)时,Δ(c,λ)>0．

证明对任意的λ∈(0,λ+)，有

Δλ(c,λ)=c-2pλ-dH′(λ)-∂2f(0,0)e-λcτ[G′(λ)-cτG(λ)],

Δλλ(c,λ)=-2p-dH″(λ)-

∂2f(0,0)e-λcτ[G″(λ)-2cτG′(λ)+c2τ2G(λ)]<0,

Δc(c,λ)=λ+cτ∂2f(0,0)e-λcτG(λ)>0,

Δ(c,0)=-∂1f(0,0)-∂2f(0,0)<0,

Δ(0,λ)=-pλ2-d[H(λ)-1]-∂1f(0,0)-∂2f(0,0)G(λ)<0,

并且

得证．

引理2假设条件(P1)和(P2)成立．对每个c≥0，方程恰有一个正实根

证明由条件(P1)、(P2)得

因此，至少有一个正的零点．由于

其中

H′(-λ)=yJ(y)(e-λy-eλy)dy<0, ∀λ>0,

G′(-λ)=yh(y)(e-λy-eλy)dy<0, ∀λ>0．

所以，对任意的即关于λ是单增的．因此恰有一个正实根．

2 比 较 原 理

方程(4)相关线性初值问题的比较原理是证明其波前解稳定性的关键[20，23，25]，下面首先证明比较原理成立．记

引理3假设p0,q0>0是任意常数，令w±(x,t):R×R+是两个连续函数且满足

1)w+(x,t)≥0且w-(x,t)≤K，其中x∈R，t≥-τ;

2)w+(x,t)≥w-(x,t)，其中

3) 对于

则有w+(x,t)≥w-(x,t)，∀x∈R,t≥0．

证明令W(x,t)=w+(x,t)-w-(x,t),∀x∈R,t≥-τ．那么W(x,t)有下界-K，且对任意的x∈R,s∈[-τ,0],W(x,s)≥0; 对任意的(x,t)∈R×[-τ,+∞)且x+ct≤L0，有W(x,t)≥0,而且对任意的有

Wt(x,t)≥pWxx(x,t)+d(J*W-W)(x,t)+p0W(x,t)+q0(h*W)(x,t)．

(10)

取

若假设不成立，则存在ε>0,t0>0使得

W(x,t)>-εe2Kt，x∈R,t∈[-τ,t0),

故存在一个正Lebesgue测度的有界集E⊂R，使得W(x,t0)≤-(15/16)εe2Kt0,x∈E．取有限区间[-a,a]，使得E⊆[-a,a]．令η(x)是一个光滑函数且满足

η(z)=1, ∀z∈[-a,a]，

对任意α∈[0,1],定义函数

z(x,t;α)

显然

∀x∈R,t∈[-τ,t0],

∀x∈R,t∈[-τ,t0],

∀x∈[-a,a]．

因此，令

α*∀t∈[-τ,t0],x∈R}．

显然

W(x,t)≥z(x,t;α*), ∀x∈R,t∈[-τ,t0]．

由于

W(x,s)≥0>z(x,s;α*), ∀x∈R,s∈[-τ,0]，

则

∀(x,t)∈R×[-τ,t0],

另外，当x+ct≤L0时,

W(x,t)≥0>z(x,t;α*), ∀(x,t)∈R×[-τ,t0]．

记Z(x,t)W(x,t)-z(x,t;α*)，不难发现，当x+ct>L0时,Z(x,t)在(x1,t1)∈R×(0,t0]处达到其下确界0．所以

Z(x1,t1)=0,Zt(x,t)|(x1,t1)≤0,

且

Zxx(x,t)|(x1,t1)≥0．

所以，对任意的y∈R,Z(x1-y,t1-τ)≥Z(x1,t1)=0,进而就有

W(x1-y,t1-τ)≥z(x1-y,t1-τ;α*)+W(x1,t1)-z(x1,t1;α*)≥

z(x1-y,t1-τ;α*)+W(x1,t1), ∀y∈R．

故由方程(10)可得

0≥Zt(x,t)|(x1,t1)=

Wt(x1,t1)-zt(x,t;α*)|(x1,t1)≥

-zt(x,t;α*)|(x1,t1)+pWxx(x1,t1)+d(J*W-W)(x1,t1)+

p0W(x1,t1)+q0(h*W)(x1,t1)≥

-2Kz(x1,t1;α*)+pzxx(x1,t1)+d(J*z-z)(x1,t1)+

p0z(x1,t1;α*)+q0(h*z)(x1,t1)≥

矛盾．因此

w+(x,t)≥w-(x,t)，∀(x,t)∈R×[-τ,+∞)．

证毕．

下面对方程(4)的特征值问题进行分析．首先考虑方程(4)在u=0处线性化的方程：

ut(x,t)=puxx(x,t)+d(J*u-u)(x,t)+

∂1f(0,0)u(x,t)+∂2f(0,0)(h*u)(x,t-τ),x∈R,

(11)

将u(x,t)=eλxv(t)代入方程(11)得

(12)

令v(t)=eμ(λ)tψμ是方程(12)的解,那么方程(11)对应的特征值问题

(13)

存在实根μ(λ)．由于u=0是不稳定平衡点，根据引理1可知μ(λ)>0．

另一方面，考虑方程(4)在u=K处线性化的方程:

ut(x,t)=puxx(x,t)+d(J*u-u)(x,t)+

∂1f(K,K)u(x,t)+∂2f(K,K)(h*u)(x,t-τ),x∈R．

(14)

方程(14)所对应的空间齐次方程为

ut(x,t)=∂1f(K,K)u(x,t)+∂2f(K,K)(h*u)(x,t-τ)，x∈R．

(15)

将u(x,t)=eλxν(t)代入方程(15),得

∂1f(K,K)ν(t)+∂2f(K,K)ν(t-τ)G(λ)．

(16)

令是方程(16)的解,则方程(15)的特征值问题

(17)

存在实根因为u=K是不稳定的平衡点，根据引理2可知

取则存在1>0使得

1w+1G(λ)w．

(18)

由于(φ(-∞),φ(+∞))=(0,K), 则对充分大的L0>0,

(∂1f(ζ1,ζ2),∂2f(ζ1,ζ2))≤(∂1f(0,0),∂2f(0,0)),x+ct≤L0，

(19)

(∂1f(ζ1,ζ2),∂2f(ζ1,ζ2))≤(∂1f(K,K)+1,∂2f(K,K)+1),

x+ct>L0．

(20)

通过构造恰当的权函数结合比较原理，讨论方程(4)的初始波动在加权最大范数空间上一致有界，满足此初值的解与行波解之间的关系．

3 稳 定 性

下面给出方程(4)非临界波前解渐近稳定性的结果．

定理2假设条件(P1)、(P2)、(P3)和(P4)成立，且在c>c*条件下，给定φ(x+ct)是方程(4)连接0和K的波前解．给定充分小的>0,令λ=λ1(c)+．定义如下权函数:

ω(x)

若初值u0满足

0≤u0(x,s)≤K，∀x∈R,s∈[-τ,0],

且

[u0(x,s)-φ(x+s)]ω(x)∈L∞(R), ∀s∈[-τ,0],

(21)

则方程(4)存在唯一解u(x,t;u0)并且满足

0≤u(x,t)≤K，∀(x,t)∈R×R+．

因此

(22)

其中μ0>0，M*>0为常数．

证明定义两个函数分别为

(23)

显然有

0≤U-(x,s)≤u0(x,s),φ(x+cs)≤U+(x,s)≤K,

∀(x,s)∈R×[-τ,0]．

令U±(x,t)分别是方程(4)具有初值U±(x,s)的解，则通过比较原理得

0≤U-(x，t)≤u(x,t;u0),φ(x+ct)≤U+(x,t)≤K, ∀x∈R,t>0．

于是

|u(x,t;u0)-φ(x+ct)|≤

max{|U+(x,t)-φ(x+ct)|, |U-(x，t)-φ(x+ct)|}．

因此，要证明此定理的结论成立，只需证明U±(x,t)指数收敛到φ(x+ct)．由对称性可知，仅证U+(x,t)收敛到φ(x+ct)即可．

记V(x,t)U+(x,t)-φ(x+ct)，显然有V(x,t)≥0,∀x∈R,t>0且

0≤V(x,s)≤|u0(x,s)-φ(x+cs)|, ∀x∈R,s∈[-τ,0],

于是，V(x,s)ω(x)在R上一致有界．接下来，分两种情况进行讨论．

情形1x+ct≤L0．由于V(x,t)≥0及U+(x,t)≤K, ∀x∈R,t≥0，则由方程(19)可得

Vt=pVxx+d(J*V-V)+f((U+),(h*U+))-f(φ,(h*φ))=

pVxx+d(J*V-V)+f((V+φ),(h*(V+φ)))-f(φ,(h*φ))≤

pVxx+d(J*V-V)+∂1f(0,0)V+∂2f(0,0)(h*V),

(24)

其中x∈R,t>0．令ψ是μ(λ)的特征向量(一维)．记μ(λ)=μ．由于V(x,s)ω(x)在R上是一致有界的，取充分大的C1>0，使得

C1ψeλ(x-L0)+μs≥V(x,s)，∀x∈R,s∈[-τ,0]．

定义

∀x∈R，t>0,

结合式(13)，容易验证

∀x∈R,t>0．

注意到∀x∈R,s∈[-τ,0],由引理3可得

∀x∈R,t∈[-τ,+∞)．

因此，对任意的(x,t)∈R×[-τ,+∞),若x+ct≤L0,则

V(x,t)≤C1ψeλ(x-L0)+μt=C1ψeλ(x+ct-L0)e-(cλ-μ)t≤

C1ψe-(cλ-μ)t．

(25)

情形2x+ct>L0．对任意的由式(20)得

Vt≤pVxx+d(J*V-V)+(∂1f(K,K)+1)V+(∂2f(K,K)+1)(h*V)．

(26)

注意到取μ0=min{cλ-μ,并选择足够大的C2>0使得

定义

∀x∈R，t∈[-τ,+∞),

(27)

则由情形1的结果以及得

∀

由式(18)得

于是

(28)

因此，函数满足

(29)

不妨令

p0=∂1f(K,K)+1,q0=∂1f(K,K)+1，

通过式(26)、(27)、(29)和引理3可以得出

∀(x,t)∈R×R+．

结合以上两种情形可得

0≤V(x,t)=U+(x,t)-φ(x+ct)≤M*e-μ0t, ∀x∈R，t≥-τ,

(30)

其中M*max{C1ψ，因此，U+(x,t)指数收敛到φ(x+ct)．

类似地，可以证明U-(x,t)指数收敛到φ(x+ct)．于是U(x,t)指数收敛到φ(x+ct)．

证毕．

4 应 用 实 例

这一节将主要通过具体的生物模型进一步验证所得稳定性的理论结果．

当方程(4)的反应项为Logistic型，即

那么，方程(4)就变为如下非局部时滞反应扩散方程：

ut(x,t)=puxx(x,t)+d(J*u-u)(x,t)+

(31)

其中x∈R,d≥0,p≥0．

对于模型(31)，条件(P3)和(P4)显然成立．容易发现f(0,0)=f(1,1)=0,f(u,u)=ru(1-u)>0, ∀u∈(0,1)，且对任意的(u,v)∈[0,1]2，∂2f(u,v)=r(1-u)≥0，即条件(P1)成立．

另一方面，总可以找到一个M>0和σ∈(0,1]使得对任意的(u,v)∈[0,1]2,

0≤∂1f(0,0)+∂2f(0,0)-f(u,v)≤M(u+v)1+σ,

并且∂1f(1,1)+∂2f(1,1)=-r<0．于是条件(P2)成立．

因此，对于c>c*，如果方程(31)的初始值u0满足

0≤u0(x,s)≤1，∀x∈R,s∈[-τ,0],

且

[u0(x,s)-φ(x+s)]ω(x)∈L∞(R), ∀s∈[-τ,0],

则方程(31)的唯一解u(x,t;u0)满足

0≤u(x,t)≤1，∀(x,t)∈R×R+．

因此

其中μ1>0，M1>0为常数．那么，定理2对于模型(31)仍然成立．

下面考虑非局部时滞Nicholson’s blowflies模型：

ut=puxx+d(J*u-u)-ru+

(32)

其中p≥0,d≥0,r>0,1<β≤e．f(u,v)=-ru+rpve-v满足先前的条件(P1)、(P2)、(P3)、(P4)，通过类似的分析可知，定理2(稳定性定理)对于模型(32)依然成立．于是，可以得到模型(32)所有非临界波前解的渐近稳定性结果．

5 总 结

本文研究一类具有非局部时滞反应项的空间非局部扩散模型．在单稳假设条件下，通过权函数结合比较原理，证明方程(4)的非临界波前解关于时间t是指数渐近稳定的．值得注意的是，该方法可以应用于更一般的非局部时滞反应扩散模型．本文仅得到c>c*时波前解的指数渐近稳定性结果，对于c=c*波前解的稳定性问题需要进一步去研究．另外,对于文献[15]中具有分布时滞的反应扩散模型行波解的稳定性也是未来研究的重点．

参考文献：

[1] FANG J, WEI J, ZHAO X Q.Spatial dynamics of a nonlocal and time-delayed reaction-diffusion system[J].JournalofDifferentialEquations, 2008,245(10): 2749-2770.

[2] GOURLEY S A, WU J.Delayed non-local diffusive systems in biological invasion and disease spread[J].NonlinearDynamics&EvolutionEquations, 2006,48: 137-200.

[3] OU C, WU J.Persistence of wavefronts in delayed nonlocal reaction-diffusion equations[J].JournalofDifferentialEquations, 2007,235(1): 219-261

[4] JOSEPH W H S, WU J, ZOU X.A reaction-diffusion model for a single species with age structure I: travelling wavefronts on unbounded domains[C]//Proceedings:Mathematical,PhysicalandEngineeringSciences, 2012,457: 1841-1853.

[5] THIEME H R, ZHAO X Q.Asymptotic speeds of spread and traveling waves for integral equations and delayed reaction-diffusion models[J].JournalofDifferentialEquations, 2003,195(2): 430-470.

[6] WANG Z C, LI W T, RUAN S.Traveling fronts in monostable equations with nonlocal delayed effects[J].JournalofDynamics&DifferentialEquations, 2008,20(3): 573-607.

[7] BRITTON N F.Spatial structures and periodic travelling waves in an integro-differential reaction-diffusion population model[J].SIAMJournalonAppliedMathematics, 1990,50(6): 1663-1688.

[8] BRITTON N F.Aggregation and the competitive exclusion principle[J].JournalofTheoreticalBiology, 1989,136(1): 57-66.

[9] BATES P W, FIFE P C, REN X, et al.Traveling waves in a convolution model for phase transitions[J].ArchiveforRationalMechanics&Analysis, 1997,138(2): 105-136.

[10] BATES P W, CHMAJ A.A discrete convolution model for phase transitions[J].ArchiveforRationalMechanics&Analysis, 1999,150(4): 281-368.

[11] COVILLE J.On uniqueness and monotonicity of solutions of non-local reaction diffusion equation[J].AnnalidiMatematicaPuraedApplicata, 2006,185(3): 461-485.

[12] COVILLE J, DUPAIGNE L.Propagation speed of travelling fronts in non local reaction-diffusion equations[J].NonlinearAnalysis:Theory,Methods&Applications, 2005,60(5): 797-819.

[13] COVILLE J, DUPAIGNE L.On a non-local eqution arising in population dynamics[J].ProceedingsoftheRoyalSocietyofEdinburghSectionA:Mathematics, 2007,137(4): 727-755.

[14] YU Z, YUAN R.Existence, asymptotic and uniqueness of traveling waves for nonlocal diffusion systems with delayed nonlocal response[J].TaiwaneseJournalofMathematics, 2013,17(6): 2163-2190.

[15] GUO S, ZIMMER J.Travelling wavefronts in nonlocal diffusion equations with nonlocal delay effects[J].BulletinoftheMalaysianMathematicalSciencesSociety, 2018,41(2): 919-943.

[16] CHENG H, YUAN R.Stability of traveling wave fronts for nonlocal diffusion equation with delayed nonlocal response[J].TaiwaneseJournalofMathematics, 2016,20(4): 801-822.

[17] GUO S, ZIMMER J.Stability of travelling wavefronts in discrete reaction-diffusion equations with nonlocal delay effects[J].Nonlinearity, 2015,28(2): 463-492.

[18] MEI M, OU C, ZHAO X Q.Global stability of monostable traveling waves for nonlocal time-delayed reaction-diffusion equations[J].SIAMJournalonMathematicalAnalysis, 2010,42(6): 2762-2790.

[19] MA S, ZHAO X Q.Global asymptotic stability of minimal fronts in monostable lattice equations[J].DiscreteandContinuousDynamicalSystems, 2008,21(1): 259-275.

[20] OUYANG Z, OU C.Global stability and convergence rate of traveling waves for a nonlocal model in periodic media[J].DiscreteandContinuousDynamicalSystems(SeriesB), 2012,17(3): 993-1007.

[21] 邢伟, 颜七笙, 杨志辉, 等.一类具有非线性传染率的SEIS传染病模型的稳定性分析[J].应用数学和力学, 2016,37(11): 1247-1254.(XING Wei, YAN Qisheng, YANG Zhihui.Stability analysis of an SEIS epidemic model with a nonlinear incidence rate[J].AppliedMathematicsandMechanics, 2016,37(11): 1247-1254.(in Chinese))

[22] 谢英超, 程燕, 贺天宇.一类具有非线性发生率的时滞传染病模型的全局稳定性[J].应用数学和力学, 2015,36(10): 1107-1116.(XIE Yingchao, CHENG Yan, HE Tianyu.Global stability of a class of delayed epidemic models with nonlinear incidence rates[J].AppliedMathematicsandMechanics, 2015,36(10): 1107-1116.(in Chinese))

[23] WU S L, CHEN G.Uniqueness and exponential stability of traveling wave fronts for a multi-type SIS nonlocal epidemic model[J].NonlinearAnalysis:RealWorldApplications, 2017,36: 267-277.

[24] CHEN X, GUO J S.Existence and asymptotic stability of traveling waves of discrete quasilinear monostable equations[J].JournalofDifferentialEquations, 2002,184(2): 549-569.

[25] WANG X S, ZHAO X Q.Pulsating waves of a partially degenerate reaction-diffusion system in a periodic habitat[J].JournalofDifferentialEquations, 2015,259(12): 7238-7259.

[26] WU Y, XING X.Stability of traveling waves with critical speeds forp-degree Fisher-type equations[J].DiscreteandContinuousDynamicalSystems, 2008,20(4): 1123-1139.

Exponential Stability of Traveling Wavefronts for Reaction-Diffusion Equations With Delayed Nonlocal Responses①

LI Panxiao

(SchoolofMathematicsandStatistics,XidianUniversity,Xi’an710071,P.R.China)

Abstract:A spatially nonlocal diffusion model with a class of delayed nonlocal responses was considered.The asymptotic stability and the convergence rate of the traveling wavefronts were mainly studied.Through construction of weighted functions and establishment of a comparison principle for the related linear equations, the conclusion that if the initial function is within a bounded distance from a certain traveling wavefront with respect to a weighted maximum norm, the solution satisfying the initial value will converge to the traveling wavefront exponentially in time, was proved, and the exponential convergence rate was also obtained.

Key words:traveling wave solution; nonlocal time delay; stability

Foundation item:The National Natural Science Foundation of China(11671315)

作者简介：李盼晓(1990—),女,硕士(E-mail: L2236135404@163.com).

基金项目：国家自然科学基金(11671315)

修订日期：2018-02-02

*收稿日期：2017-12-28；

文章编号：1000-0887(2018)11-1300-13

文献标志码：A DOI: 10.21656/1000-0887.380336

中图分类号：O175.14

① 引用本文/Citethispaper:

李盼晓.非局部时滞反应扩散方程波前解的指数稳定性[J].应用数学和力学, 2018,39(11): 1300-1312.

LI Panxiao.Exponential stability of traveling wavefronts for reaction-diffusion equations with delayed nonlocal responses[J].AppliedMathematicsandMechanics, 2018,39(11): 1300-1312.