智能优化算法之模拟退火算法

详细的代码解释及背景，过几天再进行补充

%模拟退火算法
%问题背景是tsp问题，有100个旅游景点。并且假设距离矩阵d已经给出。
%d已知
sum=inf;
for j=1:1000
	S=[1,1+randperm(100),102];
       temp=0;
      for i=1:101;
          temp=temp+d(S(i),S(i+1));
      end
      if sum<temp
          S0=S;
          sum=temp;
      end
end;
   %%%生成一千组的随机初始解，并从中挑去最小的一组随机解作为初始解。
   e=0.1^30;
   L=20000;
   at=0.999;
   T=1;
   for k=1:L
   c=2+floor(100*rand(1,2));
   c=sort(c);
   c1=c(1);
   c2=c(2);
   df=d(S0(c1-1),S0(c2))+d(S0(c1),S0(c2+1))-d(S0(c1-1),S0(c1))-d(S0(c2),S0(C2+1));
   if df<0
       S0=[S0(1:c1-1),S0(c2:-1:c1),S0(c2+1:102)];
       sum=sum+df;
   elseif rand(1)<exp(-df/T)
           
  S0=[S0(1:c1-1),S0(c2:-1:c1),S0(c2+1:102)];
       sum=sum+df;
   end
   T=T*at;
   if T<e
       break;
   end
   end