CN113552891A - 一种基于改进的蝴蝶优化算法的机器人多目标路径规划 - Google Patents
一种基于改进的蝴蝶优化算法的机器人多目标路径规划 Download PDFInfo
- Publication number
- CN113552891A CN113552891A CN202110996730.1A CN202110996730A CN113552891A CN 113552891 A CN113552891 A CN 113552891A CN 202110996730 A CN202110996730 A CN 202110996730A CN 113552891 A CN113552891 A CN 113552891A
- Authority
- CN
- China
- Prior art keywords
- butterfly
- formula
- algorithm
- function
- path
- 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
Links
- 238000005457 optimization Methods 0.000 title claims abstract description 26
- 230000006870 function Effects 0.000 claims description 55
- 239000003205 fragrance Substances 0.000 claims description 22
- 241000255777 Lepidoptera Species 0.000 claims description 12
- 230000035772 mutation Effects 0.000 claims description 12
- 238000005312 nonlinear dynamic Methods 0.000 claims description 12
- 239000006185 dispersion Substances 0.000 claims description 11
- 238000000034 method Methods 0.000 claims description 6
- 230000001131 transforming effect Effects 0.000 claims description 6
- 230000003044 adaptive effect Effects 0.000 claims description 5
- 230000004913 activation Effects 0.000 claims description 4
- 238000013528 artificial neural network Methods 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 4
- 210000002569 neuron Anatomy 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 4
- 238000012935 Averaging Methods 0.000 claims description 3
- 238000013016 damping Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 230000006872 improvement Effects 0.000 claims description 2
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 230000035943 smell Effects 0.000 description 5
- 230000009471 action Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000002349 favourable effect Effects 0.000 description 2
- 230000002028 premature Effects 0.000 description 2
- 108091026890 Coding region Proteins 0.000 description 1
- 241001414834 Ephemeroptera Species 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000019637 foraging behavior Effects 0.000 description 1
- 230000013011 mating Effects 0.000 description 1
- 230000007758 mating behavior Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/02—Control of position or course in two dimensions
- G05D1/021—Control of position or course in two dimensions specially adapted to land vehicles
- G05D1/0212—Control of position or course in two dimensions specially adapted to land vehicles with means for defining a desired trajectory
- G05D1/0214—Control of position or course in two dimensions specially adapted to land vehicles with means for defining a desired trajectory in accordance with safety or protection criteria, e.g. avoiding hazardous areas
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/02—Control of position or course in two dimensions
- G05D1/021—Control of position or course in two dimensions specially adapted to land vehicles
- G05D1/0212—Control of position or course in two dimensions specially adapted to land vehicles with means for defining a desired trajectory
- G05D1/0219—Control of position or course in two dimensions specially adapted to land vehicles with means for defining a desired trajectory ensuring the processing of the whole working surface
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/02—Control of position or course in two dimensions
- G05D1/021—Control of position or course in two dimensions specially adapted to land vehicles
- G05D1/0212—Control of position or course in two dimensions specially adapted to land vehicles with means for defining a desired trajectory
- G05D1/0221—Control of position or course in two dimensions specially adapted to land vehicles with means for defining a desired trajectory involving a learning process
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)
- Feedback Control In General (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
- Numerical Control (AREA)
Abstract
本发明涉及一种基于改进的蝴蝶优化算法的机器人多目标路径规划,可以最大程度的发挥干扰项的作用,使其具有良好的全局探索能力,有助于算法避免陷入局部最优;但在算法的运行后期,变异干扰项变为高斯分布变异,这将加快其收敛速度及提高搜索的精度,算法具有较强的局部搜索能力,在机器人多目标路径规划应用中,会提高机器人搜索的精度,使机器人具有较强得到局部搜索能力,使得机器人能够快速在起始点和终止点之间找到一条最短的、较平滑的、且避开所有障碍物的路径。
Description
技术领域
本发明具体涉及一种基于改进的蝴蝶优化算法的机器人多目标路径规划。
背景技术
近年来一种新型优化算法—蝴蝶优化算法得到了应用,它由Sankalap Arora和Satvir Singh等于2019年提出,该算法的灵感来源于蝴蝶的觅食和交配行为。蝴蝶寻找食物和伴侣主要依靠嗅觉、视觉、味觉、触觉和听觉,其中嗅觉是最重要的,它帮助蝴蝶寻找食物或者最佳的交配对象。
1算法简介
1.1基本原理
蝴蝶的嗅觉对其周围的气味非常敏感,科学研究表明,蝴蝶对这些气味来源的位置有着非常准确的判断。同时每只蝴蝶均能释放出一种具有特殊味道和浓度的香气,当蝴蝶感觉到另一只蝴蝶在这个区域散发出更多的香味时,就会去靠近,这个阶段相当于算法的全局搜索阶段。另外,当蝴蝶感知不到大于它自己的香味时,就会随机移动,相当于算法的局部搜索阶段。
将香味建模为刺激物的物理强度的函数,如下式(1);
f=cIa (1)
上式中,f是感知到的其它蝴蝶香味强度,c表示感知形态,取值范围是[0,1], I、a分别表示刺激强度和幂指数,它们取值范围都是[0,1]。参数a控制着算法的行为,当其取最大值1时,表示一只特定蝴蝶发出的香气量被其它蝴蝶以同样的能力感知到,意味着在这个区域的任何位置都可以感觉到一只散发香味的蝴蝶,算法达到了全局的最优值。若参数,取值为0,则表示任何一只蝴蝶散发的香味都不会被其它蝴蝶感觉到。
1.2算法流程
基本的蝴蝶优化算法流程如下步骤所示:
(1)初始化蝴蝶种群,明确所求问题目标函数、问题的解空间,随机初始化第一代蝴蝶;
(2)根据适应度函数计算出每只蝴蝶在其解空间中的适应度值,找出当前最优蝴蝶的位置;
(3)计算每只蝴蝶发出香味的浓度函数。利用随机函数生成数决策该蝴蝶执行全局搜索还是局部搜索,以尽可能减少外部环境的影响;
(4)全局搜索阶段,低浓度香味的蝴蝶个体受浓度最佳的蝴蝶吸引并向其移动,其全局搜索移动公式如下:
(5)局部搜索阶段,蝴蝶个体进行随机搜索,其局部搜索移动公式如下:
(6)判断是否满足迭代次数要求或精度要求,则停止计算并输出最优解及相应的参数,否则转向(2)继续循环。
蝴蝶的位置更新对种群的寻优取决定性作用,在蝴蝶算法初期,种群中的个体模式集中在适应度值较低的个体上,需要种群尽快能够在整个空间进行“搜索”,但是到了算法后期,个体模式朝高适应度的个体集中,会造成种群的“开发”能力有限,进而出现早熟收敛,无法跳出局部极值。
目前蝴蝶算法被广泛应用在机器人领域,现有技术在采用这种蝴蝶算法对机器人进行多目标路径规划,在实际应用中机器人在运动过程中也会存在上述说明的技术问题。
发明内容
为解决现有技术中存在的问题,本发明提供了一种基于改进的蝴蝶优化算法的机器人多目标路径规划,具体的技术方案如下所述:
一种基于改进的蝴蝶优化算法,其特征在于:包括有如下步骤:
S1、参数初始化设置:
初始化蝴蝶优化算法的各项基本参数;
S2、初始化种群,并计算各个解的是适应度值:
初始化种群,根据适应度函数计算出每只蝴蝶在其解空间中的适应度值;
S3、计算每只蝴蝶发出香味的浓度函数:
利用随机函数生成数决策该蝴蝶执行全局搜索还是局部搜索,以尽可能减少外部环境的影响;
S4、判断P>rand:
若P>rand,则利用式(2)进行全局搜索:
低浓度香味的蝴蝶个体受浓度最佳的蝴蝶吸引并向其移动;
若P<rand,则按照式(7)进行局部搜索:
引入动态自适应t分布变异策略:
给出进化离散度的定义,以便描述种群整体适应度的变化,设o(d)为第d代种群与第 d-1代种群的适应度值标准差的比值:
在神经网络理论中,神经元的激活函数采用Sigmoid函数来构造,其定义:
点(0,0.5)为其中心对称点,当自变量x在[-10,10]之间时函数值域变化,且变化范围为(0,1),函数在线性与非线性之间具有良好的平衡性,作为阈值函数使用,结合进化离散度o(d)和Sigmoid函数,提出一种非线性动态惯性权重因子:
依据公式(6)生成的非线性动态惯性权重因子,在常规蝴蝶算法中进行t分布变异操作,其变异程度由非线性动态惯性权重因子ω进行控制,最终通过引入这种动态自适应t分布变异策略,对算法中蝴蝶局部搜索阶段的移动公式改进如下:
其中,tdis(d)表示当迭代次数为d时,以其作为参数自由度的t分布,公式(7)在蝴蝶当前位置基础上增加了自适应t分布的变异干扰项ω·tdis(d);
S5、更新蝴蝶个体和全局最优解;
S6、越界处理;
S7、判断是否满足迭代次数要求或精度要求,则停止计算并输出最优解及适应度值,否则转向S3继续循环。
进一步地,在公式(6)中,其中ωmax为最大惯性权重,取值为0.9,ωmin为最小惯性权重,取值为0.4。
进一步地,在公式(6)中,d、Dmax分别表示当前迭代次数与最大进化代数,o(d)为种群进化离散度参数,e是阻尼因子,取值为[0,1]。
进一步地,在步骤S4中,在算法的运行初期,由于d值较小,变异干扰项ω·tdis(d)视为柯西变异,在算法的运行后期,变异干扰项ω·tdis(d)视为高斯分布变异。
一种基于改进的蝴蝶优化算法的机器人多目标路径规划,其特征在于:包括有如下步骤:
1)进行坐标系变换,且利用式(8)对起点、终点、障碍物位置的坐标进行变换:
起点B与终点E的连线为X'轴构建坐标系BX'Y',然后将坐标系OXY中的点变换到BX'Y'中,变换公式如下:
其中,(xb,yb)是坐标系OXY中起点B的坐标,(x',y')是点(x,y)在坐标系BX'Y'中对应的点,θ是X轴与直线BE的夹角;
2)建立m条平行簇l1,l2,…,lm将BE平均为m+1段,每每相邻的两条平时直线的距离设为Δl=||BE||/(m+1);
3)初始化参数;
5)利用适应度函数式(12),对蝴蝶进行评估,获得蝴蝶的适应度S;
长度指标
设起点B和终点E的坐标为B(xb,yb)、E(xe,ye),任意路径节点坐标设为Pi(xi,yi),i的取值为[1,m],设路径长度为fdistance,将其标准化表示如下式:
安全度指标
设机器人在三次样条插值节点处与障碍物无碰撞,假设已知m个路径结点的坐标(x1,y1),(x2,y2),…,(xm,ym)以及起点坐标(xb,yb)和终点坐标(xe,ye),通过三次样条差值分别得到d个插值点,其横坐标为(x1,x2,…,xd),其纵坐标为(y1,y2,…,yd),需要所生成的d个插值点是否有碰撞情况发生,将d设置为100,设路径安全度指标为fsecure,将其如下式表示:
上式中η表示安全因子,取值为100,Dj,k表示第j个插值点到第k个障碍物中心距离, R(k)表示第k个障碍物的半径,H表示路径中障碍物的数量;
光滑度指标
设路径节点为(B,P1,P2,…,Pm,E),相邻路段之间的夹角表示为ψi,建立的路径平滑度指标为:
ψi越小表示路径光滑程度越优;
将以上指标进行融合后得到机器人路径规划的综合指标S:
S=fdistance+fsecure+fflatness (12);
6)利用公式(1)计算每只蝴蝶发出香味的浓度函数,利用随机函数生成数决策该蝴蝶执行全局搜索还是局部搜索,减少外部环境的影响:
全局搜索:
利用公式(2)进行全局搜索,低浓度香味的蝴蝶个体受浓度最佳的蝴蝶吸引并向其移动;
局部搜索:
引入动态自适应t分布变异策略利用公式(7)进行局部搜索,蝴蝶个体进行随机搜索;
7)更新蝴蝶个体和全局最优解;
8)越接处理;
9)判断是否满足迭代次数要求或精度要求,若是满足进入11),否则返回5)继续循环;
10)输出最优个体适应度值。
本发明的有益效果为:
本发明在算法的运行初期,由于当前迭代次数d值较小,此时该干扰项可以视为柯西变异,使得算法保持了蝴蝶种群的多样性,因而可以最大程度的发挥干扰项的作用,使其具有良好的全局探索能力,有助于算法避免陷入局部最优;但在算法的运行后期,变异干扰项变为高斯分布变异,这将加快其收敛速度及提高搜索的精度,算法具有较强的局部搜索能力,同时变异干扰项的作用随着迭代次数的增加不断减小;
基于改进的蝴蝶优化算法应用在机器人路径规划中,会提高机器人搜索的精度,使机器人具有较强得到局部搜索能力,使得机器人能够快速在起始点和终止点之间找到一条最短的、较平滑的、且避开所有障碍物的路径。
附图说明
图1为本发明实施例中Sigmoid函数图;
图2为本发明实施例中基于改进的蝴蝶优化算法流程图;
图3为本发明实施例中的机器人在环境中的路径表示函数图。
具体实施方式:
为了加深对本发明的理解,下面结合附图对本发明的实施例做详细的说明。
须知,本说明书所附图式所绘示的结构、比例、大小等,均仅用以配合说明书所揭示的内容,以供熟悉此技术的人士了解与阅读,并非用以限定本发明可实施的限定条件,故不具技术上的实质意义,任何结构的修饰、比例关系的改变或大小的调整,在不影响本发明所能产生的功效及所能达成的目的下,均应仍落在本发明所揭示的技术内容得能涵盖的范围内。
实施例
参见图1至图3,
1自适应t分布策略
t分布是数理统计和概率论理论中一种重要的分布类型,在很多领域都得到了广泛应用。该种分布的参数自由度决定着曲线形态:当自由度越大,曲线的形态越加高耸,同时越接近标准正态分布曲线;若参数取无限大值,则该分布就是高斯分布,即当自由度越小,曲线的形态会越加低平;若参数取值为1时,则该分布就是柯西分布,即
2动态惯性权重因子
蝴蝶的位置更新对种群的寻优取决定性作用,在蝴蝶算法初期,种群中的个体模式集中在适应度值较低的个体上,需要种群尽快能够在整个空间进行“搜索”。但是到了算法后期,个体模式朝高适应度的个体集中,会造成种群的“开发”能力有限,进而出现早熟收敛,无法跳出局部极值。因此,在蝴蝶的位置更新中引入一种动态惯性权重因子,以便在算法的前期能够尽快达到较优区域,在算法的后期对较优区域进一步开发以尽快找到最优解。
首先给出进化离散度的定义,以便描述种群整体适应度的变化,设o(d)为第d代种群与第d-1代种群的适应度值标准差的比值:
在神经网络理论中,神经元的激活函数常采用Sigmoid函数来构造,其定义:
图1给出了Sigmoid函数图像,为一连续光滑且严格单调曲线,点(0,0.5)为其中心对称点,当自变量x在[-10,10]之间时函数值域变化,且变化范围为(0,1)。该函数在线性与非线性之间具有良好的平衡性,因此常作为阈值函数使用。结合进化离散度o(d)和Sigmoid函数,提出一种非线性动态惯性权重因子:
其中ωmax为最大惯性权重,一般取值为0.9,ωmin为最小惯性权重,一般取值为0.4,
d、Dmax表示当前迭代次数与最大进化代数。o(d)为种群进化离散度参数,e是阻尼因子,一般取值为[0,1]。
3动态自适应t分布变异策略
依据公式(6)生成的非线性动态惯性权重因子,在常规蝴蝶算法中进行t分布变异操作,其变异程度由非线性动态惯性权重因子ω进行控制,最终通过引入这种动态自适应t分布变异策略,对算法中蝴蝶局部搜索阶段的移动公式改进如下:
其中,tdis(d)表示当迭代次数为d时,以其作为参数自由度的t分布,公式(7)在蝴蝶当前位置基础上增加了自适应t分布的变异干扰项ω·tdis(d)。在算法的运行初期,由于d 值较小,此时该干扰项可以视为柯西变异,使得算法保持了蝴蝶种群的多样性,因而可以最大程度的发挥干扰项的作用,使其具有良好的全局探索能力,有助于算法避免陷入局部最优;但在算法的运行后期,变异干扰项变为高斯分布变异,这将加快其收敛速度及提高搜索的精度,算法具有较强的局部搜索能力,同时变异干扰项的作用随着迭代次数的增加不断减小。
4算法流程
一种基于改进的蝴蝶优化算法,包括有如下步骤:
S1、参数初始化设置:
初始化蝴蝶优化算法的各项基本参数;
S2、初始化种群,并计算各个解的是适应度值:
初始化种群,根据适应度函数计算出每只蝴蝶在其解空间中的适应度值;
S3、计算每只蝴蝶发出香味的浓度函数:
利用随机函数生成数决策该蝴蝶执行全局搜索还是局部搜索,以尽可能减少外部环境的影响;
S4、判断P>rand:
若P>rand,则利用式(2)进行全局搜索:
低浓度香味的蝴蝶个体受浓度最佳的蝴蝶吸引并向其移动;
若P<rand,则按照式(7)进行局部搜索:
引入动态自适应t分布变异策略:
给出进化离散度的定义,以便描述种群整体适应度的变化,设o(d)为第d代种群与第 d-1代种群的适应度值标准差的比值:
在神经网络理论中,神经元的激活函数采用Sigmoid函数来构造,其定义:
点(0,0.5)为其中心对称点,当自变量x在[-10,10]之间时函数值域变化,且变化范围为(0,1),函数在线性与非线性之间具有良好的平衡性,作为阈值函数使用,结合进化离散度o(d)和Sigmoid函数,提出一种非线性动态惯性权重因子:
依据公式(6)生成的非线性动态惯性权重因子,在常规蝴蝶算法中进行t分布变异操作,其变异程度由非线性动态惯性权重因子ω进行控制,最终通过引入这种动态自适应t分布变异策略,对算法中蝴蝶局部搜索阶段的移动公式改进如下:
其中,tdis(d)表示当迭代次数为d时,以其作为参数自由度的t分布,公式(7)在蝴蝶当前位置基础上增加了自适应t分布的变异干扰项ω·tdis(d);
S5、更新蝴蝶个体和全局最优解;
S6、越界处理;
S7、判断是否满足迭代次数要求或精度要求,则停止计算并输出最优解及适应度值,否则转向S3继续循环。
5机器人路径规划
假设机器人工作环境为二维空间,在空间中分布着有限数量的静态障碍物(凸多边形),机器人路径规划的任务是在起始点和终止点之间找到一条最短的、较平滑的、且避开所有障碍物的路径。本发明采用导航点模型来构建机器人工作环境,如图3所示。
6路径编码
起点B与终点E的连线为X'轴构建坐标系BX'Y',如图3所示。然后将坐标系OXY中的点变换到BX'Y'中,变换公式如下:
其中,(xb,yb)是坐标系OXY中起点B的坐标,(x',y')是点(x,y)在坐标系BX'Y'中对应的点,θ是X轴与直线BE的夹角。
如图3所示,用m条平行簇l1,l2,…,lm将BE平均为m+1段,每每相邻的两条平时直线的距离设为Δl=||BE||/(m+1)。通过每条直线上面随机产生的一个点,可以构造一条机器人移动的完整路径(B,P1,P2,…,Pm,E),将这条路径上所有的节点生成为一个蜉蝣个体的编码序列。这样就通过式(8)将一个机器人移动路径的规划问题转换为对一个集合点的变量的优化问题。
7适应度函数
路径规划主要考虑长度、安全度和光滑度3个指标.
(1)长度指标
设起点B和终点E的坐标为B(xb,yb)、E(xe,ye),任意路径节点坐标设为Pi(xi,yi),i的取值为[1,m]。设路径长度为fdistance,将其标准化表示如下式:
(2)安全度指标
为避免机器人与障碍物发生碰撞且路径更加平滑,本发明采用准则是在三次样条插值节点处与障碍物无碰撞。假设已知m个路径结点的坐标(x1,y1),(x2,y2),…,(xm,ym)以及起点坐标(xb,yb)和终点坐标(xe,ye)。通过三次样条差值分别得到d个插值点,其横坐标为 (x1,x2,…,xd),其纵坐标为(y1,y2,…,yd)。需要所生成的d个插值点是否有碰撞情况发生,本发明将d设置为100,设路径安全度指标为fsecure,将其如下式表示:
上式中η表示安全因子,取值为100,Dj,k表示第j个插值点到第k个障碍物中心距离, R(k)表示第k个障碍物的半径,H表示路径中障碍物的数量。
(3)光滑度指标
设路径节点为(B,P1,P2,…,Pm,E),相邻路段之间的夹角表示为ψi,本发明建立的路径平滑度指标为:
ψi越小表示路径光滑程度越优。
将以上指标进行融合后得到机器人路径规划的综合指标S:
S=fdistance+fsecure+fflatness (12)
8基于改进的蝴蝶优化算法的机器人多目标路径规划,包括有如下步骤:
1)进行坐标系变换,且利用式(8)对起点、终点、障碍物位置的坐标进行变换:
起点B与终点E的连线为X'轴构建坐标系BX'Y',然后将坐标系OXY中的点变换到BX'Y'中,变换公式如下:
其中,(xb,yb)是坐标系OXY中起点B的坐标,(x',y')是点(x,y)在坐标系BX'Y'中对应的点,θ是X轴与直线BE的夹角;
2)建立m条平行簇l1,l2,…,lm将BE平均为m+1段,每每相邻的两条平时直线的距离设为Δl=||BE||/(m+1);
3)初始化参数D、c、I、a等;
5)利用适应度函数式(12),对蝴蝶进行评估,获得蝴蝶的适应度S;
长度指标
设起点B和终点E的坐标为B(xb,yb)、E(xe,ye),任意路径节点坐标设为Pi(xi,yi),i的取值为[1,m],设路径长度为fdistance,将其标准化表示如下式:
安全度指标
设机器人在三次样条插值节点处与障碍物无碰撞,假设已知m个路径结点的坐标(x1,y1),(x2,y2),…,(xm,ym)以及起点坐标(xb,yb)和终点坐标(xe,ye),通过三次样条差值分别得到d个插值点,其横坐标为(x1,x2,…,xd),其纵坐标为(y1,y2,…,yd),需要所生成的d个插值点是否有碰撞情况发生,将d设置为100,设路径安全度指标为fsecure,将其如下式表示:
上式中η表示安全因子,取值为100,Dj,k表示第j个插值点到第k个障碍物中心距离, R(k)表示第k个障碍物的半径,H表示路径中障碍物的数量;
光滑度指标
设路径节点为(B,P1,P2,…,Pm,E),相邻路段之间的夹角表示为ψi,建立的路径平滑度指标为:
ψi越小表示路径光滑程度越优;
将以上指标进行融合后得到机器人路径规划的综合指标S:
S=fdistance+fsecure+fflatness (12);
6)利用公式(1)计算每只蝴蝶发出香味的浓度函数,利用随机函数生成数决策该蝴蝶执行全局搜索还是局部搜索,减少外部环境的影响:
全局搜索:
利用公式(2)进行全局搜索,低浓度香味的蝴蝶个体受浓度最佳的蝴蝶吸引并向其移动;
局部搜索:
引入动态自适应t分布变异策略利用公式(7)进行局部搜索,蝴蝶个体进行随机搜索;
7)更新蝴蝶个体和全局最优解;
8)越接处理;
9)判断是否满足迭代次数要求或精度要求,若是满足进入11),否则返回5)继续循环;
10)输出最优个体适应度值。
本发明方案所公开的技术手段不仅限于上述技术手段所公开的技术手段,还包括由以上技术特征等同替换所组成的技术方案。本发明的未尽事宜,属于本领域技术人员的公知常识。
Claims (5)
1.一种基于改进的蝴蝶优化算法,其特征在于:包括有如下步骤:
S1、参数初始化设置:
初始化蝴蝶优化算法的各项基本参数;
S2、初始化种群,并计算各个解的是适应度值:
初始化种群,根据适应度函数计算出每只蝴蝶在其解空间中的适应度值;
S3、计算每只蝴蝶发出香味的浓度函数:
利用随机函数生成数决策该蝴蝶执行全局搜索还是局部搜索,以尽可能减少外部环境的影响;
S4、判断P>rand:
若P>rand,则利用式(2)进行全局搜索:
低浓度香味的蝴蝶个体受浓度最佳的蝴蝶吸引并向其移动;
若P<rand,则按照式(7)进行局部搜索:
引入动态自适应t分布变异策略:
给出进化离散度的定义,以便描述种群整体适应度的变化,设o(d)为第d代种群与第d-1代种群的适应度值标准差的比值:
在神经网络理论中,神经元的激活函数采用Sigmoid函数来构造,其定义:
点(0,0.5)为其中心对称点,当自变量x在[-10,10]之间时函数值域变化,且变化范围为(0,1),函数在线性与非线性之间具有良好的平衡性,作为阈值函数使用,结合进化离散度o(d)和Sigmoid函数,提出一种非线性动态惯性权重因子:
依据公式(6)生成的非线性动态惯性权重因子,在常规蝴蝶算法中进行t分布变异操作,其变异程度由非线性动态惯性权重因子ω进行控制,最终通过引入这种动态自适应t分布变异策略,对算法中蝴蝶局部搜索阶段的移动公式改进如下:
其中,tdis(d)表示当迭代次数为d时,以其作为参数自由度的t分布,公式(7)在蝴蝶当前位置基础上增加了自适应t分布的变异干扰项ω·tdis(d);
S5、更新蝴蝶个体和全局最优解;
S6、越界处理;
S7、判断是否满足迭代次数要求或精度要求,则停止计算并输出最优解及适应度值,否则转向S3继续循环。
2.根据权利要求1所述的一种基于改进的蝴蝶优化算法,其特征在于:在公式(6)中,其中ωmax为最大惯性权重,取值为0.9,ωmin为最小惯性权重,取值为0.4。
3.根据权利要求2所述的一种基于改进的蝴蝶优化算法,其特征在于:在公式(6)中,d、Dmax分别表示当前迭代次数与最大进化代数,o(d)为种群进化离散度参数,e是阻尼因子,取值为[0,1]。
4.根据权利要求3所述的一种基于改进的蝴蝶优化算法,其特征在于:在步骤S4中,在算法的运行初期,由于d值较小,变异干扰项ω·tdis(d)视为柯西变异,在算法的运行后期,变异干扰项ω·tdis(d)视为高斯分布变异。
5.一种基于改进的蝴蝶优化算法的机器人多目标路径规划,其特征在于:包括有如下步骤:
1)进行坐标系变换,且利用式(8)对起点、终点、障碍物位置的坐标进行变换:
起点B与终点E的连线为X'轴构建坐标系BX'Y',然后将坐标系OXY中的点变换到BX'Y'中,变换公式如下:
其中,(xb,yb)是坐标系OXY中起点B的坐标,(x',y')是点(x,y)在坐标系BX'Y'中对应的点,θ是X轴与直线BE的夹角;
2)建立m条平行簇l1,l2,…,lm将BE平均为m+1段,每每相邻的两条平时直线的距离设为Δl=||BE||/(m+1);
3)初始化参数;
5)利用适应度函数式(12),对蝴蝶进行评估,获得蝴蝶的适应度S;
长度指标
设起点B和终点E的坐标为B(xb,yb)、E(xe,ye),任意路径节点坐标设为Pi(xi,yi),i的取值为[1,m],设路径长度为fdistance,将其标准化表示如下式:
安全度指标
设机器人在三次样条插值节点处与障碍物无碰撞,假设已知m个路径结点的坐标(x1,y1),(x2,y2),…,(xm,ym)以及起点坐标(xb,yb)和终点坐标(xe,ye),通过三次样条差值分别得到d个插值点,其横坐标为(x1,x2,…,xd),其纵坐标为(y1,y2,…,yd),需要所生成的d个插值点是否有碰撞情况发生,将d设置为100,设路径安全度指标为fsecure,将其如下式表示:
上式中η表示安全因子,取值为100,Dj,k表示第j个插值点到第k个障碍物中心距离,R(k)表示第k个障碍物的半径,H表示路径中障碍物的数量;
光滑度指标
设路径节点为(B,P1,P2,…,Pm,E),相邻路段之间的夹角表示为ψi,建立的路径平滑度指标为:
ψi越小表示路径光滑程度越优;
将以上指标进行融合后得到机器人路径规划的综合指标S:
S=fdistance+fsecure+fflatness (12);
6)利用公式(1)计算每只蝴蝶发出香味的浓度函数,利用随机函数生成数决策该蝴蝶执行全局搜索还是局部搜索,减少外部环境的影响:
全局搜索:
利用公式(2)进行全局搜索,低浓度香味的蝴蝶个体受浓度最佳的蝴蝶吸引并向其移动;
局部搜索:
引入动态自适应t分布变异策略利用公式(7)进行局部搜索,蝴蝶个体进行随机搜索;
7)更新蝴蝶个体和全局最优解;
8)越接处理;
9)判断是否满足迭代次数要求或精度要求,若是满足进入11),否则返回5)继续循环;
10)输出最优个体适应度值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110996730.1A CN113552891A (zh) | 2021-08-27 | 2021-08-27 | 一种基于改进的蝴蝶优化算法的机器人多目标路径规划 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110996730.1A CN113552891A (zh) | 2021-08-27 | 2021-08-27 | 一种基于改进的蝴蝶优化算法的机器人多目标路径规划 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113552891A true CN113552891A (zh) | 2021-10-26 |
Family
ID=78106124
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110996730.1A Withdrawn CN113552891A (zh) | 2021-08-27 | 2021-08-27 | 一种基于改进的蝴蝶优化算法的机器人多目标路径规划 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113552891A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114035586A (zh) * | 2021-11-22 | 2022-02-11 | 江苏科技大学 | 改进蚁群算法和动态窗口的车间agv小车路径规划方法 |
CN115062733A (zh) * | 2022-08-19 | 2022-09-16 | 江西派源科技有限公司 | 基于经验模态分解与蝴蝶算法的变压器声纹故障诊断方法 |
CN115357028A (zh) * | 2022-09-23 | 2022-11-18 | 桂林电子科技大学 | 一种排序野马优化算法的移动机器人路径规划方法 |
CN115435787A (zh) * | 2022-07-28 | 2022-12-06 | 南京航空航天大学 | 一种基于改进蝴蝶算法的无人机三维路径规划方法及*** |
CN115455721A (zh) * | 2022-09-23 | 2022-12-09 | 国网安徽省电力有限公司淮北供电公司 | 一种并网三相逆变器控制参数整定方法及*** |
CN115657694A (zh) * | 2022-12-29 | 2023-01-31 | 中国人民解放军陆军装甲兵学院 | 一种基于空地信息交互的无人车路径规划方法 |
CN118170033A (zh) * | 2024-05-14 | 2024-06-11 | 青岛理工大学 | 双足机器人的平地步态优化方法及*** |
-
2021
- 2021-08-27 CN CN202110996730.1A patent/CN113552891A/zh not_active Withdrawn
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114035586A (zh) * | 2021-11-22 | 2022-02-11 | 江苏科技大学 | 改进蚁群算法和动态窗口的车间agv小车路径规划方法 |
CN114035586B (zh) * | 2021-11-22 | 2024-03-29 | 江苏科技大学 | 改进蚁群算法和动态窗口的车间agv小车路径规划方法 |
CN115435787A (zh) * | 2022-07-28 | 2022-12-06 | 南京航空航天大学 | 一种基于改进蝴蝶算法的无人机三维路径规划方法及*** |
CN115435787B (zh) * | 2022-07-28 | 2024-05-28 | 南京航空航天大学 | 一种基于改进蝴蝶算法的无人机三维路径规划方法及*** |
CN115062733A (zh) * | 2022-08-19 | 2022-09-16 | 江西派源科技有限公司 | 基于经验模态分解与蝴蝶算法的变压器声纹故障诊断方法 |
CN115357028A (zh) * | 2022-09-23 | 2022-11-18 | 桂林电子科技大学 | 一种排序野马优化算法的移动机器人路径规划方法 |
CN115455721A (zh) * | 2022-09-23 | 2022-12-09 | 国网安徽省电力有限公司淮北供电公司 | 一种并网三相逆变器控制参数整定方法及*** |
CN115357028B (zh) * | 2022-09-23 | 2024-05-10 | 桂林电子科技大学 | 一种排序野马优化算法的移动机器人路径规划方法 |
CN115657694A (zh) * | 2022-12-29 | 2023-01-31 | 中国人民解放军陆军装甲兵学院 | 一种基于空地信息交互的无人车路径规划方法 |
CN118170033A (zh) * | 2024-05-14 | 2024-06-11 | 青岛理工大学 | 双足机器人的平地步态优化方法及*** |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113552891A (zh) | 一种基于改进的蝴蝶优化算法的机器人多目标路径规划 | |
Nazarahari et al. | Multi-objective multi-robot path planning in continuous environment using an enhanced genetic algorithm | |
Das et al. | Intelligent-based multi-robot path planning inspired by improved classical Q-learning and improved particle swarm optimization with perturbed velocity | |
Das et al. | A hybrid improved PSO-DV algorithm for multi-robot path planning in a clutter environment | |
Zhang et al. | Path planning of mobile robot based on hybrid multi-objective bare bones particle swarm optimization with differential evolution | |
CN109945881B (zh) | 一种蚁群算法的移动机器人路径规划方法 | |
CN113485451A (zh) | 一种基于改进的蜉蝣优化算法的机器人多目标路径规划 | |
CN111381600B (zh) | 一种基于粒子群算法的uuv路径规划方法 | |
Chakraborty et al. | Cooperative multi-robot path planning using differential evolution | |
Wang et al. | Pattern-rl: Multi-robot cooperative pattern formation via deep reinforcement learning | |
CN110530373B (zh) | 一种机器人路径规划方法、控制器及*** | |
Elmi et al. | Multi-objective grasshopper optimization algorithm for robot path planning in static environments | |
CN116834037B (zh) | 基于动态多目标优化的采摘机械臂轨迹规划方法及装置 | |
CN113791613A (zh) | 改进的黏菌优化算法及基于该算法的机器人多目标路径规划 | |
Al Dabooni et al. | Heuristic dynamic programming for mobile robot path planning based on Dyna approach | |
Belmonte-Baeza et al. | Meta reinforcement learning for optimal design of legged robots | |
Hafez et al. | Topological Q-learning with internally guided exploration for mobile robot navigation | |
CN115983130A (zh) | 一种基于改进粒子群优化算法的全局最优解搜索方法 | |
Tian et al. | Unsupervised discovery of transitional skills for deep reinforcement learning | |
Yunqiang et al. | Research on multi-objective path planning of a robot based on artificial potential field method | |
Na et al. | Gradient eigendecomposition invariance biogeography-based optimization for mobile robot path planning | |
Raiesdana | A hybrid method for industrial robot navigation | |
Yang et al. | PMDRL: Pareto-front-based multi-objective deep reinforcement learning | |
Zhang et al. | Multi-task actor-critic with knowledge transfer via a shared critic | |
Tang et al. | Reinforcement learning for robots path planning with rule-based shallow-trial |
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: 20211026 |