摘要
针对一类带有输出误差约束和外界干扰的非线性系统, 结合固定时间干扰观测器, 本文设计滤波反步控制策略. 首先, 设计预定性能函数, 处理输出误差约束问题. 其次, 利用固定时间滤波器和固定时间干扰观测器, 分别解决反步控制的“计算爆炸”问题和外界干扰问题, 干扰估计误差在固定时间内收敛到零. 再次, 采用Lyapunov理论分析闭环系统是固定时间稳定的, 且输出误差在预设时间内满足约束要求. 最后, 通过比仿真验证控制策略的有效性.
Abstract
Considering a class of nonlinear systems under output error constraints and external disturbances, a filtered backstepping control strategy is proposed by combining with a fixed-time disturbance observer. Firstly, the prescribed performance function is designed to solve the problem of output error constraint. Secondly, a fixed time filter and a fixed-time disturbance observer are designed to respectively tackle the “explosion of complexity” problem and the external disturbances problem, and to ensure that the disturbance estimation error converges to zero within a fixed time. The Lyapunov theory is utilized to analyze the fixed time stability of the closed-loop system, at the same time, the output error meets the constraint requirements in the preset time. Finally, comparative simulation is used to verify the effectiveness of the control strategy.
1 引言
近年来,非线性系统控制及其在工程系统中的应用引起了越来越多的关注,如机器人、无人机、高超声速飞行器等系统 [1-3] . 反步控制具有系统化和结构化的优势,被广泛应用于非线性系统控制问题中,并取得了大量的研究成果. 当利用传统反步控制进行高阶非线性系统的控制器设计时,会存在“计算爆炸”问题. 滤波器是一种信号处理技术,广泛应用于通信、雷达等领域. 它的主要作用是从含有多种频率成分的信号中提取出特定频率成分,有效的抑制噪声,提高系统的性能. 为了解决“计算爆炸”问题,一阶滤波器 [4-5]、指令滤波器 [6-8] 和固定时间滤波器 [9-14] 应用于反步控制中. 文献 [12] 构造了自适应固定时间滤波器,且滤波误差在固定时间内收敛到关于零的任意小邻域. 文献 [13] 采用固定时间滤波器有效处理了“计算爆炸”问题,实现了滤波误差的收敛时间不依赖于初始条件.文献 [14] 设计了快速收敛固定时间滤波器,保证滤波误差获得更快收敛速度. 本文将设计固定时间滤波器,对反步控制的虚拟控制输入信号进行滤波处理,克服反步控制的“计算爆炸”问题.
工程系统中,被控系统往往受到外界干扰的影响,可能导致控制系统性能的下降. 水下机器人在浅水区域执行任务时容易受到风浪的影响,在深水区域则容易受到水流影响 [1] . 四旋翼无人机由于体积小、质量轻,容易受到外界环境的影响,如风的干扰 [2] . 高超声速飞行器所处的复杂环境导致其容易受到外界干扰的影响. 解决外界干扰问题常用的方法之一是设计干扰观测器. 干扰观测器将观测到的干扰值用于被控系统的设计中,减小外界干扰对被控系统的影响,提高被控系统的控制性能. 文献 [15] 采用积分滑模干扰观测器估计复合干扰,实现了估计误差的有限时间收敛. 文献 [16] 针对一类受未知干扰影响的非线性系统,设计了改进的固定时间干扰观测器. 文献 [3] 提出了基于滑模理论的固定时间干扰观测器. 文献 [17] 设计了固定时间干扰观测器,实现对干扰的有效估计. 文献 [18] 利用固定时间干扰观测器对不确定性进行估计,但是文献 [16-18] 没有考虑输出误差约束问题. 本文将设计固定时间干扰观测器估计外界干扰,估计误差在固定时间内收敛到零,从而提高控制系统的性能.
高性能控制是实际工程系统对其控制系统的一般要求,这也使得它们一般都要求输出误差满足一定的约束要求. 例如,高超声速飞行器巡航段的控制系统一般要满足: 高度和速度稳定跟踪其参考信号,同时速度跟踪误差和高度跟踪误差保持在约束边界内. 另外,永磁同步电机的角速度控制系统,要实现角速度的稳定跟踪,同时要保证角速度跟踪误差始终在约束边界内. 解决输出误差约束问题的常用方法之一是构造预定性能函数. 文献 [19] 针对非线性系统输出误差的预定性能问题,构造了预定性能函数. 考虑输出误差约束影响下的严反馈非线性系统,文献 [20] 提出了预定性能反步控制策略. 文献 [21] 针对永磁同步电机,设计了有限时间性能函数,其收敛时间只与设计参数有关. 文献 [22] 引入预定性能函数,并提出固定时间预定性能控制策略. 文献 [23] 在误差转换的基础上,设计了固定时间预定性能控制策略. 文献 [24] 提出了固定时间预定性能反步控制策略,使输出误差保持在预设的范围内,但文献 [20,23-24] 没有考虑外界干扰问题.
基于以上分析,本文的主要贡献如下: 1)设计新型预定性能函数,提出固定时间反步控制器,实现闭环系统在固定时间内稳定,且输出误差在预设时间内满足约束; 2)为了处理外界干扰,结合终端滑模理论,构造干扰观测器,实现在固定时间内对外界干扰的有效估计,同时估计误差收敛到零; 3)构造固定时间滤波器处理“计算爆炸”问题,与一阶滤波器 [4-5] 不同的是,它可以保证滤波误差在固定时间内稳定,且不用选择较小的时间常数,也可以获得较好的估计精度.
2 问题描述
考虑如下非线性系统:
(1)
其中:是状态向量; di(t)表示系统受到的外界干扰; 分别表示状态、控制输入和输出; 和是已知的光滑函数.
控制目标: 设计滤波反步控制策略,保证闭环系统在固定时间内是有界稳定的; 输出y稳定跟踪参考指令,输出误差在预设时间内满足约束.
为了达到控制目标,引入如下假设和引理:
假设 1 参考指令x1d和均已知且有界.
假设 2 存在常数gia,gib >0,满足
假设 3 对于虚拟控制输入xic,存在未知常数
引理 1 [12] 考虑系统若存在正定函数Y (x),满足其中: α,β,Υ >0,0 <υ <1,ω >1,0 <δ <1,则称系统是实际固定时间稳定,收敛时间 T1 满足:
若Y(x)满足:则称系统是固定时间稳定的,且收敛时间T2满足:
引理 2 [17] 考虑不等式有
其中:
引理 3 [12] 存在有
符号说明: sigr(·)= | · |rsgn(·),sgn(·)为符号函数.分别用来表示.
3 固定时间滤波器设计
为了解决“计算爆炸”问题,构造如下固定时间滤波器:
(2)
其中:是滤波器的输入信号,且为待设计的虚拟控制输入; xid是滤波器的输出信号,用于之后的第4节观测器设计和第5节控制器设计中; τi是滤波器的时间常数,且τi >0; 滤波误差为µi = xid − xic.
滤波误差µi的稳定性可以由如下定理概括.
定理 1 考虑固定时间滤波器式(2),若输入信号则滤波误差 µi 在固定时间内收敛到零. 若则在假设3的条件下,滤波误差µi在固定时间内有界稳定.
证选取Lyapunov函数为关于时间的导数为
(3)
若上式可化为
(4)
则µi在固定时间Tf内收敛到零,且
若由Young不等式,式(3)可化为
(5)
当时,有因此
(6)
其中.
由式(6)和引理1可以得到µi有界稳定,收敛域为
当时,有因此

