卡尔曼分解、互质分解下讨论最小实现以及零极点相消

1、目录TOC o 1-5 h z HYPERLINK l bookmark2 一.卡尔曼分解下讨论零极点相消与最小实现11.1卡尔曼分解概述11.2卡尔曼分解的原理21.3卡尔曼分解与最小实现以及互质分解4零极点相消与最小实现的关系2.1概述42.2单变量系统的能控性、能观性与传递函数零极点相消之间的关系。.52.3最小实现的判据6 HYPERLINK l bookmark29 三.利用互质分解83.1互质分解与卡尔曼分解83.2matlab上的验证113.3最小实现与互质分解以及卡尔曼分解之间的关系12摘要：本文主要在卡尔曼分解以及互质分解下讨论了最小实现以及零极点相消的问题。讨论了非互质的传

2、递函数会使系统实现时出现不能控或者不能观的部分,从而引出了卡尔曼分解，卡尔曼分解后的能控能观部分的实现为最小实现，系统维数降低，说明出现了零极点相消的情况。系统的维数等于互质分解后传递函数的维数的实现时最小实现，此时的实现也是能控能观的实现。但如果传递函数中消掉的是不稳定的零极点，则不稳定的极点会导致不稳定的状态，出入输出稳定与系统渐进稳定之间是有很大差别的。卡尔曼分解下讨论零极点相消与最小实现1.1卡尔曼分解概述卡尔曼分解，即能控能观性分解，在己知系统状态方程ABCD不能控或者不能观的情况下，对其做矩阵等价变换，使其状态变量划分为能控能观、能控不能观、不能控能观、不能控不能观四个部分。状态方

3、程ABCD能够分解为不能控或者不能观部分说明了两个问题：1.对一个实际系统，并不是所有的状态都能控，也不是说所有的初始状态都能够通过系统输出反映出來，表征输入输出关系的传递函数也仅反映能控能观部分的关系，从而区分输入输出稳定以及系统渐进稳定；2.卡尔曼分解说明了系统实现时的维数是大于最小实现时的维数的（参见第二章最小实现的内容），因此在表示系统传递函数时必然存在零极点相消的现象，或者说零极点相消的现象使得系统实现时存在不能控或者不能观的部分。卡尔曼分解也提醒我们在系统实现时要注意不稳定零极点相消的问题，因为相消的极点如果为不稳定极点，则说明存在不能控或不能观的状态，而且这个状态是处于发散状态的

4、。1.2卡尔曼分解的原理能控性分解跟能观性分解是具有对偶性的，因此这里先讲能控性的分解，要说明能控性分解的原理，这里先证明一个小结论。结论1:对于能控性矩阵C=BAB.如果把B矩阵展开，即B=山U.bp,那么C矩阵可以表示为：C=b.bpAb】.辿.理飞.Ax_1bp那么结论是：如果A*与前面（左边）的向量线性相关，则同样与前面的向量线性相关。要证明这个结论非常简单，因为Abm表示以矩阵A每一列作为一个向量然后进行线性叠加，而A（Abm）,同样是矩阵A每一列的线性叠加，因此如果Abm与一组向量线性相关，则必然虑和入与同样一组向量线性相关。能控性分解：如果能控性矩阵C=BAB.不满足行满秩，设C

5、的秩为斗11，则可以构造等价变换矩阵PTqqJ其中q取自C矩阵的任意i个线性独立列，其余的缶“任意取，只要保证矩阵P非奇异就可以。则通过等价变换x=Px,状态方程可以转化为:11（式子14）利用上面的结论1,可以简单证明上面的式子。因为AP-PA,.入.qj,也就是说瓦第i列是Aq在基q%.%下面的表示乂因为禺八是与q乞线性相关的,因此A前i列为,乂因为Aq钿Ag对于基q.%.qj是线性独立的，因此瓦的斗*n列为空的形式。而B=P_1B,AJ显然,因为B与q礼线性相关，因此E的形式为的形式，至此原式得证。能观性分解以及能控能观分解由于能控性跟能观性的对偶性诫，能观性也可以参照能控的做法进行分解

