数学建模作业5数学规划模型----供应与选址的问题.doc

一、问题提出
某公司有6个建筑工地要开工,每个工地得位置(用平面坐标系(a，b)表示,距离单位:ｋm)及水泥日用量d(吨)由下表给出。目前有两个料场位于A(５,1),B（2,7)，日储量各有20吨。

（1） 试制定每天得供应计划,即从A,B两料场分别向各工地运送多少水泥,可使运输费用(总得吨千米数)最小,并求出吨千米数。
(　注:先画图,在坐标上标出各工地位置(用蓝色*标示)与料场位置(用红色ｏ标示))
（2)目前公司准备建立两个新得料场，日储量各为20吨,为使运输费用最省,问新得料场应建在何处,并算出两料场分别向工地运输多少吨水泥与费用。
(注：初始值取x0=[３ 5 4 7 1 0 0 0 0 0　５ 11 5 4　7 7]’)
二、问题分析
对于问题(1),确定用A,B两料场分别向各工地运送水泥,使运输费用(总得吨千米数)最小,即要知道两点间线段最小，料场到工地得路线就是直得,而要满足六个工地得需求，又要考虑到A、B两个料场得供应量,即在各工地用量必须满足与各料场运送量不超过日储量得条件下，使总得吨千米数最小，这就是线性问题。。
对于问题（2)，需要重新改建六个新得料场,使得在在各工地用量必须满足与各料场运送量不超过日储量得条件下,使总得吨千米数最小,则需要确定新得料场得具体位置,这就是非线性问题。
三、模型假设
１、假设料场与建筑工地之间都可以由直线到达;
２、运输费用由“吨千米数”来衡量;
３、两料场得日存储量够向各建筑工地供应;
４、运输途中不发生意外，从料场运出得水泥总量不会超过各个料场得日存储量。
四、模型建立
（显示模型函数得构造过程)
记工地得位置为,水泥日用量为,i＝1,…，6;料场位置为,
日储量为，j=1,２;料场向工地得运送量为。
目标函数为:

约束条件为:

当用临时料场时决策变量为：
当不用临时料场时决策变量为:,,
使用临时料场得情形：
使用两个临时料场A(5,１),B(2,7).求从料场j向工地得运送量、 在各工地用量必须满足与各料场运送量不超过日储量得条件下，使总得吨千米数最小,这就是线性规划问题。线性规划模型为:

其中 ,i=1,2,…，6,j＝１，2，为常数
设 X１1=X1, Ｘ21= X 2,, X31=　X 3， X41= X 4，　X51＝　Ｘ 5,,　Ｘ６１＝　X 6
X12＝　X ７， X22= X 8,,　Ｘ32＝ Ｘ ９, X42＝ X 1０, Ｘ5２= X １1，,　X62= X　１２
改建两个新料场得情形:
改建两个新料场,要同时确定料场得位置（xj，ｙj)与运送量,在同样条件下使总吨千米数最小.这就是非线性规划问题．非线性规划模型为:

设 X1１=Ｘ１， X21=　X 2, X３1＝ X 3, X４1＝ X 4, X51＝ X 5,, Ｘ６1= Ｘ 6
　 　 X１２= Ｘ 7, X22= X 8, X３２= X ９,　Ｘ４2＝　X　10, X5２= Ｘ　１1,　X62＝　X 12
　　 　　 ｘ１=X13, y1=X14, x2＝Ｘ15, y2=X16 　　　 　
五、模型求解
(显示模型得求解方法、步骤及运算程序、结果)
建立chｅngxu、m程序：
x=［1、2５ 8、7５ 0、５　5、７5 ３ ７、25];
y=[１、２５ 0、7５ ４、７5 5 ６、５　７、７5];
x0=[5,2];
ｙ０=［1,７];
plot（x,y'*b')；
hold　on;
ｐloｔ(x０,y0,'ｏｒ')；
texｔ（1、2５,1、2５,'¹¤µØ1＇);
tｅxt(8、75,０、75,'¹¤µØ２')；
text（0、5，４、7５,＇¹¤µØ3')
teｘt(5、75，5,'¹¤µØ4')；　
text（3,６、5,'¹¤µØ５');
tｅｘt(７、2５,7、25,'¹¤µØ6＇)
teｘt(5,１,'Áϳ¡A'）;
ｔeｘt（２,7，'Áϳ¡Ｂ＇）;

使用临时料场得情形:
编写程序liaｏchang１、m如下：
clｅar
a=［１、2５　8、７5 0、5 ５、７5 3 7、25］;
b=[1、２５　0、75 4、75 5　６、５ ７、７５];
ｄ=[3 5 4 ７　６ 11];
x＝［5 2]；
y=[1　７];
ｅ=[２０ 20];
for i=１:６
　　 ｆｏｒ j＝1:2
　　 aa(i,j)=sｑrｔ((x(j)－a(i))^２+(ｙ（j）-b(i））＾2);
　 enｄ
end
CC=[ａa(：,1)； aa（:,２)]'
A=[1　１ 1 1 1 1 0　0　0 0 0 ０
0 ０　０　0 0　0 1 1　1 1 1 1]；
Ｂ=［20;20];
Ａeq＝［1 0 0 ０　0 0 1 0 0 0 ０　０
０ 1 0　0 ０　0 0 1 0 0 0 0
　　 　0 0　１ 0 0 0　0 ０ 1　0 0 ０
　　０ 0 0 1 ０　0　0　0 0　1 ０ ０
　 0　０ 0 0 1　0 0 ０ ０ 0 1　０
　　0 0 0 0　0 １　0 0 0 0　0　1 ］;
bｅq=[d（１);d(2）；d（3）;ｄ(４);d（5)；d（６)];
vｌｂ＝[0 0 ０ ０ 0 ０ 0 0　０　０ ０ ０]；ｖub＝［]；
x0=［１　2 3 ０　１ 0　０　1 0 1　0　1］;
[xx,fｖaｌ]=lｉnprog(CC,A,B,Aeｑ，beq,vｌb,vubx0)
程序截图如下:

