滑模面设计及其应用实例详解
1. 滑模面的参数设计
对于线性系统
x
˙
=
A
x
+
b
u
,
x
∈
R
n
,
u
∈
R
\dot x = Ax + bu ,x \in R^n, u \in R
x˙=Ax+bu,x∈Rn,u∈R
滑模面可设计为
s
(
x
)
=
C
x
=
∑
i
=
1
n
c
i
x
=
∑
i
=
1
n
−
1
c
i
x
i
+
x
n
s(x) = Cx = \sum_{i=1}^{n}c_i x \\ = \sum_{i = 1}^{n-1}c_i x_i + x_n
s(x)=Cx=i=1∑ncix=i=1∑n−1cixi+xn
其中, x x x是状态列向量, C = [ c 1 . . . c n − 1 1 ] C = [c_1\quad . . .\quad c_{n-1} \quad 1] C=[c1...cn−11]。
在滑模控制中,参数 c 1 , c 2 , . . . , c n − 1 c_1,c_2,...,c_{n-1} c1,c2,...,cn−1应该满足多项式 p n − 1 + c n − 1 p n − 2 + . . . + c 2 p + c 1 p^{n-1}+c_{n-1}p^{n-2}+...+c_2 p+c_1 pn−1+cn−1pn−2+...+c2p+c1为 H u r w i t z Hurwitz Hurwitz(赫尔维茨)稳定, p p p为拉普拉斯算子。
例如,当 n = 2 n = 2 n=2时, s ( x ) = c 1 x 1 + x 2 s(x) = c_1x_1+x_2 s(x)=c1x1+x2,为了保证多项式 p + c 1 p+c_1 p+c1为 H u r w i t z Hurwitz Hurwitz稳定 ,需要多项式 p + c 1 = 0 p+c_1 = 0 p+c1=0的特征值实数部分小于零,即 c 1 > 0 c_1 >0 c1>0。
又如,当 n = 3 n=3 n=3时, s ( x ) = c 1 x 1 + c 2 x 2 + x 3 s(x) = c_1x_1+c_2x_2+x_3 s(x)=c1x1+c2x2+x3,为了保证多项式 p 2 + c 2 p + c 1 p^2+c_2p+c_1 p2+c2p+c1为 H u r w i t z Hurwitz Hurwitz稳定,需要多项式 p 2 + c 2 p + c 1 = 0 p^2+c_2p+c_1=0 p2+c2p+c1=0的特征值实数部分为负。
不妨取 p 2 + 2 λ + λ 2 = 0 p^2+2\lambda+\lambda^2=0 p2+2λ+λ2=0,则 ( p + λ ) 2 (p+\lambda)^2 (p+λ)2,取 λ > 0 \lambda > 0 λ>0,可满足多项式 p 2 + c 2 p + c 1 = 0 p^2+c_2p+c_1 = 0 p2+c2p+c1=0的特征值满足实数部分为负,对应得可得到 c 2 = 2 λ , c 1 = λ 2 c_2 = 2\lambda,c_1 = \lambda^2 c2=2λ,c1=λ2。
滑模问题在轨迹跟踪控制上的应用
对于跟踪问题,若设计滑模函数为
s
(
t
)
=
c
e
(
t
)
+
e
˙
(
t
)
s(t) = ce(t) + \dot e(t)
s(t)=ce(t)+e˙(t)
其中,
e
(
t
)
=
x
−
x
d
e(t) = x - x_d
e(t)=x−xd是跟踪误差,
x
x
x为实际位置,
x
d
x_d
xd为期望位置,
e
˙
(
t
)
\dot e(t)
e˙(t)为跟踪误差变化率, 则
c
c
c 必须满足
H
u
r
w
i
t
z
Hurwitz
Hurwitz条件
c
>
0
c>0
c>0。
当 s ( t ) = 0 s(t) = 0 s(t)=0时, c e ( t ) + e ˙ ( t ) = 0 , e ˙ ( t ) = − c e ( t ) ce(t)+\dot e(t) = 0, \dot e(t) = -ce(t) ce(t)+e˙(t)=0,e˙(t)=−ce(t),即 1 e ( t ) e ˙ ( t ) = − c \frac{1}{e(t)}\dot e(t) = -c e(t)1e˙(t)=−c。积分可得 ∫ 0 t 1 e ( t ) e ˙ ( t ) d t = ∫ 0 t − c d t \int_{0}^{t} \frac{1}{e(t)}\dot e(t) dt =\int_{0}^{t}-c dt ∫0te(t)1e˙(t)dt=∫0t−cdt,解得 l n e ( t ) e ( 0 ) = − c t ln \frac{e(t)}{e(0)} = -ct lne(0)e(t)=−ct
收敛结果为 e ( t ) = e ( 0 ) e − c t e(t) = e(0) e^{-ct} e(t)=e(0)e−ct,即当 t → ∞ t \rightarrow \infty t→∞,误差指数收敛于 0,收敛速度取决于 c c c值。
若通过控制律的设计,保证 s ( t ) s(t) s(t)指数收敛于零,则当 t → ∞ t \rightarrow \infty t→∞时,误差变化率也是指数趋近于零。
滑模面设计过程、滑模控制仿真、非线性系统控制器设计仿真传送门:
参考内容:《滑模变结构控制Matlab仿真》
原文地址:https://blog.csdn.net/AII_IIA/article/details/142553510
免责声明:本站文章内容转载自网络资源,如本站内容侵犯了原著者的合法权益,可联系本站删除。更多内容请关注自学内容网(zxcms.com)!