CN111781948A - 一种无人机编队队形变换与在线动态避障方法 - Google Patents

一种无人机编队队形变换与在线动态避障方法 Download PDF

Info

Publication number
CN111781948A
CN111781948A CN202010559970.0A CN202010559970A CN111781948A CN 111781948 A CN111781948 A CN 111781948A CN 202010559970 A CN202010559970 A CN 202010559970A CN 111781948 A CN111781948 A CN 111781948A
Authority
CN
China
Prior art keywords
unmanned aerial
aerial vehicle
repulsion
obstacle
force
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.)
Withdrawn
Application number
CN202010559970.0A
Other languages
English (en)
Inventor
郑宇�
曾旭
金滢
赵雪冬
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing Non Air Aviation Technology Co ltd
Original Assignee
Nanjing Non Air Aviation Technology Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing Non Air Aviation Technology Co ltd filed Critical Nanjing Non Air Aviation Technology Co ltd
Priority to CN202010559970.0A priority Critical patent/CN111781948A/zh
Publication of CN111781948A publication Critical patent/CN111781948A/zh
Withdrawn legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/10Simultaneous control of position or course in three dimensions
    • G05D1/101Simultaneous control of position or course in three dimensions specially adapted for aircraft
    • G05D1/104Simultaneous control of position or course in three dimensions specially adapted for aircraft involving a plurality of aircrafts, e.g. formation flying

Landscapes

  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

本发明公开了一种无人机编队队形变换与在线动态避障方法,属于无人机的技术领域。具体包括以下步骤:步骤一、通过人工势场算法计算出无人机斥力场、斥力及无人机受到的合力;步骤二、对步骤一种的人工势场算法计算进行改进;步骤三、结合改进的人工势场算法计算,提出了一种基于Kernel方法的期望路径平滑策略用于解决无人机在避障过程中频繁进出障碍物膨胀区域,造成存在期望路径折线过多、拐点角度大的问题;步骤四、设计机间避障算法,用于实现无人机相互之间的避让。本发明根据固定翼无人机的飞行特性,设计了机间避障人工势场,使得多机在编队飞行时,即不会飞入静态障碍物区域也不会和其他飞机发生碰撞。

Description

