灰狼优化算法求旅行商问题Matlab代码

% D-GWO（Discrete Grey Wolf Optimizer for symmetric travelling salesman problem） clc; clear all; close all; tic; %% 输入数据 %dangtzig42,最优解699，已求到675 % x=[170,166,133,140,142,126,125,119,117,99,73,72,37,6,3,21,33,4,3,27,52,57,58,88,99,95,89,83,85,104,112,112,113,125,135,147,147.5,154.5,157,158,172,174]; % y=[85,88,73,70,55,53,60,68,74,83,79,91,94,106,97,82,67,66,42,33,41,59,66,65,67,55,55,38,25,35,37,24,13,30,32,18,36,45,54,61,82,87]; %att48,最优解33522，已求到33522 % x=[6734,2233,5530,401,3082,7608,7573,7265,6898,1112,5468,5989,4706,4612,6347,6107,7611,7462,7732,5900,4483,6101,5199,1633,4307,675,7555,7541,3177,7352,7545,3245,6426,4608,23,7248,7762,7392,3484,6271,4985,1916,7280,7509,10,6807,5185,3023]; % y=[1453,10,1424,841,1644,4458,3716,1268,1885,2049,2606,2873,2674,2035,2683,669,5184,3590,4723,3561,3369,1110,2182,2809,2322,1006,4819,3981,756,4506,2801,3305,3173,1198,2216,3779,4595,2244,2829,2135,140,1569,4899,3239,2676,2993,3258,1942]; %eil51，最优解426 % x=[37.0, 49.0, 52.0, 20.0, 40.0, 21.0, 17.0, 31.0, 52.0, 51.0, 42.0, 31.0, 5.0, 12.0, 36.0, 52.0, 27.0, 17.0, 13.0, 57.0, 62.0, 42.0, 16.0, 8.0, 7.0, 27.0, 30.0, 43.0, 58.0, 58.0, 37.0, 38.0, 46.0, 61.0, 62.0, 63.0, 32.0, 45.0, 59.0, 5.0, 10.0, 21.0, 5.0, 30.0, 39.0, 32.0, 25.0, 25.0, 48.0, 56.0, 30.0]; % y=[52.0, 49.0, 64.0, 26.0, 30.0, 47.0, 63.0, 62.0, 33.0, 21.0, 41.0, 32.0, 25.0, 42.0, 16.0, 41.0, 23.0, 33.0, 13.0, 58.0, 42.0, 57.0, 57.0, 52.0, 38.0, 68.0, 48.0, 67.0, 48.0, 27.0, 69.0, 46.0, 10.0, 33.0, 63.0, 69.0, 22.0, 35.0, 15.0, 6.0, 17.0, 10.0, 64.0, 15.0, 10.0, 39.0, 32.0, 55.0, 28.0, 37.0, 40.0]; %st70,最优解675，已求到675 x=[64.0, 80.0, 69.0, 72.0, 48.0, 58.0, 81.0 , 79.0, 30.0, 42.0, 7.0, 29.0, 78.0, 64.0, 95.0, 57.0, 40.0, 68.0, 92.0, 62.0, 28.0, 76.0, 67.0, 93.0, 6.0, 87.0, 30.0, 77.0, 78.0, 55.0, 82.0, 73.0, 20.0, 27.0, 95.0, 67.0, 48.0, 75.0, 8.0, 20.0, 54.0, 63.0, 44.0, 52.0, 12.0, 25.0, 58.0, 5.0, 90.0, 41.0, 25.0, 37.0, 56.0, 10.0, 98.0, 16.0, 89.0, 48.0, 81.0, 29.0, 17.0, 5.0, 79.0, 9.0, 17.0, 74.0, 10.0, 48.0, 83.0, 84.0]; y=[96.0, 39.0, 23.0, 42.0, 67.0, 43.0, 34.0, 17.0, 23.0, 67.0, 76.0, 51.0, 92.0, 8.0, 57.0, 91.0, 35.0, 40.0, 34.0, 1.0, 43.0, 73.0, 88.0, 54.0, 8.0, 18.0, 9.0, 13.0, 94.0, 3.0, 88.0, 28.0, 55.0, 43.0, 86.0, 99.0, 83.0, 81.0, 19.0, 18.0, 38.0, 36.0, 33.0, 18.0, 13.0, 5.0, 85.0, 67.0, 9.0, 76.0, 76.0, 64.0, 63.0, 55.0, 7.0, 74.0, 60.0, 82.0, 76.0, 60.0, 22.0, 45.0, 70.0, 100.0, 82.0, 67.0, 68.0, 19.0, 86.0, 94.0]; %pr76,最优解108159，已求到108159 % x=[3600,3100,4700,5400,5608,4493,3600,3100,4700,5400,5610,4492,3600,3100,4700,5400,6650,7300,7300,6650,7300,6650,5400,8350,7850,9450,10150,10358,9243,8350,7850,9450,10150,10360,9242,8350,7850,9450,10150,11400,12050,12050,11400,12050,11400,10150,13100,12600,14200,14900,15108,13993,13100,12600,14200,14900,15110,13992,13100,12600,14200,14900,16150,16800,16800,16150,16800,16150,14900,19800,19800,19800,19800,200,200,200]; % y=[2300,3300,5750,5750,7103,7102,6950,7250,8450,8450,10053,10052,10800,10950,11650,11650,10800,10950,7250,6950,3300,2300,1600,2300,3300,5750,5750,7103,7102,6950,7250,8450,8450,10053,10052,10800,10950,11650,11650,10800,10950,7250,6950,3300,2300,1600,2300,3300,5750,5750,7103,7102,6950,7250,8450,8450,10053,10052,10800,10950,11650,11650,10800,10950,7250,6950,3300,2300,1600,800,10000,11900,12200,12200,1100,800]; %rat99,%最优解1211，已求到1212 % x=[6,15,24,33,48,57,67,77,86,6,17,23,32,43,55,65,78,87,3,12,28,33,47,55,64,71,87,4,15,22,34,42,54,66,78,87,7,17,26,32,43,57,64,78,83,5,13,25,38,46,58,67,74,88,2,17,23,36,42,53,63,72,87,2,16,25,38,42,57,66,73,86,5,13,25,35,46,54,65,73,86,2,14,28,38,46,57,63,77,85,8,12,22,34,47,58,66,78,85]; % y=[4,15,18,12,12,14,10,10,15,21,26,25,35,23,35,36,39,35,53,44,53,49,46,52,50,57,57,72,78,70,71,79,77,79,67,73,81,95,98,97,88,89,85,83,98,109,111,102,119,107,110,110,113,110,124,134,129,131,137,123,135,134,129,146,147,153,155,158,154,151,151,149,177,162,169,177,172,166,174,161,162,195,196,189,187,195,194,188,193,194,211,217,210,216,203,213,206,210,204]; %kroA100,最优解21282，已求到21282 % x=[1380,2848,3510,457,3888,984,2721,1286,2716,738,1251,2728,3815,3683,1247,123,1234,252,611,2576,928,53,1807,274,2574,178,2678,1795,3384,3520,1256,1424,3913,3085,2573,463,3875,298,3479,2542,3955,1323,3447,2936,1621,3373,1393,3874,938,3022,2482,3854,376,2519,2945,953,2628,2097,890,2139,2421,2290,1115,2588,327,241,1917,2991,2573,19,3911,872,2863,929,839,3893,2178,3822,378,1178,2599,3416,2961,611,3113,2597,2586,161,1429,742,1625,1187,1787,22,3640,3756,776,1724,198,3950]; % y=[939,96,1671,334,666,965,1482,525,1432,1325,1832,1698,169,1533,1945,862,1946,1240,673,1676,1700,857,1711,1420,946,24,1825,962,1498,1079,61,1728,192,1528,1969,1670,598,1513,821,236,1743,280,1830,337,1830,1646,1368,1318,955,474,1183,923,825,135,1622,268,1479,981,1846,1806,1007,1810,1052,302,265,341,687,792,599,674,1673,1559,558,1766,620,102,1619,899,1048,100,901,143,1605,1384,885,1830,1286,906,134,1025,1651,706,1009,987,43,882,392,1642,1810,1558]; %kroB100,最优解22141,已求到22178.5727 % x=[3140,556,3675,1182,3595,962,2030,3507,2642,3438,3858,2937,376,839,706,749,298,694,387,2801,3133,1517,1538,844,2639,3123,2489,3834,3417,2938,71,3245,731,2312,2426,380,2310,2830,3829,3684,171,627,1490,61,422,2698,2372,177,3084,1213,3,1782,3896,1829,1286,3017,2132,2000,3317,1729,2408,3292,193,782,2503,1697,3821,3370,3162,3938,2741,2330,3918,1794,2929,3453,896,399,2614,2800,2630,563,1090,2009,3876,3084,1526,1612,1423,3058,3782,347,3904,2191,3220,468,3611,3114,3515,3060]; % y=[1401,1056,1522,1853,111,1895,1186,1851,1269,901,1472,1568,1018,1355,1925,920,615,552,190,695,1143,266,224,520,1239,217,1520,1827,1808,543,1323,1828,1741,1270,1851,478,635,775,513,445,514,1261,1123,81,542,1221,127,1390,748,910,1817,995,742,812,550,108,1432,1110,1966,1498,1747,152,1210,1462,352,1924,147,791,367,516,1583,741,1088,1589,485,1998,705,850,195,653,20,1513,1652,1163,1165,774,1612,328,1322,1276,1865,252,1444,1579,1454,319,1968,1629,1892,155]; %kroC100,最优解20749，已求到20750.7625 % x=[1357,2650,1774,1307,3806,2687,43,3092,185,834,40,1183,2048,1097,1838,234,3314,737,779,2312,2576,3078,2781,705,3409,323,1660,3729,693,2361,2433,554,913,3586,2636,1000,482,3704,3635,1362,2049,2552,3939,219,812,901,2513,242,826,3278,86,14,1327,2773,2469,3835,1031,3853,1868,1544,457,3174,192,2318,2232,396,2365,2499,1410,2990,3646,3394,1779,1058,2933,3099,2178,138,2082,2302,805,22,3213,99,1533,3564,29,3808,2221,3499,3124,781,1027,3249,3297,213,721,3736,868,960]; % y=[1905,802,107,964,746,1353,1957,1668,1542,629,462,1391,1628,643,1732,1118,1881,1285,777,1949,189,1541,478,1812,1917,1714,1556,1188,1383,640,1538,1825,317,1909,727,457,1337,1082,1174,1526,417,1909,640,898,351,1552,1572,584,1226,799,1065,454,1893,1286,1838,963,428,1712,197,863,1607,1064,1004,1925,1374,828,1649,658,307,214,1018,1028,90,372,1459,173,978,1610,1753,1127,272,1617,1085,536,1780,676,6,1375,291,1885,408,671,1041,378,491,220,186,1542,731,303]; %eil101,最优解629，已求到634 % x=[41.0, 35.0, 55.0, 55.0, 15.0, 25.0, 20.0, 10.0, 55.0, 30.0, 20.0, 50.0, 30.0, 15.0, 30.0, 10.0, 5.0, 20.0, 15.0, 45.0, 45.0, 45.0, 55.0, 65.0, 65.0, 45.0, 35.0, 41.0, 64.0, 40.0, 31.0, 35.0, 53.0, 65.0, 63.0, 2.0, 20.0, 5.0, 60.0, 40.0, 42.0, 24.0, 23.0, 11.0, 6.0, 2.0, 8.0, 13.0, 6.0, 47.0, 49.0, 27.0, 37.0, 57.0, 63.0, 53.0, 32.0, 36.0, 21.0, 17.0, 12.0, 24.0, 27.0, 15.0, 62.0, 49.0, 67.0, 56.0, 37.0, 37.0, 57.0, 47.0, 44.0, 46.0, 49.0, 49.0, 53.0, 61.0, 57.0, 56.0, 55.0, 15.0, 14.0, 11.0, 16.0, 4.0, 28.0, 26.0, 26.0, 31.0, 15.0, 22.0, 18.0, 26.0, 25.0, 22.0, 25.0, 19.0, 20.0, 18.0, 35.0]; % y=[49.0, 17.0, 45.0, 20.0, 30.0, 30.0, 50.0, 43.0, 60.0, 60.0, 65.0, 35.0, 25.0, 10.0, 5.0, 20.0, 30.0, 40.0, 60.0, 65.0, 20.0, 10.0, 5.0, 35.0, 20.0, 30.0, 40.0, 37.0, 42.0, 60.0, 52.0, 69.0, 52.0, 55.0, 65.0, 60.0, 20.0, 5.0, 12.0, 25.0, 7.0, 12.0, 3.0, 14.0, 38.0, 48.0, 56.0, 52.0, 68.0, 47.0, 58.0, 43.0, 31.0, 29.0, 23.0, 12.0, 12.0, 26.0, 24.0, 34.0, 24.0, 58.0, 69.0, 77.0, 77.0, 73.0, 5.0, 39.0, 47.0, 56.0,