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