非线性ODE-PDE耦合系统边界控制的局部镇定①

何翠华， 王科， 甯懿楠

成都工业学院 大数据与人工智能学院， 成都 611730

耦合的ODE-PDE系统广泛存在于工程问题中． 近年来， 大量文献采用 Backstepping 变换方法来解决常微分方程和偏微分方程之间的耦合问题． 为满足工程需要， 我们利用边界条件设计控制律来控制这些耦合的ODE-PDE系统使其稳定在平衡状态． 边界控制律是由PDE Backstepping 控制发展而来[1-4]．

ODE-PDE耦合系统边界控制大多是线性的[5-15]， 然而非线性现象在应用数学和物理中广泛存在． 耦合的线性常微分方程和非线性偏微分方程系统具有丰富的工程应用价值， 但由于非线性在数学上是极具挑战性的问题， 故关于非线性偏微分方程系统镇定的结果较少[15-22]． 利用边界控制来稳定非线性ODE-PDE耦合系统是一个很有意义的研究领域．

本文对非线性ODE-PDE耦合系统的控制设计类似于文献[7]， 主要贡献是基于边界反馈控制的方法处理一类包含反应项是不确定的非线性偏微分方程与线性微分方程的耦合系统的局部指数稳定的问题．

因此， 本文考虑以下非线性耦合系统的边界控制

(1)

ut(x，t)=uxx(x，t)+f(u(x，t))，x∈(0， 1)

(2)

ux(0，t)=α(u(0，t)-CTX(t))

(3)

u(1，t)=U(t)

(4)

其中： 向量X(t)∈Rn是装置的一个信号； (A，B)是可控制矩阵对 (A∈Rn×n，B∈Rn) ； 标量u(x，t)∈R是装置的热量； 在x=0处， 满足Neumann边界条件，α是傅里叶常数， 它取决于装置的材料和导热性，CTX(t)是温度装置(CT∈Rn)；U(t)是边界控制输入．f(u(x，t))是非线性函数， 对于这个系统， 非线性f(u)引用文献[15]假定条件．

假定1对于系统(2)中的f(u)， 存在δ>0和β>1， 使得

f(0)=0

‖f(u)‖≤δ‖u‖β

(5)

其中‖·‖将在后面详细说明．

假定1中的第一项表示0是系统(1)-(4)在零输入下的平衡态， 第二项意味着f(u)的增长速度并不比u的幂函数快．

系统(1)-(4)代表了一类广泛的耦合非线性ODE-PDE方程， 本文主要贡献是基于 Backstepping 变换的边界控制方法使一类包含不确定反应项的非线性偏微分方程与线性微分方程的耦合系统达到局部指数稳定． 对于ODE-PDE耦合系统Backstepping边界控制律设计方法是可行的[1，6-7]． 通过引入具有理想稳定性的目标系统， PDE Backstepping 变换将系统转化为由核函数和向量值函数耦合的方程． 然后根据目标系统的边界条件得到控制律． 研究PDE Backstepping 变换解决耦合ODE-PDE的控制系统是很有价值的．

本文令

(6)

(7)

1 控制律

1.1 Backstepping变换

引入以下Backstepping变换

X(t)=X(t)

(8)

(9)

其中k(x，y)与φ(x)是待定的．

将目标系统取为

(10)

wt(x，t)=wxx(x，t)+H(w(x，t))，x∈(0， 1)

(11)

wx(0，t)=0

(12)

w(1，t)=0

(13)

其中：

选择向量KT∈Rn使得 (A+BK)是Herwitz矩阵．

利用(3)，(9)式对w(x，t)关于x和t求导

(14)

再令

k′(x，x)=0

ky(x， 0)-αk(x， 0)+φ(x)B=0

φ″(x)-αk(x， 0)CT-φ(x)A=0

kxx(x，y)-kyy(x，y)=0

(15)

接着控制律被设计为

(16)

