select_for_as.m

蚂蚁系统三种模型：蚁量模型、蚁密模型、蚁圈模型的源代码。

function  j=Select_for_as(Ak,i,Allow,A,B,Pheromone,Heuristic)
AllowCity=find(Allow(Ak,:)~=0);
BpijP=(Pheromone(i,AllowCity).^A);
BpijH=(Heuristic(i,AllowCity).^B);
Bpij=BpijP.*BpijH;
Gpij=sum(Bpij);
Pij=Bpij./Gpij;
Pij=cumsum(Pij);
Rd=rand;% 赌轮选择
[Anx Any]=size(AllowCity);
for FitIn=1:Any 
       if Rd<=Pij(FitIn)
          j=AllowCity(FitIn);
          break;
       end
end