《卡尔曼分解、互质分解下讨论最小实现以及零极点相消》由会员分享，可在线阅读，更多相关《卡尔曼分解、互质分解下讨论最小实现以及零极点相消（14页珍藏版）》请在金锄头文库上搜索。
1、目录一. 卡尔曼分解下讨论零极点相消与最小实现11.1卡尔曼分解概述11.2卡尔曼分解的原理21.3卡尔曼分解与最小实现以及互质分解4二. 零极点相消与最小实现的关系2.1概述42.2单变量系统的能控性、能观性与传递函数零极点相消之间的关系。.52.3最小实现的判据6三. 利用互质分解83.1互质分解与卡尔曼分解83.2 matlab上的验证113.3最小实现与互质分解以及卡尔曼分解之间的关系12摘要：本文主要在卡尔曼分解以及互质分解下讨论了最小实现以及零极点相消的 问题。讨论了非互质的传递函数会使系统实现时出现不能控或者不能观的部分, 从而引出了卡尔曼分解，卡尔曼分解后的能控能观部分的实现为
2、最小实现，系统 维数降低，说明出现了零极点相消的情况。系统的维数等于互质分解后传递函数 的维数的实现时最小实现，此时的实现也是能控能观的实现。但如果传递函数中 消掉的是不稳定的零极点，则不稳定的极点会导致不稳定的状态，出入输出稳定 与系统渐进稳定之间是有很大差别的。卡尔曼分解下讨论零极点相消与最小实现1.1卡尔曼分解概述卡尔曼分解，即能控能观性分解，在已知系统状态方程A B C D不 能控或者不能观的情况下，对其做矩阵等价变换，使其状态变量划分为能控能观、 能控不能观、不能控能观、不能控不能观四个部分。状态方程A B C D能够分解为不能控或者不能观部分说明了两个问 题：1.对一个实际系统，并
3、不是所有的状态都能控，也不是说所有的初始状态 都能够通过系统输出反映出来，表征输入输出关系的传递函数也仅反映能控能观 部分的关系，从而区分输入输出稳定以及系统渐进稳定；2.卡尔曼分解说明了 系统实现时的维数是大于最小实现时的维数的（参见第二章最小实现的内容）， 因此在表示系统传递函数时必然存在零极点相消的现象，或者说零极点相消的现 象使得系统实现时存在不能控或者不能观的部分。卡尔曼分解也提醒我们在系统实现时要注意不稳定零极点相消的问题，因为 相消的极点如果为不稳定极点，则说明存在不能控或不能观的状态，而且这个状 态是处于发散状态的。1.2卡尔曼分解的原理能控性分解跟能观性分解是具有对偶性的，因
4、此这里先讲能控性的分解，要 说明能控性分解的原理，这里先证明一个小结论。结论1:对于能控性矩阵C = B AB . lB如果把B矩阵展开，即8 =由 b, . bp,那么C矩阵可以表示为：C = b】. bp 叫.Abp . Ax-1bp 那么结论是：如果Abm与前面（左边）的向量线性相关，则A+tm同样与前 面的向量线性相关。要证明这个结论非常简单，因为Aibm表示以矩阵A每一列作为一个向量然后进行线性叠加，而同样是矩阵A每一列的线性叠加，因此如果Abm与一组向量线性相关，则必然A+tm与同样一组向量线性相关。能控性分解：如果能控性矩阵C = B AB .寂-电不满足行满秩，设C的秩为叫vn
5、,则可以构造等价变换矩阵pT=q qJ其中q取自C矩阵的任意叫个线性独立列，其余的任意取，只要 保证矩阵P非奇异就可以。则通过等价变换x=Px,状态方程可以转化为:11Ac Anirxcl0 Jtiejy=Cf G】*Du（式子1-1）利用上面的结论L可以简单证明上面的式子。因为APT = P/ , pT = q . % . qj,也就是说3第i列是加在基q . * . qj下面的表示，乂因为是与q线性相关的,，乂因为Aq”A%对于基q . %是线性独立的，因此M的七n列为 1 1的形式。而B = P-1B,显然，因为B与qq、线性相关，因此B的形式为| |的形式，至此原式得证。能观性分解以及
6、能控能观分解由于能控性跟能观性的对偶性福，能观性也可以参照能控的做法进行分解, 也可以得到相对应的分解形式：;: =京幻 + ?:=|巳明仕+成（式子1-2）如果对状态变量先进性能控性分解，再进行能观性分解，也就可以得到下面 的状态方程：y = CL-S 0 Cqi Oji+Dfl（式子1-3）1.3卡尔曼分解与最小实现以及互质分解再对状态方程进行能控分解时，我们还可以得到另外一个结论：原来状态方 程的传递函数与能控性分解后的能控性部分得到的传递函数相等。即：xf =+ BfUy = Cfxr 十 Du的传递函数与原状态方程相同。最直接的证明是直接对式子l-i求传递函数，然 后分块计算每一部分
7、的值以及最后得到的传递函数，由于式子l-i为上三角分块 矩阵，因此可以利用上三角矩阵的求逆公式来求，在此不再展开。卡尔曼分解后得到能控能观部分的最小实现的结果说明了一个问题，原来的 系统如果按照最小实现来构建的话，则系统的维数必然会降低，而系统的维数降 低乂说明了原来的传递函数中存在零极点相消的现象。反之亦然，按照第三章传 递函数互质与卡尔曼分解的讨论中可以看出，传递函数如果存在公共因子，那么 在传递函数实现的时候，必然会存在不能控或者不能观的部分，因此可以按照上 面的思路，构建非奇异变换矩阵来实现卡尔曼分解。因此卡尔曼分解、传递函数 的互质性以及最小实现之前是互相联系，可以互相推导的。零极点
8、相消与最小实现的关系2.1概述每个线性时不变系统都可以用输入-输出函数：y(s) = G(s)u(s)来描述，且这 系统是集中的，用状态方程描述为：x(t) = Ax(t) + Bu(t) (2-1)y(t) = Cx(t) + Du(t)如果状态方程是己知的，那么传递函数阵可以求出：G(s) = C(sl - A)】B + Do 这计算出来的传递函数阵是唯一的。相反，通过一个给定的传递函数阵求其相对 应的状态空间方程的问题，称为实现问题。如果存在有限维状态方程(2-1)或者说A, B, C, D使的G(s) = C(sl - A广】B +D,则称传递函数阵G(s)是可实现的。并且A, B,
9、C, D)称为G(s)的一个实现。 发散的线性时不变系统可以用传递函数阵描述，但是不能用有限维状态方程描述, 所以不是所有的G(s)都是可以实现的，如果G(s)是可以实现的，那么它有无线多 种实现的方法，不一定要有相同的维数，所以实现问题相当复杂，其中我们称最 小维的实现为最小实现。2.2单变量系统的能控性、能观性与传递函数零极点相消之间的关系。(2-1 )对应的传递函数为：g(s) = c(sI-A)-1b =籍告夸=黑 (2-2) QC V 买SI /I JU J其中，N(s) = c. adj (si - A)b,D(s) = det(sl - A)定理:动态方程2-1能控能观的充分必要
10、条彳牒g(s)无零极点对消，即D(s)和 N(s)无非常数的公因子。证明：首先用反证法证明条件的必要忙若有s=s既使N(So)=O,又使D(so)=O:D(s0) = det (si 一 A) = 0,N (s0) = c. adj (si - A)b = 0利用恒等式(sI-A)(sI-A)- = (sI-A)*=IB(D(s)I = (si 一 A)adj(sl - A)将S=So代入，可得Aadj(s0I - A) = soadj(soI - A)(l)将上式前乘c、后乘b后即有c. Aadj(sl 一 A)b = soc. adj(s0I - A)b=SoN(so) = 0(2)式前乘
11、cA、后乘b,并考虑到的结果后即有cA2adj(soI A)b = sc. Aadj(sI A)b = s/N(s) = 0 ，以此类推可得N(s)=c.adj(s0I A)b = 0c. A adj(s0I A)b = 0c. A2 adj(s0I A)b = 0 .这组式子乂可写成adj(s0I-A)b = 0c.AnTadj(s()I A)b = 0 c cA . . cA】T.因为假设动态方程能观的。上式中前面的能观矩阵是可逆矩阵，故 adj(s0I-A)b = 0所以我们有adj(s0I - A)b = Ek；oPk(sO)AkbPo(sO)b Ab A% Pi)=0 .Pn-l(S
12、o).但因Pn-i(s)三1故detb AbAnTb=0,这与系统可控的假设相矛盾。 矛盾表明N(s)和D(s)无相同因子，即g(s)不会出现零、极点相消的现象。 再证充分性：即若N(s)和D(s)无相同的因子，要证明(2-1)是能控能观的。用反证法。设该系统不是既能控乂能观的。不妨设系统是不能控的，这是可按能控性分解，并且可知这时传递函数g(s) = c(sl A广=c.adj(sI-A)b_ N(s) detfsI-A) =D(s)C. adj (si A1)b1det整I - Ai)Ni(s)D)在上面的式子中,D(s)是n次多项式，而Di(s)是ni次多项式，由于系统不可控, 所以i】
13、i 儿而N(s)和D(s)无相同因子可消去，显然N(s) Ni(s)-F- D(s) Di(s)这和两者应相等矛盾。同样可以证明状态方程也是不能观的。2.3最小实现的判据(A, B, C)为严格真传递函数矩阵G(s)的一个n维实现，则其为最小实现的 充分必要条件是(A,B)能控且(A,C)能观测。证明：先证必要性，即己知(A,B,C)为最小实现，欲证(A,B)能控和(A,C)能观测。 采用反证法，反设(A,B,C)不能控或不能观测，则可以通过结构分解找到能控且能 观测的(Ai,Bt,Ci),使Csl AQTBi = G(s),且有diniAi dimA(3)表明(A,B,C)不是G(s)的最小
14、实现，从而与己知条件矛盾，故反设不成立，(A,B,C) 必为能控且能观测的。必要性得证。再证充分性，即己知(A,B,C)能控且能观，欲证(A,B,C)为最小实现。也采用反证 法，反设(A, B, C)能控且能观测，但不是最小实现，这时必存在另一个最小 实现(A,B,C)使dinifA 0G(t)=G (t)式中G(t)、G(t)分别为(A,B,C)、(A, B, C)的单位脉冲响应矩阵。对G(t)求各阶导数并利用A和e役的可交换属性，得到G(i)(t) = CAB = CeABG(2)(t) = CA2 = CeB G(nT)(t) = cAn-1eAtB = CeAfiBGe a. G(nT)(t)-L(t)=CeBG1 G2(t)G(I)(t) G(n)(t)CeAB.G(n)(t). . G2(nT)(t)_. CeAtAnTBcCAeBCAeAB. CAe%nTBCA_CAn
