TSP遗传算法Matlab程序

很通用的遗传算法matlab程序,希望给大家有所帮助。

TSP遗传算法Matlab程序

function [R,Rlength]=geneticTSP(D,a,n,C,m,alpha)

[N,NN]=size(D);

farm=zeros(n,N);%用于存储种群

for i=1:n

farm(i,:)=randperm(N);%随机生成初始种群

end

R=farm(1,:);

subplot(1,3,1)

scatter(a(:,1),a(:,2),'x')

pause(1)

subplot(1,3,2)

plotaiwa(a,R)

pause(1)

farm(1,:)=R; end

maxlen=max(len);

minlen=min(len);

fitness=fit(len,m,maxlen,minlen);%计算归一化适应值

rr=find(len==minlen);

R=farm(rr(1,1),:);%更新最短路径

FARM=farm;%优胜劣汰，nn记录了复制的个数

nn=0;

for i=1:n

if fitness(i,1)>=alpha*rand

nn=nn+1;

FARM(nn,:)=farm(i,:);

end

end

FARM=FARM(1:nn,:);

[aa,bb]=size(FARM);%交叉和变异

while aa if nn<=2

nnper=randperm(2);

else

nnper=randperm(nn);

end

A=FARM(nnper(1),:);

B=FARM(nnper(2),:);

很通用的遗传算法matlab程序,希望给大家有所帮助。

[A,B]=intercross(A,B);

FARM=[FARM;A;B];

[aa,bb]=size(FARM);

end

if aa>n

FARM=FARM(1:n,:);%保持种群规模为n

end

farm=FARM;

clear FARM

counter=counter+1

end

Rlength=myLength(D,R);

subplot(1,3,3)

plotaiwa(a,R)

function D=ff01(a)

[c,d]=size(a);

D=zeros(c,c);

for j=i:c

bb=(a(i,1)-a(j,1)).^2+(a(i,2)-a(j,2)).^2;

D(i,j)=bb^(0.5);

D(j,i)=D(i,j);

end

end

function [a,b]=intercross(a,b)

L=length(a);

if L<=10%确定交叉宽度

W=9;

elseif ((L/10)-floor(L/10))>=rand&&L>10

W=ceil(L/10)+8;

else

W=floor(L/10)+8;

end

for i=1:W%交叉

x=find(a==b(1,p+i-1));

y=find(b==a(1,p+i-1));

[a(1,p+i-1),b(1,p+i-1)]=exchange(a(1,p+i-1),b(1,p+i-1));

[a(1,x),b(1,y)]=exchange(a(1,x),b(1,y));

很通用的遗传算法matlab程序,希望给大家有所帮助。

end

function [x,y]=exchange(x,y)

temp=x;

x=y;

y=temp;

function plotaiwa(a,R)

scatter(a(:,1),a(:,2),'x')

hold on

plot([a(R(1),1),a(R(31),1)],[a(R(1),2),a(R(31),2)]) hold on

for i=2:length(R)

x0=a(R(i-1),1);

y0=a(R(i-1),2);

x1=a(R(i),1);

y1=a(R(i),2);

xx=[x0,x1];

yy=[y0,y1];

plot(xx,yy)

hold on

end