细菌降解宿主组织模型波前解的全局稳定性
Global Stability of Traveling Wavefronts for a Model Describing Host Tissue Degradation by Bacteria
DOI: 10.12677/pm.2025.1512290, PDF, HTML, XML,    国家自然科学基金支持
作者: 张 敏, 张国宝*:西北师范大学数学与统计学院,甘肃 兰州
关键词: 宿主组织降解细菌感染波前解全局稳定性Host Tissue Degradation Bacterial Infection Traveling Wavefronts Global Stability
摘要: 对于细菌降解宿主组织模型,Yang-Zhang-Tian已研究了波前解的全局稳定性。然而,其稳定性结论需要较为严格的参数条件,导致参数取值范围较小。本文通过采用傅里叶变换以及具有适当加权函数的加权能量方法,证明了在其他参数条件下波前解依然稳定,扩大了参数的可取范围。
Abstract: For a model describing bacterial degradation of host tissues by bacteria, the global stability of the traveling wavefronts were proved by Yang-Zhang-Tian. However, the stability result needs more stronger conditions on the parameters, which may result in a small range of parameters. In this paper, by using Fourier’s transform and the weighted energy method with a suitably selected weighted function, we prove that the wavefronts remain stable under other conditions, which can increase the range of parameters.
文章引用:张敏, 张国宝. 细菌降解宿主组织模型波前解的全局稳定性[J]. 理论数学, 2025, 15(12): 24-33. https://doi.org/10.12677/pm.2025.1512290

1. 引言

本文研究如下的细菌降解宿主组织模型

