课件-现代控制理论-刘豹第三版-第三章
课件-现代控制理论-刘豹第三版-0绪论

1.1 控制理论的发展历程
经典控制理论 现代控制理论 新发展——大系统理论 智能控制
发表了《Cybernetics》- -控制学科诞生
经典控制理论(1935-1950)
存在的问题
1、简单对象 单输入单输出、线性、时不变系统
2、缺乏系统化方法 图形化方法,依赖于设计人员的经验
3、达到的性能要求较低,不能处理多目标性能。
面临的挑战
对象日益复杂化、控制性能要求不断提高。
3、解决的方法2---现代控制理论
1.1.3 现代控制理论的发展
70年代中期,自动化的应用开始面向大规模、复杂的系统, 如大型电力系统、交通运输系统、国民经济系统等,运用现 代控制理论方法已不能取得应有的成效,于是出现了:
——大系统理论
是指规模庞大、结构复杂、变量众多的信息与控制系统,如 生产过程、交通运输、生物工程、社会பைடு நூலகம்济和空间技术等复 杂系统。 研究对象:着重解决生物系统、社会系统这样一些众多变量 的大系统的综合自动化问题。 研究方法:时域法为主 重点:大系统多级递阶控制 核心装置:网络化的电子计算机
引言
面对未知及不断变化的世界,人 类发明了无数理论和工具,控制 论就是其中之一。
控制论是一种思想、一种方法、 一种工具、一门学科。
人类在20世纪所取得的巨大技术 成就,控制科学与技术的作用非 常显著。
引言
钱学森曾经从生产力,特别是技术革 命的进程分析了控制论的产生和发展。
现代控制理论3

3
x3
b31
b32
b3 p
u3
xn
n xn bn1 bn2 bnp u p
➢其中的对角部分应用对角形可控判据 ,即要求输入矩阵B中 不出现全零行,则系统对角部分的状态可控。
➢约当部分,展开后可得
x&1 1x1 x2 b11u1 b12u2 L b1pup
x&2 1x2 b21u1 b22u2 L b2 pup
要求约当块最后一行对应的输入矩阵B中的行不出现全零行,则 系统约当部分的状态可控。
(b)
x&1 1 1
x&2
1
x1 b11 b12 L
x2
b21
b22
L
x&3
x&4
1 1
x3 x4
b31 b41
限 时 间 间 隔 0≤t≤nT 内 , 针 对 任 意 初 态 x(0) 和 任 意 终 态
x(n),当k=n时,u(0),u(1),…,u(n-1)一定存在。
因此,k=n时,u(0),u(1),…,u(n-1)的解存在的条
件,即为系统可控时应满足的条件
令k=n ,
n1
x(n) Gn x(0) Gdetb
det Sc 0
Ab b1
b2
1b1 2b2
2b1b2 1b1b2
1 2
b1 0, b2 0
输入阵中无全零行
A
1
0
1
1
b
b1 b2
det Sc det b
Ab b1 b2
1b1 b2 1b2
1b1b2 (1b1 b2 )b2 b22
(2) 线性定常连续系统可控性判据
现代控制理论第三章PPT

