一、简介
路径规划是实现移动机器人自主导航的关键技术，是指在有障碍物的环境中，按照一定的评价标准（如距离、时间、能耗等），寻找到一条从起始点到目标点的无碰撞路径，这里选取最短距离路径规划的评价标准，即最短路径规划问题。
1．路径规划数学模型的建立
将移动机器人周围环境用一组数据进行抽象表达，建立二维或三维的环境模型，得到移动机器人能够理解分析的环境数据，是机器人路径规划的基本前提。我这里用的是栅格法，其原理是将周围环境看成一个二维平面，将平面分成一个个等面积大小的具有二值信息的栅格，每个栅格中存储着周围环境信息量，下图我给出了一个栅格法地图，方便大家更好的理解栅格地图。这里设计的栅格地图为一个20×20的地形矩阵，黑色的地方表示有障碍，白色的地方表示没有障碍。
图1栅格法地图
在用栅格法建立环境模型时，为了将环境信息转换成移动机器人可以识别的数据，一般采用序号法标记环境地图信息，即将栅格地图中一个个栅格从序号１依次累加直到标记到最后一个栅格。如图2所示。
图3八叉树搜索策略
那么，怎么判断一个栅格点是否为另一个栅格点的相邻栅格点呢，另外，又怎么判断是否为有障碍栅格呢。这就需建立矩阵D，记录每个栅格点至其相邻栅格点的代价值。本例中栅格地图有20×20个栅格点，则D的大小为400×400，其中列是起点栅格，行是局部终点栅格，各栅格点至其各相邻无障碍栅格点的代价值非零，而有障碍栅格及非相邻栅格设为0。
2．机器人最短路径规划的实现步骤
蚁周模型实现机器人最短路径规划的流程图
为了方便大家更好地理解蚁群算法的原理及实现过程，其流程图如图4所示。（流程图较长，我截图了两段。）
图4基于蚁群算法的机器人最小路径规划流程图
图中公式（3）（4）的具体表达在下边的具体步骤里。
蚁周模型实现机器人最短路径规划的具体步骤
**步骤1：**给出栅格地图的地形矩阵；初始化信息素矩阵 Tau（记录每个栅格至其他栅格的信息素量），最大迭代次数K，蚂蚁个数M，表征信息素重要程度的参数 、表征启发式信息重要程度的参数 ，信息素蒸发系数 ，信息素增加强度系数Q及启发式信息矩阵
**步骤2：**构建启发式信息矩阵。按式（1）和式（2）计算每个栅格至目标点的距离，启发式信息素取为至目标点距离的倒数，距离越短，启发式因子越大，障碍物处的启发式信息为0。建立矩阵D，用以存储每个栅格点至各自相邻无障碍栅格点的代价值。
**步骤3：**对于每一只蚂蚁，初始化蚂蚁爬行的路径及路径长度，将禁忌列表全部初始化为1；蚂蚁从起始点出发开始搜索路径，找出当前栅格点的所有无障碍相邻栅格点（即矩阵D中相应元素不为0的栅格点），再根据禁忌列表筛选出当前可选择的栅格点。
**步骤4：**如果起始点是目标点，且可选栅格点个数大于等于1，则根据式（3）计算蚂蚁从当前栅格点转移到各相邻栅格点的概率，
二、部分源代码
clear;clcclose all;ticG=[0 0 0 0 0 0 1 1 1 0 0 0 0 0 0 0 0 0 0 0; 0 0 0 0 0 0 1 1 1 0 0 0 0 0 0 0 0 0 0 0; 0 0 0 0 0 0 1 1 1 0 0 0 0 0 0 0 0 0 0 0; 0 0 0 0 0 0 1 1 1 0 0 0 0 0 0 0 0 0 0 0; 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0; 0 1 1 1 0 0 1 1 1 0 0 0 0 0 0 0 0 0 0 0; 0 1 1 1 0 0 1 1 1 0 0 0 0 0 0 0 0 0 0 0;0 1 1 1 0 0 1 1 1 0 1 1 1 1 0 0 0 0 0 0; 0 1 1 1 0 0 0 0 0 0 1 1 1 1 0 0 0 0 0 0; 0 0 0 0 0 0 0 0 0 0 1 1 1 1 0 0 0 0 0 0; 0 0 0 0 0 0 0 1 1 0 1 1 1 1 0 0 0 0 0 0; 0 0 0 0 0 0 0 1 1 0 1 1 1 1 0 0 0 0 0 0; 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 1 1 1 1 0; 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 1 1 1 1 0; 1 1 1 1 0 0 0 0 0 0 0 0 1 1 0 1 1 1 1 0; 1 1 1 1 0 0 1 1 0 0 0 1 0 0 0 0 0 0 0 0; 0 0 0 0 0 0 1 1 0 1 1 1 0 0 0 0 0 1 1 0; 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0; 0 0 1 1 0 0 0 0 0 0 1 1 0 0 1 0 0 0 0 0; 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 0;];%地图矩阵n=size(G,1);%n表示地图大小m=50; %% m 蚂蚁个数Alpha=2; %% Alpha 表征信息素重要程度的参数Beta=6; %% Beta 表征启发式因子重要程度的参数Rho=0.1; %% Rho 信息素蒸发系数NC_max=100; %%最大迭代次数Q=1; %%信息素增加强度系数Tau=ones(n,n);%Tau为信息素矩阵NC=1;%迭代计数器，记录迭代次数r_e=1; c_e=20;%地图终点在矩阵中的位置%可以通过position2rc函数产生s=n;%路径起始点在矩阵中的位置position_e=n*(n-1)+1;%路径终点在矩阵中的位置min_PL_NC_ant=inf;%%蚂蚁最短的行进距离min_ant=0;%%最短行进距离的蚂蚁坐标min_NC=0;%%最短行进距离的迭代次数% 计算邻接矩阵及启发因子%%邻接矩阵作用是计算启发因子z=1;max_generation=200;p_crossover=0.2;p_mutation=0.05;weight_length=20;weight_smooth=1;new_population={};for i=1:nfor j=1:nif G(i,j)==0 D(i,j)=((i-r_e)^2+(j-c_e)^2)^0.5;elseD(i,j)=inf;%i=j时不计算，应该为0，但后面的启发因子要取倒数，用eps（浮点相对精度）表示endendend D(r_e,c_e)=0.05;Eta=1./D;%Eta为启发因子，这里设为到终点距离的倒数Tau=10.*Eta;%%%%%创新点%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%计算移动矩阵D_move=zeros(n*n,8);%%D_move每一行代表与行标对应元素，可以前往的下一个节点的位置for point=1:n*nif G(point)==0[r,c]=position2rc(point);move=1;for k=1:nfor m1=1:nim=abs(r-k);jn=abs(c-m1); if im+jn==1||(im==1&&jn==1) if G(k,m1)==0D_move(point,move)=(m1-1)*n+k;move=move+1;endendendendend end%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%移动矩阵和邻接矩阵计算完成，检查无误%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%开始迭代routes=cell(NC_max,m);%%%%存储每次迭代每个蚂蚁的路径PL=zeros(NC_max,m); %%%%%存储每次迭代每个蚂蚁的路径长度while NC<=NC_maxfor ant=1:mcurrent_position=s;%%%当前位置为起始点path=s;%%路径初始化PL_NC_ant=0;%%长度初始化Tabu=ones(1,n*n); %%%%禁忌表，排除已经走过的位置Tabu(s)=0;%%排除已经走过的初始点D_D=D_move;%%%%D_D是D_move的中间矩阵，作用是为了不让D_move参与计算，也可不用D_D矩阵，直接用D_moveD_work=D_D(current_position,:);%%%把当前点可以前往的写一个节点的信息传送给D_worknonzeros_D_work=find(D_work);%%%找到不为0的元素的位置for i1=1:length(nonzeros_D_work)if Tabu(D_work(i1))==0D_work(nonzeros_D_work(i1))=[];%%将禁忌表中已走过的元素删除，防止走已经走过的位置D_work=[D_work,zeros(1,8-length(D_work))];%%%保证D_work向量长度为8（每个点最多能往周围的8个点走），为后面for循环做准备endend%%%%%%%%%%%%%%%%%%%%%%%%%%%%排除走过的第一点（排除起点）%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%len_D_work=length(find(D_work));while current_position~=position_e&&len_D_work>=1%%当前点是否为终点或者走进死胡同p=zeros(1,len_D_work);for j1=1:len_D_work[r1,c1]=position2rc(D_work(j1));%%利用自己编的函数把可以前进的点计算为行列表示p(j1)=(Tau(r1,c1)^Alpha)*(Eta(r1,c1)^Beta);%%%%计算每个可以前往的节点的概率endp=p/sum(p);%%%归一化pcum=cumsum(p);%%%概率累加select=find(pcum>=rand);%%%%轮盘赌法选择下个节点to_visit=D_work(select(1));%%%前往下一个节点%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%到达下一个节点%%%%%%%%%%%%%%%%%%%%%%%%path=[path,to_visit];%%%路径累加dis=distance(current_position,to_visit);%%%计算到下个节点的距离PL_NC_ant=PL_NC_ant+dis;%%距离累加current_position=to_visit;%%%当前点设为前往点D_work=D_D(current_position,:);%%%%把当前节点可以前往的下一个节点的信息传给D_workTabu(current_position)=0;%%%禁忌表中排除已经到的点for kk=1:400if Tabu(kk)==0for i3=1:8if D_work(i3)==kkD_work(i3)=[];%%%%排除禁忌表中已经走过的节点D_work=[D_work,zeros(1,8-length(D_work))];%%保证长度为8endendendendlen_D_work=length(find(D_work));%%%计算当前点可以前往的下一个节点的数量end%%%%%%%%%%%%%%%%%%%%%%%%迭代一次所有蚂蚁走完%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%routes{NC,ant}=path;%%%把蚂蚁走过的路径记录下来if path(end)==position_ePL(NC,ant)=PL_NC_ant;%%记录到达终点的蚂蚁的行进距离if PL_NC_ant<min_PL_NC_antmin_NC=NC;min_ant=ant;min_PL_NC_ant=PL_NC_ant;%%记录路径最短的蚂蚁的迭代次数和属于那一只endelsePL(NC,ant)=0;endend
三、运行结果
四、matlab版本及参考文献
1 matlab版本
a
2 参考文献
[1] 包子阳,余继周,杨杉.智能优化算法及其MATLAB实例（第2版）[M].电子工业出版社，.
[2]张岩,吴水根.MATLAB优化算法源代码[M].清华大学出版社，.
【路径规划】基于matlab蚁群优化遗传算法求解机器人栅格地图最短路径规划问题【含Matlab源码 1581期】
