上海交大杜秀华老师《现代控制理论》第四章 线性系统的能控性和能观性6

三、 多输入多输出系统的标准形

旺纳姆（Wonham）标准形和龙伯格（Luenberger）标准形。 1.多输入多输出系统的能控标准形

考虑线性定常系统

Σ:??Ax?Buxy?Cx

x为n维状态向量，u为p输入向量，y为q维输出向量

如果系统能控，则系统的能控性矩阵的秩为n，即Q中有n个线性无关列。

cQc?[b1b2?bp?Ab1Ab2?Abp???An?1b1An?1b2?An?1bp]

对多输入系统，p?1，Qc中有np列，所以，在Qc中可以找出很多种n个线性无关列的情况。这里介绍两种寻找

n个线

性无关列的方法，以构成状态变换阵，将状态空间描述形式变换为旺纳姆能控标准形和龙伯格能控标准形。

定理 [旺纳姆能控标准形]对完全能控的线性定常系统，存在线性非奇异变换

?1x?Px?Qx

使状态空间表达式转化为旺纳姆能控标准形

ΣcW:式中

??Ax?Buxccy?Ccx

?A11?0Ac?Q?1AQ??????0?

A12?A22??A1m???A2m?????0Amm??

?0?????Aii??(?i??i）???0?????1?1

1????1(?1?1)0?1??????????,i?1,2,?,m ??0???01???????11????0??ij1??ij2Aij????(?i??j)????ij?10?0?????,j?i?1,?,m???

?0?0???1(m?1)??1p??0????????0?????1?????1???Bc?QB?????

0???????????0??????1???n(m?1)np??Cc?CQ（无特殊形式）

证明：见书 例 求如下系统的旺纳姆能控标准形

?12?1??10????A???010? B??01?

???103???00??计算系统的能控性矩阵

Qc???BAB?101204??010101?A2B????? ??001042??rankQc?3，系统完全能控。按{ 1b, 1 A,2b1 } Ab和

{ b2, Ab2, A2b2}搜索，可知b1,Ab1,b2是线性独立的。得

?1?2,?2?1。

A2b1???11Ab1??12b1?1??1??0??0??1??1??, Ab??0?, A2b??0?, ?0?????0????0? b1??01111??12?????????????????0???1???4???4???1???0??4，?12解得 ?11＝－e12?b1?4

e11??11b1?Ab1得

?1??1???3??1??0????0???0? e?b??0? e?b??1?e11??4 ?0121212?????????? ????0????1????1???0???0??由Ab2?-?21b1??121e11??122e12，即

?2??0???3??1??1?????1????0????0?21??121????122??

?????0???0???1???0??解得 ?21由此可得

??1, ?121?0, ?122?2

?A11Ac?QAQ???0?1?010??00?A12??? B?Q?1B??10? ?-442?c???A22?????001???01??

??310??001?Q?式中，变换矩阵为 ?? ??100??系统的旺纳姆能控标准形为

?011??00?????442?x??10?ux???? ???001???01??龙伯格能控标准形相对于旺纳姆能控标准形在系统极点配置中应用更广，这是由于在构成状态变换阵时，相对均衡地考虑了各个输入的作用。考虑线性定常系统

Σ:??Ax?Buxy?Cx

x为n维状态向量，u为p输入向量，y为q维输出向量。

如果系统完全能控，设矩阵B的列线性无关，这在实际上通常总能成立，或可进行处理获得，bi为B的第i列，那么，系统的能控性矩阵Qc可表示为

Qc?[b1b2?bp?Ab1Ab2?Abp???An?1b1An?1b2?An?1bp]

对上式从左到右依次找出

n个线性无关列，即若某个列不

能表示成其左边各线性独立列的线性组合就为线性无关列，否则为线性相关列。考虑到Qc的组成特点，如果

Aibm,m?1,?2,p和其左边线性独立列线性相关，则,Ai?1bm也和其左边线性独立列线性相关，又由于假设B为满秩，则可将此搜索结果重新排列，并组成非奇异矩阵P：

??1??1??1P?1??bAb?Ab?bAb?Ab???bAb?Abp?11222pp?1?12p

式中，?1??2????p?n， {?1,?2,??p}为系统能控性指数集。

?1P对定义矩阵求逆，并表示成分块矩阵形式：

T?e11???????e1T?????1?1?P?(P)????eT? ?p1?????eT??p??1p式中，块的行数为?i,i?1,2,?,p。再在矩阵

12pP中，取出各

TTTe,e,?,ep?，个块矩阵的末行，即为1?2?并按如下方式构造变

换矩阵P：

cLPcLT?e1????TeA1???????T??1??e1?A???????T?ep????T?ep?A????????1T?ep?A???1111pppp

定理[龙伯格能控标准形]对完全能控的线性定常系统，存在线性非奇异变换

??PcLx?QcL?1x x使状态空间表达式转化为龙伯格能控标准形

ΣcL:??x?u??A??Bxcc ?x?y?Cc??A11??A??1??QAQ?21AccLcL????A??p1

?A12?A??A22p2??A1p???A2p????????App??

?010?????1????????Aii????(?i??i）???0???0???**???

0???????,i?1,2,?,p0? 1??*???0?????Aij?0(?i??j)??*0?0?????,i?j,i,j?1,2,?,p0?0?

?*?*??0???????0???1**????Q?1B????BccL??0? ??????0????1????CQCccL（无特殊形式）

*表示可能的零元或非零元。

证明：从略。

例 求如下系统的龙伯格能控标准形

?0?3A????1??10??0?10?31?? B???014?1???0?10??0010?0??1??0?

计算系统的能控性矩阵Qc

Qc???b1?0?1???0??0b2?Ab1Ab2?A2b1A2b2?A3b1A3b2??13?00?3?3?10?10?30?? 1144131351??00?1?1?3?3?9?0011442b,b,Ab,Ab2线性无关，式中，12所以，系统完全能控，且?1?1，1?2?3。由此构成矩阵

P?1，并求其逆阵P。

4?0?3?10??1413?

?0?1?3?01?0?12?P?1??bbAbAb2?1?12???0??0?1??1P????3??110?2?013??00?4?

?001?将

P分为两个块阵，第一个块的行数为1，第二个块的行数

T11为3。取出P中的第一行e和第四行e，组成变换矩阵PcL，并求其逆阵QcL。

T23T?e11??1?T??1ePcL??T23????e23A??1?T2???e23A??010?2?001??000?

?010??0?1?1QcL?PcL???0??00?2?30??001?

?1?10?01由此可求出系统的龙伯格能控标准形为

?0?0??Q?1AQ??AccLcL?0??10?010??001??1?34?10?1?0??Q?1B=?BccL?0??00?0??0? ?1?2.多输入多输出系统的能观标准形

利用能观性和能控性的对偶关系，可以方便地得出完全能观系统的旺纳姆能观标准形和龙伯格能观标准形。 定理 [旺纳姆能观标准形]对完全能观的线性定常系统，存在线性非奇异变换

??Po?1x?Qox x使状态空间表达式转化为旺纳姆能控标准形

ΣoW:??x?u??A??Bxoo ?x?y?Co式中

??A11????P?1AP??A21Aooo??????Al1?0??1?0?Aii?????????00?A22??Al2???? ????All????1?i??0??1???,i?1,2,?,l 0????????01??i1????ij1?ij2??ij?i??????0??0? ,i? j?1,?,l,j?1,2,?,lAij????????0??0??

??P?1B（无特殊形式） Boo?0??1?0??1???????Co?CPo??????????(l?1)1???????????p1??????0??1?

??(l?1)n?????pn??定理 [龙伯格能观标准形]对完全能观的线性定常系统，存在线性非奇异变换

?x?Po?1x?Qox

使状态空间表达式转化为龙伯格能控标准形

ΣoL:??????Ax?Buxoo??

y?Cox式中

??A11???A21?1?Ao?PoLAPoL???????Aq1?0?1???0Aii????????0?0?A12??A22???Aq2?A1q???A1q????

???Aqq???0*?0???1????,i?1,2,?,q0??????0???01*??

?0?0*??0?0*??? ,i?jAij?? ???????0?0*????1Bo?PoLB（无特殊形式）

?0??1???Co?CPoL??*???

?*0??1???能观标准形在系统观测器设计中有重要应用。 应该指出，单变量系统的传递函数与能控标准形和能观标准形的参数有明显的关系。已知传递函数可直接写出系统的能控标准形和能观标准形，不必求出变换矩阵，反之，由标准形也可以直接写出传递函数。但在多变量系统中，不存在这样的对应关系，必须求出变换矩阵，才能求出标准形。

4.7 线性系统的结构分解

介绍如何通过非奇异变换，将系统的状态空间描述按能控性和能观性进行结构分解。

通过状态空间描述的结构分解可以建立状态空间描述和系统传递函数之间的关系，为最小实现问题的提出提供理论依据。状态空间描述的结构分解在系统综合状态反馈和系统镇定问题中也有重要作用。 一、 按能控性的系统结构分解

按能控性分解是通过线性非奇异变换，将系统的完全能控状态和不能控状态，以明显的方式区分开来表示。

考虑线性定常系统

Σ:??Ax?Buxy?Cx

式中，x为输出向量。

n维状态向量，u为p输入向量，y为q维

定理 [系统按能控性结构分解] 对不完全能控的系统，其能控性判别矩阵为Qc的秩为n1

rankQc?rank??BAB?An?1B???n1?n

构造非奇异变换矩阵

P?1???q1?qn11P

qn1?1?qn??

式中，q1,?,qn为U中任意选取的n1个线性独立的列，

qn1?1,?,qn为任意的，只要使

P为非奇异的n?n个线性独立

1的列，那么线性非奇异变换现按能控性的结构分解

???A?xcc??????xc??0y???Ccx?Px或x?P?1x可使系统实

A12??xc??Bc???????uAc??xc??0?

?xc?Cc???x??c?xc为n1维系统能控状态，xc为n?n1维不能控状态。

证明：若记

T?p1???P?Q?1????T??pn??，由于PP?1?I及

P?1???q1?qn1qn1?1?qn??可以得出

piTqj?0,?i?j

由于Aqi(i?n1?1,?,n)为q1,?,qn的线性组合，由此，并利

1Tp用式iqj?0,?i?j可得

piTAqj?0,i?n1?1,?,n,j?1,?,n1

于是可得

A?PAP?1T?p1Aq1????pTAqn1??T1?pn?1Aq1?1???pTAq?n1????Tp1Aqn1Tp1Aqn1?1???TpnAqn11TpnAqn11?1?TpnAqn1?11TpnAqn1?1?1?1?TpnAqn1?TpnAqn1?1?Tp1Aqn????Tpn1Aqn??Tpn1?1Aqn???? TpnAqn???Ac????0A12??Ac??1n1同样，

B的所有列是q,?,qTi的线性组合，则有

pB?0,i?n1?1,?,n

即有

T?p1B??????T?pnB??B?cB?PB??T1?????pn1?1B???0?? ??????pTB??n?对C，无特殊形式，可以表示为

C?CP?1???Cq1?Cqn1Cqn1?1?Cqn??Cc???Cc??

由此得出系统按能控性分解的结构式。

进而证明x为能控状态。由于

crankU?rankU?rank??B?Bc???0AB?An?1B??

n?1AcBc?AcBc??0?0?n?1AcBc?AcBc???n1?rank??Bc由于Ac为n1?n1维矩阵，利用凯勒－哈密尔顿定理可

n1?1B,AB,?,A知，ABc(i?n1,?,n?1)可表示为ccccBc的线性组

ic合。基于此，可得

rank??Bcn1?1AcBc?AcBc???n1

即(Ac,Bc)完全能控，也即xc为n1维系统能控状态。同理可证其余为不能控状态。 证毕

从上述定理可以看出，不完全能控的系统通过线性非奇异变换后，系统的状态空间被显式地分成能控部分和不能控部分，其中，能控部分为n1维子系统

??Ax?Ax?Buxccc12ccy1?Ccxc

不能控部分为n?n维子系统

1??Axxcccy2?Ccxc

两者之间只存在由不能控部分到能控部分的耦合作用。

系统按能控性结构分解后的模拟结构图如图所示。

能控部分Ac?c????1x??CcuA12y1??yy2???2x??Cc不能控部分Ac

系统按能控性分解的模拟结构图

在MATLAB中函数ctrbf能对系统按能控性进行结构分解，但由于它的变换矩阵P的选取和本定理中的排列次序不

?1同，它的变换结果的形式为

?Ac??A210??0?Cc?,??,??Ac??Bc?Cc??

由于变换矩阵选取的相对任意性，因此，系统按能控性分解的结果随变换阵的不同而不同。 例 给定线性定常系统

?110??01????010?x??10?ux???????011???01??,y??111?x

由于rankB?2，只需判断

rank?B?0111???2?3AB??rank?1010??

??0111??可知系统状态不完全能控，且按能控性结构分解后能控状态

xc为二维。

按定理选取非奇异变换矩阵P

P?1?011????100??

??010??P?1的第一、第二列为能控性矩阵的第一、二列，它们线

性独立，P?1的第三列，也是最后一列，为任意选取的使P?1?1rank(P)?3?n，非奇异。 非奇异的列向量。显然上述

对系统进行线性非奇异变换x?Px，通过计算，得到

?010??111??011??100??????????1A?PAP??001??010??100???110?????10－1????111????010???001??010??01??1?????B?PB??001??10???0???10－1????01???0?011???1C?CP?1??111??100?????010??0??1?0??21?

系统按能控性分解后的状态空间描述为

?100??10???xc????xc????110?01??????x???u?xc??001??c??00?????

?x?y??121??c??xc?系统的低维能控子系统为

1010?????xc??xc??uy??12?xc??

?11??01?二、 按能观性的系统结构分解

利用对偶关系，可以得出系统按能观性结构分解的定理。

定理 [系统按能观性结构分解] 对不完全能观的系统，其能观性判别矩阵

R的秩

?C??CA???n?nrankR?rank? ???1?n?1??CA?构造非奇异变换矩阵

P

?p1???????pn?P???p?n?1? ???????pn??11式中，p1,?,pn为

1R中任意选取的n1个线性独立的行，

pn?1,?,pn为任意的，只要使

1P为非奇异的n?n1个线性独

?1立的行，那么线性非奇异变换x?Px或x?Px可使系统

实现按能观性的结构分解

???A?xoo??????xo??A21y???Co0??xo??Bo???????uAo??xo??Bo?

?xo?0???x??o?式中，xo为n1维系统能观状态，xo为n?n1维不能观状态。

从上述定理可以看出，不完全能观的系统通过线性非奇异变换后，系统的状态空间被显式地分成能观部分和不能观部分，其中，能观部分为n1维子系统

??Ax?Buxooooy1?Coxo

不能控部分为n?n1维子系统

??Ax?Ax?Buxooo21ooy2?0

两者之间只存在由能观部分到不能观部分的耦合作用。

系统按能观性结构分解后的模拟结构图如图所示。

Ao???能观部分?ou?xoyCoA21??????xo?不能观部分Ao 系统按能观性分解的模拟结构图

在MATLAB中函数obsvf能对系统按能观性进行结构分解，它是函数ctrbf的对偶形式。

由于变换矩阵选取的相对任意性，因此，系统按能控性分解的结果随变换阵的不同而不同。

两种结构分解的结果在形式上是互为对偶的。系统能观性结构分解也可用对偶方法，运算步骤类似于化标准形的对偶方法，只是将其中的化标准形运算改为能观性结构分解的运算。

例 给定线性定常系统

?00?1??1????10?3?x??1?ux???????01?3???0??,y??01?2?x

计算系统的能观性判别矩阵的秩

?c??01?2???rank?1?23??2?3rank?cA???? 2???cA????23?4??可知系统状态不完全能观，且按能观性结构分解后能观状态

xo为二维。

按定理选取非奇异变换矩阵P

?01?2??P??1?23??

??001??P的第一、第二行为能观性矩阵的第一、二行，它们线

性独立，P的第三行，也是最后一行，为任意选取的使P非奇异的行向量。显然rankP?3?n，P非奇异。

对系统进行线性非奇异变换x?Px，通过计算，得到

0??01?2??111??211??01??010??102????1?20?A?PAP?1??1?23????????????001????111????001???10?11??01?2??1??1???1????1?B?PB??1?23??????????001????0???0??211???100C?CP?1??01?2??102??????001??

系统按能观性分解后的状态空间描述为

?010??1?????x??xo???o??????1?20??????1?u?xo??10?1??xo??0?????

?x?y??100??o??xo?系统的低维能观子系统为

?01??1??xo??xo???uy??10?xo ???1?2???1?三、 按能控性和能观性的系统结构分解

如果线性定常系统不完全能控也不完全能观，若对该系统同时按能控性和能观性进行结构分解，则可以把系统结构显式地分解成能控且能观、能控不能观、不能控但能观、不能控不能观这四个部分。当然，系统是否能分解成这四个部分主要取决于系统本身的特性，并非所有系统都有这四部分组成。

系统按能控性和能观性结构分解的卡尔曼分解定理。 定理[系统按能控性和能观性结构分解] 对不完全能控和不完全能观的系统，通过引入特定的非奇异变换，可使系统实现标准形式，即有

???A 0 A 0??x??B??xcoco13co????????1??xco??A21 Aco A23 A24??xco??B2???u?x???0 0 A 0??xco?? 0?co?co???????????x??xco?? 0?????co??0 0 A43 Aco?? ?xco???xco??y???C1 0 C3 0???x?co???xco???其中，xco为能控且能观状态，xco为能控不能观状态，xco为不能控而能观状态，xco为不能控不能观状态。并且

(Aco,B1,C1)是系统能控且能观子系统。

定理给出的标准形式，在形式上是惟一的，而在具体结果上并不惟一，其结果可由选取的变换阵计算得出。

?1ΣcoC1u?A21A13?y??2ΣcoA23ΣcoC3A24A43Σco

系统卡尔曼分解的模拟结构图

反映系统输入输出特性的传递函数G(s)只能反映系统中能控且能观子系统的动力学行为，即有

G(s)?C(sI?A)B?C1(sI?Aco)B1?Gco(s)

传递函数只是对系统的一种不完全描述，如果系统中添加或去掉不能控或（和）不能观的子系统，不影响系统的传递函数。系统可由传递函数完全表征的充要条件是系统为完全能控且能观的，即系统是不可简约的。

在MATLAB中函数minreal能求出系统的完全能控且能观子系统。

例 给定线性定常系统

?1?1?1???510000??x1??1?x?x?2??0?50000??x2??5????????x?3??00?3100??x3??4???????????x4??000?300??x4??0?x?5??0000?11??x5??1????????6???00000?1????0?x???x6????x1??x??2??y1??320500??x3??y???140200??x???4??2???x5?????x6??3?7??3??u1???0???u2?6??0??

试将该系统按能控性和能观性进行结构分解。

由于系统具有约当标准形形式，利用能控性约当标准形判据和能观性约当标准形判据可知系统不完全能控、不完全能观。容易判定：系统能控且能观状态变量为x,x，系统能控

12但不能观状态变量为x,x，系统不能控但能观状态变量为x，

354系统不能控不能观状态变量为x。

6重新排列状态变量的顺序，即

?x3??x1?xco???;xco???;xco??x4?;xco??x6?

?x2??x5?则有

?x1??1????xco??x2??0?x??x??0x??co???3????xco??x5??0

??????xco??x4??0?x6????0?00000??x1??x?10000???2?01000??x3?????Px 00010??x4?00100??x5????00001??x6????由此可得系统按能控性和能观性结构分解的标准形为

??510000??1???0?50000??xco??5?xco???????00?3010????4xco??xco??????????000?101??xco??1?xco????????0000?300xx?co??????co?????00000?1???0?xco????y1??310050??xco??y???140020??x??co?2?????xco?3?7??3??u1???6???u2?0??0? ?显然，如果系统为约当标准形时，可以很方便地求得系统按能控性和能观性结构分解的标准形。其它情况是不能直接看出和分解的。

4.8 离散时间系统的能控性和能观性

一、 离散系统的能控性 1. 能控性定义

考虑离散时间线性定常系统的状态方程

x(k?1)?Ax(k)?Bu(k)

式中，x为

n维状态向量；u为p维输入向量。

定义[离散系统能控性定义] 对系统，如果对任意初始状态x(0)?x0和任意终端状态xf，存在一个有限长度输入序列

u，使系统状态由x转移到x0f，则称此系统或

(A,B)对是状态完全能控的，或简称此系统或(A,B)对

是能控的。否则，则称此系统或(A,B)对是状态不完全能控的，简称系统或(A,B)对是不能控的。

从上面的能控性定义可以看出，和连续线性系统不同，对连续线性系统而言，如果系统是能控的，则存在输入可以在很短的非零时间区间里使系统实现从任意状态到另一个任意状态的转移；而对离散时间系统来说，如果系统是能控的，则输入序列必须具有一定的长度，在下面的能控性判据证明中，我们可知输入序列长度应不小于能控性指数1. 能控性判据

定理4-18 [能控性格拉姆矩阵判据]离散时间线性定常系统，系统完全能控的充分必要条件为使如下定义的格拉姆矩阵

n?1?。

Wdc(n?1)??AkBBT(AT)k

k?0为非奇异。

定理[能控性秩判据] 离散时间线性定常系统，系统完全能 控的充分必要条件为

rank[BAB?AB]?n

n?1式中，n为矩阵

A的维数。

Qd?[BAB?An?1B]

称为系统的能控性判别矩阵。

证明：利用系统状态运动关系式，可知k?n时，系统的状态

x(n)为

x(n)?Ax(0)??An?1?iBu(i)

ni?0n?1可写为

?u(n?1)??u(n?2)???AB?An?1B????? ???u(0)?x(n)?Anx(0)???B这表明，对任意的x(0)和x(n)，当且仅当Qd满秩成立时，输入序列存在。 证毕

从上述定理的证明过程，结合能控性指数的定义可知，实现状态的任意控制，输入序列的长度为至少为?。上述定理同时给出了实现状态转移的输入序列的求取方法。 例4 离散时间线性定常系统状态方程为

?100??0??x(k)??0?u(k)x(k?1)??02?2????

????110???1??计算系统的能控性矩阵

Qd???BAB?000???A2B????0?2?4???10?2??

容易判定，rankQd?2?n。根据能控性秩判据可知，该系统不完全能控。 二、 离散系统的能观性 1. 能观性定义

考虑离散时间线性定常系统零输入时的状态空间描述

x(k?1)?Ax(k)y(k)?Cx(k)

式中，

x为n维状态向量；y为q维输出向量。

定义4-9 [能观性定义] 对系统，若能根据在有限个采样时间点上输出y(k)的观测值，唯一地确定系统的任意非零初始状态x0，则称此系统或(A,C)对是状态完全能观的，或简称此系统或(A,C)对是能观的。否则，则称此系统或(A,C)对是状态不完全能观的，或简称此系统或(A,C)对是不能观的。 1. 能观性判据

定理 [能观性格拉姆矩阵判据]离散时间线性定常系统，系统完全能观的充分必要条件为使如下定义的格拉姆矩阵

n?1Wdo(n?1)??(AT)kCTCAk

k?0为非奇异。

定理 [能观性秩判据] 离散时间线性定常系统，系统完全 能观的充分必要条件为

?C??CA???nrank???? ?n?1??CA?式中，n为矩阵A的维数。

?C??CA??Rd?????

?n?1??CA?称为系统的能观性判别矩阵。

能控性指数和能观性指数对于离散系统而言还具有特别的意义，即能控性指数决定了离散系统实现任意状态转移的最短输入序列个数，能观性指数决定了输出序列的最少需观测的长度，以便由此确定系统的任意初始状态。

i.

对原点的能控性和能达性

有关能控性的定义，在不同的文献中有所不同。其区别主要是对系统输入实现状态转移的状态有不同的限定。主要有以下三种：

1) 实现由任意状态到任意状态的转移的定义方法，如这里的能控性定义就采用了这种方法。

2) 实现由任意状态到零状态的转移的定义方法，称为对原点的能控性。

