%% 蚁群算法（ant colony optimization,ACO）
%清空变量
clear
close all
clc
 
[ graph ]  = createGraph();
 
figure 
subplot(1,3,1)
drawGraph( graph); 
%% 初始化参数
maxIter = 100;
antNo = 50;
tau0 = 10 * 1 / (  graph.n * mean( graph.edges(:)  )  );  
tau = tau0 * ones( graph.n , graph.n); % Phromone matirx 
eta = 1./ graph.edges;  % 每个边的可取性
rho = 0.5; % 蒸发率
alpha = 1;  % Phromone指数参数
beta = 1;  % Desirability exponetial paramter
%% 循环
bestFitness = inf;
bestTour = [];
for t = 1 : maxIter
    % Create Ants 
    colony = [];
    colony = createColony( graph, colony , antNo, tau, eta, alpha,  beta);
    % Calculate the fitness values of all ants 
    for i = 1 : antNo