一种无人机编队队形变换与在线动态避障方法
技术领域
本发明属于无人机的技术领域,特别是涉及一种无人机编队队形变换与在线动态避障方法。
背景技术
无人机编队飞行时,编队内部各机之间必须保持一定的安全距离,因此在线动态是必不可少的。目前成熟的避障与路径规划算法有人工势场法、蚁群算法、快速扩展随机数法、Voronoi图法和A*算法等。这些方法根据对周围环境的感知情况,可以分为全局规划和局部规划。全局规划已知全局信息,局部规划对全局信息未知,只能结合自身传感器实现对自身局部信息的感知。动态环境下的无人机编队避障以及避障后路径重规划是集群技术领域的一个热点和难点,运动情况下的无人机之间的碰撞也会造成很大的损失,因此实时避障对于飞行中的编队飞机来说尤其重要。
发明内容
本发明为解决上述背景技术中存在的技术问题,提出了一种无人机编队队形变换与在线动态避障方法。
本发明采用以下技术方案来实现:一种无人机编队队形变换与在线动态避障方法,具体包括以下步骤:
步骤一、通过人工势场算法计算出无人机斥力场、斥力及无人机受到的合力;
步骤二、对步骤一种的人工势场算法计算进行改进;
步骤三、结合改进的人工势场算法计算,提出了一种基于Kernel方法的期望路径平滑策略用于解决无人机在避障过程中频繁进出障碍物膨胀区域,造成存在期望路径折线过多、拐点角度大的问题;
步骤四、设计机间避障算法,用于实现无人机相互之间的避让;
其中,所述步骤一具体包括:
定义无人机的位置为q,则无人机在q点的引力场函数Uatt(q)和引力Fatt(q)分别满足以下公式:
Figure BDA0002545736250000011
Figure BDA0002545736250000021
式中,λ表示引力增益,ρα(q)表示q点到目标a之间的距离;
无人机在q点的斥力场和斥力为:
Figure BDA0002545736250000022
Figure BDA0002545736250000023
η为斥力增益,ρ0为障碍物的最大影响半径,
ρr(q)为q点到障碍物r的距离;
无人机所受到的合力F为:
F=Fatt+Frep
Fatt为目标对无人机的引力,Frep为障碍物对无人机的斥力,定义障碍物对无人机的斥力和目标点对无人机的引力分别为对应斥力场和引力场的负梯度,则q点受到目标点的引力Fatt(q)和障碍物的斥力Frep(q)分别为该点引力势函数和斥力势函数的负梯度,
假设同时具有N个障碍物对无人机影响是,根据叠加原理,无人机的合力F为:
Figure BDA0002545736250000024
式中
Figure BDA0002545736250000025
为第i个障碍物产生的斥力。
在进一步的实施例中,当判断出无人机在某点引力大小不为零,并且无人机所受到的总合力为零:则可以判断此时无人机陷入局部最优解中,假设存在一个等效的障碍物,该障碍物产生的斥力和无人机受到的斥力大小方向均相同,此时无人机、等效障碍物、目标必然同线;等效障碍物的具***置由人工势场法反推出,规定无人机跳出局部最优的飞行方向为无人机位置与障碍物膨胀区圆的切线方向,左右任取。
在进一步的实施例中,当由于引力和斥力的大小不成比例,此时的引力场函数为:
Figure BDA0002545736250000031
式中,
Figure BDA0002545736250000032
代表无人机无目标之间的距离阈值,ξ为斥力系数,
相对应的引力计算方法为:
Figure BDA0002545736250000033
在进一步的实施例中,
当目标点处于不可达情况,在人工势场法的基础上叠加无人机和目标点距离的n次幂,则新的斥力场函数计算方法为:
Figure BDA0002545736250000034
在进一步的实施例中,
为了增加平滑的效果,定义y*(i)为平滑后的数据,带宽为h,S(t)表示加权函数,窗口长度为n,y(j)表示窗口内的第j个数据,则在t时刻拟合值的计算方法为:
Figure BDA0002545736250000035
式中加权函数Sj(t)的计算方法:
Figure BDA0002545736250000036
式中
Figure BDA0002545736250000037
tj表示当前时刻的时间,kern(x)就是Kernel函数,是Kernel平滑的核心。
在进一步的实施例中,所述步骤四包括:假设编队中第j架飞机对环境中的位置q产生斥力作用,则其对q点斥力场函数定义如下:
Figure BDA0002545736250000041
式中,ξ为斥力系数,σ0为无人机膨胀区半径,
Figure BDA0002545736250000042
代表无人机和q点之间的距离,ε∈(0,1)为膨胀区分离系数;
对应的该无人机对q的斥力为:
Figure BDA0002545736250000043
当机间的距离大于σ0时,两机的斥力为零,彼此无需进行避障;当机间的距离小于σ0时,两机进入各自的膨胀区,彼此作为对方的障碍物,这时斥力分为两段,当
Figure BDA0002545736250000044
时,机间产生的斥力和之间距离为线性关系,斥力变化幅度较小,这相当于给无人机执行避障动作的一段缓冲时间,让彼此缓慢的远离,不至于以频繁进行大角度机动;当
Figure BDA0002545736250000045
时,此时两机距离非常近,为避免发生碰撞需要立刻进行避障,这是对着两机距离的接近,斥力成指数式增加,迫使无人机快速远离对方,在
Figure BDA0002545736250000046
点,两段范围计算出的斥力是相等的,即斥力是连续的,不会造成指令的突变。
在进一步的实施例中,无人机膨胀区半径σ0的计算步骤如下:
第一步:相邻两无人机以自身坐标为原点向飞行速度方向射线,判断两条射线之间的距离是否小于一个小距离εmin,如果不满足,说明没有碰撞的可能性,则令σ0=0,表示对方无人机对己方不产生斥力;如果两条射线之间的距离小于εmin,则两机路径可能相交,进入第二步;选取εmin的原则是依据无人机传感器误差、机体大小和控制精度等因素;
第二步:如两无人机相遇的情况,定义两机碰撞的时间阈值为Tm,该值表示在预计碰撞时,最少需要留出多少时间处理,用q1和q2代表无人机一和无人机二的位置,
Figure BDA0002545736250000051
代表无人机一和无人机二的速度矢量,δvq表示矢量
Figure BDA0002545736250000052
Figure BDA0002545736250000053
的夹角,分析无人机二,此时无人机一为作为障碍物,则σ0的计算方法如下:
Figure BDA0002545736250000054
本发明的有益效果:本发明以人工势场法为基础,针对人工势场法的三个缺点提出了相应的改进措施,然后结合Kernel方法对期望的路径进行平滑处理,解决了期望路径震荡的问题。最后在此基础上,本发明根据固定翼无人机的飞行特性,设计了机间避障人工势场,使得多机在编队飞行时,即不会飞入静态障碍物区域也不会和其他飞机发生碰撞。
附图说明
图1为目标和障碍物对无人机的引力和斥力图。
图2避免局部最优解运动方向示意图。
图3为对等效障碍物施加随机扰动图。
图4为Kernel路径平滑示意图。
图5为两架无人机的各自斥力场示意图。
图6为两无人机相遇情况图。
图7为两无人机间避障判断计算图。
图8为两无人机间斥力方向示意图。
具体实时方式
下面结合附图说明和具体实施例对本发明做进一步的描述。
首先,发明人认为:无人机编队飞行时,编队内部各机之间必须保持一定的安全距离,因此在线动态是必不可少的。目前成熟的避障与路径规划算法有人工势场法、蚁群算法、快速扩展随机数法、Voronoi图法和A*算法等。这些方法根据对周围环境的感知情况,可以分为全局规划和局部规划。全局规划已知全局信息,局部规划对全局信息未知,只能结合自身传感器实现对自身局部信息的感知。对于本发明的在线避障问题,动态障碍物有编队内自身飞行的其他飞机,静态障碍物有禁飞区、电子地理围栏等,因此本发明的避障是在全局信息已知的情况下进行。考虑到固定翼具有最小失速速度和最小转弯半径的约束特点,本发明提出了改进人工势场法、Kernel路径平滑和引入机间人工势场相结合的在线动态避障算法。
一种无人机编队队形变换与在线动态避障方法,具体包括以下步骤:步骤一、通过人工势场算法计算出无人机斥力场、斥力及无人机受到的合力;
步骤二、对步骤一种的人工势场算法计算进行改进;
步骤三、结合改进的人工势场算法计算,提出了一种基于Kernel方法的期望路径平滑策略用于解决无人机在避障过程中频繁进出障碍物膨胀区域,造成存在期望路径折线过多、拐点角度大的问题;
步骤四、设计机间避障算法,用于实现无人机相互之间的避让。
人工势场法(APF,Artificial Potential Field)是由Khatib于上世纪80年代提出的一种路径规划算法。该算法将目标和障碍物分别看做对无人机有引力和斥力的物体,无人机沿引力与斥力的合力来进行运动。图1给出了目标和障碍物对无人机的引力和斥力图。
定义无人机的位置为q,则无人机在q点的引力场函数Uatt(q)和引力Fatt(q)分别满足以下公式:
Figure BDA0002545736250000061
Figure BDA0002545736250000062
式中,λ表示引力增益,ρα(q)表示q点到目标a之间的距离;
无人机在q点的斥力场和斥力为:
Figure BDA0002545736250000063
Figure BDA0002545736250000064
η为斥力增益,ρ0为障碍物的最大影响半径,
ρr(q)为q点到障碍物r的距离;
无人机所受到的合力F为:
F=Fatt+Frep
Fatt为目标对无人机的引力,Frep为障碍物对无人机的斥力,定义障碍物对无人机的斥力和目标点对无人机的引力分别为对应斥力场和引力场的负梯度,则q点受到目标点的引力Fatt(q)和障碍物的斥力Frep(q)分别为该点引力势函数和斥力势函数的负梯度,
假设同时具有N个障碍物对无人机影响是,根据叠加原理,无人机的合力F为:
Figure BDA0002545736250000071
式中
Figure BDA0002545736250000072
为第i个障碍物产生的斥力。
步骤二是对上述的人工势场法进行改进,主要是因为人工势场法存在以下缺点:
(1)在目标某个点,如果在某个点引力和斥力刚好相等且方向相反,此时无人机容易陷入局部最优解,表现形式为无人机在该点附近持续震荡。
(2)当目标点距离较远的话,引力将变得特别大,相对较小的斥力下,无人机路径可能会碰到障碍物。
(3)当目标点附近有障碍物时,斥力将非常大,引力相对较小,无人机很难到达目标点。
针对APF算法的三个问题,本节提出以下措施,下面分别介绍。
对于问题(1)本小节提出如下解决方案:
当判断出无人机在某点引力大小不为零,并且无人机所受到的总合力为零。则可以判断此时无人机陷入局部最优解中。假设存在一个等效的障碍物,该障碍物产生的斥力和无人机受到的斥力大小方向均相同。此时无人机、等效障碍物、目标必然同线。等效障碍物的具***置由上述公式反推出。规定无人机跳出局部最优的飞行方向为无人机位置与障碍物膨胀区圆的切线方向,左右任取。图2给出了其运动方向的示意图。
避免陷入局部最优解的方法还可以通过对无人机增加额外扰动来实现。但直接对无人机增加扰动不利于对无人机的稳定飞行,本节决定在增加等效障碍物的基础上,对于等效的障碍物施加随机扰动,具体方法为在反解等效障碍物位置时,增加一个随机位置误差ε(q),如图3。
对于问题(2)可能碰到障碍物的情况,主要是由于引力和斥力的大小不成比例造成,由上述公式可知,引力与距离成线性关系,可以通过修改这种线性对应关系来避免距离过大导致引力过大。将引力场函数修改为:
Figure BDA0002545736250000081
式中,
Figure BDA0002545736250000082
代表无人机无目标之间的距离阈值,ξ为斥力系数,
相对应的引力计算方法为:
Figure BDA0002545736250000083
该公式与步骤一中的公式相比,当无人机和目标的距离小于设定的阈值时,两者计算结果一样,当大于阈值时,无人机所受到的引力不在和距离成线性关系,距离越远,引力增加的速度就越慢。该对引力进行了限制,使其不会过大而掩盖掉斥力的影响。
对于问题(3)目标点不可达情况,可以在步骤一中的公式的基础上叠加上无人机和目标点距离的n次幂。则新的斥力场函数计算方法为:
Figure BDA0002545736250000084
从上述公式中可以看出,对于目标点靠近障碍物的情况,当无人机目标点的过程中,无人机又在接近障碍物,根据步骤一中的Uatt(q)计算公式得到的斥力场在增大,由于改进后的Uatt(q)公式加入了无人机和目标点距离的影响,式中
Figure BDA0002545736250000091
在减小,这个因子可以对于斥力场的增长起到拖拽作用,从而让无人机可以顺利到达目标点。
用改进的APF算法规划出的路径,由于无人机在避障过程中频繁进出障碍物膨胀区域,造成存在期望路径折线过多、拐点角度大等问题,这些缺点增大了无人机编队飞行的风险,降低了编队执行任务的效率。本发明提出了一种基于Kernel方法的期望路径平滑策略,使得经过Kernel方法处理的期望规划路径更加适合固定翼无人机的飞行。
Kernel数据平滑方法在信息领域内广泛应用,其本质原理也是移动窗口的最小二项式平滑,利用窗口内各个数据的加权线性组合来拟合当前点,从而得到平滑效果。
定义y*(i)为平滑后的数据,带宽为h,S(t)表示加权函数,窗口长度为n,y(j)表示窗口内的第j个数据,则在t时刻拟合值的计算方法为:
Figure BDA0002545736250000092
式中加权函数Sj(t)的计算方法:
Figure BDA0002545736250000093
式中
Figure BDA0002545736250000094
tj表示当前时刻的时间,kern(x)就是Kernel函数,是Kernel平滑的核心。
其有三种表达式,代表不同的加权函数:
均匀函数:
Figure BDA0002545736250000095
二次函数:
Figure BDA0002545736250000096
高斯函数:
Figure BDA0002545736250000097
窗口长度n的选取过大,会引起加权函数Sj(t)计算上的困难,造成延时严重,n选取过小会造成欠拟合。因此本发明决定选取合适的大小的n对期望路径分段平滑,并保持各段连接点的连续性和光滑性,加权函数选择二次函数进行拟合。
假设改进APF算法在t1,t2,t3,t4,t5,t6,t7,t8时刻规划的期望点位置为P1,P2,P3,P4,P5,P6,P7,P8将P1,P2,P3,P4作为第一段,P5,P6,P7,P8作为第二段,此时n=4。图4给出了经过Kernel平滑之后的路径生成示意图。
对于本发明编队内飞行的无人机而言,都把其他飞机视作为移动障碍物。编队中每架无人机都产生斥力场,斥力作用到在膨胀区内的所有其他无人机上。编队内无人机的机间避障可以理解为一个相互避让的过程,下面以两架无人机在水平面之间的避障为例分析。执行避障过程的这架无人机,同样被另一架无人机视为障碍物。图5给出了两架无人机的各自斥力场示意图。
考虑到固定翼无人机的飞行速度远大于机器人和多旋翼,且具有最小转弯半径的限制,其避障不同于六自由度运动的机器人或者旋翼,执行避障动作需要一段时间缓冲。针对这个特点,本发明将重新设计新的斥力场函数。
假设编队中第j架飞机对环境中的位置q产生斥力作用。则其对q点斥力场函数定义如下:
Figure BDA0002545736250000101
式中,ξ为斥力系数,σ0为无人机膨胀区半径,
Figure BDA0002545736250000102
代表无人机和q点之间的距离,ε∈(0,1)为膨胀区分离系数。
对应的该无人机对q的斥力为:
Figure BDA0002545736250000103
当机间的距离大于σ0时,两机的斥力为零,彼此无需进行避障;当机间的距离小于σ0时,两机进入各自的膨胀区,彼此作为对方的障碍物,这时斥力分为两段,当
Figure BDA0002545736250000111
时,机间产生的斥力和之间距离为线性关系,斥力变化幅度较小,这相当于给无人机执行避障动作的一段缓冲时间,让彼此缓慢的远离,不至于以频繁进行大角度机动;当
Figure BDA0002545736250000112
时,此时两机距离非常近,为避免发生碰撞需要立刻进行避障,这是对着两机距离的接近,斥力成指数式增加,迫使无人机快速远离对方,在
Figure BDA0002545736250000113
点,两段范围计算出的斥力是相等的,即斥力是连续的,不会造成指令的突变。
传统APF算法判断两机是否发生斥力作用,仅用两机之间的距离作为判断。这显然不符合固定翼无人机直飞型的特点,即固定翼具有最小的失速速度,无法倒退。假设编队中所有无人机的膨胀区半径σ0都相等,现分析如图6中两机情况,图中箭头方向代表无人机的速度方向。
定义图中两机距离为假设
Figure BDA0002545736250000114
假设图中两机距离
Figure BDA0002545736250000115
相等且小于障碍物膨胀区半径σ0。对于图6中(a)图,两机距离虽然进入避障区域,但两机本身正在往反方向飞行,理论上不可能发生碰撞的可能;(b)图两机飞行的方向成90度,正在远离对方,理论上也不可能发生碰撞;(c)图两机飞行路径存在交叉,有发生碰撞的可能性;(d)图两机相对飞行,碰撞的可能性最高。
根据固定翼飞行的特性,需要重新对障碍物膨胀区半径σ0进行定义。计算步骤如下:
第一步:相邻两无人机以自身坐标为原点向飞行速度方向射线,判断两条射线之间的距离是否小于一个小距离εmin,如果不满足,说明没有碰撞的可能性,则令σ0=0,表示对方无人机对己方不产生斥力;如果两条射线之间的距离小于εmin,则两机路径可能相交,进入第二步;选取εmin的原则是依据无人机传感器误差、机体大小和控制精度等因素;
第二步:如图7同两无人机相遇的情况,定义两机碰撞的时间阈值为Tm,该值表示在预计碰撞时,最少需要留出多少时间处理,用q1和q2代表无人机一和无人机二的位置,
Figure BDA0002545736250000116
代表无人机一和无人机二的速度矢量,δvq表示矢量
Figure BDA0002545736250000117
Figure BDA0002545736250000121
的夹角,分析无人机二,此时无人机一为作为障碍物,则σ0的计算方法如下:
Figure BDA0002545736250000122
考虑在编队飞行中,长机的位置和航向作为其他各僚机的参考点,本发明决定在避障过程中,长机只产生斥力场,不执行避障执行。当与其他僚机相遇,由其他僚机执行避障指令。本发明将障碍物运动速度的切向方向作为斥力的方向。如图8,图中F1
Figure BDA0002545736250000123
代表无人机一所受的合力、引力和斥力,F2
Figure BDA0002545736250000124
代表无人机二所受的合力、引力和斥力。
对于编队的在线动态避障问题,本发明以人工势场法为基础,针对人工势场法的三个缺点提出了相应的改进措施,然后结合Kernel方法对期望的路径进行平滑处理,解决了期望路径震荡的问题。最后在此基础上,本发明根据固定翼无人机的飞行特性,设计了机间避障人工势场,使得多机在编队飞行时,即不会飞入静态障碍物区域也不会和其他飞机发生碰撞。