6、,也可以得到相对应的分解形式：如果对状态变量先进性能控性分解，再进行能观性分解，也就可以得到下面的状态方程：y二Ct-30CqiOji+Dfl（式子1-3）13卡尔曼分解与最小实现以及互质分解再对状态方程进行能控分解时，我们还可以得到另外一个结论：原来状态方程的传递函数与能控性分解后的能控性部分得到的传递函数相等。即：xf=+Bfuy=C凡十Du的传递函数与原状态方程相同。最直接的证明是直接对式子14求传递函数，然后分块计算每一部分的值以及最后得到的传递函数，由于式子1-1为上三角分块矩阵，因此可以利用上三角矩阵的求逆公式來求，在此不再展开。卡尔曼分解后得到能控能观部分的最小实现的结果说明了一

7、个问题，原來的系统如果按照最小实现來构建的话，则系统的维数必然会降低，而系统的维数降低乂说明了原來的传递函数中存在零极点相消的现象。反之亦然，按照第三章传递函数互质与卡尔曼分解的讨论中可以看出，传递函数如果存在公共因子，那么在传递函数实现的时候，必然会存在不能控或者不能观的部分，因此可以按照上面的思路，构建非奇异变换矩阵來实现卡尔曼分解。因此卡尔曼分解、传递函数的互质性以及最小实现之前是互相联系，可以互相推导的。零极点相消与最小实现的关系2.1概述每个线性时不变系统都可以用输入-输出函数：y(s)=G(s)u(s)來描述，且这系统是集中的，用状态方程描述为：x(t)=Ax(t)+Bu(t)(2

8、-1)y(t)=Cx(t)+Du(t)如果状态方程是己知的，那么传递函数阵可以求出：G(s)=C(sl-+Do这计算出來的传递函数阵是唯一的。相反，通过一个给定的传递函数阵求其相对应的状态空间方程的问题，称为实现问题。如果存在有限维状态方程(24)或者说A,B,C,D使的G(s)=C(sI-A)-1B+D则称传递函数阵G(s)是可实现的。并且(A,B,C,D称为G的一个实现。发散的线性时不变系统可以用传递函数阵描述，但是不能用有限维状态方程描述,所以不是所有的G都是可以实现的，如果G是可以实现的，那么它有无线多种实现的方法，不一定要有相同的维数，所以实现问题相当复杂，其中我们称最小维的实现为最

9、小实现。2.2单变量系统的能控性、能观性与传递函数零极点相消之间的关系。(2-1)对应的传递函数为：g(s)=c(sl-A)-Xb=器(2-2)其中fN(s)=c.adj(si-A)b,D(s)=det(sl-A)定理:动态方程2-1能控能观的充分必要条彳雄g(s)无零极点对消，即D(s)和N(s)无非常数的公因子。证明：首先用反证法证明条件的必要tto若有s二So既使N(So)=O,又使D(so)=O:D(s0)=det(si-A)=0,N(s0)=c.adj(si-A)b=0利用恒等式IA)磐豁古如(s)I=(si一A)adj(sl-A)将S二So代入，可得Aadj(s0I-A)=soad

10、j(soI-A)(l)将上式前乘c、后乘b后即有c.Aadj(sl-A)b=soc.adj(s0I一A)b=SoN(so)=0式前乘cA、后乘b,并考虑到的结果后即有cA2adj(soIA)b=soc.Aadj(solA)b=SqN(so)=0以此类推可得N(s)=c.adj(s0IA)b=0c.Aadj(s0IA)b=0c.A2adj(s0IA)b=0c.An1adj(s0I-A)b=0这组式子乂可写成ccAadj(s0I-A)b=03一1因为假设动态方程能观的。上式中前面的能观矩阵是可逆矩阵，故adj(s0I-A)b=0所以我们有adj(s0I-A)b=Sk；oPk(sO)AkbPo(sO

11、)bAbAfPE)=0Pn-l(So).但因pn_i(s)=1故detbAb一=0,这与系统可控的假设相矛盾。矛盾表明N和D无相同因子，即g(s)不会出现零、极点相消的现象。再证充分性：即若N(s)和D(s)无相同的因子，要证明(2-1)是能控能观的。用反证法。设该系统不是既能控乂能观的。不妨设系统是不能控的，这是可按能控性分解，并且可知这时传递函数1c.adi(sl-A)bN(s)g(s)=c(sI3】b=de爲打=胡zTa、-gci.adj(slAJbiNi(s)=C1(Sl-A1)S=det迤I_AQ=丽在上面的式子中,D是n次多项式，而Di(s)是ni次多项式，由于系统不可控,所以珥n

12、，而N(s)和D无相同因子可消茲显然N(s)Ni(s)-F-D(s)Di(s)这和两者应相等矛盾。同样可以证明状态方程也是不能观的。2.3最小实现的判据(A,B,C)为严格真传递函数矩阵G的一个n维实现，则其为最小实现的充分必要条件是(A,B)能控且(A,C)能观测。证明：先证必要性，即已知(A,B,C)为最小实现，欲证(A,B)能控和(A,C)能观测。釆用反证法，反设(A,BQ不能控或不能观测，则可以通过结构分解找到能控且能观测的(AbBhCJ,使Ci(slAJ-iB=G(s),且有diniAidimA(3)表明(A,B,C)不是G(s)的最小实现，从而与己知条件矛盾，故反设不成立，(A,B

13、,C)必为能控且能观测的。必要性得证。再证充分性，即己知(A,BQ能控且能观，欲证(A,B,C)为最小实现。也釆用反证法，反设(A,B,C)能控且能观测，但不是最小实现，这时必存在另一个最小实现(A；B：C)使dinifA0则G(t)=G(t)式中G(t)、G(t)分别为(A,BQ、(A,C)的单位脉冲响应矩阵。对G(t)求各阶导数并利用A和的可交换属性，得到G(t)=CA0B=CeABG(2)(t)=CA20B=G(T(t)=cAn_1eAtB=CeAn-iBG2(n-1)(0=n-lgAt-lp于是，可构造下列L(t)矩阵G(t)G1.G(n7(t)L(t)二G1G2(t).G(n)(t)

14、cBbCA0bcBabCA0AB.CeAn-iBCAeFTB=cCA_CAn-leAtBCAn-leAtAgCAn-1eAtAn1B.CA111.G(n-l)(t)G(n)(t).G2(n7(t)出上田ABAn-1B=VoUtO式中，Vo、Uc分别是(A,B,C)的能观和能控性的判别矩阵。当t=0时，有L(0)=VOUC同理，可导出L(t)=V；eA，tUc,L(O)=V；U；在上式中，V；、U；分别为(A：B：C)的能观测性和能控性的判别矩阵。由G(t)=G(t)又有L(t)=Lr(t),L(0)=L(0),故VoUc=V；U；由己知(A,B,C)能控且能观，则rankV0=n,rankUc

15、=n表示H=VOUC有rankHniinrankVOlrankUc=n又因为VjH=VVoUc从而有Uc=(VjV0)-1VjH故有n=rankUcniinrankVjVOlrankVOlrankH=rankH(8)由于式(7)和式(8)同时成立，所以必有rankH=rankVoUc=n于是，利用式VoUc=V；U；和乘积阵秩的关系式，得到n=rankV0Uc=rankV；U；n,rankU；n这表示diniA1dimAf与反设相矛盾，故反设不成立，即不存在比(A,B,C)维数更小的实现。充分性得证.所以传递函数阵无零极点相消，则实现是完全能控的且完全能观的，而这样的系统被称为是最小实现。2例

16、如,G(s)=話它的分子分母有一个公因子s-1,故其存在零极点相消，它的能控性实现为00-11x=Ax+bu=100 x+0010.0.0ri14.00Ab=10.011001,A2b=A*Ab=0.0.-100.0010010.0.1.-100所以,rank(b,Ab,A2b)=3/满秩。该实现为完全能控的。ro0-171o-3100=01_4Lo10.cA2=cA*A=1iio-?00-101-1044lo10.44-2V3,所以该实现为不能观的，综上所述系统不是最小实现。而对于系统Gi(s)=(s+l)/4(s2+s+1),它的能控性实现为:x=Ax+bu=gJx+Juy=ex=1/41

17、/4Ab七JJI=01/41/4g；=l/21/4rank(b,Ab)=2,iw秩，该实现为完全能控的,rank)=2,满秩，该实现也是完全能观的。故该实现为最小实现。利用互质分解3.1互质分解与卡尔曼分解先用一例子来说明，互质分式传递函数与卡尔曼分解之间的关系设一系统I的传递函数为：G(s)N(s)D(s)/?1S3+02S，+03S+04s4+ais3+ais2+ass+式八亠有Y(s)=G(s)*U(s)其中Y为系统输出信号的拉氏变换，U为系统输入信号的拉氏变换令D(s)V(s)=U(s)(式A2)定义变量v(t)为V(s)的反拉氏变换,有Y(s)=N(s)V(s)式A-3设状态变量为：

18、刍(t)orX(s)=X】(s)s3X?(s)s2X3(s)s_X4(s)1V(s)式a-4将式A-4代入A-2得sXl(s)=-alXl(s)-a2X2(s)-a3X3-a4X4(s)+u(s)对其进行拉式反变换有xr(t)=-a1-a2-a3-a4x(t)+u(t)同理，将式A4代入A3得y(t)=(P4p3卩2|3lx(t)综合上述各式可得5_Z4T10000 x=Ax+bu=0100 x+0u0010_0式A-5y=cx=040302创X现在开始讨论式A-1中N(s)与D(s)是否互质与A-5的能控能观性的关系由式A-5可得该系统的能控矩阵为：1-a2-+laxa2-a3TOC o 1

19、-5 h zC_01-ax-a2001-a】0001Det(C)=l,很明显系统I传递函数的该实现是一定能控的。如果A-1中N(s)与D(s)不互质的话，则必然存在一非零常数r使得N(r)=PlrA3+322+卩3r+P4=0D(r)=rA4+alrA3+a2M2+a3r+a4=0A-6令Z=1r2r1为非零向量，由A-6可知N(r)=cz=O,另外故有-a2-a2-ar3r41000r2r30100rr200101r其中0为能观矩阵，很明显0不满秩，因此当N(s)与D不互质时A-5必然不能观。假设N与D互质时，如果系统不能观则存在A的一特征值r和非零向量z使得A-1I-z=0c即Az=rz和

20、cz=0由此可得N(r)=PlrA3+P22+卩3r+P4=0应此可得r为N(s)=0的一个根，乂因为r为A的特征根即为特征方程D(s)=0的一根，从而N(s)与D(s)有公因子s-r,故而N与D不互质。与假设矛盾。综上所述，A-5当且仅当A-1互质的情况下才能控能观。对A-1取转置有01000G)=G(s)=!Y(sI-因此状态方程10-3_-4y=io也是式A-1的一种实现00此实现必定能观，由对偶律与上述描述可知，该实现当且仅当A-1为互质分式时才能控。推广至一般情况可以得出结论：对于SISO系统而言，如果系统的状态空间表达式能控且能观，则G(s)=b(sI-1c必定互质，即对于卡尔曼分

21、解得出来的能控能观状态空间表达式其传递函数必定互质O3.2matlab上的验证下面通过对A-1与其能控能观设置不同参数用matlab來对上述讨论进行验证。在matlab的command窗口中输入一下语句symss;D=(s+l)*(s+3)*(s+5)*(s+6);N=(s+2)*(s+5)*(s+6);expand(D)expand(N)可得D=sA4+15*sA3+77*sA2+153*s+90N=sA3+13*sA2+52*s+60为一对有公因子的两多项式，利用上述参数设置能控型参数输入A=-15-77-153-90;l000;0100;0010;c=l135260;O=obsv(A,c

22、);rank(O)可知O的秩为2系统不能观，可知传递函数不互质时，系统的能控标准实现不能观。利用上述参数设置能观型参数，输入A=-15100;-77010;-153001;-90000;b=l;13;52;60;C=Ctrb(A?b);rank(C)可知C的秩为2,系统不能控再在maltab的command中输入symss;D=(s+l)*(s+3)*(s+7)*(s+8);N=(s+2)*(s+5)*(s+6);expand(D)expand(N)可得D=sA4+19*sA3+119*sA2+269*s+168N=sA3+13*sA2+52*s+60为一对互质分式利用其系数设置能控标准型，输入A=-19-119-269-168;c=l13