CN110244733B - 一种基于改进蚁群算法的移动机器人路径规划方法 - Google Patents
一种基于改进蚁群算法的移动机器人路径规划方法 Download PDFInfo
- Publication number
- CN110244733B CN110244733B CN201910535562.9A CN201910535562A CN110244733B CN 110244733 B CN110244733 B CN 110244733B CN 201910535562 A CN201910535562 A CN 201910535562A CN 110244733 B CN110244733 B CN 110244733B
- Authority
- CN
- China
- Prior art keywords
- ant
- grid
- ants
- search
- representing
- 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
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 67
- 238000000034 method Methods 0.000 title claims abstract description 44
- 241000257303 Hymenoptera Species 0.000 claims abstract description 63
- 230000007704 transition Effects 0.000 claims abstract description 35
- 238000004364 calculation method Methods 0.000 claims abstract description 26
- 239000003016 pheromone Substances 0.000 claims description 31
- 230000000007 visual effect Effects 0.000 claims description 14
- 238000009826 distribution Methods 0.000 claims description 7
- 230000007613 environmental effect Effects 0.000 claims description 3
- 230000008569 process Effects 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 description 17
- 238000004458 analytical method Methods 0.000 description 5
- 230000009286 beneficial effect Effects 0.000 description 5
- 238000004088 simulation Methods 0.000 description 5
- 238000005457 optimization Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 241001465754 Metazoa Species 0.000 description 2
- 230000004888 barrier function Effects 0.000 description 2
- 230000009193 crawling Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 244000062645 predators Species 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/0217—Control of position or course in two dimensions specially adapted to land vehicles with means for defining a desired trajectory in accordance with energy consumption, time reduction or distance reduction criteria
-
- 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)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明提供了一种基于改进蚁群算法的移动机器人路径规划方法,提出了一种16个搜索方向24个邻域节点的改进蚁群算法,使蚂蚁在机器人规划的路径搜索的过程中具有更多的方向选择和更广的搜索范围,并设计了2种基于向量夹角的启发信息计算方法;另外,在转移概率部分通过引入转移概率的控制参数来控制蚂蚁在栅格地图中的有效搜索范围,可以搜索到更优的全局最短路径,可以提高算法的搜索的有效性,降低搜索的随机性。本发明通过改进的蚁群算法能够搜索到更优的全局最短路径,可以提高蚁群算法搜索的有效性,降低搜索的随机性,为机器人的路径规划提供良好的条件。
Description
技术领域
本发明属于路径规划技术领域,尤其涉及一种基于改进蚁群算法的移动机器人路径规划方法。
背景技术
路径规划是移动机器人导航研究的重要内容之一,指的是从起点到终点为机器人规划一条有效路径,并保证机器人与路径上的障碍物不发生任何碰撞。
随着智能机器人技术的发展,使智能机器人在军事、航空、智能制造、智能仓储等领域的应用备受关注。实现智能移动机器人从起点到终点的过程关键在于路径规划算法,因此,精确、有效的路径规划算法对于移动机器人导航技术的研究起着至关重要的作用。由于蚁群算法不受编码规则限制,具有较强的鲁棒性,易于与其他方法相结合的优势,目前蚁群算法已经成为了机器人路径规划领域的一个研究热点。
在栅格地图场景下,传统的蚂蚁的搜索方式主要有2种:4个搜索方向4个领域节点和8个搜索方向8个领域节点,且均存在算法易陷入局部最优,无法找到有效地最优解的问题,本发明提出了一种基于改进蚁群算法的移动机器人路径规划方法解决了上述的问题。
发明内容
针对现有技术中的上述不足,本发明提供的一种基于改进蚁群算法的移动机器人路径规划方法解决了传统蚁群算法易陷入局部最优,无法找到最优解的问题。
为了达到以上目的,本发明采用的技术方案为:
本方案提供了一种基于改进蚁群算法的移动机器人路径规划方法,包括如下步骤:
S1、根据环境障碍物分布特点建立栅格地图,初始化栅格信息素,并设置初始参数和确定移动机器人的起点和终点;
S2、初始化禁忌表,将蚁群算法的蚂蚁参数设为M只蚂蚁,并将移动机器人的起点和终点设为M只蚂蚁的起点和终点;
S3、根据蚂蚁的移动规则计算栅格地图中领域节点的启发信息值;
S4、根据所述启发信息值计算栅格的转移概率,并根据所述转移概率的控制参数确定蚂蚁下一时刻到达的位置,更新禁忌表;
S5、将所有蚂蚁完成一次路径搜索,并更新成功抵达终点的蚂蚁走过的路径上的栅格信息素;
S6、判断是否完成预设的迭代次数,若是,则输出最优路径,从而完成移动机器人的路径规划,反之,则返回步骤S2。
进一步地,所述步骤S1中设置初始参数包括设置栅格的初始信息素浓度、蚁群算法的迭代次数、蚂蚁投放数量、信息启发式因子、期望启发因子以及栅格信息素的挥发系数。
再进一步地,所述步骤S3中的蚂蚁的移动规则包括:蚂蚁的搜索方向以及蚂蚁的视野范围。
再进一步地,所述蚂蚁的搜索方向为:搜索栅格地图中的16个方向和24个领域节点,且蚂蚁的转角为22.5°。
再进一步地,所述蚂蚁的视野范围为:将所述24个领域节点分为1级视野范围栅格和2级视野范围栅格,其中,
所述1级视野范围栅格满足的条件为:
所述2级视野范围栅格满足的条件为:
其中,L表示栅格的边长,dik表示蚂蚁当前位置栅格中心到领域栅格中心的距离。
再进一步地,针对1级视野范围栅格和2级视野范围栅格的启发信息值相同时,蚂蚁选择1级视野范围内的领域栅格作为下一时刻到达的位置。
再进一步地,所述步骤S3中的启发信息值ηik的计算方法包括如下两种方法:
第一种计算方法:
其中,θik表示蚂蚁搜索方向与理论最佳搜索方向之间的夹角,Dik表示蚂蚁下一时刻的移动方向,DiE表示蚂蚁理论上的最佳搜索方向,xk表示蚂蚁下一时刻到达位置的横坐标,xi表示蚂蚁当前位置的横坐标,xE表示终点的横坐标,yk表示蚂蚁下一时刻到达位置的纵坐标,yi蚂蚁当前位置的纵坐标,yE表示终点的纵坐标;
第二种计算方法:
其中,θik表示蚂蚁搜索方向与理论最佳搜索方向之间的夹角,xk表示蚂蚁下一时刻到达位置的横坐标,xi表示蚂蚁当前位置的横坐标,xE表示终点的横坐标,yk表示蚂蚁下一时刻到达位置的纵坐标,yi蚂蚁当前位置的纵坐标,yE表示终点的纵坐标。
再进一步地,所述步骤S4包括如下步骤:
S401、根据所述启发信息值计算栅格的转移概率;
S402、生成一个随机数ω,且ω∈[0,1];
S403、判断所述转移概率的控制参数σ是否大于所述随机数ω,若是,则进入步骤S304,反之,则进入步骤S305,其中,σ∈[0,1];
S404、通过轮盘赌策略选择邻域栅格作为蚂蚁下一时刻到达位置,更新禁忌表;
S405、选择转移概率最大的邻域栅格作为蚂蚁下一时刻到达位置,更新禁忌表。
其中,τik表示t时刻路段(i,k)上的信息素量,ηik表示蚂蚁m从节点i到节点k的启发信息值,τis表示蚂蚁m从节点i到节点s的信息素浓度,α表示信息启发式因子,β为期望启发式因子,allowedm表示蚂蚁m下一步允许选择的节点集合。
再进一步地,所述步骤S5中栅格信息素τik(t+1)的表达式如下:
τik(t+1)=(1-ρ)τik(t)+Δτik(t)
本发明的有益效果:
本发明提供了一种基于改进蚁群算法的移动机器人路径规划方法,提出了一种16个搜索方向24个邻域节点的改进蚁群算法,使蚂蚁在机器人规划的路径搜索的过程中具有更多的方向选择和更广的搜索范围,并设计了2种基于向量夹角的启发信息计算方法;另外,在转移概率部分通过引入转移概率的控制参数来控制蚂蚁在栅格地图中的有效搜索范围,可以搜索到更优的全局最短路径,可以提高算法的搜索的有效性,降低搜索的随机性。本发明通过改进的蚁群算法能够搜索到更优的全局最短路径,可以提高蚁群算法搜索的有效性,降低搜索的随机性,为机器人的路径规划提供良好的条件。
附图说明
图1为本发明的方法流程图。
图2为本实施例中的栅格环境地图的示意图。
图3为本实施例中蚂蚁搜索方向的示意图。
图4为本实施例中蚂蚁搜索领域点的示意图。
图5为本实施例中蚂蚁采用不同搜索方式规划的路径的示意图。
图6为本实施例中蚂蚁搜索方向与理论最佳搜索方向之间的夹角示意图。
图7为本实施例中使用一种启发式信息计算方法的函数图像示意图。
图8为本实施例中使用另一种启发式信息计算方法的函数图像示意图。
图9为本实施例中使用两种计算方法的函数图。
图10为本实施例中启发信息值大于等于0.5时的搜索角度范围示意图。
图11为本实施例中蚂蚁的视野范围示意图。
图12为本实施例中蚂蚁查看2级视野范围内栅格分析图。
图13为本实施例中2级视野范围内栅格不可见的示意图。
图14为本实施例中栅格中心位于同一连线的情况下的视野范围内栅格分析图。
图15为本实施例中使用两种启发式信息计算方法得到的蚂蚁搜索路径轨迹示意图。
图16为本实施例中转移概率参数取不同值时的蚂蚁爬行路径轨迹示意图。
图17为本实施例中转移概率参数取不同值时的仿真结果图。
图18为本实施例中转移概率参数取不同值时10次仿真结果平均值的示意图。
图19为本实施例使用本发明和传统蚁群算法路径长度收敛对比图。
图20为本实施例使用本发明和传统蚁群算法搜索到的最短路径的示意图。
图21为本实施例中使用本发明搜索最短路径过程中的蚂蚁搜索轨迹图。
具体实施方式
下面对本发明的具体实施方式进行描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
实施例
本发明提供了一种基于改进蚁群算法的移动机器人路径规划方法,提出了一种16个搜索方向24个邻域节点的改进蚁群算法,使蚂蚁在机器人规划的路径搜索的过程中具有更多的方向选择和更广的搜索范围,并设计了2种基于向量夹角的启发信息计算方法;另外,在转移概率部分通过引入转移概率的控制参数来控制蚂蚁在栅格地图中的有效搜索范围,可以搜索到更优的全局最短路径,可以提高算法的搜索的有效性,降低搜索的随机性。
如图1所示,本发明提供了一种基于改进蚁群算法的移动机器人路径规划方法,其实现方法如下:
S1、根据环境障碍物分布特点建立栅格地图,初始化栅格信息素,并设置初始参数和确定移动机器人的起点和终点,所述设置初始参数包括设置栅格的初始信息素浓度、蚁群算法的迭代次数、蚂蚁投放数量、信息启发式因子、期望启发因子以及栅格信息素的挥发系数。在具体实施例中,使用栅格法建立环境地图可以有效表示环境布局特征,首先,对环境中的障碍物进行膨胀化处理,即只要栅格被障碍物触及,则该栅格按障碍物处理。如图2所示,图2(a)为障碍物实际形状,图2(b)表示障碍物在栅格地图中膨胀化处理之后的形状,移动机器人则在自由区域的栅格中心之间移动,另外,栅格尺寸可以根据机器人尺寸或者工作精度要求确定;
S2、初始化禁忌表,将蚁群算法的蚂蚁参数设为M只蚂蚁,并将移动机器人的起点和终点设为M只蚂蚁的起点和终点,在具体实施例中,蚂蚁为蚁群算法里面的参数,其执行的是寻路,到达指定迭代次数后会给出一条最短路径,移动机器人按照寻找到的路径移动;
S3、根据蚂蚁的移动规则计算栅格地图中领域节点的启发信息值,所述蚂蚁的移动规则包括:蚂蚁的搜索方向以及蚂蚁的视野范围,所述蚂蚁的搜索方向为:搜索栅格地图中的16个方向和24个领域节点,且蚂蚁的转角为22.5°,所述蚂蚁的视野范围为:将所述24个领域节点分为1级视野范围栅格和2级视野范围栅格,其中,
所述1级视野范围栅格满足的条件为:
所述2级视野范围栅格满足的条件为:
其中,L表示栅格的边长,dik表示蚂蚁当前位置栅格中心到领域栅格中心的距离。
现实生活中蚂蚁可以沿任意方向移动搜索食物源,但在基于栅格地图的蚁群算法中,由于栅格地图的特点限制,蚂蚁的有效搜索方向是有限的,传统的蚂蚁搜索方式有4个搜索方向的搜索方式(图3(a))和8个搜索方向的搜索方式(图3(b))。在图3(a)所示的搜索方式中蚂蚁只能沿水平或者垂直方向移动,在图3(b)所示的搜索方式中是在前者的基础之上增加了对角线方向的移动。为了增加蚂蚁移动方向的数量,本发明将蚂蚁的搜索方向扩展到16个,如图3(c)所示,将蚂蚁的最小转角细化到了22.5°,可以有效地扩大蚂蚁的搜索范围,有利于蚂蚁搜索到更优的路径。对于动物来讲,长远的视野范围可以帮助它们提前发现目标、规划路线,进行有效地捕获猎物或者躲避捕食者,同样的道理,蚂蚁在选择下一步节点前,长远的视野范围可以帮助蚂蚁做出更好地路径规划或者提前发现目标点。在4个搜索方向的搜索方式中,蚂蚁每次只能查看当前位置的4个领域节点,如图4(a)所示;在8个搜索方向的搜索方式中,蚂蚁每次可以查看当前节点的8个领域节点,如图4(b)所示;而在发明提出的16个搜索方向的搜索方式中,蚂蚁每次可以查看当前节点的24个领域节点,如图4(c)所示。细化的搜索方向和大范围的搜索邻域有利于蚂蚁找到更短的路径,图5(a)和图5(b)为传统蚂蚁搜索方式搜索到的路径,图5(c)为本发明提出的蚂蚁搜索方式搜索到的路径。假设单元栅格的边长为1,可以计算出蚂蚁采用不同的搜索方式搜索到的最优路径长度和找到最优路径的搜索次数,如表1所示,通过表1的数据可知,本发明提出的蚂蚁搜索方式可以通过更少的搜索次数找到更短的可行路径。
表1
搜索方式 | 规划的路径长度 | 搜索次数 |
4个搜索方向 | 3.00 | 3 |
8个搜索方向 | 2.41 | 2 |
16个搜索方向 | 2.24 | 1 |
在具体实施例中,标准蚁群算法中往往采用当前节点i的领域节点k到终点E之间欧式距离的倒数作为节点k的启发信息,本发明基于向量夹角提出了2种新的启发信息计算方法。已知蚂蚁当前位置i的坐标为(xi,yi),蚂蚁当前位置i的领域节点k的坐标为(xk,yk),蚂蚁的目标点E的坐标为(xE,yE),因此,蚂蚁的理论最佳搜索方向可以表示为DiE=(xE-xi,yE-yi),蚂蚁下一时刻的移动方向可以表示为Dik=(xk-xi,yk-yi),通过计算可以得到向量Dik与向量DiE之间的夹角的弧度值θ,且θ∈[0,π]。如图6所示,蚂蚁当前位置的任意一个领域节点都有一个相应的弧度值θ,因此以弧度值θ为变量设计了2种启发信息计算方法:
第一种计算方法:
其中,θik表示蚂蚁搜索方向与理论最佳搜索方向之间的夹角,Dik表示蚂蚁下一时刻的移动方向,DiE表示蚂蚁理论上的最佳搜索方向,xk表示蚂蚁下一时刻到达位置的横坐标,xi表示蚂蚁当前位置的横坐标,xE表示终点的横坐标,yk表示蚂蚁下一时刻到达位置的纵坐标,yi蚂蚁当前位置的纵坐标,yE表示终点的纵坐标。
第一种方法在区间[0,π]上的函数图像如图7所示,第一种方法可知,当节点k=E时函数无解,解决办法可以将特殊点E的启发信息设为一个较大常量;
第二种计算方法:
其中,θik表示蚂蚁搜索方向与理论最佳搜索方向之间的夹角,xk表示蚂蚁下一时刻到达位置的横坐标,xi表示蚂蚁当前位置的横坐标,xE表示终点的横坐标,yk表示蚂蚁下一时刻到达位置的纵坐标,yi蚂蚁当前位置的纵坐标,yE表示终点的纵坐标。
第二种方法在区间[0,π]上的函数图像如图8所示。启发信息值是蚂蚁选择领域节点的重要判断依据,邻域节点的启发信息值越大,该点被选择的可能性就越高。为了便于分析,用λ1和η1分别表示第一种方法的夹角和启发信息值,λ2和η2分别表示第二种方法的夹角和启发信息值,此处的λ表示度数,角度与弧度的转化关系为在图9中,两函数在A存在一个交点,在交点位置λ1=λ2=148.1°,可以发现在A之前,当η1=η2时,λ1>λ2,也就是说在相同的启发信息值的情况下,采用第一种方法搜索的角度范围更大,有利于指导蚂蚁搜索更大的范围,采用第二种方法搜索的范围则更加集中,有利于指导蚂蚁朝着最佳搜索方向移动。例如:当启发信息值η1=η2≥0.5时,采用第一种方法得到的蚂蚁搜索范围如图10(a)所示,采用第二种方法得到的蚂蚁搜索范围如图10(b)所示。因此,可以总结为:采用第一种方法体现搜索的全局性,采用第二种方法体现搜索的方向性。
在具体实施例中,可以将蚂蚁视野范围内的邻居栅格划分为:1级视野范围栅格和2级视野范围栅格,如图11所示。蚂蚁从当前位置可以直接查看1级视野范围的栅格是否存在障碍物;如果蚂蚁想查看2级视野范围的栅格是否存在障碍物,前提条件是在蚂蚁当前位置与2级视野范围内的栅格连线之间不存在障碍物。例如:如图12(a)所示,蚂蚁在节点i处若想查看2级视野区域的节点k处是否存在障碍物,必须保证1级视野区域的A栅格和B栅格处不存在障碍物。因此,若A栅格和B栅格任意一处存在障碍物均导致蚂蚁无法从节点i直接查看节点k处的情况,如图11(a)、11(b)、11(c)所示,通过上述分析可知,如果某个1级视野范围内的栅格为障碍物时会导致2级视野范围内的部分栅格不可见。图13给出了所有2级视野范围内的栅格不可见的情况,黑色代表障碍物栅格,红色代表不可见栅格。在发明提出的2种启发信息计算方法中,无论使用第一种方法还是第二种方法计算领域节点启发信息值时,会发现栅格中心位于同一连线上的邻居栅格(如图14所示)的启发信息值始终是相同的。如果在蚂蚁移动的方向上存在两栅格启发信息值相同,规定蚂蚁选择1级视野范围内的邻居栅格作为下一时刻到达的位置;
S4、根据所述启发信息值计算栅格的转移概率,并根据所述转移概率的控制参数确定蚂蚁下一时刻到达的位置,更新禁忌表,其实现方法如下:
其中,τik表示t时刻路段(i,k)上的信息素量,ηik表示蚂蚁m从节点i到节点k的启发信息值,τis表示蚂蚁m从节点i到节点s的信息素浓度,α表示信息启发式因子,β为期望启发式因子,allowedm表示蚂蚁m下一步允许选择的节点集合;
S402、生成一个随机数ω,且ω∈[0,1];
S403、判断所述转移概率的控制参数σ是否大于所述随机数ω,若是,则进入步骤S304,反之,则进入步骤S305,其中,σ∈[0,1];
S404、通过轮盘赌策略选择邻域栅格作为蚂蚁下一时刻到达位置,更新禁忌表;
S405、选择转移概率最大的邻域栅格作为蚂蚁下一时刻到达位置,更新禁忌表。
在具体实施例中,栅格的转移概率大小是蚂蚁决定是否选取某邻居栅格作为下一时刻到达位置的最终判断依据。本发明通过设置转移概率的控制参数δ,δ∈[0,1]对蚂蚁的转移策略进行调控。具体操作:在蚂蚁每次选择邻居栅格作为下一步到达位置之前,生成一个随机数ω,ω∈[0,1],如果ω<σ,则通过轮盘赌策略选择邻居栅格作为蚂蚁下一时刻到达位置;否则,选择转移概率最大的邻居栅格作为蚂蚁下一时刻到达位置。
通过分析可知,当转移概率的控制参数δ=0时,算法类似于贪心算法;当转移概率的控制参数δ=1时,算法采用的是轮盘赌策略;当转移概率的控制参数δ∈[0,1]时,算法兼具上述两种情况的特点;
S5、将所有蚂蚁完成一次路径搜索,并更新成功抵达终点的蚂蚁走过的路径上的栅格信息素,所述栅格信息素τik(t+1)的表达式如下:
τik(t+1)=(1-ρ)τik(t)+Δτik(t)
S6、判断是否完成预设的迭代次数,若是,则输出最优路径,从而完成移动机器人的路径规划,反之,则返回步骤S2。
为了进一步对本发明进行说明,采用规模为30*30的栅格地图,在栅格地图中随机生成障碍物,用ζ表示地图中障碍物的分布比例,δ表示蚂蚁的转移概率的控制参数、N表示蚁群算法的代次数、M表示蚂蚁投放数量、α表示信息启发式因子、β表示期望启发式因子、ρ表示信息素挥发系数,上述相关参数设置见表2。
表2
如图15(a)、(b)、(c)为使用第一种启发信息计算方式得到的蚂蚁搜索路径轨迹;图15(d)、(e)、(f)为使用第二种启发信息计算方式得到的蚂蚁搜索路径轨迹。对比蚂蚁的搜索路径分布特点可以发现,第一种启发信息计算方式搜索到的路径范围相对于第二种启发信息计算方式更加广泛。因此,说明了第一种启发信息计算方式相对于第一种启发信息计算方式具有搜索全局性,方法2相对于方法1具有搜索方向性。
采用规模为30*30的栅格地图,ζ=10%,ζ表示地图中障碍物的分布比例,转移概率的控制参数δ分别取值{0.2,0.4,0.6,0.8,1},相关参数见表3。此外,采用第一种启发信息计算方式:
表3
如图16所示,可以发现随着转移概率的控制参数δ取值的增大,蚂蚁的搜索范围也随着增大,因此,可以通过设置转移概率的控制参数δ的值,实现对算法搜索范围的调控。为了探寻转移概率的控制参数δ取何值时,算法的路径寻优效果最好,对转移概率的控制参数δ∈{0.2,0.4,0.6,0.8,1}的每一个取值进行10次仿真,记录每次仿真的最短路径值。图17仿真结果图,图18为转移概率的控制参数δ取不同值时10次仿真结果的平均值,从仿真结果中可以发现,当转移概率的控制参数δ=0.8时算法路径寻优效果最好。
为了进一步对本发明进行说明,分别在规模为15*15、30*30、60*60,ξ=20%的栅格地图中对本发明算法和标准蚁群算法进行分析,ζ表示地图中障碍物的分布比例,算法仿真中的相关参数见表2,采用第一种启发信息计算方式。算法搜索到的路径长度迭代结果如图19所示;本发明算法和标准蚁群算法搜索到的最短路径如图20所示,以及对应的在算法搜索过程中的蚂蚁爬行路径轨迹如图21所示。两种算法搜索到的最短路径见表4,可以发现本发明算法无论是在小、中、还是大规模场景地图中,本发明算法的路径搜索效果均比标准蚁群算法更佳。在中、小规模地图中传统蚁群算法可以搜索最优或者近似最优值;在大规模场景下,标准蚁群算法搜索具有强烈的随机性(如图21(c)所示),虽然可以找到从起点到终点的可行路径,但与最优值相差甚远。
表4
通过观察数据可以发现本发明算法搜索到的路径长度也比标准蚁群算法搜索到的最优解的路径长度短。本发明通过改进的蚁群算法能够搜索到更优的全局最短路径,可以提高蚁群算法搜索的有效性,降低搜索的随机性,为机器人的路径规划提供良好的条件。
Claims (7)
1.一种基于改进蚁群算法的移动机器人路径规划方法,其特征在于,包括如下步骤:
S1、根据环境障碍物分布特点建立栅格地图,初始化栅格信息素,并设置初始参数和确定移动机器人的起点和终点;
S2、初始化禁忌表,将蚁群算法的蚂蚁参数设为M只蚂蚁,并将移动机器人的起点和终点设为M只蚂蚁的起点和终点;
S3、根据蚂蚁的移动规则计算栅格地图中领域节点的启发信息值;
所述步骤S3中的蚂蚁的移动规则包括:蚂蚁的搜索方向以及蚂蚁的视野范围;
所述步骤S3中的启发信息值ηik的计算方法包括如下两种方法:
第一种计算方法:
其中,θik表示蚂蚁搜索方向与理论最佳搜索方向之间的夹角,Dik表示蚂蚁下一时刻的移动方向,DiE表示蚂蚁理论上的最佳搜索方向,xk表示蚂蚁下一时刻到达位置的横坐标,xi表示蚂蚁当前位置的横坐标,xE表示终点的横坐标,yk表示蚂蚁下一时刻到达位置的纵坐标,yi蚂蚁当前位置的纵坐标,yE表示终点的纵坐标;
第二种计算方法:
其中,θik表示蚂蚁搜索方向与理论最佳搜索方向之间的夹角,xk表示蚂蚁下一时刻到达位置的横坐标,xi表示蚂蚁当前位置的横坐标,xE表示终点的横坐标,yk表示蚂蚁下一时刻到达位置的纵坐标,yi蚂蚁当前位置的纵坐标,yE表示终点的纵坐标;
S4、根据所述启发信息值计算栅格的转移概率,并根据所述转移概率的控制参数确定蚂蚁下一时刻到达的位置,更新禁忌表;
其中,τik表示t时刻路段(i,k)上的信息素量,ηik表示蚂蚁m从节点i到节点k的启发信息值,τis表示蚂蚁m从节点i到节点s的信息素浓度,α表示信息启发式因子,β为期望启发式因子,allowedm表示蚂蚁m下一步允许选择的节点集合;
S5、将所有蚂蚁完成一次路径搜索,并更新成功抵达终点的蚂蚁走过的路径上的栅格信息素;
S6、判断是否完成预设的迭代次数,若是,则输出最优路径,从而完成移动机器人的路径规划,反之,则返回步骤S2。
2.根据权利要求1所述的基于改进蚁群算法的移动机器人路径规划方法,其特征在于,所述步骤S1中设置初始参数包括设置栅格的初始信息素浓度、蚁群算法的迭代次数、蚂蚁投放数量、信息启发式因子、期望启发因子以及栅格信息素的挥发系数。
3.根据权利要求1所述的基于改进蚁群算法的移动机器人路径规划方法,其特征在于,所述蚂蚁的搜索方向为:搜索栅格地图中的16个方向和24个领域节点,且蚂蚁的转角为22.5°。
5.根据权利要求4所述的基于改进蚁群算法的移动机器人路径规划方法,其特征在于,针对1级视野范围栅格和2级视野范围栅格的启发信息值相同时,蚂蚁选择1级视野范围内的领域栅格作为下一时刻到达的位置。
6.根据权利要求1所述的基于改进蚁群算法的移动机器人路径规划方法,其特征在于,所述步骤S4包括如下步骤:
S401、根据所述启发信息值计算栅格的转移概率;
S402、生成一个随机数ω,且ω∈[0,1];
S403、判断所述转移概率的控制参数σ是否大于所述随机数ω,若是,则进入步骤S404,反之,则进入步骤S405,其中,σ∈[0,1];
S404、通过轮盘赌策略选择邻域栅格作为蚂蚁下一时刻到达位置,更新禁忌表;
S405、选择转移概率最大的邻域栅格作为蚂蚁下一时刻到达位置,更新禁忌表。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910535562.9A CN110244733B (zh) | 2019-06-20 | 2019-06-20 | 一种基于改进蚁群算法的移动机器人路径规划方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910535562.9A CN110244733B (zh) | 2019-06-20 | 2019-06-20 | 一种基于改进蚁群算法的移动机器人路径规划方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110244733A CN110244733A (zh) | 2019-09-17 |
CN110244733B true CN110244733B (zh) | 2020-06-30 |
Family
ID=67888369
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910535562.9A Active CN110244733B (zh) | 2019-06-20 | 2019-06-20 | 一种基于改进蚁群算法的移动机器人路径规划方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110244733B (zh) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110647162B (zh) * | 2019-10-16 | 2022-10-14 | 厦门理工学院 | 一种导游无人机路径规划方法、终端设备及存储介质 |
CN110941267B (zh) * | 2019-11-18 | 2021-05-25 | 电子科技大学 | 蚁群算法用于拥堵环境下的动态路径规划方法 |
CN111152214B (zh) * | 2019-12-13 | 2021-08-03 | 珠海格力电器股份有限公司 | 一种四自由度码垛机器人、控制系统、码垛路径规划方法 |
CN111159523B (zh) * | 2019-12-23 | 2024-03-01 | 河北工业大学 | 一种基于Spark的并行蚁群优化社区发现方法 |
CN111738396B (zh) * | 2020-06-01 | 2023-09-26 | 北京中安智能信息科技有限公司 | 一种应用于潜艇路径规划的自适应栅格颗粒度蚁群方法 |
CN111857141B (zh) * | 2020-07-13 | 2021-07-27 | 武汉理工大学 | 一种机器人路径规划方法、装置、设备及存储介质 |
CN113159391A (zh) * | 2021-03-27 | 2021-07-23 | 桂林理工大学 | 一种求解带交通选择路程规划问题的多目标存档蚁群寻优方法 |
CN114578772B (zh) * | 2021-04-16 | 2023-08-11 | 青岛中车四方车辆物流有限公司 | 一种agv集群控制系统设计框架及方法 |
CN113110465A (zh) * | 2021-04-22 | 2021-07-13 | 哈尔滨理工大学 | 一种基于改进蚁群算法的模块路径规划方法 |
CN113110472A (zh) * | 2021-04-25 | 2021-07-13 | 深圳市跨越新科技有限公司 | 路径规划方法、装置和终端 |
CN113703391B (zh) * | 2021-08-11 | 2023-01-31 | 广东省科学院智能制造研究所 | 一种基于随机信息素优化精英蚁群算法的加工路径规划方法 |
CN113703450B (zh) * | 2021-08-23 | 2024-03-29 | 皖西学院 | 基于平滑因素改进蚁群算法的移动机器人路径规划方法 |
CN114415668B (zh) * | 2021-12-23 | 2023-07-07 | 西北工业大学 | 基于扩展视野自适应蚁群算法的移动机器人路径规划方法 |
CN116795108B (zh) * | 2023-06-09 | 2023-12-01 | 西南交通大学 | 基于多源感知信号的智能无人车配送方法 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU2009302537A1 (en) * | 2008-10-06 | 2010-04-15 | Somalogic, Inc. | Ovarian cancer biomarkers and uses thereof |
CN102288192B (zh) * | 2011-07-01 | 2013-11-27 | 重庆邮电大学 | 基于Ad Hoc网络的多机器人路径规划方法 |
CN103823466B (zh) * | 2013-05-23 | 2016-08-10 | 电子科技大学 | 一种动态环境下移动机器人路径规划方法 |
CN105760954A (zh) * | 2016-02-15 | 2016-07-13 | 南通大学 | 一种基于改进蚁群算法的泊车系统路径规划方法 |
CN106225788B (zh) * | 2016-08-16 | 2019-04-19 | 上海理工大学 | 基于路径拓展蚁群算法的机器人路径规划方法 |
CN106200650A (zh) * | 2016-09-22 | 2016-12-07 | 江苏理工学院 | 基于改进蚁群算法的移动机器人路径规划方法及系统 |
CN107272679B (zh) * | 2017-06-15 | 2020-06-16 | 东南大学 | 基于改进的蚁群算法的路径规划方法 |
CN108180914B (zh) * | 2018-01-09 | 2021-06-18 | 昆明理工大学 | 一种基于蚁群改进和尖峰平滑的移动机器人路径规划方法 |
CN109282815A (zh) * | 2018-09-13 | 2019-01-29 | 天津西青区瑞博生物科技有限公司 | 一种动态环境下基于蚁群算法的移动机器人路径规划方法 |
-
2019
- 2019-06-20 CN CN201910535562.9A patent/CN110244733B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN110244733A (zh) | 2019-09-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110244733B (zh) | 一种基于改进蚁群算法的移动机器人路径规划方法 | |
CN111024092B (zh) | 一种多约束条件下智能飞行器航迹快速规划方法 | |
CN109116841B (zh) | 一种基于蚁群算法的路径规划平滑优化方法 | |
CN113467456B (zh) | 一种未知环境下用于特定目标搜索的路径规划方法 | |
CN112985408B (zh) | 一种路径规划优化方法及系统 | |
CN111323016A (zh) | 一种基于自适应蚁群算法的移动机器人路径规划方法 | |
CN107229287A (zh) | 一种基于遗传蚂蚁算法的无人机全局路径规划方法 | |
CN106525047A (zh) | 一种基于floyd算法的无人机路径规划方法 | |
Mi et al. | Restoring connectivity of mobile robotic sensor networks while avoiding obstacles | |
CN110989352A (zh) | 一种基于蒙特卡洛树搜索算法的群体机器人协同搜索方法 | |
CN113552899B (zh) | 一种通信链路受限下的多无人机自主协同搜索方法 | |
CN110381442B (zh) | 一种基于隐式信息交互方式的群机器人目标搜索方法 | |
CN109211242B (zh) | 一种融合rrt与蚁群算法的三维空间多目标路径规划方法 | |
CN115202394A (zh) | 一种基于改进遗传算法的无人机全覆盖路径规划方法 | |
CN114167865A (zh) | 一种基于对抗生成网络与蚁群算法的机器人路径规划方法 | |
CN114675643A (zh) | 一种无线传感网络的信息传输路径规划方法及装置 | |
CN115994891B (zh) | 基于狼群算法的无人载具混凝土坝表面缺陷动态检测方法 | |
Wang et al. | UAV online path planning based on improved genetic algorithm with optimized search region | |
CN114578845B (zh) | 一种基于改进蚁群算法的无人机航迹规划方法 | |
CN115129064A (zh) | 基于改进萤火虫算法与动态窗口法融合的路径规划方法 | |
CN111896001A (zh) | 一种三维蚁群航迹优化方法 | |
CN114237282A (zh) | 面向智慧化工业园区监测的无人机飞行路径智能规划方法 | |
CN113932819B (zh) | 基于往复回填算法的无人机完整搜索路径规划方法 | |
You et al. | Multi-AUV Static Target Search Based on Improved PSO | |
CN117232512B (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 |