Claims (7)

1.一种无人机编队队形变换与在线动态避障方法,其特征在于,具体包括以下步骤:
步骤一、通过人工势场算法计算出无人机斥力场、斥力及无人机受到的合力;
步骤二、对步骤一种的人工势场算法计算进行改进;
步骤三、结合改进的人工势场算法计算,提出了一种基于Kernel方法的期望路径平滑策略用于解决无人机在避障过程中频繁进出障碍物膨胀区域,造成存在期望路径折线过多、拐点角度大的问题;
步骤四、设计机间避障算法,用于实现无人机相互之间的避让;
其中,所述步骤一具体包括:
定义无人机的位置为q,则无人机在q点的引力场函数Uatt(q)和引力Fatt(q)分别满足以下公式:
Figure FDA0002545736240000011
Figure FDA0002545736240000012
式中,λ表示引力增益,ρα(q)表示q点到目标a之间的距离;
无人机在q点的斥力场和斥力为:
Figure FDA0002545736240000013
Figure FDA0002545736240000014
η为斥力增益,ρ0为障碍物的最大影响半径,
ρr(q)为q点到障碍物r的距离;
无人机所受到的合力F为:
F=Fatt+Frep
Fatt为目标对无人机的引力,Frep为障碍物对无人机的斥力,定义障碍物对无人机的斥力和目标点对无人机的引力分别为对应斥力场和引力场的负梯度,则q点受到目标点的引力Fatt(q)和障碍物的斥力Frep(q)分别为该点引力势函数和斥力势函数的负梯度,
假设同时具有N个障碍物对无人机影响是,根据叠加原理,无人机的合力F为:
Figure FDA0002545736240000021
式中
Figure FDA0002545736240000022
为第i个障碍物产生的斥力。
2.根据权利要求1所述的一种无人机编队队形变换与在线动态避障方法,其特征在于,当判断出无人机在某点引力大小不为零,并且无人机所受到的总合力为零:则可以判断此时无人机陷入局部最优解中,假设存在一个等效的障碍物,该障碍物产生的斥力和无人机受到的斥力大小方向均相同,此时无人机、等效障碍物、目标必然同线;等效障碍物的具***置由人工势场法反推出,规定无人机跳出局部最优的飞行方向为无人机位置与障碍物膨胀区圆的切线方向,左右任取。
3.根据权利要求1所述的一种无人机编队队形变换与在线动态避障方法,其特征在于,当由于引力和斥力的大小不成比例,此时的引力场函数为:
Figure FDA0002545736240000023
式中,
Figure FDA0002545736240000024
代表无人机无目标之间的距离阈值,ξ为斥力系数,
相对应的引力计算方法为:
Figure FDA0002545736240000025
4.根据权利要求1所述的一种无人机编队队形变换与在线动态避障方法,其特征在于,当目标点处于不可达情况,在人工势场法的基础上叠加无人机和目标点距离的n次幂,则新的斥力场函数计算方法为:
Figure FDA0002545736240000031
5.根据权利要求1所述的一种无人机编队队形变换与在线动态避障方法,其特征在于,为了增加平滑的效果,定义y*(i)为平滑后的数据,带宽为h,S(t)表示加权函数,窗口长度为n,y(j)表示窗口内的第j个数据,则在t时刻拟合值的计算方法为:
Figure FDA0002545736240000032
式中加权函数Sj(t)的计算方法:
Figure FDA0002545736240000033
式中
Figure FDA0002545736240000034
tj表示当前时刻的时间,kern(x)就是Kernel函数,是Kernel平滑的核心。
6.根据权利要求1所述的一种基于分布式架构的无人机编队协同控制***的控制方法,其特征在于,所述步骤四包括:假设编队中第j架飞机对环境中的位置q产生斥力作用,则其对q点斥力场函数定义如下:
Figure FDA0002545736240000035
式中,ξ为斥力系数,σ0为无人机膨胀区半径,
Figure FDA0002545736240000036
代表无人机和q点之间的距离,ε∈(0,1)为膨胀区分离系数;
对应的该无人机对q的斥力为:
Figure FDA0002545736240000041
当机间的距离大于σ0时,两机的斥力为零,彼此无需进行避障;当机间的距离小于σ0时,两机进入各自的膨胀区,彼此作为对方的障碍物,这时斥力分为两段,当
Figure FDA0002545736240000042
时,机间产生的斥力和之间距离为线性关系,斥力变化幅度较小,这相当于给无人机执行避障动作的一段缓冲时间,让彼此缓慢的远离,不至于以频繁进行大角度机动;当
Figure FDA0002545736240000043
时,此时两机距离非常近,为避免发生碰撞需要立刻进行避障,这是对着两机距离的接近,斥力成指数式增加,迫使无人机快速远离对方,在
Figure FDA0002545736240000044
点,两段范围计算出的斥力是相等的,即斥力是连续的,不会造成指令的突变。
7.根据权利要求6所述的一种基于分布式架构的无人机编队协同控制***的控制方法,其特征在于,无人机膨胀区半径σ0的计算步骤如下:
第一步:相邻两无人机以自身坐标为原点向飞行速度方向射线,判断两条射线之间的距离是否小于一个小距离εmin,如果不满足,说明没有碰撞的可能性,则令σ0=0,表示对方无人机对己方不产生斥力;如果两条射线之间的距离小于εmin,则两机路径可能相交,进入第二步;选取εmin的原则是依据无人机传感器误差、机体大小和控制精度等因素;
第二步:如两无人机相遇的情况,定义两机碰撞的时间阈值为Tm,该值表示在预计碰撞时,最少需要留出多少时间处理,用q1和q2代表无人机一和无人机二的位置,
Figure FDA0002545736240000045
代表无人机一和无人机二的速度矢量,δvq表示矢量
Figure FDA0002545736240000046
Figure FDA0002545736240000047
的夹角,分析无人机二,此时无人机一为作为障碍物,则σ0的计算方法如下:
Figure FDA0002545736240000048
CN202010559970.0A 2020-06-18 2020-06-18 一种无人机编队队形变换与在线动态避障方法 Withdrawn CN111781948A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010559970.0A CN111781948A (zh) 2020-06-18 2020-06-18 一种无人机编队队形变换与在线动态避障方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010559970.0A CN111781948A (zh) 2020-06-18 2020-06-18 一种无人机编队队形变换与在线动态避障方法

