最佳哈密尔顿圈MATLAB程序8页

1、实用文案M文件function circle,long=modifycircle(c1,L); global a flag=1; while flag>0 flag=0; for m=1:L-3 for n=m+2:L-1 if a(c1(m),c1(n)+a(c1(m+1),c1(n+1)<. a(c1(m),c1(m+1)+a(c1(n),c1(n+1) flag=1; c1(m+1:n)=c1(n:-1:m+1); end end end end long=a(c1(1),c1(L); for i=1:L-1 long=long+a(c1(i),c1(i+1); end ci

2、rcle=c1; 命令窗口全程的：global aL=size(a,1);c1=1 2:18 20:53 19;circle,long=modifycircle(c1,L);c2=1 19 2:18 20:53;%改变初始圈，该算法的最后一个顶点不动circle2,long2=modifycircle(c2,L);if long2<longlong=long2;circle=circle2;endcircle,longcircle = Columns 1 through 19 1 16 44 45 46 47 17 48 50 49 51 53 52 3 4 21 5 22 26 Col

3、umns 20 through 38 6 25 13 37 38 39 36 10 11 31 8 29 27 28 7 30 9 32 33 Columns 39 through 53 34 35 12 40 41 42 15 43 14 24 23 20 18 2 19long =5.866999999999999e+002 1a=0.00 11.50 11.90 10.10 20.60 28.40 22.20 28.00 20.80 12.90 16.30 6.00 27.80 36.00 23.70 22.10 30.20 35.70 ;11.50 0.00 11.00 21.60 3

4、2.10 39.90 33.70 39.50 32.30 24.40 21.50 11.20 28.60 36.80 28.90 33.60 41.70 47.20 ;11.90 11.00 0.00 22.00 32.50 40.30 34.10 36.10 28.90 21.00 12.20 5.90 17.60 25.80 19.60 26.90 35.00 43.80 ;10.10 21.60 22.00 0.00 10.50 18.30 12.10 20.40 15.20 23.00 26.40 16.10 37.90 46.10 33.80 32.20 38.40 28.10 ;2

5、0.60 32.10 32.50 10.50 0.00 7.80 15.70 24.00 25.70 33.50 36.90 26.60 48.40 56.60 44.30 42.70 42.00 31.70 ;28.40 39.90 40.30 18.30 7.80 0.00 7.90 16.20 23.40 31.30 40.10 34.40 51.60 49.10 47.50 40.50 34.20 23.90 ;22.20 33.70 34.10 12.10 15.70 7.90 0.00 8.30 15.50 23.40 32.20 28.20 43.70 41.20 39.60 3

6、2.60 26.30 16.00 ;28.00 39.50 36.10 20.40 24.00 16.20 8.30 0.00 7.20 15.10 23.90 34.00 35.40 32.90 31.30 24.30 18.00 7.70 ;20.80 32.30 28.90 15.20 25.70 23.40 15.50 7.20 0.00 7.90 16.70 26.80 28.20 36.40 24.10 17.10 25.20 14.90 ;12.90 24.40 21.00 23.00 33.50 31.30 23.40 15.10 7.90 0.00 8.80 18.90 20

7、.30 28.50 16.20 9.20 17.30 22.80 ;16.30 21.50 12.20 26.40 36.90 40.10 32.20 23.90 16.70 8.80 0.00 10.30 11.50 19.70 7.40 14.70 22.80 31.60 ;6.00 11.20 5.90 16.10 26.60 34.40 28.20 34.00 26.80 18.90 10.30 0.00 21.80 30.00 17.70 25.00 33.10 41.70 ;27.80 28.60 17.60 37.90 48.40 51.60 43.70 35.40 28.20

8、20.30 11.50 21.80 0.00 8.20 18.90 26.20 23.10 33.40 ;36.00 36.80 25.80 46.10 56.60 49.10 41.20 32.90 36.40 28.50 19.70 30.00 8.20 0.00 20.30 23.00 14.90 25.20 ;23.70 28.90 19.60 33.80 44.30 47.50 39.60 31.30 24.10 16.20 7.40 17.70 18.90 20.30 0.00 7.30 15.40 25.70 ;22.10 33.60 26.90 32.20 42.70 40.5

9、0 32.60 24.30 17.10 9.20 14.70 25.00 26.20 23.00 7.30 0.00 8.10 18.40 ;30.20 41.70 35.00 38.40 42.00 34.20 26.30 18.00 25.20 17.30 22.80 33.10 23.10 14.90 15.40 8.10 0.00 10.30 ;35.70 47.20 43.80 28.10 31.70 23.90 16.00 7.70 14.90 22.80 31.60 41.70 33.40 25.20 25.70 18.40 10.30 0.00 ;global aL=size(

10、a,1);c1=1 2:11 13:18 12;circle,long=modifycircle(c1,L);c2=1 12 2:11 13:18;%改变初始圈，该算法的最后一个顶点不动circle2,long2=modifycircle(c2,L);if long2<longlong=long2;circle=circle2;endcircle,long结果circle = 1 4 5 6 7 8 9 10 11 15 16 17 18 14 13 3 2 12long = 1.808000000000000e+002global aL=size(a,1);c1=1 14 13 3 2

11、 4:11 15:18 12;circle,long=modifycircle(c1,L);c2=1 14 13 3 2 12 4:11 15:18;%改变初始圈，该算法的最后一个顶点不动circle2,long2=modifycircle(c2,L);if long2<longlong=long2;circle=circle2;endcircle,long结果circle = 1 4 5 6 7 18 8 9 10 11 15 16 17 14 13 3 2 12long =1.756000000000000e+0022a=0.00 19.80 43.80 47.60 55.70 69

12、.90 61.10 52.90 38.30 31.80 39.60 43.70 53.50 60.30 51.90 41.90 34.00 47.20 ;19.80 0.00 24.00 27.80 35.90 50.10 41.30 33.10 18.50 12.00 19.80 23.90 33.70 40.50 32.10 22.10 14.20 27.40 ;43.80 24.00 0.00 7.20 15.30 39.90 31.10 23.50 5.50 12.00 13.40 17.50 27.30 34.10 27.60 22.50 20.80 31.40 ;47.60 27.

13、80 7.20 0.00 8.10 32.70 23.90 16.30 9.30 15.80 17.20 21.30 31.10 35.70 31.40 26.30 24.60 35.20 ;55.70 35.90 15.30 8.10 0.00 24.60 15.80 8.20 17.40 23.90 21.50 17.40 27.20 27.60 27.50 30.60 32.70 39.50 ;69.90 50.10 39.90 32.70 24.60 0.00 8.80 17.00 38.20 38.10 30.30 26.20 27.40 20.60 34.10 39.40 46.9

14、0 48.30 ;61.10 41.30 31.10 23.90 15.80 8.80 0.00 8.20 29.40 29.30 21.50 17.40 18.60 11.80 25.30 30.60 38.10 39.50 ;52.90 33.10 23.50 16.30 8.20 17.00 8.20 0.00 21.20 21.10 13.30 9.20 19.00 20.00 19.30 22.40 29.90 31.30 ;38.30 18.50 5.50 9.30 17.40 38.20 29.40 21.20 0.00 6.50 7.90 12.00 21.80 28.60 2

15、2.10 17.00 15.30 25.90 ;31.80 12.00 12.00 15.80 23.90 38.10 29.30 21.10 6.50 0.00 7.80 11.90 21.70 28.50 22.00 16.70 8.80 22.00 ;39.60 19.80 13.40 17.20 21.50 30.30 21.50 13.30 7.90 7.80 0.00 4.10 13.90 20.70 14.20 9.10 16.60 18.00 ;43.70 23.90 17.50 21.30 17.40 26.20 17.40 9.20 12.00 11.90 4.10 0.0

16、0 9.80 16.60 10.10 13.20 20.70 22.10 ;53.50 33.70 27.30 31.10 27.20 27.40 18.60 19.00 21.80 21.70 13.90 9.80 0.00 6.80 6.70 16.70 24.60 25.60 ;60.30 40.50 34.10 35.70 27.60 20.60 11.80 20.00 28.60 28.50 20.70 16.60 6.80 0.00 13.50 23.50 31.40 32.40 ;51.90 32.10 27.60 31.40 27.50 34.10 25.30 19.30 22

17、.10 22.00 14.20 10.10 6.70 13.50 0.00 10.00 17.90 18.90 ;41.90 22.10 22.50 26.30 30.60 39.40 30.60 22.40 17.00 16.70 9.10 13.20 16.70 23.50 10.00 0.00 7.90 8.90 ;34.00 14.20 20.80 24.60 32.70 46.90 38.10 29.90 15.30 8.80 16.60 20.70 24.60 31.40 17.90 7.90 0.00 13.20 ;47.20 27.40 31.40 35.20 39.50 48

18、.30 39.50 31.30 25.90 22.00 18.00 22.10 25.60 32.40 18.90 8.90 13.20 0.00 ;global aL=size(a,1);c1=1 3:18 2;circle,long=modifycircle(c1,L);c2=1 2 3:18;%改变初始圈，该算法的最后一个顶点不动circle2,long2=modifycircle(c2,L);if long2<longlong=long2;circle=circle2;endcircle,longcircle = 1 10 9 3 4 5 6 7 8 11 12 13 14 15

19、 16 18 17 2long = 2.143000000000000e+002global aL=size(a,1);c1=1 10 9 11:16 18 17 3:8 2;circle,long=modifycircle(c1,L);c2=1 10 9 11:16 18 17 2 3:8;%改变初始圈，该算法的最后一个顶点不动circle2,long2=modifycircle(c2,L);if long2<longlong=long2;circle=circle2;endcircle,longcircle = 1 17 18 16 15 13 14 6 7 8 5 4 3 9 12

20、 11 10 2long = 2.019000000000000e+0023a=0.00 9.20 14.00 34.90 17.50 27.20 34.50 49.70 49.50 55.90 67.30 71.30 79.90 22.20 41.70 55.10 62.70 77.50 65.90 ;9.20 0.00 4.80 25.70 8.30 18.00 25.30 40.50 40.30 46.70 58.10 62.10 70.70 13.00 32.50 45.90 53.50 68.30 56.70 ;14.00 4.80 0.00 20.90 13.10 22.80 23

21、.30 38.50 38.30 44.70 56.10 60.10 68.70 8.20 30.50 43.90 51.50 66.30 54.70 ;34.90 25.70 20.90 0.00 24.00 33.70 27.80 20.40 36.20 42.60 54.00 58.00 66.60 12.70 28.40 41.80 49.40 64.20 52.60 ;17.50 8.30 13.10 24.00 0.00 9.70 17.00 32.20 32.00 38.40 49.80 53.80 62.40 11.30 24.20 37.60 45.20 60.00 48.40

22、 ;27.20 18.00 22.80 33.70 9.70 0.00 7.30 22.50 22.30 28.70 40.10 44.10 52.70 21.00 14.50 27.90 35.50 50.30 38.70 ;34.50 25.30 23.30 27.80 17.00 7.30 0.00 15.20 15.00 21.40 32.80 36.80 45.40 15.10 7.20 20.60 28.20 43.00 31.40 ;49.70 40.50 38.50 20.40 32.20 22.50 15.20 0.00 15.80 22.20 33.60 37.60 46.

23、20 30.30 8.00 21.40 29.00 43.80 32.20 ;49.50 40.30 38.30 36.20 32.00 22.30 15.00 15.80 0.00 22.00 17.80 34.20 37.90 30.10 7.80 5.60 25.60 28.00 16.40 ;55.90 46.70 44.70 42.60 38.40 28.70 21.40 22.20 22.00 0.00 14.60 15.40 24.00 36.50 14.20 26.80 6.80 24.80 37.60 ;67.30 58.10 56.10 54.00 49.80 40.10

24、32.80 33.60 17.80 14.60 0.00 16.40 20.10 47.90 25.60 12.20 7.80 10.20 23.00 ;71.30 62.10 60.10 58.00 53.80 44.10 36.80 37.60 34.20 15.40 16.40 0.00 8.60 51.90 29.60 28.60 8.60 18.50 39.40 ;79.90 70.70 68.70 66.60 62.40 52.70 45.40 46.20 37.90 24.00 20.10 8.60 0.00 60.50 38.20 32.30 17.20 9.90 43.10

25、;22.20 13.00 8.20 12.70 11.30 21.00 15.10 30.30 30.10 36.50 47.90 51.90 60.50 0.00 22.30 35.70 43.30 58.10 46.50 ;41.70 32.50 30.50 28.40 24.20 14.50 7.20 8.00 7.80 14.20 25.60 29.60 38.20 22.30 0.00 13.40 21.00 35.80 24.20 ;55.10 45.90 43.90 41.80 37.60 27.90 20.60 21.40 5.60 26.80 12.20 28.60 32.3

26、0 35.70 13.40 0.00 20.00 22.40 10.80 ;62.70 53.50 51.50 49.40 45.20 35.50 28.20 29.00 25.60 6.80 7.80 8.60 17.20 43.30 21.00 20.00 0.00 18.00 30.80 ;77.50 68.30 66.30 64.20 60.00 50.30 43.00 43.80 28.00 24.80 10.20 18.50 9.90 58.10 35.80 22.40 18.00 0.00 33.20 ;65.90 56.70 54.70 52.60 48.40 38.70 31.40 32.20 16.40 37.60 23.00 39.40 43.10 46.50 24.20 10.80 30.80 33.20 0.00 ;global aL=size(a,1);c1=1 3:19 2;circle,long=modifycircle(c1,L);c2=1 2 3:19;%改变初始圈，该算法的最后一个顶点不动circle2,long2=modifycircle(c2,L);if long2<longlong=long2;circle=circle2;endcircle,longcircle = Columns 1 through 1