CN104614995B - 一种二阶***有限时间滑模控制器的通用设计方法 - Google Patents
一种二阶***有限时间滑模控制器的通用设计方法 Download PDFInfo
- Publication number
- CN104614995B CN104614995B CN201510084429.8A CN201510084429A CN104614995B CN 104614995 B CN104614995 B CN 104614995B CN 201510084429 A CN201510084429 A CN 201510084429A CN 104614995 B CN104614995 B CN 104614995B
- Authority
- CN
- China
- Prior art keywords
- sliding mode
- finite time
- order system
- function
- finite
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000013461 design Methods 0.000 title claims abstract description 54
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000013507 mapping Methods 0.000 claims description 9
- 241000135164 Timea Species 0.000 claims 1
- 230000004044 response Effects 0.000 abstract description 3
- 230000000694 effects Effects 0.000 abstract 1
- 238000009415 formwork Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 17
- 230000009471 action Effects 0.000 description 15
- 238000004088 simulation Methods 0.000 description 7
- 238000004458 analytical method Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 4
- 238000012886 linear function Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 241000764238 Isis Species 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012887 quadratic function Methods 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Landscapes
- Feedback Control In General (AREA)
Abstract
本发明公开了一种二阶***有限时间滑模控制器的通用设计方法,包括:设计二阶***设计有限时间滑模变量形式;确定函数和ψ(x2)的设计条件;如果函数和ψ(x2)满足步骤二中的设计条件,则根据设计二阶***有限时间滑模控制器,从而实现二阶***的有限时间控制。考虑到已有的滑模控制中,终端滑模面和相应的有限时间控制器形式过于具体,无通用设计准则,本发明提出有限时间滑模变量的通用设计准则,获得通用形式的有限时间滑模面和相应的有限时间滑模控制器。在线性滑模控制器作用下,***状态于无限时间内渐近收敛到平衡点。本发明提出的二阶***有限时间滑模控制器的通用设计方法,相比于线性滑模控制器,有效提高了控制***的响应速度。
Description
技术领域
本发明属于滑模控制技术领域,更具体地,提出一种二阶***有限时间滑模面的通用设计方法,基于滑模控制的原理,设计出相应的有限时间滑模控制器,实现对二阶***有限时间控制的方法。
背景技术
滑模控制方法,作为一种***综合的方法,具有相当大的灵活度。因为在滑模面上,***的动态由滑模面来决定,而滑模面是预先设计的。滑模面设计的好坏,直接关系到***的动态性能。只要能够使滑模面是可到达的,并保证滑动模态的稳定性,滑模面的形式可以不拘一格。因而,对滑动模态的设计研究也成为滑模控制中的一个热点问题,从中也衍生出很多新的滑模控制方法。
对不同的滑模面,利用滑模控制器存在的充要条件,可以获得***的滑模控制器。因此,滑模控制的品质与滑模面的形式紧密相关。以往文献中设计的滑模变量的形式有:
(1)线性函数:即滑模面设计为误差和误差导数的线性函数。线性滑模面可以满足线性***的控制性能要求,其稳定性分析简单。但是,在线性滑模面,***最好的收敛是渐近收敛,当控制对象具有复杂非线性时,线性滑模面往往不能得到较好的控制性能甚至无法使***稳定。
(2)特殊二次型函数:l(x)=cxTx,其中x是***的状态,l(x)是滑模变量,c为常量。由于采用二次型滑模面的控制***,参数设计复杂,稳定性分析困难,现在已经少有研究。
(3)分段函数:文献(Man Z.H.,Palaniswami M.Decentralised three-segmentnonlinear sliding mode control for robotic manipulators[J].IEEE InternationalWorkshop on Emerging Technologies and Factory Automation,1992,607-612.)中提出了一种分段滑模面,它将滑动模态运动分为加速段、恒速段、减速段,分别对每一阶段设计滑模面。这种设计方法,滑动模态整体是不连续的,每一个分段函数是连续的,具有更加迅速的暂态响应。但是,由于滑模面是分段不连续的,稳定性分析会变得困难,控制器在滑模面切换时也有可能出现奇异的情况。
(4)终端滑模面:文献(Man Z.,Yu X.Terminal sliding mode control of MIMOsystems[J].IEEE Transactions on Systems I:Fundamental Theory andApplications,1997,44(11),1065-1070.)中提出了终端滑模面。通过将滑模面设计为非线性函数,***的状态可以在有限时间内收敛到平衡点。这种滑模面称为终端滑模面。终端滑模面方法,为***带来了比渐近收敛更好的收敛性能,稳定性分析采用Lyapunov方法即可,控制律设计简单,并且在机器人、电机等复杂非线性***中得到了成功的应用,因而逐渐形成了滑模控制的一个新的分支-终端滑模控制。为了进一步提高终端滑模的控制品质,在终端滑模面的基础上,又提出了快速终端滑模面,加快了到达阶段的趋近速度。之后,为了消除终端滑模控制的奇异性问题,先后提出了非奇异终端滑模面和快速非奇异终端滑模面。
传统滑模面通常设计为线性函数。线性滑模面可以满足***的控制性能要求。但是,在线性滑模面上,***最好的收敛是渐近收敛,即当控制对象具有较高的控制精度要求时,渐近收敛往往不能满足***控制精度的要求。终端滑模面,通过将滑模面设计为非线性函数,***的状态可以在有限时间内收敛到平衡点,相应的控制律的设计方法和线性滑模面相同。目前,有限时间滑模面的设计方法过于具体,目前尚无通用的有限时间滑模面的设计准则。
发明内容
针对现有技术中存在的上述问题,本发明提出一种二阶***有限时间滑模控制器的设计方法,通过设计通用的有限时间滑模面,解决无通用有限时间滑模面设计准则的问题,并基于有限时间滑模面的通用设计准则,给出有限时间滑模控制器的通用设计形式,符合非线性***滑模控制技术的应用需求和发展趋势。
为了解决上述技术问题,本发明提出的一种二阶***有限时间滑模控制器的通用设计方法,包括如下步骤:
步骤一、根据二阶***设计有限时间滑模变量形式:
所述二阶***的形式为
其中,x=[x1,x2]T是***状态,和分别表示状态x1和x2的导数,f(x)、b(x)是标量函数,且分别光滑可导,且|b(x)|≠0,u为二阶***的控制输入,u是标量函数,d(t)为有界干扰,d(t)为标量函数且满足|d(t)|<Ld,Ld为正常数;
有限时间滑模变量的形式为:
其中,s(x)是有限时间滑模变量,是x1的函数,ψ(x2)是x2的函数;
步骤二、函数和ψ(x2)的设计条件,包括:
条件(1):当x1=0时,当x2=0时,ψ(x2)=0;当x1≠0时,的导数满足当x2≠0时,的导数满足ψ′(x2)>0;则函数和ψ(x2)是过零点严格单调递增的函数,和ψ(x2)存在逆函数;
条件(2):当有限时间滑模变量s(x)等于0时,x1与x2一一映射,二阶***的动态在x1→0时是齐次的,记ψ-1(·)表示ψ(·)的逆函数;
条件(3):根据计算使对任意x1、x2有意义,即不会无穷大;
步骤三、如果函数和ψ(x2)满足步骤二中的设计条件,令就得到通用有限时间滑模面的表达式,执行步骤四;否则,返回步骤一重新设计有限时间滑模变量
步骤四、根据设计二阶***有限时间滑模控制器的表达式为:
其中,η>Ld是二阶***有限时间滑模控制器增益参数,sgn(s(x))是s(x)的符号函数,
步骤五:采用步骤四获得的二阶***有限时间滑模控制器,实现步骤一中所述二阶***的有限时间控制。
与现有技术相比,本发明的有益效果是:
本发明提出一种二阶***有限时间滑模控制器的通用设计方法,其中,步骤一和步骤二完成二阶***通用的有限时间滑模变量的设计,主要是通过设计有限时间滑模变量中两个独立的函数和ψ(x2),给出其应该满足的充要条件,获得通用有限时间滑模变量令从而在步骤三中就得到通用有限时间滑模面的表达式;继续通过步骤四,在二阶***通用的有限时间滑模变量设计基础上,完成二阶***通用有限时间滑模控制器的设计,该控制器能够使二阶***有限时间内收敛到原点。
通过使用本发明提供的二阶***有限时间滑模面和有限时间滑模控制器的设计方法,能够得到多种有限时间滑模面和对应的有限时间滑模控制器的形式,以往研究的终端滑模面和终端滑模控制器仅是本发明中通用滑模面和通用有限时间滑模控制器的一个具体表现形式。
附图说明
图1是二阶***有限时间滑模控制器的通用设计方法的流程图;
图2是二阶***从任意初始位置x0=(x1(0),x2(0))出发,在到达滑模面进而到达原点的状态轨迹示意图;
图3是x1∈[-5,5]时终端滑模面和快速终端滑模面的仿真图;
图4是二阶***在uFNTSM和uNTSM作用下,从初始位置(2,1)到达(0,0)的相轨迹图;
图5是x1∈[-1.2,1.2]时双曲线型非奇异终端滑模面、快速终端滑模面和终端滑模面的仿真图;
图6是二阶***在uNTSM、uFNTSM、uHNTSM作用下,从初始位置(1,3)到达(0,0)的相轨迹图;
图7是二阶***状态在uNTSM、uFNTSM、uHNTSM作用下随时间收敛曲线;
图8是有限时间控制器uNTSM、uFNTSM、uHNTSM的仿真图;
图9是有限时间滑模变量sNTSM,sFNTSM,sHNTSM在uNTSM、uFNTSM、uHNTSM作用下从初值随时间收敛到0的曲线。
具体实施方式
本发明一种二阶***有限时间滑模控制器的通用设计方法的思路是:基于终端滑模面的齐次性性质,提出有限时间滑模面的设计准则,获得通用形式的有限时间滑模面;依据李雅普诺夫稳定性理论,设计通用的有限时间滑模控制器。有限时间滑模控制器能够使二阶***有限时间到达平衡点,相比于线性滑模控制器,提高了控制***的响应速度。
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
如图1所示,本发明提供了二阶***有限时间滑模控制器的通用设计方法,具体步骤如下:
步骤一、考虑二阶***的形式为
其中x=[x1,x2]T是***状态,和分别表示状态x1和x2的导数,f(x)、b(x)是标量函数,且光滑可导,且|b(x)|≠0,u为二阶***的控制输入,u是标量函数,d(t)为有界干扰,即|d(t)|<Ld,Ld为常数。设计有限时间滑模变量的通用形式为其中s(x)是有限时间滑模变量,仅是x1的函数,ψ(x2)仅是x2的函数;
步骤二、设计函数和ψ(x2)的具体形式,有以下条件:
条件(1):使得当x1=0,x2=0时,有和ψ(x2)=0;对任何x1≠0,都有对任何x2≠0,都有ψ′(x2)>0;函数和ψ(x2)是过零点严格单调递增的函数,和ψ(x2)存在逆函数;
条件(2):使得在有限时间滑模变量s(x)等于0时,,x1与x2是一一映射,二阶***的动态在x1→0时是齐次的;
条件(3):根据计算使得对任意x1、x2有意义;
步骤三、判断:如果满足上述条件(1)、条件(2)和条件(3),令得到通用有限时间滑模面的表达式,则进入步骤四;否则,返回步骤一重新设计有限时间滑模变量
步骤四:根据满足步骤一及步骤二中的设计条件的设计二阶***有限时间滑模控制器的通用形式为其中η>Ld是控制器增益参数,sgn(s(x))是s(x)的符号函数,
步骤五:采用步骤四获得的二阶***有限时间滑模控制器,实现步骤一中所述二阶***的有限时间控制。
具体实施过程中,各步骤的详细内容如下。
步骤一中是本发明通用有限时间滑模面和有限时间控制器的适用对象,具体形式如下:
其中,x=[x1,x2]T是***状态,和分别表示状态x1和x2的导数,f(x)、b(x)是标量函数,且分别光滑可导,且|b(x)|≠0,u为二阶***的控制输入,u是标量函数,d(t)为有界干扰,d(t)为标量函数且满足|d(t)|<Ld,Ld为正常数;有限时间滑模变量函数s(x)的通用形式设计为其中s(x)是有限时间滑模变量函数,仅是x1的函数,表示的导数,ψ(x2)仅是x2的函数,ψ′(x2)表示ψ(x2)的导数,
步骤二的条件(1)中,当x1=0,x2=0时,有和ψ(x2)=0;对任何x1≠0,都有对任何x2≠0,都有ψ′(x2)>0;表明函数ψ(x2)是过零点严格单调递增的函数,因此和ψ(x2)存在逆函数。
步骤二的条件(2)中,使得在有限时间滑模面s(x)=0上,x1与x2是一一映射。这一点可以用检验函数的界和ψ(x2)的极限来保证,即计算的上确界和下确界计算ψ(x2)在x2→+∞和x2→-∞的极限和通过界和极限的关系来保证一一映射成立。特别的,当和 的上确界和下确界应该满足和
步骤二的条件(2)中,使得在有限时间滑模面s(x)=0上,二阶***的动态在x1→0时是齐次的。具体包括:
一阶***的齐次性定义是:如果标量函数f(x1)是齐次的,则***是齐次的,x1是***状态。
标量函数f(x1)的齐次性定义是:连续标量函数f(x1)是齐次的,如果它满足其中ρ>0为辅助变量,r>0,r1>0为选取的幂指数,m是齐次性系数。
由于在有限时间滑模面s(x)=0上,***的动态被描述为在x1→0时,检验微分方程的齐次性,实际上是检验函数的齐次性。令g(x1)=-γsigα(x1),γ为正的常系数。依据:g(ρx1)=-γραsigα(x1)=ρr+mh(x)=ρ1+(α-1)h(x),ρ>0,r=r1=1,m=α-1(2)
根据标量函数齐次性定义,故函数g(x1)=-γsigα(x1)是齐次性函数。因此,在x1→0时,若等价于-γsigα(x1),则***也是齐次的,满足了使得在有限时间滑模面s(x)=0上,二阶***的动态在x1→0时是齐次的条件。
由于在通用的二阶***有限时间控制器设计中存在除式,步骤二中的条件(3)中,需要计算对任意x1、x2有意义,即不会无穷大。
步骤二给出了有限时间滑模变量函数s(x)的通用设计准则。满足步骤二中的有限时间滑模变量函数s(x),可以作为有限时间滑模面s(x)=0。这是因为s(x)=0意味着故
因为函数ψ(·)是一个单调增函数,且为奇函数,因而必存在逆函数ψ-1(·)。根据公式(3),可以得到
记由于函数ψ-1(·)均为严格单调增函数,故函数也为严格单调增函数,且在x=0处满足即
也就是说,当x1≠0时,x1h(x1)<0总是成立,当x1=0时,x1h(x1)=0。
对于有限时间滑模面s(x)=0上任意一点它有限时间内收敛到平衡点(0,0),收敛时间Tx的表达式为
由于故是奇异的。公式(6)是一个广义积分,利用柯西收敛准则说明公式(6)的收敛性。公式(6)是有限的,说明有限时间滑模面s(x)=0上任意一点在有限时间Tx内收敛到平衡点(0,0)。因此考虑函数sigα(x1),根据步骤二给出的:使得在有限时间滑模面s(x)=0上,二阶***的动态在x1→0时是齐次的条件,要求在x1→0时,与sigα(x1)等价,表示为:
由于求导公式成立,故
因此满足步骤二中条件(1),(2),(3)的有限时间滑模面s(x)=0上任意一点在有限时间Tx内收敛到平衡点(0,0),s(x)是有限时间滑模变量,s(x)=0是有限时间滑模面。
步骤三进行判断,如果有限时间滑模变量满足步骤二中的各条件,令就得到通用有限时间滑模面的表达式,可以进入步骤四;否则,从步骤一开始重新设计有限时间滑模变量
步骤四中给出了通用有限时间滑模变量对应的二阶***有限时间滑模控制器,表达式为其中η>Ld>0是控制器增益参数,sgn(s(x))是s(x)的符号函数。由于在步骤二中检验了分式有意义,故二阶***通用的有限时间滑模控制器是非奇异的,即有界的。
步骤五是将由步骤四获得的二阶***有限时间滑模控制器代入步骤一中,实现二阶***的有限时间控制。加入二阶***通用的有限时间滑模控制器后,闭环***的方程为:
其中η>Ld>0为控制器增益。令两个正的系数η′和η″为η′>η-Ld,η″>η+Ld。明显η′,η″>0。由此,可以得到有限时间滑模变量s(x)的导数为
如果s(x)>0,
如果s(x)<0,
总结起来,可以有如下不等式成立:
因此,滑模控制理论提出的滑模面存在的充要条件,可以获得如下:
也就是说,如果ψ′(x2)≠0,公式(1)所示的二阶***可以有限时间内到达滑模面。
如果ψ′(x2)=0,定义x2=0附近一个无限小的ε区间Pε,Pε={x:-ε<x2<+ε},ε>0。既然在开区间Pε外,对于任意x,都有ψ′(x2)>0成立。R2表示直角坐标系中横纵坐标均为实数组成的平面,故x∈R2\{Pε}表示在平面R2除去开区间Pε的平面。令ψ′(x2)在平面x∈R2\{Pε}上的下限为ε*,即有下式成立
因此,根据公式(12),可以得到如下的不等式
根据滑模控制理论提出的滑模面存在的充要条件,由公式(14)可以获得,在除去开区间Pε的平面x∈R2\{Pε}上,公式(1)所示的二阶***在有限时间滑模控制器作用下,***状态可以有限时间内到达滑模面。
如果ψ′(x2)=0,即x2=0,由于步骤二中设计对任意x1、x2有意义,因此成立。考虑公式(1)所示***的在有限时间滑模控制器作用下的闭环***,对于s(x)>0,有
故且如果s(x)<0,有
故且由于ε是一个无限小的正数,当x2=0时因此x2能够跨越包含x2=0的无限小区间Pε,到达平面x∈R2\{Pε}上。因此,无论x2=0或者x2≠0,公式(1)所示的二阶***在有限时间滑模控制器作用下,***状态可以有限时间内到达滑模面。
当***状态有限时间到达滑模面后,***的运行动态由滑模面的形式决定。本发明中设计的通用有限时间滑模面的形式为故公式(1)所示的二阶***在滑模面上的动态方程为:
选择李雅普诺夫函数可以得到它的导数为:
因此,根据李雅普诺夫稳定性方法,***状态x1渐近收敛到原点。
根据步骤二给出的条件,当x1=0,x2=0时,有和ψ(x2)=0;对任何x1≠0,都有对任何x2≠0,都有ψ′(x2)>0。因此,对***状态运行于滑模面上,且尚未到达原点时,满足:
依据公式(9),二阶***在通用的有限时间滑模控制器作用下,***状态x2的动态为在滑模面上s(x)=0,由于***动态对于匹配性干扰具有不变性,故***状态x2的动态为当***状态运行于滑模面上,且尚未到达原点时,x2≠0总成立。故选取李雅普诺夫函数V2>0,其导数为:
结合公式(19),显然因此,根据李雅普诺夫稳定性方法,***状态x2渐近收敛到原点。
x0=(x1(0),x2(0))表示二阶***的状态任意初始位置。***状态从x0=(x1(0),x2(0))出发,经过有限时间t1后到达滑模面是表示二阶***状态在有限时间滑模面上的初始位置。设***状态从有限时间滑模面上的初始位置出发,到达原点所用的时间为Tx,则对公式(17)积分可以得到Tx的表达式为:
依据公式(6),Tx<∞。因此,结合***状态x1,x2渐近收敛到原点的结论,可以得到,在本发明所设计的通用有限时间滑模面和相应的有限时间滑模控制器下,如公式(1)所示的二阶***能够实现有限时间控制,***状态x1,x2有限时间到达滑模面并有限时间收敛到原点。图2是二阶***从任意初始位置x0=(x1(0),x2(0))出发,在到达滑模面进而到达原点的状态轨迹示意图。
为使本领域技术人员更好地理解本发明,下面结合具体实施例,对本发明的二阶***有限时间滑模控制器的通用设计方法进行详细说明。
(1)快速非奇异终端滑模变量(sFNTSM)和相应有限时间控制器(uFNTSM)
快速非奇异终端滑模变量sFNTSM的形式设计如下:
其中k,k1,α,α1是快速非奇异终端滑模变量的参数,满足k>0,k1>0,α>0,α1>0,令sFNTSM=0就可以得到快速非奇异终端滑模面的表达式。sFNTSM是按照步骤一至步骤三设计获得,这里检验公式(22)所示的快速非奇异终端滑模变量sFNTSM是否符合步骤二中所列条件。
①:步骤二中的设计条件(1),函数和满足x1=0,x2=0时,和ψ(x2)=0的条件;对任何x1≠0,都有对任何x2≠0,都有ψ′(x2)>0成立。
②:步骤二中的设计条件(2),要求在快速非奇异终端滑模面sFNTSM=0上,x1与x2是一一映射。这一点可以用检验函数的界和ψ(x2)的极限来保证。计算的上确界和下确界
计算ψ(x2)在x2→+∞和x2→-∞的极限
因此的界和ψ(x2)的极限存在关系:
因此,在快速非奇异终端滑模面sFNTSM=0上,x1与x2是一一映射。
③:步骤二中的设计条件(2),由于在快速非奇异终端滑模面sFNTSM=0上,***的动态被描述为在x1→0时,检验微分方程的齐次性,考虑到α和α1的取值,实际上是检验方程的齐次性。方程的齐次性取决于函数-(kx1)α的齐次性。-(kx1)α可以表示为-γsigα(x1)的形式,γ=kα,
由于函数g(x1)=-γsigα(x1)是齐次性函数,即-(kx1)α是齐次性函数。因此,在x1→0时,***是齐次的,满足了二阶***的动态在x1→0时是齐次的条件。
④:步骤二中的设计条件(3),计算
考虑α和α1的取值,可以得到对任意x1、x2有意义。因此,快速非奇异终端滑模变量是符合条件的。
根据步骤四,可以得到快速非奇异终端滑模面sFNTSM相应的有限时间控制器uFNTSM,表达为:
该有限时间控制器能够让公式(1)所示的二阶***有限时间到达原点。
仿真中,将设计的快速非奇异终端滑模面和已有的终端滑模面进行比较。采用的参数为k=1,k1=1。
图3是在x1∈[-5,5]非奇异终端滑模面和快速非奇异终端滑模面的仿真图,其中虚线是非奇异终端滑模面sNTSM=0,实线是快速非奇异终端滑模面sFNTSM=0。图4是二阶***在有限时间控制器作用下,从初始位置(2,1)到达(0,0)的相轨迹图,其中虚线是在非奇异终端滑模变量相应的有限时间控制器作用下,二阶***从初始位置(2,1)到达(0,0)的相轨迹图,实线是在快速非奇异终端滑模变量相应的有限时间控制器作用下,二阶***从初始位置(2,1)到达(0,0)的相轨迹图。
(2)双曲线型非奇异终端滑模变量(sHNTSM)和相应有限时间控制器(uHNTSM)
双曲线型非奇异终端滑模变量sHNTSM设计为
其中k,k1,α,α2是快速非奇异终端滑模变量的参数,满足k>0,k2>0, 令sHNTSM=0就可以得到双曲线型非奇异终端滑模面的表达式。sHNTSM是按照步骤一至步骤三设计获得。对sHNTSM是否符合设计条件进行检验。
①:步骤二中条件(1),在x1=0时,在x2=0,且对于任意x1≠0,x2≠0,都有ψ′(x2)>0成立。
②:步骤二中条件(2),计算的上确界和下确界
计算在x2→+∞和x2→-∞的极限
因此的界和ψ(x2)的极限存在关系:
因此,在双曲线型非奇异终端滑模面sHNTSM=0上,x1与x2是一一映射。
③:步骤二中条件(2),由于在双曲线型非奇异终端滑模面sHNTSM=0上,***的动态被描述为在x1→0时,微分方程与微分方程等价。由于函数g(x1)=-γsigβ(x1)是齐次性函数,故在x1→0时,***是齐次的,满足步骤二中的二阶***的动态在x1→0时是齐次的条件。
④:步骤二中条件(3),计算
考虑可以得到对任意x1、x2有意义。因此,双曲线型非奇异终端滑模变量是符合条件的。
步骤四给出了由有限时间滑模变量获得相应的有限时间器是方法。依据双曲线型非奇异终端滑模变量的形式,可以得到相应的有限时间控制器uHNTSM,表达如下:
在uHNTSM作用下,公式(1)所示的二阶***能够有限时间到达原点。
将双曲线型非奇异终端滑模面快速非奇异终端滑模面和终端滑模面一起仿真比较。仿真中采用的参数为k=1,k1=1。
图5是在x1∈[-1.2,1.2]时,双曲线型非奇异终端滑模面、快速非奇异终端滑模面和非奇异终端滑模面的仿真图,其中虚线是非奇异终端滑模面sNTSM=0,点虚线是快速非奇异终端滑模面sFNTSM=0,实线是双曲线型非奇异终端滑模面sHNTSM=0。图6是二阶***在有限时间控制器uNTSM、uFNTSM、uHNTSM作用下,从初始位置(1,3)到达(0,0)的相轨迹图,其中虚线是二阶***在有限时间控制器uNTSM作用下,从初始位置(1,3)到达(0,0)的相轨迹图,点虚线表示二阶***在有限时间控制器uFNTSM作用下,从初始位置(1,3)到达(0,0)的相轨迹图,实线是在uHNTSM作用下,二阶***从初始位置(1,3)到达(0,0)的相轨迹图。图7给出了二阶***状态x1,x2在有限时间控制器uNTSM、uFNTSM、uHNTSM作用下随时间收敛曲线,显示了从初始位置(1,3)出发,在三种有限时间控制器作用下,二阶***的状态都是有限时间收敛到原点的。图8给出了有限时间控制器uNTSM、uFNTSM、uHNTSM的值。图9描述了三种滑模变量从初始值sNTSM(x0)=7.2403,sFNTSM(x0)=8.2403,sHNTSM(x0)=4.5305在有限时间控制器uNTSM、uFNTSM、uHNTSM作用下随时间收敛到0的曲线。
尽管上面结合附图对本发明进行了描述,但是本发明并不局限于上述的具体实施方式,上述的具体实施方式仅仅是示意性的,而不是限制性的,本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨的情况下,还可以做出很多变形,这些均属于本发明的保护之内。
Claims (3)
1.一种二阶***有限时间滑模控制器的通用设计方法,其特征在于,包括如下步骤:
步骤一、根据二阶***设计有限时间滑模变量形式:
所述二阶***的形式为
其中,x=[x1,x2]T是***状态,和分别表示状态x1和x2的导数,f(x)、b(x)是标量函数,且分别光滑可导,且|b(x)|≠0,u为二阶***的控制输入,u是标量函数,d(t)为有界干扰,d(t)为标量函数且满足|d(t)|<Ld,Ld为正常数;
有限时间滑模变量的形式为:
其中,s(x)是有限时间滑模变量,是x1的函数,ψ(x2)是x2的函数;
步骤二、函数和ψ(x2)的设计条件,包括:
条件(1):当x1=0时,当x2=0时,ψ(x2)=0;当x1≠0时,的导数满足当x2≠0时,ψ(x2)的导数ψ′(x2)满足ψ′(x2)>0;则函数和ψ(x2)是过零点严格单调递增的函数,和ψ(x2)存在逆函数;
条件(2):当有限时间滑模变量s(x)等于0时,x1与x2一一映射,二阶***的动态在x1→0时是齐次的,记ψ-1(·)表示ψ(·)的逆函数;
条件(3):根据计算使对任意x1、x2有意义,即不会无穷大;
步骤三、如果函数和ψ(x2)满足步骤二中的设计条件,令得到通用有限时间滑模面的表达式,执行步骤四;否则,返回步骤一重新设计有限时间滑模变量
步骤四、根据设计二阶***有限时间滑模控制器的表达式为:
其中,η>Ld>0是二阶***有限时间滑模控制器增益参数,sgn(s(x))是s(x)的符号函数,
步骤五:采用步骤四获得的二阶***有限时间滑模控制器,实现步骤一中所述二阶***的有限时间控制。
2.根据权利要求1所述二阶***有限时间滑模控制器的通用设计方法,其特征在于,步骤二的条件(2)中,通过检验函数的界和ψ(x2)的极限来保证当有限时间滑模变量s(x)等于0时,x1与x2一一映射,即:
其中,表示x2→+∞时ψ(x2)的极限,表示x2→-∞时ψ(x2)的极限,表示x1∈R时的上确界,R表示实数域,表示x1∈R时的下确界;
如果ψ(x2)的极限则的界要满足:
3.根据权利要求1所述的二阶***有限时间滑模控制器的通用设计方法,其特征在于,步骤二的条件(2)中,所述二阶***的动态在x1→0时是齐次的,具体内容包括:
结合获得二阶***在有限时间滑模变量s(x)等于0时的运动方程
在x1→0时,检验微分方程的齐次性,也就是检验函数的齐次性;
在x1→0时,进一步检验函数是否等价于-γsigα(x1),其中,γ为正的常系数,sigα(x1)=|x1|αsgn(x1),|x1|表示x1的绝对值,幂指数α满足
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510084429.8A CN104614995B (zh) | 2015-02-16 | 2015-02-16 | 一种二阶***有限时间滑模控制器的通用设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510084429.8A CN104614995B (zh) | 2015-02-16 | 2015-02-16 | 一种二阶***有限时间滑模控制器的通用设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104614995A CN104614995A (zh) | 2015-05-13 |
CN104614995B true CN104614995B (zh) | 2017-03-29 |
Family
ID=53149497
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510084429.8A Active CN104614995B (zh) | 2015-02-16 | 2015-02-16 | 一种二阶***有限时间滑模控制器的通用设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104614995B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105179164B (zh) * | 2015-06-25 | 2018-11-09 | 江苏科技大学 | 基于t-s模糊模型的风能转换***滑模控制方法及装置 |
CN104950901B (zh) * | 2015-07-03 | 2017-07-28 | 天津大学 | 无人直升机姿态误差有限时间收敛非线性鲁棒控制方法 |
CN106444372B (zh) * | 2016-08-25 | 2019-05-31 | 浙江工业大学 | 用于电机伺服***的滑模重复控制器 |
CN106354012B (zh) * | 2016-11-02 | 2019-03-01 | 湖南工业大学 | 基于Terminal滑模的机车牵引力总量控制方法 |
CN107846019B (zh) * | 2017-11-16 | 2020-09-29 | 河海大学常州校区 | 基于分数阶高阶滑模模糊控制的有源电力滤波器控制方法 |
CN108089589A (zh) * | 2017-11-24 | 2018-05-29 | 江苏科技大学 | 一种水下机器人姿态控制方法 |
CN109213184B (zh) * | 2018-11-06 | 2021-06-08 | 哈尔滨工业大学(深圳) | 挠性航天器的有限时间多模态滑模姿态控制算法 |
CN111752157B (zh) * | 2020-07-17 | 2021-07-06 | 哈尔滨工业大学 | 一种有限时间收敛的二阶滑模控制方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011008360A (ja) * | 2009-06-23 | 2011-01-13 | Sumitomo Heavy Ind Ltd | スライディングモード制御系の設計方法及びその設計支援装置 |
JP4638827B2 (ja) * | 2006-02-28 | 2011-02-23 | ジヤトコ株式会社 | スライディングモード制御を用いた制御装置及び制御方法 |
CN102880060A (zh) * | 2012-10-25 | 2013-01-16 | 北京理工大学 | 再入飞行器自适应指数时变滑模姿态控制方法 |
CN104267605A (zh) * | 2014-10-25 | 2015-01-07 | 哈尔滨工业大学 | 适用于相对阶为1控制***的平滑非奇异终端滑模控制方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5723582B2 (ja) * | 2010-11-30 | 2015-05-27 | インターナショナル・ビジネス・マシーンズ・コーポレーションInternational Business Machines Corporation | スライディング・モード制御方式において、制御器を構成する方法、プログラム及びシステム |
-
2015
- 2015-02-16 CN CN201510084429.8A patent/CN104614995B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4638827B2 (ja) * | 2006-02-28 | 2011-02-23 | ジヤトコ株式会社 | スライディングモード制御を用いた制御装置及び制御方法 |
JP2011008360A (ja) * | 2009-06-23 | 2011-01-13 | Sumitomo Heavy Ind Ltd | スライディングモード制御系の設計方法及びその設計支援装置 |
CN102880060A (zh) * | 2012-10-25 | 2013-01-16 | 北京理工大学 | 再入飞行器自适应指数时变滑模姿态控制方法 |
CN104267605A (zh) * | 2014-10-25 | 2015-01-07 | 哈尔滨工业大学 | 适用于相对阶为1控制***的平滑非奇异终端滑模控制方法 |
Non-Patent Citations (3)
Title |
---|
有限时间收敛的滑模自适应控制器设计;赵占山 等;《山东大学学报(工学版)》;20120831;第42卷(第4期);全文 * |
超声速导弹的光滑二阶滑模控制器设计;王宏伟 等;《指挥控制与仿真》;20141031;第36卷(第5期);全文 * |
非奇异终端滑模控制***相轨迹和暂态分析;穆朝絮 等;《自动化学报》;20130630;第39卷(第6期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN104614995A (zh) | 2015-05-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104614995B (zh) | 一种二阶***有限时间滑模控制器的通用设计方法 | |
CN108303885B (zh) | 一种基于干扰观测器的电机位置伺服***自适应控制方法 | |
CN107121932B (zh) | 电机伺服***误差符号积分鲁棒自适应控制方法 | |
CN104199295B (zh) | 基于神经网络的机电伺服***摩擦补偿和变结构控制方法 | |
CN104238361B (zh) | 电机伺服***自适应鲁棒位置控制方法与*** | |
Xue et al. | On performance analysis of ADRC for a class of MIMO lower-triangular nonlinear uncertain systems | |
CN104638999B (zh) | 基于分段神经网络摩擦模型的双电机伺服***控制方法 | |
Garrido et al. | Inertia and friction estimation of a velocity-controlled servo using position measurements | |
CN104614984B (zh) | 一种电机位置伺服***的高精度控制方法 | |
CN110116409B (zh) | 一种基于扰动观测器的四通道遥操作双边控制方法 | |
CN104360635A (zh) | 一种电机位置伺服***的抗干扰控制方法 | |
Zhang et al. | Adaptive nonsingular terminal sliding mode control design for near space hypersonic vehicles | |
Zhang et al. | Dynamic feedrate optimization for parametric toolpath with data-based tracking error prediction | |
CN108762088B (zh) | 一种迟滞非线性伺服电机***滑模控制方法 | |
Su et al. | Fixed-time inverse dynamics control for robot manipulators | |
CN104007660B (zh) | 一种基于反演设计的伺服***的抗饱和控制方法 | |
CN115453890A (zh) | 一种齿轮巡检机器人***全局非奇异终端滑模控制方法 | |
CN111736472A (zh) | 一种基于rise的电机自适应预设性能渐近控制方法 | |
CN105843078A (zh) | 滑模控制方法及装置 | |
CN111752157B (zh) | 一种有限时间收敛的二阶滑模控制方法 | |
CN109995278B (zh) | 一种考虑输入受限的电机伺服***自调节控制方法 | |
CN109828586B (zh) | 基于非线性不确定性的无人船舶鲁棒h∞航向控制方法 | |
CN109194244B (zh) | 一种面向电动伺服***的控制方法及*** | |
CN108256214B (zh) | 一种航空正交加筋板刚度计算方法及装置 | |
Zhao et al. | Achieving minimum settling time subject to undershoot constraint in systems with one or two real right half plane zeros |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |