基于粒子群优化算法的TSP问题研究

function [CP BESTVALUE BESTX]=cityPositionF(sim_model) %返回城市坐标 switch sim_model case '14' %最优值为29.3405 CP=[16.47,96.10;16.47,94.44;20.09,92.54;22.39,93.37;25.23,97.24;22.00,96.05;20.47,97.02;17.20,96.29;16.30,97.38;14.05,98.12;16.53,97.38;21.52,95.59;19.41,97.13;20.09,92.55]; BESTVALUE=29.3405; BESTX=[12 6 5 4 3 14 2 1 10 9 11 8 13 7 12]; case 'bayg29' %TSPLIB提供最优解9074.1 CP=[1150.0,1760.0;630.0,1660.0;40.0,2090.0;750.0,1100.0;750.0,2030.0;1030.0,2070.0;1650.0,650.0;1490.0,1630.0;790.0,2260.0;710.0,1310.0;840.0,550.0;1170.0,2300.0;970.0,1340.0;510.0,700.0;750.0,900.0;1280.0,1200.0;230.0,590.0;460.0,860.0;1040.0,950.0;590.0,1390.0;830.0,1770.0;490.0,500.0;1840.0,1240.0;1260.0,1500.0;1280.0,790.0;490.0,2130.0;1460.0,1420.0;1260.0,1910.0;360.0,1980.0]; BESTX=[1 28 6 12 9 26 3 29 5 21 2 20 10 4 15 18 14 17 22 11 19 25 7 23 8 27 16 13 24 1]; BESTVALUE=9074.1; case 'att48' %TSPLIB提供最优解33523.7085 CP=[6734 1453;2233 10;5530 1424;401 841;3082 1644;7608 4458;7573 3716;7265 1268;6898 1885;1112 2049;5468 2606;5989 2873;4706 2674;4612 2035;6347 2683;6107 669;7611 5184;7462 3590;7732 4723;5900 3561;4483 3369;6101 1110;5199 2182;1633 2809;4307 2322;675 1006;7555 4819;7541 3981;3177 756;7352 4506;7545 2801;3245 3305;6426 3173;4608 1198;23 2216;7248 3779;7762 4595;7392 2244;3484 2829;6271 2135;4985 140;1916 1569;7280 4899;7509 3239;10 2676;6807 2993;5185 3258;3023 1942]; BESTX=[1 8 38 31 44 18 7 28 6 37 19 27 17 43 30 36 46 33 20 47 21 32 39 48 5 42 24 10 45 35 4 26 2 29 34 41 16 22 3 23 14 25 13 11 12 15 40 9 1]; BESTVALUE=33523.7085; case 'eil51' %最优解426，TSPLIB提供最优解429.9833 CP=[37 52;49 49;52 64;20 26;40 30;21 47;17 63;31 62;52 33;51 21;42 41;31 32;5 25;12 42;36 16;52 41;27 23;17 33;13 13;57 58;62 42;42 57;16 57;8 52;7 38;27 68;30 48;43 67;58 48;58 27;37 69;38 46;46 10;61 33;62 63;63 69;32 22;45 35;59 15;5 6;10 17;21 10;5 64;30 15;39 10;32 39;25 32;25 55;48 28;56 37;30 40]; BESTX=[1 22 8 26 31 28 3 36 35 20 2 29 21 16 50 34 30 9 49 10 39 33 45 15 44 42 40 19 41 13 25 14 24 43 7 23 48 6 27 51 46 12 47 18 4 17 37 5 38 11 32 1]; BESTVALUE=426; case 'berlin52' %最优解7542，TSPLIB提供最优解7544.4 CP=[565.0 575.0;25.0 185.0;345.0 750.0;945.0 685.0;845.0 655.0;880.0 660.0;25.0 230.0;525.0 1000.0;580.0 1175.0;650.0 1130.0;1605.0 620.0 ;1220.0 580.0;1465.0 200.0;1530.0 5.0;845.0 680.0;725.0 370.0;145.0 665.0;415.0 635.0;510.0 875.0 ;560.0 365.0;300.0 465.0;520.0 585.0;480.0 415.0;835.0 625.0;975.0 580.0;1215.0 245.0;1320.0 315.0;1250.0 400.0;660.0 180.0;410.0 250.0;420.0 555.0;575.0 665.0;1150.0 1160.0;700.0 580.0;685.0 595.0;685.0 610.0;770.0 610.0;795.0 645.0;720.0 635.0;760.0 650.0;475.0 960.0;95.0 260.0;875.0 920.0;700.0 500.0;555.0 815.0;830.0 485.0;1170.0 65.0;830.0 610.0;605.0 625.0;595.0 360.0;1340.0 725.0;1740.0 245.0]; BESTX=[1;49;32;45;19;41;8;9;10;43;33;51;11;52;14;13;47;26;27;28;12;25;4;6;15;5;24;48;38;37;40;39;36;35;34;44;46;16;29;50;20;23;30;2;7;42;21;17;3;18;31;22;1]'; BESTVALUE=7542; case 'st70' %最优解675，TSPLIB提供最优解678.5975 CP=[64 96;80 39;69 23;72 42;48 67;58 43;81 34;79 17;30 23;42 67;7 76;29 51;78 92;64 8;95 57;57 91;40 35;68 40;92 34;62 1;28 43;76 73;67 88;93 54;6 8;87 18;30 9;77 13;78 94;55 3;82 88;73 28;20 55;27 43;95 86;67 99;48 83;75 81;8 19;20 18;54 38;63 36;44 33;52 18;12 13;25 5;58 85;5 67;90 9;41 76;25 76;37 64;56 63;10 55;98 7;16 74;89 60;48 82;81 76;29 60;17 22;5 45;79 70;9 100;17 82;74 67;10 68;48 19;83 86;84 94]; BESTX=[1;36;29;13;70;35;31;69;38;59;22;66;63;57;15;24;19;7;2;4;18;42;32;3;8;26;55;49;28;14;20;30;44;68;27;46;25;45;39;61;40;9;17;43;41;6;53;5;10;52;60;12;34;21;33;62;54;48;67;11;64;65;56;51;50;58;37;47;16;23;1]'; BESTVALUE=675; case 'eil76' %最优解538，TSPLIB提供最优解545.3876 CP=[22 22;36 26;21 45;45 35;55 20;33 34;50 50;55 45;26 59;40 66;55 65;35 51;62 35;62 57;62 24;21 36;33 44;9 56;62 48;66 14;44 13;26 13;11 28;7 43;17 64;41 46;55 34;35 16;52 26;43 26;31 76;22 53;26 29;50 40;55 50;54 10;60 15;47 66;30 60;30 50;12 17;15 14;16 19;21 48;50 30;51 42;50 15;48 21;12 38;15 56;29 39;54 38;55 57;67 41;10 70;6 25;65 27;40 60;70 64;64 4;36 6;30 20;20 30;15 5;50 70;57 72;45 42;38 33;50 4;66 8;59 5;35 60;27 24;40 20;40 37;40 40]; BESTX=[1;33;63;16;3;44;32;9;39;72;58;10;31;55;25;50;18;24;49;23;56;41;43;42;64;22;61;21;47;36;69;71;60;70;20;37;5;15;57;13;54;19;14;59;66;65;38;11;53;7;35;8;46;34;52;27;45;29;48;30;4;75;76;67;26;12;40;17;51;6;68;2;74;28;62;73;1]'; BESTVALUE=538; case 'pr76' %最优解108159，TSPLIB提供最优解108159.4 CP=[3600 2300;3100 3300;4700 5750;5400 5750;5608 7103;4493 7102;3600 6950;3100 7250;4700 8450;5400 8450;5610 10053;4492 10052;3600 10800;3100 10950;4700 11650;5400 11650;6650 10800;7300 10950;7300 7250;6650 6950;7300 3300;6650 2300;5400 1600;8350 2300;7850 3300;9450 5750;10150 5750;10358 7103;9243 7102;8350 6950;7850 7250;9450 8450;10150 8450;10360 10053;9242 10052;8350 10800;7850 10950;9450 11650;10150 11650;11400 10800;12050 10950;12050 7250;11400 6950;12050 3300;11400 2300;10150 1600;13100 2300;12600 3300;14200 5750;14900 5750;15108 7103;13993 7102;13100 6950;12600 7250;14200 8450;14900 8450;15110 10053;13992 10052;13100 10800;12600 10950;14200 11650;14900 11650;16150 10800;16800 10950;16800 7250;16150 6950;16800 3300;16150 2300;14900 1600;19800 800;19800 10000;19800 11900;19800 12200;200 12200;200 1100;200 800]; BESTX=[1;76;75;2;3;4;5;6;7;8;9;10;11;12;13;14;74;15;16;17;18;37;36;38;39;40;34;35;33;32;29;30;31;19;20;26;27;28;43;42;54;53;52;55;56;57;58;59;60;41;61;62;63;64;73;72;71;65;66;51;49;50;67;70;68;69;47;48;44;45;46;24;25;21;22;23;1]'; BESTVALUE=108159; case 'kroa100' %最优解21282，TSPLIB提供最优解21285.4 CP=[1380 939;2848 96;3510 1671;457 334;3888 666;984 965;2721 1482;1286 525;2716 1432;738 1325;1251 1832;2728 1698;3815 169;3683 1533;1247 1945;123 862;1234 1946;252 1240;611 673;2576 1676;928 1700;53 857;1807 1711;274 1420;2574 946;178 24;2678 1825;1795 962;3384 1498;3520 1079;1256 61;1424 1728;3913 192;3085 1528;2573 1969;463 1670;3875 598;298 1513;3479 821;2542 236;3955 1743;1323 280;3447 1830;2936 337;1621 1830;3373 1646;1393 1368;3874 1318;938 955;3022 474;2482 1183;3854 923;376 825;2519 135;2945 1622;953 268;2628 1479;2097 981;890 1846;2139 1806;2421 1007;2290 1810;1115 1052;2588 302;327 265;241 341;1917 687;2991 792;2573 599;19 674;3911 1673;872 1559;2863 558;929 1766;839 620;3893 102;2178 1619;3822 899;378 1048;1178 100;2599 901;3416 143;2961 1605;611 1384;3113 885;2597 1830;2586 1286;161 906;1429 134;742 1025;1625 1651;1187 706;1787 1009;22 987;3640 43;3756 882;776 392;1724 1642;198 1810;3950 1558]; BESTX=[1;47;93;28;67;58;61;51;87;25;81;69;64;40;54;2;44;50;73;68;85;82;95;13;76;33;37;5;52;78;96;39;30;48;100;41;71;14;3;43;46;29;34;83;55;7;9;57;20;12;27;86;35;62;60;77;23;98;91;45;32;11;15;17;59;74;21;72;10;84;36;99;38;24;18;79;53;88;16;94;22;70;66;26;65;4;97;56;80;31;89;42;8;92;75;19;90;49;6;63;1]'; BESTVALUE=21282; case 'kroc100' %最优解20749，TSPLIB提供最优解20750.8 CP=[1357 1905;2650 802;1774 107;1307 964;3806 746;2687 1353;43 1957;3092 1668;185 1542;834 629;40 462;1183 1391;2048 1628;1097 643;1838 1732;234 1118;3314 1881;737 1285;779 777;2312 1949;2576 189;3078 1541;2781 478;705 1812;3409 1917;323 1714;1660 1556;3729 1188;693 1383;2361 640;2433 1538;554 1825;913 317;3586 1909;2636 727;1000 457;482 1337;3704 1082;3635 1174;1362 1526;2049 417;2552 1909;3939 640;219 898;812 351;901 1552;2513 1572;242 584;826 1226;3278 799;86 1065;14 454;1327 1893;2773 1286;2469 1838;3835 963;1031 428;3853 1712;1868 197;1544 863;457 1607;3174 1064;192 1004;2318 1925;2232 1374;396 828;2365 1649;2499 658;1410 307;299