转载请注明出处:http://blog.youkuaiyun.com/microlyc
由 http://blog.youkuaiyun.com/ibelievefly/article/details/47955473 旅行商问题改写;
function ant()
Q=400;%Q表示信息素增强的系数
%nc_now表示当前的迭代次数
nc_max=9999; %nc_max表示自主设置的最大迭代时间 表示终止条件
%first_address表示测试数据中的城市坐标
%RouteOfBest表示各代的最佳路线
%LengthOfBest表示每次迭代的最佳路线的长度
a = 3;%a=0.7表示启发因子,表示信息素的相对重要程度!!!!!!!!!!!很重要
p = 0.1;%p表示信息素的蒸发系数,(1-p)表示信息素持久性系数
%%第一步:变量初始化
SumOfant=30;%SumOfant表示蚂蚁数量
first_address = [%第一个为家,最后一个是食物源
0 3
0 2
-10 0
1 0
6 0
0 -2
0 -3
];%first_address表示测试数据中的城市坐标
target = [1,length(first_address)];
SumOfCity = size(first_address,1);
length_address =inf.*ones(SumOfCity,SumOfCity);%length_address表示两两城市间的距离
length_address(1,2)= Calculation_length(first_address,1,2);
length_address(2,3)= Calculation_length(first_address,2,3);
length_address(3,6)= Calculation_length(first_address,3,6);
length_address(6,7)= Calculation_length(first_address,6,7);
length_address(4,6)= Calculation_length(first_address,4,6);
length_address(2,4)= Calculation_length(first_address,2,4);
length_address(5,6)= Calculation_length(first_address,5,6);
length_address(2,5)= C