adapting.m

来自「此程序包是用双种群蚁群算法来求中国75个城市的最短路径问题」· M 代码 · 共 65 行

%回路长度计算
%=========================================
%第一种群计算
t=load('75城市TSP问题数据');
for i=1:ant_n_A
    dis=0;
    dis_sum=0;
    x1=0;
    x2=0;
    y1=0;
    y2=0;
    
    for j=1:(city_n-1)      %计算前段距离
        x1=t.pos(tobu_A(i,j),1);     %相邻两点的坐标
        y1=t.pos(tobu_A(i,j),2);
        x2=t.pos(tobu_A(i,j+1),1);
        y2=t.pos(tobu_A(i,j+1),2);
        dis=sqrt((x1-x2)^2+(y1-y2)^2);
        dis_sum=dis_sum+dis;    %距离累积
    end
    % 计算最后一个城市回到第一个城市的距离
    x1=t.pos(tobu_A(i,1),1);
    y1=t.pos(tobu_A(i,1),2);
    x2=t.pos(tobu_A(i,city_n),1);
    y2=t.pos(tobu_A(i,city_n),2);
    dis=sqrt((x1-x2)^2+(y1-y2)^2);
    dis_sum=dis_sum+dis;    %距离累积
    
    adapt_A(i)=dis_sum;
end

adapt_ave_A(G)=mean(adapt_A);

%============================================
%第二种群计算
for i=1:ant_n_B
    dis=0;
    dis_sum=0;
    x1=0;
    x2=0;
    y1=0;
    y2=0;
    
    for j=1:(city_n-1)      %计算前段距离
        x1=t.pos(tobu_B(i,j),1);     %相邻两点的坐标
        y1=t.pos(tobu_B(i,j),2);
        x2=t.pos(tobu_B(i,j+1),1);
        y2=t.pos(tobu_B(i,j+1),2);
        dis=sqrt((x1-x2)^2+(y1-y2)^2);
        dis_sum=dis_sum+dis;    %距离累积
    end
    % 计算最后一个城市回到第一个城市的距离
    x1=t.pos(tobu_B(i,1),1);
    y1=t.pos(tobu_B(i,1),2);
    x2=t.pos(tobu_B(i,city_n),1);
    y2=t.pos(tobu_B(i,city_n),2);
    dis=sqrt((x1-x2)^2+(y1-y2)^2);
    dis_sum=dis_sum+dis;    %距离累积
    
    adapt_B(i)=dis_sum;
end

adapt_ave_B(G)=mean(adapt_B);