CN104614995B - 一种二阶系统有限时间滑模控制器的通用设计方法 - Google Patents

一种二阶系统有限时间滑模控制器的通用设计方法 Download PDF

Info

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
finite time
sliding
function
order system
time sliding
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
Application number
CN201510084429.8A
Other languages
English (en)
Other versions
CN104614995A (zh
Inventor
穆朝絮
徐伟
贾宏杰
孙长银
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin University
Original Assignee
Tianjin University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tianjin University filed Critical Tianjin University
Priority to CN201510084429.8A priority Critical patent/CN104614995B/zh
Publication of CN104614995A publication Critical patent/CN104614995A/zh
Application granted granted Critical
Publication of CN104614995B publication Critical patent/CN104614995B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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 &CenterDot; 1 = x 2 x &CenterDot; 2 = f ( x ) + b ( x ) u + d ( t )
其中,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的绝对值,幂指数α满足
CN201510084429.8A 2015-02-16 2015-02-16 一种二阶系统有限时间滑模控制器的通用设计方法 Active CN104614995B (zh)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5723582B2 (ja) * 2010-11-30 2015-05-27 インターナショナル・ビジネス・マシーンズ・コーポレーションInternational Business Machines Corporation スライディング・モード制御方式において、制御器を構成する方法、プログラム及びシステム

Patent Citations (4)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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) 一种基于干扰观测器的电机位置伺服系统自适应控制方法
CN104238361B (zh) 电机伺服系统自适应鲁棒位置控制方法与系统
CN104614984B (zh) 一种电机位置伺服系统的高精度控制方法
Liu et al. Soft variable structure controller design for singular systems
US10581270B2 (en) Equivalent-conductance-compensated globally-linear symmetric method for obtaining power flows in DC power networks
CN107121932A (zh) 电机伺服系统误差符号积分鲁棒自适应控制方法
CN104360635A (zh) 一种电机位置伺服系统的抗干扰控制方法
Feng et al. Non-singular terminal sliding mode control and its application for robot manipulators
CN106681343A (zh) 一种航天器姿态跟踪低复杂度预设性能控制方法
CN106325075A (zh) 一类时滞线性参数变化离散系统的h∞控制方法
Bi et al. Operator-based robust control for nonlinear uncertain systems with unknown backlash-like hysteresis
CN105093936A (zh) 电液位置伺服系统的有限时间连续滑模控制方法
CN103454921B (zh) 飞行控制系统非线性跟踪控制器设计的正切线性化方法
CN104570733B (zh) 一种含磁滞补偿的电机伺服系统预设性能跟踪控制方法
CN106094526A (zh) 一种广义预测工程化应用于脱硝控制系统的方法
CN105807613A (zh) 一种基于涡扇发动机多工作点的鲁棒补偿控制方法
Liu et al. Nonlinear system modeling and application based on system equilibrium manifold and expansion model
CN101782744A (zh) 基于小波变换的线性定常分布参数系统变结构控制方法
CN109976162A (zh) 一种三阶严反馈系统的全局非线性积分滑模控制方法
Cai et al. A super-twisting algorithm combined zeroing neural network with noise tolerance and finite-time convergence for solving time-variant Sylvester equation
Kirilov et al. Global analytic hypoellipticity for a class of evolution operators on T1× S3
CN109814393B (zh) 一种悬吊漂浮物随动系统的降维观测器和控制器设计方法
Wei et al. Non-asymptotic fractional pseudo-state differentiator for a class of fractional order linear systems
Chen et al. Low Chattering Trajectory Tracking Control of Non-singular Fast Terminal Sliding Mode Based on Disturbance Observer

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