{ u t =Δuu+ωγku( 1ω ), ω t =ku( 1ω ), (1.1)

初值为 u( x,0 )= u 0 ( x ) ω( x,0 )= ω 0 ( x ) ,其中 γ,k>0 为常数, k>0 表示组织降解速率,通常很大,是关注的关键参数。变量 u 表示由细菌产生的降解性酶的浓度,而 1ω 表示健康组织的体积分数。该模型由一个反应扩散方程与一个常微分方程组成。King等人[1]在首次研究了这一模型,刻画了胞外细菌病原体对宿主组织的降解过程。细菌对组织的侵袭表现为以近似常速传播的行波。Hilhorst,King与Röger [2]证明了(1.1)的解的存在唯一性,并建立了当降解率参数 k 时,原系统的解收敛到一个类似Stefan问题的唯一极限解。

模型(1.1)对应的齐次(空间均匀)系统为

{ u =u+ωγku( 1ω ), ω =ku( 1ω ). (1.2)

通过简单计算可知,系统(1.2)有两个平衡点: ( 0,0 ) ( 1,1 ) 。此外,容易得到 ( 0,0 ) 是不稳定的。

模型(1.1)的行波解是形如 ( u( x,t ),ω( x,t ) )=( φ 1 ( ξ ), φ 2 ( ξ ) ) 的连接 ( 0,0 ) ( 1,1 ) 的特殊平移不变解,其中 ξ=x+ct 。如果行波解是单调的,则称其为波前解。将 ( φ 1 , φ 2 ) 代入(1.1),得到如下的波形方程

{ c φ 1 = φ 1 φ 1 + φ 2 γk φ 1 ( 1 φ 2 ), c φ 2 =k φ 1 ( 1 φ 2 ) (1.3)

并满足边界条件:

lim ξ ( φ 1 ( ξ ), φ 2 ( ξ ) )=( 0,0 ), lim ξ+ ( φ 1 ( ξ ), φ 2 ( ξ ) )=( 1,1 ).

Hilhorst、King与Röger [3]证明了波前解的存在性。他们指出当降解速率足够大时,最小传播速度与极限问题的最小传播速度一致。Zhang等人[4]通过构造恰当的上下解,给出了该系统最小波速的线性选择条件。彭和易[5]利用比较原理和渐近分析的方法建立了波前解的局部渐近行为。He、He与Pan [6]基于适当加权函数空间中的谱方法,证明了波前解的局部稳定性。最近,Yang、Zhang与Tian [7]在条件(H)

下(即 γ<1 1 2γ <k< 1 γ )利用加权能量方法与比较原理,建立了大波速情形下波前解的全局指数稳定性。需要指出的是,假设条件 k< 1 γ 比较严格。受文献[8]-[10]的启发,本文采用傅里叶变换与具有适当加

权函数的加权能量方法,证明在第2节所给的假设条件(A)下,波前解依然稳定,从而扩大了参数 k 的取值范围。

2. 预备知识和主要结果

为证明行波解的稳定性,我们需要作如下假设:

(A) γ<1 1 1γ <k< 2 1γ

与文献[7]中的假设(H)相比较,我们发现当 γ1/2 时,条件(A)可以扩大参数 k 的取值范围。

对于单调系统(1.1),显然成立如下的比较原理。

引理2.1 (比较原理)设 ( u ( x,t ), ω ( x,t ) ) ( u + ( x,t ), ω + ( x,t ) ) 分别为方程(1.1)在初值 ( u 0 ( x ), ω 0 ( x ) ) ( u 0 + ( x ), ω 0 + ( x ) ) 下的解。若对一切 xR 都有

( 0,0 )( u 0 ( x ), ω 0 ( x ) )( u 0 + ( x ), ω 0 + ( x ) )( 1,1 ),

则对一切 xR,t>0 都有

( 0,0 )( u ( x,t ), ω ( x,t ) )( u + ( x,t ), ω + ( x,t ) )( 1,1 ).

C>0 表示通用常数, C i >0,( i=1,2, ) 表示特定常数。 L ( R ) R 上的本性有界函数空间,其范

数定义为 | f | L = sup xR | f( x ) | L ω ( R ) 为加权 L ( R ) 空间,其权函数为 ω( x )>0 ,范数定义为 | f | L ω = sup xR | ω( x )f( x ) | l 2 为平方可和函数空间,其范数定义为 | z | 2 = ( i=1 N | z i | 2 ) 1/2 ,其中 z=( z 1 , z 2 ,, z N )

c c min ,设 ( φ 1 ( x+ct ), φ 2 ( x+ct ) ) 是系统(1.1)的一个波前解,其中波速为 c ,并连接 ( 0,0 ) ( 1,1 ) 。我们定义如下权函数:

w( x )={ e λ * ( x ζ 0 ) , x ζ 0 , 1, x> ζ 0 ,

其中 ζ 0 为充分大的常数,且 λ * = 2+( γ1 )k 。由于波的稳定性问题的困难在于不稳定平衡点 ( 0,0 ) ,因此,在利用能量方法的证明中我们引入指数权函数来克服这个困难。这里 λ * 的选取是技术性的。

定理2.2 (稳定性)假设条件(A)成立。任取系统(1.1)的一个连接 ( 0,0 ) ( 1,1 ) 的波前解

( φ 1 ( x+ct ), φ 2 ( x+ct ) ) ,其波速满足 c> c * :=max{ c min , 1+( γ+1 )k 2 2+( γ1 )k } 。若初值 ( u 0 ( x ), ω 0 ( x ) ) 满足

( 0,0 )( u 0 ( x ), ω 0 ( x ) )( 1,1 ),xR,

且扰动满足

( u 0 φ 1 , ω 0 φ 2 ) L w ( R )× L w ( R ),

则(1.1)的具有初值 u( x,0 )= u 0 ( x ),w( x,0 )= ω 0 ( x ) 的解 ( u( x,t ),v( x,t ) ) 整体存在,满足

( 0,0 )( u( x,t ),ω( x,t ) )( 1,1 ),xR,t>0,

且按时间指数收敛到该波前解 ( φ 1 ( x+ct ), φ 2 ( x+ct ) ) ,即存在常数 σ>0 C>0 ,使得

sup xR | u( x,t ) φ 1 ( x+ct ) |C e σt ,

sup xR | ω( x,t ) φ 2 ( x+ct ) |C e σt

注记:文献[7]中的 c * 为传播速度,与定理2.2中的 c * 不同。本文的 c * 与文献[7]定理2.6中的 max{ c min , c ˜ } 角色一致。在参数 γ k 的不同取值下,本文的 c * 可能会比文献[7]定理2.6中的 max{ c min , c ˜ } 小,也可能会大。

3. 行波解的稳定性

在本节中,我们致力于证明(1.1)的行波解的稳定性,即定理2.2。我们首先考虑如下线性(含时滞的)微分系统:

{ d dt z( t )=Az( t ),t>0, z( 0 )= z 0 , (3.1)

其中 A C N×N ,显然,系统(3.1)的解可表示为

z( t )= z 0 e At ,

这里 e At := m=0 ( At ) m m!

需要说明的是,在文献[9] [11]中,作者分别在 N=1 N2 的情况下研究了带时滞的系统(3.1)。我们定义

μ 1 ( A )= lim θ 0 + | I+θA |1 θ = max 1jN [ Re( a jj )+ ji N | a ij | ]

以及

μ ( A )= lim θ 0 + | I+θA | 1 θ = max 1iN [ Re( a jj )+ ji N | a ij | ].

这些定义见([9],定理3.1)的证明。

引理3.1 假设矩阵 A 的矩阵度量满足

μ( A ):= μ 1 ( A )+ μ ( A ) 2 <0 ,

其中 μ 1 ( A ) μ ( A ) 分别是由1-范数与∞-范数诱导的矩阵度量。则对系统(3.1)任意解 z( t ) ,存在一个常数 C 0 >0 ,使得

| z( t ) | C 0 e αt ,t>0, (3.2)

其中 α=| μ( A ) |>0 ,特别地

| e At | C 0 e αt ,t>0.

证明:

W( z( ),t )= | z( t ) | 2 2 = i=1 2 | z i ( t ) | 2 .

则沿着系统(3.1)的方向导数(李雅普诺夫导数) DW

DW( z( ),t )= i=1 N [ z ¯ i ( t )z( t )+ z i ( t ) z ¯ i ( t ) ] = i=1 N { z ¯ i ( t ) j=1 N a ij z j ( t )+ z i ( t ) j=1 N a ¯ ij z ¯ j ( t ) } = i=1 N { a ii z ¯ i ( t ) z i ( t )+ j = i N a ij z ¯ i ( t ) z j ( t )+ a ¯ ii z i ( t ) z ¯ i ( t )+ j = i N a ¯ ij z i ( t ) z ¯ j ( t ) } = i=1 N { 2Re( a ii ) | z i ( t ) | 2 +2 j = i N R e( a ij z ¯ i ( t ) z j ( t ) ) } i=1 N { 2Re( a ii ) | z i ( t ) | 2 + j = i N | a ij | ( | z i ( t ) | 2 + | z j ( t ) | 2 ) } ( μ 1 ( A )+ μ ( A ) ) | z( t ) | 2 2 =( μ 1 ( A )+ μ ( A ) )W( z( ),t ). (3.3)

对不等式(3.3)两边积分可得

W( z( ),t )W( z( ),0 ) e 2μ( A )t ,t0

注意到对某个 N>0 常数有 | z( t ) |N | z( t ) | 2 于是立刻推出式(3.2)成立。证明完毕。

为得到所需结论,我们需要先对系统(1.1)解的差建立先验估计。设 ( u( x,t ),ω( x,t ) ) 是系统(1.1)在初值 ( u 0 ( x ), ω 0 ( x ) ) 下的解。定义(对所有 xR ):

u 0 ( x )=min{ u 0 ( x ), φ 1 ( x ) }, ω 0 ( x )=max{ ω 0 ( x ), φ 2 ( x ) },

u 0 + ( x )=min{ u 0 ( x ), φ 1 ( x ) }, ω 0 + ( x )=max{ ω 0 ( x ), φ 2 ( x ) }.

容易得出对所有 xR ,有

0 u 0 ( x ) u 0 ( x ) u 0 + ( x )1,

0 u 0 ( x ) φ 1 ( x ) u 0 + ( x )1,

0 ω 0 ( x ) ω 0 ( x ) ω 0 + ( x )1,

0 ω 0 ( x ) φ 2 ( x ) ω 0 + ( x )1.

( u ± ( x,t ), ω ± ( x,t ) ) 分别是初值为 ( u 0 ± ( x ), ω 0 ± ( x ) ) 的系统(1.1)的解。由引理2.1的比较原理得,对于所有 xR,t>0 ,有

0 u ( x,t )u( x,t ) u + ( x,t )1,

0 u ( x,t ) φ 1 ( x+ct ) u + ( x,t )1,

0 ω ( x,t )ω( x,t ) ω + ( x,t )1,

0 ω ( x,t ) φ 2 ( x+ct ) ω + ( x,t )1.

接下来,我们将分三步完成定理2.2的证明。

步骤1:当任意给定 c c * 时, ( u + ( x,t ), ω + ( x,t ) ) 收敛到 ( φ 1 ( x+ct ), φ 2 ( x+ct ) )

ξ:=x+ct ,并设

U 1 ( ξ,t )= u + ( x,t ) φ 1 ( x+ct ), V 1 ( ξ,t )= ω + ( x,t ) φ 2 ( x+ct ) ,

其初值为 U 10 ( ξ,0 )= u 0 + ( x ) φ 1 ( x ) V 10 ( ξ,0 )= ω 0 + ( x ) φ 2 ( x )

容易验证

{ U 1t +c U 1ξ =Δ U 1 U 1 + V 1 γk U 1 +γk U 1 ( V 1 + φ 2 )+γk φ 1 V 1 , V 1t +c V 1ξ =k U 1 k U 1 ( V 1 + φ 2 )k φ 1 V 1 . (3.4)

由于 V 1 + φ 2 ( 0,1 ) φ 1 ( 0,1 ) ,有

U 1 + V 1 γk U 1 +γk U 1 ( V 1 + φ 2 )+γk φ 1 V 1 U 1 +( 1+γk ) V 1 ,

k U 1 k U 1 ( V 1 + φ 2 )k V 1 φ 1 k U 1 .

由(3.4)得

{ U 1t +c U 1ξ Δ U 1 U 1 +( 1+γk ) V 1 , V 1t +c V 1ξ k U 1 . (3.5)

( U 1 + ( ξ,t ), V 1 + ( ξ,t ) ) 为下述柯西问题的解

{ U 1t + +c U 1ξ + =Δ U 1 + U 1 + +( 1+γk ) V 1 + , V 1t + +c V 1ξ + =k U 1 + , ( U 1 + ( ξ,0 ), V 1 + ( ξ,0 ) )=( U 10 ( ξ,0 ), V 10 ( ξ,0 ) ). (3.6)

由比较原理得

( U 1 ( ξ,t ), V 1 ( ξ,t ) )( U 1 + ( ξ,t ), V 1 + ( ξ,t ) ),( ξ,t )R×[ 0,+ ). (3.7)

U ˜ ( ξ,t )= e λ * ( ξ ζ 0 ) U 1 + ( ξ,t ), V ˜ ( ξ,t )= e λ * ( ξ ζ 0 ) V 1 + ( ξ,t ) ,

其中 ζ 0 为充分大的常数(具体取值稍后给出)。由(3.6)可得 ( U ˜ ( ξ,t ), V ˜ ( ξ,t ) ) 满足

{ U ˜ t +( c2 λ * ) U ˜ ξ = U ˜ ξξ +( λ * 2 c λ * 1 ) U ˜ +( 1+γk ) V ˜ , V ˜ t +c V ˜ ξ =k U ˜ c λ * V ˜ . (3.8)

对系统(3.8)作傅里叶变换,得

{ U ^ t ( η,t )=[ η 2 iη( c2 λ * )+ λ *2 c λ * 1 ] U ^ ( η,t )+( 1+γk ) V ^ ( η,t ), V ^ t ( η,t )=k U ^ ( η,t )( c λ * +icη ) V ^ ( η,t ), (3.9)

其中

U ^ ( η,t )= U ˜ ( ξ,t ) e iηξ dξ , V ^ ( η,t )= V ˜ ( ξ,t ) e iηξ dξ .

Z( ξ,t )= ( U ˜ ( ξ,t ), V ˜ ( ξ,t ) ) T , Z ^ ( η,t )= ( U ^ ( η,t ), V ^ ( η,t ) ) T .

于是,系统(3.9)可改写为如下微分系统:

Z ^ ( η,t ) t =A( η ) Z ^ ( η,t ), (3.10)

其中,

A( η )=( η 2 iη( 2 λ * +c )+ λ *2 c λ * 1 1+γk k c λ * icη ).

(3.10)的解可表示为

Z ^ ( η,t )= e A( η )t Z ^ 0 ( η ).

因此,取逆傅里叶变换可得

Z( ξ,t )= 1 2π e iηξ e A( η )t Z ^ 0 ( η )dη .

引理3.2 成立如下结论

| U 1 ( ,t ) | L ( Ω ) + | V 1 ( ,t ) | L ( Ω ) C e αt ,

其中, Ω :=( , ζ 0 ]

证明:由于 λ * = 2+( γ1 )k ,且条件(A)成立,故经计算可得

μ 1 ( A( η ) )=max{ η 2 + λ *2 c λ * 1+k,c λ * +1+γk }=c λ * +1+γk

以及

μ ( A( η ) )=max{ η 2 + λ *2 c λ * +γk,c λ * +k }=c λ * +k.

因此

μ( A( η ) )= μ 1 ( A( η ) )+ μ ( A( η ) ) 2 =c λ * + 1+( γ+1 )k 2 <0,

由于

c> c * > 1+( γ+1 )k 2 λ * = 1+( γ+1 )k 2 2+( γ1 )k ,

故由引理3.1可得

| e A( η )t | C 0 e αt ,t>0. (3.11)

进一步,利用(3.11)得

sup ξR Z( ξ,t )= sup ξR 1 2π + e iξη e A( η )t Z ^ 0 ( η )dη C 0 e αt sup ξR 1 2π + e iξη Z ^ 0 ( η )dη C 0 e αt sup ξR Z 0 ( ξ,t ) .

由此可以推出

| U ˜ ( ,t ) | L ( R ) + | V ˜ ( ,t ) | L ( R ) C 1 e αt , (3.12)

其中 C 1 := C 0 ( | U 0 | L ( Ω ) + | V 0 | L ( Ω ) ) 。由(3.7)以及当 ξ ζ 0 时, e λ * ( ξ ζ 0 ) 1 ,得到

0 U 1 ( ξ,t ) U 1 + ( ξ,t )= e λ * ( ξ ζ 0 ) U ˜ ( ξ,t ) U ˜ ( ξ,t ),

0 V 1 ( ξ,t ) V 1 + ( ξ,t )= e λ * ( ξ ζ 0 ) V ˜ ( ξ,t ) V ˜ ( ξ,t )

对所有 ξ Ω :=( , ζ 0 ] 成立。由(3.12)立即得到

U 1 ( ,t ) L ( Ω ) + V 1 ( ,t ) L ( Ω ) C e αt ,

其中 C>0 为常数。证明完毕。

引理3.3 若 c> c * ,则存在一个足够大的数 ζ 0 >0 ,使得

sup ξ[ ζ 0 ,+ ) U 1 ( ξ,t )C e βt , sup ξ[ ζ 0 ,+ ) V 1 ( ξ,t )C e βt ,t>0, (3.13)

其中, β=min{ 1,( 1γ )k1 }

证明:注意到 ( U 1 , V 1 ) 满足

{ U 1t +c U 1ξ Δ U 1 U 1 +( 1+γk ) V 1 , V 1t +c V 1ξ k U 1 k U 1 φ 2 k φ 1 V 1 . (3.14)

对(3.14)的第一个不等式取极限 ξ ,并注意到由于 ( ξ,t )R×[ 0,+ ) U 1 ( ξ,t ) 有界,故 U 1ξ ( ,t )=0 Δ U 1 ( ,t )=0 ,于是得

d dt U 1 ( ,t ) U 1 ( ,t )+( 1+γk ) V 1 ( ,t ) (3.15)

对(3.15)在区间 [ 0,t ] 上积分,得到

U 1 ( ,t ) U 1 ( ,0 ) 0 t U 1 ( ,s )ds+( 1+γk ) 0 t V 1 ( ,s )ds. (3.16)

同理,对(3.14)的第二条不等式令 ξ ,并注意到 V 1 ( ξ,t ) ( ξ,t )R×[ 0,+ ) 有界从而 V 1ξ ( ,t )=0 ,且 φ i ( )=1,i=1,2 ,于是得

d dt V 1 ( ,t )k V 1 ( ,t ) (3.17)

对(3.17)在区间 [ 0,t ] 上积分得

V 1 ( ,t ) V 1 ( ,0 )k 0 t V 1 ( ,s )ds (3.18)

将(3.16)与(3.18)相加可得

U 1 ( ,t )+ V 1 ( ,t )[ U 1 ( ,0 )+ V 1 ( ,0 ) ] 0 t U 1 ( ,s )ds( k1γk ) 0 t V 1 ( ,s )ds. (3.19)

由条件(A)有 k> 1 1γ 。显然, k1γk=( 1γ )k1>0 ,选取 β=min{ 1,( 1γ )k1 }>0 ,由(3.19)推得

U 1 ( ,t )+ V 1 ( ,t )β[ 0 t U 1 ( ,s )ds + 0 t V 1 ( ,s )ds ]+ C 2 (3.20)

其中 C 2 = U 1 ( ,0 )+ V 1 ( ,0 ) ,由Gronwall不等式,(3.20)推出

U 1 ( ,t )+ V 1 ( ,t ) C 3 e βt . (3.21)

这意味着

lim ξ U 1 ( ξ,t )C e βt , lim ξ V 1 ( ξ,t )C e βt ,t>0.

因此,存在一个足够大的数 ζ 0 >0 ,使得(3.13)成立。证明完毕。

基于引理3.2与引理3.3,立即得到如下引理。

引理3.4 对任意,有

U 1 ( ,t ) L ( R ) + V 1 ( ,t ) L ( R ) C e σt ,t>0,

其中, 0<σ<min{ α,β }

U 1 ( ξ,t ) V 1 ( ξ,t ) 的定义,可得解 ( u + ( x,t ), ω + ( x,t ) ) 具有如下收敛性。

引理3.5 当 c> c *

sup xR | u + ( x,t ) φ 1 ( x+ct ) |C e σt ,t>0,

sup xR | ω + ( x,t ) φ 2 ( x+ct ) |C e σt ,t>0,

其中, 0<σ<min{ α,β }

步骤2: ( u ( x,t ), ω ( x,t ) ) 收敛到 ( φ 1 ( x+ct ), φ 2 ( x+ct ) )

对任意给定 c> c * ,令 ξ=x+ct ,并设

U 2 ( ξ,t )= φ 1 ( x+ct ) u ( x,t ), V 2 ( ξ,t )= φ 2 ( x+ct ) ω ( x,t ),

其初值为 U 20 ( ξ,0 )= φ 1 ( x ) u 0 ( x ), V 20 ( ξ,0 )= φ 2 ( x ) ω 0 ( x ) ,其中 t>0,ξR 。与步骤1同理,可类似证明 ( u ( x,t ), ω ( x,t ) ) 收敛到 ( φ 1 ( x+ct ), φ 2 ( x+ct ) )

引理3.6 对于任意 c> c * ,存在一个正数 σ ,使得

sup xR | u ( x,t ) φ 1 ( x+ct ) |C e σt ,t>0,

sup xR | ω ( x,t ) φ 2 ( x+ct ) |C e σt ,t>0,

其中 0<σ<min{ α,β }

步骤3: ( u( x,t ),ω( x,t ) ) 收敛到 ( φ 1 ( x+ct ), φ 2 ( x+ct ) )

定理2.2的证明:因为

( 0,0 )( u ( x,t ), ω ( x,t ) )( u( x,t ),ω( x,t ) )( u + ( x,t ), ω + ( x,t ) )( 1,1 )

( 0,0 )( u ( x,t ), ω ( x,t ) )( φ 1 ( x+ct ), φ 2 ( x+ct ) )( u + ( x,t ), ω + ( x,t ) )( 1,1 ),

其中 xR,t>0 ,所以由引理3.5与引理3.6并利用夹逼准则论证,可立即得到:对任意 c> c * ,存在常数 σ>0,C>0 使得

sup xR | u( x,t ) φ 1 ( x+ct ) |C e σt ,t>0,

sup xR | ω( x,t ) φ 2 ( x+ct ) |C e σt ,t>0.

定理2.2证毕。

基金项目

国家自然科学基金(12261081, 11861056)。

NOTES

*通讯作者。

参考文献

[1] King, J.R., Koerber, A.J., Croft, J.M., et al. (2003) Modelling Host Tissue Degradation by Extracellular Bacterial Pathogens. Mathematical Medicine and Biology, 20, 227-260. [Google Scholar] [CrossRef
[2] Hilhorst, D., King, J.R. and Röger, M. (2007) Mathematical Analysis of a Model Describing the Invasion of Bacteria in Burn Wounds. Nonlinear Analysis: Theory, Methods & Applications, 66, 1118-1140. [Google Scholar] [CrossRef
[3] Hilhorst, D., King, J.R. and Röger, M. (2007) Travelling-Wave Analysis of a Model Describing Tissue Degradation by Bacteria. European Journal of Applied Mathematics, 18, 583-605. [Google Scholar] [CrossRef
[4] Zhang, T., Chen, D., Han, Y. and Ma, M. (2021) Linear Determinacy of the Minimal Wave Speed of a Model Describing Tissue Degradation by Bacteria. Applied Mathematics Letters, 121, Article 107044. [Google Scholar] [CrossRef
[5] 彭叠, 易亚婷. 细菌降解宿主组织模型行波解的局部渐近行为[J]. 理论数学, 2022(12): 1966-1970.
[6] He, X., He, G.S. and Pan, C.H. (2022) Local Stability of Traveling Waves of a Model Describing Host Tissue Degradation by Bacteria. Mathematics, 10, Article 1805. [Google Scholar] [CrossRef
[7] Yang, X.X., Zhang, G.B. and Tian, G. (2024) The Dynamics of Traveling Wavefronts for a Model Describing Host Tissue Degradation by Bacteria. International Journal of Biomathematics, 17, Article 2350031. [Google Scholar] [CrossRef
[8] Yang, Z.J., Zhang, G.B. and Tian, G. (2024) Global Stability of Traveling Wave Solutions for a Discrete Diffusion Epidemic Model with Nonlocal Delay Effects. Journal of Mathematical Physics, 66, Article 021508. [Google Scholar] [CrossRef
[9] Ma, Z.H., Yuan, R., Wang, Y. and Wu, X. (2019) Multidimensional Stability of Planar Traveling Waves for the Delayed Nonlocal Dispersal Competitive Lotka-Volterra System. Communications on Pure & Applied Analysis, 18, 2069-2092. [Google Scholar] [CrossRef
[10] Mei, M., Zhang, K.J. and Zhang, Q.F. (2019) Global Stability of Traveling Waves with Oscillations for Nicholson’s Blowflies Equation. International Journal of Numerical Analysis and Modeling, 16, 375-397.
[11] Khusainov, D.Y., Ivanov, A.F. and Kovarzh, I.V. (2009) Solution of One Heat Equation with Delay. Nonlinear Oscillations, 12, 260-282. [Google Scholar] [CrossRef