CN113311869B - 一种多约束条件下无人机路径规划方法 - Google Patents

一种多约束条件下无人机路径规划方法 Download PDF

Info

Publication number
CN113311869B
CN113311869B CN202110592914.1A CN202110592914A CN113311869B CN 113311869 B CN113311869 B CN 113311869B CN 202110592914 A CN202110592914 A CN 202110592914A CN 113311869 B CN113311869 B CN 113311869B
Authority
CN
China
Prior art keywords
point
correction
error
aerial vehicle
unmanned aerial
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.)
Active
Application number
CN202110592914.1A
Other languages
English (en)
Other versions
CN113311869A (zh
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.)
Huaiyin Institute of Technology
Original Assignee
Huaiyin Institute of Technology
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 Huaiyin Institute of Technology filed Critical Huaiyin Institute of Technology
Priority to CN202110592914.1A priority Critical patent/CN113311869B/zh
Publication of CN113311869A publication Critical patent/CN113311869A/zh
Application granted granted Critical
Publication of CN113311869B publication Critical patent/CN113311869B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/106Change initiated in response to external conditions, e.g. avoidance of elevated terrain or of no-fly zones
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/10Internal combustion engine [ICE] based vehicles
    • Y02T10/40Engine management systems

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

本发明公开一种多约束条件下无人机路径规划方法,属于路径规划技术领域,包括:根据无人机校正点位置信息,建立多约束路径规划数学模型,将其转化为双目标规划模型,确定约束函数与约束条件;运用贪心算法对无人机最优化校正点进行求解,得到最优校正点和最短飞行距离。与现有技术相比,本发明建立最小校正次数和最短距离的双目标规划模型,优化校正点校正类型,同时利用贪心算法求解,该算法寻找局部最近三点作为无人机下一遍历节点,并在规定步长范围内,提高该算法的求解能力;本发明在二维平面上三点求解弧度曲线长度,使飞行曲线平滑过渡自然,但未过度增大路径距离,通过本发明可以获得较优的飞行路径。

Description

