物流中心選址重心法c程序設(shè)計(jì)

/單一物流中心選址重心法c++程序設(shè)計(jì)（選址）某企業(yè)有五個(gè)供應(yīng)商，他們的位置分別是（4+i,4+i），（12+i,4+i），（2+i,7+i），（11+i,11+i），（7+i,14+i），該企業(yè)從這五個(gè)供應(yīng)商處的選購 量每年為4+i、3+i、2+i、4+i、1+i，假設(shè)運(yùn)輸費(fèi)率相同，i為序號(hào)的末位數(shù)。請(qǐng)用計(jì)算機(jī)編寫一個(gè)通用程序，完成下面的題目。（1）用城市距離為該企業(yè)舉薦一個(gè)地址，使該企業(yè)的費(fèi)用最省。（2）將第1問中的結(jié)果作為一個(gè)初始解，用歐氏距離進(jìn)行重新優(yōu)化，舉薦一個(gè)最優(yōu)的位置。要求：（1）前后兩次結(jié)果之差小于0.001；（2）程序的通用性高，可以滿足任何多個(gè)供應(yīng)商；（3）程序的可讀性強(qiáng)；（4）操作界面友好；（5）在手冊(cè)中要記錄全部的程序、運(yùn)行結(jié)果及相關(guān)界面圖形。解：圖供應(yīng)地和需求地分布狀況重心法是一種模擬方法。這種方法將物流系統(tǒng)中的需求點(diǎn)和資源點(diǎn)看成是分布在某一平面范圍內(nèi)的物流系統(tǒng)，各點(diǎn)的需求量和資源量分別看成是物體的重量，物體系統(tǒng)的重心作為物流網(wǎng)點(diǎn)的最佳設(shè)置點(diǎn)，利用求物體系統(tǒng)重心的方法來確定物流網(wǎng)點(diǎn)的位置。（1）式中：Vi—i點(diǎn)運(yùn)輸量；Ri—待定物流中心到i點(diǎn)的運(yùn)輸費(fèi)率；di—待定物流中心到i點(diǎn)的距離。求解算法—數(shù)值分析法（重心法）設(shè)供應(yīng)點(diǎn)和需求點(diǎn)所在地的坐標(biāo)為（Xi,Yi）,待定物流中心的位置坐標(biāo)為（X0,Y0）則(2)將（2）式代入（1）式，然后求運(yùn)輸總費(fèi)用TC對(duì)X0和Y0的偏導(dǎo)數(shù)，并令其等于零。(3)(4)上述兩式中仍含有未知數(shù)di,因此一次不能求得X0和Y0（解析解），須要通過迭代收斂法得到數(shù)值解。迭代收斂法具體步驟：先用重心公式估算初始選址點(diǎn)（大致位置）：(5)(6)2、將X0和Y0代入公式2，計(jì)算di（i=1,2，…，m+n）；3、將di代入公式3和4，解出修正值X0和Y0；4、依據(jù)修正值X0和Y0，再重新計(jì)算di；5、重復(fù)步驟3和4，直至X0和Y0的值在連續(xù)迭代過程中不再變更，即△X0≈0，△Y0≈0,即得到精確倉庫選址位置，接著計(jì)算無意義。程序設(shè)計(jì)具體步驟：Step1:利用幾何重心公式（5）和（6）估算初始點(diǎn)X0，Y0Step2:將X0，Y0代入距離公式（2），計(jì)算di（i=1,2，…，5）Step3:將di（i=1,2，…，5）代入公式3和公式4，修正值X0和Y0Step4:重復(fù)步驟2和3，直至修正值X0和Y0的值在連續(xù)迭代過程中不再變更。（用C++程序進(jìn)行編譯），程序如下所示：#include<iostream.h>#include<math.h>#include<stdlib.h>voidmain(){intn;float*X,*Y,*V,*R;cout<<"請(qǐng)輸入需求點(diǎn)和供應(yīng)點(diǎn)總個(gè)數(shù)：";cin>>n;X=newfloat[n];Y=newfloat[n];V=newfloat[n];R=newfloat[n];cout<<"請(qǐng)輸入各結(jié)點(diǎn)的橫坐標(biāo)：";for(inti=0;i<n;i++)cin>>X[i];cout<<"請(qǐng)輸入各結(jié)點(diǎn)的縱坐標(biāo)：";for(i=0;i<n;i++)cin>>Y[i];cout<<"請(qǐng)輸入各結(jié)點(diǎn)的供貨量或需求量：";for(i=0;i<n;i++)cin>>V[i];cout<<"請(qǐng)輸入各結(jié)點(diǎn)的內(nèi)/外向運(yùn)輸費(fèi)率：";for(i=0;i<n;i++)cin>>R[i];floatsumx=0.0,sumy=0.0,sumqr=0.0;floatavex,avey;for(i=0;i<n;i++){sumx+=V[i]*R[i]*X[i];sumy+=V[i]*R[i]*Y[i];sumqr+=V[i]*R[i];}avex=sumx/sumqr;avey=sumy/sumqr;cout<<"初始物流中心地址為："<<"X="<<avex<<''<<"Y="<<avey<<endl;float*d;d=newfloat[i];floatxzsumx=0.0,xzsumy=0.0,xzsumqr=0.0;floatxzavex,xzavey;intj=0;while(1){for(i=0;i<n;i++)d[i]=sqrt((avex-X[i])*(avex-X[i])+(avey-Y[i])*(avey-Y[i]));for(i=0;i<n;i++){xzsumx+=V[i]*R[i]*X[i]/d[i];xzsumy+=V[i]*R[i]*Y[i]/d[i];xzsumqr+=V[i]*R[i]/d[i];}xzavex=xzsumx/xzsumqr;xzavey=xzsumy/xzsumqr;if(j==50)break;else{cout<<"第"<<j+1<<"次迭代："<<"tempX="<<xzavex<<''<<"tempY="<<xzavey<<endl;avex=xzavex;avex=xzavey;xzsumx=0.0;x