本文仅供学习使用
本文参考:B站:DR_CAN
控制之美(卷1)
1. 绪论
通过研究动态系统的数学模型和系统表现,可以得到在给定输入
u
(
t
)
u(t)
u(t)作用下的系统响应(Response
,即系统在输入
u
(
t
)
u(t)
u(t)作用下的输出
x
(
t
)
x(t)
x(t))。当掌握了动态系统输入与输出的关系之后,就可以设计控制器来调节动态系统的输入,使得系统的输出按照预期的目标响应。
一般而言,控制系统(Control System)
由控制器(Controller)
和动态系统
组成。
控制器会根据参考值(Reference)
r
(
t
)
r(t)
r(t)来决定控制量,即动态系统的输入
u
(
t
)
u(t)
u(t)。这种简单的控制方式称为开环(Open Loop)控制
。当系统的全部信息可知且准确时,开环控制可以完美地达成控制目标。
如果希望精确地控制系统,则需要使用闭环(Closed Loop)控制系统
,它与开环控制的最大区别是,在闭环控制中会测量系统的输出,并将其反馈(Feedback)
到输入端与参考值进行比较考值与实际系统输出的差称为误差(Error)
,控制器将根据误差调整控制量。闭环控制系统可以实现高精度的控制,同时补偿由于外界扰动及系统建模不准确而引起的偏差。
2. 状态空间表达State-Space Representation
状态空间方程(State Space Model)
描述系统数学模型的方法。状态空间方程式现代控制理论的基础,它以矩阵的形式表达系统状态变量、输入及输出之间的关系。可以描述和处理多输入输出(Multiole Input Multiple Output,MIMO)
的系统。目前流行的一些算法,如模型预测控制、卡尔曼滤波器及最优控制,都是在状态空间方程的表达形式基础上发展而来的。
对于同样的系统,在现代控制理论中则会使用状态空间方程的表达状态空间方程是一个集合,它包含了系统的输入、输出及状态变量,并把它们用一系列的一阶微分方程表达出来。对于本例中的二阶系统,为了将其写成状态空间方程,需要选取合适的状态变量(State Variables)
,才能使二阶系统转化为一系列的一阶系统。
上述形式可推广并得到状态空间方程的一般形式,即:
{
d
z
(
t
)
d
t
=
A
z
(
t
)
+
B
u
(
t
)
y
(
t
)
=
C
z
(
t
)
+
D
u
(
t
)
\begin{cases} \frac{\mathrm{d}\boldsymbol{z}\left( t \right)}{\mathrm{d}t}=A\boldsymbol{z}\left( t \right) +B\boldsymbol{u}\left( t \right)\\ \boldsymbol{y}\left( t \right) =C\boldsymbol{z}\left( t \right) +D\boldsymbol{u}\left( t \right)\\ \end{cases}
{dtdz(t)=Az(t)+Bu(t)y(t)=Cz(t)+Du(t)
其中:
z
(
t
)
\boldsymbol{z}\left( t \right)
z(t)是状态变量,是一个
n
n
n维向量:
z
(
t
)
=
[
z
1
(
t
)
,
z
2
(
t
)
,
⋯
z
n
(
t
)
]
T
\boldsymbol{z}\left( t \right) =\left[ z_1\left( t \right) ,z_2\left( t \right) ,\cdots z_{\mathrm{n}}\left( t \right) \right] ^{\mathrm{T}}
z(t)=[z1(t),z2(t),⋯zn(t)]T
y
(
t
)
\boldsymbol{y}\left( t \right)
y(t)是系统输出,是一个
m
m
m维向量:
y
(
t
)
=
[
y
1
(
t
)
,
y
2
(
t
)
,
⋯
y
m
(
t
)
]
T
\boldsymbol{y}\left( t \right) =\left[ y_1\left( t \right) ,y_2\left( t \right) ,\cdots y_{\mathrm{m}}\left( t \right) \right] ^{\mathrm{T}}
y(t)=[y1(t),y2(t),⋯ym(t)]T
u
(
t
)
\boldsymbol{u}\left( t \right)
u(t)是系统输入,是一个
p
p
p维向量:
u
(
t
)
=
[
u
1
(
t
)
,
u
2
(
t
)
,
⋯
u
p
(
t
)
]
T
\boldsymbol{u}\left( t \right) =\left[ u_1\left( t \right) ,u_2\left( t \right) ,\cdots u_{\mathrm{p}}\left( t \right) \right] ^{\mathrm{T}}
u(t)=[u1(t),u2(t),⋯up(t)]T
这说明,当使用状态空间方程来描述系统时,有
n
n
n个状态变量、
m
m
m个输出和
p
p
p个输入。它可以表达多状态多输出、多输人的系统。其中,矩阵
A
A
A是
n
×
n
n×n
n×n矩阵,表示系统状态变量之间的关系,称为状态矩阵或者系统矩阵。矩阵
B
B
B是
n
×
p
n×p
n×p矩阵,表示输入对状态变量的影响,称为输入矩阵或者控制矩阵。矩阵
C
C
C是
m
×
n
m×n
m×n矩阵,表示系统的输出与系统状态变量之间的关系,称为输出矩阵。矩阵
D
D
D是
m
×
p
m×p
m×p矩阵,表示系统的输入直接作用在系统输出的部分,称为直接传递矩阵。
上述弹簧阻尼系统,因为它只有一个输入和一个输出,所以本系统属于单输入单输出(Single Input Single Output,SISO)系统
。
多输入多输出(Multiple Inputs Multiple Outputs,MIMO)系统
见CH06
2.1 状态空间方程与传递函数的关系
如果令
G
(
s
)
G(s)
G(s)的分母部分为零,即
∣
s
I
−
A
∣
=
0
|sI-A|=0
∣sI−A∣=0,得出的
s
s
s值有两个含义:第一,从传递函数的角度考虑,它是传递函数的极点;第二,从状态矩阵的角度考虑,它是矩阵
A
A
A的特征值(令
∣
s
I
−
A
∣
=
0
|sI-A|=0
∣sI−A∣=0是求矩阵
A
A
A特征值的公式)。通过分析传递函数极点可以判断系统的表现。而当把系统写成状态空间方程之后,状态矩阵
A
A
A的特征值即为其相对应的传递函数
G
(
s
)
G(s)
G(s)的极点。因此,通过分析矩阵
A
A
A的特征值也可以判断系统的表现。
2.2 状态空间方程的解——矩阵指数函数
3. Phase Portrait相图,相轨迹
3 1. 1-D
这两个点被称为平衡点(Equilibrium Point或者Fixed Point)
。
需要注意的是,只有当初始位置在
x
f
2
x_{f2}
xf2左边,即
x
(
0
)
<
x
f
2
x(0)<x_{f2}
x(0)<xf2时,状态变量才可以回到平衡点。因此,
x
f
1
x_{f1}
xf1是一个局部(Local)稳定
的平衡点。
3 2. 2-D
图2中,相轨迹是一条曲线,但仍然是远离平衡点位置的。在这种情况下,平衡点称为不稳定节点(Unstable Node)
图3中,平衡点称为鞍点(Saddle)
,是一个不稳定的点。
3 3. General Form
II 是以原点为中心的椭圆方程,这说明相图中的点沿椭圆的轨迹运行。在这种情况下,平衡点称为中心点(Center)
,相轨迹会围绕着这个中心点做圆周运动。判断相轨迹的运动方向,可以以一个在横轴的点进行分析尝试。发现
x
1
x_1
x1与
x
2
x_2
x2此消彼长,循环往复。如果从能量的角度来分析,它们的总能量取决于初始状态
x
(
0
)
x(0)
x(0)随着时间发生改变。所以这个系统在稳定与不稳定之间,
x
(
t
)
x(t)
x(t)始终有界但始终不为0。
会持续地振荡且振幅不断地加强。在这种情况下,平衡点称为不稳定焦点(Unstable Spiral)
可以判断出它的平衡点是一个稳定焦点(Stable Spiral)
。相轨迹将沿螺旋线指向原点。其所对应的时间函数将振荡衰减,直至为0。
3 4. Summary
上述例子说明,状态矩阵
A
A
A的特征值将决定平衡点的类型及系统的表现。通过上述分析可知,状态矩阵
A
A
A的特征值实部部分决定了平衡点的稳定性,而特征值的虚部部分决定了系统是否会有振动。
3.5. 爱情中的数学-Phase Portrait 相图动态系统分析
动态系统的分析可以分为三个步骤:
- 描述系统,即通过语言来描述系统的特性;
- 数学分析,即使用数学工具对系统进行量化解析;
-
结果与讨论,即根据第二步数学分析的结果,进行深层次的思考与讨论。
3.6 连续系统离散化
3.7 Summary
-
状态空间表达形式:
用一系列的一阶微分方程表达系统的输入、输出及状态变量。其一般形式为
状态矩阵的特征值与其所对应的单输人单输出传递函数的极点相同。极点决定系统的表现,状态矩阵的特征值是研究的重点 -
相平面与相轨迹分析:
利用矩阵特征值与特征向量的性质可以将矩阵对角化,从而将系统解耦。
利用图形的方法,可以直观地分析复杂的非线性一阶系统平衡点的表现。
二阶系统平衡点的性质由状态矩阵的特征值决定。 -
动态系统的分析方法:
第一步,描述系统。通过语言来描述系统的特性。
第二步,数学分析。使用数学工具对系统进行解析。
第三步,结果与讨论。分析的结果,进行深层次的思考与讨论。
4. 系统的可控性Controllability(LTI)线性时不变
时不变性是指如果系统的输入喜好延迟了时间
T
T
T, 那么系统的输出也会延迟时间
T
T
T。一般情况下,时不变系统的数学表达式中都是常数系统(系数不是时间的函数)。线性时不变系统必须同时满足如下的两个性质。
需要说明的是,从严格意义上讲,时不变系统是不存在的,因为不能两次踏进同一条河流”,但在大部分工程情况下,在系统分析的时间区间内,参数是恒定的或者是缓慢变化的。对于非线性的系统,一般可以做线性化处理(参考附录A)。以近似为线性时不变的系统不在本书的讨论范围之内。
5. 稳定性stability-李雅普诺夫Lyapunov
CH01-2 从直观上解释了稳定性的含义,下面将介绍严谨的数学定义。在1892年,俄国数学家亚历山大·李雅普诺夫(Aleksandr Lyapunov)在其博士论文《运动稳定性的一般问题》中提出了稳定性的科学概念。本书将选用这个概念来定义系统的稳定性,它是一个通用的概念,既可以运用在线性系统中,也可以运用到非线性的系统分析中。
Stability in the sense of Lyapunov
渐近稳定性Asymptotic Stability——比李雅普诺夫意义下的稳定性更加严格的一种稳定
如果平衡点不符合上面两种条件,则不稳定。
在经典控制理论中,稳定特指渐近稳定。李雅普诺夫意义下的稳定(即极点在虚轴上的情况)会被称为临界稳定或者不稳定。在这样的定义下,从经典控制理论的角度来看,
参数 δ ( t 0 ) δ(t_0) δ(t0)和 ε ε ε的物理意义可以这样来理解: ε ε ε是一个稳定性的指标,如果状态变量始终在 ε ε ε以内,那么平衡点就符合李雅普诺夫意义下的稳定。 δ ( t 0 ) δ(t_0) δ(t0)则是平衡点稳定的前提条件,或者说是稳定所能承受的最大干扰,也可以理解为收敛域。
根据上述分析,如果
δ
(
t
0
)
δ(t_0)
δ(t0)可以任意选择(选择无限大),那么就可以推断出平衡点是全局稳定(Global Stability)
,因为不管系统的初始状态在任何位置,它都可以最终收敛到一个范围之内(渐近稳定时收敛于0)。反之,如果
δ
(
t
0
)
δ(t_0)
δ(t0)的选择是有条件的,则平衡点是局部稳定(Local Stability)
。
5.1 稳定性与状态空间方程
使用状态空间方程描述线性时不变系统的一般表达式为:
d
z
(
t
)
d
t
=
A
z
(
t
)
+
B
u
(
t
)
y
(
t
)
=
C
z
(
t
)
+
D
u
(
t
)
\frac{\mathrm{d}z\left( t \right)}{\mathrm{d}t}=Az\left( t \right) +Bu\left( t \right) \\ y\left( t \right) =Cz\left( t \right) +Du\left( t \right)
dtdz(t)=Az(t)+Bu(t)y(t)=Cz(t)+Du(t)
其中,
z
(
t
)
z(t)
z(t)是状态变量,
y
(
t
)
y(t)
y(t)是系统的输出,
u
(
t
)
u(t)
u(t)是系统的输入,矩阵
A
A
A是状态矩阵,矩阵
B
B
B是输入矩阵,矩阵
C
C
C是输出矩阵,矩阵
D
D
D是直接传递矩阵。系统的输出
y
(
t
)
y(t)
y(t)是状态变量
z
(
t
)
z(t)
z(t)与输入
u
(
t
)
u(t)
u(t)的线性组合,如果系统的状态变量和输入都有界,那么系统的输出也是有界的(符合李雅普诺夫意义下的稳定性)。
首先考虑
0
0
0输入状态,即
u
(
t
)
=
0
u(t)=0
u(t)=0。考虑状态矩阵
A
A
A的特征值与平衡点类型之间的关系,如果
A
A
A是一个二维矩阵,那么其特征值与稳定
性的关系如下所示:
上述分析方法可以推广到更高维度的状态矩阵中。对于
n
n
n维向量
z
(
t
)
z(t)
z(t),首先令
z
(
t
)
=
P
z
ˉ
(
t
)
z(t)=P\bar{z}(t)
z(t)=Pzˉ(t),其中,
P
P
P是过渡矩阵,由矩阵
A
A
A的特征向量组成。可得
d
z
ˉ
(
t
)
d
t
=
[
λ
1
⋱
λ
n
]
z
ˉ
(
t
)
\frac{\mathrm{d}\bar{z}(t)}{\mathrm{d}t}=\left[ \begin{matrix} \lambda _1& & \\ & \ddots& \\ & & \lambda _{\mathrm{n}}\\ \end{matrix} \right] \bar{z}(t)
dtdzˉ(t)=
λ1⋱λn
zˉ(t)
其中
λ
i
\lambda _{\mathrm{i}}
λi表示A的特征值,
{
z
ˉ
1
(
t
)
=
C
1
e
λ
1
t
⋮
z
ˉ
n
(
t
)
=
C
n
e
λ
n
t
\begin{cases} \bar{z}_1(t)=C_1e^{\lambda _1t}\\ \vdots\\ \bar{z}_{\mathrm{n}}(t)=C_{\mathrm{n}}e^{\lambda _{\mathrm{n}}t}\\ \end{cases}
⎩
⎨
⎧zˉ1(t)=C1eλ1t⋮zˉn(t)=Cneλnt
原向量 z ( t ) z(t) z(t)是 z ˉ ( t ) \bar{z}(t) zˉ(t)的线性组合,所以它的收敛或发散与 z ˉ ( t ) \bar{z}(t) zˉ(t)一致。此,我们可以将二维状态变量平衡点稳定判定方法沿用到更高维度的情况下,得到关于状态空间方程稳定性的两个结论(考虑零输入系统的状态空间方程为): d z ( t ) d t = A z ( t ) \frac{\mathrm{d}z(t)}{\mathrm{d}t}=Az\left( t \right) dtdz(t)=Az(t)
- 如果 A A A的特征值的实部都不大于0,它的平衡点将符合李雅普诺夫意义下的稳定性。
- 如果A的特征值的实部都小于0,它的平衡点将符合渐近稳定性。
上述结论为系统的控制器设计提供了思路。在闭环控制系统中,输入
u
(
t
)
u(t)
u(t)是状态变量的一个函数,例如
u
(
t
)
=
−
K
z
(
t
)
u(t)=-Kz(t)
u(t)=−Kz(t),其中矩阵
K
K
K是控制矩阵。
d
z
(
t
)
d
t
=
A
z
(
t
)
−
B
K
z
(
t
)
=
(
A
−
B
K
)
z
(
t
)
=
A
c
l
z
(
t
)
\frac{\mathrm{d}z(t)}{\mathrm{d}t}=Az\left( t \right) -BKz\left( t \right) =\left( A-BK \right) z\left( t \right) =A_{\mathrm{cl}}z\left( t \right)
dtdz(t)=Az(t)−BKz(t)=(A−BK)z(t)=Aclz(t)
其中,
A
c
l
=
A
−
B
K
A_{cl}=A-BK
Acl=A−BK,代表闭环控制系统的状态矩阵。根据上述分析,如果希望得到一个渐近稳定的系统,则需要设计合适的控制矩阵
K
K
K,使得
A
c
l
A_{cl}
Acl的特征值实部都为负数。
5.2 Summary
-
稳定性的定义:
李雅普诺夫意义下的稳定性:随着时间的增加,系统的状态变量始终在平衡点附近移动(有界)。
渐近稳定性:随着时间的增加,系统的状态变量会最终收敛于平衡点。 -
稳定性与传递函数:
在研究系统稳定性的时候考虑单位冲激响应,即传递函数本身。
系统稳定的条件是系统的闭环传递函数的极点均在复平面的左半部分(实部小于0)。此时系统满足有界输入有界输出稳定。 -
稳定性与状态空间方程: 状态矩阵的特征值的实部决定了系统的稳定性:
当实部都不大于0时,系统符合李雅普诺夫意义下的稳定性。
当实部都小于0时,系统符合渐近稳定性。
6. 线性控制器设计Linear Controller Design
7. LQR控制器 Linear Quadratic Regulator
线性控制器设计-轨迹跟踪(Fellow a Desired Path)
8. 状态观测器设计 Linear Observer Design
9. 可观测性与分离原理
10. 现代控制理论串讲
11. 非线性控制理论基础 Lyapunov直接方法
12. Invariance Princilpe-LaSalle;s Theorem不变性原理
13. Nonlinear Basic Feedback Stabilization 非线性系统稳定设计
14. Nonlinear Backstepping Control反馈线性化控制
Nonlinear Backstepping Controller 补充习题 - Backsteeping controller Design Review
15. NonlinearAdaptive Controller非线性自适应控制器
16. Robust Controller非线性鲁棒控制器
Robust Control : tp achieve robust performance and/or stability in percence
of bounded
modeling error
16.1 Slide Control 滑膜控制
16.2 High Gain + High Frequency
16.3 三种鲁棒控制器的比较+如何分析控制器
文章来源:https://www.toymoban.com/news/detail-800587.html
17. 三种鲁棒控制器的比较+如何分析控制器
文章来源地址https://www.toymoban.com/news/detail-800587.html
到了这里,关于[足式机器人]Part2 Dr. CAN学习笔记-Ch04 Advanced控制理论的文章就介绍完了。如果您还想了解更多内容,请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章,希望大家以后多多支持TOY模板网!