qsolve.asv

TSP问题,即旅行商问题,matlab环境下可以运行

function [inventorycost]=Qsolve(v,l)
%v,l分别是路线的集合，数组及顾客的个数
%计算路线总长度
%所有参数部分
demand=[155 98 474
        178 191 365
        50 65 200];
    depot=[38 133
           22 76
           174 193];
       UD=[3 8 6 7];
       
c=25;
cm=1;
h=0.5;
A=20;
hs=2;






length=0;
for i=2:l
length=distance(demand,v(i),v(i+1))+length;
end
%计算路线总需求
D=sum(demand(v(2:l+1),3));
%计算路线上提前期的平均需求
u=0;
for i=2:l
   u=UD(v(i))+u;
end

R0=0;
BR0=(u-R0)^2/(2*u);
Q0=sqrt((2*D*(cm+l+c+hs*BR0))/h);
R1=u-(u*h*Q0)/(hs*D);
BR1=(u-R1)^2/(2*u);
Q1=sqrt((2*D*(cm+l+c+hs*BR1))/h);
while abs(Q1-Q0)>0.01
   R0=R1;
   BR0=(u-R0)^2/(2*u);
  Q0=sqrt((2*D*(cm+l+c+hs*BR0))/h);
   R1=u-(u*h*Q0)/(hs*D);
   BR1=(u-R1)^2/(2*u);
   Q1=sqrt((2*D*(cm+l+c+hs*BR1))/h);
end
%计算库存成本
inventorycost=((c+cm*length)*D)/Q1+(Q1/2+R1-u/2)*h+(D/Q1)*A+(hs*BR1*D)/Q1;