2.845

2023影响因子

(CJCR)

  • 中文核心
  • EI
  • 中国科技核心
  • Scopus
  • CSCD
  • 英国科学文摘

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

时滞分布参数系统中和控制器设计

周笔锋 罗毅平

周笔锋, 罗毅平. 时滞分布参数系统中和控制器设计. 自动化学报, 2018, 44(12): 2222-2227. doi: 10.16383/j.aas.2018.c170084
引用本文: 周笔锋, 罗毅平. 时滞分布参数系统中和控制器设计. 自动化学报, 2018, 44(12): 2222-2227. doi: 10.16383/j.aas.2018.c170084
ZHOU Bi-Feng, LUO Yi-Ping. Neutralization Control of Distributed Parameter Systems With Delay. ACTA AUTOMATICA SINICA, 2018, 44(12): 2222-2227. doi: 10.16383/j.aas.2018.c170084
Citation: ZHOU Bi-Feng, LUO Yi-Ping. Neutralization Control of Distributed Parameter Systems With Delay. ACTA AUTOMATICA SINICA, 2018, 44(12): 2222-2227. doi: 10.16383/j.aas.2018.c170084

时滞分布参数系统中和控制器设计

doi: 10.16383/j.aas.2018.c170084
基金项目: 

国家自然科学基金 11372107

湖南省自然科学基金 2017JJ4004

湖南省自然科学基金 20 17JJ5011

风力发电机组及控制湖南省重点实验室(湖南工程学院)开放基金 2017FLFDYB03

详细信息
    作者简介:

    周笔锋  湖南电气职业技术学院教师.2015年获得湖南工程学院硕士学位.主要研究方向为复杂网络系统, 分布参数系统.E-mail:zhoubifeng99@163.com

    通讯作者:

    罗毅平  湖南工程学院教授.2006年获得华南理工大学博士学位.主要研究方向为复杂网络系统, 分布参数系统.本文通信作者.E-mail:lyp8688@sohu.com

Neutralization Control of Distributed Parameter Systems With Delay

Funds: 

National Natural Science Foundation of China 11372107

Natural Science Foundation of Hunan Province 2017JJ4004

Natural Science Foundation of Hunan Province 20 17JJ5011

Hunan Provincial Key Laboratory of Wind Generator and Its Control 2017FLFDYB03

