B=[9.10 4.20 6.10 3.60 6.90 3.95 6.55 7.80 8.25];
C=[17000 32000 12500 9500 29500 21000 41300 8600 10700]; a=0; for i=1:9 for j=1:9 if j~=i
L(j)=sqrt((A(j)-A(i))^2+(B(j)-B(i))^2); else L(j)=0; end end
aver(i)=sum(L)/8; end
for k=1:3
fc2w1(k)=aver(k)*(4819569/(1128*173115))*(173115-C(k)*(1+0.05))*230; x11=[3.60 6.90]; y11=[3.90 1.00];
l1=sqrt((x11-A(k)).^2+(y11-B(k)).^2); aver1=sum(l1)/2;
fr2w1(k)=aver1*(2162535*2/(1163+1188))/173115*(173115-C(k)*(1+0.05)); f2w1(k)=fc2w1(k)+fr2w1(k)+3.25*200000*5+5000000; end for m=1:6
fc2w2(m)=(aver(m+3))*(4819569/(1128*173115))*(173115-(C(m+3))*(1-0.1))*230; x11=[3.60 6.90]; y11=[3.90 1.00];
l1=sqrt((x11-A(m+3)).^2+(y11-A(m+3)).^2); aver1=sum(l1)/2;
fr2w2(m)=aver1*(2162535*2/(1163+1188))/173115*(173115-(C(m+3))*(1-0.1)); f2w2(m)=fc2w2(m)+fr2w2(m)+3.25*200000*5+5000000; end
f2w=[f2w1 f2w2] fc2w=[fc2w1 fc2w2] fr2w=[fr2w1 fr2w2] x=8.20; y=6.00; for n=1:9
P(n)=sqrt((A(n)-x)^2+(B(n)-y).^2); end
aver0=sum(P)/9;
fc2=aver0*(4819569/1128)*230 x11=[3.60 6.90]; y11=[3.90 1.00];
l1=sqrt((x11-x).^2+(y11).^2); aver1=sum(l1)/2;
fr2=aver1*(2162535*2/(1163+1188)) f2=fr2+fc2+2.75*200000*5+5000000 运行结果如下： 城市 堪萨斯城 西雅图 洛杉矶 丹佛 达拉斯 芝加哥 亚特兰大 纽约 多伦多 蒙特利尔 运输成本 外向运输成本（美元） 内向运输成本（美元） 5年总运输成本（美元） 4.6465e+006 8.7276e+006 6.8953e+006 5.5011e+006 5.4818e+006 4.2562e+006 5.1217e+006 5.0157e+006 5.5366e+006 6.5341e+006 0.7056e+004 1.3150e+004 0.5613e+004 0.7002e+004 1.1009e+004 1.4717e+004 1.7929e+004 2.1171e+004 2.2598e+004 2.6183e+004 1.2404e+007 1.6991e+007 1.5151e+007 1.3758e+007 1.3743e+007 1.2521e+007 1.3390e+007 1.3287e+007 1.3809e+007 1.4810e+007 由上表数据可看出最佳选址地点仍是堪萨斯城，不需改变选址决策。 方法二：利用重心坐标公式计算
A1=[3.60 6.90 0.90 1.95 5.60 7.80 10.20 11.30 14.00 12.70 14.30 8.20]; B1=[3.90 1.00 9.10 4.20 6.10 3.60 6.90 3.95 6.55 7.80 8.25 6.00];
D1=[61500*(1+0.05) 120600*(1-0.1) 17000*(1+0.05) 32000*(1+0.05) 12500*(1+0.05) 9500*(1-0.1) 29500*(1-0.1) 21000*(1-0.1) 41300*(1-0.1) 8600*(1-0.1) 10700*(1-0.1) 0]; Cx21=sum(A1.*D1); Cy21=sum(B1.*D1); V2=sum(D1); Cx2=Cx21/V2 Cy2=Cy21/V2 for j=1:12
Q(j)=sqrt((A1(j)-Cx2)^2+(B1(j)-Cx2)^2); end Q
运行结果如下：重心坐标为（7.0587，4.0895） 城市 菲尼斯特 蒙特雷 西雅图 洛杉矶 丹佛 达拉斯 注：英里=距离*230
与重心的距离 4.6840 6.0608 6.4882 5.8542 1.7456 3.5373 城市 芝加哥 亚特兰大 纽约 多伦多 蒙特利尔 堪萨斯城 与重心的距离 3.1453 5.2586 6.9599 5.6898 7.3385 1.5567
具体位置如下图所示：
由上表数据可看出堪萨斯城仍是最佳选址地点，不需改变选址决策。 3.方法一：利用求加权平均值的方法算
A=[0.90 1.95 5.60 7.80 10.20 11.30 14.00 12.70 14.30]; B=[9.10 4.20 6.10 3.60 6.90 3.95 6.55 7.80 8.25];
C=[17000 32000 12500 9500 29500 21000 41300 8600 10700]; a=0; for i=1:9 for j=1:9 if j~=i
L(j)=sqrt((A(j)-A(i))^2+(B(j)-B(i))^2); else L(j)=0; end end
aver(i)=sum(L)/8; end
for k=1:9
fc3w(k)=aver(k)*0.0235*(1+0.25)*(182100-C(k))*230; x11=[3.60 6.90]; y11=[3.90 1.00];
l1=sqrt((x11-A(k)).^2+(y11-B(k)).^2); aver1=sum(l1)/2;
fr3w(k)=aver1*((2162535*2/(1163+1188))/182100)*(1+0.05)*(182100-C(k)); f3w(k)=fc3w(k)+fr3w(k)+3.25*200000*5+5000000;
end f3w x=8.20; y=6.00; for n=1:9
P(n)=sqrt((A(n)-x)^2+(B(n)-y)^2); end
aver0=sum(P)/9;
fc3=aver0*0.0235*(1+0.25)*230*182100; x11=[3.60 6.90]; y11=[3.90 1.00];
l1=sqrt((x11-x).^2+(y11-y).^2); aver1=sum(l1)/2;
fr3=aver1*(2162535*2/(1163+1188))*(1+0.15); f3=fr3+fc3+2.75*200000*5+5000000 运行结果如下： 运输成本 城市 堪萨斯城 西雅图 洛杉矶 丹佛 达拉斯 芝加哥 亚特兰大 纽约 多伦多 蒙特利尔 外向运输成本（美元） 0.77438e+006 1.1045e+007 0.8829e+007 0.6941e+007 0.6843e+007 0.5274e+007 0.6368e+007 0.6183e+007 0.6913e+007 0.8153e+007 内向运输成本（美元） 1.0814e+004 1.3958e+004 0.6028e+004 0.7409e+004 0.6373e+004 1.1339e+004 1.1106e+004 1.4745e+004 1.7335e+004 1.9918e+004 5年总运输成本（美元） 1.5505e+007 1.9309e+007 1.7085e+007 1.5198e+007 1.5099e+007 1.3536e+007 1.4629e+007 1.4447e+007 1.5181e+007 1.6423e+007
由以上数据可看出芝加哥是最佳选址地点，
所以需改变决策。
方法二：利用重心坐标公式计算
由于外向运输资率及内向运输资率的改变只是改变了运输成本，而与年需求量和各城市坐标都没影响，所以由第一小题的运算可知堪萨斯城为最佳选址地点，故不需改变决策。
五、模型的评价与推广
上述采用了求重心坐标的方法和求加权平均值的方法，具有一定的合理性，但这两种方法均是在理想状况下建立的模型，忽略了运输时天气、地理环境等因素的影响，在实际操作中各方面因素均需考虑。在我们所采用的方法中，重心法是一种模拟方法，它将物流系统中
的需求点和资源点看成是分布在某一平面范围内的物流系统，各点的需求量和资源量分别看成是物体的重量，物体系统的重心作为物流网点的最佳设置点，利用求物体系统重心的方法来确定物流网点的位置。相关问题均可如此处理；求加权平均值的方法是一种更具普遍意义的方法，便于理解，实用性强。
参考文献
【9787810824347】 程永生 《物流系统分析》 清华大学出版社,北京交通大学出版社 2004年
【9787508426600】 刘卫国 《MAYTLAB程序设计教程》 中国水利水电出版社 2005年