3) 实现由零状态到任意状态的转移的定义方法，称为从原点出发的能控性，这种情况常称为能达性。

对于连续时间系统，状态转移矩阵总是非奇异的，这三种定义是一致的。所以一般不再讨论能达性。对于离散时间系统，如果系统矩阵A是非奇异的，那么，这三种定义对系

统来说也是一致的，但当系统矩阵

A奇异时，定义

1)和定

义3)还是等价的，而对定义2)和定义3)不再等价。下面用例子说明它们的不同。考虑离散时间线性定常系统如下

?010??0??x(k)??0?u(k)x(k?1)??001????

???000???0??由于系统的能控性矩阵秩为零，因此，按能控性定义1)系统是不能控的，按定义3)系统也是不能达的。但由于

3Ak?0,k?3,所以，有x(3)?Ax(0)?0，对任意初始状态x(0)都成立。因此，按能控性定义2)系统是能控的。所以，对离散系统当系统矩阵A奇异时，定义2)和定义3)不再等价，分别用2)定义能控性和3)定义能达性。采用1)的定义方法，可使能控性问题的讨论比较简捷。

ii.

离散化系统的能控性和能观性

对连续时间线性系统离散化为离散时间线性系统进行分析和设计是实现对系统进行计算机控制的一种常规处理方法。本小节不加证明地给出有关连续时间线性系统在时间离散化后仍能保持系统的能控性和能观性的条件。