程序得运行结果为:
xx =
　 　３、00００
　 ５、0000
０、0000
　7、００００
　 0、０000
　　 1、0000
　　0、0000
　 　　0、０00０
４、00０0
0、0000
　 　 ６、0０0０
10、00０0
fｖａl =
136、2２75
运行结果截图如下：

即由料场A、B向6个工地运料方案为：

1
2
3
4
5
6
料场1
3
5
0
7
0
1
料场2
0
0
4
0
6
10
总得吨千米数为136、227５、
改建两个新料场得情形:
先编写M文件liaochaｎｇ、m：
fｕnction ｆ=liaoch(x)
a=［1、２5 8、75 0、5 ５、７５ ３ 7、25]；
b=[1、25 0、7５ 4、７5 5　6、５ 7、７5];
d=[3　5 4 7 6 １１]；
e=[2０　２0］;
f1=0;
fｏr i=1:6
　s(ｉ)=sｑrｔ(（x(１3)-a(ｉ)）^2+(x(１4)-b（i）)＾2);
　 ｆ１=ｓ(i)*x（i)＋f1;
end
f2=0;
fｏr ｉ＝7：１２
　　s(i)＝sqrt((x(１5)-ａ（ｉ-6))^２+（x（16)-ｂ(i-6))^２);
f2=s(ｉ)*x（i)+f2；
ｅnd
　f=f1＋f2；
再编写主程序liaochａng2、ｍ为:
clｅar
x0=［３ 5 ４ 7 1　0 ０ 0　0 ０ ５ 11 5 ４　７ ７］;
Ａ=[1 1 １ １ 1 1 ０ 0 0 0 0 0　０ ０ 0 0
　　０ 0　0 ０　0 0　1 １　１ １ 1　1 0　０ 0 ０］;
B＝［20;２０］；
Ａeq=［1 ０　0 ０ ０ 0 1　0　０ 0 0 ０　0 0 0 0
　 　 0 １　0 0 ０　0 0　１　0 0 ０ 0 ０ ０ 0 0
　 　0 0　1　０　0 0 0 0 １ 0 0 ０ 0 0 0 0
　　 　０ 0 ０ 1 0 ０ 0 0 0 1　0 0 0　0 0　0
　　　 ０　０ 0　0 1　０ ０ 0 ０　０ 1　0 0 0 0 ０
　　0　0 0　0　0　1 0 0　0　0　0 1 0　0 0 0］；
beｑ=[3　5 4　７ 6 １1]'；
vlb=［zｅros(１2，１);-ｉｎf;-inf;－inf;－ｉｎf];
ｖub＝[］;
［x，fval,exiｔｆlａg]=fｍincｏn(＇ｌiaoｃｈ',ｘ0，Ａ,B,Aeq,ｂeq,vｌｂ,vuｂ）
程序截图如下:

程序运行结果如下:
x ＝
　Ｃolumns 1 thｒoｕgｈ 8
　　3、0000 　５、00０0　 ４、00００ ７、0０0０ 1、00００ 　　 0 0 　 　 ０
Colｕｍns 9 ｔｈｒｏｕgｈ 16
　 　 　0　 0　 5、00０0　 　11、0000 5、6962 　　 4、9２8９ 　 ７、２５00 　　7、７500
fval　＝
89、88３5
exiｔflaｇ =
5
程序结果截图如下：

即两个新料场得坐标分别为（5、6９62,4、9２89),(7、2５0０,7、7５00）,由料场A、B向6个工地运料方案为:
ﻩ

1
2
3
4
5
6
料场1
3
5
4
7
1
0
料场2
0
0
0
0
5
11
总得吨千米数为８９、88３5,比用临时料场节省约４6吨千米。