求解带时间窗的VRP问题

求解带时间窗的VRP问题，使用了两阶段法，内部嵌套节约算法。

clc
clear all
p1=0.9;
customer=xlsread（‘customer.xls‘）;           %需求点信息
facility=xlsread（‘facility.xls‘）;           %设施点信息
facilityposition=facility（:2:3）;               %设施坐标
customerposition=customer（:2:3）;               %需求点坐标
position=[facilityposition;customerposition];   
xlswrite（‘position.xls‘position）
position1=[position（:1） position（:2）];
distMatrix=dists（position1）;                    %计算得出的两点之间的距离
xlswrite（‘distMatrix.xls‘distMatrix）

ttimeu=fix（distMatrix）;                         %两点之间的距离
%%%%%%%%%%%%%%%%%%%%%%%%%固定数据%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Qofcar=200;                                     %车辆容量
costofallcar=5000;                               %车辆固定成本
costofunitdistance=9;                            %单位距离成本
tanpaifangyinzi=1;               

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2015-05-10 09:56  测试节约算法2\
     文件        1882  2015-05-07 06:51  测试节约算法2\ceshi.m
     文件       10752  2015-05-06 23:52  测试节约算法2\customer.xls
     文件       30452  2009-10-09 17:26  测试节约算法2\cw.m
     文件       14848  2015-05-10 09:56  测试节约算法2\distMatrix.xls
     文件         357  2009-09-24 03:29  测试节约算法2\dists.m
     文件        7680  2015-05-06 23:51  测试节约算法2\facility.xls
     文件        2274  2009-10-08 20:21  测试节约算法2\picture.m
     文件        8192  2015-05-10 09:56  测试节约算法2\position.xls
     文件        5295  2009-10-09 11:28  测试节约算法2\tabu.m