关于MATLAB用遗传算法解决TSP问题昨天的详细描述
<P>那位高手能帮小弟解决下 如何解决TSP问题最后的输出一个坐标的直观图吗<BR>前面的程序已经调试好了，， 我想在最后的输出中输出一个坐标图来描述该问题,输出的坐标图的横坐标和纵坐标分别描述的是某个城市点的具体位子，我这里包含了10个城市。。比如一个城市A的位子在坐标上的表示是（Xa，Ya）。然后依次罗列出10个城市对应的坐标上的位子```程序的输出结果是把遗传算法反复迭代后的最优的路径用线连接起来描述一个路径```<BR>对了` 我这里还有个相应的距离距阵，根据这个距离句阵来绘制坐标系```<BR>我的距离距阵是这样的：<BR>diskTSP.txt:<BR>0 0.3361 0.3141 0.3601 0.5111 0.5176 0.2982 0.4564 0.3289 0.2821<BR>0.3361 0 0.1107 0.6149 0.8407 0.8083 0.5815 0.6418 0.4378 0.3897<BR>0.3141 0.1107 0 0.5349 0.7919 0.8207 0.5941 0.6908 0.4982 0.4464<BR>0.3601 0.6149 0.5349 0 0.3397 0.6528 0.5171 0.7375 0.6710 0.6310<BR>0.5111 0.8407 0.7919 0.3397 0 0.4579 0.4529 0.6686 0.7042 0.6862<BR>0.5176 0.8083 0.8207 0.6528 0.4579 0 0.2274 0.2937 0.4494 0.4683<BR>0.2982 0.5815 0.5941 0.5171 0.4529 0.2274 0 0.2277 0.2690 0.2694<BR>0.4564 0.6418 0.6908 0.7375 0.6686 0.2937 0.2277 0 0.2100 0.2529<BR>0.3289 0.4378 0.4982 0.6710 0.7042 0.4494 0.2690 0.2100 0 0.0531<BR>0.2821 0.3897 0.4464 0.6310 0.6862 0.4683 0.2694 0.2529 0.0531 0</P><P>谢谢了 我的邮箱是<a href="mailtlaomai@cqwu.net" target="_blank" >laomai@cqwu.net</A><BR>以下是我的程序：<BR>function =ga(d,termops,num,pc,cxops,pm,alpha)<BR>clear;<BR>load diskTSP.txt;<BR>% %————————————————————————<BR>%=ga(d,termops,num,pc,cxops,pm,alpha)<BR>%d:距离矩阵<BR>%termops:种群带数<BR>%num:每带染色体的个数<BR>%pc:交叉概率<BR>%cxops:由于本程序采用单点交叉，交叉点的设置在本程序中没有很好的解决，所以本文了采用定点，即第cxops，可以随机产生。<BR>%pm:变异概率<BR>%alpha:评价函数eval(vi)=alpha*(1-alpha).^(i-1).<BR>%bestpop:返回的最优种群<BR>%trace:进化轨迹<BR>%------------------------------------------------<BR>d=diskTSP;<BR>citynum=size(d,2);<BR>n=nargin;<BR>if n<2<BR> termops=100;<BR> num=10;<BR> pc=0.25;<BR> cxops=3;<BR> pm=0.30;<BR> alpha=0.10;<BR>end<BR>if n<3<BR> num=10;<BR> pc=0.25;<BR> cxops=3;<BR> pm=0.30;<BR> alpha=0.10;<BR>end<BR>if n<4<BR> pc=0.25;<BR> cxops=3;<BR> pm=0.30;<BR> alpha=0.10;<BR>end<BR>if n<5<BR> cxops=3;<BR> pm=0.30;<BR> alpha=0.10;<BR>end<BR>if n<6<BR> pm=0.30;<BR> alpha=0.10;<BR>end<BR>if n<7<BR> alpha=0.10;<BR>end<BR>if isempty(cxops)<BR> cxops=3;<BR>end<BR>%=========================================<BR>citynum=size(d,2);<BR>for i=1:num<BR> t(i,:)=randperm(citynum);<BR>end<BR>%=========================================<BR>for i=1:termops;<BR>%=========================================<BR>=size(t);<BR>for k=1:m<BR> for i=1:n-1<BR> l(k,i)=d(t(k,i),t(k,i+1));<BR> end<BR> l(k,n)=d(t(k,n),t(k,1));<BR> l(k)=-sum(l(k,:));<BR>end<BR>%=========================================<BR>=find(max(l));<BR>trace(i)=-l(y(1));<BR>bestpop=t(y(1),:);</P>
<P>%=========================================<BR>=size(l);<BR>t1=t;<BR>=sort(l,2);%fsort from l to u<BR>for i=1:n<BR> aftersort(i)=aftersort1(n+1-i); %change <BR>end<BR>for k=1:n;<BR> t(k,:)=t1(aftersort(k),:);<BR> l1(k)=l(aftersort(k));<BR>end<BR>t1=t;<BR>l=l1;<BR>for i=1:size(aftersort,2)<BR> evalv(i)=alpha*(1-alpha).^(i-1);<BR>end<BR>m=size(t,1);<BR>q=cumsum(evalv);<BR>qmax=max(q);<BR>for k=1:m<BR> r=qmax*rand(1);<BR> for j=1:m<BR> if j==1&r<=q(1)<BR> t(k,:)=t1(1,:);<BR> elseif j~=1&r>q(j-1)&r<=q(j)<BR> t(k,:)=t1(j,:);<BR> end<BR> end<BR>end<BR>%==========================================<BR>=size(t);<BR>for k=1:m<BR> t0=1:n;<BR> for i=1:n<BR> for j=1:length(t0)<BR> if t(k,i)==t0(j)<BR> g(k,i)=j;<BR> t0(j)=[];<BR> break<BR> end<BR> end<BR> end<BR>end<BR>%============================================<BR>=size(g);<BR>ran=rand(1,m);<BR>r=cxops;<BR>=find(ran<pc);<BR>if ru>=2<BR> for k=1:2:length(ru)-1<BR> g1(ru(k),:)=),g(ru(k+1),[(r+1):n])];<BR> g(ru(k+1),:)=),g(ru(k),[(r+1):n])];<BR> g(ru(k),:)=g1(ru(k),:);<BR> end<BR>end<BR>%===============================================<BR> %均匀变异<BR>=size(g);<BR>ran=rand(1,m);<BR>r=rand(1,3); %dai gai jin<BR>rr=floor(n*rand(1,3)+1);<BR>=find(ran<pm);<BR>for k=1:length(mu)<BR> for i=1:length(r)<BR> umax(i)=n+1-rr(i);<BR> umin(i)=1;<BR> g(mu(k),rr(i))=umin(i)+floor((umax(i)-umin(i))*r(i));<BR> end<BR>end %均匀变异<BR>%===================================================<BR>=size(g);<BR>for k=1:m<BR> t0=1:n;<BR> for i=1:n<BR> t(k,i)=t0(g(k,i));<BR> t0(g(k,i))=[];<BR> end<BR>end<BR>end</P>
<P><BR> </P>
[ 本帖最后由 cdwxg 于 2006-8-8 19:41 编辑 ]
关于MATLAB用遗传算法解决TSP问题昨天的详细描述
原帖由 laomaizpy 于 2006-4-24 13:14 发表<P>那位高手能帮小弟解决下 如何解决TSP问题最后的输出一个坐标的直观图吗<BR>前面的程序已经调试好了，， 我想在最后的输出中输出一个坐标图来描述该问题,输出的坐标图的横坐标和纵坐标分别描述的是 ...
好象TSP用蚁群算法好吧?
[ 本帖最后由 cdwxg 于 2006-8-8 19:41 编辑 ]
页:
[1]