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

一、问题提出 某企业有6个建筑工地要开工，每个工地位置（用平面坐标系(a，b)表示，距离单位：km）及水泥日用量d（吨）由下表给出。现在有两个料场在A(5,1)，B(2,7)，日储量各有20吨。 试制订天天供给计划，即从A，B两料场分别向各工地运输多少水泥，可使运输费用（总吨千米数）最小，并求出吨千米数。 ( 注：先画图，在坐标上标出各工地位置（用蓝色*标示）和料场位置（用红色o标示）) （2）现在企业准备建立两个新料场，日储量各为20吨，为使运输费用最省，问新料场应建在何处，并算出两料场分别向工地运输多少吨水泥和费用。 （注：初始值取x0=[3 5 4 7 1 0 0 0 0 0 5 11 5 4 7 7]’） 二、问题分析 对于问题（1），确定用A，B两料场分别向各工地运输水泥，使运输费用（总吨千米数）最小，即要知道两点间线段最小，料场到工地路线是直，而要满足六个工地需求，又要考虑到A、B两个料场供给量，即在各工地用量必需满足和各料场运输量不超出日储量条件下，使总吨千米数最小，这是线性问题。。 对于问题（2），需要重新改建六个新料场，使得在在各工地用量必需满足和各料场运输量不超出日储量条件下，使总吨千米数最小，则需要确定新料场具体位置，这是非线性问题。 三、模型假设 1、假设料场和建筑工地之间全部能够由直线抵达； 2、运输费用由“吨千米数”来衡量； 3、两料场日存放量够向各建筑工地供给； 4、运输途中不发生意外，从料场运出水泥总量不会超出各个料场日存放量。 四、模型建立 （显示模型函数结构过程） 记工地位置为，水泥日用量为，i=1,…,6;料场位置为， 日储量为，j=1,2；料场向工地运输量为。 目标函数为： 约束条件为： 当用临时料场时决议变量为： 当不用临时料场时决议变量为：，， 使用临时料场情形： 使用两个临时料场A(5,1)，B(2,7)．求从料场j向工地运输量. 在各工地用量必需满足和各料场运输量不超出日储量条件下，使总吨千米数最小，这是线性计划问题。线性计划模型为： 其中 ，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 改建两个新料场情形： 改建两个新料场，要同时确定料场位置(xj,yj)和运输量，在一样条件下使总吨千米数最小．这是非线性计划问题．非线性计划模型为： 设 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 五、模型求解 （显示模型求解方法、步骤及运算程序、结果） 建立chengxu.m程序： x=[1.25 8.75 0.5 5.75 3 7.25]; y=[1.25 0.75 4.75 5 6.5 7.75]; x0=[5,2]; y0=[1,7]; plot(x,y,*b); hold on; plot(x0,y0,or); text(1.25,1.25,1¤μ?1); text(8.75,0.75,1¤μ?2); text(0.5,4.75,1¤μ?3) text(5.75,5,1¤μ?4); text(3,6.5,1¤μ?5); text(7.25,7.25,1¤μ?6) text(5,1,á?3?A); text(2,7,á?3?B); 使用临时料场情形： 编写程序liaochang1.m以下： 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); 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 0