乔 梁, 李琳琳, 任俊超
(1. 北京科技大学 自动化学院, 北京 100083; 2. 东北大学 理学院, 辽宁 沈阳 110819)
奇异系统是由微分方程和代数方程共同构成,既能够刻画系统的动态行为,又能够刻画系统的静态行为,区别于正常线性系统,奇异系统具有的脉冲行为,以及对初始状态的不相容性,使得有关研究变得复杂而富有新颖性.相比于正常系统,奇异系统是更具有广泛形式的一类动力系统,存在于社会生产的众多领域中,例如机器人系统、电力系统、经济系统、宇航系统、电子网络系统等,故对奇异系统进行研究具有十分重要的理论意义和实用价值[1-3].容许性(正则、无脉冲、稳定)作为奇异系统正常运行的前提,是表征系统运动行为的一类重要结构特征.目前,针对奇异系统容许性的分析与综合问题研究受到了控制学界的广泛关注,并且得到了许多优秀的理论结果[4-8].文献[4]提出一种基于干扰观测器的鲁棒控制器,并利用李雅普诺夫稳定性理论分析了闭环奇异系统的渐近稳定性.文献[5]研究一类不完全转移率的Markov奇异系统的H∞控制问题,给出Markov跳变奇异系统的新型有界实引理以及状态反馈控制方法设计.文献[6]研究时滞矩形奇异系统的镇定问题,通过引入时滞动态补偿器对系统进行反馈补偿,给出了以严格线性矩阵不等式表示的时滞矩形奇异系统的镇定条件.文献[7]基于二次型矩阵分析方法,给出了奇异凸多面体系统的容许性判据,并将理论结果推广到具有时滞情形的系统.文献[8]设计基于观测器的比例导数控制器,基于线性矩阵不等式提出模糊奇异系统容许的充分条件.
在实际系统的分析与设计时,系统的状态变量通常是从系统外部不易直接测量获得的,所以绝大部分情况下无法利用状态反馈控制器实现系统镇定与性能优化.同时在系统状态变量可以直接获取时,相比于状态反馈,对系统施加输出反馈控制更能够节约系统控制成本并提高系统可靠性.目前在输出反馈控制器设计方法中,主要是用Lyapunov理论研究系统静态输出反馈镇定问题,并结合线性不等式求解,但是其可解条件通常是一个双线性矩阵不等式,这是一个NP-hard问题[9].文献[10]中利用Lyapunov矩阵的合同变换求得系统的静态输出反馈控制器,其局限在于系统矩阵需要满足限定条件.文献[11]利用Finsler引理,通过增广矩阵不等式维数,实现静态输出反馈控制器设计,增加了计算复杂性.文献[12]利用辅助矩阵变量方法解决一类奇异系统输出反馈控制设计问题,但是系统输出矩阵要求是行满秩的.文献[13]通过将原系统输出变量增广为新系统的状态变量方法,可以解除输出反馈增益矩阵与系统输入/输出矩阵耦合关系,虽然避免出现双线性矩阵不等式,但系统维数的扩大增加了计算的复杂度.文献[14]通过输出矩阵满足矩阵等式约束,解决了一类随机奇异系统的有限时间输出反馈控制器设计中.
基于以上分析,本文利用矩阵迹不等式的理论,研究奇异系统容许性与静态输出反馈设计问题.首先给出了奇异系统容许的代数判据.并以此为基础,设计静态输出反馈控制器保证闭环奇异系统容许,并给出相应算法完成控制器参数求解.本文所提方法规避已有结果中利用Lyapunov稳定性理论结合线性矩阵不等式求解静态输出反馈增益矩阵面临双线性矩阵不等式的难点.同时该方法也可以应用到正常系统输出反馈设计中.
符号说明:
考虑如下的一类奇异系统:
(1)
其中:x(t)∈Rn是状态向量;u(t)∈Rm是控制输入;y(t)∈Rp是测量输出向量;E,A,B,C是给定的适维实数常矩阵,且
rank(E)=q≤n
定义1[1]考虑如下奇异系统:
(2)
1) 系统(2)被称为正则,如果存在s∈C,使得
det(sE-A)≠0;
2) 系统(2)被称为无脉冲,如果
deg(det(sE-A))=rank(A);
3) 系统(2)被称为稳定,如果
σ(E,A)⊂C-;
其中:
C-={s|s∈C,Re(s)<0};
σ(E,A)={s|det(sE-A)=0} .
4) 系统(2)被称为容许的,如果系统是正则、无脉冲且稳定.
引理1[1]奇异系统(2)是无脉冲的充要条件为
(3)
引理2[15]设A∈Rn×n,且满足
(4)
则A的特征值全部分布在复平面的左半部分.
引理3 如果矩阵E满足rank(E)=q, 那么存在可逆矩阵U和正交矩阵V满足:
(5)
证明 对矩阵E进行奇异值分解,有
(6)
其中:P,Q∈Rn×n是正交矩阵;E1∈Rq×q满足
且σi>0,i=1,2,…,q是矩阵非零奇异值.然后选取
引理得证.
引理4[16]设A,B∈Rn×n,则有
1) tr{A+B}=tr{A}+tr{B};
2) tr{kA}=ktr{A};
3) tr{AT}=tr{A};
4) tr{AB}=tr{BA};
5) 若矩阵U∈Rn×n为可逆矩阵,那么
tr{U-1AU}=tr{A}.
引理5 如果矩阵V是正交矩阵,那么
tr{[H(VTAU)]2}=tr{[H(A)]2} .
(7)
证明 利用引理4中矩阵迹性质以及V为正交矩阵满足VTV=I,可以推出
tr{H(VTAV)2}=
tr{VTH(A)VVTH(A)V}=
tr{VTH(A)H(A)V}=
tr{H(A)H(A)VTV}=
tr{H(A)H(A)}=tr{[H(A)]2}.
引理得证.
本文的目的是设计静态输出反馈控制器
u(t)=Ky(t)=KCx(t),
(8)
使得相应的闭环系统
(9)
是容许的.
给出自治奇异系统(2)容许的代数判据,在此基础上设计输出反馈控制器实现闭环系统容许.
定理1对于自治奇异系统(2),如果满足
(10)
tr{H(M-1AN-1T)}<
(11)
则自治奇异系统(2)是容许的.
其中:q=rank(E);M=E+A-AVUE;N=E+A-EVUA;T=M+N-2E-EVUAVUE; 矩阵U和V如式(5)中所示.
证明 根据引理3,存在可逆矩阵U和正交矩阵V满足:
(12)
(13)
由引理1,推出系统(2)是无脉冲的.再由式(10),结合式(12),(13),可得
2q+rank(A4)=
n+rank(E)=n+q.
(14)
这时rank(A4)=n-q,且A4∈R(n-q)×(n-q),即A4是可逆矩阵.
进一步,考虑到
det(s0E-A)=
det(U-1)det(V-1) .
(15)
因此系统(2)是正则的.
下面证明系统(2)是稳定的.考虑系统状态变换如下:
然后,对系统(2)进行受限等价变换,可以推出
(16)
考虑到是可逆的,进行整理可得
(17)
接下来,考虑到
UMV=UEV+UAV-UAVUEV=
以及
计算可得
VTM-1AN-1TV=
VTM-1U-1UAVVTN-1U-1UTV=
(UMV)-1UAV(UNV)-1UTV=
则有
tr{H(M-1AN-1T)}=tr{VTH(M-1AN-1T)V}= tr{H(VTM-1AN-1TV)}= tr{H(A*)}+tr{H(In-q)}= tr{H(A*)}+n-q.
(18)
再利用引理5,可得
(19)
利用式(11),可以推出
(20)
det(sE-A)=0与det (sIq-(A1-
综上,系统(2)是容许的.
注1 利用矩阵迹的不等式,定理1给出了奇异系统容许的代数判据.相对于文献[3]中的基于线性不等式的容许性判定定理需要寻求矩阵P满足不等式
(21)
定理1只需要计算系统矩阵的矩阵迹和矩阵秩来完成系统容许性判断,所提方法更加直观简单,易于计算.
由定理1,直接给出闭环系统(9)的输出反馈控制器设计定理.
定理2闭环系统(9)是容许的,如果满足如下条件:
(22)
(23)
注2 利用定理2求解输出反馈控制器时,不需要对系统输出矩阵C进行特殊设定,这就区别于文献[12]只能针对输出矩阵C是行满秩矩阵情况,以及文献[14]要求C和Lyapunov矩阵P满足等式约束
PCT=CTZ.
也可以利用定理2设计状态反馈控制器使得闭环系统容许,得到如下推论.
推论1 在状态反馈控制器u(t)=Kx(t)下,闭环系统是容许的,如果满足
当E=I,奇异系统可以转化为正常线性系统,系统(1)转化为
(24)
又因为线性系统总是正则和无脉冲的,定理2可以转为线性系统输出反馈控制设计定理.
推论2 在输出反馈控制器u(t)=Ky(t)下,正常系统(24)是容许的,如果满足
其中,
利用定理2求解奇异系统输出反馈控制器参数,关键在于求解式(22),式(23).算法的基本思路将定理条件转换为非线性规划求解问题,并通过Matlab编程实现,其中目标函数为最小化矩阵A+BKC的范数,式(22),式(23)为非线性约束条件,通过算法求解K保证闭环系统容许性.即求解最优化问题:
min‖A+BKC‖.
n-q.
例1 含管理在内的石油催化、裂化过程非常复杂.美国Profimatics公司实现了这一过程的建模和控制,其简化模型为
0=A21x1(t)+A22x2(t)+B2u(t)+F2f(t).
(25)
且输出方程为
y(t)=C1x1(t)+C2x2(t) .
(26)
其中:x1(t)为被调节量,如再生温度、滑阀位置、鼓风机能力等;x2(t)是由影响过程、企业效益和反映企业管理政策的一些变量组成的,如压力、油浆回收率、重油回收率等;u(t)为调节量;f(t)为外部干扰.
系统参数为
整理可得如下的奇异系统:
其中:
由于A22=0,系统存在脉冲行为,不满足容许条件.通过定理2,求解静态输出反馈控制器参数如下:
选取初态x1(0)=2,x2(0)=1,该闭环系统的状态响应如图1所示.
进一步计算系统的极点为(-2.254 3,-0.120 0),可见闭环奇异系统是容许的.
例2[17]考虑如图2所示的电路.此电路是由电压源、电容器、电感器和一个电流控制器组成.
图1 闭环奇异系统状态响应(输出反馈)Fig.1 State response of closed-loop singular system (output feedback)
图2 电路结构图Fig.2 The structure of circuit
按照文献[17]建立系统模型:
(26)
选取C=1 F,L=1 H,γ(t)=2,令
可以将式(26)表示成奇异系统为
其中:
计算
因此系统存在脉冲行为,不满足容许条件.通过推论1,求解得出状态反馈增益矩阵:
K=[-0.873 6 -0.109 5 0.437 9] .
进一步计算
可以得出系统是无脉冲的,同时选取初态x1(0)=2,x2(0)=1,该闭环系统的状态响应如图3所示.
图3 闭环奇异系统状态响应(状态反馈)Fig.3 State response of closed-loop singular system (state feedback)
此时系统的极点为(-0.192 5±1.138 6i),可见闭环奇异系统是容许的.
例3 考虑如下的正常线性系统:
(27)
其中:
取初态x1(0)=3,x2(0)=1,x3(0)=-1,该开环系统的状态响应如图4所示.由图4可以看出系统是不稳定的.
图4 开环系统状态响应Fig.4 State response of open-loop system
进一步由推论2可得输出反馈增益矩阵为
K=-3.221 3.
选取相同初态,该闭环线性系统的状态响应如图5所示.
此时闭环线性系统的极点为(-0.259 6±0.270 7i,-0.153 1),所以此系统是稳定的.
图5 闭环线性系统状态响应Fig.5 State response of closed-loop linear system
本文针对奇异系统的容许性和输出反馈设计问题展开研究.利用矩阵迹不等式理论,提出该类系统容许性代数判据和静态输出反馈控制器设计方法.从根本上避免了已有静态输出反馈设计方法中求解双线性不等式条件的难点问题.同时,对于奇异系统的状态反馈设计以及正常系统的输出反馈设计问题,本文所提方法仍然有效.最后,通过仿真例子证明了所提方法的可行性与有效性.