%TSP問題蟻群算法

global NC;      %迭代次數(shù)
global city_n;  %城市數(shù)量
global dis_table;   %城市距離矩陣
global G;       %記錄進化代數(shù)
global everbest;    %歷代最優(yōu)解
global hu_table;    %啟發(fā)式分布表
global adapt_best;

%=================================================================
%第一種群數(shù)據(jù)
global a_A;       %信息指數(shù)
global b_A;       %啟發(fā)式指數(shù)
global tobu_A;        %螞蟻禁忌表
global ph_table_A;    %信息素分布表
global ant_n_A;       %螞蟻數(shù)量
global temp_pool_A;   %放置備選城市
global dispose_A;     %信息素揮發(fā)率
global best_individual_A; %歷代最優(yōu)解路徑
global adapt_ave_A;   %進化數(shù)據(jù)
global adapt_best_A;  %歷代最優(yōu)解的變化
global everbest_A;    %歷代最優(yōu)解
%=================================================================
%第二種群數(shù)據(jù)
global a_B;       %信息指數(shù)
global b_B;       %啟發(fā)式指數(shù)
global tobu_B;        %螞蟻禁忌表
global ph_table_B;    %信息素分布表
global ant_n_B;       %螞蟻數(shù)量
global temp_pool_B;   %放置備選城市
global dispose_B;     %信息素揮發(fā)率
global best_individual_B; %歷代最優(yōu)解路徑
global adapt_ave_B;   %進化數(shù)據(jù)
global adapt_best_B;  %歷代最優(yōu)解的變化
global everbest_B;    %歷代最優(yōu)解
%=================================================================


initial;    %初始化數(shù)據(jù)，將螞蟻隨機放入城市位置

for G=1:NC
    pause(0.01);
    search;         %解構(gòu)成
    if mod(G,10)==0
        communication;  %每10次迭代進行一次通信
    end
    adapting;       %計算各回路長度
    keepbest;       %保優(yōu)函數(shù)
    ph_fresh;       %信息素更新
    paint;          %繪制圖形
end

result;

clear a;
clear b;
clear a_A;
clear a_B;
clear b_A;
clear b_B;
clear NC;
clear h;
clear i;
clear j;
clear k;
clear n;
clear min_pos;
clear min_dis;
clear r;
clear tobu_A;
clear tobu_B;
clear dis_sum;
clear dispose_A;
clear dispose_B;
clear dis;
clear ant_n_A;
clear ant_n_B;
clear x1;
clear x2;
clear y1;
clear y2;
clear temp_pool_A;
clear temp_pool_B;
clear city;
clear city_n;
clear adapt_A;
clear adapt_B;
clear ada_sum;
clear ada_temp;
clear G;
clear ans;
clear everbest_A;
clear everbest_B;
clear dis_table;
clear best_individual;
clear best_individual_A;
clear best_individual_B;