D0I:10.13374/1.issnl00103.2009.05.022 第31卷第5期 北京科技大学学报 Vol.31 No.5 2009年5月 Journal of University of Science and Technology Beijing May 2009 基于LMⅦ的鲁棒H∞容错动态输出反馈控制设计 陈 明 童朝南 北京科技大学信息工程学院，北京100083 摘要针对线性不确定系统，研究了其执行器失效情况下鲁棒容错H动态输出反馈控制设计问题.基于连续型执行器故 障模型，利用线性矩阵不等式(LM)方法提出了线性不确定系统动态输出反馈H©容错控制器存在的充分条件，给出了动态 输出反馈H控制器的设计方法。所获得的控制器不仅能使故障系统鲁棒稳定，并且能达到给定的H©性能指标·仿真实例 证明了所提出设计方法的有效性， 关键词线性矩阵不等式；容错控制；输出反馈；H©控制：鲁棒性 分类号TP273 Design of robust Ho fault-tolerant output feedback control based on LMI CHEN Ming.TONG Chao-nan School of Information and Engineering.University of Science and Technology Beijing.Beijing 100083,China ABSTRACT For linear uncertain systems with actuator faults,the problem of robust fault-tolerant He dynamic output feedback control was studied.Based on an actuator failure model with continuous gain.the sufficient condition for dynamic output feedback H robust fault-tolerant controllers was proposed by using a linear matrix inequality (LMI)approach.and the corresponding design method of a dynamic output feedback Ho controller was given.The proposed controller not only enables the system to keep robust stabilization but also achieves the given H performance.A simulation experiment shows that the proposed method is effective. KEY WORDS linear matrix inequality:fault-tolerant control:output feedback:Hoo control:robustness 容错控制是系统对故障的容忍技术，在设计控 控制设计的线性矩阵不等式(linear matrix inequality, 制器时，无论闭环系统控制部件是否发生故障都能 LM)方法，并把取得的理论成果应用到一个实际 保证系统稳定运行，在实际系统中，由于模型不确 F一16歼击机模型上，取得了较好的控制效果. 定性及外界扰动的广泛存在性，使得鲁棒容错控制 基于上述文献的研究成果，本文提出了一种鲁 的研究更具有实际意义]. 棒容错动态输出反馈控制的LMI设计方法，与其 近年来，关于不确定系统的鲁棒容错控制设计 他文献研究的不同之处在于：基于连续型执行器故 的研究成果很多3]，文献[5]针对线性不确定系 障模型，利用有界实引理和线性矩阵不等式，通过引 统，研究其在传感器故障情况下的鲁棒H©可靠输 入变量代换，把动态输出反馈控制器的求解条件转 出反馈控制问题，通过求解两个Riccati方程获得控 换为标准的LM1,仿真结果表明，得到的动态输出 制器参数.文献[6]研究了范数有界参数时变不确 反馈控制器不仅能使闭环系统在容许的执行器故障 定系统的鲁棒可靠H©控制问题，把执行器故障作为 模式下具有鲁棒性及容错性，并且满足预先给定的 系统的随机扰动信号，设计鲁棒可靠状态反馈控制 H∞干扰衰减指标. 器；该控制器不仅保证系统对所有容许的参数不确定 1问题描述 及执行器故障是稳定的，并且能达到给定的H∞性能 指标，文献[7]给出了连续系统的可靠鲁棒飞行跟踪 考虑如下形式的线性不确定系统： 收稿日期：2008-03-27 基金项目：北京市教委重点学科共建项目(N。,XK100080537) 作者简介：陈明(1977-)，女，博士研究生，E-mail:cm8061@ima·com:童朝南(1955一)，男，教授，博士生导师
基于 LMI 的鲁棒 H∞容错动态输出反馈控制设计 陈 明 童朝南 北京科技大学信息工程学院北京100083 摘 要 针对线性不确定系统研究了其执行器失效情况下鲁棒容错 H∞动态输出反馈控制设计问题．基于连续型执行器故 障模型利用线性矩阵不等式（LMI）方法提出了线性不确定系统动态输出反馈 H∞容错控制器存在的充分条件给出了动态 输出反馈 H∞控制器的设计方法．所获得的控制器不仅能使故障系统鲁棒稳定并且能达到给定的 H∞性能指标．仿真实例 证明了所提出设计方法的有效性． 关键词 线性矩阵不等式；容错控制；输出反馈；H∞控制；鲁棒性 分类号 TP273 Design of robust H∞ fault-tolerant output feedback control based on LMI CHEN MingT ONG Chao-nan School of Information and EngineeringUniversity of Science and Technology BeijingBeijing100083China ABSTRACT For linear uncertain systems with actuator faultsthe problem of robust fault-tolerant H∞ dynamic output feedback control was studied．Based on an actuator failure model with continuous gainthe sufficient condition for dynamic output feedback H∞ robust fault-tolerant controllers was proposed by using a linear matrix inequality （LMI） approachand the corresponding design method of a dynamic output feedback H∞ controller was given．T he proposed controller not only enables the system to keep robust stabilization but also achieves the given H∞ performance．A simulation experiment shows that the proposed method is effective． KEY WORDS linear matrix inequality；fault-tolerant control；output feedback；H∞ control；robustness 收稿日期：2008-03-27 基金项目：北京市教委重点学科共建项目（No．XK100080537） 作者简介：陈 明（1977—）女博士研究生E-mail：cm8061＠sina．com；童朝南（1955—）男教授博士生导师 容错控制是系统对故障的容忍技术．在设计控 制器时无论闭环系统控制部件是否发生故障都能 保证系统稳定运行．在实际系统中由于模型不确 定性及外界扰动的广泛存在性使得鲁棒容错控制 的研究更具有实际意义［1—2］． 近年来关于不确定系统的鲁棒容错控制设计 的研究成果很多［3—4］．文献［5］针对线性不确定系 统研究其在传感器故障情况下的鲁棒 H∞可靠输 出反馈控制问题通过求解两个 Riccati 方程获得控 制器参数．文献［6］研究了范数有界参数时变不确 定系统的鲁棒可靠 H∞控制问题把执行器故障作为 系统的随机扰动信号设计鲁棒可靠状态反馈控制 器；该控制器不仅保证系统对所有容许的参数不确定 及执行器故障是稳定的并且能达到给定的 H∞性能 指标．文献［7］给出了连续系统的可靠鲁棒飞行跟踪 控制设计的线性矩阵不等式（linear matrix inequality LMI）方法并把取得的理论成果应用到一个实际 F—16歼击机模型上取得了较好的控制效果． 基于上述文献的研究成果本文提出了一种鲁 棒容错动态输出反馈控制的 LMI 设计方法．与其 他文献研究的不同之处在于：基于连续型执行器故 障模型利用有界实引理和线性矩阵不等式通过引 入变量代换把动态输出反馈控制器的求解条件转 换为标准的 LMI．仿真结果表明得到的动态输出 反馈控制器不仅能使闭环系统在容许的执行器故障 模式下具有鲁棒性及容错性并且满足预先给定的 H∞干扰衰减指标． 1 问题描述 考虑如下形式的线性不确定系统： 第31卷 第5期 2009年 5月 北 京 科 技 大 学 学 报 Journal of University of Science and Technology Beijing Vol．31No．5 May2009 DOI:10．13374／j．issn1001－053x．2009．05．022
第5期 陈明等：基于LM的鲁棒H容错动态输出反馈控制设计 .649. (t)=(A+AA)x(t)+B1w(t)+ 其中，x(t)为控制器状态.由式(1)和式(3)构成的 (B2+△B2)uF(t) 闭环系统为： (1) z(t)=Cix(t)+D1u(t) xcl(t)=Acixel(t)+Beiw(t) (4) y(t)=C2x(t)+D21w(t) z(t)=Cexel 其中， 其中，x(t)∈R",u(t)∈Rm,z(t)∈RP,y(t)∈R9 xe(t)=[xT(t)x(t)]T, 分别为系统的状态向量、控制输入向量、控制输出向 量和量测输出向量；“(t)表示具有执行器故障的 「Acl Bel1 控制输入向量，满足uF(t)=Mu(t);w(t)∈R'能 量有限的外部干扰输入向量；A、B1、B2、C1、C2、 T「A+△A(B2+△B2)MCe1「 B D12和D21为适当维数的常值矩阵；△A、△B2分别表 L B.C2 LB.D21- 示参数的不确定性，且它们满足以下条件： [C1 D12MC] 0 定义1线性不确定系统(4)，对于系统参数不 [△A△B2]=L⊙(t)[E1E2] (2) 确定性及容许的执行器故障，给定一个常数Y>0, 其中，L、E1和E2是已知矩阵；O(t)是具有适当维 若存在动态输出反馈控制器K(s),使得下列条件 数的未知函数矩阵，其元素是Lebesgue可观测的， 满足： 且满足⊙(t)⊙(t)≤I. (1)闭环系统的所有极点均位于复平面的左半 本文采用文献[2]提出的连续型执行器故障模 平面， 式.该故障模型既包含离散故障模型，又表示离散 (2)从扰动输入w(t)到被控输出z(t)的传递 故障模型不能表示的故障情况.定义故障阵： 函数矩阵满足‖Gw(s)‖<Y,则称控制器K(s) M=diag(m1,m2,…,mm）. 为鲁棒H。容错输出反馈控制， 引理18]设E、F和△(t)是具有适当维数的 其中，0≤mh≤m;≤m(i=1,2,…,m).当m:=0 矩阵，且满足AT(t)△(t)≤I,那么对于任意实数 时，表示系统执行器完全失效；m:=1表示执行器正 a>0,有 常；0≤mh<m:<mi,m≥1表示执行器失效.引 EA(t)F+(EA(t)F)T<aEET+a-FTF. 入如下矩阵： 引理2)设R1、R2是具有适当维数的矩阵， Mo=diag(m01,mo2,…,m0m）, (t)为时变适维对角矩阵，且满足|(t)川≤U,U 为正定对角矩阵，则R:(t)R2+R(t)TRT<0 J=diag(j1,j2,…,jm）, 的充分必要条件是存在正实数B>0,使得 N=diag(n1,n2,…,nm）. BR1 URI+B1R2UR2<0. 式中， 引理310)（有界实引理）对于系统G(s)= moi=mw+mu 2 小 [AB,若矩阵A稳定，则使‖Gm(s)‖=<y的 LC:0J n=m二m0,1=1,2,,m. 充分必要条件是存在实称矩阵P>0,满足如下矩阵 moi 不等式： 由此得：M=Mo(I+N),|N≤J≤I. ATP+PA+y-2PBBTP+CTC<0. 本文目的是为故障系统(1)设计一个动态输出 2系统分析及主要结论 反馈控制器，使闭环系统在容许的执行器故障模式 下具有鲁棒性及容错性，并且使从扰动W(t)到被控输 定理1系统(4)，对于给定的任意常数Y>0, 出z(t)的传递函数矩阵G(s)满足‖Gm(s)‖m< 若存在正常数e1、e2以及对称正定矩阵X、Y和矩 Y,其中Y是预先给定的干扰衰减指标. 阵A、B、d,使得下列矩阵不等式组LMs存在可行 对于系统(1)，设其动态输出反馈控制器为： 解： sc(t)=Acxe(t)+Bcy(t) (3) (5) u(t)=C.x
第 5 期 陈 明 等 ： 基 于 L MI 的 鲁 棒 H∞ 容 错 动 态 输 出 反 馈 控 制 设 计 · 6 4 9 ·
.650, 北京科技大学学报 第31卷 元11 π21 π22 0 LTY -1。 0 π42 0 -e21I 米 r(X,Y,A,B,C,e1,e2)= BT 元52 0 0 -YI <0 π61 π62 0 0 0 元60 元71 π72 0 π7 π77 米 1 0 E1 0 0 Lπg1 0 0 0 0 0 0 0 -e2I (6) 则称故障系统(4)为鲁棒H。容错动态输出反馈控 制系统， xa-[v-w-v-tJ 其中， 由X1X=I得到： 1=AX+XAT+B2Mo+(B2MoC)T+ UVT=I-XY (8) 2eLLT+e2B2 MoJMO B7, 根据矩阵Schur补性质，Xa>0可以等价写成如下 形式： T21=AT+A+62 YB2 MoJMOBI, U-(XYX-X)U-T-VTY-1V>0 (9) 22=YA+ATY+BC2+(BC2)T, 把式(8)代入式(9)中，整理可以得到： π42=JiM6B5Y,π2=(YB1+BD21)T, U-(XYX-X)U-T-VTY-1V= U-[XYX-X-(I-XY)Y-1(I-YX)]U-T= 61=C1X+D12MoC+E2D12MoJMGB2, U-1(x-Y1)U-T>0, π62=C1+e2D12 MoJMo B7Y, 即X-Y-1>0,所以得到式(5). 66=-yI+E2D12MoJMO D12, 根据引理3，‖G2w（s)‖<y的充分必要条 =(E1+E2MoC)+2E2MoJMOB2, 件是存在一个正定对称矩阵X>0,满足如下不等 式： π2=e2E2MoJM0BiY,π76=e2E2MoJM6D, Ad Xe+XclAcl+y2XeBe BaXel+CaCel<0 =-261I+62E2MoJME o=J (10) 如果上述矩阵不等式组的可行解为X“、Y“、 也就是 X、B·和C,则故障系统(4)在H指标约束下 AdXe+XelAd< 的鲁棒容错动态输出反馈控制器各参数如下： -Y2XcBe BaXel-CaCel<0, Ac=V-1(X”-Y"AX"- 所以A日X+XA是负定的，即A1为负定的，系 统为渐近稳定的，将式(10)写成如下形式： Y'B2MoC'-B'C2X*)U-T (7) 「AXa+XeAa XelBel C B。=V-1B* BaXd -yI 0<0(11) C。=d"U-T Cel 0 -YI 证明设X具有如下形式： 令变换矩阵： vrU-l(xYx-x)U-T小 g-[=X-0] 02)
· 6 5 0 · 北 京 科 技 大 学 学 报 第 3 1 卷
第5期 陈明等：基于LMⅫ的鲁棒H~容错动态输出反馈控制设计 .651. 将式（11)两边分别左乘diag(T,0,0)和右乘 定义： diag(g1,0,0),有 TET(AaXa+XclAe)E1 EIXcBe ETCa] p11 里0= BeXelE1 0 <0 Ψ1= 21 p22 Ce31 0 -yl BT (YB+BD21)T] -YI (13) 将式(12)代入式(13)中，经推导得到： [CX+D2MC C]0-YI 「[中11+△11 (15) Lφ21+421中2+△2 其中 Ψ0= <0 [BT (YB1+VB.D21)T]-yI u1=AX+XAT+B2M+(B2MC)T L[C1X+D12MC.UT C]0-yI (14) 21=AT++YB2MC 其中 22=YA+ATY+BC2+(BC2)T (1=AX+XAT+B2MC.UT+(B2MC.UT)T 根据式(2)，里。可以改写成： 21=AT+YAX+VB.C2X+YB2 MC.UT+VA.UT 「L1 22=YA +ATY+VB.C2+(VB.C2)T 41=△AX+XAAT+△B2MCUT+(△B2MC.UT)T 里0=平1+ 0 o(t)[(E1+E2MC)E:00]+ 421=△AT+Y△AX+Y△B2MCUT T A22=Y△A+△ATY 引入如下矩阵变换： YL o(t)[(E1+E2Md)E100] <0 0 A=YAX+VB C2X+VA.UT 0 B=VBe 由引理1和矩阵Schur补性质可知，对任意满足 d=C.UT ⑧T(t)⊙(t)≤I,存在正常数e1使得： 中1+e1LLT 21+61YLLT 中2+e1 YLLTY Ψ0≤ Bi (YB1+BD21)T-yI 关 =里2 C1X+D12Md C1 0 -yI E1+E2Md E1 0 0 设 Ψ3= AX+XAT+B2Mo@+(B2MoC)T+eLLT AT+A+E1YLLT YA+ATY+BC2+(BC2)T+YLLTY BT CiX+D12MoC Ci 0 -yI E1+E2MoC E1 0 0 (16)
第 5 期 陈 明 等 ： 基 于 L MI 的 鲁 棒 H∞ 容 错 动 态 输 出 反 馈 控 制 设 计 · 6 5 1 ·
.652 北京科技大学学报 第31卷 其中A=A+YB2MC 进一步考虑执行器故障M=Mo(I+N),有 B2Mo B2Mo T YB2Mo YB2Mo Ψ2=里3+ 0 N[0000]+ 0 N[d0000] (17) D12Mo D12M0 LE2Mo E2Mo 根据引理2，有 +T+B2Mo0+(B2Mo0)T:LT+E282MoMG B2 T+t:1uT+t2B2MowG过 ATC2+C2)T1TV+:22MOMG2Y T (B1+021T CI+D12MoD2M C1*:2D12MoBY M+2D12MMD12 (E1+E2Mo)+E2E2MoM6 B2 E1+2E2MoMO B7Y 0 1252MoMD12112E2MONG E o 0 21 0 J[d0000] (18) 0 0 由矩阵Schur补性质，可将上式等价变为： +T+B2MO0+(B2Mo0TLT+E2B2MoMB2 T,i1uT+e22o6明 +ATY+C2)TTY+E22MMGB2Y 旺 (B+iDznT - G1Il2Mod+212MoG设 C+:2D12MBY 1-212063 (E1+E2Mo+M E1+2F2MoMGBZY 0 226E-1+22MM65 0 0 -2 (19) 利用矩阵Schur补性质可以把式(I9)等效变换成式 3 仿真实例验证 (6),定理1获证. 根据定理1，给出不确定系统(4)在执行器故障 某钢铁集团热连轧活套系统的标称系统数学模 下的鲁棒容错动态输出反馈控制器的设计步骤： 型为： (1)给定H。干扰指标Y>0,选取适当正常数 (t)=Ax(t)+B2u(t) e1,e2,利用LMI工具箱中的Feasp求解矩阵不等式 y(t)=C2x(t) (5)和(6)； (2)若上述不等式组有解，则从任意一个可行 名 解出发，通过SVD技术，对矩阵UVT=I-XY进 x= 行奇异值分解，确定两个满秩矩阵U、V； (3)利用式(7)求得动态输出反馈控制器各参 数
· 6 5 2 · 北 京 科 技 大 学 学 报 第 3 1 卷