Processing math: 95%

数学物理学报, 2024, 44(1): 133-139

带有时滞边界反馈的退化波动方程的镇定

白晋彦,1, 柴树根,2,*

1.晋中学院数学系 山西晋中 030619

2.山西大学数学科学学院 太原 030006

Stabilization of Degenerate Wave Equations with Delayed Boundary Feedback

Bai Jinyan,1, Chai Shugen,2,*

1. Department of Mathematics, Jinzhong University, Shanxi Jinzhong 030619

2. School of Mathematical Sciences, Shanxi University, Taiyuan 030006

通讯作者: 柴树根, E-mail:sgchai@sxu.edu.cn

收稿日期: 2023-01-25   修回日期: 2023-09-11  

基金资助: 国家自然科学基金(12271316)
国家自然科学青年基金(12001343)
山西省自然科学基金(201901D111042)
晋中学院技术创新团队建设计划(jzxyjscxtd202111)

Received: 2023-01-25   Revised: 2023-09-11  

Fund supported: NSFC(12271316)
NSFY(12001343)
Natural Science Foundation of Shanxi Province(201901D111042)
Technical Innovation Team of Jinzhong University(jzxyjscxtd202111)

作者简介 About authors

白晋彦,E-mail:bjy207@163.com;

摘要

该文在一维空间中研究了带有时滞边界反馈的退化波动方程的镇定问题. 首先运用半群理论证明了其解的适定性, 然后通过选用合适的乘子以及构造合适的 Lyapunov 泛函,证明了该退化波动方程解的指数稳定性结果.

关键词: 退化波方程; 时滞; 反馈镇定

Abstract

In this paper, we study the stabilization of degenerate wave equations with time-delay boundary feedback. Firstly, the well-posedness of the solution is proved by using semigroup theory. And then the exponential stability of the solution is proved by selecting suitable multipliers and constructing suitable Lyapunov functional.

Keywords: Degenerate wave equation; Time delay; Feedback stabilization

