@热带火原子


GA无误

以前写的,用Logistic混沌序列做了改进,讲究看了吧

好复杂，谢谢指教，我先学习一下，到时再求教一下

什么乱码？删除它，吧务～

什么乱码？删除它，吧务～

回复8楼:不可以呀

不得不说这货编程秒我n条街

回复9楼:这样会突显我的无知的。。

貌似...错过了 抱歉 这几天忙爆了今天才有空上贴吧
如果以后你还想了解我期末完了以后详细给你讲讲...我其实更喜欢BP神经网络结合遗传算法的来解决NP-hard的~~

回复12楼:已经非常感谢了，我以你给的算法为原型成功完成了作业

其实蚁群算法更适合TSP的.....PSO最简单但太早熟,标准GA只收敛到全局最优解的90%左右(所以我是用改良圈初始化染色体,不然最后结果很糟糕....不过现在我都用量子旋转门了嘿嘿..),ACA则最适合组合优化~酱紫~
@丶火贼丨浮云丶


function ACAforTSP(city)
n=size(city,1);
d=zeros(n,n);
[ii jj]=meshgrid(1:n);
d=arrayfun(@(i,j)sqrt(sum((city(i,:)-city(j,:)).^2)),ii,jj);
d=d+eye(n)*1e-4;
m=35;
alpha=1;
beta=5;
rho=.1;
q=1;
eta=1./d;
tau=ones(n,n);
table=zeros(m,n);
iter=1;
iter_max=200;
route_best=zeros(iter_max,n);
length_best=zeros(iter_max,1);
while iter<=iter_max
table(:,1)=unidrnd(n,m,1);
city_index=1:n;
for i=1:m
for j=2:n
tabu=table(i,1:(j-1));
allow=setdiff(1:n,tabu);
p=zeros(size(allow));
p=arrayfun(@(k)tau(tabu(end),allow(k))^alpha*eta(tabu(end),allow(k))^beta,1:length(allow));
p=p/sum(p);
table(i,j)=allow(find(cumsum(p)>=rand,1,'first'));
end
end
l=zeros(m,1);
for i=1:m
route=table(i,:);
for j=1:(n-1)
l(i)=l(i)+d(route(j),route(j+1));
end
l(i)=l(i)+d(route(n),route(1));
end
if iter==1
[min_length,min_index]=min(l);
length_best(iter)=min_length;
route_best(iter,:)=table(min_index,:);
else
[min_length,min_index]=min(l);
length_best(iter)=min(length_best(iter-1),min_length);
if length_best(iter)==min_length
route_best(iter,:)=table(min_index,:);
else
route_best(iter,:)=route_best(iter-1,:);
end
end
delta_tau=zeros(n,n);
for i=1:m
for j=1:(n-1)
delta_tau(table(i,j),table(i,j+1))=delta_tau(table(i,j),table(i,j+1))+q/length(i);
end
delta_tau(table(i,n),table(i,1))=delta_tau(table(i,n),table(i,1))+q/length(i);
end
tau=(1-rho)*tau+delta_tau;
iter=iter+1;
table=zeros(m,n);
end
[shortest_length index]=min(length_best);
shortest_route=route_best(index,:);
end
不送

fit=@(x,y)x.*sin(4*pi*x)+y.*sin(20*pi*y);

最优解:11.631 5.72186
最大值:17.2093