基于粒子群算法的TSP问题研究 - 图文(13)

西安工业大学毕业设计(论文)

%------进入主要循环，按照公式依次迭代------------

for t=1:M

for i=1:N

v(i,:)=w*v(i,:)+c1*rand*(y(i,:)-x(i,:))+c2*rand*(pg-x(i,:));

x(i,:)=x(i,:)+v(i,:);

if fitness(x(i,:))

p(i)=fitness(x(i,:));

y(i,:)=x(i,:);

end

if p(i)

pg=y(i,:);

end

end

Pbest(t)=fitness(pg); end

xm = pg';

fv = fitness(pg); 函数2：

function F= fitness(x) F=0;

for i=1:30

F=F+x(i)^2; end 程序2

城市的位置输入到计算机中用下面的函数： function [DLn,cityn]=tsp(n)

if n==10

35

西安工业大学毕业设计(论文)

city10=[0.4 0.4439;0.2439 0.1463;0.1707 0.2293;0.2293 0.761;0.5171 0.9414; 0.8732 0.6536;0.6878 0.5219;0.8488 0.3609;0.6683 0.2536;0.6195 0.2634]; cities d'=2.691 for i=1:10

for j=1:10

DL10(i,j)=((city10(i,1)-city10(j,1))^2+(city10(i,2)-city10(j,2))^2)^0.5; end end

DLn=DL10; cityn=city10; end

适应度函数：

function F=SumDistance(dislist,s)

DistanV=0; n=size(s,2); for i=1:(n-1)

DistanV=DistanV+dislist(s(i),s(i+1));%此函数用来计算总路程 end

DistanV=DistanV+dislist(s(n),s(1));

F=DistanV; 绘图：

function m=plotTSP10(Clist,BSF,bsf,p,f)

CityNum=size(Clist,1);

for i=1:CityNum-1 axis([0,1,0,1]); %绘制点的连线图

plot([Clist(BSF(i),1),Clist(BSF(i+1),1)],[Clist(BSF(i),2),Clist(BSF(i+1),2)],'rs-','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','g');

hold on;%'rs-','LineWidth',3,'MarkerEdgeColor','','MarkerFaceColor','g'-表示线为红色实线线宽为2，点为方形绿点边缘为黑色 end

axis([0,1,0,1]);

%绘制最后一个点和起始点的连线

plot([Clist(BSF(CityNum),1),Clist(BSF(1),1)],[Clist(BSF(CityNum),2),Clist(BSF(1),2)],'ms-','LineWidth',2,'MarkerEdgeColor','k','MarkerFaceColor','g');

36

西安工业大学毕业设计(论文)

%'ms-','LineWidth',3,'MarkerEdgeColor','','MarkerFaceColor','g'-表示线为洋红色实线线宽为2，点为方形绿点边缘为黑色 title([num2str(CityNum),'TSP']); if f==0

text(0.1,0.1,['迭代次数 ',int2str(p),' 最短距离 ',num2str(bsf)]); end hold off;

pause(0.05); 初始化和变换：

function BasePSOforTSP %初始化