传递函数矩阵的状态空间最小实现

1、传递函数矩阵最小实现方法降阶法人们在设计复杂系统时，总是希望在构造系统之前用模拟计算机或数字计算机对所设计的系统进行仿真，以检查系统性能是否达到指标要求。给定严格真传递函数矩阵G(s),为寻找一个维数最小的(A,B,C),使C(sI-A),B=G(s),则称该(A,B,C)是G(s)的最小实现，也称为不可约实现。最小实现是系统实现的一种非常重要的实现方式，关于最小实现的特性，有下列几个重要结论：(1)(A,B,C)为严格真传递函数矩阵G(s)的最小实现的充要条件是(A,B)能控且(A,C)能观测。(2)严格真传递函数矩阵G(s)的任意两个最小实现(A,B,C)与(AB,C)之问必代数等价，即两

2、个最小实现之间由非奇异线性变换阵T使得式子A=t,at,b=t,b,C=ct成立。(3)传递函数矩阵G(s)的最小实现的维数为G(s)的次数%,或G(s)的极点多项式的最高次数。为了寻求传递函数矩阵的最小实现，就意味着要把系统中不能控和不能观测的状态变量消去而不至于影响系统的传递函数。求最小实现的方法有三种：1、降阶法。根据给定的传递函数矩阵G(s),第一步先写出满足G(s)的能控型实现，第二步从中找出能观测子系统；或者第一步先写出满足G(s)的能观测型实现，第二步从中找出能控子系统，均可求得最小实现。2、直接求取约当型最小实现的方法。若G(s)诸元容易分解为部分分式形式，运用直接求取约当型最

3、小实现的方法是较为方便的。3用汉克尔矩阵法求取最小实现的方法。下面主要研究降阶法(先求能控型再求能观测子系统的方法)并举例说明。先求能控型再求能观测子系统的方法设(pxq)传递函数矩阵G(s),且p<q时，优先采用本法。取出G(s)的第j歹I，记为Gj(s),是Uj至y(s)的传递函数矩阵，有Gj(s)=gij(s).gqj(s)T=Pij(S)qij(s)Pqj(s)/qqj(s)记dj(s)为qij(s),%(s)的最小公倍式,则1TGj(s)=nij(s)nqj(s)dj(s)nini1设dj(s)=saj,njs-aj,isaj,。则nj(s)=%,njsnj'+Bkn-

4、s'/+十'j,is+%,o,i=1,.q在此dj(s)是q个子系统传递函数的公共部分，由单输入-多输出系统的实现可知，能用能控规范I型的Aj、bj实现dj(s),由n(s)的诸系数确定Cj,这时Gj(s)的实现为一0I”Aj=可。一为,。一可，,或研bj0I：Ir0IJ1I'M1j,01j,1Cj"-ppILqj，0qj,1p11J,nj4PJqj,nj1qnj令j=1,，P，便可得Gj(s)的实现为Ann二A2bnxp二Cqn=C1C2AP1Cp1bPP当p<q时，显见A、B、C的维数均较小，且有£nj=n。上述实现一定能j1控，但不一点能

5、观测，需要找出能观测部分，为此需要判别（A,C）的能观测行。若（A,C）能观测，则（A,B,C）为最小实现;若rankQo=rankCCA=n0<n_CAn则从Qo中选出no个线性无关行，记为S;在附加（n-）个任意行（通常为单位矩阵In的任意行），记为§,即MlSS1TVn构造nn的非奇异变换阵T,T=1SI5一引入变换乂=丁乂，由能观测性的结构分解可知一氏Bo-AOoAAo其中能观测子系统（Ao,Bo,Co）即为所求的最小实现。（A0,B0,C0）有如下简化求法：记丁为T1-iSnonSiUnnoU1(n_no)n=SUu/二SUIISi_SiUSU1S3nooIno有SU

6、=L>o由ct'=c!Uu/LICucu=coo,有Co=CUAUU1SAUSAU"_Ao飞iAUSiAuJ灯o,有Ao=SAUAo由TBS15一BSB1百,1sB一.Boj有Bo=SB于是由能控型化为能控能观测型的简化步骤可归结为:i.构造S阵（从Qo中选出n。个线性无关行）；2 .由SU=1%,求出U阵;3 .计算最小实现。A=SAU,B0=SB,C0=CU。由于S选择的任意性及求解U的任意性，最小实现不唯一，但最小实现的维数是唯一的，且系统都是能控能观测的。下面举例说明该法。例1、已知传递函数矩阵G(s),求最小实现。上G(s)=jS+1s+3s+1s2解：化G(

7、s)为严格真传递函数矩阵(?(s)G(s)=_1_s+1-11s3-1-C?(s)Ds+1求6(s)的最小实现gi(s)=_1_s+11l_s+1g2(s)=s+31：.s+2(s2)(s3)s2Ls-3.令d1(s)=s+1,d2(s)=s2+5s+6,其能控规范I型实现为C-?lC1-!-1jb20=，1C22小1-3-1-G?(s)的能控型实现为A-IL0-10.0011010-01_jC二匕1小；-1C=rank11cA一一121-1-11-1-6-3J62=rank=3=n&(s)的最小实现。G(s)的最小实(A,C)的能观测性判别：由于rankC=2=mCrankQ。=ra

8、nk*A1即(A,C)能观测。(A,B,C)能控且能观测，即为现为(A,B,C,D)。例2、求下列G(s)的最小实现维数及最小实现G(s)-(s1)(s2)-24s62s3(s+1)(s+3)-1_(s+1)(s+2)(s+1)(s+2).解(1)确定最小实现维数n:所有G(s)的一阶子式的最小公分母为(s+1)(s+2)；二阶子式只有一个0,其分母为任意常数。故所有子式的最小公分母仍为(s+1)(s+2),有n5=2o(2)g(s)=产；6(s+1)(s+2)1-21g2(s)=2s3(s1)(s2)1!-1令di(s)=dz(s)=(s+1)(s+2),其能控规范I型实现为I36G=1|L24103C2=2<1A/,0C-IC1C21(A,C)的能观测性判别：由于rankC=2=m一64321'C-20-10一c|CA-8-6-4-3rankQ0=rank=rankCA=rank=2<4CA10-20-1CAn口1i1210654623_(A,C)不完全可观测。从Q0中选出二行构成S阵，S=，64321-20-10由SU=12求U阵:U