蚁群算法

蚁群算法

算法概述

蚁群算法是一种用来寻找优化路径的概率型算法。它由Marco Dorigo于1992年在他的博士论文中提出，其灵感来源于蚂蚁在寻找食物过程中发现路径的行为。 这种算法具有分布计算、信息正反馈和启发式搜索的特征，本质上是进化算法中的一种启发式全局优化算法。

算法思想

当蚂蚁沿着一条路到达终点以后会立马饭回来，这样，短的路蚂蚁来回一次的时间就断短，这也意味着重复的频率就快，因而在单位时间走过的蚂蚁数目就多，洒下的信息素自然也会多，自然会有更多的蚂蚁被吸引过来，从而洒下更多的信息素….；而长的路正好相反，因此，越来越多的蚂蚁聚集到较短的路径上来，最短的路径就近似找到了。

算法解析

蚂蚁如何找到最短路径
1、信息素:信息素多的地方显然经过这里的蚂蚁多，因而会有更多的蚂蚁聚集过来。
2、正反馈现象：某一路径上走过的蚂蚁越多，则后来者选择该路径的概率就越大。
蚂蚁k（k=1,2,…,m）根据各个城市间连接路径上的信息素浓度和到邻接城市的距离决定其下一个访问城市，设Pkij(t)表示t时刻蚂蚁k从城市i转移到城市j的概率，其计算公式为：

信息素更新公式为：

算法实现

%% 旅行商问题(TSP)优化
%% 清空环境变量
clear all
clc

%% 导入数据
load citys_data.mat

%% 计算城市间相互距离
fprintf('Computing Distance Matrix... \n');
n = size(citys,1);
D = zeros(n,n);
for i = 1:n
    for j = 1:n
        if i ~= j
            D(i,j) = sqrt(sum((citys(i,:) - citys(j,:)).^2));
        else
            D(i,j) = 1e-4;      
        end
    end    
end

%% 初始化参数
fprintf('Initializing Parameters... \n');
m = 50;                              % 蚂蚁数量
alpha = 1;                           % 信息素重要程度因子
beta = 5;                            % 启发函数重要程度因子
rho = 0.1;                           % 信息素挥发因子
Q = 1;                               % 常系数
Eta = 1./D;                          % 启发函数
Tau = ones(n,n);                     % 信息素矩阵
Table = zeros(m,n);                  % 路径记录表
iter = 1;                            % 迭代次数初值
iter_max = 150;                      % 最大迭代次数 
Route_best = zeros(iter_max,n);      % 各代最佳路径       
Length_best = zeros(iter_max,1);     % 各代最佳路径的长度  
Length_ave = zeros(iter_max,1);      % 各代路径的平均长度  

%% 迭代寻找最佳路径
figure;
while iter <= iter_max
    fprintf('迭代第%d次\n',iter);
    % 随机产生各个蚂蚁的起点城市
      start = zeros