tsp问题MATLAB原程序

%GATSP.M***************************** %图论作业主程序 clear all close all clc popsize=200; %群体大小 pc=0.9; %交叉概率 pm=0.2; %变异概率 %城市坐标************************ city=[640,600;830,690;870,800;710,790; 580,60;510,670;370,840;350,680; 410,940;20,990;70,640;400,340;180,870; 870,70;910,380;830,460;720,440; 40,500;900,200;240,420;250,20;110,206;990,970;20,100;560,220 ]; ncity=size(city,1); pop=initpop(popsize,ncity); %随机产生初始群体 distcity=distance(city,ncity);%计算城市距离 for i=1:1000%遗传代数 newpop=crossrunning( pop,pc,ncity);%交叉 newpop = mutation ( newpop,ncity,pm);%变异 pop=newpop;%newpop [fitvalue,slength,length]=fit(pop,distcity,popsize,ncity);%计算群体中个体的适应度 [bestindividual,bestfit,lestlength]=best(pop,fitvalue,slength);%计算种群最大个体适应度 g(i)=lestlength;%最小值 tsp(i,:)=bestindividual; n=i; subplot(1,2,1); plot(n,g(n),'.g','LineWidth',3) hold on [selected]=maxselection(pop,fitvalue);%贪婪算子赌轮盘选择复制 pop=selected; end xlabel('\bf遗传代数'); ylabel('\bf函数最小值'); title('\slTSP问题'); [c,index]=min(g);%最小值 dy=c; hh=index;%最小值出现的代数 plot(hh,dy,'r*','LineWidth',6) dd=tsp(hh,:); subplot(1,2,2); scity=1:1:ncity; strcity=mat2str(scity); for lp=1:ncity plot(city(lp,1),city(lp,2),'o','Color','r','LineWidth',3); text(city(lp,1)+3,city(lp,2)+3,mat2str(scity(1,lp)),'FontSize',10) hold on end ylabel('\bf城市坐标'); title('\sl旅行路线'); bsetcityroad(1,:)=city(dd(ncity),:); for i=1:ncity bsetcityroad(i+1,:)=city(dd(i),:); end line(bsetcityroad(:,1),bsetcityroad(:,2),'Color','b','LineWidth',2);