众所周知,时滞问题经常在控制系统中遇到,在过去的几十年里, 稳定性分析受到了广泛关注[1, 2]. 在数学上,研究小时滞对动力 学系统的影响也是很重要的. 这个问题对于线性系统,包括有限维和 无限维的情况,都得到了很好的研究[3, 4, 5, 6, 7]. 然而对于非线性系统, 问题要复杂得多,但是也有非常好的结果[8, 9, 10, 11, 12, 13, 14, 15, 16, 17].
在本文中,我们考虑如下一般的具有小时滞的非线性梯度系统 ∂tu+Au=f(u(t),u(t−τ)). 这里梯度系统是指对应的非时滞系统 ∂tu+Au=f(u,u) 是梯度系统,这里A是Banach空间中的扇形算子, 我们的工作空间是分数次幂空间Xα.
一般来讲,小时滞通常不会改变平衡点的局部渐近稳定性,例如, 在文献[11]中,Li和Kloeden指出,如果E是非时滞系统的局部渐近稳 定平衡点,则时滞系统在平衡点E的一个小邻域内有一个紧的局部渐近稳 定集. Mao在文献[14]中也证明了具有全局Lipschitz项的时滞非线性系统, 平衡点是全局指数渐近稳定的. 而且,小时滞的梯度系统通常和无时滞的 梯度系统行为非常相似. 例如,在参考文献[4]中,作者证明了时滞微分方 程的解和非时滞的具有相同的收敛性质. 相似的结果也可以在文献[12] 中找到,它指出,对于耗散的时滞反应扩散方程,所有的解都趋于平衡点, 当时间趋于无穷大时. 在最近的工作中[9],Li和Wang得到了一 个非常漂亮的结果: 如果平衡点是孤立双曲的,一般的时滞梯度系统的 每一个有界解都会精确地走向平衡点.
在这篇文章中,我们研究无限维情形下的时滞梯度系统. 在有限个和孤立平衡点的假设下,我们证明了一定存在一个足够小的时滞 使得时滞方程的任一个有界解将会最终进入并停留在某一个平衡点的邻域 里面,我们的方法主要是利用梯度系统不变集的Morse结构. 此外我们还得到一个更加有意思的结果: 在双曲平衡点的假设下, 一定存在ε>0 和足够小的 τ>0使得任何一个落于 某个平衡点ε -邻域内的解最终收敛于该平衡点, 当时间趋于无穷大时. 由以上两个结果,我们可以得到, 如果系统的平衡点是有限的双曲孤立平衡点,则方程(1.1)每一个有界 解将会精确地走向某个平衡点.
下面我们将给出一些基本的概念和结论,首先给出以下基本假设.
(A1)~ A是Banach空间X上的正扇形算子,并且有紧的预解式, e−At是由−A生成的解析半群. 对给定的0<α<1, 分数次幂算子Aα是有定义的. 分数次幂空间Xα=D(Aα)是一个自反的Banach空间,其上的范数定义为 ‖u‖α=‖Aαu‖X. 为简单起见,我们记X 上的范数为‖⋅‖.
(A2)~ 对给定的0<α<1,非线性项f:Xα×Xα→X具有以下性质: 对Xα×Xα中任何有界集,存在常数L使得下式成立
定义2.1 设I为某个区间,函数u称为分数次幂空间Xα中方程(1.1)的 (古典)解只要满足: u:I→Xα在区间I上连续可微, 其导数∂tu∈C(I,X),并且处处满足方程(1.1).
对于方程(1.1)的存在性结果,读者可参阅文献[18]. 通过常数变易公式, 方程(1.1)的解可表示为
定理2.1 假设A是Banach空间X上的正扇形算子,e−At是由−A生成 的解析半群,则成立下面的结论
(1)~ 对任意的α≥0,一定存在常数Cα>0 使得对 所有的t>0
(2)~ 对任意的0<α≤1,存在常数Cα>0 使得当t≥0时, 对x∈D(Aα)有
(3)~ 对任意的α≥0,存在常数Cα>0 使得当t>0时, 对x∈X有
下面的引理是著名的Aubin-Lions紧性引理,我们仅给出引理的叙述, 有兴趣的读者可以参考文献[].
引理2.1 设X0,X,X1是三个Banach空间, 满足X0⊆X⊆X1. 假设X0紧嵌入到X, X连续嵌入到X1,并且假定X0和X1是自反空间. 对1<p,q<∞,记 W={u|u∈Lp([0,T;X0),u′∈Lq([0,T];X1)}, 则W紧嵌入到Lp[0,T;X).
最后我们给出一些基本的概念. 如果Au0=f(u0),i.e., ∂tu0≡0,则称u0为发展方程(1.2)的平衡点. 设u是方程(1.1)有界解,其ω -极限集定义为
定义2.2 [22] 强连续Cr -半流S(t):X→X (r≥1)称为梯度系统, 如果存在一个连续的函数V:X→R满足以下性质
(1)~ V(x)是下方有界的;
(2)~ 当‖x‖→∞时,V(x)→∞;
(3)~ 对每一个x∈X,V(S(t)x)关于时间t单调不增;
(4)~ 如果对t∈R,S(t)x有定义,并且V(S(t)x)=V(x), 则x是平衡点.
上述函数通常被称为Lyapunov函数.
注2.1 若方程(1.2)的平衡点是孤立的,u为全有界解,则
并且V(E1)>V(E2),方程(1.2)的全解是指定义在(−∞,+∞)上的解.
在证明主要结果之前,先给出一个引理,用于检验Aubin-Lions紧性引理成 立的条件.
引理3.1 假设(A1)和(A2)成立,u是方程(1.1)的解,对所有的t∈[−τ,+∞), ‖u(t)‖α≤M,则对任意的b>a≥τ,q>1, 一定存在γ∈(0,α)使得 u′(t)∈Lq(a,b;Xγ).
证 为简单起见,我们记F(t)=f(u(t),u(t−τ)),于是常数变易公式(2.3) 便可改写为
定理3.1 假设(A1)和(A2)成立,方程(1.2)的有限个平衡点是孤立的, 对任意给定的R,ε>0,一定存在一个足够小的 τ>0使得对方程(1.1)的任何满足lim supt→+∞‖u(t)‖α≤R的解最终将进入并停留在某一个平衡点的 ε -邻域内.
证 为简单起见,我们只需验证满足‖u(t)‖α≤R的解成立即可. 设{E1,⋯,En}是方程(1.2)的平衡点集,并假设按以下顺序排列
第一步,我们首先证明对任意的δ>0,存在一个足够小的τ>0 使得对所有的u∈ˉBR
我们采用反证法来证(3.15)式,如若不然,对给定的δ, 总能找到一个递减序列τk→0,其对应的解为uk,满足 d(Ej,ω(uk))≥2δ. 上式对所有的1≤j≤n 和k∈N成立. 考虑到ω(u)的定义,对每一个k我们总可以选取tk>0 使得当t≥tk,1≤j≤n
假设˜uk(t)=uk(t+tk) (t≥0), 容易验证˜uk满足下面的方程 ∂t˜uk+A˜uk=f(˜uk(t),˜uk(t−τk)). 下面我们将利用Aubin-Lions引理来证明,对b>a≥τk和 1<p<∞,在Lp(a,b;Xα)中存在一个{˜uk}∞k=1的强收敛子列. 为了得到这个结论,我们首先检验,对1<p,q<∞,0<γ<α<β<1
和
第二步,我们将证明,如果τ足够小,则对方程(1.1)的任何一个有界解, 一定存在某个平衡点Ej和足够大的T使得当t>T时
因为vi∈C([−(ηi−ti),+∞);Xα), 当然vi(t)在t=0处连续,也就是说limt→0+‖vi(t)−vi(0)‖α=0. 由连续的定义得,选取ε0<ε和t1>0使得 ‖vi(t)−vi(0)‖α<ε0,对0<t<t1. 于是选取t2∈(0,t1)使得‖vi(t)−vi(0)‖α<ε0,对t2<t<t1. 因此,我们得到 ε0>‖vi(t)−Ej1+Ej1−vi(0)‖α≥‖vi(0)−Ej1‖α−‖vi(t)−Ej1‖α. 从而‖vi(t)−Ej1‖α>‖vi(0)−Ej1‖α−ε0=ε−ε0>0,对t2<t<t1. 于是 ∫t1t2‖vi(t)−Ej1‖pαdt>(t1−t2)(ε−ε0)p. 显然limi→∞∫t1t2‖vi(t)−Ej1‖pαdt≠0. 这与在Lp(a,b;Xα)中,vi(t)强收敛于Ej1 的事实相矛盾,因此T=+∞.
假设limt→+∞v(t)=Ej, 则一定有j≥j1. 否则,如果j<j1 ∫ba‖v(t)−Ej‖pαdt=∫ba‖v(t)−vi(t)+vi(t)−Ej‖pαdt≥C1∫ba‖vi(t)−Ej‖pαdt−C2∫ba‖vi(t)−v(t)‖pαdt>C1εp−C2∫ba‖vi(t)−v(t)‖pαdt>Cεp,(C>0). 这与vi(t)在Lp(a,b;Xα)中强收敛到v(t)相矛盾.
最后我们需要验证limt→−∞v(t)=Ej1. 对于足够小的ε′>0使得−(ηi−ti)<t<−(ηi−ti)+ε′<0,于是有 ∫−(ηi−ti)+ε′−(ηi−ti)‖v(t)−Ej1‖pαdt=∫−(ηi−ti)+ε′−(ηi−ti)‖v(t)−vi(t)+vi(t)−Ej1‖pαdt≤C1∫−(ηi−ti)+ε′−(ηi−ti)‖vi(t)−v(t)‖pαdt+C2∫−(ηi−ti)+ε′−(ηi−ti)‖vi(t)−Ej1‖pαdt≤C1∫−(ηi−ti)+ε′−(ηi−ti)‖vi(t)−v(t)‖pαdt+C2ε′ε. 考虑到在Lp(−(ηi−ti),−(ηi−ti)+ε′;Xα) 空间中,vi强收敛到v,于是可以得到上述结论成立. 总之,我们有 limt→−∞v(t)=Ej1,limt→+∞v(t)=Ej(j≥j1). 根据一开始的假设(3.14),这显然与注2.1的事实想违背, 从而(3.21)式是正确的.
根据jk1i的定义,我们可以得出结论 d(Ej,ω(uk1i))≥δ1,对所有的k1i>k∗1,1≤j≤j1. 为简单起见,不妨假设上式对所有的k1i都成立. 固定 δ′2∈(0,δ1),用jk1i表示最小的满足下式成立的j ω(uk1i)⋂ˉBδ′2(Ej)≠∅. 容易看出,对所有的k1i,jk1i>j1成立. 类似的, 存在{k1i}∞i=1的子列{k2i}∞i=1 和j2∈(j1,n]使得对所有的k2i有jk2i=j2成立. 重复前面的过程,能够证明,如果j2<n,存在一个 δ2∈(0,δ′2)和k∗2>k∗1使得 d(Ej2,ω(uk2i))≥δ2,对k2i>k∗2. 根据{k2i}∞i=1 的选取,容易看出 d(Ej,ω(uk2i))≥δ2,对所有的k2i>k∗2,1≤j≤j2. 如此反复上面的过程,我们最终会得到序列 j1<j2<⋯<jm=n,ε>δ1>δ2>⋯>δm>0,k∗1<k∗2<⋯<k∗m 和{kpi}∞i=1(1≤p≤m)使得 d(Ej,ω(ukpi))≥δp,对所有的kpi>k∗p,1≤j≤jp. 特别地,当然有 d(Ej,ω(ukmi))≥δm,对所有的kmi>k∗m,1≤j≤jm=n. 这显然与(3.20)式相矛盾,定理证毕.
下面一个定理证明了方程(1.1)的解将会精确地走向平衡点, 只要平衡点是双曲的.
定理3.2 假设(A1)和(A2)成立,f∈C1. E是方程(1.2)的双曲平衡点 (也就是说,Df(E)的特征值有非零实部). 则存在ε>0和 足够小的τ>0使得当t趋于无穷大时,方程 (1.1) 的任何一个满足 lim supt→∞‖u(t)−E‖α≤ε的解最后会走向平衡点E.
证 为简单起见,我们选取E=0,以下只需验证存在ε和足 够小的τ,使得方程(1.1)的有界解满足‖u(t)‖α≤ε. 考虑方程(1.2)的线性化方程 ∂tv+Av−Df(0)v=∂tv+Lv=0, 这里L=A−Df(0),Df(0)是f在0点的Fr\'{e}chet导数, 则−L在X上能生成一个解析半群e−Lt. 因为平衡点是双曲的, 因此解析半群e−Lt在X上具有指数二分性. 也就是说, 存在一个有界的投影算子P和两个正常数β0,α0 使得对所有的u∈X,下列不等式成立
我们首先证明,如果u是方程(1.1)在[0,+∞)上的有界解,则它能够表示为
结合定理3.1和3.2,我们最终得到下面的结论.
定理3.3 假设(A1),(A2)成立,非线性项f∈C1, 方程(1.2)存在有限个孤立双曲平衡点,则对任何的R>0, 方程(1.1)的满足lim supt→+∞‖u(t)‖α≤R的解,当t→+∞时最终会走向一个平衡点, 只要时滞τ足够小.
作为例子,我们考虑下面时滞的初边值问题
定理4.1假设非线性项f满足(A2)并且是C1的, 方程(4.3)的每一个平衡点是孤立的双曲平衡点,则对任何的R>0, 方程(4.2)的满足lim supt→+∞‖u(t)‖α≤R的解,当t→+∞时最终会走向一个平衡点, 只要时滞τ足够小.