时滞不确定δ算子控制系统的鲁棒非脆弱控制

  • 格式:pdf
  • 大小:195.50 KB
  • 文档页数:4

2010年7月 第33卷第4期 四川师范大学学报(自然科学版) Journal of Sichuan Normal University(Natural Science) July,2010 V01.33.No.4 

时滞不确定I!;算子控制系统的鲁棒非脆弱控制 

汤红吉, 韩彦武, 张凤然, 于 霞 

(南通大学理学院,江苏南通226007) 

摘要:考虑了一类时滞不确定8算子控制系统的鲁棒非脆弱控制问题.基于Lyapunov稳定性理论和线 

性矩阵不等式的方法,得出了系统存在非脆弱控制器的充分条件,并以线性矩阵不等式的形式给出.最后, 

数值算例表明了方法的有效性. 

关键词:6算子控制系统;非脆弱控制;时滞相关;线性矩阵不等式 

中图分类号:0231 文献标志码:A 文章编号:1001—8395(2010)04—0483—04 

doi:10.3969/j.issn.1001—8395.2010.04.016 

时滞和不确定性是系统运行中常见的现象,也 

是影响系统性能的主要因素.目前关于时滞不确定 

系统的鲁棒控制已经有了很多的结果-】 .在系统 

的运行中,控制器的参数也会由于种种原因发生变 

化,从而可能会导致系统失稳或者性能的下降,因 

此研究被控对象和控制器同时具有不确定性的控 

制,即非脆弱控制具有重要意义.关于时滞不确定 

系统的非脆弱控制已有若干结果 -10]. 算子控制 

系统作为连续系统和离散系统的桥梁,可以把连续 

系统和离散系统统一在 算子控制系统的框架 

内_1 .目前,关于时滞不确定8算子控制系统的研 

究结果较少 13-17].本文基于LMI的方法,研究了一 

类时滞不确定 算子控制系统的非脆弱控制问题, 

得出了时滞相关的非脆弱控制器存在的充分条件. 

1 问题描述 

考虑如下同时具有时滞和不确定性的6算子控 

制系统 

8x(t)=(A+△A( )) (f)+(Ad+ 

AAd( )) (t—nT)+Bu(£), (1) 

其中, (f)∈R 是状态变量,U( )∈R 是控制输 

入变量,nT是有界的时滞量,12,是正整数, 是采样 

周期且满足nT≤h.△A(t),△A (t)是具有适当维 

数的时变矩阵并具有如下定义 

[△A( )AAd(t)]=D1 1(t)[El E2], 

(t)=F (t)[,一G F (t) , (2) 

收稿日期:2008—09—15 基金项目:国家自然科学基金(60874021)资助项目 作者简介:汤红吉(1977一),女,博士,主要从事鲁棒控制的研究 其中,D E E 、G 是具有适当维数的常数矩阵, 

F。(t)是未知的时变矩阵满足 (t)F (t)≤ 假 

设矩阵[j—G F (t)] 是可逆矩阵,其中矩阵G 

满足I—G G1>0. 

下面给出本文中需要用到的基本引理. 

引理1_1 F是对称矩阵,对于具有适当维数 

的矩阵D、 ,则 

J1+DF(t)E+[DP(f)E] ≤0, 

其中, (£)和(1)式中 (t)的定义相同,当且仅 

当存在 >0使得 

川s 圳

一 】~

~ ̄-IE] 

考虑到状态反馈控制器本身含有不确定性 

(£)=( +△ (t)) (t), (3) 

其中, 是控制器增益矩阵,△ ( )=D2 ( )E3是控 

制器增益摄动, (t)=F2(t)[I—G2F2(t) , 

(£)F2( )≤,,G2满足,一 G2>0,使得闭环系统 

8x(t)=A (t)x(t)+Ad(t)x(t—nT),(4) 

对所有的不确定性渐近稳定,其中A (t)=A+ 

( +AK(t))+△A(t),Ad(t)=Ad+△Ad(t). 

2 主要结果 

定理1 考虑时滞不确定 算子的闭环控制系 

统(4), 、b是给定的常数,如果存在正定对称矩阵 

、c,、 ,矩阵l,,纯量OL>0,/3>0使得如下的线性 

矩阵不等式成立 四川师范大学学报(自然科学版) 33卷 

= lJ l2 

22 

术 木 

:f: 木 船 

23 0 

, 

一 

 ̄wj(4)式是渐近稳定的,相应的非脆弱控制器增益 

为K=YX~,其中 

11=(A +By)+(ax+Br) +U一 V, 

l2:一a(XA 一yTB ), 

13:AdX+v—bYtB 一6 , 

吐 =XT+nVT+2oLX,哇 =bX— Ad , 

3:一 一 一bAdX一6 T. 

证明构造如下的Lyapunov—Krasovskii泛函 

V(t)=V1(t)+v2(t)+ (t), 

其中 

(t)= (t)ex(t), 

v2(t)= ∑ (£一it)(2x(t— ), 

v3(t)=∑∑ (£-jr)re(t-jr), 

(_,)= (-,)一 (_,+ ), 

P、Q、R是对称正定矩阵,则 

6V1(t)=2x (t)P6x(t)+r(sx(t)) P6x(t)= 

1 0 I 船 

一aoLD1 0 

~6 0 

G1 0 

一 O 

一8l 8BD2 

一 8BD2 

<0, (5) 

(t)e(a (t)x(t)+ 

Ad(t)x(£一nT))+r(sx(t)) P6x(t), 

6 (t)= (t)9.x(£)一 (t—n )Q (t—nT), 1 n i 6 (f)= 1∑∑ (£一(卜1)r)re(t一 

‘ =】,=j U—1) )一e'r(t—iT)re(t一 )]≤ re(t)一 

P(t—it)l R[ P( )] 

r(缸(£))Try(t)一 [ (£一n )一 

(t)] [ (t—nT)一 (t)]. 

对于任意的具有适当维数的矩阵M、Ⅳ有 

2[g ̄x(t)+Nx(t—nT)] Esx(t)一A (t)x(t)一 

Ad(t)x(t—nT)]=0. (6) 

