大学生数学建模:规划的基础实验-答案_第1页
大学生数学建模:规划的基础实验-答案_第2页
大学生数学建模:规划的基础实验-答案_第3页
大学生数学建模:规划的基础实验-答案_第4页
大学生数学建模:规划的基础实验-答案_第5页
已阅读5页,还剩9页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

实验作业——规划基础练习格式要求——写出:程序、结果、解释、进一步(一)线性规划问题1.用matlab及lingo求解下列线性规划问题:程序:lindoMax3x1-x2-x3StX1-2x2+x3<11-4x1+x2+2x3>3-2x1+x3=1Endmax=3*x1-x2-x3

;x1-2*x2+x3<=11

;-4*x1+x2+2*x3>=3

;-2*x1+x3=1

;结果:z=4,x1=4,x2=1,x3=92.某班男同学30人、女同学20人,植树。工作效率(个/人、天)如下表。如何安排,植树最多?解:假设分别有x1、x2、x3个男同学挖坑、栽树、浇水,y1、y2、y3个女同学挖坑、栽树、浇水。Maxf=20x1+10y1程序MAX20x1+10y1STx1+x2+x3=30y1+y2+y3=2020x1+10y1-30x2-20y2=030x2+20y2-25x3-15y3=0ENDGIN6所求最优解为fmax=340棵,X1=13(男13人全天挖坑),X2=4(男4人全天栽树),X3=13(男13人全天浇水);Y1=8,(女8人全天挖坑),Y2=11(女11人全天栽树),Y3=1(女1人全天浇水)其实可以取消整数的限制MAX20x1+10y1STx1+x2+x3=30y1+y2+y3=2020x1+10y1-30x2-20y2=030x2+20y2-25x3-15y3=0END所求最优解为fmax=350棵,X1=35/2(男17人全天挖坑,1个人挖半天坑)X2=0,X3=25/2(男12人全天浇水,1个人浇水半天);Y1=0,Y2=35/2(女17人全天栽树,1个人栽树半天),Y3=5/2(女2人全天浇水,1个人浇水半天)。结果更优,一半可以解释为一天中的一半时间若用lingo语言编程,则程序如下:model:sets:nan/wa,zai,jiao/:nanxl,x;nv/wa,zai,jiao/:nvxl,y;endsetsdata:nanxl=20,30,25;nvxl=10,20,15;nannum=30;nvnum=20;enddatamax=nanxl(1)*x(1)+nvxl(1)*y(1);@sum(nan:x)=nannum;@sum(nv:y)=nvnum;!@for(nan:@gin(x));!@for(nv:@gin(y));nanxl(1)*x(1)+nvxl(1)*y(1)=nanxl(2)*x(2)+nvxl(2)*y(2);nanxl(1)*x(1)+nvxl(1)*y(1)=nanxl(3)*x(3)+nvxl(3)*y(3);end3.某牧场饲养一批动物,平均每头动物至少需要700g蛋白质、30g矿物质和100g维生素。现有A、B、C、D、E五种饲料可供选用,每千克饲料的营养成分(单位:g)与价格(单位:元/kg)如下表所示:蛋白质矿物质维生素价格A31.00.50.4B20.51.01.4C10.21.20.8D62.02.01.6E120.50.81.6试求能满足动物生长营养需求又最经济的选用饲料方案。模型:变量命名—设A饲料选用X1kg,B饲料选用X2kg,C饲料选用X3kg,D饲料选用X4kg,E饲料选用X5kg.目标函数:设买饲料的总成本为Z。MinZ=0.4X1+1.4X2+0.8X3+1.6X4+1.6X5.约束条件:蛋白质:3X1+2X2+X3+6X4+12X5>700;矿物质:X1+0.5X2+0.2X3+2X4+0.5X5>30;维生素:0.5X1+X2+1.2X3+2X4+0.8X5>100;X1,X2,X3,X4,X5均为非负。程序:MIN0.4X1+1.4X2+0.8X3+1.6X4+1.6X5ST3X1+2X2+X3+6X4+12X5>700X1+0.5X2+0.2X3+2X4+0.5X5>300.5X1+X2+1.2X3+2X4+0.8X5>100END解释:经济即总成本最小,满足需求为约束条件,变量为非负约束93.33334VARIABLEVALUEREDUCEDCOSTX1233.3333280.000000X20.0000001.133333X30.0000000.666667X40.0000000.800000X50.0000000.000000ROWSLACKORSURPLUSDUALPRICES2)0.000000-0.1333333)203.3333280.0000004)16.6666660.000000进一步:三个需求的约束条件在报告中显示影子价格分别为-0.133333,0,0,则说明调整需求量对优化结果影响很小。Lingo程序model:sets:yingyang/dbz,kwz,wss/:xuqiu;ciliao/a,b,c,d,e/:jiage,x;chengfen(ciliao,yingyang):hanliang;endsetsdata:xuqiu=700,30,100;jiage=0.4,1.4,0.8,1.6,1.6;hanliang=3 1.0 0.52 0.5 1.01 0.2 1.26 2.0 2.012 0.5 0.8;enddatamin=@sum(ciliao:jiage*x);@for(yingyang(i):@sum(ciliao(j):hanliang(j,i)*x(j))>=xuqiu(i));end4.在以色列,为分享农业技术服务和协调农业生产,常常由几个农庄组成一个公共农业社区。在本课题中的这个公共农业社区由三个农庄组成,我们称之为南方农庄联盟。南方农庄联盟的全部种植计划都由技术协调办公室制订。当前,该办公室正在制订来年的农业生产计划。南方农庄联盟的农业收成受到两种资源的制约。一是可灌溉土地的面积,二是灌溉用水量。这些数据由下表给出。表1南方农庄联盟可用资源数据农庄可灌溉耕地面积(英亩)分配的用水量(英亩-英尺)ABC400600300600800375注:英亩-英尺是水容积单位,1英亩-英尺就是面积为1英亩,深度为1英尺的体积;1英亩-英尺≈1233.48立方米。南方农庄联盟种植的作物是甜菜、棉花和高粱,这三种作物的纯利润及耗水量不同。农业管理部门根据本地区资源的具体情况,对本联盟农田种植规划制定的最高限额数据由下表给出。表2南方农庄联盟农作物种植数据农作物最大种植面积(英亩)耗水量(英亩-英尺/英亩)纯利润(美元/英亩)甜菜棉花高粱600500325321400300100三家农庄达成协议:各家农庄的播种面积与其可灌溉耕地面积之比相等;各家农庄种植何种作物并无限制。所以,技术协调办公室面对的任务是:根据现有的条件,制定适当的种植计划帮助南方农庄联盟获得最大的总利润,现请你替技术协调办公室完成这一决策。对于技术协调办公室的上述安排,你觉得有何缺陷,请提出建议并制定新的种植计划。变量:设A农庄种植甜菜面积为X11,种植棉花面积为X12,种植高粱面积为X13.B农庄种植甜菜面积X21,种植棉花面积为X22,种植高粱面积为X23.C农庄种植甜菜面积为X31,种植棉花面积为X32,种植高粱面积为X33.模型&程序:max400x11+400x21+400x31+300x12+300x22+300x32+100x13+100x23+100x33st3x11+2x12+x13<6003x21+2x22+x23<8003x31+2x32+x33<375x11+x12+x13<400x21+x22+x23<600x31+x32+x33<300x11+x21+x31<600x12+x22+x32<500x13+x23+x33<325400x21+400x22+400x23-600x11-600x12-600x13=0300x21+300x22+300x23-600x31-600x32-600x33=0endLingo程序为:model:sets:nzh/A,B,C/:kymj,kysl;nzw/tian,mian,gao/:maxmj,hsl,clr;renwu(nzh,nzw):x;endsetsdata:kymj=400,600,300;kysl=600,800,375;maxmj=600,500,325;hsl=3,2,1;clr=400,300,100;enddata!目标是纯利润最大化;max=@sum(renwu(I,J):clr(J)*x(I,J));!可灌溉面积约束;@for(nzh(I): @sum(nzw(J):x(I,J))<=kymj(I));!可用水量约束;@for(nzh(I): @sum(nzw(J):hsl(J)*x(I,J))<=kysl(I));!最大种植面积约束;@for(nzw(J): @sum(nzh(I):x(I,J))<=maxmj(J));!播种面积与可灌溉面积之比相等的约束;@sum(nzw(J):x(1,J))*kymj(2)=@sum(nzw(J):x(2,J))*kymj(1);@sum(nzw(J):x(2,J))*kymj(3)=@sum(nzw(J):x(3,J))*kymj(2);end可见运行结果中,约束条件1和约束条件3中的值对优化结果影响很大,即分配给A和C农场的用水量不合理,需重新分配,故设A用水量为a,B用水量为b,C用水量为c,得程序:max400x11+400x21+400x31+300x12+300x22+300x32+100x13+100x23+100x33st3x11+2x12+x13-a<03x21+2x22+x23-b<03x31+2x32+x33-c<0a+b+c=1775x11+x12+x13<400x21+x22+x23<600x31+x32+x33<300400x21+400x22+400x23-600x11+600x12+600x13=0300x21+300x22+300x23-600x31+600x32+600x33=0x11+x21+x31<600x12+x22+x32<500x13+x23+x33<325end结果:比原方案更优。5.有一艘货轮,分前、中、后三个舱位,它们的容积与最大允许载重量如下表所示:前舱中舱后舱最大允许载重量(t)200030001000容积(m3)400054001000现有三种货物待运,已知有关数据如下表所示:商品数量(件)每件体积(m3/件)每件重量(t/件)运价(元/件)A6001081000B100056700C80075600又为了航运安全,要求前、中、后舱在实际载重量上大体保持各舱最大允许载重量的比例关系。具体要求前、后舱分别与中舱之间载重理比例上偏差不超过15%,前、后舱之间不超过10%。问该货轮应装载A,B,C各多少件,其运费收入为最大?解:设前舱运AX11件,运BX12件,运CX13,中舱运AX21件,运BX22件,运C23件,后舱运AX31件,运BX32件,运CX33件。程序max1000x11+1000x12+1000x13+700x21+700x22+700x23+600x31+600x32+600x33stx11+x12+x13<600x21+x22+x23<1000x31+x32+x33<8008x11+6x21+5x31<20008x12+6x22+5x32<30008x13+6x23+5x33<100010x11+5x21+7x31<400010x12+5x22+7x32<540010x13+5x23+7x33<1000x12+x22+x32-3x13-3x23-3x33<450x12+x22+x32-3x13-3x23-3x33>-450x11+x21+x31-2x13-2x23-2x33<200x11+x21+x31-2x13-2x23-2x33>-200end结果OBJECTIVEFUNCTIONVALUE1)742000.0VARIABLEVALUEREDUCEDCOSTX11250.0000000.000000X12300.0000000.000000X1350.0000000.000000X210.00000020.000000X220.00000020.000000X23100.0000000.000000X310.0000000.000000X32120.0000000.000000X330.0000006.000000ROWSLACKORSURPLUSDUALPRICES2)0.00000040.0000003)900.0000000.0000004)680.0000000.0000005)0.000000120.0000006)0.000000120.0000007)0.000000110.0000008)1500.0000000.0000009)1560.0000000.00000010)0.0000008.00000011)480.0000000.00000012)420.0000000.00000013)250.0000000.00000014)150.0000000.000000这个线性规划的最优解为x11=250,x12=300,x13=50,x23=100,x32=120,最优值为z=742000即前舱运A250件,中仓运A300件,运B100件后,运C120件,后舱运A50件,运费收入最大为742000。注意:由于此问题不加整数限制时已经得到整数最优解,故不用整数规划。6.某战略轰炸机群奉命摧毁敌人军事目标。已知该目标有四个要害部位,只要摧毁其中之一即可达到目的。为完成此项任务的汽油消耗量限制为48000升、重型炸弹48枚、轻型炸弹32枚。飞机携带重型炸弹时每升汽油可飞行2千米,带轻型炸弹时每升汽油可飞行3千米。又知每架飞机每次只能装载一枚炸弹,每出发轰炸一次除来回路程汽油消耗(空载时每升汽油可飞行4千米)外,起飞和降落每次各消耗100升。有关数据如表所示。为了使摧毁敌方军事目标的可能性最大,应如何确定飞机轰炸的方案,要求建立这个问题的线性规划模型,并用lingo求解。解:max=1-(0.9^x11)*(0.92^x12)*(0.8^x21)*(0.84^x22)*(0.85^x31)*(0.88^x32)*(0.75^x41)*(0.8^x42);537.5*x11+687.5*x12+560*x21+480*x22+605*x31+515*x32+650*x41+550*x42<=48000;x11+x21+x31+x41<=48;x12+x22+x32+x42<=32;@gin(x11);@gin(x12);@gin(x21);@gin(x22);@gin(x31);@gin(x32);@gin(x41);@gin(x42);改进:讲目标函数的乘变成对数化后的加能避免计算误差带来的解的计算问题model:sets: yaohai/1234/:jl; danyao/qingzhong/:yh,ssl; fangan(yaohai,danyao):logpp,pp,cishu;!储存每个要害,每种弹药打击可能性与次数;endsetsmax=@sum(yaohai(I): @sum(danyao(J):logpp(I,J)*cishu(I,J)));!失败率最低;@sum(yaohai(I): @sum(danyao(J):cishu(I,J)*(yh(J)*jl(I)+200)))<=48000;!油耗限制;@for(danyao(J): @sum(yaohai(I):cishu(I,J))<=ssl(J));!弹药限制;@for(fangan:@gin(cishu));@for(fangan:logpp=-@log(1-pp));!每次的失败率的自然对数;data: jl=450480540600; ssl=3248; pp=.1 .08 .2 .16 .15 .12 .25 .20; yh=.583333.75;enddataend7.有四个工人,要指派他们分别完成4项工作,每人做各项工作所消耗的时间如表。问指派哪个人去完成哪项工作,可使所有工作被完成时所用的时间(假设所有工作并行完成)(用lingo求解)?程序:model:!4个工人,4个工作的分配问题;sets:workers/w1..w4/;jo

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论