定理 4-4 若 (A,B)可控,则对 B值域中的任一非零向量 b,均存在一个状态反馈增益阵 K1,使得 (A+BK1,b)可控。
证明
1,引理,(见习题 4-3)
若 (A,B)可控,可选取向量 u1,u2,…,un-1,使得由下式定义的 n
个向量 x1,x2,…,xn线性无关 。
x1=b,xk+1=Axk+Buk (k=1,2,…,n-1) (S—1)
习题 4-3的证明用归纳法 。 k=1 x1=b
若存在 使 线性无关,
121,,,?kuuu? )(,,,21 nkxxx k
kxxx,,,21?要证:存在,使 与 线性无关,
即 不属于,此处 是所张成的空间 。
kkk BuAxx 1ku
kLkL1?kx
状态反馈所以 包含了,
用反证法。
,ku?
kkk BuAxx 1
kL
ku
kk Axx1 kL kkk BuAxx 1
,ku?
kkk LAxx 1 kLB?Im
kL
)( kixi?
kiLBuxAxx kiiik,,2,111
kL
kL BA|
BA|
均属于,特别的取 =0,
属于 。由于 ( ),因为
,所以 。又因为 中的任一向量有 。 是 A的不变子且包含了 B的值域,空间,这与的维数为 n相矛盾。
1K
2,定义矩阵
0121211 nn uuuxxxK
1211211 0 nn xxxuuuK (S—2)
nxxx,,,21? 0,,,,121?nuuu?1K
映射为 p维向量将 n维向量是 p× n矩阵
,
111221111 nn uxKuxKuxK?
)( 1 bBKA?3,证明 可控
bBKAxBKABuAxx
bBKAxBKABuAxx
bBKAxBKAxBKAxBuAxx
bBKAxBKAbBuAbx
bx
n
nnnn
n
nnnn
1
11111
2
121221
2
121212223
11112
1
)()(
)()(
)()(
)(
nbBKAbBKAbBKAbr a n k n 11211 )()()(?
nxxx,,,21?因为 线性无关,即所以取 K= K1+Lk,即可证明定理 4-5。
A+BK1+bk=A+BK1+BLk= A+B(K1+Lk)
定理 4-5 若系统 (4-1)可控,则存在状态反馈增益阵 K,使得
A+BK的 n个特征值配置到复平面上 n个任意给定的位置 ( 复数共軛成对出现 ) 。
证明 首先选取非零向量 L,可得 b=BL,由定理 4-4可知存在
K1,使 (A+BK1 b)可控 。 由单变量极点配置定理可知存在 n维行向量 k,使得 A+BK1+bk 的特征值可任意配置,由于
TLuxx 11
10
01
00
100
010
011
例题 1 系统方程为试构造 K1,使 (A+BK1,b=BL)可控。
解 取 x1 =BL=b,由
x1=b,xk+1=Axk+Buk (k=1,2,…,n-1)
因为 Ax1与 x1线性无关,故取 x2= Ax1,可得 u1=[00]T。 又因为 A x2
与 x1 x2构成线性相关组,u2不能取 [00]T,可取 u2=[-1 1]T,这样可得 x3= Ax2+B u2= [2 0 2] T。 由的计算式 (S—2)可得
101
101
110
)()(
111
111
011
111
111
211
011
210
010
010
2
11
1
1
1
bBKAbBKAb
BKA
K
不难验证 (A+BK1 b)可控。
例题 2 系统方程为
uxx
10
01
00
00
1000
0100
0100
0010
欲使闭环系统 (A+BK)具有特征值 -2,-2,-1+j,-1-j,试确定状态反馈增益阵 K。
解 取 L=[1 0]T,x1=b1; 取 u1=[-10]T,可得 x2=[0100]T ; u2=[00] T,
可得 x3=[1000] T ;
取 u3=[01] T,可得 x4=[0001] T ; 于是由 的计算式可得
1K
0001
0100
1000
0001
0010
0100
0100
0001
1K
显然,(A+BK1,b1)可控。令 k=[k1 k2 k3 k4],直接计算
1001
0100
0010
4321
11 kkkkkbBKA
它的特征式为 s4-(1+k3)s3+(k3-k2)s2+(k2-k1)s+k1-k4,
期望特征式为 s4+6s3+14s2+16s+8,
比较上述两多项式的系数,可得
k1 = -37,k2= -21,k3= -7,k4= - 45
状态反馈阵可取为
0001
4582137
1 LkKK
在上面的做法中,在 L和 ui取定后,k就唯一的确定了 。 但 L
和 ui是非唯一的,这一事实至少可以说明达到同样极点配置的 K值有许多,K的这种非唯一性是多输入系统与单输入系统极点配置问题主要区别之一 。 如何充分利用 K的自由参数,
以满足系统其它性能的要求,是多输入系统状态反馈设计的一个活跃的研究领域 。
多输入系统状态反馈配置极点问题的另一特点是,非线性方程,,说明如下:如将 K阵的元素用待定系数 kij表示,闭环的多项式可以写为
)()()()()](de t [ 12211 KfsKfsKfsKfsBKAsI nnnnn
式中 f i (K)表示某一个以 K的元素 kij为变量的非线性函数 。 如果将期望多项式表成
nnnnn ssss 12211?
比较两式的系数,可知应有上式在单输入情况始终是线性方程组,在多输入时,一般是非线性方程。定理 4-4所提供的事实表明:当系统可控时,可以通过牺牲 K的自由参数,可以使 (S—3)简化为一组能解出的线性方程组。对例题 2,也可以用求解上述方程来做,通过 K
中自由参数的适当选取,往往可以方便地求出需要的 K阵。
),,2,1()( niKf ii (S—3)
例题 3 对例题 2中的系统,用直接求解 (S—3)式的方法,计算达到极点配置的 K阵 。
解 因为这时
8765
4321
1
1
0100
0010
kkkk
kkkk
BKA
方案1:取 k4=k5=k6=k7=0,1+k8= -2,由
464]1)1) [ (2( 232 sssss
易得 416,4
321 kkk
即有
3000 0564K
方案 2,取 k1=k2= 0,k3= -1,k4=1 由
816146]1)1[()2( 23422 ssssss
可得 k5= -8,k6= -16,k7= -14,k8= -7,即有
714168
1100K
以上的做法中,充分利用了将矩阵分块和相伴标准形的有关知识,从而方便了计算 。
如同单输入系统一样,定理 4-4中可控条件对于任意配置极点是充分必要条件,但对于某一组指定的特征值进行配置时,系统可控只是充分条件,而不是必要条件。给定极点组可用状态反馈达到配置的充分必要条件是给定极点组需包含系统的不可控模态。因此判别原来系统的模态可控性就成了关键。
例题 4 系统动态方程为
uxx
0
1
1
210
201
100
xy 210
3,2,14,3,2给定两组极点,分别为,和,问哪组极点可用状态反馈进行配置。
35.05.0 j解 计算出 A阵的特征值,分别为 –1,,可验证
–1 是不可控的,其它两个特征值是可控的。极点组 {-1,-2,-3}包含了不可控模态 -1,所以可用状态反馈进行配置;极点组
{-2,-3,-4}则不能达到配置。
如前所述,积分器解耦系统是不稳定的 。 还需在此基础上进一步改善控制律,这一问题由于比较复杂,现省略 。 下面介绍一种简单情况 。
npddd p21
定理 2-4 若系统可用状态反馈解耦,且则采用状态反馈
,
)(
)(
)(
,
01
1
20212
1
2
10111
1
1
11
2
2
2
1
1
1
pp
d
pd
d
p
d
d
d
d
d
d
kAkAkAc
IkAkAkAc
IkAkAkAc
EKEH
p
p
p?
可以将闭环传函矩阵化为
]1[)(
01
1
1
1
ii
d
id
d
id
df ksksksksd i a gsG i
i
i
i
i
(2-24)
解耦其中 kij是可调参数,可用来对闭环传递函数矩阵的对角元进行极点配置 。 由 (2-24)式可知 Gf(s)的麦克米伦阶为 n,说明这时解耦状态反馈律 ( 2-23) 未改变系统的可观测性 。
例题 2-8 系统动态方程为
xy
uxx
101
001
10
11
10
121
130
102
问可否用状态反馈律 u=Kx+Hv,将闭环传递函数阵变为
13
10
0
1
1
)(
2 ss
ssG
f
如有可能求出 K和 H。
解 c1B=[0 -1],d1=0; c2B=[0 0],c2AB=[2 1],d2=1,
12 10
2
1
E
EE d1+d2 +p=0+1+2=3=n
由 ( 2-2 4) 可知状态反馈律应按 ( 2-23 ) 选取
103
5.785.10
)3(
)(
02
11
2
1
2
2
11
1
IAAc
IAc
EK
EH
稳态误差定理 2-5 系统 (2-28)可控的的充分必要条件为系统 (2-25)可控,且注意:条件 (2-29)式只有当输出数至多等于输入数 (p ≥q)
且 rankC=q时才有可能 。
qnC BAr a n k
0
(2-29)
0
0
sIC
BAsI (2-30)
当 s不等于零时,(2-30)式中的矩阵,由于 ( A B) 可控,它的从上往下数的 n行是线性无关的,并且下面的 q行由于 s不是零,
它和前述 n行也是线性无关的,这时 (2-30)矩阵的秩为 n+q。 当
s=0时,由于 (2-29)式它的秩也是 n+q。 因此,对任意的 s,(2-30)
证明 考虑矩阵式矩阵的秩均为 n+q,故 (2-28)式的系统可控。反之,同样易于说明。
条件 (2-29) 式的说明,若 (2-29)式不满足,即有
qnsIC BAsIr a n k
s
00
0
这表明增广系统 (2-28)式有 s=0这一特征值不可控,状态反馈 (2-
31)不能改变这个特征值,它在闭环动态方程中仍存在,从而导致闭环系统不稳定。但在闭环传递函数形成时,产生了零、极对消,闭环传函无此 s=0的极点,这一零、极对消的原因是正向通道中引入了积分环节 (参考本节插图 ),也就是引入了一个
s=0 的极点,它一定是与对象的零点发生了相消。因此 (2-29)
式不成立意味着对象有 s=0的零点。在单变量的情况,就是对象传递函数的零点,但是在系统迥路中不稳定的零、极对消是不允许的。以上说明的一个例子如下
xyduxx 10,1032 10
试问是否可以用本节的方法设计控制器,使得在 yr(t)=1(t)作用下无稳态误差。
例题 2-9 系统方程为解 验证 (2-29)式
32
010
132
010
qnr a n k
计算对象部分的传递函数
23)( 2 ss
ssG
对这一系统无法用本节的方法进行控制器设计以消去稳态误差。
证明
1,引理,(见习题 4-3)
若 (A,B)可控,可选取向量 u1,u2,…,un-1,使得由下式定义的 n
个向量 x1,x2,…,xn线性无关 。
x1=b,xk+1=Axk+Buk (k=1,2,…,n-1) (S—1)
习题 4-3的证明用归纳法 。 k=1 x1=b
若存在 使 线性无关,
121,,,?kuuu? )(,,,21 nkxxx k
kxxx,,,21?要证:存在,使 与 线性无关,
即 不属于,此处 是所张成的空间 。
kkk BuAxx 1ku
kLkL1?kx
状态反馈所以 包含了,
用反证法。
,ku?
kkk BuAxx 1
kL
ku
kk Axx1 kL kkk BuAxx 1
,ku?
kkk LAxx 1 kLB?Im
kL
)( kixi?
kiLBuxAxx kiiik,,2,111
kL
kL BA|
BA|
均属于,特别的取 =0,
属于 。由于 ( ),因为
,所以 。又因为 中的任一向量有 。 是 A的不变子且包含了 B的值域,空间,这与的维数为 n相矛盾。
1K
2,定义矩阵
0121211 nn uuuxxxK
1211211 0 nn xxxuuuK (S—2)
nxxx,,,21? 0,,,,121?nuuu?1K
映射为 p维向量将 n维向量是 p× n矩阵
,
111221111 nn uxKuxKuxK?
)( 1 bBKA?3,证明 可控
bBKAxBKABuAxx
bBKAxBKABuAxx
bBKAxBKAxBKAxBuAxx
bBKAxBKAbBuAbx
bx
n
nnnn
n
nnnn
1
11111
2
121221
2
121212223
11112
1
)()(
)()(
)()(
)(
nbBKAbBKAbBKAbr a n k n 11211 )()()(?
nxxx,,,21?因为 线性无关,即所以取 K= K1+Lk,即可证明定理 4-5。
A+BK1+bk=A+BK1+BLk= A+B(K1+Lk)
定理 4-5 若系统 (4-1)可控,则存在状态反馈增益阵 K,使得
A+BK的 n个特征值配置到复平面上 n个任意给定的位置 ( 复数共軛成对出现 ) 。
证明 首先选取非零向量 L,可得 b=BL,由定理 4-4可知存在
K1,使 (A+BK1 b)可控 。 由单变量极点配置定理可知存在 n维行向量 k,使得 A+BK1+bk 的特征值可任意配置,由于
TLuxx 11
10
01
00
100
010
011
例题 1 系统方程为试构造 K1,使 (A+BK1,b=BL)可控。
解 取 x1 =BL=b,由
x1=b,xk+1=Axk+Buk (k=1,2,…,n-1)
因为 Ax1与 x1线性无关,故取 x2= Ax1,可得 u1=[00]T。 又因为 A x2
与 x1 x2构成线性相关组,u2不能取 [00]T,可取 u2=[-1 1]T,这样可得 x3= Ax2+B u2= [2 0 2] T。 由的计算式 (S—2)可得
101
101
110
)()(
111
111
011
111
111
211
011
210
010
010
2
11
1
1
1
bBKAbBKAb
BKA
K
不难验证 (A+BK1 b)可控。
例题 2 系统方程为
uxx
10
01
00
00
1000
0100
0100
0010
欲使闭环系统 (A+BK)具有特征值 -2,-2,-1+j,-1-j,试确定状态反馈增益阵 K。
解 取 L=[1 0]T,x1=b1; 取 u1=[-10]T,可得 x2=[0100]T ; u2=[00] T,
可得 x3=[1000] T ;
取 u3=[01] T,可得 x4=[0001] T ; 于是由 的计算式可得
1K
0001
0100
1000
0001
0010
0100
0100
0001
1K
显然,(A+BK1,b1)可控。令 k=[k1 k2 k3 k4],直接计算
1001
0100
0010
4321
11 kkkkkbBKA
它的特征式为 s4-(1+k3)s3+(k3-k2)s2+(k2-k1)s+k1-k4,
期望特征式为 s4+6s3+14s2+16s+8,
比较上述两多项式的系数,可得
k1 = -37,k2= -21,k3= -7,k4= - 45
状态反馈阵可取为
0001
4582137
1 LkKK
在上面的做法中,在 L和 ui取定后,k就唯一的确定了 。 但 L
和 ui是非唯一的,这一事实至少可以说明达到同样极点配置的 K值有许多,K的这种非唯一性是多输入系统与单输入系统极点配置问题主要区别之一 。 如何充分利用 K的自由参数,
以满足系统其它性能的要求,是多输入系统状态反馈设计的一个活跃的研究领域 。
多输入系统状态反馈配置极点问题的另一特点是,非线性方程,,说明如下:如将 K阵的元素用待定系数 kij表示,闭环的多项式可以写为
)()()()()](de t [ 12211 KfsKfsKfsKfsBKAsI nnnnn
式中 f i (K)表示某一个以 K的元素 kij为变量的非线性函数 。 如果将期望多项式表成
nnnnn ssss 12211?
比较两式的系数,可知应有上式在单输入情况始终是线性方程组,在多输入时,一般是非线性方程。定理 4-4所提供的事实表明:当系统可控时,可以通过牺牲 K的自由参数,可以使 (S—3)简化为一组能解出的线性方程组。对例题 2,也可以用求解上述方程来做,通过 K
中自由参数的适当选取,往往可以方便地求出需要的 K阵。
),,2,1()( niKf ii (S—3)
例题 3 对例题 2中的系统,用直接求解 (S—3)式的方法,计算达到极点配置的 K阵 。
解 因为这时
8765
4321
1
1
0100
0010
kkkk
kkkk
BKA
方案1:取 k4=k5=k6=k7=0,1+k8= -2,由
464]1)1) [ (2( 232 sssss
易得 416,4
321 kkk
即有
3000 0564K
方案 2,取 k1=k2= 0,k3= -1,k4=1 由
816146]1)1[()2( 23422 ssssss
可得 k5= -8,k6= -16,k7= -14,k8= -7,即有
714168
1100K
以上的做法中,充分利用了将矩阵分块和相伴标准形的有关知识,从而方便了计算 。
如同单输入系统一样,定理 4-4中可控条件对于任意配置极点是充分必要条件,但对于某一组指定的特征值进行配置时,系统可控只是充分条件,而不是必要条件。给定极点组可用状态反馈达到配置的充分必要条件是给定极点组需包含系统的不可控模态。因此判别原来系统的模态可控性就成了关键。
例题 4 系统动态方程为
uxx
0
1
1
210
201
100
xy 210
3,2,14,3,2给定两组极点,分别为,和,问哪组极点可用状态反馈进行配置。
35.05.0 j解 计算出 A阵的特征值,分别为 –1,,可验证
–1 是不可控的,其它两个特征值是可控的。极点组 {-1,-2,-3}包含了不可控模态 -1,所以可用状态反馈进行配置;极点组
{-2,-3,-4}则不能达到配置。
如前所述,积分器解耦系统是不稳定的 。 还需在此基础上进一步改善控制律,这一问题由于比较复杂,现省略 。 下面介绍一种简单情况 。
npddd p21
定理 2-4 若系统可用状态反馈解耦,且则采用状态反馈
,
)(
)(
)(
,
01
1
20212
1
2
10111
1
1
11
2
2
2
1
1
1
pp
d
pd
d
p
d
d
d
d
d
d
kAkAkAc
IkAkAkAc
IkAkAkAc
EKEH
p
p
p?
可以将闭环传函矩阵化为
]1[)(
01
1
1
1
ii
d
id
d
id
df ksksksksd i a gsG i
i
i
i
i
(2-24)
解耦其中 kij是可调参数,可用来对闭环传递函数矩阵的对角元进行极点配置 。 由 (2-24)式可知 Gf(s)的麦克米伦阶为 n,说明这时解耦状态反馈律 ( 2-23) 未改变系统的可观测性 。
例题 2-8 系统动态方程为
xy
uxx
101
001
10
11
10
121
130
102
问可否用状态反馈律 u=Kx+Hv,将闭环传递函数阵变为
13
10
0
1
1
)(
2 ss
ssG
f
如有可能求出 K和 H。
解 c1B=[0 -1],d1=0; c2B=[0 0],c2AB=[2 1],d2=1,
12 10
2
1
E
EE d1+d2 +p=0+1+2=3=n
由 ( 2-2 4) 可知状态反馈律应按 ( 2-23 ) 选取
103
5.785.10
)3(
)(
02
11
2
1
2
2
11
1
IAAc
IAc
EK
EH
稳态误差定理 2-5 系统 (2-28)可控的的充分必要条件为系统 (2-25)可控,且注意:条件 (2-29)式只有当输出数至多等于输入数 (p ≥q)
且 rankC=q时才有可能 。
qnC BAr a n k
0
(2-29)
0
0
sIC
BAsI (2-30)
当 s不等于零时,(2-30)式中的矩阵,由于 ( A B) 可控,它的从上往下数的 n行是线性无关的,并且下面的 q行由于 s不是零,
它和前述 n行也是线性无关的,这时 (2-30)矩阵的秩为 n+q。 当
s=0时,由于 (2-29)式它的秩也是 n+q。 因此,对任意的 s,(2-30)
证明 考虑矩阵式矩阵的秩均为 n+q,故 (2-28)式的系统可控。反之,同样易于说明。
条件 (2-29) 式的说明,若 (2-29)式不满足,即有
qnsIC BAsIr a n k
s
00
0
这表明增广系统 (2-28)式有 s=0这一特征值不可控,状态反馈 (2-
31)不能改变这个特征值,它在闭环动态方程中仍存在,从而导致闭环系统不稳定。但在闭环传递函数形成时,产生了零、极对消,闭环传函无此 s=0的极点,这一零、极对消的原因是正向通道中引入了积分环节 (参考本节插图 ),也就是引入了一个
s=0 的极点,它一定是与对象的零点发生了相消。因此 (2-29)
式不成立意味着对象有 s=0的零点。在单变量的情况,就是对象传递函数的零点,但是在系统迥路中不稳定的零、极对消是不允许的。以上说明的一个例子如下
xyduxx 10,1032 10
试问是否可以用本节的方法设计控制器,使得在 yr(t)=1(t)作用下无稳态误差。
例题 2-9 系统方程为解 验证 (2-29)式
32
010
132
010
qnr a n k
计算对象部分的传递函数
23)( 2 ss
ssG
对这一系统无法用本节的方法进行控制器设计以消去稳态误差。