物资配送matlab

matlab 数学建模 物资配送模型求解，最优配送方式选择

function [cost]=calcost（gpool） 
load distvrptw.txt; 
distance=distvrptw; 
time=distance/50;  
load1=0;load2=0;load3=0; 
g=[2 1.5 4.5 3 1.5 4 2.5 3];  
t=[1 4;4 6;1 2;4 7;3 5.5;2 5;5 8;1.5 4]; 
xt=[1 2 1 3 2 2.5 3 0.8]; 
[mn]=size（gpool）; 
cost=zeros（m1）; 
for ngpooli=1:m      
    r=gpool（ngpooli:）;     
    center=find（r==0）;     
    rk=r（r~=0）;     
    distcost=0; 
for i=1:n-1 
    distcost=distcost+distance（r（i）+1r（i+1）+1）; 
end
   timecost=0;
   dt=[]; 
   dt（1）=time（r（1）+1r（2）+1）; 
if center（2）>3
for i=2:center（2）-2 
    dt（i）=dt（i-1）+xt（rk（i-1））+time（rk（i-1）+1rk（i）+1）; 
end
end
    dt（center（2）-1）=time（r（center（2））+1r（center（2）+1）+1）; 
if center（3）-center（2）>=3 
for i=center（2）:center（3）-3  
    dt（i）=dt（i-1）+xt（rk（i-1））+time（rk（i-1）+1rk（i）+1）; 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2017-05-29 17:55  物资配送\
     文件        1766  2017-05-28 00:08  物资配送\calcost.m
     文件         139  2017-05-28 00:10  物资配送\calfitvalue.m
     文件         266  2017-05-28 00:15  物资配送\crossover.m
     文件         269  2017-05-28 00:27  物资配送\distvrptw.txt
     文件         218  2017-05-28 00:16  物资配送\eliminate.m
     文件         576  2017-05-28 00:52  物资配送\gavrptw.m
     文件         683  2017-05-28 00:20  物资配送\mutation.m
     文件        1020  2017-05-28 00:22  物资配送\rcrossover.m
     文件         903  2017-05-28 00:48  物资配送\vrpmain.m