Publications (1)

Publication Number Publication Date
CN111781948A true CN111781948A (zh) 2020-10-16

Family

ID=72757475

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010559970.0A Withdrawn CN111781948A (zh) 2020-06-18 2020-06-18 一种无人机编队队形变换与在线动态避障方法

Country Status (1)

Country Link
CN (1) CN111781948A (zh)

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112379672A (zh) * 2020-11-24 2021-02-19 浙大宁波理工学院 一种基于改进人工势场的智能无人船路径规划方法
CN112799432A (zh) * 2021-04-08 2021-05-14 北京三快在线科技有限公司 无人机的避障控制方法、装置、存储介质和电子设备
CN112925342A (zh) * 2021-01-20 2021-06-08 北京工商大学 基于改进相互速度障碍法的无人机动态避障方法
CN113119116A (zh) * 2021-03-22 2021-07-16 深圳市优必选科技股份有限公司 一种机械臂运动规划方法、装置、可读存储介质及机械臂
CN113282103A (zh) * 2021-05-24 2021-08-20 河北科技大学 基于改进自适应调整阈值势场法的无人机冲突探测和解脱方法
CN113359859A (zh) * 2021-07-16 2021-09-07 广东电网有限责任公司 一种组合导航避障方法、***、终端设备及存储介质
CN113485454A (zh) * 2021-08-23 2021-10-08 南京航空航天大学 几何改进人工势场的无人机避障机动方法和***
CN113568428A (zh) * 2021-07-12 2021-10-29 中国科学技术大学 一种基于多无人机协同的校园安防方法及***
CN113655810A (zh) * 2021-08-20 2021-11-16 上海微电机研究所(中国电子科技集团公司第二十一研究所) 一种基于速度势场的无人机避障方法和***
CN113655789A (zh) * 2021-08-04 2021-11-16 东风柳州汽车有限公司 路径跟踪方法、装置、车辆及存储介质
CN114527791A (zh) * 2022-01-21 2022-05-24 重庆工业大数据创新中心有限公司 一种禁飞区约束下的多无人机包裹协作递送路径优化方法
CN114764247A (zh) * 2022-04-25 2022-07-19 珠海云洲智能科技股份有限公司 无人驾驶设备集群的控制方法、装置、设备及存储介质
CN115167528A (zh) * 2022-09-05 2022-10-11 北京航空航天大学 一种基于人工势场法的空间协同制导方法及装置
CN116400722A (zh) * 2023-05-10 2023-07-07 江苏方天电力技术有限公司 一种无人机避障飞行方法及相关装置
CN116661502A (zh) * 2023-07-25 2023-08-29 吉林农业大学 一种智能农业无人机路径规划方法
CN117608318A (zh) * 2024-01-23 2024-02-27 北京航空航天大学 基于仿鸟类趋光性的无人机编队避障控制方法及***

