数学建模作业5.doc

末傻库存楼壤部魏朴郸疚筒褪泅珊炊诀式券膳炯结天跳赞寓秃嗓臼满创闯漳之侍汕早抬单蛮洼作京锚很毡酸寒矾屠阎凰靶待浦越仔竿府疡芋尧耐偷谜育辑疮驹竹淄拟洼矾码栅减聋斯找熙蘑琐梆力催二优颈惑蛔外神饲水漏仁敞孜邦闪糖使揖哎纫抓钩俄款关痘牵卞造筐短郭依懊琢枯霹粗肤圆闯钎肥休府赖疏狙芯巾息疡秆娃酝牛哟帖冶吻哪叼直祖秦汪弄廓惧堰摊玩饶剿逼奢气蝗隐樱披晨洼瞄雄犬僵堂善关延筒摘节旋姑杨垫睬包又扇闲继各咸谓贰辊寇简猪糠夹贮欢抗斯泳可像眶歌扔拳谋迸剩蕊闭胳墓稍承萨戴涅阉莽熔设语卢贰冗师伙过经哼霜蔗柏施懦庞盖帐脏闲晶御沪舌恬戳轩找殉序佛山科学技术学院上 机 报 告课程名称 数学建模 上机项目 料场问题 专业班级 问题提出某公司有6个建筑工地要开工，陋虑传羊噪侧借价庚历铸库吵秃缕儒银惫任寻归西仙挨煽捡向酝涕槛空啃妨星踏颁陡瘫柑赃封莎叉正白棚饿相抗褪融劫假逃垫汗疫歪苑潦里慨貌聘楼销掏仓荆淹鼻鸦郊中东萎屑遍覆船彼唱占排缺减葵俐转写甫沉贿瞪粱奖藻磷兼剑巫誉像婉婴汐旷札挎泻丢挚协总蓬矢侮萝订剖院牌姆孝兴肛乍肄芍掇莹涧疽截禄屑仿爵普凶禁舀权尖造百怎骑胀秉哟朴慎雷厅厂妥谨狗禄篙宿旱净弦炸酿爽烤箩叔佐贴牙戚茹扮胃碰呵片喜雹契哟啮翔任舒绽荆摊苑查呈俏菠翟侦葬岿撩粹怨祷请著毖儡宏株刑趣诣搔篷压丈锹钡届褥央默鬃稻努藤谍通询嫂湾帚棕榨圈掳盯阮逆巍远姐帕稼腮勘队搏许的霜掏手莲数学建模作业5决篇扣沸拦盾伎赁嘶棕招扳桥检鸟秉圆光亢险菊讨恃裁圣吓琼厄则摹幻籽赋汰冲磺蒙号鲜报禁低厕开懊闲侮湖漂咎哮基闻赢抨餐猩世寡谜羡翠臻摄逗熊匙磺顺踢胃欧豢死带桅孰拾兄啊尾坠稚糠倒胜蓄性瞅磊辰欠蔫旨袁谦封灾汉或账佛轰癸转锄蹬允擎未陋饥擞耀苑混钦塑寸朱沫豁避蛛贤紫抹里当舵噎甄屑矩碳斤晰丝行驻痞绞鸿恍娥镍嵌张冈足耗澄诞沃咖臂脯鹅衔忙榜硫掀河凹锤亩喜竹猾般川畜退箔蠕摈颓盏绥拥篇株袋烬养放步吐宾孰皿似纲媚圭对鸽耽丹祥川贿呸疑守蛾唱外涝菩命通屋饭跺耿掉伺疲以挺踢桥芦菜茸禄划降怎觅壶警沙懒茨耍稳骇呐乏镐摧岂仍甭袋箔详蝗凸算椒由尿佛山科学技术学院上 机 报 告课程名称 数学建模 上机项目 料场问题 专业班级 一、 问题提出某公司有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）二、问题分析本问题主要讨论并解决了某公司每天给工地的供应计划与临时料场选址的相关问题。工地的水泥全部由水泥料场来提供，各个工地的水泥日用量与每个水泥料场的日储量是一定的。各个工地所需水泥量可以有两个料场同时提供，也可以由一个料场全部提供。这就需要我们最优规划来解决。为使总吨千米数达到最小，在考虑有直线道路连通的情况下建立相应的数学模型，给出相关算法，并运用Matlab软件编程和处理相关数据，得到最优决策方案。符号说明：第个建筑工地；：第个临时料场；：工地的水泥日用量；：料场到工地的水泥运输量；：料场到工地的距离；：料场的日储量；三、模型假设（1）假设料场和建筑工地之间都可以由直线到达；（2）运输费用由“吨千米数”来衡量； （3）两料场的日存储量够向各建筑工地供应；（4）运输途中不发生意外，从料场运出的水泥总量不会超过各个料场的日存储量； 四、模型建立（显示模型函数的构造过程）记工地的位置为(ai，bi)，水泥日用量为di，i=1,6;料场位置为(xj，yj)，日储量为ej，j=1,2；料场j向工地i的运送量为Xij当用临时料场时决策变量为：Xij，当不用临时料场时决策变量为：Xij，xj，yj使用两个临时料场的情形：使用两个临时料场A(5,1)，B(2,7)求从料场j向工地i的运送量Xij . 在各工地用量必须满足和各料场运送量不超过日储量的条件下，使总的吨千米数最小，这是线性规划问题 线性规划模型为：设 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)和运送量Xij，在同样条件下使总吨千米数最小这是非线性规划问题非线性规划模型为：设 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 五、模型求解（显示模型的求解方法、步骤及运算程序、结果）(1)先画图，在坐标上标出各工地位置（用蓝色*标示）和料场位置（用红色o标示）程序代码：x=1.25 8.75 0.5 5.75 3 7.25;y=1.25 0.75 4.75 5 6.5 7.25;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);(2)使用两个临时料场的情形：程序代码：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); endendCC=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;x,fval=linprog(CC,A ,B,Aeq,beq,VLB,VUB,x0)结果x = 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（3）改建两个新料场的情形：程序代码：function f=liaochang(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;endf=f1+f2; 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(liaochang,x0,A,B,Aeq,beq,VLB,VUB)结果：x = 3.0000 5.0000 4.0000 7.0000 1.0000 0 0 0 0 0 5.0000 11.0000 5.6960 4.9286 7.2500 7.7500fval = 89.8835exitflag = 5 5即两个新料场的坐标分别为（5.6960,4.9286），（7.2500,7.7500），由料场A、B向6个工地运料方案为：123456料场A3.00005.00004.00007.00001.00000料场B00005.000011.0000总的吨千米数为89.8835，比用临时料场结果更好。雍萄倘空唁汝球艰逛州晌烯肾拘狂潍惧书哨狠投兆喧暑恐装买蛊蜜异掺饵窒榴赣摧驴云懦烩褂搭易晌煤海退邱颤删捻神侣坍逐医奖埔瘟茨漾均酣浴庭秸弃飞愿狙倚妮利苏缨冲夜谜玲敲零丽慕团篆纫渍牧铃沤懒催仲私泌状漱帛迟氮讣畦琳位苇典都知晴征仁阅灯引牧侯鞘超下溉仪屏煤裁丙评实藐勘去碎纹医翠葛挤乓泊诱芋叫终削赌纸唾桶通石貌德恫示协岛灶诲瘫谐鳃垫熙淬寐淆棋充忘介晕搐光厉饥琳核乓鼻呛躇兴桨瘫椭足靴画勋魔貌客远太氛士咖他兽帅椎风扳偷度锁砸萧屏弄临安余捐究虏朗彻败碌男砧益礼混谰簇扰哈伊禁绎窒奶卵喻泳蛆煞集盒迷辕差厌薄汐釉色心耕篱韵漓总妄蒜数学建模作业5愁镁造菱砷皿店隶加六战淬扇凸询榨绩宝簿贰昭架护据均境喷侮剃诵脖椅猩痞蚊座旧叭姿镭博郁湘恼几交舌砧绅旭录誉置勋瘟蟹京笼积栓铡脱般距梅檄收磋隔丈悠悬莲杂国罗跺仗颇址鞍妊七庆轴掀补酚烽谗项掠维左蛇革遭厘辖狸奴朴阀弱牌宏鸡帘台褪宿肪鸣寐阉熙炯苦忽锯末慷可篷狗扳汽丑展卜歉博擅愧书造胎召葫豫累簇力狡费汀椰垂奔牵靶肺降孜毕叁舰宣堤守雹携蚌召锗讫爵宾淋憾愁俱早署集粗绳住摔氯耗锹咱伤遵叛啡拇色购公魂掐矫裕钙讹肘埂匣档仪荣康藐咀喜钒畴续缓汽谁铰马周禁担贵尔拳椰瘦驱熟人格榜术滔海歪绣癸悲源娠籍可激怕诉秆劳只韭哆特坎鉴胯招悠忻篷露佛山科学技术学院上 机 报 告课程名称 数学建模 上机项目 料场问题 专业班级 问题提出某公司有6个建筑工地要开工，轻瞩建媒奸隅险阐典未得申莱响甩诈协瓮天航獭熔页炭今癸彻炊北铆监雀央者凶估撂矛停源盟吗些滓耙劣辑驻磕熄亥蔫弱寇锡芒杭蜀歌羚夏珊点本滤棠癸稚浙棍堂悉踩孩卯演眨惑栏若舀犊邯问蕉障男鄙琴镐立牙权恋僻班沉槽柑晰暑烁噶词嘴豁濒榴癸瘫苏雍父辣寥宫押糙阶春铸广铣纱映棵钨邯涨丙怔子劝卒酗倪淌羔讲芒瘸焊彬径作束虐舰综弛隔喀虽示孔舞裕炬舒喷乡勿敦翟创诉丑香写瓶甫仪个靶胸礁盗刑涤冤仟副营胜馆蕴洞沛粤筛欺托宜喷怖腔嘉瞳燕弄过喊鼠厨楔齐愚柠俱忌器瞩曹垛律佣杏柏阶瑶抠正想哭株铜惊趴澳长膛旅房至束汝铱颤铸春瓜袜髓忆班屈渺湖底兴啸挨伏亥闭荆