( A
c1
,bc1 ) 的能控性,其中
1 0 0 0 A c1 0 0 2 5
解:
0 0 1 0 0 1 1 10
0 0 b c1 0 1
0 1 0 0 0 0 1 10 A3 c1b c1 0 1 10 101 1 10 101 1025
若取
u( t ) B( t )T ΦT ( t0 ,t )Wc1( t0 ,t f )x( t0 )
tf t0
x( t f ) Φ( t f ,t0 )[ x( t0 )
Φ( t0 ,t )B( t )B( t )T ΦT ( t0 ,t )Wc1( t0 ,t f )x( t0 )dt ]
( k 1,2, , n 1 )
假设 F( t ) Φ( t0 ,t )B( t ) 对上式关于时间t求一阶、二阶、直至n-1阶导数 ,可得
(t ) Φ (t , t )B(t ) Φ(t , t )B (t ) F 0 0
(t ) Φ(t0 , t )A(t )B(t ) Φ(t0 , t )B
实现最优控制和最优估值及其它系统综合
与校正的必要条件。
4.1 系统的能控性
[定义]设系统的状态方程为
(t ) A(t )x(t ) B(t )u(t ) x
对于任意非零初始状态 x(t0 ) ,如果存在容许控制u(t ) ,在有限时区
t [t0 , t f ] 将其转移到状态空间原点,即 x(t f ) 0 ,则称系统在
(t )] Φ(t0 , t )[A(t )B(t ) B
Φ(t0 , t )B1 (t )
现代控制理论(刘豹版)第三章部分习题及解答

例3-3 有系统如下,判断其是否能控。
u x a a a X ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=100100010321解:由A 阵为友阵,设其特征根分别为321,,λλλ,有 (1)若321,,λλλ互异,其变换阵为一个范德蒙德矩阵:⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=232221321111λλλλλλT 由3210λλλ==⇒=T ,与题设矛盾,故0≠T 从而得:⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=-1231231******λλλλλλT ,所以⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=-1231231λλλλλλb T 各行均不为0因此系统是能控的。
(2)若3121,λλλλ≠=,则有:⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=232213121101λλλλλT 由0)(213≠-=λλT ,得⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=-1****1**1311λλT T , 故⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=-11311λλb T 各行均不为0,因此系统是能控的。
(3)若321λλλ==,则有:⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=22010012211λλλT ,从而得⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=-120100122111λλλT 故⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=-1001b T 最后一行元素不为0,因此系统是可控的。
3-6 已知系统的微分方程试写出其对偶系统的状态空间表达式,及其传递函数。
u y y yy 66116=+++解:系统的传递函数为状态空间表达式为∑1其对偶系统的状态空间表达式为∑2由系统∑2中,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+++=--**6******61161)(2312S s s A sI故其传递函数为=-=-21222)()(b A sI c s W 6116623+++S s s3-7 【习题3-7 】已知能控系统的状态方程A ,b 为试将该状态空间表达式化为能控标准型 解:,所以系统是能控的。
特征多项式: 即:系统能控Ⅰ型:从而得 []x y u x x 0061006116100010=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---= 61166)()()(23+++==s s s s U s Y s W []xy ux x 1000066101101600=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---= ⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡-=113421b A []n rankM Ab b M==⎥⎦⎤⎢⎣⎡-==2711110543212+-=---=-λλλλλA I []⎥⎦⎤⎢⎣⎡-=⎥⎦⎤⎢⎣⎡-⎥⎦⎤⎢⎣⎡-=⎥⎦⎤⎢⎣⎡=12161501171110111a b Ab T c ⎥⎦⎤⎢⎣⎡----=⎥⎦⎤⎢⎣⎡-=--6211811216111c T 10,501=-=a a系统能控Ⅱ型:3-8 已知能观系统的A ,b ,c 为试将该状态空间表达式化为能观标准型 解:,所以系统是能观的。
现代控制理论--刘豹优秀PPT

1.3.1 从系统框图出发建立状态空间表达式 该法是首先将系统的各个环节,变换成相应的模拟结构图,并把每个积
1.4.2 传递函数中有零点时的实现
22
此时,系统的微分方程为: 相应地,系统传递函数为:
设待实现的系统传递函数为:
因为
上式可变换为
(26)
23
令 则 对上式求拉氏反变换,可得: 每个积分器的输出为一个状态变量,可得系统的状态空问表达式:
24
或表示为:
推广到 阶系统,式(26)的实现可以为:
25
(4)
6
在经典控制理论中,用指定某个输出量的高阶微分方程来描述系统的 动态过程。如上图一所示的系统,在以 作输出时,从式(1)消去中间变量 i,得到二阶微分方程为:
(5)
其相应的传递函数为:
(6)
回到式(5)或式(6)的二阶系统,若改选 和 作为两个状态变量,
即令
则得一阶微分方程组为:
7
设单输入一单输出定常系统,其状态变量为 状态方程的一般形式为:
2
1.1.3 状态方程 以状态变量
为坐标轴所构成的 维空间,称为
状态空间。
1.1.4 状态方程
由系统的状态变量构成的一阶微分方程组称为系统的状态方程。 用图下所示的 网络,说明如何用状态变量描述这一系统。
图一
3
根据电学原理,容易写出两个含有状态变量的一阶微分方程组:
亦即
(1)
式(1)就是图1.1系统的状态方程,式中若将状态变量用一般符号 ,
现代控制理论3章

3.2.2 状态能观性的判据 已知系统的动态方程
x’(t ) A(t ) x(t )
y C (t ) x(t )
x(0) x0
1、线性定常连续系统能观测性的格拉姆矩阵判据
线性定常连续系统{A,C} 状态完全能观测的充要条件是:存在 时刻t1>t0,使如下定义的能观性格拉姆(Gram)矩阵为非奇异。
16
例:给定线性定常系统的状态方程,判断能控性。
1 x1 1 x1 0 x1 x 2.5 1.5 x 1 u y 1 0 x 2 2 2
解:
1 1 1 0 AB 1 1 2.5 1.5 1 1 rank[ B AB] rank =1 2 系统不能控 1 1 s 1.5 1 1 1 0 s 1 Cadj ( sI A) B s 2.5 1 2.5 g ( s) s 1 sI A ( s 2.5)( s 1) ( s 1) 2.5 s 1.5
0 t1 0
t1
T x0 Wc (t0 , t1 ) x0 [ BT T (t0 , ) x0 ]T BT T (t0 , ) x0 d 0
T x0 Wc
(t0 , t1 ) x0 B (t0 , ) x0 d (t0 , ) x0 0
At k 0 n 1
y (t ) k (t )CAk x(0)
k 0
n 1
y (t ) 0 (t )Cx(0) 1 (t )CAx(0) n 1 (t )CAn 1 x(0) C CA x(0) y (t ) 0 (t ) 1 (t ) ... n 1 (t ) n 1 CA C CA =n rank n 1 CA
《现代控制理论》第三版(刘豹

bpppp K x x x x J xx K xJx K K K x x u K K K xJ J K x x p x K xK xxx Ax+Buy=CxA=BCM Mf f M MM y f K c c B y yM y f K c B y K c c B y yffx Ax+Buy=CxA=By y y y u usG ss s sa a ab b bAB n n n b nba a ab a b aC bba b a a b a aabb a b a T TBCa ab a aa aa aaby y y y u u us sG ss s sb b b bCBTW s C SI A B d SI ASSS S S S SSSSSAA P P P A P P P T P P TPPP P T ATJ T ATT BCTW s W s I W s W s2 4.At e10140 A5 15!4!Atett 11A P P112T jP222A P P212T jP1122Tj j412cos 2sin 214sin 22cos 22At t tt e Te T t t 112111()444S S SI A S SS22t t 22))jt jt e21141A223323232233112!3!13715,2!66281344,1562!At e I AtA t A t t t t t t t t t t t t t 6e At e 3031()31()4t tt t t e e t e e 0Ate I 1T2.010001(1,0)xx u y x 1(0),()1x u t 22.1...2!A I Att ()00201()122111222tAt A t t t td t td t x x = =01te e Bu d)(10t t 110110110t27.231.00()()(0)(),()At At t t x x x t e x e t u K BK Ku ()00()t At A t e e B d x x K ()00()t At A t e e d x B ()00|At A t e e0At At e e B x 2t t u K 0(0)x x ()01()t d ()00t At A t e e d B x K 1()00|At A t t e A e B x K 10()At At e A e I B x K 1x x K Kx K B K ()1()A e e B t 0At t3()1()t t u k 0(0)x x ()00(t)1()tAt A t e eB t d x x K 00tAt A At e e d e B x K ()'A A e e A e A 0t 000|t t A tA A e e d A e A 0|t110|A A t Ate A e 11()A At A te A e I 210()()At At At At t e A e I A tee B x x 210()At At e A e I A t B x KK d 00t 1t A A e d Ae e d2-9.2-21s,0.1T s 11u 2u 111212x x Ku x x u 12u u y210120;1.11110010I A P P 1101T P 22200011I A P P 112P 0111T 1T 11101100()0AT T T G T eT T e e 000()0101001(1)T TAt t T T K H T e dtBdt e e K e K T e T T 111Te t e 11K T e TT 0(1)011221212(1)()()()()(1)()()()0()11()(1)0()(1)T TT T12x k x k G T H T u k x k x k x k ex k eu k K e u k K T e T12()()21()u k x k y k x k 0.1Ts0.110.1210.12(1)0(1)()()x k ex k e K e u k K e12()()21()x k y k x k 1Ts10.1)(1x k ()x k ()u k 12109).9)0.9)0.1(1111221112(1)()0(1)()11()(1)0()1x k x k e x k x k eu k K e u k Ke12()()21()x k y k x k 112222(1)()()()(1)()()()()100()()10x k x k u k TA I TBx k x k u k x k u k T KT x k u k TT20.1T ss12112212(1)()0.90(1)()0.11()0.10()0.1x k x k x k x k u k K u k 1T 1122(1)()()000(1)()()111x k x k u k K x k x k u k 1111At e t2AB BA3-11x u2x4 ya, b, c, d1122334412340001000110000100010x x a x x b u x x c x x dx x yx x 02M b AbA bA b 22100001001a a a a c a ac c a c d,a, b, c, d4rankM 323223222230011100C CA c NCA a c b c cCA a ac c b bc c c 4rankN 00,a, b, c, d 222121110x x a b u x x c dx yx 11b a M b Abc d11b a Ma b c c dda b c d1110C NCA a bb3abcdabcd331121[]1M b Ab 210M211211C NCA 21N212121234A 11b10c 123411M bAbrank 2M 3412det 0M12341TC N C CACA2det 0N20432233231341322232831441333232223232332232323det 1410142510486624M M 012123b Ab A b1220010141413C NCA CA det 0N 32226825()14343s s s W s s s s s I 0134A 01b52Tc1II01TA2b 01T 1d 311200011012MbA bdet 0M 3rankM 3425d c Ab010001100cR det 0c3R 1cR 0011000101014ˆ1202c cA R AR 1ˆ100Tc b R b cc cR 3()W s 231()1s s W s s s2111231001001001001sss12101,,10010000r r rrrI A I II200100100r crrI0010000000000011010000000000000000002r00000000,001001r cr rB I 012cC 00:20000000000c cc cC NC A C A 100011001000001000()36Rank N c 1010100100000110011o R 1000110011000000001000010100000010000o R 00000100110010000010001000oR 100000010121212121ˆˆˆˆ0000000000000000ˆ00110000100011001000o oA AR AR AAA 1ˆˆ1000ˆ00o BBB R B10000112ˆ0ˆˆoCCR C C 100100000001211111111ˆˆˆ00110ˆˆ100,00,00001100ˆ010m mA B C W AA B B C C W s 23()11001001001001ss W s s s s s s120012000110,,011000312()01200000000000000100000001000000010000100mm m om m mmmmIA I I I I 000,0100B 00001000000001o m m mC I :1000101122000000010000010000100100100100001001o o o ooMB A B A B 36rankM 1cR 010010100000100100001000cR 0010000100001000100010000000001000110011000100000cR 10110ˆc o cAR A R 000010100000010100ˆ00000c oBR B 00000000000000000001100100ˆ001000o cCC R 11111111ˆˆˆ010ˆˆ000,100100ˆ001m mA B C W AA B B C C 1n 1001,00M b Ab A b1nA 233-33234-11()TQ P x x 121322233233P P P P P P P P P x 110112131P =0220317.750i221101003()Q x 2TQ x x 131P11233160()Q x ()117.75ix P 11434-22112211221221A a a a a a a1211221211221221a aa a a aQ IPTA12221222122221221001 a a P P P P a a11112112()1a P a P211211221112122()a a P a P a P121221222()1a P a PP=11P1221222P PAP PA QP Pa a P P a a1121111211121112112r P t A A222r P t A A11222111122r a a a a P t A A r t AA1122r t A a a 11221221A a a a a P11P 112212244r r P P P t A At A A22111212212()()04r a a a a t A A 112212210A a a a a 11P 0A0r t A 1122a a 112212210a a a a 1122a a 4-3121()v x22x x 0ex ()0v x22121()266v x x x x x 2336PP 2i0iiiP ()0v x x()v x TA P PAQQI302x xP T111P120P22212()v xx x11x 22x5-21120010,,01011010100990Mb Ab A b 3rankM)100k 2010,,01011010100990M b Ab A b 3rankM2I321111A1112112212100100010010110100111019901001011111 cIT A b Ab b10cIA20120100011111A bKb T0.13221()(11)(11)f *32()122440f 40131K314 1.20.c KKT 15-512Mb Ab A b 24001015det 0M 5-8(1);1100200c A B20121010det detdet10c A B c A B Ec A B;RankM 13kMc A B 05-10(1)cNcA(2)II:01y cTx x(3)12,gGg :12221()det()1g f I A Gc g g g (4):*2()()(2)322f r r rr(*222113ˆ;2002rrx u y rr2ˆˆˆ()0113ˆˆ()002Ax bu G y y r;x u y yr(1)(2)22101- =10000g g A Gc g g 2121=detf I A Gcg g g :*2()()(2)322f r r r r232r Gr12g 11)(det[(-)]。
《现代控制理论》刘豹著(第3版)课后习题答案(最完整版)

第一章习题答案1-1 试求图1-27系统的模拟结构图,并建立其状态空间表达式。
11K s K K p +sK s K p 1+s J 11sK n 22s J K b -++-+-)(s θ)(s U 图1-27系统方块结构图解:系统的模拟结构图如下:)(s U )(s θ---+++图1-30双输入--双输出系统模拟结构图1K pK K 1pK K 1+++pK n K ⎰⎰⎰11J ⎰2J K b ⎰⎰-1x 2x 3x 4x 5x 6x系统的状态方程如下:uKK x KK x KK x X K x K x x x x J Kx J x J K x J Kx x J K x x x ppppn pb 1611166131534615141313322211+--=+-==++--===∙∙∙∙∙∙令y s =)(θ,则1x y =所以,系统的状态空间表达式及输出方程表达式为[]⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡+⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡-----=⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡∙∙∙∙∙∙65432116543211111111265432100000100000000000000010010000000000010x x x x x x y uK K x x x xx x K K K K K K J K J J K J KJ K x x x x x xp p pp n pb1-2有电路如图1-28所示。
以电压)(t u 为输入量,求以电感中的电流和电容上的电压作为状态变量的状态方程,和以电阻2R 上的电压作为输出量的输出方程。
R1L1R2L2CU---------Uc---------i1i2图1-28 电路图解:由图,令32211,,x u x i x i c ===,输出量22x R y =有电路原理可知:∙∙∙+==+=++3213222231111x C x x x x R x L ux x L x R 既得22213322222131111111111x R y x Cx Cx x L x L R x uL x L x L R x =+-=+-=+--=∙∙∙写成矢量矩阵形式为:[]⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡----=⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡32121321222111321000010111010x x x R y u L x x x CCL L R L L R x x x 。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
3.2 线性定常系统的能控性判别
线性定常系统能控性判别准则有两种形式,一种是先将系统进行状态变
换,把状态方程化为约旦标准型 ,再根据 阵,确定系统的能控性; 另一种方法是直接根据状态方程的 A 阵和 B 阵,确定其能控性。 3.2.1 具有约旦标准型系统的能控性判别 1.单输入系统 具有约旦标准型系统矩阵的单输入系统,状态方程为: (1) 或 (2) 式中
3.1 3.2 3.3 3.4 3.5 3.6 3.7 3.8 3.9 3.10
能控性的定义 线性定常系统的能控性判别 线性连续定常系统的能观性 离散时间系统的能控性与能观性 时变系统的能控性与能观性 能控性与能观性的对偶关系 状态空间表达式的能控标准型与能观标准型 线性系统的结构分解 传递函数阵的实现问题 传递函数中零极点对消与状态能控性和能观 性之间的关系
(12)
3.7.2 单输出系统的能观标准型 与变换为能控标准型的条件相似,只有当系统是状态完全能观时,即
有:
系统的状态空间表达式才可能导出能观标准型。
状态空间表达式的能观标准型也有两种形式,能观标准 型和能观标 准 型,它们分别与能控标准 1.能观标准 型 若线性定常系统: (13) 是能观的,则存在非奇异变换: (14) 型和能控标准 型相对偶。
5)如果 6)如果
是能控状态,则 和
也是能控状态, 是任意非零实数。 也必定是能控状态。
是能控状态,则
7)由线性代数关于线性空间的定义可知,系统中所有的能控状态构成
状态空间中的一个子空间。此子空间称为系统的能控子空间, 记为 。
2.线性连续时变系统的能控性判别 时变系统的状态方程如下: (1)
系统在
对多输入系统,其状态方程为: 式中,B 为 阶矩阵; 为 r 维列矢量。
(15)
其能控的充分必要条件是矩阵: 的秩为 。
3.3 线性连续定常系统的能观性
3.3.1能观性定义 能观性所表示的是输出 出方程出发,即 (1) 如果对任意给定的输入 期间的输出 ,在有限观测时间 ,使得根据 ,则称状态 能唯一地确定系统在初始时刻的状态 反映状态矢量 的能力,与控制
3.4 离散时间系统的能控性与能观性
3.4.1 能控性矩阵 M 离散时间系统的状态方程如下: (1) 当系统为单输入系统时,式中 列矢量;G为系统矩阵 3.4.2 能观性矩阵N 离散时间系统的能观性,是从下述两个方程出发的。 (2) 式中, 为 维列矢量;C 为 输出矩阵,其余同式(6)。 ; 为标量控制作用.控制阵 为状态矢量 。 为 维
能观性矩阵。仿连续时间系统,记为N。即
(5)
3.5 时变系统的能控性与能观性
3.5.1 能控性判别 1.有关线性时变系统能控性的几点说明 1)定义中的允许控制 绝对平方可积的,即 ,在数学上要求其元在 区间是
这个限制条件是为了保证系统状态方程的解存在且唯一。 2)定义中的 ,是系统在允许控制作用下,由初始状态 转移到 目标状态(原点)的时刻。 3)根据能控性定义,可以导出能控状态和控制作用之问的关系式。 4)非奇异变换不改变系统的能控性。
3.6.1 线性系统的对偶关系
有两个系统,一个系统 为:
另一个系统
:为:
若满足下述条件,则称
与
是互为对偶的。
式中,
为 维状态矢量; 为 各为 与
各为r与m维控制矢量; 系统矩阵; 各为, 维输出矩阵。
各为 与 ,
与 维输出矢量; 维控制矩阵;
3.6.2 对偶原理 系统 则 的能控性等价于 者说,若 全能控的)。 和 的能观性, 是状态完全能控的(完全能观的),则 是互为对偶的两个系统, 的能观性等价于 的能控性。或 是状态完全能观的(完
根据3.3节中能观性定义,如果知道有限采样周期内的输出 一地确定任意初始状态矢量 导能观性条件。从式(1),有:
,就能唯
,则系统是完全能观的,现根据此定义推
(3) 若系统能观,那么在知道 出 , 时,应能确定 ,现从式(7)可得:
写成矩阵形式:
(4)
有唯一解的充要条件是其系数矩阵的秩等于 。这个系数矩阵称为
判别法之间的关系
众所周知,一个矩阵:
式中,
为列矢量,当且仅当由 为非奇异时,
构成的格拉姆矩阵 列矢量是线性无关的。现在
因此,有 奇异等价。
这个矩阵的列矢量线性无关与
非
3.6 能控性与能观性的对偶关系
能控性与能观性有其内在关系,这种关系是由卡尔曼提出的对偶原理确
定的,利用对偶关系可以把对系统能控性分析转化为对其对偶系统能观性的 分析。从而也沟通了最优控制问题和最优估计问题之间的关系。
3.6.3 时变系统的对偶原理
时变系统的对偶关系和定常系统稍有不同,且其对偶原理的证明也复 杂得多。 对偶原理是现代控制理论中一个十分重要的概念,利用对偶原理可以
把系统能控性分析方面所得到的结论用于其对偶系统,从而很容易地得到
其对偶系统能观性方面的结论。
3.7 状态空间表达式的能控标准型与能观标准型
2)变换矩阵R确定之后.只需经讨一次变换便可对系统同时按能控性
和能观性进行结构分解.但是R阵的构造需要涉及较多的线性空间概念。 3)结构分解的另一种方法:先把待分解的系统化约旦标准型,然后按 能空判别法则和能管判别个状态变量的能控型和能观性,最后按能控能观、 能控不能观、不能控能观、不能控不能观四种类型分类排列,即可组成 相应的子系统。
的各项系数。
2.能控标准
型 (6)
若线性定常单输入系统:
是能控的,则存在线性非奇异变换: (7) 相应的状态空间表达式(6)转换成: (8) 其中
(9)
(10) 并称形如式(8)的状态空间表达式为能控标准 式(9)中的 型。 是系统特征多项式:
(11)
的各项系数,亦即系统的不变量。 式(11)中的是 相乘的结果,即:
型,∑的状态空间
式中,
为系统
为系统
的能控标准II型对应的系数阵;
的能控标准I型对应的系数阵; 的能控标准 型对应的系数阵。
为系统 的对偶系统
2.能观标准 型 若线性定常单输出系统:
(20) 是能观的,则存在非奇异变换
(21)
使其状态空问表达式(20)变换为:
其中
(22)
(23)
(24) 称形如式(22)的状态空间表达式为能观标准 型。
这时式(2)用房承租形式表示,可有:
(3)
(4)
从而可得结构图如图所示。将式(3)带入输出方程式(4),得:
(2)A 为约旦标准型矩阵 以三阶为例:
这时,状态方程的解为:
从而
(5)
由式(5)可知,当且仅当输出.矩阵C中第一列元素不全为零时,y( t ) 中总包含着系统的全部自由分量而为完全能观。 2.直接从A、C阵判断系统的能观性 约旦标准型 系统具有串联型 的结构,如图所示:
测的。
5)如果 和 都是不能观的,则 也是不能观的。
6)根据前面分析可以看出,系统的不能观测状态构成状态空间的一个子
空间,称为不能观子空间,记为
。只有当系统的不能观子空问
。在状
态空间中是零空间,则该系统才是完全能观的。
2.线性连续时变系统能观性判别 时变系统 (4) 在 上状态完全能观测的充分必要条件是格拉姆矩阵 (5) 为非奇异的。 3.5.3 连续时变系统可控性和可观性判别法则和连续定常系统的
是能控的,而
维子系统:
是不能控的。对于这种状态结构的分解情况如图所示,因为 对 作用, 仅作无控的自由运动。显然,若不考虑 可得到一个低维的能控系统。
不起
维子系统,便
至于非奇异变换阵: (7) 其中 个列矢量可以按如下方法构成,前 个列矢量 能控性矩阵M中的 确保 个线性无关的列,另外的 个列 为非奇异的条件下,完全是任意的。 是 在
3.7.1 单输入系统的能控标准艰 对于一般的 维定常系统:
如果系统是状态完全能控的,即满足: 1.能控标准 型
若线性定常单输入系统: (1) 是能控的,则存在线性非奇异变换:
(2)
ቤተ መጻሕፍቲ ባይዱ
(3)
使其状态空间表达式(1)化成: (4) 其中
(5)
称形如式(4)的状态空间表达式为能控标准 型。其中 , 为特征多项式:
使其状态空间表达式(13)化成: (15) 其中
(16)
(17) 称形如式(15)的状态空间表达式为能观标准 型。其中 是矩阵A的特征多项式的各项系数。
(18)
取变换阵
:
(19)
直接验证,或者用对偶原理来证明。证明过程如下: 首先构造 的对偶系统
然后写出对偶系统 表达式的能观标准 型即是
的能控标准 的能控标准 型,即
维不能观测的子系统,便得到一个 。
非奇异变换阵
是这样构成的,取
(14)
3.8.3
按能控性和能观性进行分解
1)如果线性系统是不完全能控和不完全能观的,若对该系统同时按能
控性和能观性进行分解,则可以把系统分解成能控且能观、能控不能观、
不能控能观、不能控不能观四部分。当然,并非所有系统都能分解成有这 四个部分的。
为简明起见,下面列举三个具有上述类型的二阶系统,对其能控性加以 剖析。 (3)
(4)
(5)
1)对于式(3)的系统,系统矩阵A为对角线型,其标量微分方程形式为: (6) (7) 2)对于式(4)的系统,系统矩阵A为约旦型,微分方程组为: (8) (9) 3)对于式(5)的系统,系统矩阵虽也为约旦型,但控制矩阵第二行的元素
却为0,其微分子方程组为:
(10) (11) 2.具有一般系统矩阵的多输入系统 系统的状态方程为:
(12)
3.2.2 直接从A与B判别系统的能控性
1.单输入系统
线性连续定常单输入系统:
其能控的充分必要条件是由 A、b 构成的能控性矩阵: (14)