在前面各章中,我们总是假设已经知道了受控对象的模型,但由于实际中存在种种不确定因素,如:
参数变化;
未建模动态特性;
平衡点的变化;
传感器噪声;
不可预测的干扰输入;
等等,所以我们所建立的对象模型只能是实际物理系统的不精确的表示。鲁棒系统设计的目标就是要在模型不精确和存在其他变化因素的条件下,使系统仍能保持预期的性能。如果模型的变化和模型的不精确不影响系统的稳定性和其它动态性能,这样的系统我们称它为 鲁棒控制系统 。
鲁棒性 (Robustness)
所谓鲁棒性,是指标称系统所具有的某一种性能品质对于具有不确定性的系统集的所有成员均成立,如果所关心的是系统的稳定性,那么就称该系统具有鲁棒稳定性;如果所关心的是用干扰抑制性能或用其他性能准则来描述的品质,那么就称该系统具有鲁棒性能。
系统的不确定性
参数不确定性,如二阶系统:
可以代表带阻尼的弹簧装置,RLC电路等。这种不确定性通常不会改变系统的结构和阶次。
动态不确定性也称未建模动态,我们通常并不知道它的结构、
阶次,但可以通过频响实验测出其幅值界限:
],[,11)( 2 aaaasssG
为确定函数)(,,)()( jWRjWj
)(s?
加性不确定性:
乘性不确定性:
)()(),( 0 ssGsG
)())((),( 0 sGsIsG
一个例子设汽车质量为 M,路面摩擦系数为?,汽车的力学模型如下图所示:
其运动方程为:
如果考虑到汽车的质量 M随车载负荷发生变化,且?也随路面状况不同而变化,则方程的系数就具有一定的不确定性,即,无法得到 M和?的精确值。假设 M和?的取值范围给定如下:
为给定常数i
MMM


,2020
1010


M
fvdtdvM
v
f?v
那么实际的被控对象就可以描述为如果用 f 到 v的传递函数来描述,则有其中可以找到适当的界函数
2100,,)() Mfvdt
dvMM +(
)()()( 1)( 0
00
ssGsMMsG
)]()[()(
)(
,
1
)(
0000
00
0




sMMsM
Ms
s
sM
sG
)()(),( jWjjW有鲁棒控制理论是分析和处理具有不确定性系统的控制理论,包括两大类问题:鲁棒性分析及鲁棒性综合问题。 鲁棒性分析 是根据给定的标称系统和不确定性集合,找出保证系统鲁棒性所需的条件;而 鲁棒性综合 (鲁棒控制器设计问题)就是根据给定的标称模型和不确定性集合,基于鲁棒性分析得到的结果来设计一个控制器,使得闭环系统满足期望的性能要求。
主要的鲁棒控制理论有:
Kharitonov区间理论;
H?控制理论;
结构奇异值理论 (?理论 );
等。
Kharitonov定理具有不确定参数的系统假设系统的特征多项式为其系数满足我们称 (1)为区间多项式,为了判定系统的稳定性,应该研究所有可能的参数组合,这是个无穷检验问题。
前苏联数学家 Kharitonov于 1978年给出了关于判断区间多项式族鲁棒稳定性的四多项式定理,为研究参数不确定系统的鲁棒性分析奠定了基础。
( 1 ) )( 0111 asasasasf nnnn
],[0,,,1,0, iiiii aaniaaa?
Kharitonov定理,(1)中的每一个多项式均稳定当且仅当下面的四个多项式稳定注,定理中的四个多项式通常被称作 Kharitonov顶点多项式。 Kharitonov定理的意义在于它将区间多项式中无穷多个多项式的稳定性与四个定点的稳定性等价起来,
将无穷检验变为有限检验(顶点检验)。








5
5
4
4
3
3
2
2104
5
5
4
4
3
3
2
2103
5
5
4
4
3
3
2
2102
5
5
4
4
3
3
2
2101
)(
)(
)(
)(
sasasasasaasP
sasasasasaasP
sasasasasaasP
sasasasasaasP
考虑下图所示的闭环系统其中


m
i
iii
in
i
i
i
i rrrsrrsDsqsNrsD
sNrsG
0 0
,,),(,)(,),( )(),(
闭环传递函数为
),(1
),(),(
rskG
rsGrsG
CL
Gcl(s)的分母为 )(),( skNrsD?
- G(s)
k
u y
例:
1
122),(
1
2
2
3
3
4
23


srsrsrs
sssrsG
]3,2[],4,3[,54 321 rrr,
取 k= 1,此时闭环传递函数的分母为
21221 12233423122334 spspspsssssrsrsrs
其中 ]4,3[],6,5[],3,2[ 321 ppp
此时上面的闭环系统稳定当且仅当下面的四个多项式稳定
432
4
432
3
432
2
432
1
4622)(
4522)(
3632)(
3532)(
sssssF
sssssF
sssssF
sssssF




H?控制理论
H?控制理论提出的背景现代控制理论的许多成果在理论上很漂亮,但实际应用并不成功。主要原因是忽略了对象的不确定性,并对系统所存在的干扰信号作了苛刻的要求。
加拿大学者 Zames在 1981年提出了著名的 H?控制思想,考虑如下一个单输入单输出系统的设计问题:对于属于一个有限能量的干扰信号,设计一个控制器使得闭环系统稳定且干扰对系统期望输出影响最小。由于传递函数的 H?范数可描述有限输入能量到输出能量的最大增益,所以用表示上述影响的传递函数的 H?范数作为目标函数对系统进行优化设计,就可使具有有限功率谱的干扰对系统期望输出的影响最小。
对于反馈系统如果 P(s)具有误差,那么相应地开环和闭环频率特性也具有误差
)()(1
)()( ),()()(


jKjP
jKjPGjKjPjG
BK
其中 K(s)为控制器,w为干扰信号,r为参考输入,u
为控制输入,e为控制误差信号,y为输出信号。系统的开环和闭环频率特性为
)()()( 0 sPsPsP
-
r y
P(s)kK(s)
e
w
u
)()()(
)()()(
0
0


jGjGjG
jGjGjG
BBB
KKK


其中体现了开环特性的相对偏差 到闭环频率特性的增益,因此,如果我们在设计控制器 K时,
能够使 S的增益足够小,即
)()(1
)()( ),()()(
0
0
000

jKjP
jKjPGjKjPjG
BK
分别为开环和闭环频率特性的标称函数,简单的推导可得而传递函数
)(
)(
)()(1
1
)(
)(
0?

jG
jG
jKjPjG
jG
K
K
B
B?

)()(1
1 )(
0 sKsP
sS
KK GG?
BB GG?
为充分小正数,)(?jS
那么闭环特性的偏差将会抑制在工程允许的范围内。
传递函数 S(s)称为系统的 灵敏度函数 。实际上 S(s)还等于干扰 w到输出的闭环传递函数,因此减小 S(s)的增益就等价于减小干扰对控制误差的影响。引入定义
)]([s u p)(
jSsS
R
其中 表示最大奇异值,即)(,)}({)( 21*m a x AAA
H?控制问题即为对于给定的?> 0,设计控制器 K
使得闭环系统稳定且满足
)(sS
为最大特征值。的共轭转置阵,为 m a x*?AA
H?理论中考虑干扰信号是不确定的,而是属于一个可描述集
L2中包含的是能量有限的信号。考虑抑制干扰 w?L2对系统性能的影响,为此引入表示干扰抑制水准的标量?,
求控制器 K使得满足
z为输出信号。定义
0 22 )(|)( dttwtwL
( 1 ),222222 Lwwz
2
2
0
s u p)( wzjT
w
zw


其中 Tzw(s)为由 w至 z的闭环传递函数,则 (1)等价于求使?最小的控制器 K就是 H?最优设计问题。
)( jT zw
传递函数的 H?范数对于系统的传递函数 G (s),若其在右半平面无极点,定义下面的范数为 H?范数


djuu
u
Gu
sG
2
2
2
2
)(
2
1
,s u p)(
定理:
)]([s u p)(
jGsG
R
其中闭环系统鲁棒稳定性分析其中?(s)为任意稳定的真有理分式且满足 ||?(s)||1
定理,上图所示的闭环系统对任意的?(s)均稳定当且仅当
1))()()(( 1sKsGIsK
-
u y
G(s)kK(s)
G(s)?
加性不确定性考虑下图所示系统闭环系统鲁棒稳定性分析其中?(s)为任意稳定的真有理分式且满足 ||?(s)||1
定理,上图所示的闭环系统对任意的?(s)均稳定当且仅当
1)()())()(( 1 sKsGsKsGI
-
u y
G(s)kK(s)
G(s)?
乘性不确定性考虑下图所示系统时域模型的鲁棒性考虑系统其中,?为任意满足 的实矩阵,E,F为已知矩阵。
定理,对任意的?,,A+E?F 稳定当且仅当
xAAx )(
1)( 1 EAsIF控制系统的鲁棒控制问题即为设计 K使得 A+BK+ E?F稳定,也即
KxuBuxFEAx,)(?
1)( 1 EBKAsIF
FEA IT
IT
实 验
Furuta摆实验三自由度直升机系统