1、完全图哈密尔顿圈的遗传模拟退火算法 matlab 通用源程序% maxpop 给定群体规模% pop 群体% newpop 种群%t0 初始温度function codmin,finmin=fc0(cc,v0,t0)N=length(cc(1,:);%定群体规模if N50maxpop=2*N-20;endif N0%用模拟退火产生新的群体 pop=fc1(maxpop,pop,N,cc,v0,t);%转轮赌选择种群 f=zeros(1,maxpop);for i=1:maxpopfor j=1:N-1x=pop(i,j);y=pop(i,j+1);fo1=cc(pop(i,j),pop(i,
2、j+1);f(i)=f(i)+fo1;endf(i)=f(i)+cc(pop(i,1),pop(i,N);endfmin=min(f);for i=1:maxpopif fmin=infendif fmin=inf|f(i)=infdd=fmin-f(i);endftk(i)=exp(dd/t);endfin1,cod=sort(-ftk);fin=abs(fin1);%f(cod(1)if f(cod(1)=RR);% codnewpop(i,:)=pop(cod(cod2(end),:);end%单亲繁殖if N32jmax=round(N/9);endif NR1for j=1:2:jm
3、ax+2nn=randperm(N);x=nn(j);y=nn(j+1);if newpop(i,x)=v0|newpop(i,y)=v0pop(i,:)=newpop(i,:);continue;endbox1=newpop(i,x);newpop(i,x)=newpop(i,y);newpop(i,y)=box1;pop(i,:)=newpop(i,:); endendend%温度下降t=t-0.1;endfunction pop=fc1(maxpop,pop,N,cc,v0,t) ff(N-1)=0;f=0;pop1=zeros(maxpop,N);for i=1:maxpopfor j
4、=1:N-1x=pop(i,j);y=pop(i,j+1);ff(j)=cc(pop(i,j),pop(i,j+1);pop1(i,:)=pop(i,:);nn=randperm(N);x=nn(1);y=nn(2);pop1=pop;if pop(i,x)=v0|pop(i,x)=v0continuebox1=pop(i,x);pop1(i,x)=pop1(i,y);pop1(i,y)=box1;endff1(j)=cc(pop1(i,j),pop1(i,j+1);endf=sum(ff);f1=sum(ff1);if f=infendif f=inf|f1=infdd=f-f1;endAij=min(1,exp(dd/t);Pacept=rand(1);if AijPaceptpop(i,:)=pop1(i,:);endend
