- 1、本文档共28页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
利用LINGO软件优化数学模型 邓 磊 西南大学 数学与统计学院 2011年4月15日 提 纲 问题1 合建水厂问题的优化 问题2 选址问题的优化 问题1 合建水厂问题的优化 在河边(X轴)有两座城市A、B, 其坐标分别是A(1,5), B(6,4), 现要合建一个自来水厂, 向它们供水.试设计自来水厂的位置以及铺设水管的方法, 使所铺设的水管总长度为最短. 这个问题实际上也是一个选址问题, 通过对几种选址方案的比较而得出最佳方案. 1. 特殊点法建模 图1 特殊点法 首先考虑在几个特殊位置建立水厂,如图1,(1)在离A最近的河边A1处建水厂, 水管从A1分别铺设到A、B, 得出水管总长度11.40; (2)在离B最近的B1处建水厂, 水管从B1分别铺设到A、B, 得出水管总长度11.07; (3) 仍在离A最近的A1处建水厂, 但B厂的水管经过A, 得出水管总长度10.40; (4)仍在离B最近的B1处建水厂, 但到A厂的水管经过B, 得出水管总长度9.10. 上述四个模型中,水管总长度每后一个比前一个小, 同时表明水管总长与水厂位置及水管铺设方法均有关. 2. 对称点连线法建模 如果水管必须分别从水厂铺设, 则1中的(3)、(4)不合要求, 但(1)、(2)中的铺设方法也不是最佳的. 利用光行最速原理求解:如图2,先作点B关于X轴的对称点B(6,-4), 再作线段AB交X轴于点W, 以点W为自来水厂的位置, 水管按折线AWB铺设成V型, 所用水管总长为10.30, 该结果优于(1)、(2)、(3)方法.但还不如(4)中方法好. 应该说,许多中学数学教师认为,2中的模型已经是最优了,未参加过数学建模培训的大学生基本上都是给出这个答案. 3. 利用费尔马点建模 在2中W点建立水厂, 改进水管铺设方法. 图3 费尔马点法 当然, 也可在W点取水, 而在N点建立水厂进行净化处理后再分别输送到A、B, 水管总长度不变. 4. 改进费尔马点建模 即便是参加过数学建模培训的大学生,通常也会认为3中的解就是最优解了. 但只要注意到NW并不与x轴垂直就会想到还可通过改变水厂位置进一步优化.移动W点使NW与x轴垂直, 则水管总长还会变小, 但随着W点的移动,相应地N又不是新位置的△ABW的费尔马点了, 进一步改变N点的位置, 又改变W点的位置, …, 如此下去, 何时才是最优解呢? 如图4, 设水厂的位置建在点M处, 显然M应在A1与B1之间, 设N(x,y)为△ABM的费尔马点, 则只有当MN⊥X轴时, 才可能使AN+BN+MN达到最短, 于是设点M的坐标为M(x,0).过点N作直线PQ∥X轴交AA1于P、交BB1于Q. 由费尔马点性质可知, ANM=∠BNM= 120度, 从而∠ANP=∠BNQ=30度, 由 列出方程组得 求解得 由此得N(4.37,3.06), M(4.36, 0), 这时的水管总长度为8.83. 5. 利用Lingo求解 利用Lingo建模求解, 设N(x,y),M(m,0), 则只需求x,y,m使NA+NB+NM最小, 将AB的坐标代入, 得如下程序. MODEL: MIN=@SQRT((x-1)^2+(y-5)^2)+@SQRT((x-6)^2+(y-4)^2)+@SQRT((x-m)^2+(y-0)^2); END 运行得出的结果与4一致, 表明4中的结果已经是优化到最佳状态了. 6. 优化效果比较 1. 料场位置确定时的最少运量 设第i个工地的坐标为(ai,bi), 水泥用量为di, 第j个两个临时料场的坐标为(xj,yj), 水泥日储量为ej, 从第j个料场运到第i个工地的水泥为c(i,j), 则有如下数学模型 目标函数: c(i,j)为决策变量, 可编制如下程序: MODEL: SETS: demand/1..6/:a,b,d; supply/1..2/:x,y,e; link(demand,supply):c; ENDSETS DATA: 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; e=20,20; !
文档评论(0)