蟻群算法最短路徑MATLAB直接調(diào)用

1、function R_best,L_best,L_ave,Shortest_Route,Shortest_Length=ACATSP(C,NC_ max,m,Alpha,Beta,Rho,Q) %=% ACATSP.m% Ant Colony Algorithm for Traveling Salesman Problem% ChengAihua,PLA Information Engineering University,ZhengZhou,China% Email:aihuacheng% All rights reserved% 主要符號說明%Cn個城市的坐標，nx 2的矩陣% NC_ma

2、x 最大迭代次數(shù)%m螞蟻個數(shù)%Alpha表征信息素重要程度的參數(shù)%Beta表征啟發(fā)式因子重要程度的參數(shù)%Rho信息素蒸發(fā)系數(shù)%Q信息素增加強度系數(shù)% R_best 各代最佳路線% L_best 各代最佳路線的長度%第一步：變量初始化n=size(C,1);%*表示問題的規(guī)模(城市個數(shù))D=zeros( n,n );%D表示完全圖的賦權(quán)鄰接矩陣for i=1:nfor j=1:nif i=jD(i,j)=(C(i,1)-C(j,1)八2+(C(i,2)-C(j,2)八2)95;elseD(i,j)=eps;endD(j,i)=D(i,j);endendEta=1./D;%Eta為啟發(fā)因子，這里設(shè)

3、為距離的倒數(shù)Tau=o nes( n,n );%Tau為信息素矩陣Tabu=zeros(m, n);%存儲并記錄路徑的生成NC=1;%迭代計數(shù)器R_best=zeros(NC_max, n);%各代最佳路線L_best=i nf.*o nes(NC_max,1);% 各代最佳路線的長度L_ave二zeros(NC_max,1);%各代路線的平均長度 while NC=rand);to_visit=J(Select(1);Tabu(i,j)=to_visit;endendif NC=2Tabu(1,:)=R_best(NC-1,:);end%第四步：記錄本次迭代最佳路線L=zeros(m,1);

4、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 i=1:mfor j=1:(n-1)Delta_Tau(Tabu(i,j),Tabu(i,j+1)=Delta_Tau(Tabu(i,j),Tabu(i,j+1)+Q/L

5、(i); endDelta_Tau(Tabu(i,n),Tabu(i,1)=Delta_Tau(Tabu(i,n),Tabu(i,1)+Q/L(i);endTau=(1-Rho).*Tau+Delta_Tau;%第六步：禁忌表清零Tabu=zeros(m,n);end%第七步：輸出結(jié)果Pos=find(L_best=min(L_best);Shortest_Route=R_best(Pos(1),:);Shortest_Length=L_best(Pos(1);disp( 最短距離 : num2str(Shortest_Length);disp( 最短路徑 : num2str(Shortest

6、_RouteShortest_Route(1);figureDrawRoute(C,Shortest_Route)for i = 1:size(C,1)text(C(i,1),C(i,2), num2str(i);endxlabel( 城市位置橫坐標 )ylabel( 城市位置縱坐標 ) grid onfigure plot(L_best)xlabel( 迭代次數(shù) )ylabel(距離)grid onhold onplot(L_ave)title( 平均距離和最短距離對比 )function DrawRoute(C,R)%=%DrawRoute.m%畫路線圖的子函數(shù)%-%CCoordinate節(jié)點坐標，由一個NX 2的矩陣存儲%=RRoute路線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