CN111752158B - 一种有限时间收敛的二阶滑模控制方法 - Google Patents
一种有限时间收敛的二阶滑模控制方法 Download PDFInfo
- Publication number
- CN111752158B CN111752158B CN202010693051.2A CN202010693051A CN111752158B CN 111752158 B CN111752158 B CN 111752158B CN 202010693051 A CN202010693051 A CN 202010693051A CN 111752158 B CN111752158 B CN 111752158B
- Authority
- CN
- China
- Prior art keywords
- sliding mode
- order
- spacecraft
- mode control
- attitude
- 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
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
- Feedback Control In General (AREA)
Abstract
Description
技术领域
本发明涉及一种控制方法。属于自动化控制领域。
背景技术
滑模控制方法目前已经被广泛应用于航天器控制、飞行器控制、机电***控制等领域。按照***相对阶数,滑模控制方法可以分为一阶滑模控制方法、二阶滑模控制方法和高阶滑模控制方法。二阶滑模及高阶滑模相对于一阶滑模控制的显著优点是能够削弱抖振现象。
目前,基于终端滑模面,非奇异终端滑模面及快速非奇异终端滑模面设计的一阶滑模控制方法能够使***状态在有限时间内收敛到平衡点附近的邻域,但由于控制策略中存在切换函数,因此所得到的一阶滑模控制策略将使***出现明显的抖振现象。二阶滑模控制能够显著消除***中存在的抖振现象,但是由于终端滑模面和非奇异终端滑模面在进行微分运算时将产生奇异项,因此二阶滑模控制无法基于终端滑模面和非奇异终端滑模面进行设计,这就意味着,目前二阶滑模控制策略只能使被控对象(***)的状态渐进稳定到平衡点。
发明内容
本发明是为了解决目前存在的二阶滑模控制无法基于终端滑模面和非奇异终端滑模面进行设计问题。
一种有限时间收敛的二阶滑模控制方法,包括以下步骤:
针对被控对象建立含有不确定性的非线性二阶***:
其中,x1和x2是***的状态,且x1和x2及其各自的一阶导数及二阶导数均可测,为x1的一阶导数;b(x1,x2)和f(x1,x2)是已知的非线性函数,d(t,x1,x2)是未知的非线性扰动,u是***的输入控制信号,y是***的输出;t为时间;
设计对数双曲正切终端滑模面为
基于滑模面(2)设计二阶对数双曲正切滑模控制策略并进行控制。
进一步地,所述未知的非线性扰动d(t,x1,x2)满足约束|d(t,x1,x2)|<η,η是约束上限。
进一步地,所述基于滑模面(2)设计的二阶对数双曲正切滑模控制策略为:
进一步地,当被控对象为航天器姿态控制***时,含有不确定性的非线性二阶***的具体形式如下:
其中:表示航天器姿态的四元数,x1=q, D(t,x1,x2)=E(q)J-1d,B(x1,x2)=E(q)J-1;I3×3为单位矩阵;是航天器姿态四元数的标量部分,q=[q1 q2 q3]T是航天器的姿态四元数的矢量部分,和分别是q的一阶导数和二阶导数,ω是航天器的转动角速度是其一阶导数,ω×的定义方式与q ×相同;d是航天器受到的外部有界扰动;u是作用于航天器上的控制信号;J为转动惯量矩阵。
有益效果:
本发明解决了二阶滑模控制无法基于终端滑模面和非奇异终端滑模面进行设计问题。而且可以使***的状态在有限时间内收敛到平衡点附近的微小邻域内,且本发明的二阶控制策略具有更小的抖振。
附图说明
图1为控制策略(3)作用下,实施例一中的***状态到达滑模面的过程;图1的横坐标为时间,纵坐标为x1;
图2为控制策略(3)作用下,实施例一中的***状态沿滑模面滑动到平衡点附近的微小邻域的过程;图2的横坐标为时间,纵坐标为s;
图3为实施例一中的一阶控制策略的曲线;图3的横坐标为时间,纵坐标为u1;
图4为控制策略(4)作用下,实施例一中的***状态到达滑模面的过程;图4的横坐标为时间,纵坐标为x1;
图5为控制策略(4)作用下,实施例一中的***状态沿滑模面滑动到平衡点附近的微小邻域的过程;图5的横坐标为时间,纵坐标为s;
图6为实施例一中的二阶控制策略的曲线;图3的横坐标为时间,纵坐标为u2;
图7为实施例二的***状态收敛到滑模面的过程;
图8为实施例二中姿态跟踪误差收敛过程;
图9为实施例二中角速度跟踪误差的收敛过程;
图10为实施例二中控制策略(10)的输入信号。
具体实施方式
具体实施方式一:
本实施方式所述的一种有限时间收敛的二阶滑模控制方法,包括以下步骤:
针对被控对象建立含有不确定性的非线性二阶***:
其中,x1和x2是***的状态,且x1和x2及其各自的一阶导数及二阶导数均可测,为x1的一阶导数;对于一个二阶***而言,x1可以为位置或角度等信息,对应x2可以为速度或角速度等信息;b(x1,x2)和f(x1,x2)是已知的非线性函数,d(t,x1,x2)是未知的非线性扰动,d(t,x1,x2)满足约束|d(t,x1,x2)|<η,η是约束上限;u是***的输入控制信号,y是***的输出;t为时间;
设计对数双曲正切终端滑模面为
其中,k>0是由设计者根据性能要求设置的常值参数,p和q都是正奇数,且0<p/q<1;sgn(·)为符号函数;
当s=0时,公式(2)可以改写为
选取Lyapunov函数为
对V1求导可得
对于一个任意的、有限小的、包含原点的微小邻域|x1|<ε,都有下式成立:
那么,公式
可以改写为
证毕。
基于该滑模面(2)设计得到的一阶对数双曲正切滑模控制策略为:
要证明控制策略(3)作用下的***(1)的可在有限时间内到达s=0只需选定公式Lyapunov函数为如下公式即可:
基于滑模面(2)设计得到的二阶对数双曲正切滑模控制策略为:
要证明控制策略(4)作用下的***(1)的可在有限时间内到达s=0只需选定公式Lyapunov函数为如下公式即可:
这里的证明过程就是参照经典的超螺旋滑模,即Super-Twisting方法。
证毕。
一阶对数双曲正切终端滑模控制策略(3)和二阶对数双曲正切终端滑模控制策略(4)都可以使***(1)的状态x1和x2在有限时间内收敛到平衡点附近的微小邻域内,且控制策略(4)相比于控制策略(3)具有更小的抖振。
实施例一:
将具体实施方式一的方案对用于广义的二阶***,进行仿真,具体如下:
对数双曲正切终端滑模面(2)中的参数取为:k=5,p/q=3/5。
一阶对数双曲正切终端滑模控制策略(3)中的参数取为:η=0.1。在一阶对数双曲正切终端滑模控制策略(3)的作用下,***状态x1和x2到达滑模面的过程如图1所示,***状态x1和x2沿滑模面滑动到平衡点的过程如图2所示。一阶对数双曲正切终端滑模控制策略(3)的曲线如图3所示。结果表明,一阶对数双曲正切终端滑模控制策略可以使***(1)的状态x1和x2在有限时间内收敛到平衡点附近的微小邻域内;一阶对数双曲正切终端滑模控制策略(3)曲线平缓,无突变到无穷大的趋势(即无奇异项)。
二阶对数双曲正切终端滑模控制策略(4)中的参数取为:K1=0.1,K2=0.0014。在二阶对数双曲正切终端滑模控制策略(4)的作用下,***状态x1和x2到达滑模面的过程如图4所示,***状态x1和x2沿滑模面滑动到平衡点的过程如图5所示。一阶对数双曲正切终端滑模控制策略(4)的曲线如图6所示。结果表明,二阶对数双曲正切终端滑模控制策略可以使***(1)的状态x1和x2在有限时间内收敛到平衡点附近的微小邻域内;二阶对数双曲正切终端滑模控制策略(3)曲线平缓,无突变到无穷大的趋势(即无奇异项);控制策略(4)相比于控制策略(3)具有更小的抖振。
以传统的二阶线性超扭滑模控制策略作为对比,进行说明:
传统的二阶线性超扭滑模控制策略的形式为
其中
1、一阶对数双曲正切终端滑模控制策略(3)和二阶对数双曲正切终端滑模控制策略(4)都可以使***(1)的状态x1和x2在有限时间内收敛到平衡点附近的微小邻域内(即有限时间内到达某一精度,而非无穷时间);而控制策略(5)只能保证当时间无穷大的时候,***状态才能收敛到平衡点。
2、一阶对数双曲正切终端滑模控制策略(3)和二阶对数双曲正切终端滑模控制策略(4)均无奇异项。
3、控制策略(4)相比于控制策略(3)具有更小的抖振。
实施例二:
将具体实施方式一的方案对用于航天器姿态控制***,具体如下:
用姿态四元数描述姿态控制模型:
其中,表示航天器姿态的四元数,是航天器姿态四元数的标量部分,q=[q1 q2 q3]T是航天器的姿态四元数的矢量部分,I3×3为单位矩阵;和分别是姿态四元数的一阶导数和二阶导数,本实施例中,粗体的矢量形式q用于表示航天器的姿态四元数,区别于细体形式表示的正奇数q,也正是为了区分于正奇数q,航天器姿态四元数的标量部分用表示;ω是航天器的转动角速度,是其一阶导数,ω×的定义方式与q ×相同;d是航天器受到的外部有界扰动;u是作用于航天器上的控制信号;J为转动惯量矩阵;
显然,***(7)可以变化为如下形式:
由此可得用于航天器姿态控制的二阶对数双曲正切滑模控制策略具体为:
不失一般性,给定***状态的初值为转动惯量 ω=[0 0 0]Trad/s,控制策略中的参数选取为与实施例1中的参数一致,将(10)所得的控制策略作用于***(7)后,***状态收敛到滑模面的过程如图7所示,姿态跟踪误差和角速度跟踪误差的收敛过程分别如图8和图9所示,控制策略(10)的输入信号如图10所示。显然,***状态在有限时间内收敛到平衡点附近的邻域,且控制信号无高频抖振,也无奇异现象。
需要注意的是,具体实施方式和实施例仅仅是对本发明技术方案的解释和说明,不能以此限定权利保护范围。凡根据本发明权利要求书和说明书所做的仅仅是局部改变的,仍应落入本发明的保护范围内。
Claims (4)
1.一种有限时间收敛的二阶滑模控制方法,其特征在于,包括以下步骤:
针对被控对象建立含有不确定性的非线性二阶***:
其中,x1和x2是***的状态,且x1和x2及其各自的一阶导数及二阶导数均可测,为x1的一阶导数;b(x1,x2)和f(x1,x2)是已知的非线性函数,d(t,x1,x2)是未知的非线性扰动,u是***的输入控制信号,y是***的输出;t为时间;
设计对数双曲正切终端滑模面为
其中,k>0是常值参数,p和q都是正奇数,且0<p/q<1;sgn(·)为符号函数;
基于滑模面(2)设计二阶对数双曲正切滑模控制策略并进行控制;
所述基于滑模面(2)设计的二阶对数双曲正切滑模控制策略为:
2.根据权利要求1所述的一种有限时间收敛的二阶滑模控制方法,其特征在于,所述未知的非线性扰动d(t,x1,x2)满足约束|d(t,x1,x2)|<η。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010693051.2A CN111752158B (zh) | 2020-07-17 | 2020-07-17 | 一种有限时间收敛的二阶滑模控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010693051.2A CN111752158B (zh) | 2020-07-17 | 2020-07-17 | 一种有限时间收敛的二阶滑模控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111752158A CN111752158A (zh) | 2020-10-09 |
CN111752158B true CN111752158B (zh) | 2021-08-03 |
Family
ID=72711761
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010693051.2A Active CN111752158B (zh) | 2020-07-17 | 2020-07-17 | 一种有限时间收敛的二阶滑模控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111752158B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113064349B (zh) * | 2021-03-22 | 2022-06-07 | 中国人民解放军国防科技大学 | 固定时间收敛的非线性平台滑模控制方法、装置及*** |
CN115407788B (zh) * | 2022-08-26 | 2023-09-15 | 哈尔滨工业大学 | 一种水下高速航行体固定时间收敛二阶滑模控制方法及控制*** |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104467595A (zh) * | 2014-12-05 | 2015-03-25 | 沈阳工业大学 | 直接驱动伺服***的二阶滑模控制***及其控制方法 |
CN104485864A (zh) * | 2014-11-24 | 2015-04-01 | 沈阳工业大学 | 直接驱动伺服***的二阶滑模控制***及其控制方法 |
JP2017157121A (ja) * | 2016-03-04 | 2017-09-07 | 日産自動車株式会社 | スライディングモード制御方法及びスライディングモード制御装置 |
CN109557814A (zh) * | 2018-12-04 | 2019-04-02 | 四川航天***工程研究所 | 一种有限时间积分滑模末制导律 |
CN110161852A (zh) * | 2019-05-08 | 2019-08-23 | 杭州电子科技大学 | 一种基于二阶滑模算法的移动机械臂运动控制方法 |
CN110487132A (zh) * | 2019-08-30 | 2019-11-22 | 哈尔滨工业大学 | 一种基于非奇异快速终端滑模控制的角度约束制导方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8436283B1 (en) * | 2008-07-11 | 2013-05-07 | Davidson Technologies Inc. | System and method for guiding and controlling a missile using high order sliding mode control |
-
2020
- 2020-07-17 CN CN202010693051.2A patent/CN111752158B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104485864A (zh) * | 2014-11-24 | 2015-04-01 | 沈阳工业大学 | 直接驱动伺服***的二阶滑模控制***及其控制方法 |
CN104467595A (zh) * | 2014-12-05 | 2015-03-25 | 沈阳工业大学 | 直接驱动伺服***的二阶滑模控制***及其控制方法 |
JP2017157121A (ja) * | 2016-03-04 | 2017-09-07 | 日産自動車株式会社 | スライディングモード制御方法及びスライディングモード制御装置 |
CN109557814A (zh) * | 2018-12-04 | 2019-04-02 | 四川航天***工程研究所 | 一种有限时间积分滑模末制导律 |
CN110161852A (zh) * | 2019-05-08 | 2019-08-23 | 杭州电子科技大学 | 一种基于二阶滑模算法的移动机械臂运动控制方法 |
CN110487132A (zh) * | 2019-08-30 | 2019-11-22 | 哈尔滨工业大学 | 一种基于非奇异快速终端滑模控制的角度约束制导方法 |
Non-Patent Citations (3)
Title |
---|
Adaptive Nonsingular Fast Terminal Sliding Mode Guidance Law with Impact Angle Constraints;Junhong Song,等;《International Journal of Control , Automation and Systems》;20161231;第14卷(第1期);第99-114页 * |
Non-singularTerminal Sliding Mode Control for Fast Circumnavigate Mission Centered on Non-cooperative Spacecraft;Xuebo Yang,等;《IEEE》;20191231;第3658-3663页 * |
近圆轨道下航天器相对运动的滑模控制;柴鑫彤,等;《电子设计工程》;20170131;第25卷(第1期);第93-96页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111752158A (zh) | 2020-10-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107577144B (zh) | 一种基于增强型指数趋近律的飞行器有限时间自适应姿态控制方法 | |
Fang et al. | A novel sliding mode controller for small-scale unmanned helicopters with mismatched disturbance | |
CN103728882B (zh) | 微陀螺仪的自适应反演非奇异终端滑模控制方法 | |
CN111752158B (zh) | 一种有限时间收敛的二阶滑模控制方法 | |
CN110018637B (zh) | 一种考虑完成时间约束的航天器姿态跟踪保性能控制方法 | |
Wu et al. | Quantized fixed-time fault-tolerant attitude control for hypersonic reentry vehicles | |
CN111752157B (zh) | 一种有限时间收敛的二阶滑模控制方法 | |
CN110543184B (zh) | 一种刚性飞行器的固定时间神经网络控制方法 | |
CN109634291B (zh) | 一种基于改进型障碍李雅普诺夫函数的刚性飞行器姿态约束跟踪控制方法 | |
CN110286595B (zh) | 一类受饱和非线性输入影响的分数阶***自适应控制方法 | |
CN104122794A (zh) | 微陀螺仪的自适应模糊神经补偿非奇异终端滑模控制方法 | |
CN109062240B (zh) | 一种基于神经网络估计的刚性飞行器固定时间自适应姿态跟踪控制方法 | |
CN113589689B (zh) | 一种基于多参数自适应神经网络的滑模控制器设计方法 | |
CN110501911A (zh) | 一种考虑执行器受限问题的刚性飞行器自适应固定时间姿态跟踪控制方法 | |
CN106113046B (zh) | 基于死区和摩擦补偿的机械臂伺服***动态面瞬态控制方法 | |
CN112631316B (zh) | 变负载四旋翼无人机的有限时间控制方法 | |
Hu et al. | Attitude output feedback control for rigid spacecraft with finite-time convergence | |
CN108762088B (zh) | 一种迟滞非线性伺服电机***滑模控制方法 | |
CN104503246A (zh) | 微陀螺仪***的间接自适应神经网络滑模控制方法 | |
CN115256386B (zh) | 考虑跟踪误差约束的不确定机械臂神经自适应控制方法 | |
CN113110048A (zh) | 采用hosm观测器的非线性***输出反馈自适应控制***和方法 | |
CN110488854B (zh) | 一种基于神经网络估计的刚性飞行器固定时间姿态跟踪控制方法 | |
CN115981162A (zh) | 一种基于新型扰动观测器的机器人***滑模控制轨迹跟踪方法 | |
Gu et al. | Model free adaptive control design for a tilt trirotor unmanned aerial vehicle with quaternion feedback: Theory and implementation | |
CN105116905A (zh) | 一种飞行器姿态控制方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |