CN111207752B - 一种基于动态切点调整的无人机航迹规划方法 - Google Patents

一种基于动态切点调整的无人机航迹规划方法 Download PDF

Info

Publication number
CN111207752B
CN111207752B CN202010050196.0A CN202010050196A CN111207752B CN 111207752 B CN111207752 B CN 111207752B CN 202010050196 A CN202010050196 A CN 202010050196A CN 111207752 B CN111207752 B CN 111207752B
Authority
CN
China
Prior art keywords
point
correction
aerial vehicle
unmanned aerial
vertical
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
CN202010050196.0A
Other languages
English (en)
Other versions
CN111207752A (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.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
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 Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN202010050196.0A priority Critical patent/CN111207752B/zh
Publication of CN111207752A publication Critical patent/CN111207752A/zh
Application granted granted Critical
Publication of CN111207752B publication Critical patent/CN111207752B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/20Instruments for performing navigational calculations
    • 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)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Automation & Control Theory (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
  • Traffic Control Systems (AREA)

Abstract

本发明属于无人机航迹规划技术领域,公开了一种基于动态切点调整的无人机航迹规划方法,基于动态切点调整的航迹规划方法,针对现有技术存在的问题,将航迹规划问题形式化为一个与定位误差校正次数、航迹长度和转弯半径相关的多目标优化问题,在保证定位误差得到有效校正的前提下优化航迹长度,并利用动态切点调整的方式来制定无人机的转弯策略。本发明可有效降低无人机航迹长度与算法复杂度并能符合无人机动力结构,相较于其他无人机航迹规划方法更能适用于实际场景应用,实现多约束条件下的无人机航迹快速规划。

Description

一种基于动态切点调整的无人机航迹规划方法
技术领域
本发明属于无人机航迹规划技术领域,尤其涉及一种基于动态切点调整的无人机航迹规划方法。
背景技术
无人机是一种有动力、可控制、能携带多种任务设备、执行多种任务并能重复使用的无人驾驶航空器。由于具备无人驾驶、尺寸小、成本低以及机动性强等优点,被广泛应用于各大军用和民用领域当中,如无人机通信、农业、火灾监控和检测和物联网等领域。无人机航迹规划技术是无人机应用的重要研究方向。自20世纪中期以来,一些西方国家已相继开始了无人机航迹规划技术的研究,并在20世纪70年代逐步构建了无人机飞行方面的一系列基础理论。而伴随着计算机技术的普遍应用和快速发展,从20世纪80年代中期以及后期开始,相关学者已将无人机航迹规划技术应用于遥控、遥感等领域,并在20世纪90年代,无人机航迹规划技术逐步趋于成熟。但现有的无人机航迹规划技术仍难以满足多样化的应用需求,无人机所用的航迹规划通常需要人工经验操纵处理,同时由于无人机系统结构限制,无人机的定位系统无法对自身进行精准定位,一旦定位误差积累到一定程度可能导致任务失败,并且无人机在转弯时受到结构和控制系统的限制,无法完成即时转弯。因此,在飞行过程中对定位误差进行校正并设计符合无人机动力结构的飞行轨迹是无人机航迹规划中的一项重要任务。
现有技术一提出运用几何数学的手段对无人机航迹进行规划,并在航迹规划过程中对转弯半径、转弯方向予以约束,该技术能快速规划出最短的无人机飞行路径,但该技术假设装载在无人机设备上的定位系统具有精准定位功能,不会产生定位误差。而在实际的无人机飞行过程中,定位误差是会逐渐累积的,一旦超过定位误差阈值则可能会导致无人机偏离预定轨迹,从而影响预期任务目标的实现。
现有技术二提出了与三维空间相结合的快速A*搜索算法,通过无源定位的方式对无人机航迹进行规划,该技术可以通过无源定位的方式来降低定位误差所带来的影响,但该技术所规划的航迹没有考虑无人机动力结构的限制。受限于无人机的动力结构,无人机在转弯时需有一定的转弯半径,无法实现即时转弯,因此该技术不能直接应用于实际的无人机航迹规划应用中。
现有技术三针对目前无人机在飞行任务中无法以最小路径代价实现航带间转弯飞行的问题,提出一种基于最小转弯半径的无人机转弯航迹规划算法。该技术依据无人机情报处理对于航片旁向重叠率的要求,确定航带间距,同时依据无人机飞行速度及姿态数据,确定无人机的最小转弯半径。该技术所规划的航迹不仅短而且符合无人机的动力结构,但该技术复杂程度较高,需根据不同的航带间距来制定不同的转弯策略,而现有无人机的计算能力有限,所规划的无人机航迹不能过于复杂,因此所规划的航迹不仅需符合无人机的动力结构还需尽可能的复杂度低。
综上所述,现有技术存在的问题是:
(1)现有技术一运用几何数学的手段对无人机航迹进行规划,没有考虑无人机定位存在误差的问题,而误差会逐步积累可能导致任务失败的情况发生。
(2)现有技术二通过无源定位的方式对无人机航迹进行规划,在航迹规划过程中没有考虑无人机动力结构的问题,无法满足无人机转弯半径所带来的约束。
(3)现有技术三基于最小转弯半径来对无人机航迹进行规划,但该技术复杂程度较高,针对不同的场景需制定不同的转弯策略,无法满足无人机航迹规划对低程度复杂度的要求。
解决上述技术问题的难度:现有技术一:该技术运用几何数学的手段对无人机航迹进行规划时需假设无人机时刻按预定路线进行飞行,而在飞行过程中产生定位误差是必然的,因此在规划航迹过程中需考虑定位误差所带来的影响。现有技术二:该技术采用快速A*搜索算法来规划航迹,并通过无源定位的方式来降低定位误差,但采用即时转弯的转弯策略不符合无人机的动力结构。现有技术三:该技术考虑无人机转弯半径,并根据不同的场景设计不同的转弯策略,但该技术复杂度较高。
解决上述技术问题的意义:本发明采用基于动态切点调整的航迹规划方法,针对现有技术存在的问题,将航迹规划问题形式化为一个与定位误差校正次数、航迹长度和转弯半径相关的多目标优化问题,在保证定位误差得到有效校正的前提下优化航迹长度,并利用动态切点调整的方式来制定无人机的转弯策略。本发明可有效降低无人机航迹长度与算法复杂度并能符合无人机动力结构,相较于其他无人机航迹规划方法更能适用于实际场景应用,实现多约束条件下的无人机航迹快速规划。
发明内容
针对现有技术存在的问题,本发明提供了一种基于动态切点调整的无人机航迹规划方法。
本发明是这样实现的,一种基于动态切点调整的无人机航迹规划方法,所述基于动态切点调整的无人机航迹规划问题形式化为一个与定位误差校正次数、航迹长度和转弯半径相关的多目标优化问题,在保证定位误差得到有效校正的前提下优化航迹长度,并利用动态切点调整的方式来制定无人机的转弯策略。
进一步,所述基于动态切点调整的无人机航迹规划方法包括以下步骤:
第一步,构建由1个出发点A、1个目的地B、M个水平校正点、N个垂直校正点组成的无人机航迹规划区域;
第二步,基于无人机航迹规划区域,构建校正点选择策略;
第三步,基于校正点选择策略,构建航迹规划方案;
第四步,基于航迹规划方案,构建基于动态切点调整的无人机航迹规划方法。
进一步,所述第一步构建由1个出发点A、1个目的点B、R个水平校正点、L个垂直校正点组成的无人机航迹规划区域,具体包括:
步骤(1.1),构建包含2+R+L点的无人机航迹规划区域,无人机在空间飞行过程中需要实时定位,其定位误差包括垂直误差和水平误差,无人机每飞行1m,垂直误差和水平误差将各增加δ个专用单位,以下简称单位,其中δ为无人机飞行1m,所各增加的垂直误差和水平误差,到达目的点时垂直误差和水平误差均应小于θ个单位,当垂直误差和水平误差均小于θ个单位时,无人机仍能够按照规划航迹飞行;
步骤(1.2),无人机在飞行过程中需要对定位误差进行校正,航迹规划区域内存在校正点可用于误差校正,当无人机到达校正点即能够根据校正点的误差校正类型进行误差校正,校正垂直和水平误差的位置可根据地形在航迹规划前确定,若垂直误差、水平误差都能得到及时校正,则无人机可以按照预定航线飞行,并通过若干个校正点进行误差校正后最终到达目的地;
步骤(1.3),无人机在出发点的垂直和水平误差均为0,在垂直误差校正点进行垂直误差校正后,其垂直误差将变为0,水平误差保持不变,在水平误差校正点进行水平误差校正后,其水平误差将变为0,垂直误差保持不变;
步骤(1.4),当无人机的垂直误差不大于α1个单位,水平误差不大于α2个单位时才能进行垂直误差校正,其中α1为无人机到达垂直校正点时所允许的最大垂直误差,α2为无人机到达垂直校正点时所允许的最大水平误差;
步骤(1.5),当无人机的垂直误差不大于β1个单位,水平误差不大于β2个单位时才能进行水平误差校正,其中β1为无人机到达水平校正点时所允许的最大垂直误差,β2为无人机到达水平校正点时所允许的最大水平误差;
步骤(1.6),无人机在转弯时受到结构和控制系统的限制,无法完成即时转弯,定义无人机的最小转弯半径为rmin
进一步,所述第二步基于无人机航迹规划区域,构建校正点选择策略具体包括:
步骤(2.1),经垂直/水平校正点交替校正后,无人机能按照规划路径飞行的最远距离D1应满足以下关系:
Figure GDA0003744021270000051
应在该距离内寻找新的垂直校正点进行垂直校正,否则根据第一步所述约束条件无法按规划航迹飞行;
步骤(2.2),经水平/垂直校正点交替校正后,无人机能按照规划路径飞行的最远距离D2应满足以下关系:
Figure GDA0003744021270000052
应在该距离内寻找新的水平校正点进行水平校正,否则根据第一步所述约束条件无法按规划航迹飞行;
步骤(2.3),若经垂直校正点校正后,无人机不寻找水平校正点进行校正,仍寻找垂直校正点进行校正,那么无人机能按照规划路径飞行的最远距离D3应满足以下关系:
Figure GDA0003744021270000053
若飞行超过该距离未进行校正,则在未来将无法进行垂直校正;
步骤(2.4),若经水平校正点校正后,无人机不寻找垂直校正点进行校正,仍寻找水平校正点进行校正,那么无人机能按照规划路径飞行的最远距离D4应满足以下关系:
Figure GDA0003744021270000054
若飞行超过该距离未进行校正,则在未来将无法进行水平校正;
步骤(2.5),综上所述,无人机选择垂直/水平交替校正方案可在飞行相同距离的条件下,取得最少校正次数,并且该方案在选择下一校正点时,仅需对不同组别的点进行计算即可,无需计算同组别的校正点,因此可以降低算法计算量,减少搜索次数;
进一步,所述第三步基于校正点选择策略,构建航迹规划方案具体包括代价函数构建、校正点的范围构建、循环体构建、终止条件构建和优化方案比较五部分:
步骤(3.1),采用改进的A*算法进行航迹规划,首先构建代价函数,代价函数可表示为:
f(n)=g(n)+h(n) (5)
其中,g(n)为备选校正点到出发点A的欧式距离,h(n)为备选校正点到目的点B的欧式距离,假设起点A周围存在m个同类型校正点可作为备选,将这些校正点编号为N1,...,Nm,根据式(5)可计算每个点的代价函数,可表示为:
Figure GDA0003744021270000061
其中,
Figure GDA0003744021270000062
为起点A到校正点N1之间的欧式距离,DN1,B为校正点N1到目的点B之间的欧式距离,
Figure GDA0003744021270000063
为起点A到校正点Nm之间的欧式距离,
Figure GDA0003744021270000064
为校正点Nm到目的点B之间的欧式距离;
根据式(6)的计算结果,可选出最优的校正点Ni,选择方案如下所示:
f(Ni)=min{f(N1),...,f(Nm)}i∈(1,m) (7)
步骤(3.2),基于构建的代价函数和校正点选择策略,无人机单次飞行且不进行航迹校正的飞行距离是有限的,因此本方法通过确定搜索范围、计算梯度方向和预测未来轨迹来减少备选校正点个数;
步骤(3.2.1),从出发点A出发先选择垂直校正点,无人机最远飞行距离DAmax应满足以下条件:
DAmax≤D3 (8)
其中,D3为无人机经垂直校正点校正后,不寻找水平校正点进行校正,仍寻找垂直校正点进行校正,能按照规划路径飞行的最远距离;
若DAmax>D3,即使找到垂直校正点也无法进行垂直校正,因此以出发点A为球心,以半径为D3做球,在球内所有水平校正点可作为备选校正点,则在球外的校正点则舍去;
步骤(3.2.2),在球内所有的点也非全是有效点,通过计算梯度方向来判断该点是否为有效点,点NC、ND、NE均为出发点A的待选校正点,分别计算
Figure GDA0003744021270000071
Figure GDA0003744021270000072
之间梯度,此处梯度具体为向量夹角,若夹角>90°则将该点舍去,夹角可采用如下的向量夹角公式进行计算:
Figure GDA0003744021270000073
其中,ω为向量夹角,
Figure GDA0003744021270000074
Figure GDA0003744021270000075
表示向量,具体指代为向量
Figure GDA0003744021270000076
Figure GDA0003744021270000077
即参数
Figure GDA0003744021270000078
可为向量
Figure GDA0003744021270000079
Figure GDA00037440212700000710
之一,参数
Figure GDA00037440212700000711
可为向量
Figure GDA00037440212700000712
Figure GDA00037440212700000713
Figure GDA00037440212700000714
之一,参数
Figure GDA00037440212700000715
Figure GDA00037440212700000716
不可相等,依次求解各向量之间夹角;
步骤(3.2.3),采用垂直/水平交替校正的方案,但数据集中的校正点位置是随机布置的,若校正点NC为垂直校正点,无人机经点NC垂直校正后根据本方法应选择水平校正点,但若水平校正点附近无可达的垂直校正点时,应选择择次优的水平校正点NE,并将点ND舍去;
步骤(3.2.4),通过确定搜索范围、计算梯度方向和预测未来轨迹可以减少备选校正点的个数,从而在下一步循环体构建中减少搜索次数;
步骤(3.3),对目标函数和校正点的范围完成构建后,对循环体予以构建,假设无人机从出发点A出发后先选择的校正点为垂直校正点NC,无人机到达NC时的误差为:
Figure GDA00037440212700000717
其中,
Figure GDA0003744021270000081
为无人机到达校正点NC后的垂直误差,
Figure GDA0003744021270000082
为无人机到达校正点NC后的水平误差,
Figure GDA0003744021270000083
为起点A到校正点NC之间的欧式距离;
经点NC进行垂直校正后,无人机的垂直/水平误差为:
Figure GDA0003744021270000084
其中,
Figure GDA0003744021270000085
为无人机离开校正点NC后的垂直误差,
Figure GDA0003744021270000086
为无人机离开校正点NC后的水平误差;
之后将选择水平校正点进行水平校正,待选的水平校正点处于以点NC为球心,半径为r′球内,r′满足以下约束条件:
Figure GDA0003744021270000087
再根据构建的校正点范围和目标函数,即可选出最优的下一校正点;
步骤(3.4),因采用循环体结构,需对终止条件予以构建,由第一步可知当无人机飞到目的点B时的垂直/水平误差均小于预设值θ,因此在循环体部分选出校正点后,需计算该点是否能直接到达目的点B,若可以则终止循环,否则继续循环;
设点NN为水平校正点,点NM为垂直校正点,无人机由点NM飞向点NN,此时经点NN校正后的无人机误差为:
Figure GDA0003744021270000088
其中,
Figure GDA0003744021270000089
为无人机离开校正点NN后的垂直误差,
Figure GDA00037440212700000810
为无人机离开校正点NN后的水平误差;
若无人机由点NN直接飞至目的点B,则飞至目的点B时的误差为:
Figure GDA0003744021270000091
其中,
Figure GDA0003744021270000092
为校正点NM与校正点NN之间的欧式距离,
Figure GDA0003744021270000093
为校正点NN与目的点B之间的欧式距离,
Figure GDA0003744021270000094
为无人机到达目的点B的水平误差,
Figure GDA0003744021270000095
为无人机到达目的点B的垂直误差;
Figure GDA0003744021270000096
则允许直接飞至目的点B,完成循环,否则在校正点NN继续循环,直至达到终止条件;
步骤(3.5),由于采用垂直/水平交替校正的方案,因此从出发点A出发的情况将有2种,分别是:
(a)先选择垂直校正点进行校正;
(b)先选择水平校正点进行校正;
因此将产生2种航迹规划方案,需对这2种方案进行比较,选择航迹更短的方案。
进一步,所述第四步基于航迹规划方案,构建基于动态切点调整的无人机航迹规划方法具体包括起点飞向校正点、校正点飞向校正点、校正点飞向终点三部分:
步骤(4.1),由于无人机的初始飞行方向可自行设定,因此本方法使无人机的飞行方向dirA与出发点A,校正点Na构成的向量
Figure GDA0003744021270000097
方向一致即
Figure GDA0003744021270000098
步骤(4.2),无人机经过校正点Na
Figure GDA0003744021270000099
与出发点A,校正点Na构成的向量
Figure GDA00037440212700000910
方向应具有如下关系:
Figure GDA00037440212700000911
以Na为切点做圆
Figure GDA00037440212700000912
Figure GDA00037440212700000913
半径为,
Figure GDA00037440212700000914
无人机沿圆弧飞行,飞行的圆心角ξ应满足如下关系:
Figure GDA00037440212700000915
根据圆心角可计算无人机在Na点沿圆弧飞行的航迹长度
Figure GDA0003744021270000101
为:
Figure GDA0003744021270000102
之后在N′a点飞出,飞向方向
Figure GDA0003744021270000103
应与校正点Na和校正点Nb构成的向量
Figure GDA0003744021270000104
方向一致,即
Figure GDA0003744021270000105
再沿
Figure GDA0003744021270000106
方向飞行l长度,l可由下式进行计算:
Figure GDA0003744021270000107
此时无人机的位置为点N′b,以点N′b为切点做圆
Figure GDA0003744021270000108
Figure GDA0003744021270000109
半径仍为
Figure GDA00037440212700001010
无人机沿圆
Figure GDA00037440212700001011
飞行即可到达下一校正点Nb
因此点Na飞向点Nb采用航迹平滑方案所航行的航迹长度
Figure GDA00037440212700001012
Figure GDA00037440212700001013
之间的关系如下所示:
Figure GDA00037440212700001014
在上式中
Figure GDA00037440212700001015
是一个定值,为点Na与Nb之间的欧式距离,圆心角ξ也是一个定值,由
Figure GDA00037440212700001016
Figure GDA00037440212700001017
决定,因此若要使航迹最短,即使
Figure GDA00037440212700001018
最小;因此,无人机按最小转弯半径进行转弯,即满足以下关系:
Figure GDA00037440212700001019
其中,
Figure GDA00037440212700001020
为无人机在校正点Ni时的转弯半径;
无人机此时离开点Na时的飞向方向
Figure GDA00037440212700001021
与校正点Na和校正点Nb构成的向量
Figure GDA00037440212700001022
垂直;
步骤(4.3),在圆
Figure GDA00037440212700001023
上重新寻找一个切点N″m使之满足以下条件:
Figure GDA00037440212700001024
其中Om为圆
Figure GDA00037440212700001025
的圆心,
Figure GDA00037440212700001026
为圆
Figure GDA00037440212700001027
的圆心Om与切点N″m所组成的向量;
由点Nm″直接飞向目的点B,到达终点;该段航迹长度可表示为:
Figure GDA00037440212700001028
其中,
Figure GDA0003744021270000111
为切点N″m与目的点B之间的欧式距离,φ为无人机在圆
Figure GDA0003744021270000112
处转过的圆心角,可由下式进行计算:
Figure GDA0003744021270000113
其中,
Figure GDA0003744021270000114
为圆
Figure GDA0003744021270000115
的圆心Om与无人机经过的最后一个校正点Nm所组成的向量;
因此规划后的航迹长度可由式(24)进行表示;在式(24)中,点Na为第一个校正点,点Nm为最后一个校正点,无人机共经过k个校正点;
Figure GDA0003744021270000116
其中
Figure GDA0003744021270000117
为起点A到校正点Na之间的欧式距离,
Figure GDA0003744021270000118
为校正点Ni与校正点Nj之间的欧式距离。
综上所述,本发明的优点及积极效果为:本发明采用基于动态切点调整的方法,针对现有技术存在的问题,将无人机航迹规划问题形式化为一个与定位误差校正次数、航迹长度和转弯半径相关的多目标优化问题,在保证定位误差得到有效校正的前提下优化航迹长度,并利用动态切点调整的方式来制定无人机的转弯策略。本发明可在保证无人机定位精度的前提下缩短航迹长度和校正次数,同时满足无人机的动力结构。
本发明与现有技术相比,可有效减少定位误差所带来的影响,同时满足无人机的动力结构,平滑航迹,并且本发明方法简单,易于实现,使之能实际运用于无人机通信、农业、火灾监控和检测和物联网等领域。
下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。
附图说明
图1是本发明实施例提供的一种基于动态切点调整的无人机航迹规划方法的流程图。
图2是本发明实施例提供的由出发点飞向校正点航迹平滑方案。
图3是本发明实施例提供的由校正点飞向校正点航迹平滑方案。
图4是本发明实施例提供的由校正点飞向目的点航迹平滑方案。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
针对现有技术存在的问题,本发明提供了一种基于动态切点调整的无人机航迹规划方法,下面结合附图对本发明作详细的描述。
如图1所示,本发明实施例提供的一种基于动态切点调整的无人机航迹规划方法包括以下步骤:
S101:构建由1个出发点A、1个目的地B、M个水平校正点、N个垂直校正点组成的无人机航迹规划区域;
S102:基于无人机航迹规划区域,构建校正点选择策略;
S103:基于校正点选择策略,构建航迹规划方案;
S104:基于航迹规划方案,构建基于动态切点调整的无人机航迹规划方法。
下面结合具体实施例对本发明的技术方案作进一步的描述。
本发明实施例提供的一种基于动态切点调整的无人机航迹规划方法包括以下步骤:
第一步,构建由1个出发点A、1个目的点B、R个水平校正点、L个垂直校正点组成的无人机航迹规划区域,具体包括:
步骤(1.1),构建包含2+R+L点的无人机航迹规划区域,无人机在空间飞行过程中需要实时定位,其定位误差包括垂直误差和水平误差,无人机每飞行1m,垂直误差和水平误差将各增加δ个专用单位,以下简称单位,其中δ为无人机飞行1m,所各增加的垂直误差和水平误差,到达目的点时垂直误差和水平误差均应小于θ个单位,当垂直误差和水平误差均小于θ个单位时,无人机仍能够按照规划航迹飞行;
步骤(1.2),无人机在飞行过程中需要对定位误差进行校正,航迹规划区域内存在校正点可用于误差校正,当无人机到达校正点即能够根据校正点的误差校正类型进行误差校正,校正垂直和水平误差的位置可根据地形在航迹规划前确定,若垂直误差、水平误差都能得到及时校正,则无人机可以按照预定航线飞行,并通过若干个校正点进行误差校正后最终到达目的地;
步骤(1.3),无人机在出发点的垂直和水平误差均为0,在垂直误差校正点进行垂直误差校正后,其垂直误差将变为0,水平误差保持不变,在水平误差校正点进行水平误差校正后,其水平误差将变为0,垂直误差保持不变;
步骤(1.4),当无人机的垂直误差不大于α1个单位,水平误差不大于α2个单位时才能进行垂直误差校正,其中α1为无人机到达垂直校正点时所允许的最大垂直误差,α2为无人机到达垂直校正点时所允许的最大水平误差;
步骤(1.5),当无人机的垂直误差不大于β1个单位,水平误差不大于β2个单位时才能进行水平误差校正,其中β1为无人机到达水平校正点时所允许的最大垂直误差,β2为无人机到达水平校正点时所允许的最大水平误差;
步骤(1.6),无人机在转弯时受到结构和控制系统的限制,无法完成即时转弯,定义无人机的最小转弯半径为rmin
第二步,基于无人机航迹规划区域,构建校正点选择策略具体包括:
步骤(2.1),经垂直/水平校正点交替校正后,无人机能按照规划路径飞行的最远距离D1应满足以下关系:
Figure GDA0003744021270000131
应在该距离内寻找新的垂直校正点进行垂直校正,否则根据第一步所述约束条件无法按规划航迹飞行;
步骤(2.2),经水平/垂直校正点交替校正后,无人机能按照规划路径飞行的最远距离D2应满足以下关系:
Figure GDA0003744021270000141
应在该距离内寻找新的水平校正点进行水平校正,否则根据第一步所述约束条件无法按规划航迹飞行;
步骤(2.3),若经垂直校正点校正后,无人机不寻找水平校正点进行校正,仍寻找垂直校正点进行校正,那么无人机能按照规划路径飞行的最远距离D3应满足以下关系:
Figure GDA0003744021270000142
若飞行超过该距离未进行校正,则在未来将无法进行垂直校正;
步骤(2.4),若经水平校正点校正后,无人机不寻找垂直校正点进行校正,仍寻找水平校正点进行校正,那么无人机能按照规划路径飞行的最远距离D4应满足以下关系:
Figure GDA0003744021270000143
若飞行超过该距离未进行校正,则在未来将无法进行水平校正;
步骤(2.5),综上所述,无人机选择垂直/水平交替校正方案可在飞行相同距离的条件下,取得最少校正次数,并且该方案在选择下一校正点时,仅需对不同组别的点进行计算即可,无需计算同组别的校正点,因此可以降低算法计算量,减少搜索次数;
第三步,基于校正点选择策略,构建航迹规划方案具体包括代价函数构建、校正点的范围构建、循环体构建、终止条件构建和优化方案比较五部分:
步骤(3.1),采用改进的A*算法进行航迹规划,首先构建代价函数,代价函数可表示为:
f(n)=g(n)+h(n) (5)
其中,g(n)为备选校正点到出发点A的欧式距离,h(n)为备选校正点到目的点B的欧式距离,假设起点A周围存在m个同类型校正点可作为备选,将这些校正点编号为N1,...,Nm,根据式(5)可计算每个点的代价函数,可表示为:
Figure GDA0003744021270000151
其中,
Figure GDA0003744021270000152
为起点A到校正点N1之间的欧式距离,
Figure GDA0003744021270000157
为校正点N1到目的点B之间的欧式距离,
Figure GDA0003744021270000153
为起点A到校正点Nm之间的欧式距离,
Figure GDA0003744021270000154
为校正点Nm到目的点B之间的欧式距离;
根据式(6)的计算结果,可选出最优的校正点Ni,选择方案如下所示:
f(Ni)=min{f(N1),...,f(Nm)}i∈(1,m) (7)
步骤(3.2),基于构建的代价函数和校正点选择策略,无人机单次飞行且不进行航迹校正的飞行距离是有限的,因此本方法通过确定搜索范围、计算梯度方向和预测未来轨迹来减少备选校正点个数;
步骤(3.2.1),从出发点A出发先选择垂直校正点,无人机最远飞行距离DAmax应满足以下条件:
DAmax≤D3 (8)
其中,D3为无人机经垂直校正点校正后,不寻找水平校正点进行校正,仍寻找垂直校正点进行校正,能按照规划路径飞行的最远距离;
若DAmax>D3,即使找到垂直校正点也无法进行垂直校正,因此以出发点A为球心,以半径为D3做球,在球内所有水平校正点可作为备选校正点,则在球外的校正点则舍去;
步骤(3.2.2),在球内所有的点也非全是有效点,通过计算梯度方向来判断该点是否为有效点,点NC、ND、NE均为出发点A的待选校正点,分别计算
Figure GDA0003744021270000155
Figure GDA0003744021270000156
之间梯度,此处梯度具体为向量夹角,若夹角>90°则将该点舍去,夹角可采用如下的向量夹角公式进行计算:
Figure GDA0003744021270000161
其中,ω为向量夹角,
Figure GDA0003744021270000162
Figure GDA0003744021270000163
表示向量,具体指代为向量
Figure GDA0003744021270000164
Figure GDA0003744021270000165
即参数
Figure GDA0003744021270000166
可为向量
Figure GDA0003744021270000167
Figure GDA0003744021270000168
之一,参数
Figure GDA0003744021270000169
可为向量
Figure GDA00037440212700001610
Figure GDA00037440212700001611
Figure GDA00037440212700001612
之一,参数
Figure GDA00037440212700001613
Figure GDA00037440212700001614
不可相等,依次求解各向量之间夹角;
步骤(3.2.3),采用垂直/水平交替校正的方案,但数据集中的校正点位置是随机布置的,若校正点NC为垂直校正点,无人机经点NC垂直校正后根据本方法应选择水平校正点,但若水平校正点附近无可达的垂直校正点时,应选择择次优的水平校正点NE,并将点ND舍去;
步骤(3.2.4),通过确定搜索范围、计算梯度方向和预测未来轨迹可以减少备选校正点的个数,从而在下一步循环体构建中减少搜索次数;
步骤(3.3),对目标函数和校正点的范围完成构建后,对循环体予以构建,假设无人机从出发点A出发后先选择的校正点为垂直校正点NC,无人机到达NC时的误差为:
Figure GDA00037440212700001615
其中,
Figure GDA00037440212700001616
为无人机到达校正点NC后的垂直误差,
Figure GDA00037440212700001617
为无人机到达校正点NC后的水平误差,
Figure GDA00037440212700001618
为起点A到校正点NC之间的欧式距离;
经点NC进行垂直校正后,无人机的垂直/水平误差为:
Figure GDA00037440212700001619
其中,
Figure GDA00037440212700001620
为无人机离开校正点NC后的垂直误差,
Figure GDA00037440212700001621
为无人机离开校正点NC后的水平误差;
之后将选择水平校正点进行水平校正,待选的水平校正点处于以点NC为球心,半径为r′球内,r′满足以下约束条件:
Figure GDA0003744021270000171
再根据构建的校正点范围和目标函数,即可选出最优的下一校正点;
步骤(3.4),因采用循环体结构,需对终止条件予以构建,由第一步可知当无人机飞到目的点B时的垂直/水平误差均小于预设值θ,因此在循环体部分选出校正点后,需计算该点是否能直接到达目的点B,若可以则终止循环,否则继续循环;
设点NN为水平校正点,点NM为垂直校正点,无人机由点NM飞向点NN,此时经点NN校正后的无人机误差为:
Figure GDA0003744021270000172
其中,
Figure GDA0003744021270000173
为无人机离开校正点NN后的垂直误差,
Figure GDA0003744021270000174
为无人机离开校正点NN后的水平误差;
若无人机由点NN直接飞至目的点B,则飞至目的点B时的误差为:
Figure GDA0003744021270000175
其中,
Figure GDA0003744021270000176
为校正点NM与校正点NN之间的欧式距离,
Figure GDA0003744021270000177
为校正点NN与目的点B之间的欧式距离,
Figure GDA0003744021270000178
为无人机到达目的点B的水平误差,
Figure GDA0003744021270000179
为无人机到达目的点B的垂直误差;
Figure GDA00037440212700001710
则允许直接飞至目的点B,完成循环,否则在校正点NN继续循环,直至达到终止条件;
步骤(3.5),由于采用垂直/水平交替校正的方案,因此从出发点A出发的情况将有2种,分别是:
(a)先选择垂直校正点进行校正;
(b)先选择水平校正点进行校正;
因此将产生2种航迹规划方案,需对这2种方案进行比较,选择航迹更短的方案;
第四步,基于航迹规划方案,构建基于动态切点调整的无人机航迹规划方法具体包括起点飞向校正点、校正点飞向校正点、校正点飞向终点三部分:
步骤(4.1),由于无人机的初始飞行方向可自行设定,因此本方法使无人机的飞行方向dirA与出发点A,校正点Na构成的向量
Figure GDA0003744021270000181
方向一致即
Figure GDA0003744021270000182
步骤(4.2),无人机经过校正点Na
Figure GDA00037440212700001822
与出发点A,校正点Na构成的向量
Figure GDA0003744021270000183
方向应具有如下关系:
Figure GDA0003744021270000184
以Na为切点做圆
Figure GDA0003744021270000185
Figure GDA0003744021270000186
半径为,
Figure GDA0003744021270000187
无人机沿圆弧飞行,飞行的圆心角ξ应满足如下关系:
Figure GDA0003744021270000188
根据圆心角可计算无人机在Na点沿圆弧飞行的航迹长度
Figure GDA0003744021270000189
为:
Figure GDA00037440212700001810
之后在N′a点飞出,飞向方向
Figure GDA00037440212700001811
应与校正点Na和校正点Nb构成的向量
Figure GDA00037440212700001812
方向一致,即
Figure GDA00037440212700001813
再沿
Figure GDA00037440212700001814
方向飞行l长度,l可由下式进行计算:
Figure GDA00037440212700001815
此时无人机的位置为点N′b,以点N′b为切点做圆
Figure GDA00037440212700001816
Figure GDA00037440212700001817
半径仍为
Figure GDA00037440212700001818
无人机沿圆
Figure GDA00037440212700001819
飞行即可到达下一校正点Nb
因此点Na飞向点Nb采用航迹平滑方案所航行的航迹长度
Figure GDA00037440212700001820
Figure GDA00037440212700001821
之间的关系如下所示:
Figure GDA0003744021270000191
在上式中
Figure GDA0003744021270000192
是一个定值,为点Na与Nb之间的欧式距离,圆心角ξ也是一个定值,由
Figure GDA0003744021270000193
Figure GDA0003744021270000194
决定,因此若要使航迹最短,即使
Figure GDA0003744021270000195
最小;因此,无人机按最小转弯半径进行转弯,即满足以下关系:
Figure GDA0003744021270000196
其中,
Figure GDA0003744021270000197
为无人机在校正点Ni时的转弯半径;
无人机此时离开点Na时的飞向方向
Figure GDA0003744021270000198
与校正点Na和校正点Nb构成的向量
Figure GDA0003744021270000199
垂直;
步骤(4.3),在圆
Figure GDA00037440212700001910
上重新寻找一个切点N″m使之满足以下条件:
Figure GDA00037440212700001911
其中Om为圆
Figure GDA00037440212700001912
的圆心,
Figure GDA00037440212700001913
为圆
Figure GDA00037440212700001914
的圆心Om与切点N″m所组成的向量;
由点Nm″直接飞向目的点B,到达终点;该段航迹长度可表示为:
Figure GDA00037440212700001915
其中,
Figure GDA00037440212700001916
为切点N″m与目的点B之间的欧式距离,φ为无人机在圆
Figure GDA00037440212700001917
处转过的圆心角,可由下式进行计算:
Figure GDA00037440212700001918
其中,
Figure GDA00037440212700001921
为圆
Figure GDA00037440212700001919
的圆心Om与无人机经过的最后一个校正点Nm所组成的向量;
因此规划后的航迹长度可由式(24)进行表示;在式(24)中,点Na为第一个校正点,点Nm为最后一个校正点,无人机共经过k个校正点;
Figure GDA00037440212700001920
其中
Figure GDA0003744021270000201
为起点A到校正点Na之间的欧式距离,
Figure GDA0003744021270000202
为校正点Ni与校正点Nj之间的欧式距离。
综上,本发明实施例提供的一种基于动态切点调整的无人机航迹规划方法,针对现有技术存在的问题,将航迹规划问题形式化为一个与定位误差校正次数、航迹长度和转弯半径相关的多目标优化问题,在保证定位误差得到有效校正的前提下优化航迹长度,并利用动态切点调整的方式来制定无人机的转弯策略。本发明可有效降低无人机航迹长度与算法复杂度并能符合无人机动力结构,相较于其他无人机航迹规划方法更能适用于实际场景应用,实现多约束条件下的无人机航迹快速规划。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (1)

1.一种基于动态切点调整的无人机航迹规划方法,其特征在于,所述一种基于动态切点调整的无人机航迹规划方法将航迹规划问题形式化为一个与定位误差校正次数、航迹长度和转弯半径相关的多目标优化问题,在保证定位误差得到有效校正的前提下优化航迹长度,并利用动态切点调整的方式来制定无人机的转弯策略;本方法可有效降低无人机航迹长度与算法复杂度并能符合无人机动力结构,相较于其他无人机航迹规划方法更能适用于实际场景应用,实现多约束条件下的无人机航迹快速规划;
所述一种基于动态切点调整的无人机航迹规划方法包括以下步骤:
第一步,构建由1个出发点A、1个目的地B、M个水平校正点、N个垂直校正点组成的无人机航迹规划区域;
第二步,基于无人机航迹规划区域,构建校正点选择策略;
第三步,基于校正点选择策略,构建航迹规划方案;
第四步,基于航迹规划方案,构建基于动态切点调整的无人机航迹规划方法;
所述第一步构建由1个出发点A、1个目的点B、R个水平校正点、L个垂直校正点组成的无人机航迹规划区域,具体包括:
步骤(1.1),构建包含2+R+L点的无人机航迹规划区域,无人机在空间飞行过程中需要实时定位,其定位误差包括垂直误差和水平误差,无人机每飞行1m,垂直误差和水平误差将各增加δ个专用单位,以下简称单位,其中δ为无人机飞行1m,所各增加的垂直误差和水平误差,到达目的点时垂直误差和水平误差均应小于等于θ个单位,当垂直误差和水平误差均小于等于θ个单位时,无人机仍能够按照规划航迹飞行;
步骤(1.2),无人机在飞行过程中需要对定位误差进行校正,航迹规划区域内存在校正点可用于误差校正,当无人机到达校正点即能够根据校正点的误差校正类型进行误差校正,校正垂直和水平误差的位置可根据地形在航迹规划前确定,若垂直误差、水平误差都能得到及时校正,则无人机可以按照预定航线飞行,并通过若干个校正点进行误差校正后最终到达目的地;
步骤(1.3),无人机在出发点的垂直和水平误差均为0,在垂直误差校正点进行垂直误差校正后,其垂直误差将变为0,水平误差保持不变,在水平误差校正点进行水平误差校正后,其水平误差将变为0,垂直误差保持不变;
步骤(1.4),当无人机的垂直误差不大于α1个单位,水平误差不大于α2个单位时才能进行垂直误差校正,其中α1为无人机到达垂直校正点时所允许的最大垂直误差,α2为无人机到达垂直校正点时所允许的最大水平误差;
步骤(1.5),当无人机的垂直误差不大于β1个单位,水平误差不大于β2个单位时才能进行水平误差校正,其中β1为无人机到达水平校正点时所允许的最大垂直误差,β2为无人机到达水平校正点时所允许的最大水平误差;
步骤(1.6),无人机在转弯时受到结构和控制系统的限制,无法完成即时转弯,定义无人机的最小转弯半径为rmin=200m;
所述第二步基于无人机航迹规划区域,构建校正点选择策略具体包括:
步骤(2.1),经垂直/水平校正点交替校正后,无人机能按照规划路径飞行的最远距离D1应满足以下关系:
Figure FDA0003797168690000021
应在该距离内寻找新的垂直校正点进行垂直校正,否则根据第一步约束条件无法按规划航迹飞行;
步骤(2.2),经水平/垂直校正点交替校正后,无人机能按照规划路径飞行的最远距离D2应满足以下关系:
Figure FDA0003797168690000022
应在该距离内寻找新的水平校正点进行水平校正,否则根据第一步约束条件无法按规划航迹飞行;
步骤(2.3),若经垂直校正点校正后,无人机不寻找水平校正点进行校正,仍寻找垂直校正点进行校正,那么无人机能按照规划路径飞行的最远距离D3应满足以下关系:
Figure FDA0003797168690000031
若飞行超过该距离未进行校正,则在未来将无法进行垂直校正;
步骤(2.4),若经水平校正点校正后,无人机不寻找垂直校正点进行校正,仍寻找水平校正点进行校正,那么无人机能按照规划路径飞行的最远距离D4应满足以下关系:
Figure FDA0003797168690000032
若飞行超过该距离未进行校正,则在未来将无法进行水平校正;
步骤(2.5),综上所述,无人机选择垂直/水平交替校正方案可在飞行相同距离的条件下,取得最少校正次数,并且该方案在选择下一校正点时,仅需对不同组别的点进行计算即可,无需计算同组别的校正点,因此可以降低算法计算量,减少搜索次数;
所述第三步基于校正点选择策略,构建航迹规划方案具体包括代价函数构建、校正点的范围构建、循环体构建、终止条件构建和优化方案比较五部分:
步骤(3.1),采用改进的A*算法进行航迹规划,首先构建代价函数,代价函数可表示为:
f(n)=g(n)+h(n) (5)
其中,g(n)为备选校正点到出发点A的欧式距离,h(n)为备选校正点到目的点B的欧式距离,假设起点A周围存在m个同类型校正点可作为备选,将这些校正点编号为N1,...,Nm,根据式(5)可计算每个点的代价函数,可表示为:
Figure FDA0003797168690000033
其中,
Figure FDA0003797168690000041
为起点A到校正点N1之间的欧式距离,
Figure FDA0003797168690000042
为校正点N1到目的点B之间的欧式距离,
Figure FDA0003797168690000043
为起点A到校正点Nm之间的欧式距离,
Figure FDA0003797168690000044
为校正点Nm到目的点B之间的欧式距离;
根据式(6)的计算结果,可选出最优的校正点Ni,选择方案如下所示:
f(Ni)=min{f(N1),...,f(Nm)}i∈(1,m) (7)
步骤(3.2),基于构建的代价函数和校正点选择策略,无人机单次飞行且不进行航迹校正的飞行距离是有限的,因此本方法通过确定搜索范围、计算梯度方向和预测未来轨迹来减少备选校正点个数;
步骤(3.2.1),从出发点A出发先选择垂直校正点,无人机最远飞行距离DAmax应满足以下条件:
DAmax≤D3 (8)
其中,D3为无人机经垂直校正点校正后,不寻找水平校正点进行校正,仍寻找垂直校正点进行校正,能按照规划路径飞行的最远距离;
若DAmax>D3,即使找到垂直校正点也无法进行垂直校正,因此以出发点A为球心,以半径为D3做球,在球内所有水平校正点可作为备选校正点,则在球外的校正点则舍去;
步骤(3.2.2),在球内所有的点也非全是有效点,通过计算梯度方向来判断该点是否为有效点,点NC、ND、NE均为出发点A的待选校正点,分别计算
Figure FDA0003797168690000045
Figure FDA0003797168690000046
之间梯度,此处梯度具体为向量夹角,若夹角>90°则将该点舍去,夹角可采用如下的向量夹角公式进行计算:
Figure FDA0003797168690000047
其中,ω为向量夹角,
Figure FDA0003797168690000048
Figure FDA0003797168690000049
表示向量,具体指代为向量
Figure FDA00037971686900000410
Figure FDA00037971686900000411
即参数
Figure FDA00037971686900000412
可为向量
Figure FDA00037971686900000413
Figure FDA00037971686900000414
之一,参数
Figure FDA00037971686900000415
可为向量
Figure FDA00037971686900000416
Figure FDA00037971686900000417
Figure FDA00037971686900000418
之一,参数
Figure FDA00037971686900000419
Figure FDA00037971686900000420
不可相等,依次求解各向量之间夹角;
步骤(3.2.3),采用垂直/水平交替校正的方案,但数据集中的校正点位置是随机布置的,若校正点NC为垂直校正点,无人机经点NC垂直校正后根据本方法应选择水平校正点,但若水平校正点附近无可达的垂直校正点时,应选择择次优的水平校正点NE,并将点ND舍去;
步骤(3.2.4),通过确定搜索范围、计算梯度方向和预测未来轨迹可以减少备选校正点的个数,从而在下一步循环体构建中减少搜索次数;
步骤(3.3),对目标函数和校正点的范围完成构建后,对循环体予以构建,假设无人机从出发点A出发后先选择的校正点为垂直校正点NC,无人机到达NC时的误差为:
Figure FDA0003797168690000051
其中,
Figure FDA0003797168690000052
为无人机到达校正点NC后的垂直误差,
Figure FDA0003797168690000053
为无人机到达校正点NC后的水平误差,
Figure FDA0003797168690000054
为起点A到校正点NC之间的欧式距离;
经点NC进行垂直校正后,无人机的垂直/水平误差为:
Figure FDA0003797168690000055
其中,
Figure FDA0003797168690000056
为无人机离开校正点NC后的垂直误差,
Figure FDA0003797168690000057
为无人机离开校正点NC后的水平误差;
之后将选择水平校正点进行水平校正,待选的水平校正点处于以点NC为球心,半径为r′球内,r′满足以下约束条件:
Figure FDA0003797168690000058
再根据构建的校正点范围和目标函数,即可选出最优的下一校正点;
步骤(3.4),因采用循环体结构,需对终止条件予以构建,由第一步可知当无人机飞到目的点B时的垂直/水平误差均小于预设值θ,因此在循环体部分选出校正点后,需计算该点是否能直接到达目的点B,若可以则终止循环,否则继续循环;
设点NN为水平校正点,点NM为垂直校正点,无人机由点NM飞向点NN,此时经点NN校正后的无人机误差为:
Figure FDA0003797168690000061
其中,
Figure FDA0003797168690000062
为无人机离开校正点NN后的垂直误差,
Figure FDA0003797168690000063
为无人机离开校正点NN后的水平误差;
若无人机由点NN直接飞至目的点B,则飞至目的点B时的误差为:
Figure FDA0003797168690000064
其中,
Figure FDA0003797168690000065
为校正点NM与校正点NN之间的欧式距离,
Figure FDA0003797168690000066
为校正点NN与目的点B之间的欧式距离,
Figure FDA0003797168690000067
为无人机到达目的点B的水平误差,
Figure FDA0003797168690000068
为无人机到达目的点B的垂直误差;
Figure FDA0003797168690000069
则允许直接飞至目的点B,完成循环,否则在校正点NN继续循环,直至达到终止条件;
步骤(3.5),由于采用垂直/水平交替校正的方案,因此从出发点A出发的情况将有2种,分别是:
(a)先选择垂直校正点进行校正;
(b)先选择水平校正点进行校正;
因此将产生2种航迹规划方案,需对这2种方案进行比较,选择航迹更短的方案;
所述第四步基于航迹规划方案,构建基于动态切点调整的无人机航迹规划方法具体包括起点飞向校正点、校正点飞向校正点、校正点飞向终点三部分:
步骤(4.1),由于无人机的初始飞行方向可自行设定,因此本方法使无人机的飞行方向dirA与出发点A,校正点Na构成的向量
Figure FDA0003797168690000071
方向一致即
Figure FDA0003797168690000072
步骤(4.2),无人机经过校正点Na
Figure FDA0003797168690000073
与出发点A、校正点Na构成的向量
Figure FDA0003797168690000074
方向应具有如下关系:
Figure FDA0003797168690000075
以Na为切点做圆
Figure FDA0003797168690000076
Figure FDA0003797168690000077
半径为,
Figure FDA0003797168690000078
无人机沿圆弧飞行,飞行的圆心角ξ应满足如下关系:
Figure FDA0003797168690000079
根据圆心角可计算无人机在Na点沿圆弧飞行的航迹长度
Figure FDA00037971686900000710
为:
Figure FDA00037971686900000711
之后在N′a点飞出,飞向方向
Figure FDA00037971686900000712
应与校正点Na和校正点Nb构成的向量
Figure FDA00037971686900000713
方向一致,即
Figure FDA00037971686900000714
再沿
Figure FDA00037971686900000715
方向飞行l长度,l可由下式进行计算:
Figure FDA00037971686900000716
此时无人机的位置为点N′b,以点N′b为切点做圆
Figure FDA00037971686900000717
Figure FDA00037971686900000718
半径仍为
Figure FDA00037971686900000719
无人机沿圆
Figure FDA00037971686900000720
飞行即可到达下一校正点Nb
因此点Na飞向点Nb采用航迹平滑方案所航行的航迹长度
Figure FDA00037971686900000721
Figure FDA00037971686900000722
之间的关系如下所示:
Figure FDA00037971686900000723
在式(19)中,
Figure FDA00037971686900000724
中ξ为弧度,
Figure FDA00037971686900000725
中ξ为角度;
在上式中
Figure FDA00037971686900000726
是一个定值,为点Na与Nb之间的欧式距离,圆心角ξ也是一个定值,由
Figure FDA00037971686900000727
Figure FDA00037971686900000728
决定,因此若要使航迹最短,即使
Figure FDA00037971686900000729
最小;因此,无人机按最小转弯半径进行转弯,即满足以下关系:
Figure FDA0003797168690000081
其中,
Figure FDA0003797168690000082
为无人机在校正点Ni时的转弯半径;
无人机此时离开点Na时的飞向方向
Figure FDA0003797168690000083
与校正点Na和校正点Nb构成的向量
Figure FDA0003797168690000084
垂直;
步骤(4.3),在圆
Figure FDA0003797168690000085
上重新寻找一个切点N″m使之满足以下条件:
Figure FDA0003797168690000086
其中Om为圆
Figure FDA0003797168690000087
的圆心,
Figure FDA0003797168690000088
为圆
Figure FDA0003797168690000089
的圆心Om与切点N″m所组成的向量;
由点Nm″直接飞向目的点B,到达终点;该段航迹长度可表示为:
Figure FDA00037971686900000810
其中,
Figure FDA00037971686900000811
为切点N″m与目的点B之间的欧式距离,φ为无人机在圆
Figure FDA00037971686900000812
处转过的圆心角,可由下式进行计算:
Figure FDA00037971686900000813
其中,
Figure FDA00037971686900000814
为圆
Figure FDA00037971686900000815
的圆心Om与无人机经过的最后一个校正点Nm所组成的向量;
因此规划后的航迹长度可由式(24)进行表示;在式(24)中,点Na为第一个校正点,点Nm为最后一个校正点,无人机共经过k个校正点;
Figure FDA00037971686900000816
其中
Figure FDA00037971686900000817
为起点A到校正点Na之间的欧式距离,
Figure FDA00037971686900000818
为校正点Ni与校正点Nj之间的欧式距离。
CN202010050196.0A 2020-01-17 2020-01-17 一种基于动态切点调整的无人机航迹规划方法 Active CN111207752B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010050196.0A CN111207752B (zh) 2020-01-17 2020-01-17 一种基于动态切点调整的无人机航迹规划方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010050196.0A CN111207752B (zh) 2020-01-17 2020-01-17 一种基于动态切点调整的无人机航迹规划方法

Publications (2)

Publication Number Publication Date
CN111207752A CN111207752A (zh) 2020-05-29
CN111207752B true CN111207752B (zh) 2022-10-04

Family

ID=70784849

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010050196.0A Active CN111207752B (zh) 2020-01-17 2020-01-17 一种基于动态切点调整的无人机航迹规划方法

Country Status (1)

Country Link
CN (1) CN111207752B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111928838B (zh) * 2020-07-06 2021-10-22 西安羚控电子科技有限公司 一种基于舰载无人机前轮转弯技术的路径规划方法
CN112066991A (zh) * 2020-09-07 2020-12-11 常州微亿智造科技有限公司 多约束条件下智能飞行器航迹快速规划方法、装置
CN111998858B (zh) * 2020-09-15 2024-01-19 长春工业大学 一种基于改进a*算法的无人机航路规划方法
CN112346482B (zh) * 2020-11-25 2023-03-03 中国工程物理研究院总体工程研究所 飞行航线管理方法
CN112729308B (zh) * 2020-12-24 2024-05-03 广州融赋数智技术服务有限公司 多约束条件下无人机快速航迹规划方法
CN113904743B (zh) * 2021-10-11 2023-11-21 重庆邮电大学 一种无人机中继系统的安全通信资源优化设计方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110262548A (zh) * 2019-06-21 2019-09-20 北京理工大学 一种考虑抵达时间约束的无人机航迹规划方法
CN110320933A (zh) * 2019-07-29 2019-10-11 南京航空航天大学 一种巡航任务下无人机避障运动规划方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107504972B (zh) * 2017-07-27 2018-08-07 北京航空航天大学 一种基于鸽群算法的飞行器航迹规划方法及装置
CN107607120B (zh) * 2017-09-06 2020-07-07 北京理工大学 基于改进修复式Anytime稀疏A*算法的无人机动态航迹规划方法
CN108592925B (zh) * 2018-07-26 2020-09-08 中国人民解放军陆军工程大学 基于最小转弯半径的无人机转弯航迹规划算法
CN109947136B (zh) * 2019-02-28 2020-07-10 清华大学 一种面向无人机群体快速目标搜索的协同主动感知方法
CN110617818A (zh) * 2019-10-15 2019-12-27 中国人民解放军军事科学院国防科技创新研究院 一种无人机航迹生成方法
CN110658850B (zh) * 2019-11-12 2022-07-12 重庆大学 一种基于贪心策略的无人飞行器的航迹规划方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110262548A (zh) * 2019-06-21 2019-09-20 北京理工大学 一种考虑抵达时间约束的无人机航迹规划方法
CN110320933A (zh) * 2019-07-29 2019-10-11 南京航空航天大学 一种巡航任务下无人机避障运动规划方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
满足直飞限制二维航迹规划方法研究;涂吉林等;《电子学报》;20000325(第03期);全文 *
突发威胁下的无人机航迹规划算法;孙健等;《飞行力学》;20180301(第03期);全文 *

Also Published As

Publication number Publication date
CN111207752A (zh) 2020-05-29

Similar Documents

Publication Publication Date Title
CN111207752B (zh) 一种基于动态切点调整的无人机航迹规划方法
US11727812B2 (en) Airplane flight path planning method and device based on the pigeon-inspired optimization
CN111595343B (zh) 一种基于定位误差校正的无人机航迹规划方法
CN110470301B (zh) 多动态任务目标点下的无人机路径规划方法
Wu et al. Bi-directional adaptive A* algorithm toward optimal path planning for large-scale UAV under multi-constraints
CN111024092B (zh) 一种多约束条件下智能飞行器航迹快速规划方法
CN109947136B (zh) 一种面向无人机群体快速目标搜索的协同主动感知方法
CN110715664B (zh) 多约束条件下智能无人飞行器航迹快速规划方法
CN111580556A (zh) 一种时空约束下多无人机协同路径规划与制导方法
CN111006693B (zh) 智能飞行器航迹规划系统及其方法
CN103700286A (zh) 一种舰载无人机自动着舰引导方法
CN113311869B (zh) 一种多约束条件下无人机路径规划方法
CN112363539B (zh) 多无人机协同目标搜索方法
CN110849355B (zh) 一种地磁多参量多目标快速收敛的仿生导航方法
CN113342030B (zh) 基于强化学习的多无人机协同自组织控制方法及系统
CN110032214B (zh) 基于矢量场的快速Standoff目标跟踪方法
CN111273686A (zh) 一种三维环境下多无人机同时到达指定地的路径规划方法
CN116088576A (zh) 一种基于改进鲸鱼算法的无人机三维路径规划方法
CN111897362A (zh) 一种复杂环境下的翼伞组合式航迹规划方法
CN115494866A (zh) 一种多无人机全局及局部路径智能规划方法及系统
CN116880561A (zh) 基于改进粒子群无人机路径规划安全增强的优化方法及系统
CN110794874A (zh) 一种定位误差约束下的飞行器航迹快速规划方法
CN117724524A (zh) 一种基于改进球面向量粒子群算法的无人机航线规划方法
Liu et al. Improved ant colony algorithm for path planning of mobile robots based on compound prediction mechanism
CN115933736A (zh) 一种航迹快速规划精准控制飞行器的方法

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