一、机器人路径规划介绍
移动机器人(Mobile robot,MR)的路径规划是 移动机器人研究的重要分支之,是对其进行控制的基础。根据环境信息的已知程度不同,路径规划分为基于环境信息已知的全局路径规划和基于环境信息未知或局部已知的局部路径规划。随着科技的快速发展以及机器人的大量应用,人们对机器人的要求也越来越高,尤其表现在对机器人的智能化方面的要求,而机器人自主路径规划是实现机器人智能化的重要步骤,路径规划是指规划机器人从起点位置出发,无碰撞、安全到达指定目标位置的最优路径。目前,常用的移动机器人全局路径规划方法很多,如栅格法和人工势场法。对于栅格法,当空间增大时,所需存储空间剧增,决策速度下降;而人工势场法容易产生局部最优解问题和死锁现象。随着智能控制技术的发展,出现了如遗传算法算法、粒子群优化算法、麻雀搜索算法、灰狼优化算法、鲸鱼优化算法等。
参考文献:
[1]史恩秀,陈敏敏,李俊,等.基于蚁群算法的移动机器人全局路径规划方法研究[J].农业机械学报, 2014, 45(6):5.DOI:CNKI:SUN:NYJX.0.2014-06-009.
[2]朱庆保,张玉兰.基于栅格法的机器人路径规划蚁群算法[J].机器人, 2005, 27(2):5.DOI:10.3321/j.issn:1002-0446.2005.02.008.
[3]曹新亮,王智文,冯晶,等.基于改进蚁群算法的机器人全局路径规划研究[J].计算机工程与科学, 2020, 42(3):7.DOI:CNKI:SUN:JSJK.0.2020-03-027.
二、栅格地图环境搭建
首先建立移动机器人工作环境,设移动机器人的工作空间为二维空间(记为RS),工作环境中的障碍物即为机床。在机器人运动过程中,障碍物为静止且大小不发生变化。按栅格法划分RS,移动机器人在栅格间行走。无障碍物的栅格为可行栅格,有障碍物的栅格为不可行栅格。栅格集包含所有栅格。栅格标识有:直角坐标法和序号法。本文采用序号标识法。
在移动机器人工作空间下按从左到右,从上到下的顺序,依次标记为序号1,2,3,···,n,每一个序号代表一个栅格。为了避免移动机器人与障碍物发生碰撞,可以将障碍物膨胀,障碍物在占原有栅格的同时,再占多个栅格,按 个栅格算。这种划分方法简单实用,能够满足环境模型与真实情况相符。从而使移动机器人在路径规划时畅通无阻。令S={1,2,3,···,N}为栅格序号集。根据上述对应关系,可知g(0,0)的序号为1,g(1,0)序号为2,直至g(X,Y)的序号为n。规划起始位置、目标位置均为任意且都属于S(但不在同一栅格内)。
在实际工作环境中,移动机器人工作环境是复杂多变的,且为三维空间。为了便于研究,本文对环境进行简化建模。栅格法是一种常用的环境表示方法,因其简单方便(二维环境),环境建模的复杂性小,因而本文环境建模采用栅格法。在栅格地图中,工作环境被划分为很多栅格,其中包括有障碍物和无障碍的栅格,在仿真程序中用0表示此栅格无障碍物,机器人可以通过此栅格,用1表示栅格有障碍物,机器人无法通过,需选择其他栅格。栅格的尺寸大小可根据工作环境中的障碍物尺寸以及安全距离进行设置。为了实现程序仿真,需要对栅格进行标识,如下图所示,以20x20的栅格环境为例来说明。
如上图所示,白色栅格表示无障碍物的栅格,黑色栅格则表示有障碍物的栅格,在地图中对每个栅格编号,不同序号的栅格在坐标系中的坐标可用下式来表示:
x=mod(Ni/N)-0.5
y=N-ceil(Ni/N)+0.5
其中,mod为取余运算,ceil表示向后取整,Ni是对应栅格的标号,N表示每 列的栅格数量,取栅格中心位置作为栅格在坐标系中的坐标。这样机器人全局路径规划的问题就转变成了利用算法在栅格地图上寻找由起始点到目标点的有序的栅格子集,这些栅格子集的中心连线便是算法寻找的路径。
参考文献:
[1]史恩秀,陈敏敏,李俊,等.基于蚁群算法的移动机器人全局路径规划方法研究[J].农业机械学报, 2014, 45(6):5.DOI:CNKI:SUN:NYJX.0.2014-06-009.
[2]曹新亮,王智文,冯晶,等.基于改进蚁群算法的机器人全局路径规划研究[J].计算机工程与科学, 2020, 42(3):7.DOI:CNKI:SUN:JSJK.0.2020-03-027.
三、机器人路径规划多目标模型
3.1路径成本
当机器人从起点向目标点移动时,通常选择最短的一条路径。将一条路径上每 条线段的长度累加求和即可得到路径总长度。任何两点形成的线段都是根据欧氏距 离来计算的,其中,P = [P0, P1, …, Pi, Pi + 1, … Pn, Pn + 1]代表路径 P,S = P0代表起 点,T = Pn + 1代表目标点。路径长度目标计算方法如下式
其中,Pi =(xi, yi)和 Pi + 1 = (xi + 1, yi + 1)是路径中的两个连续点,d(Pi, Pi + 1)是 路径中的线段距离;Length(P)表示由所有线段相加得到的总路径长度;n 表示路 径中的点数量。
3.2平滑成本
路径平滑度表示路径的弯曲程度,换句话说,只有当路径是平滑的,机器人 在移动时才会使用较少的能量。为了衡量可行路径的平滑度,使用两个连续路段之 间的夹角 Ang [Pi, Pi + 1, Pi + 2]表示。路径平滑度的计算方法如下式
其中 Pi,Pi + 1和 Pi + 2是路径上的三个相邻节点。
3.3目标函数
移动机器人(Mobile robot,MR)的路径规划的目标函数f1和f2分别是路径成本最小和平滑成本最小。
参考文献:
[1]于振翱. 面向多目标优化的移动机器人路径规划方法研究[D]. 山东:聊城大学,2023.
[2]杨嘉. 基于改进NSGA-Ⅱ算法的移动机器人路径规划研究[D]. 陕西:长安大学,2021.
四、MMO-DE-CSCD求解移动机器人路径规划
基于聚类技术的差分进化算法(MMODE-CSCD)在解决多模式多目标优化问题(MMOPs)上效果显著。以下是该算法的介绍和描述:
问题定义:MMOPs 是指在目标空间中具有相同或相似Pareto最优前沿,而在决策空间中存在多个Pareto最优解集的问题。解决这类问题需要优化算法能够定位多个Pareto集(PSs)。
算法设计:MMO-DE-CSCD算法设计了一种基于聚类的特殊拥挤距离(CSCD)方法,用于在决策和目标空间中计算解的综合拥挤度。此外,引入了基于距离的精英选择机制(DBESM),以确定各个个体的学习样本,并在样本周围生成新的个体,以获得在决策和目标空间中分布良好的种群。
CSCD方法:CSCD方法利用k-means算法将同一非支配层的个体划分为多个类别,使得同一类别中的个体来自相同的Pareto集,从而更准确地计算拥挤距离,并在精英选择和环境选择过程中将拥挤距离作为指标。
DBESM机制:DBESM机制旨在改善所有Pareto集上的种群分布,综合考虑种群的多样性和收敛性,为每个个体确定学习样本。
算法框架:MMODE-CSCD算法的框架包括差分进化的基本操作,如变异、交叉和选择,以及将CSCD和DBESM嵌入到非支配排序方案中,以获得在决策空间和目标空间中分布良好的解。
参考文献:
[1]Liang, J.J., Qiao, K., Yue, C., Yu, K., Qu, B., Xu, R., Li, Z., & Hu, Y. (2021). A clustering-based differential evolution algorithm for solving multimodal multi-objective optimization problems. Swarm Evol. Comput., 60, 100788.
4.1部分代码
close all
clear
clc
dbstop if all error
addpath("./ /")
global G S E
MultiObj= fun_info();%获取无人机模型信息
params.maxgen=100; % 最大迭代次数
params.Np=50; % 种群大小
params.Nr=100; %外部存档大小(不得小于种群大小)
[Xbest,Fbest] = (params,MultiObj);
%% 获取算法得到的所有路径存于Result中
for i=1:size(Xbest,1)
global_best = round(Xbest(i,:));
route = [S(1) global_best E(1)];
path=generateContinuousRoute(route,G);
% path=shortenRoute(path);
path=GenerateSmoothPath(path,G);
Result(i).BestPosition= path;
Result(i).BestFit=Fbest(i,:);
end
%% 获取 路径成本最小 与 平滑成本最小
if size(Fbest,1)>1
idx=min(Fbest);
else
idx=Fbest;
end
KK=find(Fbest(:,1)==idx(1));
IDX(1)=KK(1);
KK=find(Fbest(:,2)==idx(2));
IDX(2)=KK(1);
https://mbd.pub/o/bread/mbd-ZpqYmplw
https://mbd.pub/o/bread/mbd-ZpqYmpdx
https://mbd.pub/o/bread/mbd-ZpqYmZ5s
https://mbd.pub/o/bread/mbd-ZpqYmZxr
https://mbd.pub/o/bread/mbd-ZpqYmZpu
%% 画图pareto前沿图
dbclear all
figure
plot(Fbest(:,1),Fbest(:,2),'r*');
xlabel('路径成本')
ylabel('平滑成本')
legend('')
4.2部分结果
五、完整MATLAB代码
见下方联系方式