More Information
    Author Bio:

       Teacher at Hunan Electrical College of Technology. He received his master degree from Hunan Institute of Engineering in 2015. His research interest covers complex networks and distributed parameter systems

    Corresponding author: LUO Yi-Ping   Professor at Hunan Institute of Engineering. He received his Ph. D. degree from South China University of Technology in 2006. His research interest covers complex networks and distributed parameter systems. Corresponding author of this paper
  • 摘要: 首先给出中和控制器设计思路,针对具有时滞特性的分布参数系统,设计中和控制器,讨论此类系统的稳定问题.利用Lyapunov稳定性理论并结合矩阵不等式处理方法,得出具有时滞特性分布参数系统稳定中和控制器存在的充分条件.最后结合所给条件,给出一个数值仿真说明其有效性.
  • 实际生活中许多物理系统都具有时空特性, 例如热扩散、流体换热器、化学工程、旋转梁、可变几何形状静电微致动器、集成和消防神经元等[1-7], 其行为必须依赖于时间和空间位置, 这些系统的时空过程被称为分布参数系统(Distributed parameter systems, DPSs).为了更好地研究此类系统, 根据能量守恒定律, 通常构建拟线性抛物型偏微分方程(Quasi-linear parabolic partial differential equation)或拟线性抛物型偏微分-积分方程(Quasi-linear partial differential-integral equations).构建拟线性抛物型偏微分方程研究分布参数系统稳定性一直是国内外相关领域学者的重点研究课题[8-15].

    文献[8]针对拟线性抛物型偏微分方程设计控制器, 利用Lyapunov稳定性定理结合线性矩阵不等式(Linear matrix inequality, LMI)计算方法, 得出了分布参数系统指数稳定控制器存在的充分条件.文献[9]针对一类同时具有变时滞和连续分布时滞的分布参数系统的状态反馈控制问题进行了研究, 通过选择适当的Lyapunov-Krasovskii函数, 采用LMI方法, 得到了变时滞闭环系统渐近稳定的一个充分条件.文献[10]针对具有时滞性的抛物型偏微分方程, 利用滑模控制(Sliding mode control, SMC), 研究了系统的稳定性.文献[11]针对一类具有耦合性质的时滞分布参数系统, 构建拟线性抛物型偏微分方程并结合Lyapunov稳定性定理, 研究了系统节点的同步问题.文献[12]针对工艺参数未知的化学分布参数系统, 构建半线性偏微分方程, 利用Lyapunov稳定性定理, 研究了其自适应输出反馈控制问题.文献[13]基于边界观测和边界控制, 针对扩散抛物型偏微分方程, 运用自适应控制方案, 研究了系统的同步稳定性.

    分布参数系统分布式控制是针对分布参数系统中的状态变量设计控制器, 达到控制的目的.这种方法的好处在于控制器设计简单, 但实际应用中有一定难度.考虑现实生活中存在另一类控制方法, 例如工厂的工业废水处理, 工业废水中污染物的种类较多, 较常见的有酸碱污染物, 对于工业废水, 在进入河道之前需进行一系列处理.对于工业废水中的酸碱污染物, 比较常用的方法是用相对应的化学物质进行中和.例如对于酸性污水需用碱性物质来中和, 首先, 溶液中酸性离子的分布可以看成一个分布参数系统, 因为其浓度不是均匀分布的; 其次, 在中和过程中, 渗入工业废水中的碱性离子的扩散也可以看作是一个分布参数系统, 对于不同的时间, 不同的位置, 碱性离子的密度不尽相同, 若添加了过量的碱性物质, 污水又将呈现碱性, 污染环境; 若碱性物质的量过少, 则存在部分酸性离子不能中和, 同样对环境存在危害.基于此, 若将碱性物质在工业废水中的扩散现象构建成一个合适的分布参数系统, 用偏微分方程描述, 在构建的系统中同时考虑酸性离子对碱性物质的影响, 将碱性物质作为控制酸性离子的控制器, 这种控制方法称为中和控制.

    在实际系统中, 时滞现象普遍存在[9-11, 16-18], 在系统的信息传递和信号检测过程中通常具有滞后现象.因此, 从时滞分布参数系统出发, 设计中和控制器, 研究时滞分布参数系统的中和控制就显得尤为重要.基于此, 本文针对具有变时滞特性的分布参数系统, 设计中和控制器, 讨论此类系统的稳定问题.利用Lyapunov稳定性理论并结合LMI处理方法, 得出了具有时滞特性分布参数系统稳定中和控制器存在的充分条件.最后结合所给条件, 给出一个数值仿真说明其有效性.

    考虑下列具有时滞特性的分布参数系统

    $ \begin{align} \frac{\partial w_{i}(x, t)}{\partial t}=&\ D_{a}\sum\limits_{k=1}^{m}\frac{\partial^2 w_{i}(x, t)}{\partial x_k^2} +\nonumber\\ &\ \sum\limits_{j=1}^{n}a_{ij}^{0}w_{j}(x, t) +\nonumber\\ &\ \sum\limits_{j=1}^{n}a_{ij}w_{j}(x, t-\tau)+u_{i}, \nonumber\\ &\qquad\qquad\qquad\quad i=1, 2, \cdots, n \end{align} $

    (1)

    将系统变为矩阵形式, 即为

    $ \begin{align} \frac{\partial W(x, t)}{\partial t}=&\ D_{a}\Delta W(x, t)+A_{0} W(x, t) +\nonumber\\ &\ A W(x, t-\tau)+U(x, t) \end{align} $

    (2)

    其中, $(x, t)\in \Omega \times {\bf R}_{+}$, $D_{a}>0$, $\tau$为时滞, $A_{0}=$ $(a_{ij}^{0})_{n \times n}$和$A=(a_{ij})_{n \times n}$为常数矩阵, $\Omega = \{ x, | x |$ $ <$ 为具有光滑边界$\Theta$的有界区域, 且$mes \Omega>0$ (表示$\Omega$的测度), $U(x, t)=$ $(u_{1}(x, t)$, 为控制输入, 状态变量$W(x, t)=(w_{1}(x, t), \cdots, w_{n}(x, t))\in {\bf R}^{n}$, $\Delta=$ $\sum_{k=1}^{m}\frac{\partial^2}{\partial x_k^2}$为Laplace扩散-吸收算子, 系统初始边界条件为

    $ \begin{align} W(x, t)=0, \quad (x, t)\in \Theta \times [0, +\infty ) \end{align} $

    (3)

    $ \begin{align} & \frac{\partial W(x, t)}{\partial t}=0, && (x, t)\in \Theta \times [0, +\infty ) \end{align} $

    (4)

    $ \begin{align} & W(x, t)=\psi (x, t), && (x, t)\in \Theta \times [-\tau, 0 ) \end{align} $

    (5)

    其中, $n$为$\Theta$的单位外法向量, $\psi (x, t)$为适当光滑的函数.

    为使系统(1)达到稳定, 设计如下中和控制器

    $ \begin{align} \begin{cases}\dfrac{\partial \eta_{i}(x, t)}{\partial t}=D_{c}\sum\limits_{k=1}^{m}\dfrac{\partial^2 \eta_{i}(x, t)}{\partial x_k^2}+\sum\limits_{j=1}^{n}b_{cij}\eta_{j}(x, t) +\\[3mm] \qquad\qquad\quad \sum\limits_{j=1}^{n}b_{cij}^{1}w_{j}(x, t)+ \sum\limits_{j=1}^{n}b_{cij}^{2}w_{j}(x, t-\tau)\\[5mm] u_{i}(x, t)=k_{i}\eta_{i}(x, t), \quad i=1, \cdots, n \end{cases} \end{align} $

    (6)

    将控制系统(6)变为矩阵形式, 即为

    $ \begin{align} \begin{cases} \dfrac{\partial \chi (x, t)}{\partial t}=D_{c}\Delta \chi(x, t)+ B_{c} \chi(x, t) +\\\qquad\qquad\ \ B_{c1} W(x, t)+B_{c2} W(x, t-\tau)\\[2mm] U(x, t)=K\chi(x, t) \end{cases} \end{align} $

    (7)

    其中, 为状态变量, $D_{c}>0$, $B_{c}=(b_{cij})_{n\times n}$, $B_{c1}=(b_{cij}^{1})_{n\times n}$, $B_{c2}$ $=$ $(b_{cij}^{2})_{n\times n}$为控制器参数, 控制律$U(x, t)$为系统(2)输入控制, 控制系统(7)初始边界条件为

    $ \begin{align} \chi(x, t)=0, \quad (x, t)\in \Theta \times [0, +\infty ) \end{align} $

    (8)

    $ \begin{align} & \frac{\partial \chi(x, t)}{\partial t}=0, &&(x, t)\in \Theta \times [0, +\infty ) \end{align} $

    (9)

    $ \begin{align} &\chi(x, t)=\varphi (x, t), && (x, t)\in \Theta \times [-\tau, 0) \end{align} $

    (10)

    其中, $n$为$\Theta$的单位外法向量, $\varphi (x, t)$为适当光滑的函数.

    引理1 (Schur complements)[19].对于给定的对称矩阵, 以下三个条件是等价的:

    1) $S <0$;

    2) $S_{11} < 0$, $S_{22}-S_{12}^{{\rm T}}S_{11}^{-1}S_{12} < 0$;

    3) $S_{22} < 0$, $S_{11}-S_{12}S_{22}^{-1}S_{12}^{{\rm T}} < 0$.

    引理2[20].设是边界$\Theta$内的光滑有界区域, $n$为$\Theta$的单位外法向量, $G\subset \Omega $为一光滑子域, 若$u, v \in C^{2}(\bar{G})$, 则

    $ \begin{align} \int_{G}u\Delta v{\rm d}x=\int_{\partial G}u\frac{\partial v}{\partial n}{\rm d}s-\int_{\Omega }\bigtriangledown u \bigtriangledown v{\rm d}x \end{align} $

    (11)

    本文的目标是设计控制器使系统达到稳定状态.将设计的中和控制器(7)作用于系统(2), 得到如下闭环系统.

    $ \begin{align} \begin{aligned} \frac{\partial z}{\partial t}=D_{d}\Delta z(x, t)+Hz(x, t)+Gz(x, t-\tau)\\ \end{aligned} \end{align} $

    (12)

    其中,

    $ \begin{align*} &z(x, t)=\begin{bmatrix} W(x, t)\\ \chi (x, t) \end{bmatrix}\in {\bf R}^{2n} \\&D_{d}=\begin{bmatrix} D_{a}I_{n}&0\\ 0&D_{c}I_{n} \end{bmatrix}\in {\bf R}^{2n\times 2n} \\&H=\begin{bmatrix} A_{0}&K\\ B_{c1}&B_{c} \end{bmatrix}\in {\bf R}^{2n\times 2n} \\&G=\begin{bmatrix} A&0\\ B_{c2}&0 \end{bmatrix}\in {\bf R}^{2n\times 2n} \end{align*} $

    定理1. 在给定参数条件下, 若存在矩阵$Y$, $W$, $S_{1n}$, $S_{2n}\in {\bf R}^{n \times n}$为正定对称矩阵, 矩阵$M$, $D$, $X_{1}$, $X_{2}$, $X_{3}$, $X_{4}\in {\bf R}^{n \times n}$满足以下矩阵不等式

    $ \begin{align} &\begin{bmatrix} \Pi_{1}&X_{3}& AS_{1n}& 0&W& I\\ *&\Pi_{2}&X_{4}&0&M^{{\rm T}}&0\\ *&*&-S_{1n}&0&0&0\\ *&*&*&-S_{2n}&0&0\\ *&*&*&*&-S_{1n}&0\\ *&*&*&*&*&-S_{2n}\\ \end{bmatrix} <0 \end{align} $

    (13)

    $ \begin{align} &\begin{bmatrix} W& I\\ I&Y\\ \end{bmatrix}>0, \quad \begin{bmatrix} S_{1n}&0\\ 0&S_{2n}\\ \end{bmatrix}>0 \end{align} $

    (14)

    其中

    $ \begin{align*}\Pi_{1}= &W^{{\rm T}}A_{0}^{{\rm T}}+A_{0}W+X_{2} +X_{2}^{{\rm T}} \\ \Pi_{2}= &A_{0}^{{\rm T}}Y+Y^{{\rm T}}A_{0}+X_{1}+X_{1}^{{\rm T}} \\ X_{1}= &DB_{c1}\\X_{2}= &KM^{{\rm T}}\\ X_{3}= &W^{{\rm T}}A_{0}^{{\rm T}}Y+X_{2}^{T}Y+W^{{\rm T}}X_{1}^{{\rm T}} +\\ & MB_{c}^{{\rm T}}D^{{\rm T}}+A_{0} \end{align*} $

    $ \begin{align*}X_{4}=&Y^{{\rm T}}AS_{1n}+DB_{c2}S_{1n} \end{align*} $

    则构建的中和控制器使系统镇定.

    证明.构造李雅普诺夫函数为

    $ \begin{align} &V(t, z(x, t))=\int_{\Omega }z^{{\rm T}}(x, t)Pz(x, t){\rm d}x +\nonumber\\ &\qquad \int_{\Omega }\int_{t-\tau}^{t}z^{{\rm T}}(x, \theta )Sz(x, \theta ){\rm d}\theta {\rm d}x \end{align} $

    (15)

    $ \begin{align} &\dot{V}(t, z(x, t))=\int_{\Omega }z^{{\rm T}}(x, t)P\dot{z}(x, t) +\notag\\&\qquad z^{{\rm T}}(x, t)P\dot{z}(x, t){\rm d}x+ \int_{\Omega }z^{{\rm T}}(x, t)S z(x, t) -\notag\\&\qquad z^{{\rm T}}(x, t-\tau)Sz(x, t-\tau){\rm d}x \notag\\[2mm] &\dot{V}(t, z(x, t))=\int_{\Omega }z^{{\rm T}}(x, t)(PH+H^{{\rm T}}P)z(x, t) +\notag\\ &\qquad z^{{\rm T}}(x, t)PG z(x, t-\tau) +\notag\\ &\qquad 2z^{{\rm T}}(x, t)PD_{d}\Delta z(x, t) +\notag\\ &\qquad z^{{\rm T}}(x, t-\tau)G^{{\rm T}}P z(x, t) +\notag \\&\qquad z^{{\rm T}}(x, t)Sz(x, t)- z^{{\rm T}}(x, t-\tau)Sz(x, t-\tau){\rm d}x\notag \end{align} $

    令$I(t)=\int_{\Omega }z^{{\rm T}}(x, t)\Gamma \Delta z(x, t){\rm d}x$, $\Gamma=(\Gamma_{ij})=2PD_{d}$, 利用引理2和边界条件

    $ \begin{align*} I(t)= &\sum\limits_{i=1}^{2n}\sum\limits_{j=1}^{2n} \Gamma_{ij}\int_{\Omega }z_{i}(x, t)\Delta z_{i}(x, t){\rm d}x =\\ &\sum\limits_{i=1}^{2n}\sum\limits_{j=1}^{2n} \Gamma_{ij}[\int_{\Theta }z_{i}(x, t)\frac{\partial z_{i}(x, t) }{\partial n}{\rm d}x -\\ &\int_{\Omega }\bigtriangledown z_{i}(x, t) \bigtriangledown z_{i}(x, t){\rm d}x]=\\ &-\int_{\Omega }\bigtriangledown z^{{\rm T}}(x, t)\Gamma \bigtriangledown z(x, t){\rm d}x\end{align*} $

    $ \begin{align*} &\dot{V}(t, z(x, t))=\int_{\Omega }z^{{\rm T}}(x, t)(PH+H^{{\rm T}}P)z(x, t) +\\ &\qquad z^{{\rm T}}(x, t)PG z(x, t-\tau) -\\ &\qquad \bigtriangledown z^{{\rm T}}(x, t)2PD_{d} \bigtriangledown z(x, t) +\\ &\qquad z^{{\rm T}}(x, t-\tau)G^{{\rm T}}P z(x, t) +\\ &\qquad z^{{\rm T}}(x, t)Sz(x, t)- z^{{\rm T}}(x, t-\tau)Sz(x, t-\tau){\rm d}x \\ &\dot{V}(t, z(x, t))\leq\int_{\Omega } \begin{bmatrix} z(x, t) \\ z(x, t-\tau) \end{bmatrix}^{{\rm T}}\times\\&\qquad \begin{bmatrix} PH+H^{{\rm T}}P+S& PG\\ G^{{\rm T}}P&-S \end{bmatrix} \begin{bmatrix} z(x, t) \\ z(x, t-\tau) \end{bmatrix}{\rm d}x \end{align*} $

    所以若存在不等式

    $ \begin{align} \begin{bmatrix} PH+H^{{\rm T}}P+S& PG\\ G^{{\rm T}}P&-S \end{bmatrix} <0 \end{align} $

    (16)

    成立, 则

    $ \begin{align} \dot{V}(t, z(x, t)) <0 \end{align} $

    (17)

    由矩阵不等式(16)求解控制器时, 控制器参数$K$, $B_{c}$, $B_{c1}$, $B_{c2}$包含在矩阵$H$和矩阵$G$中不能被求出, 所以接下来的目的是求解控制器参数.

    首先, 将矩阵$P$, $S$分区, 定义:

    $ \begin{align} &P= \left[ \begin{array}{ccc} Y&D\\ D^{{\rm T}}&T_{c}\\ \end{array} \right]\notag\\[1mm] &P^{-1}= \left[ \begin{array}{ccc} W&M\\ M^{{\rm T}}&N\\ \end{array} \right]\notag\\[1mm] &S^{-1}= \left[ \begin{array}{ccc} S_{1n}&0\\ 0&S_{2n}\\ \end{array} \right] \end{align} $

    (18)

    其中, $Y$, $W$是正定对称矩阵, $D$, 是可逆矩阵, 且$PP^{-1}=I$

    $ \begin{align} MD^{{\rm T}}=I-WY \end{align} $

    (19)

    定义两个矩阵

    $ \begin{align} \begin{array}{cc} E= \left[ \begin{array}{cc} W&I\\ M^{{\rm T}}&0\\ \end{array} \right], ~~~ F= \left[ \begin{array}{cc} I&Y\\ 0&D^{{\rm T}}\\ \end{array} \right]\\ \end{array} \end{align} $

    (20)

    $ \begin{align} %20式 PE=F, ~~~ E^{{\rm T}}PE=E^{{\rm T}}F= \left[%第一个矩阵 \begin{array}{cc} W&I\\ I&Y\\ \end{array} \right]>0 \end{align} $

    (21)

    对式(16)两边左乘, 右乘, 有

    $ \begin{align} %22式 \begin{bmatrix} E^{{\rm T}}H^{{\rm T}}F+F^{{\rm T}}HE+E^{{\rm T}}SE&F^{{\rm T}}GS^{-1}\\ S^{-1}G^{{\rm T}}F&-S^{-1} \end{bmatrix} <0 \end{align} $

    (22)

    由式(18) $\sim$ (22), 根据Schur补定理及所给条件可得定理成立.

    推论1.对于系统(1), 若取时滞$\tau=0$, 则系统变为

    $ \begin{align*} %2式 \frac{\partial W(x, t)}{\partial t}= &D_{a}\Delta W(x, t) +\\&(A_{0}+A) W(x, t)+U(x, t) \end{align*} $

    则可构造控制器为

    $ \begin{align*} \begin{cases} \dfrac{\partial \chi (x, t)}{\partial t}=D_{c}\Delta \chi(x, t) +\\\qquad\qquad \quad B_{c} \chi(x, t)+ B_{c1} W(x, t)\\[1mm] U(x, t)=K\chi(x, t) \end{cases} \end{align*} $

    在给定参数条件下, 若存在矩阵$Y$, 为正定对称矩阵, 矩阵$M$, $D$, $X_{1}$, $X_{2}$, 满足以下矩阵不等式

    $ \begin{align*} %12式 &\begin{bmatrix} \Phi_{1}&X_{3}\\ *&\Phi_{2}\\ \end{bmatrix} <0 \\[1mm]& %13式 \begin{bmatrix} W& I\\ I&Y\\ \end{bmatrix}>0 \end{align*} $

    其中

    $ \begin{align*}\Phi_{1}= &W^{{\rm T}}(A_{0}^{{\rm T}}+A^{{\rm T}})+(A_{0}+A)W+X_{2}+X_{2}^{{\rm T}}\\[1mm] \Phi_{2}= &(A_{0}^{{\rm T}}+A^{{\rm T}})Y+Y^{{\rm T}}(A_{0}+A)+X_{1}+X_{1}^{{\rm T}}\\[1mm] X_{1}= &DB_{c1}\\[1mm] X_{2}= &KM^{{\rm T}}\\[1mm] X_{3}= &W^{{\rm T}}(A_{0}^{{\rm T}}+A^{{\rm T}})Y+X_{2}^{{\rm T}}Y+W^{{\rm T}}X_{1}^{{\rm T}} + \\&MB_{c}^{{\rm T}}D^{{\rm T}}+A_{0}+A\end{align*} $

    则构建的中和控制器使系统镇定.

    证明. 构造李雅普诺夫函数为

    $ \begin{align*} %11式 V(t, z(x, t))=\int_{\Omega }z^{{\rm T}}(x, t)Pz(x, t){\rm d}x \end{align*} $

    后续参考定理1的证明.

    注1.本文针对一类具有时空特性的分布参数系统提出了中和控制方法.这类系统的控制许多学者进行了研究, 主要是用分布式控制设计控制器, 例如文献[8-11]等.我们采用完全不同的设计思想, 在相应的时空上设计控制器, 控制器形如式(6)和式(7), 利用设计的分布参数系统控制目标分布参数系统, 控制过程为设计的分布参数系统自行运动的过程.当知道控制系统的模型时, 这种控制方法在实际应用中容易实现, 此控制方法最大优点在于不用对状态点进行直接控制.

    为了说明问题, 考虑如下分布参数系统及控制系统

    $ \begin{align*} \begin{cases} \dfrac{\partial W(x, t)}{\partial t}=D_{a}\Delta W(x, t)+A_{0} W(x, t) +\\ \qquad\qquad \quad\ A W(x, t-\tau)+U(x, t)\\[1mm] \dfrac{\partial \chi(x, t)}{\partial t}=D_{c}\Delta \chi(x, t)+B_{c} \chi(x, t) +\\ \qquad\qquad \quad B_{c1} W(x, t)+B_{c2} W(x, t-\tau)\\[1mm] U(x, t)=K\chi (x, t) \end{cases} \end{align*} $

    对分布参数系统, $m=1$, $n=2$, 系统参数$D_{a}$ $=$ $1.9$, , ; 取时滞$\tau$ $=$ $1$, 应用定理1提出的方法, $D_{c}=1.8$, 通过MAT- LAB软件中的LMI工具箱, 可以得到控制系统参数.

    $ \begin{align*}&K=10^{4}\times\begin{bmatrix} 2.4122&0\\ 0&2.3778 \end{bmatrix}\\[1mm]&B_{c}=10^{4}\times\begin{bmatrix} -1.6755&-0.6569\\ -0.6664&-1.2168 \end{bmatrix}\\[1mm]&B_{c1}=\begin{bmatrix} -1.5127&0.3684\\ 0.3684&-1.2435 \end{bmatrix} \\[1mm]&B_{c2}=\begin{bmatrix} -0.6946&-0.5525\\ -0.2763&-1.0235 \end{bmatrix}\end{align*} $

    给定系统的初始条件

    $ \begin{align*} &w_{1}(x, 0)=-\exp(0.7 (-x+5))\\[1mm] &w_{2}(x, 0)=0.5\sin(2x)\\[1mm] &\eta_{1}(x, 0)=0\\[1mm] &\eta_{2}(x, 0)=0\end{align*} $

    图 1图 2是系统状态和控制系统的状态图.从图 1图 2可以看出, 系统在控制器的作用下, 经过一段时间后可达到稳定状态.

    图 1  分布参数系统$W(x, t)$状态图
    Fig. 1  The state of distributed parameter systems $W(x, t)$
    图 2  控制系统$\chi (x, t)$状态图
    Fig. 2  The state of controller $\chi(x, t)$

    本文首先给出中和控制器设计思路, 针对具有时滞特性的分布参数系统, 设计中和控制器, 讨论此类系统的稳定问题.利用Lyapunov稳定性理论并结合LMI处理方法, 得出了具有时滞特性分布参数系统稳定中和控制器存在的充分条件.最后结合所给条件, 给出一个数值仿真说明其有效性.


  • 本文责任编委 姚鹏飞
  • 图  1  分布参数系统$W(x, t)$状态图

    Fig.  1  The state of distributed parameter systems $W(x, t)$

    图  2  控制系统$\chi (x, t)$状态图

    Fig.  2  The state of controller $\chi(x, t)$

  • [1] Ray W H. Advanced Process Control. New York: McGraw-Hill, 1981.
    [2] Christofides P D. Nonlinear and Robust Control of PDE Systems. Boston: Birkhäuser Boston, 2001.
    [3] Deng H, Li H X, Chen G R. Spectral-approximation-based intelligent modeling for distributed thermal processes. IEEE Transactions on Control Systems Technology, 2005, 13:686-700 doi: 10.1109/TCST.2005.847329
    [4] Padhi R, Ali S F. An account of chronological developments in control of distributed parameter systems. Annual Reviews in Control, 2009, 33:59-68 doi: 10.1016/j.arcontrol.2009.01.003
    [5] Baillieul J. Linearized models for the control of rotating beams. In: Proceedings of the 27th Conference on Decision and Control. Austin, Texas, USA: IEEE, 1989, 3: 1726-1731
    [6] Najar F, Choura S, Abdel-Rahman E M, El-Borgi S, Nayfeh A. Dynamic analysis of variable-geometry electrostatic microactuators. Journal of Micromechanics and Microengineering, 2006, 16(11):2449-2457 doi: 10.1088/0960-1317/16/11/028
    [7] Ly C, Doiron B. Correction:divisive gain modulation with dynamic stimuli in integrate-and-fire neurons. PLoS Computational Biology, 2009, 5(6), DOI: 10.1371/annota-tion/9c4474f6-e58e-41ea-af4b-090bbb86f474
    [8] Luo Y P, Xia W H, Liu G R, Deng F Q. LMI approach to exponential stabilization of distributed parameter control systems with delay. Acta Automatica Sinica, 2009, 35(3):299-304 http://www.ams.org/mathscinet-getitem?mr=2535951
    [9] 高存臣, 刘振, 徐瑞萍.一类具有连续分布时滞的分布参数系统的反馈控制.控制与决策, 2013, 28(3):445-450 http://d.old.wanfangdata.com.cn/Periodical/kzyjc201303022

    Gao Cun-Chen, Liu Zhen, Xu Rui-Ping. Feedback control for a class of distributed Parameter systems with continuous distributed time-delay. Control and Decision, 2013, 28(3):445-450 http://d.old.wanfangdata.com.cn/Periodical/kzyjc201303022
    [10] Xing H L, Li D H, Gao C C, Kao Y G. Delay-independent sliding mode control for a class of quasi-linear parabolic distributed parameter systems with time-varying delay. Journal of the Franklin Institute, 2013, 350(2):397-418 doi: 10.1016/j.jfranklin.2012.12.007
    [11] Yuan K, Alofi A, Cao J D, Al-Mazrooei A, Elaiw A. Synchronization of the coupled distributed parameter system with time delay via proportional-spatial derivative control. Discrete Dynamics in Nature and Society, 2014, 2014: Article No. 418258
    [12] Babaei D, Armaou P A. Adaptive control of chemical distributed parameter systems. IFAC-PapersOnLine, 2015, 48(8):681-686 doi: 10.1016/j.ifacol.2015.09.047
    [13] Demetriou M A. Boundary adaptive synchronization of networked PDEs with adaptive parameter estimators. IFAC-PapersOnLine, 2016, 49(8):242-247 doi: 10.1016/j.ifacol.2016.07.448
    [14] Jiang Z X, Cui B T, Wu W, Zhuang B. Event-driven observer-based control for distributed parameter systems using mobile sensor and actuator. Computers and Mathematics with Applications, 2016, 72(12):2854-2864 doi: 10.1016/j.camwa.2016.10.009
    [15] Wang Z P, Wu H N. Finite dimensional guaranteed cost sampled-data fuzzy control for a class of nonlinear distributed parameter systems. Information Sciences, 2016, 327:21-39 doi: 10.1016/j.ins.2015.08.009
    [16] 郭涛, 梁燕军.不确定非线性时滞关联大系统自适应分散容错控制.自动化学报, 2017, 43(3):486-492 http://www.aas.net.cn/CN/abstract/abstract19027.shtml

    Guo Tao, Liang Yan-Jun. Adaptive decentralized fault-tolerant control for uncertain nonlinear time-delay large scale systems. Acta Automatica Sinica, 2017, 43(3):486-492 http://www.aas.net.cn/CN/abstract/abstract19027.shtml
    [17] Chen W H, Zhong J C, Zheng W X. Delay-independent stabilization of a class of time-delay systems via periodically intermittent control. Automatica, 2016, 71:89-97 doi: 10.1016/j.automatica.2016.04.031
    [18] Xie T T, Liao X F, Li H Q. Leader-following consensus in second-order multi-agent systems with input time delay:an event-triggered sampling approach. Neurocomputing, 2016, 177:130-135 doi: 10.1016/j.neucom.2015.11.013
    [19] Boyd S, El Ghaoui L, Feron E, Balakrishnan V. Linear Matrix Inequalities in System and Control Theory, Philadelphia, PA, USA: Society for Industrial and Applied Mathematics, 1994.
    [20] 崔宝同, 楼旭阳.时滞分布参数系统理论及其应用.北京:国防工业出版社, 2009. 8-12

    Cui Bao-Tong, Lou Xu-Yang. Theory and Application of Time-Delay Distributed Parameter System. Beijing:National Defense Industry Press, 2009. 8-12
  • 期刊类型引用(9)

    1. 李力厚,李杰,王书义. 基于间歇单点测量的PDE系统观测器设计及应用. 控制工程. 2023(03): 536-542 . 百度学术
    2. 李成群,李延波,张晶华,苏志伟. 混合时滞随机中立型分布参数系统的稳定性. 数学的实践与认识. 2022(08): 186-192 . 百度学术
    3. 周笔锋,罗毅平,唐果宁. 分布参数系统源控制系统设计. 自动化学报. 2022(12): 3062-3066 . 本站查看
    4. 李延波,陈超洋,张晶华,李成群. 分布时滞随机偏微分系统的均方指数稳定性. 控制理论与应用. 2022(11): 2185-2192 . 百度学术
    5. 刘青松. 基于嵌套-伪预估器反馈的时滞控制系统输入时滞补偿. 自动化学报. 2021(10): 2464-2471 . 本站查看
    6. 周笔锋,罗毅平. 时滞分布参数系统间歇非周期中和控制器设计. 计算机工程与应用. 2020(08): 220-225 . 百度学术
    7. 吴鸿娟,熊江,冯玉明,涂正文. 一种带时间窗口的双脉冲控制器下的忆阻振荡系统滞同步. 自动化学报. 2020(07): 1507-1516 . 本站查看
    8. 潘月斗,王国防. 基于中立型系统理论的感应电机高精度磁链观测器研究. 自动化学报. 2020(10): 2109-2120 . 本站查看
    9. 刘秀翀,王占山. 系统H_∞范数计算:Lyapunov函数的直接优化方法. 自动化学报. 2019(08): 1606-1610 . 本站查看

    其他类型引用(5)

  • 加载中
  • 图(2)
    计量
    • 文章访问数:  1802
    • HTML全文浏览量:  241
    • PDF下载量:  587
    • 被引次数: 14
    出版历程
    • 收稿日期:  2017-02-17
    • 录用日期:  2017-05-26
    • 刊出日期:  2018-12-20

    目录

    /

    返回文章
    返回