matlab 模拟退火求解tsp.txt

MATLAB绘制出权值wv和阀值bv确定的误差曲面

模拟退火求解tsp2008-01-21 10:05d = Import["E:\\nb\\模拟退火.xls"][[1]] // Floor(*来自w2008z空间*);

exch[pi0_, d_] := Module[{n, u, v, pi, k, i, pi1, df},
n = Length[d]; u = Random[Integer, {2, n - 2}];
v = Random[Integer, {1, n - 1}];
v = u + v;
If[v > n, v = n
   ](*if*);

pi1 = pi0;
{pi1[[u]], pi1[[v]]} = {pi0[[v]], pi0[[u]]};

Do[pi1[[u + k]] = pi0[[v - k]], {k, 1, v - u - 1}];

If[v < n,
   df = Sum[
     d[[pi1[[i - 1]], pi1[[i]]]] - d[[pi0[[i - 1]], pi0[[i]]]], {i, u, 
      v + 1}],
   df = Sum[
      d[[pi1[[i - 1]], pi1[[i]]]] - d[[pi0[[i - 1]], pi0[[i]]]], {i, u, n}] + 
     d[[pi1[[n]], pi1[[1]]]] - d[[pi0[[n]], pi0[[1]]]]
   ](*if*);
{df, pi1}
]

trp[d_, t0_, tf_] := 
Module[{n, Lk, t, pi0, minf, pmin, flag, k, ftemp, i, df, pi1}(*来自w2008z空间*),
n = Length[d]; Lk = 100*n;
t = t0; pi0 = Range[n];
minf = Sum[d[[pi0[[i]], pi0[[i + 1]]]], {i, n - 1}] + 
    d[[pi0[[n]], pi0[[1]]]];
pmin = pi0; flag = 1;

While[t > tf && flag <= 1000,
   flag = flag + 1;
   For[k = 1, k <= Lk, k++(*来自w2008z空间*),
    {df, pi1} = exch[pi0, d];
    If[df < 0 || E^(df/t) > Random[],
     pi0 = pi1;
     ](*if*);
    ](*for*);
   ftemp = 
    Sum[d[[pi0[[i]], pi0[[i + 1]]]], {i, n - 1}] + d[[pi0[[n]], pi0[[1]]]];
   If[minf > ftemp,
    {minf, pmin} = {ftemp, pi0}(*来自w2008z空间*);
    ](*if*);
   t = 0.87*t
   ](*While*);
{minf, pmin}
]

Table[trp[d, 10, 0.5], {i, 5}] // Union

{{43, {1, 3, 2, 5, 4}}, {43, {1, 5, 4, 2, 3}}}