MATLAB蚁群算法求解带时间窗的旅行商TSPTW问题代码实例

MATLAB蚁群算法求解带时间窗的旅行商TSPTW问题代码实例

蚁群算法编程求解TSPTW问题实例：

在经纬度范围为(121, 43)到(123, 45)的矩形区域内，散布着1个商家（编号1）和25个顾客点（编号为2_{26），各个点的坐标见表1。节点1表示商家，2}26表示与商家对应的顾客点。假定车辆行驶速度为30km/h，配送时间从配送员从商家位置出发算起，然后送到各个顾客点，最后返回商家位置。每个顾客点有到达的时间窗范围，不在规定的时间窗范围内到达会有一定的时间惩罚成本。早到的时间惩罚系数为0.2，迟到的时间惩罚系数为0.5，单位距离成本为1。一旦车辆出发则处于配送状态，直到所有配送任务完成停止计算，用蚁群算法求总成本最小（距离成本+时间惩罚成本）的最优路径。
编号 X坐标
（经度） Y坐标
（纬度） 最早时间窗
（分钟） 最晚时间窗
（分钟）
1 122.044 44.4703 150 1860
2 121.2915 44.5052 50 2810
3 122.9038 44.6702 450 2810
4 121.082 44.1601 400 1290
5 122.8507 44.2793 300 1100
6 122.9272 43.3611 400 2430
7 121.1113 44.6827 400 2810
8 121.801 43.3733 450 2810
9 121.9528 44.5705 100 2240
10 122.7934 44.3812 450 2240
11 121.9571 44.0387 300 2240
12 121.6636 44.5933 0 2430
13 121.5183 44.3164 100 2620
14 122.2505 43.0642 400 1480
15 122.671 44.243 300 1480
16 122.8145 44.6338 0 2430
17 121.9659 43.9129 150 1480
18 122.0954 43.4043 400 2240
19 121.1637 44.4004 350 1860
20 121.6539 43.2019 500 1670
21 121.241 43.6942 300 1860
22 122.227 44.7439 500 2050
23 121.0417 43.9196 350 1290
24 121.5933 44.2417 200 3000
25 121.2904 44.1993 250 1670
26 122.5063 44.7392 400 1100

先看下求解结果！
运行结果:

优化结束！
最优总成本为：1230.2245
其中：
距离成本为：978.1766
早到时间惩罚成本为：203.8713
晚到时间惩罚成本为：48.1766
最优配送路径为：
1->9->22->26->16->3->10->5->15->6->14->18->8->20->21->23->4->25->17->11->24->13->2->19->7->12->1
到达各点的时间为：
0->30->102->164->236->258->327->352->393->597->761->844->910->960->1103->1169->1224->1271->1433->1461->1554->1577->1643->1680->1743->1868->1956
下面进行程序演示！