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

1、单一物流中心选址重心法程序设计 图 供应地和需求地分布情况重心法是一种模拟方法。这种方法将物流系统中的需求点和资源点看成是分布在某一平面范围内的物流系统，各点的需求量和资源量分别看成是物体的重量，物体系统的重心作为物流网点的最佳设置点，利用求物体系统重心的方法来确定物流网点的位置。（1）式中：Vii点运输量；Ri待定物流中心到i点的运输费率；di待定物流中心到i点的距离。求解算法数值分析法（重心法）1) 设供应点和需求点所在地的坐标为（Xi,Yi）,待定物流中心的位置坐标为（X0,Y0）则 (2)2) 将（2）式代入（1）式，然后求运输总费用TC对X0和Y0的偏导数，并令其等于零。(3)(4)

2、上述两式中仍含有未知数di,因此一次不能求得X0和Y0（解析解），需要通过迭代收敛法得到数值解。迭代收敛法具体步骤：1、 先用重心公式估算初始选址点（大致位置）：(5)(6)2、将X0和Y0代入公式2，计算di（i=1,2，m+n）；3、将di代入公式3和4，解出修正值X0和Y0；4、根据修正值X0和Y0，再重新计算di；5、重复步骤3和4，直至X0和Y0的值在连续迭代过程中不再变化，即X00，Y00,即得到精确仓库选址位置，继续计算无意义。程序设计具体步骤：Step1:利用几何重心公式（5）和（6）估算初始点X0，Y0Step2:将X0，Y0代入距离公式（2），计算di（i=1,2，5）St

3、ep3:将di（i=1,2，5）代入公式3和公式4，修正值X0和Y0Step4:重复步骤2和3，直至修正值X0和Y0的值在连续迭代过程中不再变化。（用C+程序进行编译），程序如下所示：#include <iostream.h>#include <math.h>#include <stdlib.h>void main() int n; float *X,*Y,*V,*R; cout<<"请输入需求点和供应点总个数：" cin>>n; X=new floatn; Y=new floatn; V=new floatn;

4、R=new floatn; cout<<"请输入各结点的横坐标：" for (int i=0;i<n;i+) cin>>Xi; cout<<"请输入各结点的纵坐标：" 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&

5、gt;>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; avey=sumy/sumqr; cout<<"初始物流中心地址为："<<"X="<<avex<<' '<<"Y="<<avey<<endl; fl

6、oat *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-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

7、 (j=50) break; else cout<<"第"<<j+1<<"次迭代："<<"tempX="<<xzavex<<' '<<"tempY="<<xzavey<<endl; avex=xzavex; avex=xzavey; xzsumx=0.0; xzsumy=0.0; xzsumqr=0.0; j+; cout<<"物流中心最佳地址为："<<"