PDF (491KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

白晋彦, 柴树根. 带有时滞边界反馈的退化波动方程的镇定[J]. 数学物理学报, 2024, 44(1): 133-139

Bai Jinyan, Chai Shugen. Stabilization of Degenerate Wave Equations with Delayed Boundary Feedback[J]. Acta Mathematica Scientia, 2024, 44(1): 133-139

1 引言

近年来, 退化偏微分方程的能控性和反馈镇定受到了人们越来越多的关注, 并且出现在了很多科学领域, 比如影像学[1]、气候学[2]、人类遗传学[3]等. 各科学领域的发展对退化偏微分方程的数学问题提出了越来越大的挑战. 这些问题的一个共同特征是涉及到一个变系数的椭圆算子, 在所考虑的空间区域上失去了一致椭圆性, 但在空间区域的紧子集上是保持一致椭圆性的. 这样就出现退化发生在区域的边界或内部小子集上的情形. 对退化波动方程而言, 一致椭圆性的缺失不仅影响解的适定性所属的泛函空间, 更对退化椭圆算子的各种估计提出了新的挑战, 因而对退化发生在边界和内部时退化波动方程的能控性和反馈镇定就成为了人们研究的热点问题[4-9].

众所周知, 时滞效应在日常生活的实际问题中也经常出现, 这些延迟效果的影响可将本来稳定的系统变得不稳定[10-13]. 从而要使具有时滞的系统稳定, 克服由时滞带来的负面影响, 控制的设计较无时滞时要困难很多. 具有时滞边界和内部反馈的非退化波动方程解的稳定性问题, 目前已经有不少相关的结论[14-19]. 特别地, 文献[14]指出, 当边界输入控制同时包括有时滞和无时滞两部分时, 非退化波动方程解的稳定性取决于这两部分的权重, 当无时滞部分的权重大于有时滞部分的权重时, 系统才能稳定.

具有边界控制的退化波动方程的一般形式为

{utt(xαux)x=0,(x,t)(0,1)×(0,+),u(0,t)=0,t(0,+),ux(1,t)=y(t),t(0,+),u(x,0)=u0(x),  ut(x,0)=u1(x),x(0,1),w(t)=ut(1,t),t(0,+),
(1.1)

其中α>0是退化程度的度量,y(t)是一个输入控制,w(t)是观测器. 如果输入控制y(t)没有时滞, 运用最简单的反馈控制y(t)=w(t), 即可使得闭环系统耗散且指数稳定 (见文献[4]). 那么, 当输入控制有了时滞, 上述闭环系统是否稳定? 本文将研究下面具有时滞边界反馈的退化波方程解的镇定问题

{utt(xαux)x=0,(x,t)(0,1)×(0,+),u(0,t)=0,t(0,+),ux(1,t)+μ1ut(1,t)+μ2ut(1,tτ)=0,t(0,+),u(x,0)=u0(x),  ut(x,0)=u1(x),x(0,1),ut(1,tτ)=f0(tτ),t(0,τ),
(1.2)

其中0α<1,u是状态变量,(u0,u1,f0)为任意属于合适空间的初值, 常数τ>0是时滞,μ1,μ2为正实数. 本文将证明μ1>μ2时, 上述系统的解是指数稳定的.

定义系统(1.2)的能量泛函为

E(t)=1210(u2t+xαu2x)dx+ξ210u2t(1,tτρ)dρ,t>0,
(1.3)

其中ξ是正实数, 满足τμ2<ξ<τ(2μ1u2). 注意到本文考虑μ1>μ2时的情形, 故满足条件的ξ是存在的.

本文的主要结果如下.

定理 1.1 假设μ1>μ2, 对任意的0α<1, 都存在正常数C1,C2, 使得

E(t)C1eC2tE(0),
(1.4)

即系统(1.2)的解是指数稳定的.

本文将在第 2 部分证明系统 (1.2) 解的适定性, 第 3 部分证明系统 (1.2) 解的指数稳定性 (定理 1.1).

2 退化波动方程解的适定性

为了研究系统 (1.2) 解的适定性, 首先引入一些加权的 Sobolev 空间.

W1α(0,1)={uL2(0,1)|u (0,1] 上是局部绝对连续的, xα2uxL2(0,1)  u(0)=0},
W2α(0,1)={uW1α(0,1)|xαuxH1(0,1)},

其中H1(0,1)是经典的 Sobolev 空间.

z(ρ,t)=ut(1,tτρ),ρ(0,1), t>0,则系统 (1.2) 可转化为

{utt(xαux)x=0,(x,t)(0,1)×(0,+),zρ(ρ,t)+τzt(ρ,t)=0,(ρ,t)(0,1)×(0,+),u(0,t)=0,t(0,+),ux(1,t)+μ1ut(1,t)+μ2z(1,t)=0,t(0,+),z(0,t)=ut(1,t),t(0,+),u(x,0)=u0(x),  ut(x,0)=u1(x),x(0,1),z(ρ,0)=f0(τρ),ρ(0,1).
(2.1)

U:=(u,ut,z), 则

U:=(ut,utt,zt)=(ut,(xαux)x,τ1zρ).

因此, 系统 (2.1) 可写为抽象发展方程

{U=AU,U(0)=(u0,u1,f0(τ)),
(2.2)

其中算子A被定义为

A(u,v,z)=(v,(xαux)x,τ1zρ),

D(A):={(u,v,z)W2α(0,1)×W1α(0,1)×H1ρ(0,1),ux(1)+μ1v(1)+μ2z(1)=0,z(0)=v(1)},

这里的H1ρ(0,1)以及后边的L2ρ(0,1)均是关于变量ρ的函数所属的经典的 Sobolev 空间.考虑具有如下内积的 Hilbert 空间H=W1α(0,1)×L2(0,1)×L2ρ(0,1)

(u,v,z),(˜u,˜v,˜z)H=10(xαux˜ux+v˜v)dx+ξ10z(ρ)˜z(ρ)dρ,(u,v,z),(˜u,˜v,˜z)H.

可以验证, 算子AH上是极大耗散的. 从而利用半群理论, 可以得到系统 (2.1) 解的适定性结果. 具体如下.

命题 2.1 算子AH上是极大耗散的.

U=(u,v,z)D(A), 则

(AU,U)=(v,(xαux)x,τ1zρ),(u,v,z)H=10[xαvxux+(xαux)xv]dxξτ10zρ(ρ)z(ρ)dρ=xαuxv|x=1x=0ξ2τz2(1)+ξ2τz2(0)=ux(1)v(1)ξ2τz2(1)+ξ2τz2(0)=μ1v2(1)μ2z(1)v(1)ξ2τz2(1)+ξ2τz2(0)(μ1+μ22+ξ2τ)v2(1)+(μ22ξ2τ)z2(1)0.
(2.3)

因此算子AH上是耗散的, 下面证算子AH上是极大的, 即对任意给定的λ>0,(λIA)是满射. 接下来, 对任给的(f,g,h)H, 寻找下面问题的一个解U=(u,v,z)D(A)

{λuv=f,λv(xαux)x=g,λz+1τzρ=h.
(2.4)

首先可以求得z. 由问题(2.4)的第三个方程可知,zρ+τλz=hτ,结合关系z(0)=v(1)以及 (2.4) 式的第一个方程, 有

z(ρ)=v(1)eλρτ+eλρτρ0hτeλστdσ=λu(1)eλρτf(1)eλρτ+eλρτρ0hτeλστdσ.
(2.5)

特别地,

z(1)=λu(1)eλτf(1)eλτ+eλτ10hτeλστdσ.
(2.6)

另外, 结合(2.4)式的第一个和第二个方程, 可知需要寻找uW2α(0,1)且满足

λ2u(xαux)x=λf+g.
(2.7)

现在假设uW2α(0,1)存在, 且是(2.7)式的解, 则有

10(λ2u(xαux)x)wdx=10(λf+g)wdx,  wW1α(0,1).
(2.8)

分部积分后整理可得

10(λ2uw+xαuxwx)dx+(μ1+μ2eλτ)λu(1)w(1)=10(λf+g)wdx+μ1f(1)w(1)+μ2eλτ(f(1)10hτeλστdσ)w(1).
(2.9)

Λ(u,w)=10(λ2uw+xαuxwx)dx+(μ1+μ2eλτ)λu(1)w(1),
L(w)=10(λf+g)wdx+μ1f(1)w(1)+μ2eλτ(f(1)10hτeλστdσ)w(1).

由于ΛW1α(0,1)×W1α(0,1)上是连续双线性泛函且满足强制性,LW1α(0,1)上是连续的, 从而根据 Lax-Milgram 定理可知, (2.9) 式存在唯一解uW1α(0,1). 下面将证(u,v,z)D(A)且是 (2.4) 式的解.

D(0,1)(0,1)上具有紧支集的C函数构成的集合, 由于D(0,1)W1α(0,1), 所以

10(λ2u(xαux)x)wdx=10(λf+g)wdx,  wD(0,1).

根据对偶可知在分布意义下u是 (2.7) 式的解, 因此uW2α(0,1)且在(0,1)上几乎处处有

λ2u(xαux)x=λf+g.

将上面方程在(0,1)上积分并利用分部积分公式, 结合(2.9)式有

ux(1)+(μ1+μ2eλτ)λu(1)=μ1f(1)+μ2eλτ(f(1)10hτeλστdσ),

运用(2.4)和(2.6)式可得ux(1)+μ1v(1)+μ2z(1)=0.综上所知,(u,v,z)D(A)且是问题(2.4)的解. 因此有如下结果. 证毕.

定理 2.1 对任意初值U(0)H, 系统(2.2)存在唯一解UC([0,+),H). 进一步, 若U(0)D(A), 那么UC([0,+),D(A))C1([0,+),H).

3 退化波动方程解的指数稳定性

命题 3.1 对系统 (1.2) 的任意解, 对应解的能量是递减的, 且存在正常数C使得下面的估计式成立

E(t)C(u2t(1,t)+u2t(1,tτ)),t>0.
(3.1)

由于

E(t)=1210(u2t(x,t)+xαu2x(x,t))dx+ξ210u2t(1,tτρ)dρ=1210(u2t(x,t)+xαu2x(x,t))dx+ξ2τttτu2s(1,s)ds,
(3.2)

从而

E(t)=10(ututt+xαuxuxt)dx+ξ2τu2t(1,t)ξ2τu2t(1,tτ).

对上述等式进行分部积分并利用系统 (1.2) 的边界条件以及 Hölder 不等式, 可得

E(t)=μ1u2t(1,t)μ2ut(1,t)ut(1,tτ)+ξ2τu2t(1,t)ξ2τu2t(1,tτ)(μ1+μ22+ξ2τ)u2t(1,t)+(μ22ξ2τ)u2t(1,tτ).

最后根据假设τμ2<ξ<τ(2μ1u2)可知, 存在常数C>0, 使得估计式(3.1)成立. 证毕.

命题 3.20α<1时, 对系统 (1.2) 的任意解有下面的不等式成立

ddt102xuxutdx(1+2μ21)u2t(1,t)+2μ22u2t(1,tτ)10u2tdx(1α)10xαu2xdx.
(3.3)

根据分部积分公式和 Cauchy-Schwarz 不等式, 有

ddt102xuxutdx=10(2xuxtut+2xux(xαux)x)dx=u2t(1,t)10u2tdx+2u2x(1,t)210xαu2xdx10xα+1(u2x)xdx=u2t(1,t)10u2tdx+u2x(1,t)(1α)10xαu2xdx(1+2μ21)u2t(1,t)+2μ22u2t(1,tτ)10u2tdx(1α)10xαu2xdx.

得证.

现在, 引入函数

δ(t)=1τttτestu2t(1,s)ds,

可以很容易估计到

δ(t)=1τ(u2t(1,t)eτu2t(1,tτ))1τttτestu2t(1,s)ds1τ[u2t(1,t)eτu2t(1,tτ)eτttτu2t(1,s)ds]=1τ[u2t(1,t)eτu2t(1,tτ)]eτ10u2t(1,tτρ)dρ.
(3.4)

构造 Lyapunov 泛函

ε(t):=E(t)+γ1102xuxutdx+γ2δ(t),
(3.5)

其中γ1,γ2是足够小的正常数, 将在后边给出.

由于

|γ1102xuxutdx+γ2δ(t)|γ110(xαu2x+u2t)dx+γ2τttτu2t(1,s)dsmax{2γ1,2γ2ξ}E(t),
(3.6)

所以通过选择足够小的γ1γ2, 使得1max{2γ1,2γ2ξ}=a1>0, 且记1+max{2γ1,2γ2ξ}=a2, 则有

a1E(t)ε(t)a2E(t).
(3.7)

ε(t)和能量E(t)是等价的.

定理 3.10α<1, 则存在正常数c1,c2使得系统 (1.2) 的任意解有

ε(t)c1ec2tε(0),t>0.
(3.8)

对已构造的 Lyapunov 泛函ε(t)关于t求导, 并结合命题 3.1, 命题 3.2 可得

ε(t)=E(t)+γ1ddt102xuxutdx+γ2δ(t)[C+γ1(1+2μ21)+γ2τ]u2t(1,t)+[C+2γ1μ22γ2τeτ]u2t(1,tτ) min{2γ1(1α),2γ2eτξ}E(t).
(3.9)

对固定的C>0, 通过选择足够小的γ1,γ2, 使得

C+γ1(1+2μ21)+γ2τ0,C+2γ1μ22γ2τeτ0.

显然满足条件的γ1,γ2是存在的, 从而

ε(t)cE(t),
(3.10)

其中c>0, 结合(3.7)式的第二个不等式, 即可证得(3.8) 式. 证毕.

由定理3.1和能量的等价性(3.7) 式, 即可证得定理1.1.

参考文献

Citti G, Manfredini M.

A degenerate parabolic equation arising in image processing

Commun Appl Anal, 2004, 8(1): 125-141

[本文引用: 1]

Ghil M.

Climate stability for a Sellers type mode

J Atmo Sci, 1976, 29(5): 483-493

[本文引用: 1]

Ethier S N, Kurtz T G.

Fleming-Viot processes in population genetics

SIAM J Control Optim, 1993, 31(2): 345-386

DOI:10.1137/0331019      URL     [本文引用: 1]

Alabau-Boussouira F, Cannarsa P, Leugering G.

Control and stabilization of degenerate wave equation

SIAM J Control Optim, 2017, 55(3): 2052-2087

DOI:10.1137/15M1020538      URL     [本文引用: 2]

Zhang M M, Gao H.

Null controllability of some degenerate wave equations

J Syst Sci Complex, 2017, 30(5): 1027-1041

DOI:10.1007/s11424-016-5281-3      [本文引用: 1]

This paper is devoted to a study of the null controllability problems for one-dimensional linear degenerate wave equations through a boundary controller. First, the well-posedness of linear degenerate wave equations is discussed. Then the null controllability of some degenerate wave equations is established, when a control acts on the non-degenerate boundary. Different from the known controllability results in the case that a control acts on the degenerate boundary, any initial value in state space is controllable in this case. Also, an explicit expression for the controllability time is given. Furthermore, a counterexample on the controllability is given for some other degenerate wave equations.

Zhang M M, Gao H.

Persistent regional null controllability of some degenerate wave equations

Math Meth Appl Sci, 2017, 40(16): 5821-5830

DOI:10.1002/mma.v40.16      URL     [本文引用: 1]

Zhang M M, Gao H.

Interior controllability of semi-linear degenerate wave equations

J Math Anal Appl, 2018, 457(1): 10-22

DOI:10.1016/j.jmaa.2017.07.057      URL     [本文引用: 1]

Bai J Y, Chai S G.

Exact controllability for some degenerate wave equations

Math Meth Appl Sci, 2020, 43(12): 7292-7302

DOI:10.1002/mma.v43.12      URL     [本文引用: 1]

Bai J Y, Chai S G.

Exact controllability for a one-dimensional degenerate wave equation in domains with moving boundary

Appl Math Lett, 2021, 119: Article ID 107235

[本文引用: 1]

Gumowski I, Mira C. Optimization in Control Theory and Practice. Cambridge: Cambridge University Press, 1968

[本文引用: 1]

Datko R, Lagness J, Poilis M P.

An example on the effect of time delays in boundary feedback stabilization of wave equations

SIAM J Control Optim, 1986, 24(1): 152-156

DOI:10.1137/0324007      URL     [本文引用: 1]

Datko R.

Not all feedback stabilized hyperbolic systems are robust with respect to small time delay in their feedbacks

SIAM J Control Optim, 1988, 26(3): 697-713

DOI:10.1137/0326040      URL     [本文引用: 1]

Wang J N, Yang D Z.

Stability and bifurcation of a pathogen-immune model with delay and diffusion effects

Acta Math Sci, 2021, 41A(4): 1204-1217

[本文引用: 1]

Xu G Q, Yung S P, Li L K.

Stabilization of wave systems with input delay in the boundary control

ESAIM Control Optim Calc Var, 2006, 12(4): 770-785

DOI:10.1051/cocv:2006021      URL     [本文引用: 2]

Nicaise S, Pignotti C.

Stability and instability results of the wave equation with a delay term in the boundary or internal feedbacks

SIAM J Control Optim, 2006, 45(5): 1561-1585

DOI:10.1137/060648891      URL     [本文引用: 1]

Nicaise S, Pignotti C.

Stabilization of the wave equation with boundary or internal distributed delay

Differential Integral Equ, 2008, 21(9/10): 935-958

[本文引用: 1]

Ait Benhassi E M, Ammari K, Boulite S, Maniar L.

Feedback stabilization of a class of evolution equations with delay

J Evol Equ, 2009, 9(1): 103-121

DOI:10.1007/s00028-009-0004-z      URL     [本文引用: 1]

Nicaise S, Valein J.

Stabilization of second order evolution equations with unbounded feedback with delay

ESAIM Control Optim Calc Var, 2010, 16(2): 420-456

DOI:10.1051/cocv/2009007      URL     [本文引用: 1]

Ammari K, Nicaise S, Pignotti C.

Feedback boundary stabilization of wave equations with interior delay

Syst Control Lett, 2010, 59(10): 623-628

DOI:10.1016/j.sysconle.2010.07.007      URL     [本文引用: 1]

/