传递函数矩阵最小实现方法
降阶法人们在设计复杂系统时,总是希望在构造系统之前用模拟计算机或数字计算机对所设计的系统进行仿真,以检查系统性能是否达到指标要求。给定严格真传递函数矩阵G(s),为寻找一个维数最小的(A,B,C),使C(sl - 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)与(A,B,C5之
间必代数等价,即两个最小实现之间由非奇异线性变换阵T使得式子
A =T」AT,
B =T J B,
C =CT 成立。
(3)传递函数矩阵G(s)的最小实现的维数为G(s)的次数n.,或G(s)的极点多项式的最高次数。
为了寻求传递函数矩阵的最小实现,就意味着要把系统中不能控和不能观测的状态变量消去而不至于影响系统的传递函数。求最小实现的方法有三种:
1、降阶法。根据给定的传递函数矩阵G(s),第一步先写出满足G(s)的能控型实现,第二步从中找出能观测子系统;或者第一步先写出满足G(s)的能观测型实现,第二步从中找出能控子系统,均可求得最小实现。
2、直接求取约当型最小实现的方法。若G(s)诸元容易分解为部分分式形式,运用直接求取约当型最小实现的方法是较为方便的。
3用汉克尔矩阵法求取最小实现的方法。
下面主要研究降阶法(先求能控型再求能观测子系统的方法)并举例说明。
先求能控型再求能观测子系统的方法设(px q)传递函数矩阵G(s),且p v q时,
优先采用本法。取出G(s)的第j列,记为G j(s),是u j至y(s)的传递函
数矩阵,有
由n ij (s)的诸系数确定C j ,这时
A j
「a j,0
-6,° III Pj i
n j 1
1j,0
"I" j ,
b p -
G j (s)=[g 1j (s)....g qj (s)]T
=[
( )
l qu (s)
/((sy q qj
(s)
记d j (s)为q ij (s), q qj (s)的最小公倍式,则
1
Gj(s)
=帀g ⑼叭创
n | n , _1 ,,,
设d j (s) = s a j ,n_i S a“s a j,。
则 m(s) = '-ij,n j
- 1ij,n j ^
s^^ J|| - ■- j,i^ ■- ij,o , i = 1,...q
在此d j (s)是q 个子系统传递函数的公共部分,由单输入-多输出系统的实现 可知,能用能控规范I 型的 A j 、b j 实现d j (s),
G j (s)的实现为
C j
i q
询
令j =1^1, p ,便可得G j (s)的实现为
b 2
n ;.p
A P -
C
q n = C1 C2 川C P 1
当p v q时,显见A、B、C的维数均较小,且有n: =n。上述实现一定能
|.Si (nqQ
)n J
U nn o
1(n-OQ
)n
」
j Su SU 」L 。
由 TT- S U U} SU SU i
_Si
n o
|
Q
,
1
n-n o
有 SU = I n
o
由 CT 4 =C U uj-lcu cu 」-c o o , 由 TA 宀鬥A U Ui]=[SAU SAU i]=B
[S i 」
:S i AU S i AU- ]A 21
B J SB L 可,
〔SB 一旦」’
于是由能控型化为能控能观测型的简化步骤可归结为: 有 C o = CU U i
,有 A 。= SAU
A Q
有 B o 二
SB
控,但不一点能观测,需要找出能观测部分,为此需要判别(A,C )的能观测行<
C
CA 若(A,C )能观测,贝U ( A,B,C )为最小实现;若rankQ ° = rank
.
= n 0
v n
+
k 一
则从Q o 中选出n o 个线性无关行,记为S;在附加(n - n °)个任意行(通常为单 位矩阵I n 的任意行),记为S i ,即
构造W 的非奇异变换阵T ,€ 引入变换x =Tx ,由能观测性的结构分解可知
其中能观测子系统(A,B o ,C o )即为所求的最小实现。
(瓦,瓦,C 0)有如下简化求法:记T‘为
I .构造S 阵(从Q o 中选出n o 个线性无关行);
Vi
+ +
=S , Vn°
卑
+ T
T
■
V
n
=S
_
AQ
B = T |B° 1
C = c r = PC
o~\
'A 21 AQ j
I BQ J
L
J
— _1
A =TAT
A」
_0
■-1
-6
0〕
1
-5
r 一01
000
〔0b2
一0
1
一
C>_-1
1 1
一
1
g^s)= 丄
S +1 1 | 1
—1 s +1 —1
_s+1」
g2(s)二
s+3 _ 1 s+2
| -1 _(s+2)(s+3) s_3_
_s+2」
d=1
b2
「1」
C2屮
_-3
1
-1
2?由SU =%,求出U阵;
3. 计算最小实现。A00 =SAU , 瓦二SB , C0 =CU。
由于S选择的任意性及求解U的任意性,最小实现不唯一,但最小实现的维数是唯一的,且系统都是能控能观测的。下面举例说明该法。
例1、已知传递函数矩阵G(s),求最小实现。
s 2 1
G(s)j s+1 s+3
s s+1
l-S+1 s + 2_j
解:化G(s)为严格真传递函数矩阵G(s)
J ___
G(s)=卜:1 s:3+£ 0L(?(S)+D
-1 -1 1 1
?s + 1 s + 2」
求G(s)的最小实现
令dMs) =s 1,d2(s) =s2 5s 6,其能控规范I型实现为
A 一-1
(?(s)的能控型实现为
(A,C)的能观测性判别:由于ran kC =2=m
「C l
~C 1 rankQ ° =rank | n m =rank |
=rank
]CA _
]CA
-1 -1
2 1
七
-1 =3 = n
-6 -3
即(A,C)能观测。(A,B,C )能控且能观测,即为G(s)的最小实现 G(s)的最小实
现为(A,B,C,D )
例2、求下列G(s)的最小实现维数及最小实现
-4s + 6
2s+3 ]
(s+1)(s + 2) (s+1)(s+3)
G(s
)=
c ” -2 -1 ] (s+1)(s + 2) (s + 1)(s + 2) _
解(1)确定最小实现维数n,..:所有G(s)的一阶子式的最小公分母为(s ,1)(s ?2); 二阶子式只有一个
0,其分母为任意常数。故所有子式的最小公分母仍为
(s 1)(s 2),有 n .=2。
(2)
g(s)二
1 :4s + 6]
(s+1)(s+2)] -2 一
心乔E 2:3
]
C 』3 2
_-1
A 0
「 th
C
b =
1
IL 0
A 2
_0 b 2
令d 1 (s) = d 2(s) = (s 1)(s 2),其能控规范I 型实现为
(A,C)的能观测性判别:由于rankC = 2二m 1 1
-3
7 4]
A
■ C C 21
_C 1 CA +
E-
1
=rank CA
+
CA njm
ran kQ 0 = ra nk =rank
■6 4
3 2〕 -2
-1 0 -8 -6 -4 -3 0 -2 0 -1 12 10
6 5
6
2
3
=2 <4
(A,C)不完全可观测 从Q o 中选出二行构成
S 阵亠 6
4
3 2
-2 0-10
由SU “2求U 阵: 6 4 IL-2 0 3
-1
U
11
2
] | U
21
U 31 U 41
1 0 「0 1」
四个方程含8个未知数,设任意规定
u
31 = U32 = u
41
=U 42 = 0,可解得
1 U = 一
4 0 2
3
4 0
故最小实现为 _3 一 2
A)= SAU = 1 2
4
_0
C 。=cu