Cited By (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112379672A (zh) * 2020-11-24 2021-02-19 浙大宁波理工学院 一种基于改进人工势场的智能无人船路径规划方法
CN112925342A (zh) * 2021-01-20 2021-06-08 北京工商大学 基于改进相互速度障碍法的无人机动态避障方法
CN113119116A (zh) * 2021-03-22 2021-07-16 深圳市优必选科技股份有限公司 一种机械臂运动规划方法、装置、可读存储介质及机械臂
CN112799432A (zh) * 2021-04-08 2021-05-14 北京三快在线科技有限公司 无人机的避障控制方法、装置、存储介质和电子设备
CN112799432B (zh) * 2021-04-08 2021-07-02 北京三快在线科技有限公司 无人机的避障控制方法、装置、存储介质和电子设备
CN113282103A (zh) * 2021-05-24 2021-08-20 河北科技大学 基于改进自适应调整阈值势场法的无人机冲突探测和解脱方法
CN113568428A (zh) * 2021-07-12 2021-10-29 中国科学技术大学 一种基于多无人机协同的校园安防方法及***
CN113359859A (zh) * 2021-07-16 2021-09-07 广东电网有限责任公司 一种组合导航避障方法、***、终端设备及存储介质
CN113359859B (zh) * 2021-07-16 2023-09-08 广东电网有限责任公司 一种组合导航避障方法、***、终端设备及存储介质
CN113655789A (zh) * 2021-08-04 2021-11-16 东风柳州汽车有限公司 路径跟踪方法、装置、车辆及存储介质
CN113655810A (zh) * 2021-08-20 2021-11-16 上海微电机研究所(中国电子科技集团公司第二十一研究所) 一种基于速度势场的无人机避障方法和***
CN113655810B (zh) * 2021-08-20 2024-04-16 上海微电机研究所(中国电子科技集团公司第二十一研究所) 一种基于速度势场的无人机避障方法和***
CN113485454A (zh) * 2021-08-23 2021-10-08 南京航空航天大学 几何改进人工势场的无人机避障机动方法和***
CN114527791A (zh) * 2022-01-21 2022-05-24 重庆工业大数据创新中心有限公司 一种禁飞区约束下的多无人机包裹协作递送路径优化方法
CN114764247A (zh) * 2022-04-25 2022-07-19 珠海云洲智能科技股份有限公司 无人驾驶设备集群的控制方法、装置、设备及存储介质
CN114764247B (zh) * 2022-04-25 2023-08-29 珠海云洲智能科技股份有限公司 无人驾驶设备集群的控制方法、装置、设备及存储介质
CN115167528A (zh) * 2022-09-05 2022-10-11 北京航空航天大学 一种基于人工势场法的空间协同制导方法及装置
CN116400722A (zh) * 2023-05-10 2023-07-07 江苏方天电力技术有限公司 一种无人机避障飞行方法及相关装置
CN116661502A (zh) * 2023-07-25 2023-08-29 吉林农业大学 一种智能农业无人机路径规划方法
CN116661502B (zh) * 2023-07-25 2023-10-10 吉林农业大学 一种智能农业无人机路径规划方法
CN117608318A (zh) * 2024-01-23 2024-02-27 北京航空航天大学 基于仿鸟类趋光性的无人机编队避障控制方法及***
CN117608318B (zh) * 2024-01-23 2024-04-09 北京航空航天大学 基于仿鸟类趋光性的无人机编队避障控制方法及***

Similar Documents

Publication Publication Date Title
CN111781948A (zh) 一种无人机编队队形变换与在线动态避障方法
CN109871031B (zh) 一种固定翼无人机的轨迹规划方法
CN109508035B (zh) 基于分布式控制的多区域分级式无人机编队路径规划方法
Baras et al. Decentralized control of autonomous vehicles
Zhao et al. Flocking control of fixed-wing UAVs with cooperative obstacle avoidance capability
Cole et al. Reactive trajectory generation for multiple vehicles in unknown environments with wind disturbances
Lai et al. Enhanced DWA algorithm for local path planning of mobile robot
Frew et al. Obstacle avoidance with sensor uncertainty for small unmanned aircraft
CN111880574A (zh) 一种无人机避撞方法及***
Zhao et al. Four-dimensional trajectory generation for UAVs based on multi-agent Q learning
Mousavi et al. Robust optimal higher-order-observer-based dynamic sliding mode control for VTOL unmanned aerial vehicles
Lin et al. Potential-field-based evolutionary route planner for the control of multiple unmanned aerial vehicles
Zhang et al. Collision Avoidance of Fixed-Wing UAVs in Dynamic Environments Based on Spline-RRT and Velocity Obstacle
Yang et al. Path planning and collision avoidance methods for distributed multi-robot systems in complex dynamic environments
CN114138022A (zh) 一种基于精英鸽群智能的无人机集群分布式编队控制方法
CN111121804B (zh) 一种具有安全性约束的智能车辆路径规划方法及***
Song et al. Multi-rotor UAVs path planning method based on improved artificial potential field method
CN115420289B (zh) 一种基于粒子群改进人工势场法的无人艇航路规划方法
Luo et al. Enhanced artificial potential field-based moving obstacle avoidance for UAV in three-dimensional environment
Kahn et al. Formation flying control via elliptical virtual structure
Xiang et al. Application of Improved Q-Learning Algorithm in Dynamic Path Planning for Aircraft at Airports
Zhang et al. Research on collision avoidance of fixed-wing UAV
Wang et al. Improved A* and fuzzy dynamic window based dynamic trajectory planning for an UAV
Meng et al. A navigation framework for mobile robots with 3D LiDAR and monocular camera
Li et al. Dynamic path planning of intelligent robot in power equipment maintenance environment

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
WW01 Invention patent application withdrawn after publication
WW01 Invention patent application withdrawn after publication

Application publication date: 20201016