(7)
则滤波误差µi有界稳定,收敛域为
证毕.
注 2 文献 [4-5] 的一阶滤波器通过减小时间常数τi加快滤波误差的收敛速度,但滤波误差无法在固定时间内实现稳定,且较小的τi会受到采样周期和高频噪声的限制. 相较于文献[9]提出的固定时间滤波器,本文滤波器式(2)可以实现滤波误差在固定时间内收敛的,同时,在一定程度上,估计精度比文献 [9] 的滤波器的估计精度高.
4 固定时间干扰观测器设计
为了处理外界干扰,基于非奇异快速终端滑模理论,设计如下固定时间干扰观测器:

(8)
式中:是di的估计,估计误差为待设计参数. si表示如下的非奇异快速终端滑模面:
(9)
其中:由滤波器式(2)估计得到,滤波误差为.
估计误差的稳定性可以由如下定理概括.
定理 2 对于固定时间干扰观测器式(8)的干扰估计误差在时收敛到零,Td满足
其中Tzf是滤波误差zf的收敛时间,
证根据式(1)(8),可得
(10)
由定理1可知,当时,式(9)可化为
(11)
将式(10)代入式(11),则对si求导得
(12)
选取Lyapunov函数为,对Ys求导得
(13)
由引理1,在时收敛到零,Td满足
证毕.
5 预设时间反步控制器设计
本节首先设计预定性能函数,然后基于固定时间干扰观测器,设计固定时间预定性能反步控制器.
5.1 预设时间预定性能函数设计
为了确保输出误差z1在预设时间内满足约束,设计如下性能函数:

(14)
其中:是预设时间,上式关于时间的导数为

(15)
可知因此,性能函数有如下两个特性:
1)是连续函数,即P(t)在预设时间Tc内收敛到
2)P(t)是单调递减函数,时,P(t)保持在不变.
1)文献 [19] 的性能函数为
p0 >p∞ >0,a2 >0,P2(t)在t → +∞时趋于p∞,收敛时间不能预先确定.
2)文献 [26] 的性能函数为

(16)
其中: a1,b1,c1 >0,PT1 = P(T1)>0. 其导数为

(17)
观察式(14)(16)可知,式(16)的预设时间T1由参数b1,c1设定,当b1和c1满足时,输出误差在预设时间T1内收敛到约束区域.
考虑如下3种性能函数: 情形1: 文献 [19] 的传统性能函数; 情形2: 文献 [26] 的性能函数式(16); 情形3: 本文的性能函数式(14). 在相同的初始条件 P(0)= e + 3,PTc = 4,Tc = 0.5 s和参数a = a1 = c = c1 = 1,b = b1 = a2 = 2下,对比曲线如图1所示. 由图1可知,在相同参数和预设时间下,0 <t <Tc内,本文的性能函数式(14)表现出更快的收敛速度.
图13种性能函数对比曲线
Fig.1Compared curves of three performance functions
5.2 控制器设计
步骤 1 输出误差为z1 =x1−x1d,z1的导数为
(18)
设计如下误差变换:
(19)
其中ξ1 = z1/P,对e1求导得
(20)
其中
设计如下虚拟控制输入:
(21)
其中:是d1的估计,由固定时间干扰观测器式(8)估计得到. 将式(21)代入式(20),可得
(22)
步骤 2 定义z2 =x2−x2d,x2d由滤波器式(2)提供. z2的导数为
(23)
设计如下虚拟控制输入:
(24)
其中:是d2的估计,由固定时间干扰观测器式(8)估计得到. 将式(24)代入式(23),可得
(25)
步骤定义zi =xi−xid,xid由滤波器式(2)提供. zi的导数为
(26)
设计如下虚拟控制输入:
(27)
其中:是di的估计,由固定时间干扰观测器式(8)估计得到. 将式(27)代入式(26),可得
(28)
步骤 n 定义zn =xn−xnd,xnd由滤波器式(2)提供. zn的导数为
(29)
设计实际控制输入u
(30)
其中:是dn的估计,由固定时间干扰观测器式(8)估计得到. 将式(30)代入式(29),有
(31)
6 稳定性分析
根据Lyapunov理论,证明闭环系统的稳定性.
定理 3 考虑由非线性系统(1)和实际控制输入式(30)组成的闭环系统,满足: 所有信号在固定时间内有界稳定,输出误差z1在预设时间内进入到预设区域(−kLP(t),kUP(t)).
证选取Lyapunov函数为
(32)
其中:
由定理2可知,当时,估计误差收敛到零,即.
在假设3的条件下,由Young不等式, Y 的导数为
根据引理2知,满足
由引理3可知
进一步,可得