考虑连续时间线性定常系统，其状态空间描述为

??Ax?Bu,Σ:xy?Cxt?0

取采样周期为T，并采用零阶保持器，则对应的时间离散化

系统为

ΣT:x(k?1)?Ax(k)?Bu(k),y(k)?Cx(k)k?0,1,2?

ATAt式中A?e,B?(?0edt)B。

T定理4-22 [离散化系统保持能控能观性条件] 对连续时间线性定常系统

Σ经离散化后得到系统ΣT，若系统Σ是不能

控的（不能观的），则其离散化系统ΣT也是不能控的（不能观的）。若系统

Σ是能控的（能观的），则其离散化系统ΣT为能控的（能观的）的充分条件是系统采样周期T的值应满足，当有

Re(?i??j)?0

时，成立

T?2?lIm(?i??j),l?1,2,?

式中，?i为系统矩阵A的特征值。

上述定理在系统是单输入时，此条件也是必要的。 例 线性定常系统的状态空间描述为

??3???1x???0y??01?7?5??1??0?u00?x?????10???0??

2?x可判定上述系统完全能控、能观，系统的特征值为

?1,?1?j2，三个特征值具有相同的实部，按照定理，离散

化系统完全能控、能观的条件是当且仅当T?T?T?2?l??l，且22?l?0.5?l，l?1,2,?。即系统完全能控、能观的条件为42?l?0.5?l，l4?1,2,?。