一种多约束条件下无人机路径规划方法
技术领域:
本发明属于无人机路径规划技术领域,具体涉及一种多约束条件下无人机路径规划方法。
背景技术
近年来,得益于人工智能的飞速发展及计算机技术的显著提升,无人机控制已成为军事作战任务的核心内容。无人机以其操作灵活性,在电力巡航、区域勘测、无损检测、抢险救灾、物资投放等作业中发挥着越来重要的作用,因此无人机路径规划也成为研究热点。在无人机执行作战任务前,根据不同目标位置有效且快速规划出最佳航迹路线能够缩短任务执行时间、提高能源利用率和扩展活动空间范围。然而,无人机飞行环境复杂多变,存在着诸如建筑物、山体障碍物及突如其来的不可抗因素,极大影响着无人机路径规划的合理性与飞行安全性。
在现实环境应用中,由于地理环境的复杂性和无人机自身结构等诸多因素,使得无人机面对复杂多变的外界环境快速路径规划便成为重中之重。无人机在固定的规划路径中,避免导航定位误差造成检测失败。因此,在无人机飞行过程中需要考虑外界环境,如何快速而最优的规划出水平误差、垂直误差、路径约束和无人机转弯半径条件约束下的空间节点,并使这些空间节点形成一条有效的路径规划尤为重要。
发明内容
技术目的:针对背景技术中存在的问题,本发明提出一种多约束条件下无人机路径规划方法,建立最小校正次数和最短距离的双目标规划模型,利用贪心算法求解,提高该算法的求解能力;在二维平面上三点求解弧度曲线长度,使飞行过渡自然,曲线平滑,但未过度增大路径距离,最终获得较优的飞行路径。
技术方案:本发明公开了一种多约束条件下无人机路径规划方法,包括以下步骤:
步骤1:首先分析无人机路径规划特点,建立三维数学模型,忽略无人机大小,将其视为质点,用(xi,yi,zi)表示第i个校正点的坐标i=1,2,3,…,n,起点A为(x0,y0,z0),终点B点(xn+1,yn+1,zn+1),无人机在校正过程中,从A出发,经过k个校正点到达终点,轨迹可表示为A点与k个校正点到B点依次相连的轨迹,建立所述三维数学模型的目标函数并确定其约束条件:
Figure BDA0003089874640000021
min k (2)
其中,公式(1)表示飞行轨迹尽可能小,公式(2)表示矫正次数尽可能少,
Figure BDA0003089874640000022
表示最小轨迹为起点A(x0,y0,z0)到第一个校正点i1的距离;
Figure BDA0003089874640000023
表示第1个到第k个校正点之间相邻校正点累积距离之和;
Figure BDA0003089874640000024
表示第k个校正点到目的点B(xn+1,yn+1,zn+1)的距离;if表示为第f个校正点;
Figure BDA0003089874640000025
表示为第if-1个校正点到第if个校正点的距离;
其约束条件为:
1)用ξi表示Ci校正的类型,无人机飞行过程中的定位误差有垂直误差和水平误差:
Figure BDA0003089874640000026
2)无人机在出发点处不存在定位误差,即水平误差和垂直误差为0;
3)到达终点B时垂直误差和水平误差均应小于θ个单位:
Figure BDA0003089874640000027
其中,
Figure BDA0003089874640000028
Figure BDA0003089874640000029
分别表示经过校正点
Figure BDA00030898746400000210
处调整后的水平误差和垂直误差;无人机每飞行1m,垂直误差和水平误差增加δ个专用单位;
4)飞行器在转弯时受到结构和控制系统的限制,无法完成即时转弯,利用转弯模型算法计算无人机的转弯半径,其约束半径r为:
(a-x1)2+(b-y1)2+(c-z1)2=r2
其中,(a,b,c)为圆心坐标O1坐标,(x1,y1,z1)为第一个校正点的坐标;
5)无人机在部分校正点进行误差校正时存在无法达到理想校正点,即将某个误差精确校正为0,用ηi表示校正点Ci是否为理想校正点:
Figure BDA00030898746400000211
根据实时位置,获取校正点坐标和校正误差,实时判断是否在设置范围内,并判断误差校正的类型;
步骤2:确定双目标规划模型,运用贪心算法进行全局搜索最优解对无人机最优化校正点进行求解,得到最优校正点和最短飞行距离。
进一步地,无人机在出发点A原始误差为0,到达第一个校正点i1前无需校正;同时由最后一个校正点ik到达终点B时定位误差小于θ也无需矫正。
进一步地,所述最小轨迹为A(x0,y0,z0)到第一个校正点i1的距离:
Figure BDA0003089874640000031
第一个校正点到第k个校正点的距离之和为:
Figure BDA0003089874640000032
其中,if表示为第f个校正点;
Figure BDA0003089874640000033
表示为第if-1个校正点到第if个校正点的距离;
第k个校正点到目的点B(xn+1,yn+1,zn+1)的距离:
Figure BDA0003089874640000034
进一步地,在出发点A,无人机的水平误差和垂直误差均为0,无人机从第if-1个校正点到if校正点的距离为
Figure BDA0003089874640000035
产生的积累误差为
Figure BDA0003089874640000036
同时在到达校正点if时含有if-1未校正的原始误差
Figure BDA0003089874640000037
无人机在校正点if校正前的垂直误差
Figure BDA0003089874640000038
为积累误差
Figure BDA0003089874640000039
与原始误差rif-1之和,水平误差
Figure BDA00030898746400000310
与垂直误差
Figure BDA00030898746400000311
′可表示为:
Figure BDA00030898746400000312
情况一,若校正点
Figure BDA00030898746400000313
是垂直误差校正点,需满足:
Figure BDA00030898746400000314
经过垂直误差校正点
Figure BDA00030898746400000315
校正后的误差为:
Figure BDA00030898746400000316
其中,α1为垂直误差设定值,α2为水平误差设定值,且都为常数;
情况二,若校正点
Figure BDA0003089874640000044
是水平误差校正点,需满足:
Figure BDA0003089874640000041
其中,β1为垂直误差设定值,β2为水平误差设定值,且都为常数;
经过水平误差校正点
Figure BDA0003089874640000042
校正后的误差为:
Figure BDA0003089874640000043
进一步地,所述步骤3中求解具体过程为:
第一步:赋初值;设置起点A、终点B,无人机垂直和水平误差均为0,s0=0,r0=0;
第二步:计算任意的一个校正点飞到离自身最近的一个校正点飞行积累下来的误差,判断累积下来的误差是否在允许范围内,若水平误差在允许范围内,累计误差为垂直误差则进行垂直校正,若在垂直误差允许范围内,累计误差为水平误差则进行水平校正,否则,该点为不可校正点;
第三步:从当前点出发,找到所有垂直校正点集合为Conset,判断Conset是否为空集,如果Conset是空集,则显示“路径寻找失败”,并且退出循环;如果Conset不是空集,则在集合Conset寻找距离终点B最近的一个校正点,记为Ck
第四步:设置起点为Ck,记录调整前后的误差,若Ck不是理想校正点,用均匀分布随机数ζ来仿真纠正是否成功校正后水平剩余误差均为0,记录校正前后的误差,并设置校正后的误差值;
第五步:判断当前点能否到达终点B,如果能到达B则输出最终路线,并退出循环;如果不能到达B,转到下一步;
第六步:从当前点出发,找到所有水平校正点集合为Conset′,判断Conset′是否为空集,如果Conset′是空集,则显示“路径寻找失败”,并且退出循环;如果Conset′不是空集,则在集合Conset′寻找距离终点B最近的一个校正点,记为Ck′;
第七步:设置起点为Ck′,记录调整前后的误差,若Ck′不是理想校正点,用均匀分布随机数ζ来仿真纠正是否成功,校正后水平剩余误差均为0,记录校正前后的误差,并设置校正后的误差值;
第八步:判断当前点能否到达终点B,如果能到达B则输出最终路线,并退出循环;如果不能到达B,转到第二步。
进一步地,所述第四步与第七步中用服从(0,1)均匀分布的随机数ζ纠正具体为:若ζ∈[0,0.2]则视为校正失败,校正后的剩余误差为设为min(error,5),其中error为校正前误差;反之,ζ∈(0.2,1]视为校正成功。
有益效果:
1、本发明针对无人机在不同实际飞行环境中有效路径规划方法,建立最优路径和最小校正次数为目标函数,定位误差约束,转弯控制受限约束和无法达到理想校正点约束条件下的多种约束构建数学模型,利用贪心算法对无人机路径进行快速规划,垂直校正点和水平校正点可校正的误差允许范围不相同,且达到终点时也在允许范围内,路径规划更加符合实际飞行过程。
2、寻找局部最近三点作为求解圆弧长度,将三维空间转化为二维平面上的圆弧曲线长度求解,使飞行过渡平滑更加自然,但并未过多增大飞行路径,可获得最优路径。
附图说明
图1为本发明实施例无人机航迹(以C标记的为校正点);
图2为本发明实施例求取无人机转弯轨迹的切点与转弯半径示意图;
图3为本发明实施例无人机选择的航迹路径示意图;
图4为本发明实施例无人机转弯放大后的轨迹示意图;
图5为本发明模型求解算法运行1000次的路径图。
具体实施方式
下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
本发明公开了一种多约束条件下无人机路径规划方法,具体包括以下步骤:
步骤1:先分析无人机路径规划特点,建立多约束条件的数学模型;
建立三维数学模型,用(xi,yi,zi)表示第i个校正点Ci的坐标i=1,2,3,…,n。A点坐标为(x0,y0,z0),B点坐标(xn+1,yn+1,zn+1),用ξi表示Ci校正的类型:
Figure BDA0003089874640000051
无人机在校正过程中,从A点出发,经过k个校正点到达终点B,轨迹可表示为A点与k个校正点到B点依次相连的轨迹,表示为
Figure BDA0003089874640000061
其中,
Figure BDA0003089874640000062
表示校正点Ci中的第k个子集元素。
最小轨迹为A(x0,y0,z0)到第一个校正点i1的距离:
Figure BDA0003089874640000063
第一个校正点
Figure BDA0003089874640000064
到第k个校正点
Figure BDA0003089874640000065
的距离之和为:
Figure BDA0003089874640000066
其中,if表示为第f个校正点;
Figure BDA0003089874640000067
表示为第if-1个采样校正点到第if个采样校正点的距离。
第k个校正点到目的点B(xn+1,yn+1,zn+1)的距离:
Figure BDA0003089874640000068
无人机路径规划过程中要考虑两个条件:
1)第一个是:飞行轨迹尽可能小,即:
Figure BDA0003089874640000069
2)第二个条件是:矫正次数尽可能少,即
min k。
无人机每飞行1m,垂直误差和水平误差增加δ个专用单位。si和ri分别表示经过校正点Ci处调整后的水平误差和垂直误差。水平误差si′和垂直误差ri′表示在经过校正点Ci处调整前的误差。在出发A点,无人机的水平误差和垂直误差均为0,所以有s0=0和r0=0。
无人机从第if-1个校正点到if校正点的距离为
Figure BDA00030898746400000610
产生的积累误差为
Figure BDA00030898746400000611
同时在到达校正点if时含有if-1未校正的原始误差
Figure BDA00030898746400000612
无人机在校正点if校正前的垂直误差
Figure BDA00030898746400000613
为积累误差
Figure BDA00030898746400000614
与原始误差
Figure BDA00030898746400000615
之和。水平误差
Figure BDA00030898746400000616
与垂直误差
Figure BDA00030898746400000617
可表示为:
Figure BDA0003089874640000071
情况一:若校正点
Figure BDA0003089874640000072
是垂直误差校正点,需满足:
Figure BDA0003089874640000073
其中,α1为垂直误差设定值,α2为水平误差设定值,且都为常数;
经过垂直误差校正点
Figure BDA0003089874640000074
校正后的误差为:
Figure BDA0003089874640000075
情况二:若校正点
Figure BDA0003089874640000076
是水平误差校正点
Figure BDA0003089874640000077
其中,β1为垂直误差设定值,β2为水平误差设定值,且都为常数;
经过水平误差校正点
Figure BDA0003089874640000078
校正后的误差为:
Figure BDA0003089874640000079
到达终点时垂直误差和水平误差均应小于θ个单位即可,所以有:
Figure BDA00030898746400000710
其中,
Figure BDA00030898746400000711
Figure BDA00030898746400000712
分别表示经过校正点
Figure BDA00030898746400000713
处调整后的水平误差和垂直误差。
其次无人机在转弯时受到结构和控制系统的限制,无法完成即时转弯(飞行器前进方向无法突然改变),且转弯时不能直角转弯,故需增加约束条件,转弯半径不能小于200米,故对飞行轨迹进行优化。
利用转弯模型算法计算无人机的转弯半径,其约束半径r为:
(a-x1)2+(b-y1)2+(c-z1)2=r2
其中,(a,b,c)为圆心坐标O1坐标,(x1,y1,z1)为第一个校正点的坐标。
下面需要对空间上k个校正点轨迹进行优化,将三维空间飞行转变为二维平面求解。参见附图1:
1)从出发点A出发时无误差,无人机可直接飞到校正点C1,线段AC1为第一段飞行轨迹。
2)不在同一直线上的三点A、C1和C2,确定一个平面AC1C2;在平面AC1C2上,以直线AC1为切线,C1为切点做半径为200米的圆O1(圆心),过C2点的做圆01的切线,切点为F1(在圆01上的属于非校正点切点),圆弧C1F1即为旋转时轨迹(第二段轨迹)。
3)线段C2F1为第三段轨迹。
4)接着以C2F1为切线做半径为200米的圆O2,再过C3点的做圆O2的切线,切点为F2(在圆O2上的属于非校正点切点),圆弧C3F2即为旋转时轨迹(第四段轨迹)。剩下校正点以此类推,直至B点。
参见附图2,具体过程为:
1.确定平面:
参见附图2,设A、C1、C2三点坐标为(x0,y0,z0)、(x1,y1,z1)、(x2,y2,z2),根据三点坐标AC1C2的方程为
Figure BDA0003089874640000081
圆心坐标O1坐标为(a,b,c),由3个条件:
a)点在平面AC1C2
Figure BDA0003089874640000082
b)线段C1O1长度等于200米
(a-x1)2+(b-y1)2+(c-z1)2=2002
c)向量
Figure BDA0003089874640000083
与向量
Figure BDA0003089874640000084
垂直
(a-x1)(x1-x0)+(b-y1)(y1-y0)+(c-z1)(z1-z0)=0
可知,存在两个符合条件的圆心O1和01′,需要选择一个合适初值点作为参考点。有效的圆心点所在位置可由
Figure BDA0003089874640000085
设R1坐标为(XR1,YR1,ZR1),则:
Figure BDA0003089874640000086
2.求切点F1坐标
设F1坐标为(d,e,f),由下面的3个条件:
a)F1点在平面AC1C2
Figure BDA0003089874640000091
b)线段O1F1长度等于200米
(a-d)2+(b-e)2+(c-f)2=2002
c)向量
Figure BDA0003089874640000092
与向量
Figure BDA0003089874640000093
垂直
(a-x2)(d-a)+(b-y2)(e-b)+(c-z2)(f-c)=0
可知,存在两个符合条件的切点F1与F1′,但只有一个为有效切点。有效的切点所在位置可由向量
Figure BDA0003089874640000094
设R2坐标为(XR2,YR2,ZR2),则:
Figure BDA0003089874640000095
参见附图4,看出无人机路径在直线飞行与转弯飞行时过渡自然,曲线平滑,效果更好。
无人机在部分校正点进行误差校正时存在无法达到理想校正点(即将某个误差精确校正为0)。假设无人机在部分校正点能够成功将某个误差校正为0的概率为80%,如果校正失败,则校正后的剩余误差为min(error,5)个单位。用ηi表示校正点Ci的是否为理想校正点,即:
Figure BDA0003089874640000096
综上所述,多约束条件转化为双目标规划模型如下:
Figure BDA0003089874640000101
α1=25,α2=15,β1=20,β2=25,θ=30,δ=0.001,当无人机的垂直误差不大于α1个单位,水平误差不大于α2个单位时才能进行垂直误差校正。当无人机的垂直误差不大于β1个单位,水平误差不大于β2个单位时才能进行水平误差校正。
下面以垂直校正再水平校正误差进行对上述的双目标规划模型进行求解,其算法思想如下:
第一步:赋初值;设置起点A,设置终点B,无人机垂直和水平误差均为0,s0=0,r0=0;
第二步:计算校正点之间距离,判断飞到自身最近的一个校正点累积下来的误差是否在允许范围内,如果水平误差在允许范围内,累计误差为垂直误差则进行垂直校正,如果在垂直误差允许范围内,累计误差为水平误差则进行水平校正,否则,该点为不可校正点;
第三步:从当前点出发,找到所有垂直校正点集合为Conset,判断Conset是否为空集,如果Conset是空集,则显示“路径寻找失败”,并且退出循环;如果Conset不是空集,则在集合Conset寻找距离终点B最近的一个校正点,记为Ck
第四步:设置起点为Ck,记录调整前后的误差,若Ck不是理想校正点,用均匀分布随机数ζ来仿真纠正是否成功,若ζ∈[0,0.2]则视为校正失败,校正后的剩余误差为设为min(error,5),其中error为校正前误差;反之,ζ∈(0.2,1]视为校正成功,校正后水平剩余误差均为0,记录校正前后的误差,并设置校正后的误差值;
第五步:判断当前点能否到达终点B,如果能到达B则输出最终路线,并退出循环;如果不能到达B,转到下一步;
第六步:从当前点出发,找到所有水平校正点集合为Conset′,判断Conset′是否为空集,如果Conset′是空集,则显示“路径寻找失败”,并且退出循环;如果Conset′不是空集,则在集合Conset′寻找距离终点B最近的一个校正点,记为Ck′;
第七步:设置起点为Ck′,记录调整前后的误差,若Ck′不是理想校正点,用均匀分布随机数ζ来仿真纠正是否成功,若ζ∈[0,0.2]则视为校正失败,校正后的剩余误差为设为min(error,5),其中error为校正前误差;反之,ζ∈(0.2,1]视为校正成功,校正后水平剩余误差均为0,记录校正前后的误差,并设置校正后的误差值;
第八步:判断当前点能否到达终点B,如果能到达B则输出最终路线,并退出循环;如果不能到达B,转到第二步。
第九步:无人机校正点集合为
Figure BDA0003089874640000111
设置飞行航迹上的切点集合
Figure BDA0003089874640000112
集合中
Figure BDA0003089874640000118
第十步:设置三点为
Figure BDA0003089874640000113
用A1、A2和A3确定圆心O,再确定切点
Figure BDA0003089874640000114
求出圆弧长度
Figure BDA0003089874640000115
第十一步:判断A3是否等于
Figure BDA0003089874640000116
若相等则A3=B,退出循环;否则转到第八步,直至结束。
程序循环运行1000次,参见附图5所示的路径图,统计成功路径规划的次数,从而得到该算法的成功率96.70%。
经过校正前后的误差结果对比表
Figure BDA0003089874640000117
Figure BDA0003089874640000121
基于优化转弯时的算法,仍使得校正点的个数为8,与前面的校正点个数相等,且飞行轨迹为110185.13米,只比未优化的飞行轨迹多107.89米,比对结果优良,因此该算法的性能较好。
上述实施方式只为说明本发明的技术构思及特点,其目的在于让熟悉此项技术的人能够了解本发明的内容并据以实施,并不能以此限制本发明的保护范围。凡根据本发明精神实质所做的等效变换或修饰,都应涵盖在本发明的保护范围之内。