(33)
其中:
A1 =min{q11−qi1−qn1},B1 =min{q12 −,qi2−qn2},C =K1 =−1,B2 =min{B1,K1},A2 = 2 min{},且A1,B1,C,K1,A2,B2 >0.
存在常数使可得
(34)
令则
由引理1可知, Y 在固定时间内收敛,收敛区域为{Y :}. 收敛时间T满足
即闭环系统在固定时间内有界稳定. 由式(19)可知,转换误差e1满足
进一步可以得到
因此,输出误差满足由 P(t)的表达式可知,z1在预设时间内进入到约束区域(−kLP(t),kUP(t)). 证毕.
注 5 文献 [20-21] 考虑了输出误差约束问题,但是没有考虑输出误差的固定时间收敛性. 文献 [23-24] 考虑了输出误差约束问题,但没有考虑外界干扰问题. 本文综合考虑输出误差约束和外界干扰问题,设计基于终端滑模理论的固定时间干扰观测器(8)估计外界干扰; 构造预设时间预定性能函数(14),解决了输出误差约束问题.
7 仿真验证
7.1 滤波器的对比仿真
考虑如下两种滤波器的输入信号:
1) xic = 3; 2) xic = sin t + sin (0.5t) .
对比仿真中,选取相同的初始条件和滤波器参数,第1种输入信号下,滤波器参数为li = 10,τi = 0.5,第2种输入信号下,滤波器的参数为图2和图3为滤波器的对比结果,由图可知,本文的滤波器式(2)在3个滤波器中波动幅度最小,估计精度相对较好.
图2情形1的对比结果
Fig.2Compared results in case1
图3情形2的对比结果
Fig.3Compared results in case2
注 6 对不同结构的滤波器,文献 [10-12] 通过选取相同的时间常数τi进行比较,以说明所设计滤波器的优越性. 本文在对不同结构的滤波器进行对比仿真时,针对两种类型的滤波器输入信号进行仿真,以说明本文设计滤波器的优越性. 另外,若选取相同的较小时间常数τi = 0.01,本文设计的滤波器的估计精度较高; 若选取不同的时间常数τi ,且一阶滤波器的时间常数选取的较小,即τi = 0.01,在第1种输入信号下固定时间滤波器 [9] 和本文的滤波器的时间常数选取为τi = 0.5时,一阶滤波器的估计精度较高,但在第2种输入信号下固定时间滤波器 [9] 和本文的滤波器的时间常数选取为τi = 0.2 时,本文设计的滤波器的估计精度较高. 再者,进行了多组时间常数的仿真后,结果表明本文设计的滤波器的收敛速度和精度不仅受到时间常数τi的影响,还受到vi1,vi2的影响,放松了对较小的时间常数的要求,且滤波器误差可以在固定时间内实现稳定. 因此,本文的滤波器式(2)不仅可以保证滤波误差在固定时间内收敛,在一定程度上,还具有较高的估计精度.
7.2 仿真验证
考虑如下范德堡尔电路系统 [25] :
(35)
其中:选取 x1(0)= 0.2,x2(0)= 0,参考指令为 x1d = sin(2t),外界干扰为 di(t)= 2 sin(0.1πt),i = 1,2.
控制器参数如表1所示,在相同的初始条件和控制器参数下,考虑如下5种对比情形:
情形1: 基于本文滤波器式(2),采用性能函数式(14)的反步控制策略; 情形2: 基于自适应固定时间滤波器 [9],采用性能函数式(14)的反步控制策略; 情形3: 基于一阶滤波器 [4-5],采用性能函数式(14)的反步控制策略; 情形4: 基于本文滤波器式(2)和文献 [19] 的传统性能函数的反步控制策略; 情形5: 基于本文滤波器式(2)和文献 [26] 的性能函数(16)的反步控制策略.
表1控制器参数
Table1Controller parameters
图4–5为情形1–3的对比结果. 图4描述了3种情形下的输出x1的跟踪曲线和x2的变化曲线,可以看出情形1的输出可以以更快的速度实现对参考指令的跟踪; 图5为输出误差z1的变化曲线和控制输入u的变化曲线,由图可知,z1收敛到零附近的小邻域; 情形1的u 表现出较快的稳定速度,波动幅度较小.
图4输出x1的跟踪曲线和x2的曲线
Fig.4Tracking curves of output x1 and curves of x2
图5输出误差z1和控制输入u的曲线
Fig.5Curves of output error z1 and control input u
图6–7为情形1、情形4和情形5的对比结果. 图6为输出x1的跟踪曲线和x2的变化曲线. 可以看出,输出 x1稳定跟踪参考指令x1d. 图7为输出误差z1的曲线和控制输入u的变化曲线. 在本文所设计的控制策略下,情形1中的跟踪精度相较于情形4和情形5的跟踪精度更优.
8 结论
本文综合考虑外界干扰和输出误差约束对非线性系统的影响,提出滤波反步控制策略. 基于终端滑模理论,设计固定时间干扰观测器估计外界干扰,同时利用固定时间滤波器处理“计算爆炸”问题,然后设计固定时间反步控制器,证明了闭环系统的固定时间稳定性,滤波器对比仿真和范德堡尔电路系统的对比仿真体现了所设计控制策略的有效性.
图6输出x1的跟踪曲线和x2的曲线
Fig.6Tracking curves of output x1 and curves of x2
图7输出误差z1和控制输入u的曲线
Fig.7Curves of output error z1 and control input u