0.?。 用MATLAB检验如下，若取l?1，即T?5在MATLAB

中输入

a=[-3 -7 -5;1 0 0;0 1 0]；b=[1;0;0]；cd=[0 1 2]; [ad,bd]=c2d(a,b,pi/2)； 可得离散化系统状态方程为

? -0.1039 0.2079 0.5197??-0.1039??x(k)?? 0.1039?u(k)x(k?1)?? -0.1039 -0.4158 -0.5197????

??? 0.1039 0.2079 0.3118??? 0.1376??y(k)??012?x(k)

离散化系统的能控性矩阵和能观性矩阵可以通过函数ctrb和obsv产生。在MATLAB中输入

ud=ctrb(ad,bd);vd=(ad,cd); 可得

? -0.1039 0.1039 -0.0045?? 0.0000 1.0000 2.0000??,V?? 0.1039 0.0000 0.1039?Ud?? 0.1039 -0.1039 0.0045??d???? 0.0000 0.0432 0.0864?? 0.1376 0.0537 0.0059????

用rank(ud)和rank(vd)可算出秩都为2，因此，在取采样周期为T?0.5?时，原来能控能观的系统离散化后，既不能控，也不能观。由此，验证了定理。

本文来源：https://www.bwwdw.com/article/ywp5.html