改建两个新料场的情形：
改建两个新料场，要同时确定料场的位置(xj,yj)和运送量 Xij，在同样条件下使总吨千米数
最小．这是非线性规划问题．非线性规划模型为：
26
min f
X ij (x j ai )2 ( y j bi )2
j1 i1
2
s.t. Xij di , i 1, 2,, 6 j 1
五、模型求解
（显示模型的求解方法、步骤及运算程序、结果）
(1)先画图，在坐标上标出各工地位置（用蓝色*标示）和料场位置（用红色 o 标示）
程序代码： x=[1.25 8.75 0.5 5.75 3 7.25]; y=[1.25 0.75 4.75 5 6.5 7.25]; x0=[5,2]; y0=[1,7]; plot(x,y,'*b'); hold on; plot(x0,y0,'or'); text(1.25,1.25,'¹¤µØ1'); text(8.75,0.75,'¹¤µØ2'); text(0.5,4.75,'¹¤µØ3'); text(5.75,5,'¹¤µØ4'); text(3.6,5,'¹¤µØ5'); text(7.25,7.25,'¹¤µØ6'); text(5,1,'Áϳ¡A'); text(2,7,'Áϳ¡B');
四、模型建立
（显示模型函数的构造过程）
记工地的位置为(ai，bi)，水泥日用量为 di，i=1,…,6;料场位置为(xj，yj)，日储量为 ej，j=1,2； 料场 j 向工地 i 的运送量为 Xij．
26
目标函数为： min f
X ij (x j ai )2 ( y j bi )2
j1 i1
26
min f aa(i, j) X ij j1 i1 2
s.t. X ij di , i 1,2,,6 j 1 6 X ij e j , j 1,2 i 1
其中 aa(i, j) (x j ai )2 ( y j bi )2 ，i=1,2,…,6,j=1,2,为常数.
设 X11=X1, X21= X 2,, X31= X 3, X41= X 4, X51= X 5,, X61= X 6 X12= X 7, X22= X 8,, X32= X 9, X42= X 10, X52= X 11,, X62= X 12
（2）目前公司准备建立两个新的料场，日储量各为 20 吨，为使运输费用最省，问新的料场 应建在何处，并算出两料场分别向工地运输多少吨水泥和费用。
（注：初始值取 x0=[3 5 4 7 1 0 0 0 0 0 5 11 5 4 7 7]’）
二、问题分析
本问题主要讨论并解决了某公司每天给工地的供应计划与临时料场选址的相关问题。工地的水 泥全部由水泥料场来提供，各个工地的水泥日用量与每个水泥料场的日储量是一定的。各个工 地所需水泥量可以有两个料场同时提供，也可以由一个料场全部提供。这就需要我们最优规划 来解决。为使总吨千米数达到最小，在考虑有直线道路连通的情况下建立相应的数学模型，给 出相关算法，并运用 Matlab 软件编程和处理相关数据，得到最优决策方案。 符号说明：
(2)使用两个临时料场的情形：
程序代码：
clear a=[1.25 8.75 0.5 5.75 3 7.25]; b=[1.25 0.75 4.75 5 6.5 7.75]; d=[3 5 4 7 6 11]; x=[5 2]; y=[1 7]; e=[20 20]; for i=1:6
for j=1:2 aa(i,j)=sqrt((x(j)-a(i))^2+(y(j)-b(i))^2);
i ：第 i 个建筑工地； j ：第 j 个临时料场；
di ：工地 i 的水泥日用量；
c ：料场 j 到工地 i 的水泥运输量； ij
aa(i, j) ：料场 j 到工地 i 的距离；
e ：料场 j 的日储量； j
三、模型假设
（1）假设料场和建筑工地之间都可以由直线到达； （2）运输费用由“吨千米数”来衡量； （3）两料场的日存储量够向各建筑工地供应； （4）运输途中不发生意外，从料场运出的水泥总量不会超过各个料场的日存储量；
2
X ij di ,
约束条件为： j1 6 X ij e j , i 1
i 1,2,,6 j 1,2
当用临时料场时决策变量为：Xij， 当不用临时料场时决策变量为：Xij，xj，yj．
使用两个临时料场的情形：
使用两个临时料场 A(5,1)，B(2,7)．求从料场 j 向工地 i 的运送量 Xij . 在各工地用量必须满足和 各料场运送量不超过日储量的条件下，使总的吨千米数最小，这是线性规划问题． 线性规划 模型为：
1
2
3
4
5
6
a
1.25
8.75
0.5
5.75

3
7.25
b
1.25
0.75
4.75
5
6.5
7.75
d
3
5
4
7
6
11
（1）试制定每天的供应计划，即从 A，B 两料场分别向各工地运送多少水泥，可使运输费用（总 的吨千米数）最小，并求出吨千米数。 ( 注：先画图，在坐标上标出各工地位置（用蓝色*标示）和料场位置（用红色 o 标示）)
课程名称 上机项目 专业班级
佛山科学技术学院
上机报告
数学建模 料场问题
一、 问题提出
某公司有 6 个建筑工地要开工，每个工地的位置（用平面坐标系(a，b)表示，距离单位： km）及水泥日用量 d（吨）由下表给出。目前有两个料场位于 A(5,1)，B(2,7)，日储量各有 20 吨。
工地位置（a，b）及水泥日用量 d
end end CC=[aa(:,1);aa(:,2)]'; A=[1 1 1 1 1 1 0 0 0 0 0 0
0 0 0 0 0 0 1 1 1 1 1 1]; B=[20;20]; Aeq=[1 0 0 0 0 0 1 0 0 0 0 0
6
Xij ej , j 1, 2
i 1
设 X11=X1, X21= X 2, X31= X 3, X41= X 4, X51= X 5,, X61= X 6
X12= X 7, X22= X 8, X32= X 9, X42= X 10, X52= X 11, X62= X 12
x1=X13, y1=X14, x2=X15, y2=X16