基于人工蜂群算法的旅行商问题（TSP）【matlab代码】

function [map] = create_tsp_map() %% berlin 52 % x = [565 25 345 945 845 880 25 525 580 650 1605 1220 1465 1530 845 725 145 415 510 560 300 520 480 835 975 1215 1320 1250 660 410 420 575 1150 700 685 685 770 795 720 760 475 95 875 700 555 830 1170 830 605 595 1340 1740]; % y = [575 185 750 685 655 660 230 1000 1175 1130 620 580 200 5 680 370 665 635 875 365 465 585 415 625 580 245 315 400 180 250 555 665 1160 580 595 610 610 645 635 650 960 260 920 500 815 485 65 610 625 360 725 245]; %% kroA100 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 % 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 % 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]; %% kroD100 % x = [2995 202 981 1346 781 1009 2927 2982 555 464 3452 571 2656 1623 2067 1725 3600 1109 366 778 386 3918 3332 2597 811 241 2658 394 3786 264 2050 3538 1646 2993 547 3373 460 3060 1828 1021 2347 3535 1529 1203 1787 2740 555 47 3935 3062 387 2901 931 1766 401 149 2214 3805 1179 1017 2834 634 1819 1393 1768 3023 3248 1632 2223 3868 1541 2374 1962 3007 3220 2356 1604 2028 2581 2221 2944 1082 997 2334 1264 1699 235 2592 3642 3599 1766 240 1272 3503 80 1677 3766 3946 1994 278]; % y = [264 233 848 408 670 1001 1777 949 1121 1302 637 1982 128 1723 694 927 459 1196 339 1282 1616 1217 1049 349 1295 1069 360 1944 1862 36 1833 125 1817 624 25 1902 267 781 456 962 388 1112 581 385 1902 1101 1753 363 540 329 199 920 512 692 980 1629 1977 1619 969 333 1512 294 814 859 1578 871 1906 1742 990 697 354 1944 389 1524 1945 1568 706 1736 121 1578 632 1561 942 523 1090 1294 1059 248 699 514 678 619 246 301 1533 1238 154 459 1852 165]; %% kroE100 % x = [3477 91 3972 198 1806 538 3430 2186 1513 2143 53 3404 1034 2823 3104 3232 2790 374 741 3083 3502 1280 3326 217 2503 3527 739 3548 48 1419 1689 3468 1628 382 3029 3646 285 1782 1067 2849 920 1741 876 2753 2609 3941 3613 1754 2916 2445 3825 2779 201 2502 765 3105 1937 3364 3702 2164 3019 3098 3239 3359 2081 1398 618 1878 3803 397 3035 2502 3230 3479 958 3423 78 96 3431 2053 3048 571 3393 2835 144 923 989 3061 2977 1668 878 678 1086 640 3551 106 2243 3796 2643 48]; % y = [949 1732 329 1632 733 1023 1088 766 1646 1611 1657 1307 1344 376 1931 324 1457 9 146 1938 1067 237 1846 38 1172 41 1850 1999 154 872 1223 1404 253 872 1242 1758 1029 93 371 1214 1835 712 220 283 1286 258 523 559 1724 1820 1101 435 693 1274 833 1823 1400 1498 1624 1874 189 1594 1376 1693 1011 1100 1953 59 886 1217 152 146 380 1023 1670 1241 1066 691 78 1461 1 1711 782 1472 1185 108 1997 1211 39 658 715 1599 868 110 1673 1267 1332 1401 1320 267]; %% kroA150 % 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 3477 91 3972 198 1806 538 3430 2186 1513 2143 53 3404 1034 2823 3104 3232 2790 374 741 3083 3502 1280 3326 217 2503 3527 739 3548 48 1419 1689 3468 1628 382 3029 3646 285 1782 1067 2849 920 1741 876 2753 2609 3941 3613 1754 2916 2445]; % 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 949 1732 329 1632 733 1023 1088 766 1646 1611 1657 1307 1344 376 1931 324 1457 9 146 1938 1067 237 1846 38 1172 41 1850 1999 154 872 1223 1404 253 872 1242 1758 1029 93 371 1214 1835 712 220 283 1286 258 523 559 1724 1820]; %% kroB150 % 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 3825 2779 201 2502 765 3105 1937 3364 3702 2164 3019 3098 3239 3359 2081 1398 618 1878 3803 397 3035 2502 3230 3479 958 3423 78 96 3431 2053 3048 571 3393 2835 144 923 989 3061 2977 1668 878 678 1086 640 3551 106 2243 3796 2643 48]; % 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 1101 435 693 1274 833 1823 1400 1498 1624 1874 189 1594 1376 1693 1011 1100 1953 59 886 1217 152 146 380 1023 1670 1241 1066 691 78 1461 1 1711 782 1472 1185 108 1997 1211 3