综上可得 

6V(t)=6V1(t)+6 (t)+ (t)≤ 

(t) 。 (t), (7) 其中 

(t)=[ (t)(瓤(t)) (t—nT)】 , 

A (£)P+PA (£)+Q一 R —A (£) PA (£)+ R—A (£)Ⅳ 

P+nTR+^ +M N—MtAd(t) 

一Q一 一NTA (£)一A (z)Ⅳ n』 

根据Lyapunov稳定性理论知,如果 <0,则 

闭环系统(4)渐近稳定. 

下证 <0.注意到 可以表示为如下形式 

=F+ l (f)豆,+[ 。 ( ) ] + 

( )豆:+[ ( ) :] , 

= A 

Ⅳ , [: 

=(A+ ) P+P(A+ )+Q一 R, 

3=PAd+R/nT一(A+BK) N, 

1"22=TP+nTR+MT+M, r3,=一a一 R—NtA 一A , D 0 o

 第4期 汤红吉等:时滞不确定艿算子控制系统的鲁棒非脆弱控制 

E1=[E1 0 E2],E2=[E3 0 0】. 

根据引理1,不等式 <0的充要条件是存在实数 

1>0, 2>0使得 

口 =f+ [一 一 】一aT+ 

T一 aT 【一G, j‘J , 

其中 

= (1,2) 

(2,2) 

球 

出 (1,3) 

N—MTAd 

(3,3) 

采 

术 

(1,1)=F11, (1,2)=一(A+BK) M, (1,3)= ,, (2,2)= :, (3,3)= ,. 在矩阵 的左右两边分别同时乘以矩阵 

diag{P~,P~,P~, J, j, J, J},并且令 

6"1~= , = ,X=P~,Y=KX,M=aP,N= 

,U=P Qe一,V=P P~,可得矩阵 又由 定理条件知 <0,从而可得 <0.证毕. 

注1 在定理的证明过程中,通过“牛顿一莱 

布尼兹”公式引入了矩阵变量M、 但是为了得出 

线性矩阵不等式的结果,人为的令“M=aP、N= 

6P”,引入了自由参数a,b.与文献[12]相比,增强 

了线性矩阵不等式的可解性.另一方面, 、Ⅳ本是 

自由矩阵,定义了“M=aP、N=bP”以后,限制了 

矩阵 、Ⅳ的形式,在一定程度上使所得的结果具 

有保守性. 

注2 如果控制器增益不存在摄动,即 (t) 

=0,取a=一1,b=0,则定理1就是文献[12]的主 

要结果;如果系统(1)中不存在不确定性,且控制 

器增不存在摄动,即 (t)=0, (t)=0,取。= 

一1,6=0,则定理1退化为文献[12]中的推论. 

3 数值算例 

考虑时滞不确定 算子控制系统(1),其中 

r Q 500 1 0 lq 厂 Q 200 2 0 0q 

A=I—Q 500 1 0 0 l,Ad=I—Q 200 2 0 0 l, 

Q 500 1 3 0a L Q 200 2 0 01 : 

[二 

1E PD1 

0 —61 MTD1 1 一s Ⅳ D1 

一, G1 

一I 

;f= 术 

木 术 82E PBD2 

0 一s MTBD2 

0 一 N BD2 

0 0 

0 0 

一J G2 

一I 

=r 

厂Q 5 0 0] r 0.500 1 0 0q Gl=1 0 Q 5 0 l,El=I—Q 500 1 0 0 I, 

0 0 Q 5a L Q 500 1 0 0 

r Q 200 2 0 0q r 0.106 2 0 0q 

=I—Q 200 2 0 0 l, =I—Q 106 2 0 0 I, 

Q 200 2 0 0J L—Q 106 2 0 0 

G2=G1,D2=[0.5 0.5 0.5]. 

取a=一3,b:0.1,n=1 000,T=0.O1,利用 

MATLAB线性矩阵不等式(LMI)工具箱求解线性 

矩阵不等式(LMI)(5),可得相应的鲁棒非脆弱控 

制器增益为 

K=[0.951 7 —0.028 3 0.703 9]. 

本文研究了一类时滞不确定算子控制系统的 

鲁棒非脆弱控制问题,针对控制器增益存在加性扰 

动的情形,以LMI的形式给出了系统存在非脆弱状 

态反馈控制器的充分条件,最后,数值算例表明本 

文方法的可行性. 

致谢感谢南通市应用研究计划(K2009035)、南 

通大学自然科学&,/ ̄(oszoo9)和南通大学引进人才项 

目& ̄(o3o8oo55)对本文的资助.