无人水面艇(Unmanned Surface Vehicle,USV)作为一种小型水面智能任务平台,具有体积小、成本低、高速智能、雷达反射面积小、无人员伤亡等优势,在军用和民用领域具有广泛的应用前景[1]。同时,USV还是整合低空无人机和水下机器人跨网络的关键节点,备受世界各国科研人员的关注[2]。路径跟踪是USV完成各项使命任务的关键技术之一,定义为USV在控制系统的驱动下,不考虑时间约束,从任意初始位置出发,驶入预先设定的期望航线,并沿此航线最终到达指定位置[3]。不具有侧向推进器或侧向推进器在高速航行时失效的USV是一类典型的欠驱动系统,由于控制输入个数少于其所需控制的自由度个数,故该类系统存在不可积的加速度约束[3]。目前国内外对于欠驱动USV路径跟踪控制已有较多的研究成果。
Caharija等设计了基于坐标变换的全局k-指数直线路径跟踪控制器和海流观测器,但并未考虑更为一般的曲线路径跟踪控制问题[4]。田勇等基于视线导引策略和抗饱和PID控制算法设计了直线航迹跟踪控制器,并进行湖上试验验证,但并未考虑外界扰动的影响[5]。Fossen和Lekkas考虑海流等外界干扰对USV运动学模型的扰动, 进行了路径跟踪控制研究,但并未考虑USV模型的非对称特性[6]。Børhaug等提出了用于补偿外界环境干扰的积分视线(Integral Line-of-Sight,ILOS)导引策略,并给出系统的稳定性证明[7]。文献[8]设计了基于相对速度模型的路径跟踪ILOS导引策略,并证明了控制系统是全局k-指数稳定(Global k-Exponential Stable, GKES)的,但只实现直线路径跟踪。
本文在研究国内外现有成果的基础上,提出了2种改进ILOS导引策略,之后综合考虑USV的欠驱动性和不对称性以及外界扰动的影响,基于改进ILOS导引策略和反馈控制思想实现了USV水平面内的直线和曲线路径跟踪控制,并利用级联系统理论证明了当所有控制目标实现时,控制系统为全局k-指数稳定的。仿真对比实验表明了算法的有效性和先进性。
1 USV运动数学模型及问题描述 1.1 USV运动建模欠驱动USV的水平面三自由运动学和动力学模型[3]如下:
(1) |
式中:η=[x, y, ψ]T为大地坐标系i下USV在水平面的位置(x, y)和航向角ψ;J(ψ)为旋转矩阵;v=[u, v, r]T为USV的纵向、横向线速度(u、v)及航向角速度(r);f=[Tu, Tr]T为控制输入矩阵,Tu和Tr分别为前进推力和转艏力矩;B为控制输入配置矩阵;M为惯性参数矩阵;C(v)=CRB(v)+CA(v), CRB(v)和CA(v)分别为刚体自身和附加的科氏力和向心力矩阵; D(v)为阻尼参数矩阵。各矩阵定义如下:
(2) |
(3) |
(4) |
(5) |
(6) |
上述USV模型中没有考虑作用在艇体上的外界干扰,因此以上模型中的速度均为绝对运动速度。考虑存在海流等外界干扰的情况,若满足以下2个条件:
条件1 在大地坐标系i中,海流速度Vc=[Vx, Vy, 0]T为常值无旋,合速度
条件2 CRB(v)与速度u、v无关,即CRB(v)=CRB(vr), 其中vr为相对运动速度。
此时可通过在USV运动数学模型中加入海流速度Vc=[Vx, Vy, 0]T的方式,直接用相对运动速度vr=[ur, vr, r]T替代式(1)中的绝对运动速度v=[u, v, r]T,模型如下:
(7) |
模型式(7)不仅可避免同时使用绝对前进速度u和相对前进速度ur给控制器设计带来的复杂公式推导,而且由于Tu作用于船体产生的速度为ur,而不是u,因此可通过控制总能量消耗直接控制USV的相对前进速度ur。
1.2 USV模型坐标变换为便于2.3节欠驱动USV航向、航速控制律的设计,本文将船体坐标系的原点变换到船体枢心位置处,去除偏航力矩对横向运动的直接影响[9]。
定义:
(8) |
式中:ε=-(m33b22-m23b32)/(m22b32-m23b22)。坐标变换后,式(7)可写为
(9) |
(10) |
(11) |
(12) |
(13) |
(14) |
其中:
(15) |
(16) |
(17) |
(18) |
由式(13)可以看出,经过坐标变换,
考虑计划航线为由一系列航路点(xp(θ), yp(θ))相连组成期望路径P(θ),其中θ为大于0的路径参数变量,且对于非闭合曲线路径θ满足
(19) |
为便于USV路径跟踪误差的描述,基于Serret-Frenet坐标系来描述路径跟踪误差,如图 1所示。该坐标系以期望路径P(θ)上的自由点为原点,横轴XP沿P(θ)的切线方向,纵轴YP沿P(θ)的法线方向并指向右侧;xe(t)为USV的纵向位置误差;Plos(t)为视线导引矢量,t为时间;v为坐标变换后USV的横向速度;Δ=nL为USV的前视距离,n=2~5, L为船长。
在大地坐标系i中,USV当前位置为P(t)=[x, y]T,由图 1得横向位置误差
(20) |
对ye(t)微分得
(21) |
由式(19)可得
(22) |
由图 1所示的几何关系可知:
(23) |
将式(7)、式(22)及式(23)代入式(21)中得
(24) |
式中:β=arctan(u, v)为USV转向时造成的横向速度与外界干扰环境海流等共同作用产生的漂角,这将导致USV的航迹角χ与航向角ψ不相等,且χ=ψ+β。在船舶路径跟踪应用中,χ一般可通过导航系统测量得到,那么可将漂角对路径跟踪的影响直接映射到航迹角中,设计相应的航迹控制器去进行路径跟踪。考虑到实际USV路径跟踪应用中,未安装导航系统或导航系统失灵情况下,仍可根据USV上安装的电子罗盘量测到的航向以完成路径跟踪任务,此时应考虑设计航向控制器,故控制器设计过程中对漂角的处理方式就显得尤为重要。一般USV航行时都满足条件3。
条件3 漂角β很小(一般小于5°)且为常值。即存在βmax>0,使得|β|≤βmax,且
那么,式(24)可改写为
(25) |
本文设计的路径跟踪控制系统采用ILOS导引子系统与航向控制子系统级联的方式。导引子系统根据USV的位置、姿态信息和计划航线的信息,求解出期望的USV航向,而航向控制子系统则根据导引子系统给定的期望信号和当前USV的运动信息计算出期望的舵机操纵信号,以此使得USV路径跟踪的位置和航向都收敛到平衡位置,同时航速控制子系统实现USV对期望速度的跟踪。
因此本文的3自由度欠驱动USV路径跟踪控制目标为
(26) |
(27) |
(28) |
式中:ψd为导引子系统给出的航向参考信号,即期望航向角;ud为期望的相对前进速度。
2.1 传统ILOS导引策略ILOS导引策略是一种经典有效的导航算法,文献[3]中定义了船舶运动控制领域中视线导引的概念:ILOS矢量是船舶当前位置在计划航线投影点的切线上一点与船舶当前位置的连线。通过控制船舶合速度的方向始终对准ILOS矢量,引导船舶逼近计划航线。ILOS导引原理独立于动力学控制,不依赖于任何模型,对高频白噪声敏感度低,且需要设计的参数少,期望航向的获取只与船舶的实时位置和给定的期望航线有关,能够实时高效地计算期望航向并传递给控制层。ILOS导引策略可分为基于包围圈和前视距离2种,本文采用基于前视距离的ILOS导引策略来完成路径跟踪目标,传统ILOS导引策略[7]定义如下:
(29) |
(30) |
式中: yint为虚拟控制输入,是一个积分环节。
2.2 改进ILOS导引策略本节基于自适应观测技术,设计自适应观测器来估计漂角β,再选择恰当的虚拟控制输入对漂角进行补偿,最终实现对路径的精确跟踪。基于以上思想,本文提出2种改进ILOS导引策略。
第1种改进ILOS导引策略设计如下:
(31) |
(32) |
(33) |
式中:kp=1/Δ、ki和ρ为要设计的严格大于零的常值参数。
记
由式(31)~式(33)可以看出,当USV距离计划航线较远时,首要任务是快速减小ye(t),此时积分增益
第2种改进ILOS导引策略在前者的基础上将常值增益kp用位置误差函数来替代:
(34) |
式中:λ为严格大于零的常值参数;Δmax和Δmin分别为前视距离Δ的最大值和最小值。设计如下:
(35) |
(36) |
(37) |
(38) |
时变的Δ可使得船舶操纵更加灵活[3]。由式(37)和式(38)可知,当USV距离计划航线很远时,
下面给出以上导引子系统的稳定性证明。
证明 定义
(39) |
又因为
(40) |
(41) |
将式(40)和式(41)代入式(39)中,可得
(42) |
式中:
(43) |
(44) |
显然,
(45) |
当航向控制器精确跟踪期望航向角ψd,可选择虚拟控制输入
(46) |
整理得
(47) |
设计李雅普诺夫函数V1,定义如下:
(48) |
式中:ye(t)≠0,
对其微分得
(49) |
因为
(50) |
故导引子系统在平衡点ye(t)=0处是一致全局渐近稳定的(Uniform Global Asymptotic Stability,UGAS)。
此外,存在μ>0,使得ye(t)满足D={ye(t)∈R‖ye(t)≤μ}时,
(51) |
式中:
证毕
2.3 航向、航速控制子系统设计本节采用反馈线性化比例微分控制和比例控制分别设计航向、航速控制器,实现对导引子系统给出的参考信号ψd和期望速度的跟踪[11]。控制律设计如下:
(52) |
(53) |
式中:urd为坐标变换后USV的参考相对期望速度;k1、k2、kur为严格大于零的常值增益,参考信号ψd和
证明 定义速度误差为:
(54) |
(55) |
系统式(54)和式(55)为线性时不变的,此外k1、k2、kur和
证毕
3 系统稳定性证明导引子系统与航向控制子系统一起构成级联系统,完成对期望航线的跟踪。航向控制子系统通过航向跟踪误差
证明 将位置误差动态与航向跟踪误差动态写成如下级联形式:
(56) |
(57) |
式中:
由2.2节可知,取虚拟控制输入
(58) |
(59) |
式中:
(60) |
(61) |
证毕
下面给出级联系统式(60)和式(61)的稳定性及证明过程。
定理1 系统F(t, z)在平衡点(ye(t), h)处是全局k-指数稳定的。
证明 由式(58)可知,当
(62) |
(63) |
式中:
(64) |
(65) |
由本节可知子系统式(63)在平衡点
(66) |
由于U>0,
(67) |
式中:
由式(45)可知函数ϕ1(ye(t),
(68) |
式中:β0 max为漂角的最大值。
由式(68)知|g(t, ye(t), h)|不恒为零,因此级联系统式(62)、式(63)是全局k-指数稳定的,即系统式(58)的标称系统F(t, z)是全局k-指数稳定的。
证毕
定理2 级联系统式(58)、式(59)是全局k-指数稳定的。
证明 首先证明式(59)是全局k-指数稳定的,取李雅普诺夫函数:
(69) |
对其微分得
(70) |
由定理1可知式(59)是全局k-指数稳定的,又由于|G(t, z,
证毕
最后,讨论时变参数kp和积分增益k对导引子系统稳定性的影响:根据式(36)可知,较大的前视距离Δ和积分增益k会限制导引系统的指数稳定域。在第1种改进导引策略中,只有参数k影响系统的指数稳定域,而在第2种改进ILOS导引策略中,参数Δ和k一起限制了系统的指数稳定域[13-14]。
4 仿真实验为验证本文提出的改进ILOS导引策略的有效性和优越性,以文献[15]中的模型为例,在MATLAB/Simulink中进行路径跟踪对比仿真实验。USV模型具体参数如表 1所示。
参数 | 数值 |
m11 | 25.8 |
m22 | 33.8 |
m23 | -11.748 |
m32 | -11.748 |
m33 | 6.813 |
b11 | 1 |
b22 | 0 |
b32 | 1 |
d11 | 2.0 |
d22 | 7.0 |
d23 | -2.542 5 |
d32 | -2.542 5 |
d33 | 1.422 |
通常情况下USV的期望航线分为直线和曲线2种,故将本文提出的2种改进ILOS导引新策略下和传统ILOS导引策略下的路径跟踪控制算法在直线和曲线2种工况下进行仿真实验对比。2种工况下航向、航速控制子系统的参数选择为:k1=5,k2=3,kur=5;导引子系统参数选择如下:传统的ILOS导引策略中,γ=0.005,Δ=8 m;第1种改进ILOS导引策略中,选择ki=0.001,ρ=5,Δ=8 m;第2种改进ILOS导引策略中,选择ki=0.000 9,ρ=5,λ=3,Δmin=6 m,Δmax=14 m。
1) 工况1。USV的初始位置为(x0, y0)=(50, -100)m,初始航向ψ0=π/4 rad,初始航速ur=0.5 m/s;大地坐标系下的海流设置为:Uc=0.2 m/s,ψc=30°;期望航线为y=x,期望航速ud=3 m/s。仿真结果如图 2所示。其中,传统ILOS导引策略下(以下简称为算法1)和2种改进ILOS导引策略下(以下简称为算法2、3)的路径跟踪控制算法的相对纵向速度分别对应变量下标为1、2、3的仿真结果。
2) 工况2。USV的初始位置为(x0, y0)=(0, -10) m,初始航向ψ0=π/4 rad,初始航速ur=0.5 m/s;大地坐标系下的海流设置为:Uc=0.2 m/s,ψc=30°;期望航线为由(0, 0) m、(100, 50) m、(200, 100) m、(300, 50) m、(600, 200) m经过3次样条差值生成的连续曲线,期望航速ud=3 m/s。仿真结果如图 3所示。
分析图 2(a)和图 3(a)可知,在仅有纵向控制力和转艏控制力矩、没有侧向推进器的欠驱动情况下,本文设计的3种控制算法均能够保证欠驱动USV快速地跟踪期望的直线和曲线路径,并沿该路径航行;曲线路径跟踪过程中,在曲率变化较大时,传统ILOS导引策略超调和振荡较大,而第1和第2种改进ILOS导引策略超调和振荡较小,且第2种ILOS改进导引策略最优。
从图 2(b)和图 3(b)可以看出,直线和曲线路径跟踪过程中,3种控制算法均能够保证位置误差ye最终收敛到0附近;传统ILOS导引策略下,ye的振荡大、持续时间长,第1和第2种改进ILOS导引策略下,ye的振荡小、持续时间较短。
分析图 2(c)、(f)和图 3(c)、(f)可以看出,直线和曲线路径跟踪过程中,3种控制算法均能够保证航速控制在ud=3 m/s附近,稳态误差较小,但算法1较其他2种改进算法作用下的前进推力在初始加速阶段所需推力较大。
从图 2(d)可以看出,直线路径跟踪过程中,3种算法都能很好的跟踪期望航向角,但算法1给出的期望航向角振荡较大,而算法2、3给出的期望航向角相对较为平缓,这也导致图 2(e)中算法1作用下的转艏力矩振荡幅度较大,而算法2、3则作用下的转艏力矩变化相对较为平滑。
从图 3(d)、(e)可以看出,曲线路径跟踪过程中,3种算法都能很好的跟踪期望航向角,且算法3给出的期望航向角最为平滑,同时转艏力矩也最小。
5 结论1) 本文研究了一类非对称欠驱动USV的路径跟踪问题,综合考虑了USV的欠驱动性、非对称性以及外界干扰如海流等因素,提出了2种改进ILOS导引策略和反馈控制思想的路径跟踪控制算法,实现了USV水平面内的路径跟踪。
2) 基于级联系统和李雅普诺夫理论证明了当所有控制目标实现时,控制系统是全局k-指数稳定的。
3) 通过与传统ILOS导引策略的路径跟踪算法进行仿真实验对比分析可以看出,所提出的2种改进引导策略算法具有更好的路径跟踪效果和动态性能,具有一定的先进性。
下一步将研究风、浪、流等复杂海况下以及模型不确定情况下欠驱动USV的路径跟踪控制问题。
[1] |
张树凯, 刘正江, 张显库, 等. 无人船艇的发展及展望[J].
世界海运, 2015, 38 (9): 29–36.
ZHANG S K, LIU Z J, ZHANG X K, et al. The development and outlook of unmanned vessel[J]. World Shipping, 2015, 38 (9): 29–36. (in Chinese) |
[2] |
薛春祥, 黄孝鹏, 朱咸军, 等. 外军无人系统现状与发展趋势[J].
雷达与对抗, 2016 (36): 1–5.
XUE C X, HUANG X P, ZHU X J, et al. Status quo and development trends of foreign military's unmanned systems[J]. Radar & ECM, 2016 (36): 1–5. (in Chinese) |
[3] | FOSSEN T I. Handbook of marine craft hydrodynamics and motion control[M]. Hudson County: John Wiley & Sons, 2011: 60-89. |
[4] | CAHARIJA W, PETTERSEN K Y, SØRENSEN A J, et al. Relative velocity control and integral LOS for path following of autonomous surface vessels:Merging intuition with theory[J]. Engineering for the Marine Environment, 2014, 228 (2): 180–191. |
[5] |
田勇, 王丹, 彭周华, 等. 无人水面艇直线航迹跟踪控制器的设计与验证[J].
大连海事大学学报, 2015, 41 (4): 14–18.
TIAN Y, WANG D, PENG Z H, et al. Design and validation of path tracking controller for USV along straight-lines[J]. Journal of Dalian Maritime University, 2015, 41 (4): 14–18. (in Chinese) |
[6] | FOSSEN T I, LEKKAS A M. Direct and indirect adaptive integral line-of-sight path-following controllers for marine craft exposed to ocean currents[J]. International Journal of Adaptive Control and Signal Processing, 2015, 28 (3): 20–35. |
[7] | BØRHAUG E, PAVLOV A, PETTERSEN K Y. Integral LOS control for path following of underactuated marine surface vessels in the presence of constant ocean currents[C]//Proceedings of the 47th IEEE Conference on Decision and Control. Piscataway, NJ: IEEE Press, 2008: 4984-4991. http://ieeexplore.ieee.org/xpls/icp.jsp?arnumber=4739352 |
[8] | WANG H, WANG D, PENG Z H. Adaptive neural control for cooperative path following of marine surface vehicles:State and output feedback[J]. International Journal of Systems Science, 2016, 47 (2): 343–359. DOI:10.1080/00207721.2015.1056274 |
[9] | FREDRIKSEN E, PETTERSEN K Y. Global-exponential waypoint maneuvering of ships:Theory and experiments[J]. Automatica, 2006, 42 (4): 677–687. DOI:10.1016/j.automatica.2005.12.020 |
[10] | FOSSEN T I, PETTERSEN K Y, GALEAZZI R. Line-of-sight path following for Dubins paths with adaptive sideslip compensation of drift forces[J]. IEEE Transactions on Control Systems Technology, 2015, 23 (2): 820–827. DOI:10.1109/TCST.2014.2338354 |
[11] | LEKKAS A M, FOSSEN T I. Integral LOS path following for curved paths based on a monotone cubic Hermite spline parametrization[J]. IEEE Transactions on Control Systems Technology, 2014, 22 (6): 2287–2301. DOI:10.1109/TCST.2014.2306774 |
[12] | ZHENG Z W, HUO W, WU Z. Direct-adaptive fuzzy path following control for an autonomous airship[J]. Control and Decision, 2014, 29 (3): 1418–1424. |
[13] | WⅡG M S, PETTERSEN K Y, KROGSTAD T R. Uniform semiglobal exponential stability of integral line-of-sight guidance laws[J]. IFAC Papers Online, 2015, 48 (16): 61–68. DOI:10.1016/j.ifacol.2015.10.259 |
[14] | FOSSEN T I, PETTERSEN K Y. On uniform semiglobal exponential stability (USGES) of proportional line-of-sight guidance laws[J]. Automatica, 2014, 50 (11): 2912–2917. DOI:10.1016/j.automatica.2014.10.018 |
[15] | BREIVIK M. Nonlinear maneuvering control of underactuated ships[D]. Trondheim: Norwegian University of Science and Technology, 2003: 10-30. https://www.researchgate.net/publication/215523100_Nonlinear_Maneuvering_Control_of_Underactuated_Ships |