function [ distances ] = calculateDistance( city )
%計算距離
[~, col] = size(city);
distances = zeros(col);
for i=1:col
for j=1:col
distances(i,j)= distances(i,j)+ sqrt( (city(1,i)-city(1,j))^2 + (city(2,i)-city(2,j))^2 );
end
end
end
交叉函數：crossover.m
function [childPath] = crossover(parent1Path, parent2Path, prob)
% 交叉
random = rand();
if prob >= random %行：随機生成一個01之前的數，和機率值比較，相當于一個輪盤賭的選擇，來執行機率選擇執行操作
[l, length] = size(parent1Path);
childPath = zeros(l,length);
setSize = floor(length/2) -1; %行：floor 朝負無窮大方向取整
offset = randi(setSize); %行：随機 産生1-setSize之間的整數
parent1OXstart=offset; %行：0X交叉，父代1遺傳給子代的片段起點
parent1OXend=setSize+offset-1; %行：0X交叉，父代1遺傳給子代的片段終點
for i=parent1OXstart:parent1OXend
childPath(1,i) = parent1Path(1,i);
end
childPath(1,1) = parent1Path(1,1);%行：将路徑的起點和終點直接賦給子代，因為起點和終點不變
childPath(1,length) = parent1Path(1,length);
%行：parent2依順序放入子代位坐标
parent2Point=2;
for x=2:length-1
if childPath(1,x)==0
for y=parent2Point:length-1
if ~any(childPath == parent2Path(1,y))%行：如果子代路徑相應元素不等于父代元素
childPath(1,x)=parent2Path(1,y);
parent2Point=y+1;
break;
end
end
end
end
else
childPath = parent1Path;
end
end
變異函數：mutate.m
function [ mutatedPath ] = mutate( path, prob )
%對指定的路徑利用指定的機率進行更新
%行：由于路徑不予許存在相同的數，是以變異不能随意變，隻能用随機互換位置的方法
random = rand();
if random <= prob
[l,length] = size(path);
index1 = randi([2,length-1]);
index2 = randi([2,length-1]);
%交換
temp = path(l,index1);
path(l,index1) = path(l,index2);
path(l,index2)=temp;
end
mutatedPath = path;
結果路徑圖檔顯示函數：paint.m
function [ bestPopPath ] = paint( cities, pop, minPath, totalDistances,gen)
gNumber=gen;
[~, length] = size(cities);
xDots = cities(1,:);
yDots = cities(2,:);
%figure(1);
%行：先畫城市的固定坐标
plot(xDots,yDots, 'p', 'MarkerSize', 10, 'MarkerFaceColor', 'green');
xlabel('經度');
ylabel('緯度');
axis equal
[num,txt,raw]=xlsread('provincial_ capital_coordinate12.xlsx'); %經緯度資料%
text(xDots(1)+0.2,yDots(1)+0.2,txt(1+1,2),'FontSize',9.5,'Color','red','FontWeight','Bold') ; %先畫起點，突出起點
for i=2:length-1
text(xDots(i)+0.2,yDots(i)+0.2,txt(i+1,2),'FontSize',8,'Color','blue','FontWeight','Bold') ; %加上0.01使标号和點不重合，可以自己調整
end
hold on
[minPathX,~] = find(totalDistances==minPath,1, 'first');
bestPopPath = pop(minPathX, :);
bestX = zeros(1,length);
bestY = zeros(1,length);
for j=1:length
bestX(1,j) = cities(1,bestPopPath(1,j));
bestY(1,j) = cities(2,bestPopPath(1,j));
end
% title('目前代數%d 的路徑圖');
title(['第', num2str(gNumber),'代種群的最優路徑圖']) ;
%行：再畫城市的路徑連線
plot(bestX(1,:),bestY(1,:), 'red', 'LineWidth', 1.25);
legend('城市', '路徑');
axis equal%行：axis equal/将橫軸縱軸的定标系數設成相同值。
grid on
%text(5,0,sprintf('疊代次數: %d 總路徑長度: %.2f',gNumber, minPath),'FontSize',10);
drawnow
hold off
end
