Processing math: 43%

数学物理学报, 2023, 43(5): 1471-1482

一类 2×2 双曲偏微分系统的 PDP 边界控制

庞玉婷,1, 赵东霞,1,*, 赵鑫,2, 高彩霞,1

1中北大学数学学院 太原 030051

2太原理工大学电气与动力工程学院 太原 030024

The PDP Boundary Control for a Class of 2×2 Hyperbolic Partial Differential System

Pang Yuting,1, Zhao Dongxia,1,*, Zhao Xin,2, Gao Caixia,1

1School of Mathematics, North University of China, Taiyuan 030051

2School of Electrical and Power Engineering, Taiyuan University of Technology, Taiyuan 030024

通讯作者: * 赵东霞,Email: zhaodongxia6@sina.com

收稿日期: 2022-03-18   修回日期: 2023-03-8  

基金资助: 山西省基础研究计划资助项目(20210302123046)

Received: 2022-03-18   Revised: 2023-03-8  

Fund supported: Fundamental Research Program of Shanxi Province(20210302123046)

作者简介 About authors

庞玉婷,Email:2116786325@qq.com;

赵鑫,Email:1808642517@qq.com;

高彩霞,Email:1519546532@qq.com

摘要

该文研究了具有恒定坡度和底部摩擦的单段明渠系统的指数稳定性, 该系统由2×2双曲偏微分方程描述. 通过设计位置反馈和延迟位置反馈(简称PDP)边界控制器来解决系统的反馈镇定问题. 首先, 利用算子半群理论给出系统适定性的证明. 然后, 通过构造合适的Lyapunov函数分析闭环系统的指数稳定性, 并获得了反馈参数和时滞需要满足的充分条件. 此外, 还利用谱分析方法建立了系统算子特征根和特征向量的渐近表达式. 最后, 通过一个数值仿真示例来评估PDP控制器性能.

关键词: 双曲PDE系统; Lyapunov函数; PDP边界控制器; 指数稳定性

Abstract

This paper studies the exponential stability of a single open-channel system with constant slope and bottom friction, which is described by a 2×2 hyperbolic partial differential equation. The position feedback and delayed position feedback (PDP for short) boundary controller is designed to solve the problem of feedback stabilization. Firstly, the well-posedness of the system is proved by using operator semigroup theory. Then, the exponential stability of the closed-loop system is analyzed by constructing an appropriate Lyapunov function, and sufficient conditions for feedback parameters and time-delay are obtained. In addition, the asymptotic expressions of the eigenvalues and the eigenfunctions of the system operator are given by spectral analysis method. Finally, a numerical example is used to evaluate the performance of the PDP controller.

Keywords: Hyperbolic PDE system; Lyapunov function; PDP boundary controller; Exponential stability

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

本文引用格式

庞玉婷, 赵东霞, 赵鑫, 高彩霞. 一类 2×2 双曲偏微分系统的 PDP 边界控制[J]. 数学物理学报, 2023, 43(5): 1471-1482

Pang Yuting, Zhao Dongxia, Zhao Xin, Gao Caixia. The PDP Boundary Control for a Class of 2×2 Hyperbolic Partial Differential System[J]. Acta Mathematica Scientia, 2023, 43(5): 1471-1482

1 引言

1871年, 法国科学家Saint-Venant建立了著名的非线性Saint-Venant偏微分方程, 该方程由反映质量守恒的连续方程和反映动量守恒的运动方程组成, 可以描述渠道和其他具有自由表面的浅水体中非恒定渐变水流的运动规律. 到目前为止, 已有大量关于该方程的研究成果, 被广泛应用于明渠灌溉[1], 运河航运[2], 管道运输[3]和交通流系统[4].

Saint-Venant方程属于一阶拟线性2×2双曲偏微分方程. 在过去的几十年里, 许多形式的2×2双曲型偏微分系统得到了广泛的研究, 并取得了一系列兼具理论意义和应用价值的成果. 文献[5]致力于有界闭区间上一维双曲守恒律系统和平衡律系统的控制器设计和指数稳定性的证明. 在文献[6]中, 针对一类具有变系数的一维2×2线性双曲型偏微分方程组, 基于反步法研究了边界稳定和状态估计问题. 在文献[7]中, 利用频域法分析双曲守恒律系统在比例反馈下的稳定性问题, 并利用Nyquist定理给出了闭环系统稳定性的充要条件. 为了控制双曲守恒律系统达到指定状态, 文献[8]提出一种基于Riemann不变量的边界控制方法. 文献[9]则是利用Lyapunov函数研究了该系统在比例积分(PI)边界条件下的反馈镇定和扰动抑制问题.

对于双曲系统的边界控制问题, 关键是控制器设计和稳定性分析. 从已有文献来看, 稳定性分析方法通常有Lyapunov函数法, 频域法, Riemann不变量法, 反步法以及谱分析方法等[10,11,12,13,14]. 在控制工程实践中, PI控制器能够改善系统的稳态性能, 受到广大学者的广泛关注. 在文献[15][16]中, 相继研究了星形网络结构和级联网络结构的双曲偏微分系统的PI边界控制器的设计. 此外, 控制器的更详细设计可参见文献[17](第5章, 第7章). 在文献[18,19]中, 对于一类双曲守恒律偏微分系统, 提出了一种前馈边界控制方法来抑制量测扰动. 当受控系统受到显著的输入扰动时, 这种控制技术可以在影响系统输出之前进行测量和补偿.

