旅行商模型利用LINGO规划算法求最佳哈密顿圈的算法

lingo部分
model:
sets:
cities/1..19/:u;
link(cities,cities):
distance,x;endsets
data:
distance=
    ;在此输入matlab计算的边权值矩阵
enddata
n=@size(cities);
min=@sum(link(i,j)|i#ne#j:distance(i,j)*x(i,j));
@for(cities(i):
@sum(cities(j)|j#ne#i:x(j,i))=1;
@sum(cities(j)|j#ne#i:x(i,j))=1;
@for(cities(j)|j#gt#1#and#j#ne#i:
u(j)>=u(i)+x(i,j)
-(n-2)*(1-x(i,j))+(n-3)*x(j,i);
);
);
@for(link:@bin(x));
@for(cities(i)|i#gt#1:
u(i)<=n-1-(n-2)*x(1,i);
u(i)>=1+(n-2)*x(i,1););
end

matlab对图片的复原代码：
circle=%输入lingo计算出的排序
sum2;
for k=1:19
  
    circle2(k)=circle(20-k);
end
circle2

for i=1:19
begin=72*(i-1)+1;
end_f=72*i;
num=circle2(i);
all_file_new(:,begin:end_f)=I{num};
end
figure
show_file_new=all_file_new(1:1980,:);
imshow(show_file_new);