原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
最佳哈密尔顿圈MATLAB程序.doc
M文件
function [circle,long]=modifycircle(c1,L);
global a
flag=1;
while flag0
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
circle=c1;
命令窗口
全程的：
global a
L=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 long2long
long=long2;
circle=circle2;
end
circle,long
circle =
Columns 1 through 19
1 16 44 45 46 47 17 48 50 49 51 53 52 3 4 21 5 22 26
Columns 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 19
long =
5.866999999999999e+002
1
a=[
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 32.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 ;
20.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 32.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.2
文档评论（0）