众所周知, 在实际应用中, 系统控制器和传感器之间的通信不可避免地会产生时间延迟, 这有时可以提高受控系统的性能[20], 亦有可能导致控制性能下降, 甚至使原本稳定的系统变得不稳定[21,22]. 在文献[23,24]中,Suh和Bien首先设计了比例负延迟(PMD, Proportional Minus Delay)控制器, 并证明在传统控制器上增加时间延迟反而可以改善系统的性能. 然而, 文献[23,24]中给出的结果是基于近似分析和数值模拟的. 之后, Atay[25]将PMD控制器应用于单摆和倒立摆系统, 采用特征根分析方法给出了系统参数的稳定域, 并得出无论时滞取值如何, 合适选取反馈参数均可使系统达到渐近稳定. Hu[26]使用延迟位置反馈或延迟速度反馈或两者结合来镇定线性无阻尼振荡器的周期振动. 随后, Liu和Hu[27]提出了一种结合位置反馈和延迟位置反馈(PDP, Position feedback and Delayed Position feedback)的控制器, 以稳定一类多自由度线性无阻尼系统. 近年来, 时滞反馈亦被用于镇定PDE系统. 在文献[28]中, 通过使用反步法和Lyapunov函数法, 研究了输入端具有任意长延迟的反应扩散偏微分方程的指数稳定性. 在文献[29]中, 基于简化的Hayami扩散波模型, 建立了系统的PDP 反馈控制策略, 并利用Lyapunov函数法建立了闭环系统的指数稳定性. 然而, 据我们所知, 对于2×2双曲型PDE 系统, 目前还没有关于时滞反馈控制器的研究结果.

本文致力于设计PDP边界控制器来镇定双曲型PDE系统, 该系统描述的是具有底部坡度和底部摩擦的渠道水流动力学. 本文的贡献是在输入端设计PDP 控制器, 并利用算子半群理论证明了系统的适定性. 然后通过构造合适的Lyapunov函数, 建立反馈参数和时滞值所满足的充分条件, 以保证闭环系统的指数稳定性. 本文的结构如下: 在第2节中, 利用一阶运输方程初值问题的解来描述PDP控制器中的时滞项, 得到了PDE-PDE耦合闭环系统, 并将其改写为抽象发展方程的形式. 第3节利用算子半群理论证明了系统的适定性. 在第4节中, 结合Lyapunov函数方法给出控制参数和时滞值的充分条件, 证明了闭环系统在L2范数下的指数稳定性. 此外, 在第5节中, 采用谱分析方法给出了系统特征值和特征向量的渐近表达式. 最后, 在第6节中采用数值仿真来评估PDP控制器性能.

2 模型建立和PDP控制器设计

本文所进行的研究是基于文献[30]中的一个经典模型, 该模型描述了一个具有矩形横截面和恒定底部坡度的棱柱形渠道池, 如图1所示. 具体地, 该双曲系统的动力学方程如下:

t(HV)+x(HVV22+gH)+(0g[Sf(H,V)Sb])=0,
(2.1)

其中, H(t,x)V(t,x)分别表示t时刻x位置处的水面高度和水流速度, t>0,x(0,L), L是渠道长度, g为重力加速度, Sb为底部坡度, Sf=CfV2H表示底部摩擦, 且Cf为常数, 代表摩擦系数.

图1

图1   具有底部坡度的单渠道系统侧视图


假设系统(2.1)的稳态为(H,V)且满足Sf(H,V)=Sb. 在平衡状态附近将系统(2.1)作线性化处理

