logo资料库

数学建模中的供应与选址问题.doc

第1页 / 共6页
第2页 / 共6页
第3页 / 共6页
第4页 / 共6页
第5页 / 共6页
第6页 / 共6页
资料共6页,全文预览结束
供应与选址 摘要:本论文主要讨论并解决了某公司每天的供应计划与临时料场选址的相关问 题。为使总吨千米数达到最小,在考虑有直线道路连通的情况下建立相应 的数学模型,给出了相关算法。并运用 Lingo9.0 等软件编程和处理相关数 据,得到最优决策方案。 问题一是一个线性规划问题,我们首先建立单目标的优化模型,也即 模型一.借助 Lingo 软件得到了该公司每天向六个建筑工地运输水泥的供应 计划如表 1,从而可使得总的吨千米数最小. 问题二是在问题一的基础上建立一个非线性规划模型,保持供应计划 不变的情况下,改变临时料场的位置以使吨千米数进一步减少。用 lingo 求解可知当新建的临时料场位于 C(6,4),D(7,8)两位置时,节省的吨千米 数可达到 30 . 表 1 j 1 jc (A 料场) 2 jc (B 料场) A,B 两料场每天分别向各工地运送的水泥量 1 4 0 2 6 0 3 0 6 4 0 7 5 0 8 6 2 9 关键字:线性规划 吨千米数 供应计划 一、 问题重述 某公司有 6 个建筑工地要开工,每个工地的位置(用平面坐标系 a,b 表示, 距离单位:千米)及水泥日用量 d(吨)由下表给出. 目前有两个临时料场位于 A(5,1),B(2,7),日储量各有 30 吨. (1)试制定每天的供应计划,即从 A,B 两料场分别向各工地运送多少吨水 泥,使总的吨千米数最小? (2)为了进一步减少吨千米数,打算舍弃两个临时料场,改建两个新的, 日储量各为 20 吨,问应建在何处,节省的吨千米数会多大? a b d 1 1 1 4 2 8 0 6 3 0 4 6 4 5 6 7 5 3 6 8 6 8 7 11 二、基本假设 1、料场与工地之间有直线道路; 2、两料场供应量应与工厂日用量达到平衡; 3、改建后供应计划保持原计划不变; 4、每个工地的位置用平面坐标的形式表示; 三、基本符号说明 i :第i 个临时料场;
j :第 j 个建筑工地; jd :工地 j 的水泥日用量; ijc :料场i 到工地 j 的水泥运输量; ijr :料场i 到工地 j 的距离; ie :料场i 的日储量; 四、问题的分析,模型的建立及求解 4.1.1 问题一的分析 某公司有 6 个建筑工地,位置坐标为( ja , jb ) (单位:千米),水泥日用量 jd (单位:吨) 现有 A(5,1),B(2,7) 两料场,记( ix , iy ),i =1,2,日储量 ie 各有 30 吨. 已知每个工地的位置及水泥日用量如下表 a b d 1 1 1 4 2 8 0 6 3 0 4 6 4 5 6 7 5 3 6 8 6 8 7 11 4.1.2 模型一的建立 由题知,问题一是一个线性规划模型,确定分配量求最小值,即使 总的吨千米数最小. 则目标函数: min 其中 r  ij 约束条件: s.t 6 2  i 1  j 1  c ij  6  j 1  d j 6 2  i 1  ij ij 1  j c r ( a x  j 2 6  j c 1  1  i  e i ij 其中 ie 为 30 吨 4.1.3 模型一的求解 2 ) i  ( b y  j i 2 )
将已知数据代入模型中,用 lingo 软件求解(程序见附录 1),得到结果(程 序运行结果见附录 2)如下表: A,B 两料场每天分别向各工地运送的水泥量 j 1 jc (A 料场) 2 jc (B 料场) 1 4 0 4.2.1 问题二的分析 2 6 0 3 0 6 4 0 7 5 0 8 6 2 9 问题二是在问题一的基础上,进一步减少吨千米数,舍弃两个临时料场, 改建两个新的临时料场,日储量各为 20 吨,求新建的料场的位置,在其它条件 不变下使总吨公里数最小,此时节省的吨千米数最大.为此,需建立一个非线形 规划模型. 4.2.2 模型二的建立 问题二是一个非线性模型,求解取最小值时需满足的最优条件. 目标函数: min 6 2  i 1  j 1  c a x  ( ij j 2 ) i  ( b j  2 y ) i 约束条件: s.t 6 2  i 1  j 1  c ij  6  j 1  d j 2 6 i 1  1   j c ix  , 0  e i ij 其中 ie 为 20 吨 iy  0 4.2.2 模型二的求解 将模型一求得的供用计划数据代入模型二中,用 lingo(其程序见附录 3) 求解得(4,6)(7,8)此时节省的最大吨千米数最大 89.88349. 五、模型的评价 本文优点是建立了规划模型,通过lingo软件进行线性求解,得出各种供应 计划方案的最优解;不足之处,在处理供应计划与选址的关系上比较含糊,没有 深入讨论. 六、参考文献 [1] 姜启源、谢金星等,数学模型 ,北京:高等教育出版社. 2007.8.
[10]席少霖等,最优化计算方法,上海:上海科学技术出版社,2003. [8] 谢金星等,优化建模与 LINDO/LINGO 软件,北京:清华大学出版社,2005.7. 附录七 附录1 MODEL: Title Location Problem; sets: demand/1..6/:a,b,d; supply/1..2/:x,y,e; link(supply,demand):c; endsets data: !需求点的位置; a=1,8,0,5,3,8; b=1,0,4,6,6,7; ! 供需 量; d=4,6,6,7,8,11; e=30,30; enddata init: ! 初始点; x,y=5,1,2,7; endinit ! 目标函数; [OBJ] min= @sum(link(i,j): c(i,j)*((x(i)-a(j))^2+(y(i)-b(j))^2)^(1/2) ); ! 需求约束; @for(demand(j):[DEMAND_CON] @sum(supply(i):c(i,j)) =d(j);); ! 供应约束; @for(supply(i):[SUPPLY_CON] @sum(demand(j):c(i,j)) <=e(i); ); @for(supply: @bnd(0,X,8); @bnd(0,Y,7); ); END 附录 2 Feasible solution found. Objective value: Total solver iterations: Model Title: Location Problem 119.0558 53 Variable A( 1) A( 2) A( 3) A( 4) Value 1.000000 8.000000 0.000000 5.000000 Reduced Cost 0.000000 0.000000 0.000000 0.000000
A( 5) A( 6) B( 1) B( 2) B( 3) B( 4) B( 5) B( 6) D( 1) D( 2) D( 3) D( 4) D( 5) D( 6) X( 1) X( 2) Y( 1) Y( 2) E( 1) E( 2) C( 1, 1) C( 1, 2) C( 1, 3) C( 1, 4) C( 1, 5) C( 1, 6) C( 2, 1) C( 2, 2) C( 2, 3) C( 2, 4) C( 2, 5) C( 2, 6) 3.000000 8.000000 1.000000 0.000000 4.000000 6.000000 6.000000 7.000000 4.000000 6.000000 6.000000 7.000000 8.000000 11.00000 7.999999 4.999999 0.000000 6.000000 30.00000 30.00000 4.000000 6.000000 0.000000 0.000000 0.000000 2.000000 0.000000 0.000000 6.000000 7.000000 8.000000 9.000000 Row OBJ DEMAND_CON( 1) DEMAND_CON( 2) DEMAND_CON( 3) DEMAND_CON( 4) DEMAND_CON( 5) DEMAND_CON( 6) SUPPLY_CON( 1) SUPPLY_CON( 2) Slack or Surplus 119.0558 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 18.00000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 -2.040203 2.328466 -2.565686 5.838852 0.000000 0.000000 0.000000 0.000000 -0.2786167 2.870480 1.972527 0.000000 3.169781 10.54593 0.000000 0.000000 0.000000 0.000000 Dual Price -1.000000 -7.071066 -0.1727034E-05 -9.222887 -3.837723 -5.837722 -7.000000 0.000000 3.837722
附录 3 MODEL: Title Location Problem; sets: demand/1..6/:a,b,d; supply/1..2/:x,y,e; link(supply,demand):c; endsets data: !需求点的位置; a=1,8,0,5,3,8; b=1,0,4,6,6,7; ! 供需 量; d=4,6,6,7,8,11; c=4 6 0 0 0 2 0 0 6 7 8 9 ; e=20,20; enddata init: ! 初始点; x,y=5,1,2,7; endinit ! 目标函数; [OBJ] min= @sum(link(i,j): c(i,j)*((x(i)-a(j))^2+(y(i)-b(j))^2)^(1/2) ); ! 需求约束; @for(demand(j):[DEMAND_CON] @sum(supply(i):c(i,j)) =d(j);); ! 供应约束; @for(supply(i):[SUPPLY_CON] @sum(demand(j):c(i,j)) <=e(i); ); @for(supply: @bnd(0,x,8); @bnd(0,y,7); ); @for(supply: @bnd(0,x,8);@gin(x); @bnd(0,y,7); @gin(y);); END
分享到:
收藏