摘要
针对传统无人机航迹规划算法应用在突发威胁场景下存在搜索点冗余、路径规划实时性较差等问题,提出了一种基于多因素Dubins路径的无人机动态航迹规划算法。该算法首先根据无人机自身性能约束及突发威胁区域的位置,并且考虑无人机的起始和最终位置,利用传统的Dubins路径找到有效的路径扩展点。然后结合启发式搜索思想建立基于路径长度和威胁的路径扩展点评估函数。最后通过路径评估函数计算,比较路径点的代价值,选取每一步的路径扩展点,规划出较优路径。仿真结果表明,在突发威胁场景下利用该算法进行航迹规划时路径长度较短、路径扩展点较少,并且符合无人机实际飞行过程中航向角变化,可有效保障无人机的安全性和航迹规划的实时性。
无人机航迹规
在突发威胁环境下,无人机需要实时躲避突发威胁区域。目前,研究人员已经提出了多种航迹规划算法,然而文献[
综上,面对突发威胁情况,上述传统航迹规划算法都存在一些局限性,不能应用于真实场景中无人机的航迹规划。因此,本文提出一种适合突发威胁场景的多因素Dubins算法(Multiple factors Dubins algorithm, MFDA)。该算法通过将无人机所处位置与突发威胁区域圆作Dubins曲线,并建立路径扩展点评估函数进行选择,不需要对威胁外的区域进行大量计算和重复搜索,大大地降低了路径搜索点的冗余,同时也节省了路径规划时间。论文的主要创新点和贡献如下:(1)突发威胁场景中,在满足曲率约束和规定的起始点和终点的切线方向的条件下引入Dubins路径规划出适合无人机真实环境下的航迹曲线。(2)建立扩展点评估函数,引入路径长度代价及威胁代价作为路径评估因素,提高路径搜索精度。(3)以Dubins路径与威胁区域相切的点作为待选路径扩展点,在待选扩展点集合中选择最佳扩展点从而降低搜索点冗余性。
本文提出的MFDA算法是基于Dubins路径实现的,因此本节对Dubins路径的类型进行了介绍,并在此基础上对路径规划过程中路径扩展点的选择模型进行了分析。
Dubins路径是在满足曲率约束和规定的起始点和终点的切线方向的条件下,连接两个二维平面(即X⁃Y平面)的最短路

图1 Dubins最短路径
Fig.1 The shortest path of Dubins
本文中节点表示为(xi,yi,ψi),其中xi表示无人机在X平面映射的横坐标,yi表示在Y平面映射的纵坐标,ψi表示无人机的航向。通过对上述4种Dubins路径的分析,可以发现,在航迹规划过程中,除了考虑无人机的起始和目标位置信息,在添加起始位置航向以及目标位置航向的基础上,当无人机跨越威胁区域飞行时,Dubins路径无疑是最短路
本文提出的MFDA算法将直线航迹与威胁区域圆的切点作为算法的下一步扩展节点,如

图2 扩展点选择
Fig.2 Extension point selection
(1) |
式中:函数g(n)定义为从已知开始节点到任意节点n的一条最优路径的代价,函数h(n)定义为从节点n到目标节点的整个目标节点集合内的最小代价路径的代价。节点n到目标节点的随机路径为h(n),即节点n到目标节点的最佳路径。
人工势场法在1986年由Khatib提出,主要思想是通过建立目标点的引力场以及障碍物的斥力场,在两个势力场合力的作用下沿着势场函数下降的方向移动,由此形成一条避开障碍物的最优路径。由

图3 人工势场法基本原理
Fig.3 Basic principles of artificial potential field method
在人工势场法中合力即为势能函数的梯度,联合合力和无人机的运动状态方程即可得到其飞行轨迹,由于人工势场法避免了搜索最优解或非线性参数优化的过程,因此路径规划的速度相比于
本文提出一种多因素Dubins路径算法,重新建立基于路径长度代价及威胁代价作为路径评估因素的评估函数,通过1.2节路径扩展点的选取可以看出,在当前路径扩展点选择下一步路径扩展点的过程中,根据Dubins曲线的特性,待选路径扩展点中最多只会存在两个选择,并且Dubins路径只作与威胁区域相切的曲线,那么在从起始点到目标点的整个路径规划过程中路径扩展点的数量会大幅度降低,另外Dubins曲线是在满足曲率约束和规定的起始点和终点的切线方向的条件下,连接两个二维平面(即X⁃Y平面)的最短路径,那么最终通过MFDA算法得到的规划路径长度也会减小。MFDA算法路径规划过程如

图4 MFDA算法路径规划过程
Fig.4 Path planning process of MFDA
算法路径规划过程步骤如下:
(1)首先将起始点当作当前路径扩展点,生成到目标点的路径。如果在此方向上不存在威胁区域,则将目标点Pf加入到CLOSE表中,构建成一条从起始点到目标点的路径曲线,算法运行结束。如果在此方向上存在威胁区域,则生成当成当前路径扩展点到威胁区域的多条Dubins曲线。
(2)通过MFDA算法估价函数计算并选出代价值最小的点作为下一步的路径扩展点,将此路径扩展点加入到CLOSE表中。
(3)将步骤(2)中选取的路径扩展点作为当前的位置,然后重新执行步骤(1)的操作。之后重复上述过程,不断产生新的路径扩展点,直到到达目的点,构建出完整路径。
由上述MFDA算法路径规划过程可知,路径搜索的过程也就是通过计算不同Dubins路径的代价值从而选取路径扩展点的过程,因此下节着重分析步骤(2)中路径扩展点的代价计算和选取。
MFDA算法同时考虑路径代价和威胁代价对无人机轨迹的影响,估价函数建立如下
(2) |
式中:g(n)=w1·Distance+w2·Hazard,w1表示距离代价的影响因子,w2表示威胁代价的影响因子,w1+w2=1。Distance距离由两端弧长和一段直线长度构成,由

图5 LSR路径类型
Fig.5 Path type of LSR
起始节点(xi,yi,ψi),求解圆Oi和圆Othreat的切线QiPn的坐标。
步骤1 计算圆Oi的圆心坐标
(3) |
(4) |
威胁区域Othreat的圆心坐标已知。
步骤2 生成直线航迹QiPn
(5) |
(6) |
(7) |
因此,逆时针旋转β1角度与重合,的模L已知,之后按照比例进行缩放可以得到,从而计算出切点Pn坐标,同理可得到切点Qi坐标。
当无人机在位置为Pi,航向为ψi准备绕过威胁区域时,选择切点Pn或者作为
(8) |
式中节点Qi、Pn的位置已有前述步骤计算可知,因此和OiPn的长度也可计算得到。接下来计算Hazard,距离威胁区域越近,所受到的威胁就越大,假设把不经过威胁区域的航迹路径所受到的威胁设置为0,因此只有绕威胁区域飞行的圆弧段航迹才会受到威胁。从当前节点到下一个扩展节点的威胁代价值为
(9) |
因此,,其中
(10) |
因此
(11) |
计算h(n),同样将距离代价和威胁代价的预评估影响因子分别用w1、w2表示。当无人机以Pn或点切入威胁区域,则无人机分别以PnQnPf路径或路径以指定位姿到达目标点,则和分别预示了无人机绕过威胁区域的飞行轨迹,所以h(Pn)表示为
(12) |
式中
(13) |
(14) |
通过上述代价值的计算与比较,选择代价值较小的节点作为MADA算法的下一步路径扩展点。
因此,根据2.3节和2.4节的MFDA算法的路径扩展点选取过程,通过不断的循环计算,可以规划出一条在突发威胁场景下从起始点到目标点的无人机飞行航迹。
根据前述分析,对本文设计的航迹规划算法进行了仿真验证,并与传统的
图

图6 单威胁区域场景下航迹规划结果
Fig.6 Track planning results under single threat area scenario

图7 双威胁区域场景下航迹规划结果
Fig.7 Track planning results under two threat zone scenarios
在此场景下,MFDA算法所得路径由2段Dubins曲线组成,长度为270.914 6 m,并且路径光滑。在相同条件下,由
区域中共有6个威胁区域,其位置、半径和威胁等级表示为:,,,,,。无人机起始位姿和最终位姿分别表示为:,。

图8 多个威胁区域场景下航迹规划结果
Fig.8 Track planning results under multiple threat area scenarios
分别对上述不同场景下航迹规划算法的性能指标进行统计并分析。图

图9 3个场景下路径规划长度
Fig.9 Path planning length under three scenarios

图10 3个场景下路径扩展点数量
Fig.10 Number of path extension points under three scenarios

图11 3个场景下路径规划时间
Fig.11 Path planning time under three scenarios

图12 路径航向变化
Fig.12 Change of course heading
从
从
针对无人机飞行过程中遭遇突发威胁区域的场景,提出了一种多因素Dubins算法,该算法在Dubins路径的基础上建立路径扩展点评估函数,通过引入路径长度评估因子和威胁评估因子对路径扩展点进行选择,从而可以有效地降低路径搜索点的数量。同时结合启发式搜索的思想,对可能出现的路径长度代价和威胁代价进行评估,达到了缩短路径长度的目的。并且增加无人机起始点方向和到达目的点方向限制条件,构建了适用于真实环境中的无人机飞行航迹。仿真表明,MFDA算法在突发威胁区域场景下,可以规划出较短的路径,同时相比于传统的航迹规划算法具有较少的路径搜索点,并且得到的路径符合无人机实际飞行时的航向变化。
本文仅考虑了单架无人机的航迹规划,在下一步的研究中将研究在无人机集群情况下的航迹安全和威胁区域规避等基本问题。
参考文献
Bortoff S A. Path planning for UAVs[C]//Proceedings of the 2000 American Control Conference. [S.l.]: IEEE, 2000: 364-368. [百度学术]
Edsger W D. A note on two problems in connexion with graphs[J]. Numerische mathematik, 1959,1: 269-271. [百度学术]
LaValle S M. Planning algorithms[M]. Cambridge, UK: Cambridge University Press, 2006. [百度学术]
姚远, 周兴社, 张凯龙, 等. 基于稀疏
YAO Yuan, ZHOU Xingshe, ZHANG Kailong, et al. Dynamic trajectory planning for unmanned aerial vehicle based on sparse
Wang C, Wang L, Qin J, et al. Path planning of automated guided vehicles based on improved A-star algorithm[C]//Proceedings of IEEE International Conference on Information and Automation. [S.l.]: IEEE, 2015: 2071-2076. [百度学术]
Khatib O. Real-time obstacle avoidance for manipulators and mobile robots[M]. New York, USA: Springer, 1986: 396-404. [百度学术]
Weiguang L, Xia S. AGV path planning based on improved A* algorithm[J]. Modern Manufacturing Engineering, 2015 (10): 9. [百度学术]
Cui Shigang, Wang Hui, Yang Li. A simulation study of A-star algorithm for robot path planning[C]//Proceedings of 16th International Conference on Mechatronics Technology. [S.l.]: IEEE,2012: 506-510. [百度学术]
Wang C, Wang L, Qin J, et al. Path planning of automated guided vehicles based on improved A-Star algorithm[C]//Proceedings of 2015 IEEE International Conference on Information and Automation. [S.l.]: IEEE, 2015: 2071-2076. [百度学术]
许凯波, 鲁海燕, 黄洋, 等. 基于双层蚁群算法和动态环境的机器人路径规划方法[J]. 电子学报, 2019, 47(10): 2166-2176. [百度学术]
XU Kaibo, LU Haiyan, HUANG Yang, et al. Robot path planning based on double-layer ant colony optimization algorithm and dynamic environment[J]. Acta Electronica Sinica, 2019, 47(10): 2166-2176. [百度学术]
Bounini F, Gingras D, Pollart H, et al. Modified artificial potential field method for online path planning applications[C]//Proceedings of 2017 IEEE Intelligent Vehicles Symposium (IV). [S.l.]: IEEE, 2017: 180-185. [百度学术]
Wu L, Zha H, Xiu C, et al. Local path planning for intelligent vehicle obstacle avoidance based on dubins curve and tentacle algorithm[C]//Intelligent & Connected Vehicles Symposium. [S. l.]: IEEE, 2017: 148-191. [百度学术]
Tuqan M, Daher N, Shammas E. A simplified path planning algorithm for surveillance missions of unmanned aerial vehicles[C]//Proceedings of 2019 IEEE/ASME International Conference on Advanced Intelligent Mechatronics. [S.l.]: IEEE, 2019: 1341-1346. [百度学术]
韩庆田. 基于 Dubins 曲线的多 UAV 协同航迹规划研究[J]. 航空计算技术, 2019, 49(4): 13-16. [百度学术]
HAN Qingtian. Research on cooperate path planning of multiple UAVs based on Dubins curve[J]. Aeronautical Computing Technique, 2019, 49(4): 13-16. [百度学术]
Živojević D, Velagić J. Path planning for mobile robot using Dubins-curve based RRT algorithm with differential constraints[C]//Proceedings of 2019 International Symposium ELMAR. [S.l.]: IEEE, 2019: 139-142. [百度学术]
Guimaraes M D, Alves N A, FIUZA D C N, et al. An evolutionary approach for the dubins’ traveling salesman problem with neighborhoods[C]//Proceedings of the 14th Annual Conference on Genetic and Evolutionary Computation. New York: [s.n.], 2012: 377-384. [百度学术]
Duchoň F, Babinec A, Kajan M, et al. Path planning with modified A-star algorithm for a mobile robot[J]. Procedia Engineering, 2014, 96: 59-69. [百度学术]
Chen Y, Luo G, Mei Y, et al. UAV path planning using artificial potential field method updated by optimal control theory[J]. International Journal of Systems Science, 2016, 47(6): 1407-1420. [百度学术]