{th+Vxh+Hxv=0,tv+gxh+Vxv[gHSb]h+[2gVSb]v=0,
(2.2)

其中h(t,x)=H(t,x)H,v(t,x)=V(t,x)V. 作Riemann坐标变换

{ξ1(t,x)v(t,x)+h(t,x)gH,ξ2(t,x)v(t,x)h(t,x)gH,
(2.3)

可得线性化Saint-Venant方程(2.2)的特征形式(详见文献[5]第一章或文献[30]):

{tξ1(t,x)+λ1xξ1(t,x)+γξ1(t,x)+δξ2(t,x)=0,tξ2(t,x)λ2xξ2(t,x)+γξ1(t,x)+δξ2(t,x)=0,
(2.4)

其中

λ1=gH+V,λ2=gHV,γ=gSb(1V12gH),δ=gSb(1V+12gH).
(2.5)

本文仅考虑亚临界流动状态的情形, 即gH(V)2>0. 因此, λ1>λ2>00<γ<δ.

不妨设系统(2.4)的初始条件为

(ξ1(0,x)ξ01(x),ξ2(0,x)ξ02(x))L2((0,L);R2),
(2.6)

边界条件及控制器U(t)的设置为

ξ1(t,0)=R1ξ2(t,0)+U(t),ξ2(t,L)=R2ξ1(t,L),
(2.7)

其中, 参数R1,R2为实常数. 在控制器的设计中, 充分考虑到传输过程中不可避免地存在时间延迟, 即

U(t)=Kp(H(t,0)H)+Kd(H(tτ,0)H)=Kph(t,0)+Kdh(tτ,0)=Kp2Hg(ξ1(t,0)ξ2(t,0))+Kd2Hg(ξ1(tτ,0)ξ2(tτ,0)),
(2.8)

利用位置反馈和时滞位置反馈的线性组合设计控制器以保证闭环系统(2.4)-(2.8)的指数稳定性, 其中, 反馈参数Kp,KdR, 时滞τ>0.

考虑到一阶运输方程初值问题

{tu(t,x)+Lτxu(t,x)=0,u(0,x)=ξ1(τxL,0)
(2.9)

的解可以表示为

u(t,x)=ξ1(tτxL,0).
(2.10)

同理, 令

v(t,x)=ξ2(tτxL,0),
(2.11)

则闭环系统(2.4)-(2.8)可写为如下PDE-PDE耦合系统的形式

{tξ1(t,x)+λ1xξ1(t,x)+γξ1(t,x)+δξ2(t,x)=0,tξ2(t,x)λ2xξ2(t,x)+γξ1(t,x)+δξ2(t,x)=0,tu(t,x)+Lτxu(t,x)=0,tv(t,x)+Lτxv(t,x)=0,ξ1(t,0)=k0ξ2(t,0)+k1(u(t,L)v(t,L)),ξ2(t,L)=k2ξ1(t,L),u(t,0)=ξ1(t,0),v(t,0)=ξ2(t,0),
(2.12)

其中

k0=2R1gKpH2gKpH,k1=KdH2gKpH,k2=R2.
(2.13)

在Hilbert状态空间

H=L2(0,L)×L2(0,L)×L2(0,L)×L2(0,L)

中, 定义内积

X1,X2=L0[f1(x)¯f2(x)+g1(x)¯g2(x)]dx+τL0[h1(x)¯h2(x)+q1(x)¯q2(x)]dx,
(2.14)

其中Xi=(fi,gi,hi,qi)H (i=1,2). 定义线性算子A:D(A)HH

AX=(λ1xγδ00γλ2xδ0000Lτx0000Lτx)(fghq),
(2.15)
D(A)={(f,g,h,q)(H1(0,L))4f(0)=k0g(0)+k1(h(L)q(L)),g(L)=k2f(L),h(0)=f(0),q(0)=g(0)}.
(2.16)

那么, 在Hilbert状态空间H上, 系统(2.12)可改写为如下抽象发展方程的形式

{˙X(t)=AX(t),t>0,X(0)=X0,
(2.17)

其中X(t)=(ξ1(,t),ξ2(,t),u(,t),v(,t)).

3 系统(2.17)的适定性

为了后续计算方便, 不失一般性, 在本节令L=1.

定理3.1 设算子A由(2.15)和(2.16)式给出. 如果参数满足k0k1,k11,k2λ2+λ10, 那么A1存在并且是紧的. 从而A的谱集σ(A)仅由有穷代数重数的孤立特征值组成.

对于 X1=(f1,g1,h1,q1)H, 求解AX=X1, 其中 X=(f,g,h,q)D(A), 可得

A(fghq)=(λ1fγfδgλ2gγfδg1τh1τq)=(f1g1h1q1),
(3.1)

{λ1f+γf+δg+f1=0,λ2gγfδgg1=0,1τh=h1,1τq=q1,f(0)=k0g(0)+k1(h(1)q(1)),g(1)=k2f(1),h(0)=f(0),q(0)=g(0).
(3.2)

由(3.2)式的前两个方程可得

λ1f(x)+λ2g(x)λ1f(0)λ2g(0)=x0(g1(s)f1(s))ds.
(3.3)

进而

g(x)=λ1λ2f(x)+λ1(k0k1)+λ2(1k1)λ2(1k1)g(0)λ1k1τλ2(1k1)10(h1(s)q1(s))ds+1λ2x0(g1(s)f1(s))ds.
(3.4)

再将(3.4)式代入(3.2)式的第一个方程可得

f(x)=(δλ2γλ1)f(x)λ1(k0k1)δ+λ2(1k1)δλ1λ2(1k1)g(0)f1(x)λ1+δk1τλ2(1k1)10(h1(s)q1(s))dsδλ1λ2x0(g1(s)f1(s))ds.
(3.5)

因此, 由一阶线性常微分方程的通解公式得

f(x)= cec1xx0ec1(xν)[λ1(k0k1)δ+λ2(1k1)δλ1λ2(1k1)g(0)δk1τλ2(1k1)10(h1(s)q1(s))ds+δλ1λ2ν0(g1(s)f1(s))ds+f1(ν)λ1]dν,
(3.6)

其中c=f(0), c1=δλ2γλ10. 结合(3.2)式的第六个方程和(3.4)式得

k2f(1)=λ1λ2f(1)+[λ1(k0k1)λ2(1k1)+1]g(0)λ1k1τλ2(1k1)10(h1(s)q1(s))ds+1λ210(g1(s)f1(s))ds.
(3.7)

另外, 由(3.6)式可得

f(1)=f(0)ec110ec1(1ν)[λ1(k0k1)δ+λ2(1k1)δλ1λ2(1k1)g(0)δk1τλ2(1k1)10(h1(s)q1(s))ds+δλ1λ2ν0(g1(s)f1(s))ds+1λ1f1(ν)]dν.
(3.8)

故结合(3.7)和(3.8)式, 有

[ec1+δ(1ec1)c1λ1(k0k1)+λ2(1k1)λ1λ2(k0k1)λ1(k0k1)+λ2(1k1)(k2λ2+λ1)(k0k1)]f(0)=τk11k1[δ(ec11)c1λ2(λ1(k0k1)+λ2(1k1)(k0k1)λ11)+λ1(k0k1)+λ2(1k1)(k0k1)(k2λ2+λ1)λ1k2λ2+λ1]10(h1(s)q1(s))dsδλ1λ2c110(1ec1(1s))(g1(s)f1(s))ds+1k2λ2+λ110(g1(s)f1(s))ds+1λ110ec1(1ν)f1(ν)dν.
(3.9)

不妨记

M1=τk11k1[δ(ec11)c1λ2(λ1(k0k1)+λ2(1k1)(k0k1)λ11)+λ1(k0k1)+λ2(1k1)(k0k1)(k2λ2+λ1)λ1k2λ2+λ1]10(h1(s)q1(s))dsδλ1λ2c110(1ec1(1s))(g1(s)f1(s))ds+1k2λ2+λ110(g1(s)f1(s))ds+1λ110ec1(1ν)f1(ν)dν,M2=ec1+δ(1ec1)c1λ1(k0k1)+λ2(1k1)λ1λ2(k0k1)λ1(k0k1)+λ2(1k1)(k2λ2+λ1)(k0k1),
(3.10)

于是当M20时可解得f(0)=M1M2. 进而, 结合(3.2)式的第三个方程和(3.6)式可确定f(x), 再由(3.4)式可得g(x), 具体如下

{f(x)=M1M2ec1xx0ec1(xν)[λ1(k0k1)+λ2(1k1)λ1λ2(k0k1)δM1M2+δk1τλ1(k0k1)10(h1(s)q1(s))ds+1λ1f1(ν)]dνδλ1λ2c1x0(ec1(xν)1)(g1(ν)f1(ν))dν,g(x)=λ1λ2f(x)+λ1(k0k1)+λ2(1k1)λ2(k0k1)M1M2+1λ2x0(g1(s)f1(s))ds+k1τk0k1x0(h1(s)q1(s))ds,h(x)=M1M2τx0h1(s)ds,q(x)=(1k1)M1(k0k1)M2+k1τk0k1x0(h1(s)q1(s))dsτx0q1(s)ds.
(3.11)

因此, 由Sobolve嵌入定理可知, A1存在并且在H上是紧的. 从而, σ(A)仅由有穷代数重数的孤立特征值构成. 证毕.

4 系统(2.12)的指数稳定性

本节致力于构建一个合适的Lyapunov函数来分析系统(2.12)的指数稳定性.

定义4.1 如果存在正常数ϑC, 使得柯西问题(2.12)的解对任意初始条件(ξ01(x),ξ02(x))L2((0,L);R2) 都满足

则闭环系统(2.12)在L^2范数意义下是指数稳定的.

定理4.1 假定参数k_i\neq0\ (i=0,1,2). 如果存在正常数p_i\ (i=1,2,3,4)和充分小的\mu>0使得系统参数以及调控参数k_i和时滞\tau满足不等式

\begin{equation} \begin{cases} \lambda_2p_2k_2^2{\rm e}^{\mu L}\leq\lambda_1p_1{\rm e}^{-\mu L},\; (\frac{1}{\tau}p_3L+\lambda_1p_1)(k_0^2+2\left|k_0k_1\right|)+\frac{1}{\tau}p_4L\leq\lambda_2p_2,\\ (\frac{1}{\tau}p_3L+\lambda_1p_1)(2k_1^2+\left|k_0k_1\right|)\leq\frac{1}{\tau}L{\rm e}^{-\mu L}\min\{p_3,p_4\},\; \delta p_1=\gamma p_2, \end{cases} \end{equation}
(4.1)

则闭环系统(2.12)是指数稳定的.

构造Lyapunov候选函数

\begin{equation} V(t)=\displaystyle\int_0^L(\xi_1^2p_1{\rm e}^{-\mu x}+\xi_2^2p_2{\rm e}^{\mu x}){\rm d}x+\displaystyle\int_0^L{\rm e}^{-\mu x}(p_3u^2+p_4v^2){\rm d}x, \end{equation}
(4.2)

其中, p_1,\;p_2,\;p_3,\;p_4,\;\mu均为正常数. 沿系统(2.12)对该Lyapunov函数求导得

\begin{equation} \dot{V}(t)=\dot{V_1}(t)+\dot{V_2}(t), \end{equation}
(4.3)

\begin{equation} \dot{V_1}(t)=-\big[\lambda_1p_1\xi_1^2{\rm e}^{-\mu x}-\lambda_2p_2\xi_2^2{\rm e}^{\mu x}+\frac{1}{\tau}p_3 Lu^2{\rm e}^{-\mu x} +\frac{1}{\tau}p_4 Lv^2{\rm e}^{-\mu x}\big]_0^L,\\ \end{equation}
(4.4)
\begin{equation} \dot{V_2}(t)=-\displaystyle\int_0^L M^T \Lambda M{\rm d}x, \end{equation}
(4.5)

其中

\begin{equation*} M= \left( \begin{array}{c} \xi_1 \\ \xi_2\\ u\\ v \end{array} \right),\; \Lambda= \left( \begin{array}{cccc} (\lambda_1\mu+2\gamma)p_1{\rm e}^{-\mu x} & \delta p_1{\rm e}^{-\mu x}+\gamma p_2{\rm e}^{\mu x} &0 & 0\\ \delta p_1{\rm e}^{-\mu x}+\gamma p_2{\rm e}^{\mu x} & (\lambda_2\mu+2\delta)p_2{\rm e}^{\mu x} & 0 & 0\\ 0 & 0 & \frac{1}{\tau}p_3\mu L {\rm e}^{-\mu x} &0\\ 0 & 0 & 0 & \frac{1}{\tau}p_4\mu L{\rm e}^{-\mu x} \end{array} \right). \end{equation*}

将(2.12)式中的边界条件代入\dot{V_1}(t)中得

\begin{matrix} \dot{V_1}(t) &=&(-\lambda_1p_1{\rm e}^{-\mu L}+\lambda_2p_2{\rm e}^{\mu L} k_2^2)\xi_1^2(t,L)+\big[(\lambda_1p_1+\frac{1}{\tau}p_3L)k_0^2-\lambda_2p_2+\frac{1}{\tau}p_4L\big]\xi_2^2(t,0) \\ &&+2(\lambda_1p_1+\frac{1}{\tau}p_3L)k_0k_1\xi_2(t,0)(u(t,L)-v(t,L))+(\lambda_1p_1+\frac{1}{\tau}p_3L)k_1^2(u(t,L) \\ &&-v(t,L))^2-\frac{1}{\tau}p_3Lu^2(t,L){\rm e}^{-\mu L}-\frac{1}{\tau}p_4Lv^2(t,L){\rm e}^{-\mu L} \\ &\leq&(-\lambda_1p_1{\rm e}^{-\mu L}+\lambda_2p_2{\rm e}^{\mu L} k_2^2)\xi_1^2(t,L) \\ &&+\big[(\lambda_1p_1+\frac{1}{\tau}p_3L)k_0^2-\lambda_2p_2+\frac{1}{\tau}Lp_4+2(\lambda_1p_1+\frac{1}{\tau}p_3L)\left|k_0k_1\right|\big]\xi_2^2(t,0) \\ &&+\big[(\lambda_1p_1+\frac{1}{\tau}p_3L)\left|k_0k_1\right|+2(\lambda_1p_1+\frac{1}{\tau}p_3L)k_1^2-\frac{1}{\tau}p_3L{\rm e}^{-\mu L}\big]u^2(t,L) \\ &&+\big[(\lambda_1p_1+\frac{1}{\tau}p_3L)\left|k_0k_1\right|+2(\lambda_1p_1+\frac{1}{\tau}p_3L)k_1^2-\frac{1}{\tau}p_4L{\rm e}^{-\mu L}\big]v^2(t,L). \end{matrix}
(4.6)

显然, 当参数满足条件

\begin{equation} \begin{cases} \lambda_2p_2k_2^2{\rm e}^{\mu L}\leq\lambda_1p_1{\rm e}^{-\mu L}, (\frac{1}{\tau}p_3L+\lambda_1p_1)(k_0^2+2\left|k_0k_1\right|)+\frac{1}{\tau}p_4L\leq\lambda_2p_2,\\ (\frac{1}{\tau}p_3L+\lambda_1p_1)(2k_1^2+\left|k_0k_1\right|)\leq \frac{1}{\tau}L{\rm e}^{-\mu L}\min\{p_3,p_4\} \end{cases} \end{equation}
(4.7)

时, 必有\dot{V_1}(t)\leq0.

为使\dot{V_2}(t)<0, 当且仅当

\begin{equation} \left( \begin{array}{cc} (\lambda_1\mu+2\gamma)p_1{\rm e}^{-\mu x}\quad & \delta p_1{\rm e}^{-\mu x}+\gamma p_2{\rm e}^{\mu x} \\ \delta p_1{\rm e}^{-\mu x}+\gamma p_2{\rm e}^{\mu x}\quad & (\lambda_2\mu+2\delta)p_2{\rm e}^{\mu x} \end{array} \right)>0, \end{equation}
(4.8)

即各阶顺序主子式均大于0

\begin{equation} \begin{array}{ll} {\rm (a)}\quad (\lambda_1\mu+2\gamma)p_1{\rm e}^{-\mu x}>0,\\ {\rm (b)}\quad p_1p_2(\lambda_1\mu+2\gamma)(\lambda_2{\mu}+2\delta)-(\delta p_1{\rm e}^{-\mu x}+\gamma p_2{\rm e}^{\mu x})^2>0. \end{array} \end{equation}
(4.9)

显然条件(4.9)(a)成立. 对于条件(4.9)(b), 若选择参数p_1,\;p_2使得\delta p_1=\gamma p_2, 则函数(\delta p_1{\rm e}^{-\mu x}+\gamma p_2{\rm e}^{\mu x})^2x=L处取得最大值, 于是,

\begin{equation*} p_1p_2(\lambda_1\mu+2\gamma)(\lambda_2\mu+2\delta)-(\delta p_1{\rm e}^{-\mu x}+\gamma p_2{\rm e}^{\mu x})^2 \\ > p_1p_2(\lambda_1\mu+2\gamma)(\lambda_2\mu+2\delta)-(\delta p_1{\rm e}^{-\mu L}+\gamma p_2{\rm e}^{\mu L})^2\\ \\ p_1p_2\big[\mu^2\lambda_1\lambda_2+2\mu(\lambda_1\delta+\lambda_2\gamma)\big]-\delta^2 p_1^2({\rm e}^{-\mu L}-{\rm e}^{\mu L})^2. \end{equation*}

考虑到

{\rm e}^{-\mu L}=1-\mu L+\frac{\mu^2L^2}{2}-\frac{\mu^3L^3}{3!}+\cdots,
{\rm e}^{\mu L}=1+\mu L+\frac{\mu^2L^2}{2}+\frac{\mu^3L^3}{3!}+\cdots,

因此当\mu>0充分小时, ({\rm e}^{-\mu L}-{\rm e}^{\mu L})^2是关于\mu^2的同阶无穷小, 从而必有

p_1p_2\big[\mu^2\lambda_1\lambda_2+2\mu(\lambda_1\delta+\lambda_2\gamma)\big]-\delta^2 p_1^2({\rm e}^{-\mu L}-{\rm e}^{\mu L})^2>0.

于是(4.9)(b)成立.

综上可得, 必存在正常数\nu使得\dot{V}(t)\leq-\nu V(t). 故结论成立.

注4.1 若选择如下Lyapunov候选函数

\begin{equation} V(t)=\displaystyle\int_0^L(\xi_1^2p_1{\rm e}^{-\mu x}+\xi_2^2p_2{\rm e}^{\mu x}){\rm d}x+\frac{p_5\tau}{L}\displaystyle\int_0^L{\rm e}^{-\mu x}(u^2+v^2){\rm d}x, \end{equation}
(4.10)

即(4.2)式中p_3=p_4的情形, 则参数的稳定性条件(4.1)便转化为

\begin{equation} k_2^2\leq \frac{\lambda_1\gamma}{\lambda_2\delta}{\rm e}^{-2\mu L},\; p_5+\lambda_1p_1\leq \min\big\{\frac{p_5{\rm e}^{-\mu L}}{2k_1^2+\left|k_0k_1\right|}, \frac{\delta\lambda_2p_1-p_5\gamma}{\gamma(k_0^2+2\left|k_0k_1\right|)}\big\}. \end{equation}
(4.11)

5 系统算子的谱分析

在这一节, 旨在分析系统算子 \mathcal{A} 的谱分布. 考虑特征值问题

\begin{equation*} \mathcal{A} X=\mu X,\;\;X=(f,g,h,q)\in D(\mathcal{A}), \end{equation*}

f,g,h,q 满足

\begin{equation} \begin{cases} -\lambda_1f'-\gamma f-\delta g=\mu f,\;\lambda_2g'-\gamma f-\delta g=\mu g,\;-\frac{1}{\tau}h'=\mu h,\;-\frac{1}{\tau}q'=\mu q, \\ f(0)=k_0g(0)+k_1(h(1)-q(1)),\;g(1)=k_2f(1),\;h(0)=f(0),\;q(0)=g(0). \end{cases} \end{equation}
(5.1)

类似于文献[定理3]的证明, 易得如下结论成立.

定理5.1\mathcal{A} 由(2.15)和(2.16)式给出. 则 \sigma(\mathcal{A})=\sigma_p(\mathcal{A})=\{\mu_n,n\in Z^+\}, 且 \mu_n 的渐近表达式为

\begin{equation} \mu_n= \begin{cases} -\frac{B}{A}+\frac{\ln\sqrt{k_0k_2}}{A}+\frac{n\pi {\rm i}}{A}+\mathcal{O}(n^{-1}), &\;\;n\rightarrow\infty, \hbox{若}\;\; k_0k_2>0;\\ -\frac{B}{A}+\frac{\ln\sqrt{-k_0k_2}}{A}+\frac{(n+\frac{1}{2})\pi {\rm i}}{A}+\mathcal{O}(n^{-1}),&\;\;n\rightarrow\infty, \hbox{若}\;\; k_0k_2<0, \end{cases} \end{equation}
(5.2)

其中

\begin{equation} \begin{aligned} A=\frac{\lambda_1+\lambda_2}{2\lambda_1\lambda_2},\;\;\;\;B=\frac{\lambda_1\delta+\lambda_2\gamma}{2\lambda_1\lambda_2}. \end{aligned} \end{equation}
(5.3)

此外, 对应的特征函数 X_n=(f_n,g_n,h_n,q_n) 的渐近表达式为

\begin{equation} \begin{aligned} \begin{cases} f_n(x)=P(x)(k_0k_2)^{-\frac{x}{2}}{\rm e}^{-\frac{2\lambda_2}{\lambda_1+\lambda_2}\cdot n\pi {\rm i}x}+\mathcal{O}(n^{-1}),\\ g_n(x)=\frac{1}{k_0}P(x)(k_0k_2)^{\frac{x}{2}}{\rm e}^{\frac{2\lambda_1}{\lambda_1+\lambda_2}\cdot n\pi {\rm i}x}+\mathcal{O}(n^{-1}),\\ h_n(x)=\mathcal{O}(n^{-1}),\;\; q_n(x)=\mathcal{O}(n^{-1}), \end{cases} \end{aligned} \end{equation}
(5.4)

其中

\begin{equation} P(x)=\exp\Big\{-\frac{(\lambda_2-\lambda_1)\ln\sqrt{k_0k_2}+\gamma-\delta}{\lambda_1+\lambda_2}x\Big\}. \end{equation}
(5.5)

注5.1 特征值的渐近表达式(5.2)表明, 特征值的实部趋于某一常数, 虚轴不是特征值的渐近线. 此外, 基于特征函数的渐近表达式(5.4), 可知Riesz基性质和谱确定性增长条件成立, 进而结合系统的耗散性也可建立系统(2.17)的指数稳定性, 可参考文献[31].

6 数值模拟

在这一节, 对于由 Saint-Venant 方程 (2.1) 描述的具有底部坡度和摩擦的单渠道系统进行数值仿真, 参数如下: 渠道长度 L=1000 m, 宽度 W=80 m, 底部坡度 S_b=0.0002, 摩擦系数 C_f=0.001 s^2/{\rm m}, 稳态流率 Q^*=400\;{\rm m}^3/{\rm s}. 根据

\begin{equation*} S_b H^*=C_f V^{*2},\;Q^*=WH^*V^*, \end{equation*}

计算得稳态水深和流速分别是 H^*=5 m, V^*=1\;{\rm m}/{\rm sec}. 取系统参数为

\begin{equation*} \begin{aligned} \lambda_1=8\;{\rm m}/{\rm s},\;\lambda_2=6\;{\rm m}/{\rm s},\;\gamma=0.0018\;{\rm s}^{-1},\;\delta=0.0021\;{\rm s}^{-1}. \end{aligned} \end{equation*}

\tau=2, Lyapunov 函数中的参数值分别为

\begin{equation*} p_1=0.18,\;p_2=0.21,\;p_3=0.02,\;p_4=0.001,\;\mu=0.001, \end{equation*}

满足 \delta p_1=\gamma p_2. 根据定理 4.1 中的条件 (4.1) 可知, 如果反馈参数满足条件

\begin{equation*} |k_2|\leq0.395, k_0^2+2|k_0k_1|\leq0.07, 2k_1^2+|k_0k_1|\leq0.016, \end{equation*}

则闭环系统 (2.12) 的稳定性可以保证. 图2图3 展示了调控参数 k_0=0.1,\;k_1=0.05,\;k_2=0.3, 初始条件为 \xi_1(0,x)=\cos(0.002\pi x)-\frac{1}{2},\;\xi_2(0,x)=\sin(0.002\pi x)+\frac{1}{2} 时系统 (2.12) 的稳定性. 图4 则展示了 (2.8) 中 PDP 控制器 U(t) 的稳定性, 其中 K_p=-0.155,\;K_d=0.147.

图2

图2   系统 (2.12) 中状态 \xi_1(t,x) 的收敛性


图3

图3   系统 (2.12) 中状态 \xi_2(t,x) 的收敛性


图4

图4   (2.8) 式控制器 U(t) 的收敛性


7 结论

本文研究了由Saint-Venant方程描述的具有非零底部坡度和底部摩擦的单渠道系统的PDP边界控制器设计. 首先通过线性化手段和Riemann坐标变换将系统重构, 得到一类2\times2双曲型偏微分方程, 进而将其改写为抽象发展方程的形式. 然后, 利用算子半群理论证明了在一定条件下闭环系统解的适定性. 接着, 构造了严格Lyapunov函数求得系统反馈参数和时滞参数的一个充分条件, 保证了在亚临界流动条件下闭环系统的指数稳定性, 无需对底部摩擦以及坡度大小做额外的假设. 更进一步, 利用谱分析方法, 给出了系统特征值和特征函数的渐近表达式. 下一步的研究工作是将单渠道系统扩展为具有n段级联渠道, 星形网络渠道或树形网络的情形, 并考虑控制器中的时滞因素对系统稳定性的影响.

参考文献

Litrico X, Fromion V, Baume J P, et al.

Experimental validation of a methodology to control irrigation canals based on Saint-Venant equations

Control Engineering Practice, 2005, 13(11): 1425-1437

DOI:10.1016/j.conengprac.2004.12.010      URL     [本文引用: 1]

Santos V D, Prieur C.

Boundary control of open channels with numerical and experimental validations

IEEE Transactions on Control Systems Technology, 2008, 16(6): 1252-1264

DOI:10.1109/TCST.2008.919418      URL     [本文引用: 1]

Bastin G, Coron J M, Tamasoiu S O.

Stability of linear density-flow hyperbolic systems under PI boundary control

Automatica, 2015, 53: 37-42

DOI:10.1016/j.automatica.2014.12.025      URL     [本文引用: 1]

Zhang L G, Prieur C, Qiao J F.

PI boundary control of linear hyperbolic balance laws with stabilization of ARZ traffic flow models

Systems & Control Letters, 2019, 123: 85-91

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

Bastin G, Coron J M.

Stability and Boundary Stabilization of 1-D Hyperbolic Systems

Switzerland: Birkhäuser, 2016

[本文引用: 2]

Vazquez R, Krstic M, Coron J M.

Backstepping boundary stabilization and state estimation of a 2\times2 linear hyperbolic system

Orlando: IEEE Conference on Decision and Control and European Control Conference, 2011

[本文引用: 1]

Litrico X, Fromion V.

Boundary control of hyperbolic conservation laws using a frequency domain approach

Automatica, 2009, 45(3): 647-656

DOI:10.1016/j.automatica.2008.09.022      URL     [本文引用: 1]

Bastin G, Coron J M, D'Andrea-Novel B, et al. Boundary control for exact cancellation of boundary disturbances in hyperbolic systems of conservation laws. Seville: IEEE Conference on Decision and Control, and the European Control Conference, 2005

[本文引用: 1]

Bastin G, Coron J M.

Exponential stability of PI control for Saint-Venant equations with a friction term

Methods and Applications of Analysis, 2019, 26(2): 101-112

DOI:10.4310/MAA.2019.v26.n2.a1      URL     [本文引用: 1]

Coron J M, D'Andrea-Novel B, Bastin G.

A strict Lyapunov function for boundary control of hyperbolic systems of conservation laws

IEEE Transactions on Automatic Control, 2007, 52(1): 2-11

DOI:10.1109/TAC.2006.887903      URL     [本文引用: 1]

Pazy A. Semigroup of Linear Operators and Applications to Partial Differential Equations. New York: Springer-Verlag, 1983

[本文引用: 1]

Hayat A, Shang P P.

A quadratic Lyapunov function for Saint-Venant equations with arbitrary friction and space-varying slope

Automatica, 2019, 100: 52-60

DOI:10.1016/j.automatica.2018.10.035      URL     [本文引用: 1]

Coron J M, Tamasoiu S O.

Feedback stabilization for a scalar conservation law with PID boundary control

Chinese Annals of Mathematics, 2015, 36B(5): 763-776

[本文引用: 1]

Chentouf B, Wang J M.

Boundary feedback stabilization and Riesz basis property of a 1-d first order hyperbolic linear system with L^\infty-coefficients

Journal of Differential Equations, 2009, 246(3): 1119-1138

DOI:10.1016/j.jde.2008.08.010      URL     [本文引用: 1]

Trinh N T, Andrieu V, Xu C Z.

Boundary PI controllers for a star-shaped network of 2\times2 systems governed by hyperbolic partial differential equations

International Federation of Automatic Control, 2017, 50(1): 7070-7075

[本文引用: 1]

Trinh N T, Andrieu V, Xu C Z.

Output regulation for a cascaded network of 2\times2 hyperbolic systems with PI controller

Automatica, 2018, 91: 270-278

DOI:10.1016/j.automatica.2018.01.010      URL     [本文引用: 1]

Aström K J, Hägglund T.

PID Controllers: Theory, Design, and Tunning

Research Triangle Park, NC: ISA-The Instrumentation, Systems and Automation Society, 1995

[本文引用: 1]

Litrico X, Fromion V, Scorletti G.

Robust feedforward boundary control of hyperbolic conservation laws

Networks and Heterogeneous Media, 2017, 2(4): 717-731

DOI:10.3934/nhm.2007.2.717      URL     [本文引用: 1]

Bastin G, Coron J M, Hayat A.

Feedforward boundary control of 2\times2 nonlinear hyperbolic systems with application to Saint-Venant equations

European Journal of Control, 2021, 57: 41-53

DOI:10.1016/j.ejcon.2020.11.002      URL     [本文引用: 1]

Tallman G, Smith O.

Analog study of dead-beat posicast control

Ire Transactions on Automatic Control, 2003, 4(1): 14-21

DOI:10.1109/TAC.1958.1104844      URL     [本文引用: 1]

Silva G J, Datta A, Bhattacharyya S P. PID Controllers for Time-Delay Systems. Boston: Birkhäuser, 2005

[本文引用: 1]

Cai G, Chen L.

Some problems of delayed feedback control

Advances in Mechanics, 2013, 43(1): 21-28

[本文引用: 1]

Suh I H, Bien Z.

Proportional minus delay controller

IEEE Transactions on Automatic Control, 1979, 24(2): 370-372

DOI:10.1109/TAC.1979.1102024      URL     [本文引用: 2]

Suh I H, Bien Z.

Use of time-delay actions in the controller design

IEEE Transactions on Automatic Control, 1980, 25: 600-603

DOI:10.1109/TAC.1980.1102347      URL     [本文引用: 2]

Atay F M.

Balancing the inverted pendulum using position feedback

Applied Mathematics Letters, 1999, 12(5): 51-56

[本文引用: 1]

Hu H Y.

Using delayed state feedback to stabilize periodic motions of an oscillator

Journal of Sound and Vibration, 2004, 275(3-5): 1009-1025

DOI:10.1016/j.jsv.2003.07.006      URL     [本文引用: 1]

Liu B, Hu H.

Stabilization of linear undamped systems via position and delayed position feedbacks

Journal of Sound and Vibration, 2008, 312(3): 509-525

DOI:10.1016/j.jsv.2007.11.001      URL     [本文引用: 1]

Krstic M.

Control of an unstable reaction-diffusion PDE with long input delay

Systems & Control Letters, 2009, 58(10/11): 773-782

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

范东霞, 赵东霞, 史娜, .

一类扩散波方程的PDP反馈控制和稳定性分析

数学物理学报, 2021, 41A(4): 1088-1096

[本文引用: 1]

Fan D X, Zhao D X, Shi N, et al.

The PDP feedback control and stability analysis of a diffusive wave equation

Acta Mathematica Scientia, 2021, 41A(4): 1088-1096

[本文引用: 1]

Bastin G, Coron J M, D'Andrea-Novel B.

On Lyapunov stability of linearised Saint-Venant equations for a sloping channel

Networks and Heterogeneous Media, 2009, 4(2): 177-187

DOI:10.3934/nhm.2009.4.177      URL     [本文引用: 2]

Zhao D X, Fan D X, Guo Y P.

The spectral analysis and exponential stability of a 1-d 2\times2 hyperbolic system with proportional feedback control

International Journal of Control, Automation and Systems, 2022, 20(8): 2633-2640

DOI:10.1007/s12555-021-0507-0      [本文引用: 1]

/