Claims (5)

1.一种多约束条件下无人机路径规划方法,其特征在于,包括以下步骤:
步骤1:首先分析无人机路径规划特点,建立三维数学模型,忽略无人机大小,将其视为质点,用(xi,yi,zi)表示第i个校正点的坐标i=1,2,3,…,n,起点A为(x0,y0,z0),终点B为(xn+1,yn+1,zn+1),无人机在校正过程中,从A出发,经过k个校正点到达终点,轨迹可表示为A点与k个校正点到B点依次相连的轨迹,建立所述三维数学模型的目标函数并确定其约束条件:
Figure FDA0003887104620000011
mink (2)
其中,公式(1)表示飞行轨迹尽可能小,公式(2)表示矫正次数尽可能少,
Figure FDA0003887104620000012
表示最小轨迹为起点A(x0,y0,z0)到第一个校正点i1的距离;
Figure FDA0003887104620000013
表示第1个到第k个校正点之间相邻校正点累积距离之和;
Figure FDA0003887104620000014
表示第k个校正点到目的点B(xn+1,yn+1,zn+1)的距离;if表示为第f个校正点;
Figure FDA0003887104620000015
表示为第if-1个校正点到第if个校正点的距离;
所述最小轨迹为A(x0,y0,z0)到第一个校正点i1的距离:
Figure FDA0003887104620000016
第一个校正点到第k个校正点的距离之和为:
Figure FDA0003887104620000017
其中,if表示为第f个校正点;
Figure FDA0003887104620000018
表示为第if-1个校正点到第if个校正点的距离;
第k个校正点到目的点B(xn+1,yn+1,zn+1)的距离:
Figure FDA0003887104620000019
在出发点A,无人机的水平误差和垂直误差均为0,无人机从第if-1个校正点到if校正点的距离为
Figure FDA00038871046200000110
产生的积累误差为
Figure FDA00038871046200000111
同时在到达校正点if时含有if-1未校正的原始误差
Figure FDA00038871046200000112
无人机在校正点if校正前的垂直误差
Figure FDA00038871046200000113
为积累误差
Figure FDA00038871046200000114
与原始误差
Figure FDA00038871046200000115
之和,水平误差
Figure FDA00038871046200000116
与垂直误差
Figure FDA00038871046200000117
可表示为:
Figure FDA0003887104620000021
其约束条件为:
1)用ξi表示Ci校正的类型,无人机飞行过程中的定位误差有垂直误差和水平误差:
Figure FDA0003887104620000022
2)无人机在出发点处不存在定位误差,即水平误差和垂直误差为0;
3)到达终点B时垂直误差和水平误差均应小于θ个单位:
Figure FDA0003887104620000023
其中,
Figure FDA0003887104620000024
Figure FDA0003887104620000025
分别表示经过校正点
Figure FDA0003887104620000026
处调整后的水平误差和垂直误差;无人机每飞行1m,垂直误差和水平误差增加δ个专用单位;
4)飞行器在转弯时受到结构和控制系统的限制,无法完成即时转弯,利用转弯模型算法计算无人机的转弯半径,其约束半径r为:
(a-x1)2+(b-y1)2+(c-z1)2=r2
其中,(a,b,c)为圆心坐标O1坐标,(x1,y1,z1)为第一个校正点的坐标;
转弯半径不能小于200米 ,对飞行轨迹进行优化,具体优化步骤为:
(1)从出发点A出发时无误差,无人机可直接飞到校正点C1,线段AC1为第一段飞行轨迹;
(2)不在同一直线上的三点A、C1和C2,确定一个平面AC1C2;在平面AC1C2上,以直线AC1为切线,C1为切点做半径为200米的圆O1,O1为圆心,过C2点的做圆O1的切线,切点为F1,在圆O1上的属于非校正点切点,圆弧C1F1即为旋转时轨迹,即为第二段轨迹;
(3)线段C2F1为第三段轨迹;
(4)接着以C2F1为切线做半径为200米的圆O2,再过C3点的做圆O2的切线,切点为F2,在圆O2上的属于非校正点切点,圆弧C3F2即为旋转时轨迹,即第四段轨迹,剩下校正点以此类推,直至B点;
设A、C1、C2三点坐标为(x0,y0,z0)、(x1,y1,z1)、(x2,y2,z2),根据三点坐标AC1C2的方程为
Figure FDA0003887104620000031
圆心坐标O1坐标为(a,b,c),由3个条件:
a)点在平面AC1C2
Figure FDA0003887104620000032
b)线段C1O1长度等于200米
(a-x1)2+(b-y1)2+(c-z1)2=2002
c)向量
Figure FDA0003887104620000033
与向量
Figure FDA0003887104620000034
垂直
(a-x1)(x1-x0)+(b-y1)(y1-y0)+(c-z1)(z1-z0)=0
可知,存在两个符合条件的圆心O1和O1′,需要选择一个合适初值点作为参考点,有效的圆心点所在位置可由
Figure FDA0003887104620000035
设R1坐标为(XR1,YR1,ZR1),则:
Figure FDA0003887104620000036
求切点F1坐标:
设F1坐标为(d,e,f),由下面的3个条件:
a)F1点在平面AC1C2
Figure FDA0003887104620000037
b)线段O1F1长度等于200米
(a-d)2+(b-e)2+(c-f)2=2002
c)向量
Figure FDA0003887104620000038
与向量
Figure FDA0003887104620000039
垂直
(a-x2)(d-a)+(b-y2)(e-b)+(c-z2)(f-c)=0
可知,存在两个符合条件的切点F1与F1′,但只有一个为有效切点,有效的切点所在位置可由向量
Figure FDA00038871046200000310
设R2坐标为(XR2,YR2,ZR2),则:
Figure FDA0003887104620000041
5)无人机在部分校正点进行误差校正时存在无法达到理想校正点,即将某个误差精确校正为0,用ηi表示校正点Ci是否为理想校正点:
Figure FDA0003887104620000042
根据实时位置,获取校正点坐标和校正误差,实时判断是否在设置范围内,并判断误差校正的类型;
步骤2:确定双目标规划模型,运用贪心算法进行全局搜索最优解对无人机最优化校正点进行求解,得到最优校正点和最短飞行距离
多约束条件转化为双目标规划模型如下:
Figure FDA0003887104620000043
2.根据权利要求1所述的多约束条件下无人机路径规划方法,其特征在于,无人机在出发点A原始误差为0,到达第一个校正点i1前无需校正;同时由最后一个校正点ik到达终点B时定位误差小于θ也无需矫正。
3.根据权利要求1所述的多约束条件下无人机路径规划方法,其特征在于,
情况一,若校正点
Figure FDA0003887104620000044
是垂直误差校正点,需满足:
Figure FDA0003887104620000051
经过垂直误差校正点
Figure FDA0003887104620000052
校正后的误差为:
Figure FDA0003887104620000053
其中,α1为垂直误差设定值,α2为水平误差设定值,且都为常数;
情况二,若校正点
Figure FDA0003887104620000054
是水平误差校正点,需满足:
Figure FDA0003887104620000055
其中,β1为垂直误差设定值,β2为水平误差设定值,且都为常数;
经过水平误差校正点
Figure FDA0003887104620000056
校正后的误差为:
Figure FDA0003887104620000057
4.根据权利要求1所述的多约束条件下无人机路径规划方法,其特征在于,所述步骤2中求解具体过程为:
第一步:赋初值;设置起点A、终点B,无人机垂直和水平误差均为0,s0=0,r0=0;
第二步:计算任意的一个校正点飞到离自身最近的一个校正点飞行积累下来的误差,判断累积下来的误差是否在允许范围内,若水平误差在允许范围内,累计误差为垂直误差则进行垂直校正,若在垂直误差允许范围内,累计误差为水平误差则进行水平校正,否则,该点为不可校正点;
第三步:从当前点出发,找到所有垂直校正点集合为Conset,判断Conset是否为空集,如果Conset是空集,则显示“路径寻找失败”,并且退出循环;如果Conset不是空集,则在集合Conset寻找距离终点B最近的一个校正点,记为Ck
第四步:设置起点为Ck,记录调整前后的误差,若Ck不是理想校正点,用均匀分布随机数ζ来仿真纠正是否成功校正后水平剩余误差均为0,记录校正前后的误差,并设置校正后的误差值;
第五步:判断当前点能否到达终点B,如果能到达B则输出最终路线,并退出循环;如果不能到达B,转到下一步;
第六步:从当前点出发,找到所有水平校正点集合为Conset′,判断Conset′是否为空集,如果Conset′是空集,则显示“路径寻找失败”,并且退出循环;如果Conset′不是空集,则在集合Conset′寻找距离终点B最近的一个校正点,记为Ck′;
第七步:设置起点为Ck′,记录调整前后的误差,若Ck′不是理想校正点,用均匀分布随机数ζ来仿真纠正是否成功,校正后水平剩余误差均为0,记录校正前后的误差,并设置校正后的误差值;
第八步:判断当前点能否到达终点B,如果能到达B则输出最终路线,并退出循环;如果不能到达B,转到第二步。
5.根据权利要求4所述的多约束条件下无人机路径规划方法,其特征在于,所述第四步与第七步中用服从(0,1)均匀分布的随机数ζ纠正具体为:若ζ∈[0,0.2]则视为校正失败,校正后的剩余误差为设为min(error,5),其中error为校正前误差;反之,ζ∈(0.2,1]视为校正成功。
CN202110592914.1A 2021-05-28 2021-05-28 一种多约束条件下无人机路径规划方法 Active CN113311869B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110592914.1A CN113311869B (zh) 2021-05-28 2021-05-28 一种多约束条件下无人机路径规划方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110592914.1A CN113311869B (zh) 2021-05-28 2021-05-28 一种多约束条件下无人机路径规划方法

