![数学建模作业5数学规划模型----供应与选址的问题.doc_第1页](http://file1.renrendoc.com/fileroot_temp2/2020-3/15/fabbc0b9-d3d4-4215-90ce-06c76465f4ae/fabbc0b9-d3d4-4215-90ce-06c76465f4ae1.gif)
![数学建模作业5数学规划模型----供应与选址的问题.doc_第2页](http://file1.renrendoc.com/fileroot_temp2/2020-3/15/fabbc0b9-d3d4-4215-90ce-06c76465f4ae/fabbc0b9-d3d4-4215-90ce-06c76465f4ae2.gif)
![数学建模作业5数学规划模型----供应与选址的问题.doc_第3页](http://file1.renrendoc.com/fileroot_temp2/2020-3/15/fabbc0b9-d3d4-4215-90ce-06c76465f4ae/fabbc0b9-d3d4-4215-90ce-06c76465f4ae3.gif)
![数学建模作业5数学规划模型----供应与选址的问题.doc_第4页](http://file1.renrendoc.com/fileroot_temp2/2020-3/15/fabbc0b9-d3d4-4215-90ce-06c76465f4ae/fabbc0b9-d3d4-4215-90ce-06c76465f4ae4.gif)
![数学建模作业5数学规划模型----供应与选址的问题.doc_第5页](http://file1.renrendoc.com/fileroot_temp2/2020-3/15/fabbc0b9-d3d4-4215-90ce-06c76465f4ae/fabbc0b9-d3d4-4215-90ce-06c76465f4ae5.gif)
已阅读5页,还剩7页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
一、问题提出某公司有6个建筑工地要开工,每个工地的位置(用平面坐标系(a,b)表示,距离单位:km)及水泥日用量d(吨)由下表给出。目前有两个料场位于A(5,1),B(2,7),日储量各有20吨。(1) 试制定每天的供应计划,即从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 6X12= 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); 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);使用临时料场的情形:编写程序liaochang1.m如下:cleara=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); endend 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 1 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 1 ;beq=d(1);d(2);d(3);d(4);d(5);d(6);vlb=0 0 0 0 0 0 0 0 0 0 0 0;vub=;x0=1 2 3 0 1 0 0 1 0 1 0 1;xx,fval=linprog(CC,A,B,Aeq,beq,vlb,vub,x0)程序截图如下:程序的运行结果为:xx = 3.0000 5.0000 0.0000 7.0000 0.0000 1.0000 0.0000 0.0000 4.0000 0.0000 6.0000 10.0000fval = 136.2275运行结果截图如下:即由料场A、B向6个工地运料方案为:123456料场1350701料场20040610总的吨千米数为136.2275.改建两个新料场的情形:先编写M文件liaochang.m:function f=liaoch(x)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;f1=0;for i=1:6 s(i)=sqrt(x(13)-a(i)2+(x(14)-b(i)2); f1=s(i)*x(i)+f1;endf2=0;for i=7:12 s(i)=sqrt(x(15)-a(i-6)2+(x(16)-b(i-6)2); f2=s(i)*x(i)+f2;end f=f1+f2;再编写主程序liaochang2.m为:clearx0=3 5 4 7 1 0 0 0 0 0 5 11 5 4 7 7;A=1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 0 0 0 0;B=20;20;Aeq=1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0;beq=3 5 4 7 6 11;vlb=zeros(12,1);-inf;-inf;-inf;-inf;vub=;x,fval,exitflag=fmincon(liaoch,x0,A,B,Aeq,beq,vlb,vub) 程序截图如下:程序运行结果如下:x = Columns 1 through 8 3.0000 5.0000 4.0000 7.0000 1.0000 0 0 0 Columns 9 through 16 0 0 5.0000 11.0000 5.6962 4.9289
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 【正版授权】 IEC 60636:1979 EN-FR Flexible waveguide assembly performance
- 【正版授权】 IEC 60601-2-46:2016 RLV EN Medical electrical equipment - Part 2-46: Particular requirements for the basic safety and essential performance of operating tables
- 【正版授权】 IEC 60588-3:1977 EN-FR Askarels for transformers and capacitors. Part 3: Specifications for new askarels
- 【正版授权】 IEC 60512-19-1:2010 EN-FR Connectors for electronic equipment - Tests and measurements - Part 19-1: Chemical resistance tests - Test 19a: Fluid resistance of pre-insulated crimp barre
- 【正版授权】 IEC 60371-3-7:1995/AMD1:2006 EN-D Amendment 1 - Insulating materials based on mica - Part 3: Specifications for individual materials - Sheet 7: Polyester film mica paper with an epoxy
- 签橱柜合同范本
- 租客非法分租合同范本
- 戒烟合同范本模板
- 传媒工作室与主播合同范本
- 【正版授权】 IEC 60335-2-84:2019 EN-FR Household and similar electrical appliances - Safety - Part 2-84: Particular requirements for toilet appliances
- 2024年工艺品定制合同
- 20230802基于SAP锂电池数字化转型解决方案参考
- 酒店客人投诉处理课件
- 2022重庆中考语文A卷试卷分析
- 肝胆脾胃科护士出科小结
- 学历(学籍)认证报告
- 四人股东合作协议书完整版doc(二篇)
- 多年采购经验总结汇报
- 2024年春季小学三年级下册语文部编版“快乐读书吧”课时练习
- 关节功能28天复原方案
- 铁路劳动安全教育课件
评论
0/150
提交评论