物流中心选址重心法c程序设计

1、单一物流中心选址重心法C+程序设计(选址)某企业有五个供应商，他们的位置分别是(4+i,4+i ), (12+i,4+i ), ( 2+i,7+i ),(11+i,11+i ), (7+i,14+i ),该企业从这五个供应商处的采购量每年为4+i、3+i、2+i、4+i、1+i ,假设运输费率相同，i为序号的末位数。请用计算机编写一个通用程序，完成下面的题目。(1)用城市距离为该企业推荐一个地址，使该企业的费用最省。(2)将第1问中的结果作为一个初始解，用欧氏距离进行重新优化，推荐一个最优的位置。要求：(1)前后两次结果之差小于0.001 ;(2)程序的通用性高，可以满足任何多个供应商；(3)

2、程序的可读性强；(4)操作界面友好；(5)在手册中要记录所有的程序、运行结果及相关界面图形。解：F Kxl.y.)图供应地和需求地分布情况重心法是一种模拟方法。这种方法将物流系统中的需求点和资源点看成是分布在 某一平面范围内的物流系统，各点的需求量和资源量分别看成是物体的重量，物体系统 的重心作为物流网点的最佳设置点，利用求物体系统重心的方法来确定物流网点的位 置。运输总费用min TC£ ViRidi (1)二1式中：V i点运输量;R一待定物流中心到i点的运输费率;di待定物流中心到i点的距离。求解算法一数值分析法（重心法）1） 设供应点和需求点所在地的坐标为 （X,Yi）,待定

3、物流中心的位置坐标为（Xo,Y。）则di = , Xi -X0 2 Yi -Y0 2 2） 将（2）式代入（1）式，然后求运输总费用TC对Xo和Yo的偏导数，并令其等于 零。TCMR%CXi -X。）2 +（Y Yo）2,:TC元一 MR (Xi -Xo) =0di、MRXidiVi Ri -1 =0diXoYo' (ViRiYi di)“(ViR di)x (ViRiXi di)(3) (MR di)上述两式中仍含有未知数di,因此一次不能求得Xo和Yo （解析解），需要通过迭代收敛 法得到数值解。迭代收敛法具体步骤：1、先用重心公式估算初始选址点（大致位置）：Xo工(ViRXi)%

4、 (MR)YoZ (ViRY) (6)x (ViRi)2、将Xo和Yo代入公式2,计算di （i=1,2,，m+n）;3、将di代入公式3和4,解出修正值Xo和Yo；4、根据修正值Xo和Yo,再重新计算di；5、重复步骤3和4,直至Xo和Yo的值在连续迭代过程中不再变化，即 XoQo,Yo= o,即得到精确仓库选址位置，继续计算无意义程序设计具体步骤：Stepl:利用几何重心公式（5）和（6）估算初始点Xo, Yo float *X,*Y ,*V,*R;X0£ (ViRXi)9 (ViRi)v " (MRY)Y0 =,、(Vi Ri)Step2:将Xo, Yo代入距离公式(

5、2),计算di (i=1,2,，5)di =、(Xi 二Xo)2一(Yi 二 Y0)2d2 = ；(X2 -X0)2 (Y2 -Yo)2d3 = .,(X3 =Xo)2一(Y3 =Yo)2d4 = .(X4 =Xo)L(Y4 二Yo)2 199d5 = . (X5 -Xo)2 (Y5 -Yo)2Step3:将di (i=1,2 ,，5)代入公式3和公式4,修正值Xo和Y0' (ViRiXi di)X 0 =z“ (MR di)7” (ViRiY di)YoZ (ViR/di)Step4:重复步骤2和3,直至修正值Xo和Yo的值在连续迭代过程中不再变化 行编译)，程序如下所示：#incl

6、ude <iostream.h>#include <math.h>#include <stdlib.h>void main()（用C+利序进int n;II.cout<<"请输入需求点和供应点总个数：cin>>n;X=new floatn;Y=new floatn;V=new floatn;R=new floatn;cout<<"请输入各结点的横坐标："for (int i=0;i<n;i+) cin>>Xi;cout<<"请输入各结点的纵坐标：&quo

7、t;for (i=0;i<n;i+) cin>>Yi;cout<<"请输入各结点的供货量或需求量：for (i=0;i<n;i+) cin>>Vi;cout<<”请输入各结点的内/外向运输费率："for (i=0;i<n;i+)cin>>Ri;float sumx=0.0,sumy=0.0,sumqr=0.0;float avex,avey;for (i=0;i<n;i+)sumx+=Vi*Ri*Xi;sumy+=Vi*Ri*Yi;sumqr+=Vi*Ri;avex=sumx/sumqr;a

8、vey=sumy/sumqr;cout<<"初始物流中心地址为："<<"X="<<avex<<' '<<"Y="<<avey<<endl;float *d;d=new floati;float xzsumx=0.0,xzsumy=0.0,xzsumqr=0.0;float xzavex,xzavey;int j=0;while (1)for (i=0;i<n;i+)di=sqrt(avex-Xi)*(avex-Xi)+(avey-

9、Yi)*(avey-Yi);for (i=0;i<n;i+)xzsumx+=Vi*Ri*Xi/di;xzsumy+=Vi*Ri*Yi/di;xzsumqr+=Vi*Ri/di;xzavex=xzsumx/xzsumqr;xzavey=xzsumy/xzsumqr;if (j=50)break;elsecout<<"第"<<j+1<<"次迭代："<<"tempX="<<xzavex<<' '<<"tempY="<<xzavey<<endl;avex=xzavex;avex=xzavey;xzsumx=0.0;xzsumy=0.0;x