Publications (2)

Publication Number Publication Date
CN113311869A CN113311869A (zh) 2021-08-27
CN113311869B true CN113311869B (zh) 2023-01-13

Family

ID=77376062

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110592914.1A Active CN113311869B (zh) 2021-05-28 2021-05-28 一种多约束条件下无人机路径规划方法

Country Status (1)

Country Link
CN (1) CN113311869B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114460969A (zh) * 2022-02-16 2022-05-10 中银金融科技有限公司 无人机路径规划方法、装置、电子设备及计算机存储介质
CN114545974B (zh) * 2022-03-02 2023-12-19 广东翼景信息科技有限公司 一种基于智能机场的有能力约束无人机弧路径规划方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110794874B (zh) * 2019-10-11 2022-12-30 东南大学 一种定位误差约束下的飞行器航迹快速规划方法
CN110715664B (zh) * 2019-11-05 2021-04-20 大连理工大学 多约束条件下智能无人飞行器航迹快速规划方法
CN111024092B (zh) * 2019-12-31 2020-10-30 西南交通大学 一种多约束条件下智能飞行器航迹快速规划方法
CN111595343B (zh) * 2020-04-29 2022-10-04 西北工业大学 一种基于定位误差校正的无人机航迹规划方法
CN112066991A (zh) * 2020-09-07 2020-12-11 常州微亿智造科技有限公司 多约束条件下智能飞行器航迹快速规划方法、装置
CN112729308B (zh) * 2020-12-24 2024-05-03 广州融赋数智技术服务有限公司 多约束条件下无人机快速航迹规划方法

