【路径规划】基于matlab人工蜂群优化粒子群算法避障最短路径规划【含Matlab源码 124期】

举报
海神之光 发表于 2022/05/29 00:03:44 2022/05/29
【摘要】 一、人工蜂群算法简介 人工蜂群算法从提出开始, 由于其参数少, 自适应强, 受到很多研究者的追捧, 广泛应用于各大领域, 其中就包括移动机器人的路径规划领域。 人工蜂群算法的数学模型 初始化会形成一个...

一、人工蜂群算法简介

人工蜂群算法从提出开始, 由于其参数少, 自适应强, 受到很多研究者的追捧, 广泛应用于各大领域, 其中就包括移动机器人的路径规划领域。

人工蜂群算法的数学模型
初始化会形成一个初始食物源。即初始解 (x1, x2, …, xn, …, xN) :
在这里插入图片描述
式中:lbi为可行解的最小值;ubi是可行解的最大值;rand (0, 1) 生成 (0, 1) 之间的随机数;n∈{1, 2, 3, …Fn}, Fn为食物源数目。
食物源生成之后, 评价函数会对生成的食物源进行评价, 评价函数如式 (2) 所示:
在这里插入图片描述
式中:f (xn) 为对应的解向量xn的目标函数值。

雇佣蜂在当前食物源附近搜索新食物源的方法如式 (3) 所示:
在这里插入图片描述
式中:xni′为新食物源;xni为当前食物源;ε为搜索系数, 其值介于[-1, 1]之间;xki是食物源中一个不同于xni的食物源, k∈{1, 2, 3, …Fn}。

当搜索到新食物源, 雇佣蜂会对这个新食物源进行评价, 并与之前的食物源进行对比, 来决定放弃哪个食物源。

观察蜂参照雇佣蜂返回时所带来的食物源信息进行选择觅食。它们会利用食物源的评价函数来计算他们选择某个食物源的可开采的概率probabilityn, 计算方法如式 (4) 所示:
在这里插入图片描述
观察蜂会因个体差异, 存在各自对食物源的期望。然后将期望和可开采概率比较, 来决定观察蜂是否选择此食物源进行开采。

侦察蜂的作用是搜索新的食物源, 并不参照当前食物源, 而是类似初始化一样随机生成。在此之前会对选择次数超过限定值的食物源进行放弃, 如式 (5) 所示:
在这里插入图片描述
式中:trial是当前蜜蜂对食物源优化的次数。

二、部分源代码


clc;
clear;
close all;
%% ABC + PSO Path Planning Problem
% 这里,系统试图找到起点之间的最优路径
%借助人工蜂群 (ABC) 算法的百分比和目的地点
% 和粒子群优化算法相结合。 这是个好策略
% 在机器人寻路中。 在每次运行中新位置的新障碍
% 定义和一条曲线试图找到最佳路径。 运行多次
% 以找到最佳结果。

%% 清除
clc;
clear;
warning('off');

%% 启动 ABC + PSO 最优路径查找器
model=Basics();
model.n=6;  % 处理点数
CostFunction=@(x) Cost(x,model);    % 成本函数
nVar=model.n;       % 决策变量的数量
VarSize=[1 nVar];   % 决策变量矩阵的大小
VarMin.x=model.xmin;           % 变量的下界
VarMax.x=model.xmax;           % 变量的上限
VarMin.y=model.ymin;           % 变量的下界
VarMax.y=model.ymax;           % 变量的上限

%% PSO + ABC 参数
MaxIt=150;          % 最大迭代次数
nPop=20;           % 种群规模(群体规模)
w=1;                % 惯性重量
wdamp=0.98;         % 惯性重量阻尼比
c1=1.5;             % 个人学习系数
c2=1.5;             % 全局学习系数
nOnlooker = nPop;         % 围观蜜蜂的数量
L = round(0.6*nVar*nPop); % 放弃限制参数(试用限制)
a = 1;
alpha=0.1;
VelMax.x=alpha*(VarMax.x-VarMin.x);    % 最大速度
VelMin.x=-VelMax.x;                    %最小速度
VelMax.y=alpha*(VarMax.y-VarMin.y);    % 最大速度
VelMin.y=-VelMax.y;                    % 最小速度

