滑模变结构理论一、引言滑模变结构控制本质上是一类特殊的非线性控制,其非线性表现为控制的不连续性,这种控制策略与其它控制的不同之处在于系统的“结构”并不固定,而是可以在动态过程中根据系统当前的状态(如偏差及其各阶导数等)有目的地不断变化,迫使系统按照预定“滑动模态”的状态轨迹运动。

由于滑动模态可以进行设计且与对象参数及扰动无关,这就使得变结构控制具有快速响应、对参数变化及扰动不灵敏、无需系统在线辩识,物理实现简单等优点。

该方法的缺点在于当状态轨迹到达滑模面后,难于严格地沿着滑模面向着平衡点滑动,而是在滑模面两侧来回穿越, 从而产生颤动。

滑模变结构控制出现于20世纪50年代,经历了 50余年的发展,已形成了一个相对独立的研究分支,成为自动控制系统的一种一般的设计方法。

以滑模为基础的变结构控制系统理论经历了 3个发展阶段.第1阶段为以误差及其导数为状态变量研究单输入单输出线性对象的变结构控制; 20世纪60年代末开始了变结构控制理论研究的第2阶段, 研究的对象扩大到多输入多输出系统和非线性系统;进入80年代以来, 随着计算机、大功率电子切换器件、机器人及电机等技术的迅速发展, 变结构控制的理论和应用研究开始进入了一个新的阶段, 所研究的对象已涉及到离散系统、分布参数系统、滞后系统、非线性大系统及非完整力学系统等众多复杂系统, 同时,自适应控制、神经网络、模糊控制及遗传算法等先进方法也被应用于滑模变结构控制系统的设计中。

二、基本原理带有滑动模态的变结构控制叫做滑模变结构控制(滑模控制)。

所谓滑动模态是指系统的状态被限制在某一子流形上运动。

通常情况下,系统的初始状态未必在该子流形上,变结构控制器的作用在于将系统的状态轨迹于有限时间内趋使到并维持在该子流形上,这个过程称为可达性。

系统的状态轨迹在滑动模态上运动并最终趋于原点,这个过程称为滑模运动。

滑模运动的优点在于,系统对不确定参数和匹配干扰完全不敏感。

下图简要地描述了滑模变结构控制系统的运动过程,其中S(t)为构造的切换函数(滑模函数), S(t)=0为滑模面。

图1为了更好的解释滑模变结构控制,考虑如下单输入线性系统的状态调节问题：(t) = A (t) + Bu(t)x x其中, (t)R (t)R x u ∈∈和，分别表示系统的状态和输入。

在线性系统状态调节器设计中,状态反馈控制器设计为:u(t) = Kx(t)其中,状态反馈阵K 可以通过极点配置方法或者线性二次调节器方法设计。

可以看到,上面设计的控制器是固定不变的,但是在滑模变结构控制系统中,控制器结构根据切换函数而变化。

滑模变结构控制器通常设计为如下形式:(t)S(t)>0(t)(t)S(t)<0u u u +-⎧⎪=⎨⎪⎩其中, S(t)为切换函数, S(t)= 0为滑模面。

因此,变结构控制主要体现为滑模面两侧所设计的控制器-(t)(t)u u +≠。

从上式可以看出，当系统状态从区域}{+=(t)(t)0x S ϕ>进入}{=(t)(t)0x S ϕ-<时，系统由(t) = A (t) + Bu (x(t),t)x x +变成(t) = A (t) + Bu (x(t),t)x x -，即滑模变结构控制系统在滑模面S(t)=0附近不连续。

因此,滑模变结构控制的本质是将具有不同结构旳反馈控制系统按照一定的逻辑规则进行切换,并且使得闭环控制系统具备渐近稳定等良好的动态品质。

滑模变结构控制系统的响应由到达(或趋近)阶段、滑动阶段和稳态阶段组成,因此,滑模变结构控制需要满足以下三个条件:(1)滑模面满足可达条件,即系统状态轨迹于有限时间内到达切换面S(t)=0上;(2)滑模面上存在滑动模态;(3)滑动模态具有渐近稳定等良好的动态品质;三、控制的设计滑模控制的设计一般包含以下两步:第一,设计适当的滑模面,使得系统的状态轨迹进入滑动模态后具有渐近稳定等良好的动态特性;第二,设计滑模控制律,使得系统的状态轨迹于有限时间内被趋使到滑模面上并维持在其上运动。

(1)滑模面的设计目前,滑模面主要有线性滑模面、分段滑模面、移动滑模面、积分型滑模面和模糊滑模面等。

滑模面的设计方法有基于标准型(正则型)的设计方法、基于李雅普洛夫的设计方法、基于频率整形的设计方法以及基于LMI 方法等。

面时,不仅要选择合适的设计方法,还要考虑被控系统的结构特点,设计形式恰当的滑模面。

本文工作釆用的主要是线性滑模面和积分型滑模面,下面我们将针对这两种滑模面进行阐述。

1)线性滑模面1)线性滑模面研究表明,线性滑模面是建立滑模变结构控制系统的一种最好结构,因此,我们首先介绍常用的基于标准型的线性滑模面设计方法。

考虑线性不确定系统:(t) = A (t) + Bu(t)(t)x x Df + （1-1） 其中, ()n x t R ∈为系统状态，(t)m u R ∈，输入矩阵n m B R ⨯∈为列满秩，(t)f 为外界干扰。