Also Published As

Publication number Publication date
CN113311869A (zh) 2021-08-27

Similar Documents

Publication Publication Date Title
CN113311869B (zh) 一种多约束条件下无人机路径规划方法
CN110308740B (zh) 一种面向移动目标追踪的无人机群动态任务分配方法
CN111024092B (zh) 一种多约束条件下智能飞行器航迹快速规划方法
CN110244750B (zh) 一种无人机巡视路径规划方法及装置
CN111595343B (zh) 一种基于定位误差校正的无人机航迹规划方法
CN111207752B (zh) 一种基于动态切点调整的无人机航迹规划方法
CN111006693B (zh) 智能飞行器航迹规划系统及其方法
CN111580556A (zh) 一种时空约束下多无人机协同路径规划与制导方法
CN111352417A (zh) 异构多无人机协同路径的快速生成方法
CN110597286B (zh) 一种利用智慧机库实现输电线路无人机自主巡检的方法
CN103700286A (zh) 一种舰载无人机自动着舰引导方法
CN110032209A (zh) 一种多无人机任务区集结分布式控制方法
CN107656523A (zh) 仓储agv导航方法及装置、计算机设备及存储介质
CN109405830B (zh) 基于线路坐标序列的无人机自动巡检方法
CN112729308B (zh) 多约束条件下无人机快速航迹规划方法
CN111862200B (zh) 一种煤棚内无人机定位方法
CN113093801B (zh) 一种多无人机协同编队的航路文件制作方法
CN110794874B (zh) 一种定位误差约束下的飞行器航迹快速规划方法
CN112066992A (zh) 一种基于视场约束的反辐射无人机搜索航迹规划方法
CN111381605A (zh) 一种应用于多无人机大范围海域水下多目标协同搜索方法
CN111121779B (zh) 一种无人机所处飞行区域的实时检测方法
CN110109481B (zh) 基于轨道变换法的多无人飞行器跟踪目标相角修正方法
CN115933736A (zh) 一种航迹快速规划精准控制飞行器的方法
Khuralay et al. Computer simulation of intelligent control systems for high-precision cruise missiles
CN114237304B (zh) 一种基于Dubins动态路径规划的多无人机一致性集结方法

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
GR01 Patent grant
GR01 Patent grant
EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20210827

Assignee: Shanghai Yanqiao Information Technology Co.,Ltd.

Assignor: HUAIYIN INSTITUTE OF TECHNOLOGY

Contract record no.: X2023980047013

Denomination of invention: A Path Planning Method for Unmanned Aerial Vehicles under Multiple Constraints

Granted publication date: 20230113

License type: Common License

Record date: 20231115