W公司配送案例浅分析

W公司配送案例浅分析PresentationT第2小组案例一背景介绍本案例仅选取公司的四个配送中心进行研究。这四个配送中心具有以下特点：（1）配送网络比较广阔。（2）客户数目不多，但比较固定，几乎需要每日配送。（3）每个客户需求量波动不大。（4）货物种类较少，需求量较少，配送周期短。（5）客户与公司虽然联系较多，但仍处于关系建立阶段，客户问题时有发生。W公司配送系统现状企业的配送系统中有四个配送中心（M=1,2,3,4），30个客户点（N=1,2，，30），单个配送中心基本建设费用R=600万元。单位产品单位路程运费C1为固定的1.4元/(t.km)，车辆单位路程的费用C2为固定的0.6元/km。配送中心网点布局

1、配送中心坐标配送中心编号XY119352758132870470242,配送中心容量、费用情况配送中心设施位置容量库存费用/(元/天。T)服务顾客数111010072110100831101204411010011配送中心经营状况

（每个配送中心有车辆6辆，每辆车的承载能力为12t。配送中心总人员30人左右。）配送中心编号总成本/元总收入/元客户总数满意客户数总收益/元123080003020000757120002368000040440008736400031653000177000044117000439890004935000119946000配送中心的行车路线配送中心行车路线配送中心行车路线1配送中心1--5--28--23--配送中心1配送中心1--14--11--配送中心1配送中心1--18--3--配送中心13配送中心3--30--2--9--配送中心3配送中心3--27--配送中心32配送中心2--6--22--配送中心2配送中心2--16--17--配送中心2配送中心2--10--20--配送中心2配送中心2--7--8--配送中心24配送中心4--26--12--配送中心4配送中心4--25--19--15--配送中心4配送中心4--4--24--配送中心4配送中心4--13--29--21--配送中心4配送中心最初配送情况（注意：客户1的遗漏）配送中心最初配送成本详细情况成本项目配送中心1配送中心2配送中心3配送中心4建设成本/(元/年)600000600000600000600000库存总成本100800015120006480001728000运输总成本69000015530004050001642000顾客不满意所造成的成本1000015000019000总成本2308000368000016530003989000目录分析步骤发现问题确认目标解决方法优化措施问题的发现配送路线规划问题车辆配置问题系统整体优化车载问题服务质量问题遗传算法解决配送路线规划车辆重调度解决车辆配置问题解决车载问题提升服务质量建立关系综合调整优化系统认确标目遗传算法在MATLAB的操作输入数据编程过程相关代码functionD=distance(P)fori=1:length(P)forj=1:length(P)D(i,j)=sqrt((P(i,1)-P(j,1))^2+(P(i,2)-P(j,2))^2);endendfunction[s,q]=TSP(D,n)%D距离矩阵，n配送点个数ticm=input('种群大小，偶数最佳m=');p=input('遗传代数=');%产生初始种群fori=1:me=1:n;forj=1:nt=unidrnd(length(e));q(i,j)=e(t);e(t)=[];endends=len(q,D,n,m);%计算路程长度fori=1;pq=newq(s,q,m);%产生新种群q=change(q,m,n);%交叉q=vf(q,n,m);%变异s=len(q,D,n,m);%计算路程长度end[s0,b]=min(s);fprintf('总路程为s=%4.2f\n',s0);fprintf('站点顺序为\n');fprintf('%d',q(b,:));toc%计算路程函数functions=len(q,D,n,m)fori=1:ms(i)=0;forj=1:n-1s(i)=s(i)+D(q(i,j),q(i,j+1));endend%选择目标函数小的个体，产生新种群functionq=newq(s,q,m)t=0.2;%淘汰的比例fori=1:m-1forj=i:mifs(i)<s(j)k=s(j);s(j)=s(i);s(i)=k;Q=q(j,:);q(j,:)=q(i,:);q(i,:)=Q;endendendf=fix(m*t);%淘汰个数%将目标函数值大的前f个个体用目标函数值小的后f个个体替换fori=1:fq(i,:)=q(m-f+i,:);end%打乱顺序fori=1:2*mN1=unidrnd(m);N2=unidrnd(m);Q=q(N1,:);q(N1,:)=q(N2,:)q(N2,:)=Q;end%两两交叉操作functionq=change(q,m,n)while1N1=unidrnd(n-1);N2=unidrnd(n-1);ifN2>N1+1break;endend%采用部分匹配交叉，且一对一对交叉fori=1:2:m-1%中间部分交叉Q=q(i,N1+1:N2);q(i,N1+1:N2)=q(i+1,N1+1:N2);q(i+1,N1+1:N2)=Q;%选出可以保留的部分t1=1;t2=1;forj=1:nif(q(i,j)~=q(i+1,N1+1:N2))e1(t1)=j;t1=t1+1;endif(j>N1&j<=N2)e1(t1)=j;t1=t1+1;endif(q(i+1,j)~=q(i+1,N1+1:N2))e2(t2)=j;t2=t2+1;endif(j>N1&j<=N2)e2(t2)=j;t2=t2+1;endend%替换重复部分forj=1:nifj~=e1fork=N1+1:N2ifq(i+1,k)~=q(i,e1)q(i,j)=q(i+1,k);e1(t1)=j;t1=t1+1;break;endendendifj~=e2fork=N1+1:N2ifq(i,k)~=q(i+1,e2)q(i+1,j)=q(i,k);e2(t2)=j;t2=t2+1;break;endendendendforj=1:nifj~=e1fork=1:nifk~=q(i,e1)q(i,j)=k;e1(t1)=j;t1=t1+1;break;endendendifj~=e2fork=1:nifk~=q(i+1,e2)q(i+1,j)=k;e2(t2)=j;t2=t2+1;break;endendendendcleare1e2end%变异functionq=vf(q,n,m)F=0.01;%变异概率为0.01N1=unidrnd(m);N2=1+unidrnd(n-1);k=1+unidrnd(n-1);p0=unifrnd(0,1);%产生随机变异概率ifp0<Ffori=1:nifq(N1,i)==kq(N1,i)=q(N1,N2);break;endendq(N1,N2)=k;end第一个配送中心及配送点的坐标，分别是3、5、11、14、18、23、28、及p145 2010 4315 105 3024 1033 372 4819 35第二个配送中心及配送点的坐标，分别是6、7、8、10、16、17、20、22及p255 7876 6080 6580 9058 6065 6988 8968 8875 81第三个配送中心及配送点的坐标，分别是2、9、27、30及p336 7940 90