构造线性切换函数：S(x) = Cx(t) （1-2） 当系统（1-1）中的干扰(t)f 满足匹配条件时，存在矩阵M ，使得D=BM ，对系统(1-1进行相似变换z = Tx 得到标准型为:1111122(t)(t)A Z (t)Z A Z =+22112222(t)(t)A Z (t)B (u(t)Mf(t))Z A Z =+++ （1-3） 其中，11121112212220,,A A TAT TB C CT C C A A B --⎡⎤⎡⎤⎡⎤====⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦，此时滑模函数（1-2）变为：1122(t)()(t)(t)Cz C z C t z S =+= 当系统处于滑动模态阶段时,系统状态轨迹到达滑模面S(t)=0,即 1122(t)(t)0C z C z +=假设CB 非奇异，由于B 2是列满秩的，因此,矩阵C 2非奇异，结合上式可得： 12211(t)(t)z C C z -=将其带入系统（1-3）)的前n-m 维子系统得到滑动模态: 111112211(t)()Z (t)Z A A C C -=- （1-4）利用线性系统理论的传统设计方法,如极点配置设计切换函数中的12C C 和可以保证滑动模态(1-4)的稳定性以及其它动态品质。

近年来,基于标准型设计滑模面的方法得到了进一步的推广和应用。

这种方法应用到线性多变量系统的滑动超平面设计中,同时将二次型范函指标优化方法引入到滑模面的设计中。

值得一提的是Slotine J.-J.E.在1983提出的对于相变量表示的非线性系统给出一个很好的线性滑模面形式： 11()n d s x dtλ-=+ 其中，λ>0。

对于这种线性滑模面形式下面将详细叙述其原理：考查单输入动态系统：(x)b(x)u n x f =+u 是控制输入，(n 1)[,,...,]T x x x x -=是状态向量，f(x)不是精确已知，但是不精确性范围的上界是x 的一个已知连续函数，类似的b(x)不精确已知，但其符号已知并且其范围受x 的连续函数界定。

记d x x x =-令：(1)[,,...,]n T d x x x x x x -==-用标量方程S （x,t ）=0定义状态空间n R 中的时变曲线S （t ）1(x;t)()n d s x dtλ-=+ （1-5） 当给定初始状态：(0)x(0)d x =时，跟踪d x x ≡的问题就等价于当t>0时，轨线必须停留在曲面S （t ）上，事实上0s ≡代表了一个线性微分方程，怎唯一的解是0x ≡。

因此跟踪n 维向量d x 的问题可简化成使标量S 恒为零的问题。

更确切的说，跟踪n 维向量d x 的问题能够有效的被用S 表示的一阶镇定问题取代.事实上（1-5）包含(1)n x -，则对S 微分一次就可使得u 出现。

进一步的，S 的界可直接转换成跟踪误差向量x 的界限，因此标量S 是跟踪性能的真实度量，特别的，假定(0)0x =，则有：0,(t)0,(t)(2)0,...,1i i t S t x i n φλε∀≥≤⇒∀≥≤=- （1-6） 其中1/n εφλ-=，根据（1-5）定义，跟踪误差x 可由S 通过一系列的一阶低通滤波器获得，记y 1为第一个滤波器的输出，可得：(t T)10(t)(T)dT ty e s λ--=⎰ 从s φ≤可得：(t T)10y (t)(/)(1e )/t t e dT λλφφλφλ---≤=-≤⎰对于第二个滤波器运用相同的推理，并继续，直到1y n x -=。

同理可得1/n x φλε-≤=类似的，()i x 可以看成是S 通过图2中的步骤得到，由前面的结论得 11/n i z φλ--≤其中z 1是（n-i-1）阶滤波器的输出。

此外，注意到：1p p p p p λλλλλλ+-==-+++ 可见图2中的步骤意味着：(i)1()(1)(2)i i n i x φλλελλ--≤+= 这就是界（1-6），最后当0x ≠时，界（1-6）可以渐进的得到，即在一小段恒定时间λ（n-1）/内得到，因此这就由一个一阶镇定问题有效地取代了n 阶的跟踪问题，并且用（1-6）式量化了性能度量的相应变换。

图2可以选择(x)b(x)u n x f =+中的控制u ，使得曲面S （t ）之外满足：21d 2s s dtη≤- （1-7） 以得到使得S 恒为零的简化问题，其中η是正常数，（1-7）表达的是以2s 为度量到曲面的平方距离沿所有系统轨线减小，因此使得轨线趋近于曲面S(t)，如图3所示，轨线一旦进入曲面就将一直停留在曲面上，换句话说，系统轨线满足（1-7）即滑动条件，使曲面成为一个不变集。

此外，（1-7）表明在有一些干扰和系统不确定时，仍保持曲面是个不变集，从图3看出，不在曲面上的轨线仍能指向曲面运动。

满足（1-7）的曲面S （t ）称为滑动曲面，并且系统的性态一旦在曲面就被称为滑动形态或者滑动模。

图3不变集S(t)另一个有趣的方面是，一旦系统轨线在曲面上，系统的轨线由不变集自身的方程定义，即： 1()=0n d x dtλ-+ 换句话说，S(t)既是一个曲面也是一个动态。

最后，满足（1-7）的条件保证了即使初始条件不严格成立（即x(t=0)偏离了(t 0)d x =）时，系统轨线仍能在小于(t 0)/s η=的有限时间内到达曲面S(t)。

事实上，如果假定S （t=0）>0,记t reach 为到达曲面S=0所需要的时间，对于（1-7）从t=0到t=t reach 积分可得到0-S （t=0）=S(t=t reach )-S(t=0) ≤ (0)reach t η--这表明由S （t=0）<0可以得到类似结果，因此s(t 0)/reach t η≤= 进而，1(x;t)()n d s x dtλ-=+一旦在曲面上，跟踪误差以时间常数（n-1）/λ指数趋于零（因设计中共（n-1）个时间常数等于1/λ的滤波器）。