%% Initialization PSO + ABC
% 创建空粒子结构
empty_particle.Position=[];
empty_particle.Velocity=[];
empty_particle.Cost=[];
empty_particle.Sol=[];
empty_particle.Best.Position=[];
empty_particle.Best.Cost=[];
empty_particle.Best.Sol=[];
% 空蜂结构
empty_bee.Position = [];
empty_bee.Cost = [];
% 初始化全局最佳
GlobalBest.Cost=inf;
%初始化人口数组
pop = repmat(empty_bee, nPop, 1);
bee = repmat(empty_bee, nPop, 1);
% 初始化有史以来最好的解决方案
BestSol.Cost = inf;
GlobalBest.Cost=inf;
% 创建粒子矩阵
particle=repmat(empty_particle,nPop,1);
% 初始化循环
for i=1:nPop
% 初始化位置
if i > 1
particle(i).Position=CRSolution(model);
else
% 从源头到目的地的直线
xx = linspace(model.xs, model.xt, model.n+2);
yy = linspace(model.ys, model.yt, model.n+2);
particle(i).Position.x = xx(2:end-1);
particle(i).Position.y = yy(2:end-1);
end
%初始化速度
particle(i).Velocity.x=zeros(VarSize);
particle(i).Velocity.y=zeros(VarSize);
%评估
[particle(i).Cost, particle(i).Sol]=CostFunction(particle(i).Position);
% 更新个人最佳
particle(i).Best.Position=particle(i).Position;
particle(i).Best.Cost=particle(i).Cost;
particle(i).Best.Sol=particle(i).Sol;
%更新全球最佳
if particle(i).Best.Cost<GlobalBest.Cost
GlobalBest=particle(i).Best;
end
end
% 初始化循环
for i=1:nPop
% 初始化位置
if i > 1
pop(i).Position=CRSolution(model);
else
% 从源头到目的地的直线
xx = linspace(model.xs, model.xt, model.n+2);
yy = linspace(model.ys, model.yt, model.n+2);
pop(i).Position.x = xx(2:end-1);
pop(i).Position.y = yy(2:end-1);
end



  
 
  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25
  • 26
  • 27
  • 28
  • 29
  • 30
  • 31
  • 32
  • 33
  • 34
  • 35
  • 36
  • 37
  • 38
  • 39
  • 40
  • 41
  • 42
  • 43
  • 44
  • 45
  • 46
  • 47
  • 48
  • 49
  • 50
  • 51
  • 52
  • 53
  • 54
  • 55
  • 56
  • 57
  • 58
  • 59
  • 60
  • 61
  • 62
  • 63
  • 64
  • 65
  • 66
  • 67
  • 68
  • 69
  • 70
  • 71
  • 72
  • 73
  • 74
  • 75
  • 76
  • 77
  • 78
  • 79
  • 80
  • 81
  • 82
  • 83
  • 84
  • 85
  • 86
  • 87
  • 88
  • 89
  • 90
  • 91
  • 92
  • 93
  • 94
  • 95
  • 96
  • 97
  • 98
  • 99
  • 100
  • 101
  • 102
  • 103
  • 104
  • 105
  • 106

三、运行结果

在这里插入图片描述
在这里插入图片描述

四、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1] 包子阳,余继周,杨杉.智能优化算法及其MATLAB实例(第2版)[M].电子工业出版社,2016.
[2]张岩,吴水根.MATLAB优化算法源代码[M].清华大学出版社,2017.
[3]巫茜,罗金彪,顾晓群,曾青.基于改进PSO的无人机三维航迹规划优化算法[J].兵器装备工程学报. 2021,42(08)
[4]方群,徐青.基于改进粒子群算法的无人机三维航迹规划[J].西北工业大学学报. 2017,35(01)
[5]周国春,肖本贤.基于人工蜂群算法的机器人路径规划研究[J].自动化与仪表. 2018,33(05)

3 备注
简介此部分摘自互联网,仅供参考,若侵权,联系删除

文章来源: qq912100926.blog.csdn.net,作者:海神之光,版权归原作者所有,如需转载,请联系作者。

原文链接:qq912100926.blog.csdn.net/article/details/123996080

【版权声明】本文为华为云社区用户转载文章,如果您发现本社区中有涉嫌抄袭的内容,欢迎发送邮件进行举报,并提供相关证据,一经查实,本社区将立刻删除涉嫌侵权内容,举报邮箱: cloudbbs@huaweicloud.com
  • 点赞
  • 收藏
  • 关注作者

评论(0

0/1000
抱歉,系统识别当前为高风险访问,暂不支持该操作

全部回复

上滑加载中

设置昵称

在此一键设置昵称,即可参与社区互动!

*长度不超过10个汉字或20个英文字符,设置后3个月内不可修改。

*长度不超过10个汉字或20个英文字符,设置后3个月内不可修改。