通过边界条件(13)获得．

1.2 核方程的解

为了方便， 引用文献[14]直接给出(15)式的解． 令

其中I是单位矩阵．

Γ(0)=(K， -αCT，KA+α2CT)

核函数k(x，y)和φ(x)分别为

φ(x)=Γ(0)eDxJ

1.3 逆变换

同样我们给出逆变换

X(t)=X(t)

(17)

(18)

令

其中I是单位矩阵．

Θ(0)=(K，α(K-CT))

核函数Ψ(x)和n(x，y)分别为

Ψ(x)=Θ(0)eΩxF

2 稳定性

为得到稳定性定理， 首先给出几个引理．

引理1函数w(x，t)由(6)式定义， (PB)T和X(t)由(7)式定义， 有以下不等式成立

(19)

证现在利用Schwartz’s不等式和Young’s不等式来估计， 可得

接着根据w(1，t)=0与Agmon’s不等式， 有

然后利用Poincare不等式， 可得

因此， 不等式(19)成立．

引理2根据变换(9)与逆变换(18)， 可得

‖w(x，t)‖2≤(1+s)‖u(t)‖2+‖φ‖2‖X(t)‖

‖u(x，t)‖2≤(1+t)‖w(t)‖2+‖Ψ‖2‖X(t)‖

(20)

其中：

(21)

(22)

证由变换(9)， 可以获得

由Holder’s不等式， 可以得到

其中：

(23)

接着根据Schwartz’s不等式， 有

其中：

因此

‖w(x，t)‖2≤(1+s)‖u(t)‖2+‖φ‖2‖X(t)‖

同理可估计

‖u(x，t)‖2≤(1+t)‖w(t)‖2+‖Ψ‖2‖X(t)‖

其中t和‖Ψ‖2由(22)式定义， 则(20)式成立．

引理3H(w(x，t))由(11)式定义， 可得

(24)

其中：

(25)

证

我们令

其中：

利用Holder不等式 与Young’s不等式， 可得

其中：

由此可得

同理可估计

其中s来自(23)式．

然后

其中ξ是如(25)式所示的常数．

接下来证明目标系统(10)-(13)是局部稳定的．

定理1存在一个正常数γ且任意的初始条件满足‖X(0)，w(0)‖2≤γ， 有下列不等式成立

(26)

其中

(27)

则目标系统(10)-(13)在‖·‖2意义下是局部指数稳定的．

证

考虑以下Lyapunov函数

(28)

其中矩阵P=PT>0满足

P(A+BK)+(A+BK)TP=-I

(29)

且参数a>0是待定的．

矩阵P是(29)式的解， 因此

(30)

其中λmin(P)与λmax(P)是矩阵P的最小和最大的特征值．

对(28)式关于t求导， 则

代入(19)，(20)和(24)式可以得到

令a=8‖PB‖2+3， 可得

其中ξ由(25)式定义．

(31)

选取‖X(0)，w(0)‖2≤γ， 则V(0)≤σγ2， 因此不等式(31)成立． 利用(30)和(31)式， 可以推出

其中μ由(27)式定义， 因此不等式(26)成立． 证毕．

3 仿真设计

系统(1)-(4)的模拟仿真

其中：

f(u(x，t))=12u2(x，t)

同时

我们令

K=(-0.7， -9.2)

可以验证矩阵A+BK是Hurwitz矩阵． 因此， 系统满足(5)式， 容易得到

其中：

Γ(0)=(-0.7， -9.2， 0， -5， -18.75， 13.82)

则控制律为

(32)

仿真结果见图1-4． 图1开环系统是不稳定的． 图2和图3显示在控制(32)输入下， 闭环系统的信号在时间趋近于无穷时趋近于零． 图4显示控制(32)是有界的．

图1 开环系统下的X(t)

图2 闭环系统下的X(t)

图3 闭环系统下的u(x， t)

图4 U(t)的输入