【蚁群系统】解决TSP旅行商问题(matlab实现)
目录
- 1. 蚁群系统(Ant Colony System)(ACS)
- 2. 旅行商问题
- 3. 蚁群系统解决TSP问题(Matlab)
1. 蚁群系统(Ant Colony System)(ACS)
1997年,蚁群算法的创始人Dorigo在“Ant colony system: a cooperative learning approach to the traveling salesman problem”一文中提出了一种具有全新机制的ACO算法——蚁群系统(Ant Colony System,ACS),进一步提高了ACO算法的性能。
(1)使用一种伪随机比例规则(pseudorandom proportional)选择下城市节点,建立开发当前路径与探索新路径之间的平衡。
- q0是一个[0, 1]区间内的参数,当产生的随机数q≤q0时,蚂蚁直接选择使启发式信息与信息素量的指数乘积最大的下城市节点,我们通常称之为开发(exploitation);反之,当产生的随机数q>q0时ACS将和各种AS算法一样使用轮盘赌选择策略,我们称之为偏向探索(bias exploration)。通过调整q0,我们能有效调节“开发”与“探索”之间的平衡,以决定算法是集中开发最优路径附近的区域,还是探索其它的区域。
(2)使用信息素全局更新规则,每轮迭代中所有蚂蚁都已构建完路径后,在属于至今最优路径的边上蒸发和释放信息素。
- 不论是信息素的蒸发还是释放,都只在属于至今最优路径的边上进行,这里与AS有很大的区别。因为AS算法将信息素的更新应用到了系统的所有边上,信息素更新的计算复杂度为O(n^2),而ACS算法的信息素更新计算复杂度降低为O(n)。参数ρ 代表信息素蒸发的速率,新增加的信息素 被乘上系数ρ 后,更新后的信息素浓度被控制在旧信息素量与新释放的信息素量之间,用一种隐含的又更简单的方式实现了MMAS算法中对信息素量取值范围的限制。
(3)引入信息素局部更新规则,在路径构建过程中,对每一只蚂蚁,每当其经过一条边(i, j)时,它将立刻对这条边进行信息素的更新。
- 信息素局部更新规则作用于某条边上会使得这条边被其他蚂蚁选中的概率减少。这种机制大大增加了算法的探索能力,后续蚂蚁倾向于探索未被使用过的边,有效地避免了算法进入停滞状态。
- 顺序构建和并行构建。顺序构建是指当一只蚂蚁完成一轮完整的构建并返回到初始城市之后,下一只蚂蚁才开始构建;并行构建是指所有蚂蚁同时开始构建,每次所有蚂蚁各走一步(从当前城市移动到下一个城市)。对于ACS,要注意到两种路径构建方式会造成算法行为的区别。在ACS中通常我们选择让所有蚂蚁并行地工作。
总结:
- ACS解决了基本蚁群算法在构造解过程汇总,随机选择策略造成的算法进化速度慢的缺点。该算法在每一次循环中仅让最短路径上的信息量作更新,且以较大的概率让信息量最大的路径被选中,充分利用学习机制,强化最优信息的反馈。ACS的核心思想是:蚂蚁在寻找最佳路径的过程中只能使用局部信息,即采用局部信息对路径上的信息量进行调整;在所有进行寻优的蚂蚁结束路径的搜索后,路径上的信息量会再一次调整,这次采用的是全局信息,而且只对过程中发现的最好路径上的信息量进行加强。
ACS模型与AS模型的主要区别有三点:
(1)蚂蚁的状态转移规则不同
(2)全局更新规则不同
(3)新增了对各条路径信息量调整的局部更新规则
2. 旅行商问题
假设有一个旅行商人要拜访n个城市,他必须选择所要走的路径,路径的限制是每个城市只能拜访一次,而且最后要回到原来出发的城市。路径的选择目标是要求得的路径路程为所有路径之中的最小值。这就是TSP,属于NPC问题。
3. 蚁群系统解决TSP问题(Matlab)
function [R_best,L_best,L_ave,Shortest_Route,Shortest_Length]=ACS(NC_max,m,Alpha,Beta,Rho,Q)
clear all;
close all;
clc;
% 重要参数说明:
% R_best 各代最佳路线
% L_best 各代最佳路线的长度%% 第一步:变量初始化
C=[5.294,1.558;4.286,3.622;4.719,2.774;4.185,2.230;0.915,3.821;4.771,6.041;1.524,2.871;3.447,2.111;3.718,3.665;2.649,2.556];%n个城市的坐标,n×2的矩阵
n=size(C,1); % 表示城市个数n=10
D=zeros(n,n); % D表示完全图的赋权邻接矩阵,产生n*n的值全为0的矩阵
for i=1:nfor j=1:nD(i,j)=((C(i,1)-C(j,1))^2+(C(i,2)-C(j,2))^2)^0.5;end
endm = 50; % 蚂蚁个数 50
q0=0.9; % 参数q0
NC_max = 100; % 最大迭代次数 100
Alpha = 2; % 表征信息素重要程度的参数2
Beta = 4; % 表征启发式因子重要程度的参数4
Rho = 0.1; % 信息素蒸发系数0.1
Q = 0.01; % 信息素增加强度系数 0.01
Eta=1./D; % Eta为启发因子,为距离的倒数
Lnn=10;
TauStart=(1/(n*Lnn)).*ones(n,n); %TauStart为信息素有个初始值 Tau:产生10*10行矩阵
Tabu=zeros(m,n); % 存储并记录路径 Tabu:50*10矩阵
NC=1; % 迭代计数器
R_best=zeros(NC_max,n); % 各代最佳路线 100*10
L_best=inf.*ones(NC_max,1); % 各代最佳路线的长度 100*1
L_ave=zeros(NC_max,1); % 各代路线的平均长度 100*1%% 第二步: 将m=50只蚂蚁放到n=10个城市上
while NC<=NC_max Randpos=[];Tau=TauStart;for i=1:(ceil(m/n)) % m为50,n为10Randpos=[Randpos,randperm(n)]; % 将每只蚂蚁放到随机的城市endTabu(:,1)=(Randpos(1,:)); % 初始化禁忌表中第一个元素 矩阵的第一行赋值给另外一个矩阵的第一列%% 第三步:对每只蚂蚁按概率函数选择下一座城市,完成各自的周游for j=2:nfor i=1:m% 选择下一个访问的城市visited=Tabu(i,1:(j-1)); % 在访问第j个城市的时候,第i个蚂蚁访问过的城市, visited也是矩阵unvisited=zeros(1,(n-j+1)); % 待访问的城市 P=unvisited; % 待访问城市的选择概率分布 J、P都是一行值为0的矩阵Jc=1;% 找出未访问的城市 for k=1:n if length(find(visited==k))==0unvisited(Jc)=k;Jc=Jc+1;endend% 计算待选择城市的概率for k=1:length(unvisited) % Tau(visited(end),unvisited(k))访问过的城市的最后一个与所有未访问的城市之间的信息素P(k)=(Tau(visited(end),unvisited(k))^Alpha)*(Eta(visited(end),unvisited(k))^Beta);end% 状态转移规则q=rand();if q<=q0 Select=find(P==max(P));else % 轮盘赌P=P/(sum(P));Pcum=cumsum(P);Select=find(Pcum>=rand);end% 添加到禁忌表to_visit=unvisited(Select(1));Tabu(i,j)=to_visit;end% 信息素局部更新Tau=(1-Rho).*Tau+Rho.*TauStart; endif NC>=2Tabu(1,:)=R_best(NC-1,:); % 上一代最短路径作为本代第一条路径end%% 第四步:记录本次迭代最佳路线L=zeros(m,1);for i=1:mR=Tabu(i,:);for j=1:(n-1)L(i)=L(i)+D(R(j),R(j+1)); % 要走一周回到原来的地点 endL(i)=L(i)+D(R(1),R(n));endL_best(NC)=min(L); % 记录每一代中路径的最短值pos=find(L==L_best(NC));R_best(NC,:)=Tabu(pos(1),:); % 最优的路径L_ave(NC)=mean(L); % 每一代中平均路径长度NC=NC+1;%% 第五步:更新全局信息素(增加了全局最短路径的信息素)Delta_Tau=zeros(n,n);for j=1:(n-1)Delta_Tau(Tabu(pos(1),j),Tabu(pos(1),j+1))=Delta_Tau(Tabu(pos(1),j),Tabu(pos(1),j+1))+1/min(L);endDelta_Tau(Tabu(pos(1),n),Tabu(pos(1),1))=Delta_Tau(Tabu(pos(1),n),Tabu(pos(1),1))+1/min(L);Tau=(1-Rho).*Tau+Rho.*Delta_Tau; % 更新路径上的信息素含量%%第六步:禁忌表清零Tabu=zeros(m,n);
end%% 第七步:输出结果
Pos=find(L_best==min(L_best));
Shortest_Route=R_best(Pos(1),:);
Shortest_Length=L_best(Pos(1));
% 画图
subplot(1,2,1)
DrawRoute(C,Shortest_Route)
subplot(1,2,2)
plot(L_best,'b')
hold on
plot(L_ave,'r')
title('平均距离和最短距离')disp('最短距离为:')
disp(min(L_best))num=find(min(L_best)==L_best)
disp('最短路径是:')
disp(R_best(num(1),:))%% 画路线图
function DrawRoute(C,R)% C 节点坐标,由一个N×2的矩阵存储
% R 路线N=length(R);
scatter(C(:,1),C(:,2));
hold on
plot([C(R(1),1),C(R(N),1)],[C(R(1),2),C(R(N),2)],'g')
hold on
for ii=2:Nplot([C(R(ii-1),1),C(R(ii),1)],[C(R(ii-1),2),C(R(ii),2)],'g')hold on
end
title('旅行商问题优化结果');
【蚁群系统】解决TSP旅行商问题(matlab实现)相关推荐
- 蚁群算法解决tsp问题c语言,蚁群算法解决TSP问题程序.doc
蚁群算法解决TSP问题程序 蚁群算法用于求解TSP问题,经过仿真测试,发现此程序的优化效率和鲁棒性都非常好. 这与在无线多媒体传感器网络路由算法应用到的寻找最佳路径的蚁群算法非常相似. functio ...
- 蚁群算法解决tsp问题python_蚁群算法在解决TSP问题中的应用
陈灵佳 文章首先对蚁群算法与TSP问题进行简要介绍,在此基础上对蚁群算法在解决TSP问题中的应用进行论述.期望通过本文的研究能够对TSP问题的解决有所帮助. [关键词]蚁群算法 TSP问题 最优解 1 ...
- 蚁群算法解决 TSP 问题
蚁群算法解决 TSP 问题 数据集 Tools.py Ant.py ACO_G.py 运行效果 数据集 json 形式(c.json)的中国各省市经纬度数据集,一共 2241 个市的数据,为后来的 T ...
- 机器学习(MACHINE LEARNING)MATLAB蚁群算法解决TSP问题
文章目录 1 蚁群算法 2 蚁群算法与TSP问题的关系 3 代码实现 1 蚁群算法 基本原理: (1)蚂蚁在携带等量的信息素一路释放 (2)信息素浓度会和路径的长度成反比 (3)下次蚂蚁来到该路口会选 ...
- 蚂蚁算法求解tsp问题matlab,蚁群算法解决TSP问题的MATLAB程序
蚁群算法TSP(旅行商问题)通用matlab程序 function [R_best,L_best,L_ave,Shortest_Route,Shortest_Length]=ACATSP(C,NC_m ...
- 蚁群算法求解TSP问题 matlab程序
https://blog.csdn.net/wayjj/article/details/72809344 蚁群算法,单单学习算法还是不够深入了解,得实际编程实现了,理解才能更加透彻,本文根据这篇博文贴 ...
- 蚁群算法解决TSP问题
原文转载于https://www.cnblogs.com/biaoyu/archive/2012/09/26/2704456.html 蚁群算法简介 蚁群算法(Ant Clony Optimizati ...
- 蚁群算法解决TSP问题(2#JAVA代码+详细注释+对比动态规划【JAVA】)
第一部分:原理 TSP10cities.txt 1 2066 2333 2 935 1304 3 1270 200 4 1389 700 5 984 2810 6 2253 478 7 949 302 ...
- 蚁群算法解决TSP问题(matlab代码实现)
最新文章
- 编程自动化,未来机器人将超越人类?
- python程序员最害怕的东西,面试!恭喜你收到了offer,来看下吧
- 下一代脑电图可以帮助恢复失去的大脑功能
- java poi设置单元格格式为数值,Apache POI 如何读取Excel中数值类型单元格所规定的保留小数位?...
- BZOJ 3160 FFT+Manacher
- SpringSecurity 案例之创建资源服务器准备工作
- js 技巧杂引(转)
- [导入]C#好书盘点【月儿原创】
- 线程间通信之eventfd
- 高级java技术web组件_(重温)JavaWeb--Servlet技术(二)(JavaWeb 的三大组件之一)...
- python自动化常用数字_Python全栈自动化系列之Python编程基础(基本数据类型)
- ServerGuide U盘安装Windows Server 2008 R2安装
- 神经网络浅讲:从神经元到深度学习
- 写给零基础小白的网站开发入门
- 重学JavaWeb —— JSP,简单全面一发入魂
- warning: TCG doesn‘t support requested feature: CPUID.01H:ECX.vmx [bit 5]Could not initialize SDL
- 3.7 使用吸管和颜色取样器工具 [Ps教程]
- The Meaning of Life
- SAP 生产订单工序外协对应的工序报工后,对应的采购订单自动入库
- 计算机中word的快捷键大全,Word中全选的快捷键是什么 word快捷键大全【图解】...