CN112880688B - 基于混沌自适应麻雀搜索算法的无人机三维航迹规划方法 - Google Patents
基于混沌自适应麻雀搜索算法的无人机三维航迹规划方法 Download PDFInfo
- Publication number
- CN112880688B CN112880688B CN202110110630.4A CN202110110630A CN112880688B CN 112880688 B CN112880688 B CN 112880688B CN 202110110630 A CN202110110630 A CN 202110110630A CN 112880688 B CN112880688 B CN 112880688B
- Authority
- CN
- China
- Prior art keywords
- sparrow
- algorithm
- track
- adaptive
- aerial vehicle
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/20—Instruments for performing navigational calculations
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine 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)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了基于混沌自适应麻雀搜索算法的无人机三维航迹规划方法,包括以下步骤:根据飞行环境建立飞行环境模型;建立无人机飞行代价函数,对无人机航迹性能进行评价;采用混沌初始化种群策略,自适应权重策略和柯西‑高斯混合变异策略对麻雀搜索算法进行改进,提出混沌自适应麻雀搜索算法;采用混沌自适应麻雀搜索算法对三维环境中的无人机航迹进行规划,得到无人机航迹规划最优解,得到规划结果;提出改进麻雀算法在求解质量上有着明显的优势,其混沌策略和自适应策略让算法拥有快速的收敛速度,和优秀的收敛精度,变异策略则让算法拥有很强的跳出局部最优的能力,进而可快速得到一条优秀的无人机飞行路线。
Description
技术领域
本发明涉及无人机航迹优化的研究领域,特别涉及基于混沌自适应麻雀搜索算法的无人机三维航迹规划方法。
背景技术
众所周知,无人机航迹规划已成为定义无人机任务的重要内容之一。无人机飞行任务所需的三维路径规划可以定义为非确定性多项式(NP)问题,其主要目的是对出发点到目标点之间的路径进行优化。在此过程中需要处理多个约束条件,如无人机自身的性能约束、威胁约束和环境约束等等。为了使得规划出的航迹满足实际可飞的要求,还必须对航迹进行平滑处理。综上所述,无人机三维航迹规化问题可以被认为是一个多约束优化问题。
随着规划问题的复杂程度的不断增加,对应求解的难度和计算量也迅速的增长。根据没有免费的午餐定理(NFL)可知,一个优化算法可能在一系列问题中表现良好,而在另一系列问题中表现较差。因此现有的算法之中很难找到一种可以有效应对所有复杂环境的无人机航迹规划方法。
麻雀搜索算法(Sparrow Search Algorithm,SSA)是一种新型的群智能优化算法。该算法是通过麻雀的觅食行为和反捕食行为的启发而提出的。标准SSA算法具有调整参数少,收敛速度快,计算简单等优点。但在缺点方面,SSA算法在求解复杂工程优化问题时,会开始出现“早熟”现象导致收敛精度不高,易于陷入局部最优,因此以标准SSA算法直接求解航迹规划问题时得到的航迹可飞性不高。
发明内容
本发明的主要目的在于克服现有技术的缺点与不足,提供基于混沌自适应麻雀搜索算法的无人机三维航迹规划方法,首先考虑避障约束与航迹的长度代价,高度代价,威胁代价,路径平滑代价,建立三维航迹规划模型,然后提出一种改进的混沌自适应麻雀搜索算法,并应用于求解上述模型,得到规划结果。
本发明的目的通过以下的技术方案实现:
基于混沌自适应麻雀搜索算法的无人机三维航迹规划方法,其特征在于,包括以下步骤:
根据飞行环境建立飞行环境模型;
建立无人机飞行代价函数,对无人机航迹性能进行评价;
采用混沌初始化种群策略,自适应权重策略和柯西-高斯混合变异策略对麻雀搜索算法进行改进,提出混沌自适应麻雀搜索算法;
采用混沌自适应麻雀搜索算法对三维环境中的无人机航迹进行规划,得到无人机航迹规划最优解,得到规划结果。
进一步地,所述根据飞行环境建立飞行环境模型,具体为:采集目标飞行区域的数字高程模型,通过MATLAB处理,得到三维地图模型。
进一步地,所述三维地图模型包括威胁区域,在三维地图模型中设定飞行起始点和目标坐标点以及威胁分布坐标。
进一步地,所述建立三维航迹规划模型,即建立无人机飞行代价函数,对无人机航迹性能进行评价,具体为:
设定无人机三维航迹规划的性能指标:航迹长度,飞行高度,最大转角;
通过航迹长度代价函数获取航迹长度;通过飞行高度模型获取飞行高度;通过转角代价模型获取最大转角;
通过航迹长度代价函数、飞行高度模型、转角代价模型,得到航迹代价函数:
Jcost=w1Lpath+w2Hheight+w3Jturn,
其中,Jturn是总的航迹代价函数,Lpath为航迹长度代价函数,hheight为高度的标准差代价函数,Jturn是转角代价函数,参数wi,i=1,2,3表示每个代价函数的权重,且满足以下条件:
通过对代价函数的求解处理,得到由线段组成的航迹,再使用插值的方式调用B样条曲线函数,对已获取的航迹再进行平滑处理,最终得到无人机实际可飞的航迹。
进一步地,所述通过航迹长度代价函数获取航迹长度,如下所示:设定完整的航迹有n个节点,第i个航迹点和第i+1个航迹点之间的距离为li,两个航迹点的坐标分别表示为g(i)=(xi,yi,zi),g(i+1)=(xi+1,yi+1,zi+1),航迹需要满足以下条件:
其中,Lpath为航迹长度代价函数,说
所述通过飞行高度模型获取飞行高度,如下所示:飞行高度模型为:
所述通过转角代价模型获取最大转角,如下所示:转角代价模型如下:
其中,Jturn是转角代价函数,φ为最大转角,θ是当前转角,ai是第i段航迹向量,|ai|表示矢量a的长度。
进一步地,所述提出改进的混沌自适应麻雀搜索算法,在算法初始化阶段采用混沌初始化种群策略增强算法稳定性,在发现者位置更新阶段引入了自适应权重因子,对精英麻雀个体采用柯西-高斯混合变异策略;如下所示:
采用标准SSA算法,将麻雀群体模拟为发现者-加入者模型,同时叠加侦察预警机制;发现者是容易找到食物的个体,其他个体为加入者,同时种群选取一定比例个体作为预警者,若探测到威胁则放弃在该位置觅食;基于标准SSA算法,提出改进的混沌自适应麻雀搜索算法,在算法初始化阶段采用混沌初始化种群策略增强算法稳定性,在发现者位置更新阶段引入了自适应权重因子,对精英麻雀个体采用柯西-高斯混合变异策略;
进一步地,所述标准SSA算法,其过程如下:
初始化:设置麻雀种群规模N,发现者个数Pd,预警者个数Sd,目标函数的维度D,初始值的上界ub,下界lb,最大迭代次数Tmax;
在D维搜索空间中,存在N只麻雀,则第i只麻雀在D维空间的位置为:
Xi=[x1,…xD],
其中i=1,2,...,N,
第i只麻雀的适应度值可以表示为:
其中f表示适应度值;
具有适应度好的发现者在觅食过程中会优先获取食物,并且为所有加入者提供觅食方向,因此发现者比加入者拥有更大的搜索范围;其位置更新公式如下:
式中,Xi表示第i只麻雀在D维空间的位置,t代表算法当前迭代的次数;α∈(0,1]代表一个随机数;R2∈(0,1]代表预警值,ST∈[0.5,1)代表安全值;Q代表服从正态分布的随机数;L代表所有元素都为1的1行d列的矩阵;当R2<ST时,意味觅食环境安全,发现者可以执行大范围的搜索操作,若R2≥ST,则表示部分麻雀发现了捕食者并发出警报信号,此时所有麻雀都需要转移至安全位置;
种群其余麻雀皆为加入者,加入者位置更新公式为:
其中,Xi表示第i只麻雀在D维空间的位置,t代表算法当前迭代的次数,Tmax代表算法最大迭代次数,Xbest是当前种群已占的最佳位置,即精英麻雀个***置,Xworst是最差位置,Q代表服从正态分布的随机数,A为一d×d矩阵,该矩阵每个元素被随机赋值1或-1;L代表所有元素都为1的1行D列的矩阵,当则在最佳位置附近觅食,/>时则第i加入者没有获取食物,需要飞往它处觅食;/>
种群所有麻雀都具有侦察预警机制,一般探测到危险的麻雀占种群的10%~20%,侦察预警麻雀的位置更新公式为:
式中,Xi表示第i只麻雀在D维空间的位置,t代表算法当前迭代的次数,Xbest是当前种群已占的最佳位置,β是步长参数,为均值=0,方差=1且服从正态分布的随机数;K∈[-1,1]为随机数,fi为第i个麻雀的适应度,fb是当前最佳适应度值,fw是当前最差适应度值;ε为为一常数,取值在(1.00E-10,1.00E-9)之间,可让式子避免出现分母等于零的情况。
进一步地,所述改进的混沌自适应麻雀搜索算法,具体如下:
采用立方映射混沌算子对麻雀算法进行种群初始化:
pi+1=4pi 3-3pi
-1<pi<1,pi≠0,i=0,1,...,N
Xi=Xlb+(Xlb-Xub)×(pi+1)×0.5
其中,p表示立方映射混沌算子,是一个每维取值均为(-1,1)的D维向量,式中Xi为麻雀个体变量值,Xlb,Xub为麻雀个体在各个维度的上下界;
采用立方映射混沌算子初始化种群大小为N的麻雀种群:先随机产生一个d维向量,每维取值均为(-1,1),作为第一个算子,然后对第一个算子的每一维做迭代运算,得到剩余的(N-1)个算子,最后采用(10)式将立方映射产生的算子值映射到麻雀个体上,式中Xi为麻雀个体变量值,Xlb,Xub为麻雀个体在各个维度的上下界;
引入自适应权重因子策略:
在发现者的位置更新公式中引入非线性时变的自适应权重因子w,通过引入自适应权重因子,在迭代初期,小的适应度值使算法有大的寻优范围,而在后期,大的适应度值则利于提高算法的收敛精度;标准SSA发现者位置更新公式中的α是一个随机数,引入了一个动态非线性变化的权重因子来控制α的取值范围,其定义如下:
式中fi是当前麻雀个体的适应度值,fg是全局最佳适应度值,t是当前迭代次数;
引入柯西-高斯混合变异策略:
在麻雀算法的求解后期,麻雀种群逐渐都会向最优个体聚集,这会导致种群多样性不足,令算法陷入停滞。为了解决这一问题,本发明再引入柯西-高斯混合变异策略,对当前适应度最好的个体进行变异操作,并从变异前后的个体选择最优个体进入下一次迭代,其定义如下:
式中,是当前种群适应度最好的精英麻雀个***置,/>是精英个体变异后的位置,cauchy(0,1)是满足柯西分布的随机变量,Gauss(0,1)是满足高斯分布的随机变量,/>为随迭代次数自适应调整的动态参数。
本发明与现有技术相比,具有如下优点和有益效果:
本发明所提出改进麻雀算法在求解质量上有着明显的优势,其混沌策略和自适应策略让算法拥有快速的收敛速度,和优秀的收敛精度,变异策略则让算法拥有很强的跳出局部最优的能力,进而可快速得到一条优秀的无人机飞行路线。
附图说明
图1是本发明所述基于混沌自适应麻雀搜索算法的无人机三维航迹规划方法流程图;
图2是本发明所述实施例中三维地图模型示意图;
图3是本发明所述实施例中威胁分布模型示意图;
图4是本发明所述实施例中混沌自适应麻雀搜索算法流程图;
图5是本发明所述实施例中基于CASSA算法的无人机三维航迹规划仿真图;
图6是本发明所述实施例中二维航迹规划仿真图;
图7是本发明所述实施例中CASSA算法收敛曲线图。
具体实施方式
下面结合实施例及附图对本发明作进一步详细的描述,但本发明的实施方式不限于此。
实施例:
基于混沌自适应麻雀搜索算法的无人机三维航迹规划方法,如图1所示,包括以下步骤:
根据飞行环境建立飞行环境模型;
建立三维航迹规划模型,即建立无人机飞行代价函数,对无人机航迹性能进行评价;
提出改进的混沌自适应麻雀搜索算法,具体为:在算法初始化阶段采用混沌初始化种群策略增强算法稳定性,在发现者位置更新阶段引入了自适应权重因子,对精英麻雀个体采用柯西-高斯混合变异策略;
通过混沌自适应麻雀搜索算法应用于航迹规划模型,得到无人机航迹规划最优解,即规划结果。
具体如下:
第一步:建立飞行环境模型:首先采集目标飞行区域的数字高程模型(DEM),通过MATLAB处理,得到三维地图模型,如图2所示。另外,无人机在执行飞行任务时经常会遇到具有威胁飞行安全的区域,称之为威胁区域,这些区域可以是敌人的雷达或者导弹***等等,无人机一旦进入威胁区域很有可能坠毁。为了简化模型,本发明采用半径为r的圆柱形区域表示威胁区域。设定无人机飞行起始点(10,90)与目标点坐标(140,10),以及威胁分布坐标(10,60),(30,40),(60,20),(70,60),(120,70),(110,30),威胁分布模型如图3所示。
第二步:建立三维航迹规划模型,建立合理的无人机飞行代价函数,对无人机航迹性能进行评价。无人机三维航迹规划的性能指标主要包括航迹长度,飞行高度,最大转角,具体如下:。
航迹长度对于规划任务来说是非常重要的,较短的航迹可以节省更多的燃料和更多的时间,并且碰上其它未知威胁的概率也会更低。航迹的定义为无人机从起点到终点所经历路程的值。
一条完整的航迹有n个节点,其中第i个航迹点和第i+1个航迹点之间的距离为li,两个航迹点的坐标分别表示为g(i)=(xi,yi,zi),g(i+1)=(xi+1,yi+1,zi+1),航迹需要满足以下条件:
其中Lpath为航迹长度代价函数,若无人机陷入危险区域或碰撞上障碍物面临坠毁危险,记为Lpath=∞,在仿真中采用惩罚的方式处理。
飞行高度需要足够的稳定性,太高的高度不利于避免未知威胁,太低的高度会加大无人机同山体障碍碰撞几率。本发明使用以下飞行高度模型:
转角代价函数影响无人机的可操作性和飞行稳定性,航迹规划过程中,航迹的转角应不大于预设的最大转角。本发明使用以下转角代价模型:
其中,Jturn是转角代价函数,φ为最大转角,θ是当前转角,ai是第i段航迹向量,|ai|表示矢量a的长度。
通过对航迹长度,飞行高度,最大转角建立代价函数,得到的航迹代价函数如下:
Jcost=w1Lpath+w2Hheight+w3Jturn (4)
其中,Jturn是总的航迹代价函数,参数wi,i=1,2,3表示每个代价函数的权重,且满足以下条件:
通过对代价函数的求解处理,可以得到由线段组成的航迹,再使用插值的方式调用B样条曲线函数,对已获取的航迹再进行平滑处理,最终得到无人机实际可飞的航迹。
第三步:基于混沌自适应麻雀搜索算法的航迹规划:提出混沌自适应麻雀搜索算法,在算法初始化阶段采用混沌初始化种群策略以增强算法稳定性,在发现者位置更新阶段引入了自适应权重因子,再对精英麻雀个体采用柯西-高斯混合变异策略提升算法跳出局部最优能力。
标准SSA算法是通过模拟麻雀的觅食过程来寻找优化问题的解。其算法原理为:将麻雀群体模拟为发现者-加入者模型,同时叠加侦察预警机制。发现者是容易找到食物的个体,其他个体为加入者,同时种群选取一定比例个体作为预警者,若探测到威胁则放弃在该位置觅食。
step1:初始化:设置麻雀种群规模N,发现者个数Pd,预警者个数Sd,目标函数的维度D,初始值的上界ub,下界lb,最大迭代次数Tmax。
step3:具有较好适应度的发现者在觅食过程中会优先获取食物,并且为所有加入者提供觅食方向,因此发现者比加入者拥有更大的搜索范围。其位置更新公式如下:
式中,t代表算法当前迭代的次数。α∈(0,1]代表一个随机数。R2∈(0,1]代表预警值,ST∈[0.5,1)代表安全值。Q代表服从正态分布的随机数。L代表所有元素都为1的1行d列的矩阵。当R2<ST时,意味觅食环境安全,发现者可以执行大范围的搜索操作,若R2≥ST,则表示部分麻雀发现了捕食者并发出警报信号,此时所有麻雀都需要转移至安全位置。
step4:种群其余麻雀皆为加入者,加入者位置更新公式为:
其中,Xbest是当前种群已占的最佳位置,即精英麻雀个***置,Xworst是最差位置,A为一d×d矩阵,该矩阵每个元素被随机赋值1或-1。当则在最佳位置附近觅食,时则第i加入者没有获取食物,需要飞往它处觅食。
Step5:种群所有麻雀都具有侦察预警机制,一般探测到危险的麻雀占种群的10%~20%,侦察预警麻雀的位置更新公式为:
式中,β是步长参数,为均值=0,方差=1且服从正态分布的随机数。K∈[-1,1]为随机数,fi为第i个麻雀的适应度,fb是当前最佳适应度值,fw是当前最差适应度值。ε为一常数,取值在(1.00E-10,1.00E-9)之间,可让式子避免出现分母等于零的情况。
第四步:混沌自适应麻雀搜索算法的具体改进及应用如下:
混沌自适应麻雀搜索算法流程图如图4所示:
引入混沌初始化种群策略:
麻雀算法在求解复杂优化问题时,在迭代后期种群存在着种群多样性降低的缺点,易于早熟收敛,陷入局部最优。因此本发明采用立方映射混沌算子对麻雀算法进行种群初始化,根据混沌算子具有的随机性和规律性优点,可以提高算法的种群多样性。
Xi=Xlb+(Xlb-Xub)×(pi+1)×0.5 (10)
采用立方映射混沌算子初始化种群大小为N的麻雀种群:先随机产生一个d维向量,每维取值均为(-1,1),作为第一个算子,然后采用(9)式对第一个算子的每一维做迭代运算,得到剩余的(N-1)个算子,最后采用(10)式将立方映射产生的算子值映射到麻雀个体上,式中Xi为麻雀个体变量值,Xlb,Xub为麻雀个体在各个维度的上下界。
引入自适应权重因子策略:
在step 3的发现者的位置更新公式中引入非线性时变的自适应权重因子w,通过引入自适应权重因子,在迭代初期,较小的适应度值使算法有较大的寻优范围,而在后期,较大的适应度值则利于提高算法的收敛精度。标准SSA发现者位置更新公式中的α是一个随机数,本发明引入了一个动态非线性变化的权重因子来控制α的取值范围,其定义如下:
式中fi是当前麻雀个体的适应度值,fg是全局最佳适应度值,t是当前迭代次数。
引入柯西-高斯混合变异策略:
在麻雀算法的求解后期,麻雀种群逐渐都会向最优个体聚集,这会导致种群多样性不足,令算法陷入停滞。为了解决这一问题,本发明再引入柯西-高斯混合变异策略,对当前适应度最好的个体进行变异操作,并从变异前后的个体选择最优个体进入下一次迭代,其定义如下:
式中,是当前种群适应度最好的精英麻雀个***置,/>是精英个体变异后的位置,cauchy(0,1)是满足柯西分布的随机变量,Gauss(0,1)是满足高斯分布的随机变量,/> 为随迭代次数自适应调整的动态参数。
仿真验证:
运用CASSA算法求解无人机三维航迹规划问题,实际上就是用CASSA对目标函数进行寻优求解。在算法中,每一只麻雀代表一条航迹,每只麻雀的维数代表航迹点的个数,麻雀个体在各维度的取值代表航迹点的坐标值,麻雀个体适应度即为无人机航迹的目标函数值,是越低越好的。麻雀种群发现者是当前目标函数值较优个体,跟随者是较差个体,通过对发现者和跟随者的位置更新来得到麻雀种群的最优位置,算法达到最大迭代次数后返回的最优解即是所寻找到的最优航迹。基于CASSA算法的无人机三维航迹规划仿真图如图5所示,图6为二维航迹规划仿真图;图7为CASSA算法收敛曲线图。
具体实现流程如下:
第一步:参数初始化,对CASSA算法的各参数初始化,如最大迭代次数,种群规模,发现者和预警者数量等,对无人机航迹规划的三种代价函数的权重进行设定。
第二步:设置无人机的飞行起点坐标和终点坐标,设定威胁区域的坐标和威胁范围大小。
第三步:运用CASSA算法对无人机航迹评价函数进行优化求解,不断更新并存储已获取的最优航迹。
第四步:判断算法是否达到最大迭代次数,若无则继续执行循环,若达到则算法结束,输出已获取的最优航迹。
上述实施例为本发明较佳的实施方式,但本发明的实施方式并不受上述实施例的限制,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。
Claims (5)
1.基于混沌自适应麻雀搜索算法的无人机三维航迹规划方法,其特征在于,包括以下步骤:
根据飞行环境建立飞行环境模型;
建立无人机飞行代价函数,对无人机航迹性能进行评价;
采用混沌初始化种群策略,自适应权重策略和柯西-高斯混合变异策略对麻雀搜索算法进行改进,提出混沌自适应麻雀搜索算法;
采用混沌自适应麻雀搜索算法对三维环境中的无人机航迹进行规划,得到无人机航迹规划最优解,得到规划结果;
所述采用混沌初始化种群策略,自适应权重策略和柯西-高斯混合变异策略对麻雀搜索算法进行改进,提出混沌自适应麻雀搜索算法,如下所示:
采用标准SSA算法,将麻雀群体模拟为发现者-加入者模型,同时叠加侦察预警机制;发现者是容易找到食物的个体,其他个体为加入者,同时种群选取一定比例个体作为预警者,若探测到威胁则放弃在该位置觅食;基于标准SSA算法,提出改进的混沌自适应麻雀搜索算法,在算法初始化阶段采用混沌初始化种群策略增强算法稳定性,在发现者位置更新阶段引入了自适应权重因子,对精英麻雀个体采用柯西-高斯混合变异策略;
所述标准SSA算法,其过程如下:
初始化:设置麻雀种群规模N,发现者个数Pd,预警者个数Sd,目标函数的维度D,初始值的上界ub,下界lb,最大迭代次数ss;
在D维搜索空间中,存在N只麻雀,则第i只麻雀在D维空间的位置为:
Xi=[x1,...xD],
其中i=1,2,...,N,
第i只麻雀的适应度值可以表示为:
其中f表示适应度值;
具有适应度好的发现者在觅食过程中会优先获取食物,并且为所有加入者提供觅食方向,因此发现者比加入者拥有更大的搜索范围;其位置更新公式如下:
式中,Xi表示第i只麻雀在D维空间的位置,t代表算法当前迭代的次数;α∈(0,1]代表一个随机数;R2∈(0,1]代表预警值,ST∈[0.5,1)代表安全值;Q代表服从正态分布的随机数;L代表所有元素都为1的1行d列的矩阵;当R2<ST时,意味觅食环境安全,发现者可以执行大范围的搜索操作,若R2≥ST,则表示部分麻雀发现了捕食者并发出警报信号,此时所有麻雀都需要转移至安全位置;
种群其余麻雀皆为加入者,加入者位置更新公式为:
其中,Xi表示第i只麻雀在D维空间的位置,t代表算法当前迭代的次数,Tmax代表算法最大迭代次数,Xbest是当前种群已占的最佳位置,即精英麻雀个***置,Xworst是最差位置,Q代表服从正态分布的随机数,A为一d×d矩阵,该矩阵每个元素被随机赋值1或-1;L代表所有元素都为1的1行D列的矩阵,当则在最佳位置附近觅食,/>时则第i加入者没有获取食物,需要飞往它处觅食;
种群所有麻雀都具有侦察预警机制,一般探测到危险的麻雀占种群的10%~20%,侦察预警麻雀的位置更新公式为:
式中,Xi表示第i只麻雀在D维空间的位置,t代表算法当前迭代的次数,Xbest是当前种群已占的最佳位置,β是步长参数,为均值=0,方差=1且服从正态分布的随机数;K∈[-1,1]为随机数,fi为第i个麻雀的适应度,fb是当前最佳适应度值,fw是当前最差适应度值;ε为常数,取值在1.00E-10,1.00E-9之间,使得公式避免出现分母等于零的情况;
所述改进的混沌自适应麻雀搜索算法,具体如下:
采用立方映射混沌算子对麻雀算法进行种群初始化:
pi+1=4pi 3-3pi
-1<pi<1,pi≠0,i=0,1,K,N
Xi=Xlb+(Xlb-Xub)×(pi+1)×0.5
其中,p表示立方映射混沌算子,是一个每维取值均为(-1,1)的D维向量,,式中Xi为麻雀个体变量值,Xlb,Xub为麻雀个体在各个维度的上下界;
采用立方映射混沌算子初始化种群大小为N的麻雀种群:先随机产生一个d维向量,每维取值均为(-1,1),作为第一个算子,然后对第一个算子的每一维做迭代运算,得到剩余的(N-1)个算子,最后将立方映射产生的算子值映射到麻雀个体上,式中Xi为麻雀个体变量值,Xlb,Xub为麻雀个体在各个维度的上下界;
引入自适应权重因子策略:
在发现者的位置更新公式中引入非线性时变的自适应权重因子w,通过引入自适应权重因子,在迭代初期,小的适应度值使算法有大的寻优范围,而在后期,大的适应度值则利于提高算法的收敛精度;标准SSA发现者位置更新公式中的α是一个随机数,引入了一个动态非线性变化的权重因子来控制α的取值范围,其定义如下:
式中fi是当前麻雀个体的适应度值,fg是全局最佳适应度值,t是当前迭代次数;
引入柯西-高斯混合变异策略:
在麻雀算法的求解后期,麻雀种群逐渐都会向最优个体聚集,这会导致种群多样性不足,令算法陷入停滞;为了解决这一问题,引入柯西-高斯混合变异策略,对当前适应度最好的个体进行变异操作,并从变异前后的个体选择最优个体进入下一次迭代,其定义如下:
2.根据权利要求1所述的基于混沌自适应麻雀搜索算法的无人机三维航迹规划方法,其特征在于,所述根据飞行环境建立飞行环境模型,具体为:采集目标飞行区域的数字高程模型,通过MATLAB处理,得到三维地图模型。
3.根据权利要求2所述的基于混沌自适应麻雀搜索算法的无人机三维航迹规划方法,其特征在于,所述三维地图模型包括威胁区域,在三维地图模型中设定飞行起始点和目标坐标点以及威胁分布坐标。
4.根据权利要求1所述的基于混沌自适应麻雀搜索算法的无人机三维航迹规划方法,其特征在于,所述建立无人机飞行代价函数,对无人机航迹性能进行评价,具体为:
设定无人机三维航迹规划的性能指标:航迹长度,飞行高度,最大转角;
通过航迹长度代价函数获取航迹长度;通过飞行高度模型获取飞行高度;通过转角代价模型获取最大转角;
通过航迹长度代价函数、飞行高度模型、转角代价模型,得到航迹代价函数:
Jcost=w1Lpath+w2Hheight+w3Jturn,
其中,Jcost是总的航迹代价函数,Lpath为航迹长度代价函数,Hheight为高度的标准差代价函数,Jturn是转角代价函数,参数wi,i=1,2,3表示每个代价函数的权重,且满足以下条件:
通过对代价函数的求解处理,得到由线段组成的航迹,再使用插值的方式调用B样条曲线函数,对已获取的航迹再进行平滑处理,最终得到无人机实际可飞的航迹。
5.根据权利要求4所述的基于混沌自适应麻雀搜索算法的无人机三维航迹规划方法,其特征在于,所述通过航迹长度代价函数获取航迹长度,如下所示:设定完整的航迹有n个节点,第i个航迹点和第i+1个航迹点之间的距离为li,两个航迹点的坐标分别表示为g(i)=(xi,yi,zi),g(i+1)=(xi+1,yi+1,zi+1),航迹需要满足以下条件:
其中,Lpath为航迹长度代价函数,
所述通过飞行高度模型获取飞行高度,如下所示:飞行高度模型为:
所述通过转角代价模型获取最大转角,如下所示:转角代价模型如下:
其中,Jturn是转角代价函数,φ为最大转角,θ是当前转角,ai是第i段航迹向量,|ai|表示矢量ai的长度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110110630.4A CN112880688B (zh) | 2021-01-27 | 2021-01-27 | 基于混沌自适应麻雀搜索算法的无人机三维航迹规划方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110110630.4A CN112880688B (zh) | 2021-01-27 | 2021-01-27 | 基于混沌自适应麻雀搜索算法的无人机三维航迹规划方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112880688A CN112880688A (zh) | 2021-06-01 |
CN112880688B true CN112880688B (zh) | 2023-05-23 |
Family
ID=76053375
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110110630.4A Active CN112880688B (zh) | 2021-01-27 | 2021-01-27 | 基于混沌自适应麻雀搜索算法的无人机三维航迹规划方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112880688B (zh) |
Families Citing this family (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113344828B (zh) * | 2021-06-30 | 2023-05-30 | 太原理工大学 | 一种基于混沌麻雀搜索算法的煤岩分界图像增强方法 |
CN113759313B (zh) * | 2021-07-23 | 2023-09-29 | 哈尔滨工程大学 | 一种基于混沌麻雀算法的时差/频差定位方法 |
CN113762603B (zh) * | 2021-08-13 | 2023-08-22 | 广西大学 | 基于改进麻雀算法优化的光伏基站短期光伏功率预测方法 |
CN113805568B (zh) * | 2021-08-17 | 2024-04-09 | 北京理工大学 | 基于多智能体时空建模与决策的人机协同感知方法 |
CN113947009B (zh) * | 2021-08-31 | 2024-03-08 | 西北大学 | 基于柯西扰动麻雀优化的稀疏lstm滑坡动态预测方法 |
CN114115329B (zh) * | 2021-10-28 | 2023-06-16 | 中国人民解放军陆军工程大学 | 中继协同无人机任务规划方法及装置 |
CN114020031B (zh) * | 2021-11-16 | 2023-09-22 | 南京航空航天大学 | 基于改进鸽群优化的无人机集群协同动态目标搜索方法 |
CN114237282B (zh) * | 2021-11-29 | 2023-10-27 | 南京邮电大学 | 面向智慧化工业园区监测的无人机飞行路径智能规划方法 |
CN113867369B (zh) * | 2021-12-03 | 2022-03-22 | 中国人民解放军陆军装甲兵学院 | 一种基于交流学习的海鸥算法的机器人路径规划方法 |
CN114485665A (zh) * | 2021-12-31 | 2022-05-13 | 东华大学 | 一种基于麻雀搜索算法的无人机航迹规划方法 |
CN114417712B (zh) * | 2022-01-01 | 2023-07-04 | 西北工业大学 | 基于混沌初始化ssa-bp神经网络的飞艇螺旋桨可靠性估计方法 |
CN114460941B (zh) * | 2022-02-09 | 2023-10-24 | 云南大学 | 一种基于改进麻雀搜索算法的机器人路径规划方法及*** |
CN114440893B (zh) * | 2022-02-16 | 2024-07-09 | 北京邮电大学 | 用于tdoa信号解算的协同定位方法、***及存储介质 |
CN115018149B (zh) * | 2022-06-01 | 2023-04-28 | 南京林业大学 | 一种直升机航空施药飞机起降点位置规划方法 |
CN114995149B (zh) * | 2022-06-13 | 2023-06-13 | 南通大学 | 液压位置伺服***改进混沌变权麻雀搜索参数辨识方法 |
CN115277770B (zh) * | 2022-07-20 | 2023-04-25 | 华北电力大学(保定) | 一种节点接入和飞行策略联合优化的无人机信息收集方法 |
CN115278774A (zh) * | 2022-07-20 | 2022-11-01 | 云南电网有限责任公司电力科学研究院 | 一种北斗短报文缺失数据补录方法及*** |
CN115542944B (zh) * | 2022-10-24 | 2024-06-28 | 广东电网有限责任公司云浮供电局 | 一种基于配电网环境的多无人机路径规划方法及相关装置 |
CN115542936B (zh) * | 2022-10-28 | 2024-04-16 | 中国地质大学(武汉) | 一种航磁无人机舵机布局智能优化方法、装置及电子设备 |
CN115657693B (zh) | 2022-12-28 | 2023-03-17 | 安徽省交通航务工程有限公司 | 一种船舶路径优化方法、电子设备及存储介质 |
CN116562331B (zh) * | 2023-05-19 | 2023-11-21 | 石家庄铁道大学 | 一种改进爬行动物搜索算法优化svm的方法及其应用 |
CN117034459B (zh) * | 2023-08-18 | 2024-05-31 | 华东交通大学 | 基于改进蜣螂优化算法的磁浮列车运行优化方法及*** |
CN117056847B (zh) * | 2023-10-10 | 2024-01-30 | 中南大学 | 一种流式数据的异常检测方法、***、设备及存储介质 |
CN117170413B (zh) * | 2023-11-03 | 2024-01-30 | 北京卓翼智能科技有限公司 | 基于改进正余弦算法的无人机路径规划方法和装置 |
CN117807356B (zh) * | 2024-02-29 | 2024-05-10 | 齐鲁工业大学(山东省科学院) | 基于改进麻雀算法优化粒子滤波的双矢量水听器定位方法 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1835370A2 (fr) * | 2006-03-17 | 2007-09-19 | Thales | Planificateur automatique de trajectoire |
CN108413959A (zh) * | 2017-12-13 | 2018-08-17 | 南京航空航天大学 | 基于改进混沌蚁群算法的无人机航路规划方法 |
US10140875B1 (en) * | 2017-05-27 | 2018-11-27 | Hefei University Of Technology | Method and apparatus for joint optimization of multi-UAV task assignment and path planning |
CN109144102A (zh) * | 2018-09-19 | 2019-01-04 | 沈阳航空航天大学 | 一种基于改进蝙蝠算法的无人机航路规划方法 |
CN109631900A (zh) * | 2018-12-24 | 2019-04-16 | 中国矿业大学 | 一种无人机三维航迹多目标粒子群全局规划方法 |
CN110608743A (zh) * | 2019-10-18 | 2019-12-24 | 南京航空航天大学 | 基于多种群混沌灰狼算法的多无人机协同航路规划方法 |
JP2020019336A (ja) * | 2018-07-31 | 2020-02-06 | 旻新科技股▲分▼有限公司 | 鳥の飛行ルートを使用する、無人航空機の飛行経路の計画方法 |
CN111815055A (zh) * | 2020-07-11 | 2020-10-23 | 哈尔滨理工大学 | 一种基于改进莱维飞行蚁狮优化算法的无人机航迹规划方法 |
CN112033416A (zh) * | 2020-09-18 | 2020-12-04 | 中原工学院 | 一种无人机静态航迹规划方法和装置 |
CN112130581A (zh) * | 2020-08-19 | 2020-12-25 | 昆明理工大学 | 一种面向空中机动作战的无人机集群协同任务规划方法 |
CN112230678A (zh) * | 2020-10-29 | 2021-01-15 | 皖江工学院 | 基于粒子群算法的三维无人机路径规划方法及规划*** |
-
2021
- 2021-01-27 CN CN202110110630.4A patent/CN112880688B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1835370A2 (fr) * | 2006-03-17 | 2007-09-19 | Thales | Planificateur automatique de trajectoire |
US10140875B1 (en) * | 2017-05-27 | 2018-11-27 | Hefei University Of Technology | Method and apparatus for joint optimization of multi-UAV task assignment and path planning |
CN108413959A (zh) * | 2017-12-13 | 2018-08-17 | 南京航空航天大学 | 基于改进混沌蚁群算法的无人机航路规划方法 |
JP2020019336A (ja) * | 2018-07-31 | 2020-02-06 | 旻新科技股▲分▼有限公司 | 鳥の飛行ルートを使用する、無人航空機の飛行経路の計画方法 |
CN109144102A (zh) * | 2018-09-19 | 2019-01-04 | 沈阳航空航天大学 | 一种基于改进蝙蝠算法的无人机航路规划方法 |
CN109631900A (zh) * | 2018-12-24 | 2019-04-16 | 中国矿业大学 | 一种无人机三维航迹多目标粒子群全局规划方法 |
CN110608743A (zh) * | 2019-10-18 | 2019-12-24 | 南京航空航天大学 | 基于多种群混沌灰狼算法的多无人机协同航路规划方法 |
CN111815055A (zh) * | 2020-07-11 | 2020-10-23 | 哈尔滨理工大学 | 一种基于改进莱维飞行蚁狮优化算法的无人机航迹规划方法 |
CN112130581A (zh) * | 2020-08-19 | 2020-12-25 | 昆明理工大学 | 一种面向空中机动作战的无人机集群协同任务规划方法 |
CN112033416A (zh) * | 2020-09-18 | 2020-12-04 | 中原工学院 | 一种无人机静态航迹规划方法和装置 |
CN112230678A (zh) * | 2020-10-29 | 2021-01-15 | 皖江工学院 | 基于粒子群算法的三维无人机路径规划方法及规划*** |
Non-Patent Citations (5)
Title |
---|
A novel swarm intelligence optimization approach:sparrow search algorithm;XUE J等;《Systems Science and Control Engineering》;第8卷(第1期);第22-34页 * |
Obstacle Avoidance of UAV Based on Neural Networks and Interfered Fluid Dynamical System;Wang Y等;《2020 3rd International Conference on Unmanned Systems》;第1066-1071页 * |
一种新型的群智能优化技术的研究与应用 ——麻雀搜索算法;薛建凯;《中国优秀硕士学位论文全文数据库信息科技辑(月刊)》(第1期);第I140-134页 * |
基于动态路径规划的车载自组织网络贪婪路由协议;赵普凡等;《国外电子测量技术》;第39卷(第11期);第46-51页 * |
基于混沌麻雀搜索算法的无人机航迹规划方法;汤安迪等;《计算机应用》;第41卷(第7期);第1-11页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112880688A (zh) | 2021-06-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112880688B (zh) | 基于混沌自适应麻雀搜索算法的无人机三维航迹规划方法 | |
CN110608743B (zh) | 基于多种群混沌灰狼算法的多无人机协同航路规划方法 | |
CN113791634B (zh) | 一种基于多智能体强化学习的多机空战决策方法 | |
Yang et al. | UAV air combat autonomous maneuver decision based on DDPG algorithm | |
CN108153328B (zh) | 一种基于分段贝塞尔曲线的多导弹协同航迹规划方法 | |
CN111240353B (zh) | 基于遗传模糊树的无人机协同空战决策方法 | |
CN110456823B (zh) | 一种针对无人机计算与存储能力受限的双层路径规划方法 | |
CN108549233B (zh) | 一种带有直觉模糊信息的无人机空战机动博弈方法 | |
CN110928329B (zh) | 一种基于深度q学习算法的多飞行器航迹规划方法 | |
CN110412874B (zh) | 针对机动目标和时延通信的多导弹协同制导律设计方法 | |
CN111859541B (zh) | 一种基于迁移学习改进的pmaddpg多无人机任务决策方法 | |
CN114330115B (zh) | 一种基于粒子群搜索的神经网络空战机动决策方法 | |
CN110187640B (zh) | 针对机动目标和允许通信时滞的多导弹协同作战制导律设计方法 | |
CN112666981B (zh) | 基于原鸽群动态群组学习的无人机集群动态航路规划方法 | |
CN113962012B (zh) | 无人机对抗策略优化方法及装置 | |
CN108983800B (zh) | 一种基于深度学习的飞机姿态控制方法 | |
CN114063644B (zh) | 基于鸽群反向对抗学习的无人作战飞机空战自主决策方法 | |
CN110986958B (zh) | 基于多种群协作果蝇优化的多无人机协同路径规划方法 | |
CN113282061A (zh) | 一种基于课程学习的无人机空中博弈对抗的解决方法 | |
CN113159266B (zh) | 基于麻雀搜索神经网络的空战机动决策方法 | |
CN113741186B (zh) | 一种基于近端策略优化的双机空战决策方法 | |
Leuenberger et al. | Actor-critic reinforcement learning with neural networks in continuous games | |
CN117313561B (zh) | 无人机智能决策模型训练方法及无人机智能决策方法 | |
Kaneshige et al. | Artificial immune system approach for air combat maneuvering | |
CN116661496B (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 |