张春燕 盛安冬 戚国庆 李银伢
环航是指机器人保持一定距离环绕目标的运动,其控制器需要使机器人运行在以目标为圆心、指定半径的期望轨迹上.环航跟踪在军用侦察监视、民用搜索救护等领域均有着广泛的应用前景.目前关于完整机器人对目标环航跟踪问题的研究已相对完善[1−6],为进一步便于工程应用,一些学者研究了非完整机器人的环航问题.根据其控制量的不同,主要分为以下两类:
1)将机器人速率视为常数,仅以角速率为控制量[7−21].其中Park[7]在环航研究中,令速率恒定并以向心加速度作为控制量,与角速率输入类似.文献[8−13]针对纯方位量测和纯距离量测,给出了环航跟踪系统的观测器与控制器设计方法,并应用到多目标的跟踪中.Lawrence[14]提出了一种基于Lyapunov矢量场引导法的环航控制器;另外,国内外研究学者将多种现代控制策略应用到环航跟踪研究中,并取得了一系列成果[15−21].
2)部分学者在设计非完整机器人角速率控制器的同时,增加了速率测量,即以速率与角速率作为机器人输入量[22−27],以满足更复杂的跟踪要求,但也增加控制器设计的难度.Summers等[22]在多机协同环航跟踪中考虑了速率变化,但在单机环航研究中仍视速率为常数;Oh等分别利用滑模控制[23]和模型预测控制[24]提出的环航控制方法在保证机器人与目标距离固定的同时,实现了以额定速率对目标进行环航跟踪.文献[25−27]则实现了环航系统中,机器人轨迹渐近收敛到相对距离和环航速率固定的期望轨迹上.目标静止时,将机器人速率固定为期望环航速率,即可利用上述控制方案达到跟踪要求,Yu等[25]亦给出了多机定距定速环航跟踪静止目标的分布式控制方法.对运动目标,Arranz等[26]提出的控制方案受目标与机器人初始位置限制;Miao等[27]针对系统直角坐标系和极坐标系下运动模型分别提出了控制器的设计方法,但未考虑机器人速度受限问题,会出现速率输入幅值偏大的问题.
上述研究中均假设机器人的运动能不受限制地满足控制器要求,然而在实际应用中,机器人速度普遍会受到限制.如何保证机器人在速度受限情形下达到期望的控制要求,并能在有限时间内运行在期望轨迹上,是工程应用中存在的两个现实问题.针对这些问题,本文研究了非完整机器人在速度受限下,以规定环绕速率沿固定半径跟踪已知目标的控制器设计问题.首先,针对机器人速度受限的情形,根据机器人运动模型及跟踪要求,利用反步法给出了使环航系统渐近稳定的控制器;为避免有限时间控制器设计时会出现控制量幅值过大的问题[28−30],进一步给出了一种在满足输入约束情形下的有限时间饱和控制器,使机器人轨迹在有限时间内收敛到期望轨迹.
本文结构安排如下:第1节描述了所研究的问题;第2节给出了两种控制器的设计方案;第3节给出了仿真比较结果;第4节为文章的结论.在文中,R表示实数集,R+表示正实数集,AT是向量或矩阵A的转置,表示向量x的欧几里得范数,sup表示上确界,AB表示集合{x|x∈A,x∈/B},sgn(x)表示符号函数(sgn(0)=0).
本文主要研究已知目标的速度信息情况下,非完整机器人对平面目标的环航控制器设计问题.且基于工程应用,机器人的速度受限.
将机器人视为质点,在二维笛卡尔坐标系下,机器人的运动模型为[26]:
目标在二维笛卡尔坐标系下的运动模型为:
根据环航跟踪的控制要求,对机器人和目标的相对位置建模.为简便书写,定义一个角度函数,
注1.当x=y=0时,y/x是不确定的,目前数学上常用的处理方法是令atan(0,0)为某一常值或在函数定义域中排除点(0,0)[31].在本文中,atan(y,x)函数不会出现x,y均为0的情况.以式(4)中atan(yo−yt,xo−xt)为例,yo−yt=xo−xt=0表示机器人与目标位置相同,在环航应用中可以排除这种情形.因此,考虑数学完备性和工程实用性,文中定义函数atan(·)时,直接令atan(0,0)=0.
如图1所示,在以目标为原点的极坐标下,系统相对位置可以表示为
图1 环航跟踪系统模型Fig.1 Model of the standofftracking system
其中,d∈R+表示机器人与目标之间的相对距离,θ∈(−π,π]表示机器人与目标之间的视线角.
问题描述:针对非完整机器人,设计速率控制器u1和角速率控制器u2,使其实现在速度受限条件下,以固定环绕速率沿指定半径跟踪目标.即设计u=[u1u2]T,实现
其中,dc∈R+表示规定距离,ωc∈R表示规定环绕角速率.
注2.环航跟踪系统中,固定环绕速率等同于环绕角速率固定,考虑到机器人顺时针或逆时针环绕目标,以环绕角速率代替环绕速率作为跟踪要求.ωc>0表示机器人逆时针环绕目标;ωc<0,机器人顺时针环绕目标.
假设1.基于环航跟踪要求,机器人速率vo与目标速度νt应满足以下条件:
注3.根据机器人环航要求可知,机器人速度与目标速度需满足. 由于实际应用中机器人速率受限制,并不能实现对任意运动目标的环航跟踪.在假设中给出了目标速度的条件,保证速率受限的机器人仍能完成对目标的环航跟踪,满足实际应用需求.
定义1[32].考虑如下连续系统:
引理1[33].若存在一个连续可微的Lyapunov方程V(xxx):D→R满足下以下两个条件:
2)存在k∈R+和α∈(0,1),使得对任意.则系统在平衡点有限时间稳定.且收敛时间.
引理2[34].对于x∈R,y∈R与常数p≥1,下列不等式成立
如图1所示,定义φ:=ψo−θ+π ∈(−π,π]表示机器人航向与视线之间的角度,即ψo=θ+φ−π;定义距离误差为实际距离与期望距离的差:ed=d−dc,ed∈(−dc,∞).
对目标静止情况分析,在极坐标下环航跟踪系统可以表示为
2.1.1 反步法控制器
如式(7)所述,系统状态方程中存在非微分同胚映射,本节考虑了如何在此情况下利用反步法设计控制器.另外,如何选取控制器参数使机器人速率满足约束条件的同时,提高系统收敛速度也是本节研究重点之一.
为方便控制器设计,构造一个非线性函数f(·):R→R,函数f(x)都需满足:
下文中的fs(·)和fm(·)均为满足上式条件的非线性函数,且函数绝对值的上界分别为,即,.
注4.根据式(8)所示条件,函数f(·)为有界的递增函数,且f(0)=0,存在满足条件的函数,如:等.
受文献[27]启发,用反步法按如下步骤设计控制器.
步骤1.将φd作为系统(7)中ed和θ子系统的虚拟控制量,即
如图1,由于相对距离变化率与相对速度νr在视线上的分量相关,因此令
为保证虚拟控制量φd到fs(ed)的映射为微分同胚映射,即u1cos(·)存在连续可微逆映射,需sinφd0.考虑到机器人环绕目标时的期望环绕角速率为ωc,可令
构造Lyapunov函数
其导数沿闭环系统(13)为
步骤2.定义相位误差,则.令,当时,ς→0.
构造Lyapunov函数
其中,k1>0,V2是正定的,其导数为
令机器人角速率控制器为
其中,k2>0.
定理 1.机器人速率与角速率控制器分别为式(11)和式(14)时,系统(7)中距离和相位误差(ed,eφ)渐近收敛到(0,0),且环绕角速率固定为ωc,达到式(5)所示的跟踪要求,即机器人稳定的环绕跟踪静止目标.
证明.由式(11)和式(14)可得:
Lyapunov函数V2的导数沿上述闭环系统为:
根据LaSalle不变集原理[35],距离和相位误差ed,eφ渐近收敛到不变集,即.当ωc>0时,φ=−π/2,机器人逆时针环绕目标;当ωc<0时,φ=π/2,机器人顺时针环绕目标.□
定理1给出了非完整机器人环航跟踪静止目标的渐近稳定控制器.与文献[27]中控制器相比,本文控制器考虑了机器人速率限制.以静止目标为例,文献[27]中环航跟踪控制器如下所示:
其中,k,k1,k2∈R+,φd=atan(−ωcdc,ked). 在此控制器下,相对速度νr在视线上的分量为u1cosφd=ked.当机器人速率受限时,参数k应满足
注5.为满足式(2)中机器人速率限制条件,参数k选取时,需假设一个相对距离误差的上界,使得.考虑到实际应用中距离d必然存在上下界的情况,即|ed|存在一个具有一定保守性的上界,上述假设是合理的.因此,参数k选取时应满足.定理1中控制器以fs(ed)替代了控制器(15)中参数ked,可以选取适当的函数fs(·),使得.由式(13)可知,ed的收敛速度与ked或fs(ed)有直接关系,在此情况下,收敛速度方面本文算法优于文献[27]方法.
2.1.2 有限时间稳定控制器
渐近稳定的系统状态随着时间趋向于平衡点,但不能保证在规定时间内达到期望值.进一步考虑设计机器人控制器,使相对距离在有限时间内收敛到期望距离dc,同时环绕角速率达到期望值ωc,即
其中,0 由于在进行有限时间控制器设计时通常需利用不等式进行放缩,所以控制器增益取值一般偏大,进而导致控制量幅值偏大[30].若直接利用有限时间控制器,可能会影响系统的性能,且无法保证满足式(2)中机器人速率限制条件,因此设计了输入受限下的有限时间饱和控制器(Finite-time saturated controller,FTSC). 定义函数 其中,x∈R,α∈(0,1),则sigα(·)为连续的递增函数,且 sigα(0)=0. 在函数f(·)的基础上,构造一个非线性函数h(·):R→R, 其中,c2>0,且. 函数h(·)同样满足式(8)中条件,即h(·)为有界连续的递增函数,且h(0)=0.下文中的hs(·)和hm(·)均为上述形式的函数,且函数绝对值的上界分别为. 设计控制器 定理 2.机器人控制器为式(18)和式(19)时,环航系统(7)中距离和相位误差(ed,eφ)在有限时间内收敛到(0,0),满足式(17)中跟踪要求,即机器人轨迹在有限时间内收敛到期望轨迹. 证明.定义一个非线性函数χ(·):R→R+,使得h(x)=χ(x)k3sigα(x),∀x∈R,则 当x有界时,必然存在,使得.考虑到距离误差ed有界,因此存在,使得对任意,都有. 环航系统(7)中控制器为式(18)和式(19)时, 根据LaSalle不变集原理[35],系统(7)在控制器(18)和(19)作用下在距离和相位误差为(0,0)处是渐近稳定的,进一步证明其有限时间稳定. 由于0<α<1,则(α+1)/2∈(0,1).根据引理2,可得 由引理1,距离和相位误差(ed,eφ)在有限时间内收敛到(0,0),且收敛时间满足. 综上所述,对任意ed∈(−dc,∞),φ∈(−π,π],在控制器(18)和(19)作用下,在有限时间内收敛到,即. □ 注 6.根据函数χ(·)的定义,当x>c2时,,因此,.可以由相对距离误差上界,给出一个具有一定保守性的下界.当工程应用中给出规定的收敛时间Tf时,通过选取适当的控制参数使得,即可保证机器人轨迹在规定收敛时间Tf内收敛到期望轨迹. 针对运动目标,极坐标下环航跟踪系统可以表示为 利用反步法设计非完整机器人控制器,实现对运动目标的环航跟踪.并针对不同的目标运动状态,给出参数选取的方法,使机器人速率满足约束条件的同时,提高收敛速度. 步骤1.将φd作为系统(20)中ed和θ子系统的虚拟控制量,令 由式(22)可得,|fm(·)|的上界时,;当时,;考虑如式(2)所示机器人速率限制,对于满足假设条件的运动目标,应满足 步骤2.对运动目标和静止目标的环航系统,均有子系统,因此令机器人角速率输入控制器为式(14),其中变量由式(22)和式(23)求得. 推论 1.对于满足假设条件的运动目标,针对环航系统(20)设计控制器为式(22)和式(14)时,系统在距离和相位误差为(0,0)处渐近稳定,满足式(5)中跟踪要求,实现对运动目标的环绕跟踪. 证明部分与定理1类似,此处不再赘述. 在系统渐近稳定的基础上,结合有限时间稳定理论设计控制器,使运动目标的环航跟踪系统在有限时间内稳定.令 推论 2.对于满足假设条件的运动目标,针对环航系统(20)设计控制器为式(25)和式(19)时,系统在距离和相位误差为(0,0)处有限时间稳定,满足式(17)中要求,机器人轨迹有限时间内收敛到期望轨迹. 证明部分与定理2类似,此处不再赘述. 本节针对静止目标、运动目标的跟踪环航问题进行了仿真实验.并通过与文献[27]中控制方法进行比较,说明本文所提方法的有效性和优越性. 为确保仿真实验的公平性,实验中选取的机器人仿真参数均一致:最大线速度,最小线速度;机器人初始位置均为,初始航向;目标初始位置均设为;规定跟踪距离dc=6 m,环绕角速率ωc=−0.5 rad/s.为方便选取控制器参数,假设相对距离误差上界为. 目标静止时,为使文献[27]中控制器速率输入满足式(2)所示的限制条件,速率控制器(15)中参数.而本文提出的考虑速率约束的反步控制器(11)和有限时间饱和控制器(18)中,c1=0,参数. 以环航跟踪静止目标为例,如注6所述,针对文中所提的有限时间饱和控制器,给定一个合理的收敛时间Tf,可以设计控制参数,使系统在规定时间内收敛到平衡点.根据仿真数据,令Tf=20 s.由,可得,.结合初始相对位置,可令k1=0.01,k2=5,则.为公平地比较分析三种控制器的不同点,本文所提考虑速度受限的反步法角速率控制器(14)和文献[27]角速率控制器中参数均设为k1=0.01,k2=5. 图2 文献[27]方法对静止目标的环航跟踪Fig.2 Standofftracking a static target using control strategy in[27] 图2为利用文献[27]中控制器时,机器人环航跟踪静止目标的仿真结果,其中图2(a)为目标位置和机器人轨迹,图2(b)和图2(c)分别为机器人速率和角速率输入量u1,u2.图3和图4分别为本文所提考虑速度受限时的反步控制方法和有限时间饱和控制方法下,机器人对静止目标的环航跟踪仿真结果.图5为三种控制器下对静止目标环航跟踪的相对距离d和视线角变化率θ的仿真结果对比图. 图3 反步法对静止目标的环航跟踪Fig.3 Standofftracking a static target using backstepping 由图5中对比文献[27]控制方法和本文考虑速率约束的反步控制方法,可看出,,即两种控制方法均能实现稳定跟踪.文献[27]中参数k需满足不等式(16)中条件,当较大时,必然导致k取值较小.图2(b)所对应的速率控制器(15)中参数k已使机器人初始速率非常逼近,即所选取参数k已逼近其可选范围的上界.本文提出的反步控制方法用非线性函数fs(ed)代替ked,在保证速率符合约束条件的同时,始终满足|ked|<|fs(ed)|.因此,如图5所示,本文所提的反步法较文献[27]中方法,系统收敛速度有明显的提高.另外,本文反步法控制器参数选取时,不需要考虑相对距离误差上界. 图4 有限时间饱和控制方法对静止目标的环航跟踪Fig.4 Standofftracking a static target using FTSC 由图5仿真结果对比本文反步控制方法和有限时间饱和控制方法,可以看出,反步法控制器能使系统渐近收敛到平衡点,而在有限时间饱和控制器下,,即机器人在有限时间饱和控制器下,能在规定收敛时间Tf内达到跟踪要求(17).另外,图4(c)中,角速率输入u2在t=14 s左右时,出现较大变化,主要是由于距离误差收敛到ed=0 m时,变化较大,进而导致u2变化. 图5 三种控制方法对静止目标的环航跟踪结果对比Fig.5 Standofftracking a static target with different control strategies 运动目标分为匀速目标和变速目标进行仿真分析: 图6为利用文献[27]中控制器时,机器人环航跟踪匀速目标的仿真结果,其中图6(a)为目标位置和机器人轨迹,图6(b)为相对距离d和视线角θ的变化曲线,图6(c)和图6(d)分别为机器人速率和角速率输入量u1,u2.图7和图8分别为机器人在本文所提考虑速度约束反步法和FTSC方法下,针对匀速目标的环航跟踪仿真结果.图9~11则为三种方法下对变速目标的环航跟踪仿真结果. 对于匀速目标vt= 1m/s,令文献[27]控制方法中参数为k= 0.075以满足,进而使机器人速率输入满足式(2)所示的限制条件p; 令本文所提反步√法和FTSC方法中,满足式(24)中不等式条件.而变速目标速率,令文献[27]控制方法中参数为k=0.06,本文所提反步法和FTSC方法中.控制器中其他参数与静止目标环航控制器参数一致. 图6 文献[27]方法对匀速目标的环航跟踪Fig.6 Standofftracking a constant velocity target using control strategy in[27] 图7 反步法对匀速目标的环航跟踪Fig.7 Standofftracking a constant velocity target using backstepping 由图7和10可以看出,对任意满足假设条件的运动目标,机器人在具有速度约束的反步法控制方案下均能达到跟踪要求(5),实现稳定跟踪.与图6和图9对比,机器人收敛速度明显较快.本文所提考虑速度约束反步法,较文献[27]所提方法,具有显著的优势. 由仿真参数计算可得,有限时间饱和控制器下,对匀速目标的环航跟踪系统,收敛时间应T≤30 s;对变速目标的环航跟踪系统,收敛时间T≤40 s.由图8和图11可以看出,机器人与目标相对距离d和环绕角速率均能在规定时间内达到期望值.即,对任意满足假设条件的运动目标,在本文所提有限时间饱和控制器下,机器人轨迹均能在一定时间内收敛到期望轨迹. 图8 有限时间饱和控制方法对匀速目标的环航跟踪Fig.8 Standofftracking a constant velocity target using FTSC 图9 文献[27]方法对变速目标的环航跟踪Fig.9 Standofftracking a time-varying velocity target using control strategy in[27] 图10 反步法对变速目标的环航跟踪Fig.10 Standofftracking a time-varying velocity target using backstepping 图11 有限时间饱和控制方法对变速目标的环航跟踪Fig.11 Standofftracking a time-varying velocity target using FTSC 本文研究了非完整机器人在二维平面内静止和运动目标的环航跟踪问题.针对静止目标,首先给出了一种考虑机器人输入约束的反步法控制器,解决了传统控制方法中存在速率输入过大的问题,使机器人以规定环绕速率沿期望距离跟踪目标.在此基础上,设计了有限时间饱和控制器,可根据应用中要求的收敛时间,设计合理的控制参数,实现速度受限情形下,机器人达到期望的跟踪要求.最后,将控制方法推广到对运动目标的跟踪中,并证明了所提方法的可行性.未来将考虑与协同控制结合,将文中控制方法运用到多机器人协同环航跟踪目标问题中.2.2 情形2:运动目标
3 仿真结果及分析
3.1 静止目标环航跟踪
3.2 运动目标环航跟踪
4 结论