大家好,我是小鸭酱,博客地址为:http://www.cnblogs.com/xiaoyajiang
%解决完备图中的较优H圈
clc
clear
w = [ inf 6 1 8 3 1 ;... %权赋矩阵
6 inf 1 1 6 9 ;...
1 1 inf 3 6 7 ;...
8 1 3 inf 1 2 ;...
3 6 6 1 inf 1 ;...
1 9 7 2 1 inf ] ;
v = [ 11 12 13 14 15 16 ;... %两点之间的路径,如25表示从点2到点5
21 22 23 24 25 26 ;...
31 32 33 34 35 36 ;...
41 42 43 44 45 46 ;...
51 52 53 54 55 56 ;...
61 62 63 64 65 66 ] ;
p = [ 12 23 34 45 56 61 ] %保存某个H圈的每两点间的路径,H圈便是由这些小的路径构成的
W = zeros(1,6) ;
S=0;
for h = 1 : 6
d= rem(p(h),10) ;
c= (p(h) - d)/10 ;
W(h) = w(c,d) ;
S= S + W(h) ;
end
S
for i=1:1000
x = randperm(6) ;
m = x(1) ;
n = x(2) ;
j = rem(p(m),10) ;
i = (p(m) - j)/10 ;
z = rem(p(n),10) ;
y = (p(n) - z)/10 ;
if ( w(i,j) + w(y,z) > w(i,y) + w(j,z) )
s = 0;
for b = 1 : 6
if ( p(b) == v(i,y) | p(b) == v(j,z) | p(b) == v(y,i)| p(b) == v(z,j) )
s = 1 ;
break ;
end
end
if s
i = i + 1 ;
break ;
end
p(m) = v(i,y) ;
p(n) = v(j,z) ;
end
end
W = zeros(1,6) ;
S=0;
for h = 1 : 6
d= rem(p(h),10) ;
c= (p(h) - d)/10 ;
W(h) = w(c,d) ;
S= S + W(h) ;
end
p
S