13.建模作业_优化问题.doc

《13.建模作业_优化问题.doc》由会员分享，可在线阅读，更多相关《13.建模作业_优化问题.doc（13页珍藏版）》请在三一文库上搜索。

1、数学建模课程作业题-13第五章 优化模型-优化问题1.已知某工厂计划生产I,II,III三种产品，各产品需要在A,B,C设备上加工，有关数据如下： 产品设备IIIIII设备有效台数（每月）A8210300B1058400C21310420单位产品利润 (千元)322.9试问:（1）如何发挥生产能力，使生产盈利最大？模型的建立及求解：设生产I,II,III产品x1，x2，x3件z为所获得的利润。于是数学模型如下：利用matlab求解（附录一）得到最优值Z =135.2667（千元），生产方案如下表。产品IIIIII数量23237生产I,II,III产品分别为23,23,7利润最大为125.266

2、7千元。（2）若为了增加产量，可租用别的工厂设备B,每月可租用60台，租金1.8万元，租用B设备是否划算？模型的建立及求解：租用别的工厂设备B以后模型为：利用matlab求解（附录二）得到最优值Z =129（千元），生产方案如下表。产品IIIIII数量31280生产I,II,III产品分别为31,28,0利润最大为129千元。（3）若另有俩种新产品、,其中新产品需用设备A为12台时，B为5台时，C为10台时，单位产品盈利2.1千元；新产品需用设备A为4台时，B为4台时，C为12台时，单位产品盈利1.87千元，如A,B,C的设备台时不增加，这两种新产品投产在经济上是否划算？ 模型的建立及求解：添

3、加两个新产品、后，、对应的产品数分别为x4，x5，建立模型如下：利用matlab求解（附录三）得到最优值Z =136.9625（千元），生产方案如下表。产品IIIIII数量27160014生产I,II,III，产品分别为27,16,0,0,14利润最大为136.9625千元。（4）对产品工艺重新进行设计，改进结构.改进后生产每件产品I需用设备A为9台时，设备B为12台时，设备C为4台时，单位盈利4.5千元，这时对原计划有何影响？ 模型的建立及求解：改进结构后，建立的模型如下：利用matlab求解（附录四）得到最优值Z =153.1618（千元），生产方案如下表。产品IIIIII数量23250生

4、产I,II,III产品分别为23,25,0利润最大为153.1618千元。 2. 有一个大型的冶金矿山公司，共有14个出矿点，已知其年产量及各矿点矿石的平均品位（含铁量的百分比）如下表所示：各矿点信息矿点号出矿量（万吨）平均铁品位（%）矿点号出矿量（万吨）平均铁品位（%）17037.16815.448.342751.2592.749.0831740.00107.640.2242347.001113.552.715342.00122.756.9269.549.96131.240.727151.41147.250.20按照炼铁生产要求，在矿石产出后，需按要求指定的品位值T进行不同品位矿石的混合配料

5、，然后进入烧结工序.最后，将小球状的烧结球团矿送入高炉进行高温炼铁，生产出生铁.该企业要求：将这14个出矿点的矿石进行混合配矿.依据生产设备及生产工艺要求，混合矿石的平均品位T规定为45% .问：应如何配矿才能获得最佳效益？模型的建立及求解： 设从第一矿点到第十四个矿点，每个矿点的配矿量分别为万吨（i表示矿点数），每个矿点铁的平均品味为。由题目给点条件，可得如下线性规划模型：（1）将（1）展开约束条件为混矿后的平均品味限制和各矿点的含矿量限制：（2）将（2）展开简化得： 得到最终模型： S.T利用matlab求解（附录五）得到最佳效益Max=63.8991，具体分配方案见下表。矿点123456

6、7891011121314产量317172339.5115.42.77.613.52.71.27.23. 三个家具商店购买办公桌：A需要30张，B需要50张，C需要45张.这些办公桌由两个工厂供应：工厂1生产70张，工厂2生产80张.下表给出了工厂和商店的距离（单位公里），假设每张每公里运费0.5元.寻求一个运送方案使运费最少？工厂和商店的距离工厂家具店ABC11053027205模型的建立及求解：设工厂1运给A x1a张，给B x1b，给C x1c张。工厂2运给A x2a张，给B x2b，给C x2c张，z表示最小费用。利用matlab求解（附录六）得到A B C分别在工厂1和工厂2的购买张

7、数，如下表：工厂家具店ABC10500230045最优方案为：工厂一运给A 店铺0张，给B 店铺50张，给C 店铺0张。工厂二运给A 店铺30张，给B 店铺0张，给C 店铺45张。总运费为342.5元 4. 某车间有一批长度为180公分的钢管（数量充分多），今为制造零件，要将其截成三种不同长度的管料，70公分，52公分，35公分.生产任务规定，这三种料的需要量分别不少于100根，150根，100根.所有截法如下表所示.我们知道，截钢管时不免要产生“边角料”，从节约原料的观点来考虑，应该采取怎样的截法，才能在完成任务的前提下，使总的边角料达到最小限度？所有可能的截法截法（1）（2）（3）（4）（

8、5）（6）（7）（8）需要量长度702111000010052021032101503510130235100边料（cm）56235246235模型的建立及求解：设表示第i种方法截的数量，Z表示剩余边料的总和，为了节约材料，Z越小越好，而且还得满足各个长度的数量要求。建立模型如下：利用matlab求解（附录七）得到剩余边料最小值为600cm，具体截取方案如下表。方案12345678数量（根）1560011015055. 某人有一笔50万元的资金可用于长期投资,可供选择的投资机会包括购买国库券、公司债券、投资房地产、购买股票或银行保值储蓄等.不同的投资方式的具体参数如下表所示.投资者希望投资组合

9、的平均年限不超过5年，平均的期望收益率不低于13%，风险系数不超过4，收益的增长潜力不低于10%.问在满足上述要求的前提下投资者该如何选择投资组合使平均年收益率最高？投资参数序号投资方式（）投资期限（年）年收益率%风险系数增长潜力%1国库券311102公司债券10153153房地产6258304股票2206205短期存款110156长期储蓄5122107现金存款0300模型的建立及求解：设国库券、公司债券、房地产、股票、短期存款、长期储蓄、现金存款分别存，Z表示平均年收益，由题意可建立模型如下：利用matlab求解（附录八）得到最优年收益为17%，投资方案如下表格。存款方案国库券公司债券房地产

10、股票短期存款长期储蓄现金存款存款额（万）26.5782.790720.63120000 6. 设有M=400万元资金，要求4年内使用完，若在一年内使用资金万元，则可获得效益万元，效益不能再使用，当年不用的资金可存入银行，年利率为=10%，试制定出这笔资金的使用方案，使4年的经济效益总和最大.模型的建立及求解：设前四年使用的资金分别为万元，总的经济效益为Z，第一年使用了万元，则可剩余400-万元，则第一年末的时候得到的效益为万元，第二年可使用的资金为，第二年末得到经济效益为，第三年可使用的资金为万元，第三年末经济效益为，第四年可使用的资金为万元，第四年末总效益为因此可以建立模型如下：利用matl

11、ab求解（附录九）得到四年可获得最大效益为Z=43.0858万元，投资方案如下表所示。年份第一年投资额第二年投资额第三年投资额第四年投资额投资额（万元）86.3027103.7161127.0870152.2390 7. 某个中型的百货商场要求售货人员每周工作5天，连续休息2天，工资200元/周，已知对售货人员的需求经过统计分析如下表所示，问如何安排可使配备销售人员的总费用最少？销售人员调查表星期一二三四五六日所需售货员人数18151216191412开始休息的人数模型的建立及求解：设星期一到星期天每天休息的人数分别为由于要求售货人员每周工作5天，连续休息2天，工资200元/周，则可建立如下模

12、型：利用matlab求解（附录十）得到最少费用为46400元，每天安排休息的人数如下表。星期一二三四五六日所需售货员人数18151216191412开始休息的人数0000135148. 某工厂生产甲、乙两种产品，已知有关数据见下表。工厂在做决策时，要考虑如下的问题：（1）根据市场信息，产品甲的需求有所上升，故产品甲的产量大于乙的2倍；（2）超过计划供应的原材料时，需高价采购，这就增加成本；（3）不要使设备超负荷运行；（4）应尽可能达到并超过计划利润指标48；试问问如何安排生产？给出数学模型和计算结果。生产数据甲乙拥有量原材料51060设备4436利润68模型的建立及求解：设生产甲产品为，生产乙

13、产品为，获得利润为Z。根据问题一，应有；根据问题二，可知；根据问题三可知；根据问题四，要尽可能实现利润最大化。由此可建立如下模型：利用matlab求解（附录十一）得到最大利润Z=400，生产甲乙产品的数量如下表。 生产类别甲产品乙产品数量40.000020.0000附录一：%zs13_1_1.mc=-3,-2,-2.9;a=8,2,10;10,5,8;2,13,10;b=300,400,420;vlb=0,0,0;vub=;x,z=linprog(c,a,b,vlb,vub);xi=round(x)Z=-z附录二：%zs13_1_2.mc=-3,-2,-2.9;a=8,2,10;10,5,8;

14、2,13,10;b=300,460,420;vlb=0,0,0;vub=;x,z=linprog(c,a,b,vlb,vub);xi=round(x)Z=-z-18附录三：%zs13_1_3.mc=-3,-2,-2.9,-2.1,-1.87;a=8,2,10,12,4;10,5,8,5,4;2,13,10,10,12;b=300,400,420;vlb=0,0,0,0,0;vub=;x,z=linprog(c,a,b,vlb,vub);xi=round(x)Z=-z附录四：%zs13_1_4.mc=-4.5,-2,-2.9;a=9,2,10;12,5,8;4,13,10;b=300,400,4

15、20;vlb=0,0,0;vub=;x,z=linprog(c,a,b,vlb,vub);xi=round(x)Z=-z附录五：%zs13_2.mc=-0.3716 -0.5125 -0.4 -0.47 -0.42 -0.4996 -0.5141 -0.4838 -0.4908 -0.4022 -0.5271 -0.5692 -0.4072 -0.5020;a=0.0784,-0.0625,0.05, -0.02,0.03,-0.0496,-0.0641,-0.0338,-0.0408,0.0478,-0.0771,-0.1192,0.0428,-0.052;b=0;vlb=0;0;0;0;0

16、;0;0;0;0;0;0;0;0;0;vub=70;7;17;23;3;9.5;1;15.4;2.7;7.6;13.5;2.7;1.2;7.2;x,z=linprog(c,a,b,vlb,vub);disp(取最优值对应的各个矿点的产值)xi=xdisp(最优值)Z=-z附录六：%zs13_3.mc=10 5 30 7 10 5;a=1 1 1 0 0 0;0 0 0 1 1 1;b=70 80;aeq=1 0 0 1 0 0; 0 1 0 0 1 0;0 0 1 0 0 1;beq=30 50 45;vlb=0 0 0 0 0 0;vub=30 50 45 30 50 45;x,z=linp

17、rog(c,a,b,aeq,beq,vlb,vub);xi=round(x)Z=0.5*z附录七：%zs13_4.mc=5,6,23,5,24,6,23,5;a=-2,-1,-1,-1,0,0,0,0;0,-2,-1,0,-3,-2,-1,0;-1,0,-1,-3,0,-2,-3,-5;b=-100,-150,-100;vlb=0;0;0;0;0;0;0;0;vub=;x,z=linprog(c,a,b,vlb,vub);disp(取最优值所对应各个方案截的数量)xi=round(x)disp(最优值)z附录八：%zs13_5.mc=11,15,25,20,10,12,3./(-50);a=(

18、1/50).*3 10 6 2 1 5 0;1 3 8 6 1 2 0;0 -15 -30 -20 -5 -10 0;b=5 4 -10;aeq=1 1 1 1 1 1 1;beq=50;vlb=0;0;0;0;0;0;0;0;vub=;x,z=linprog(c,a,b,aeq,beq,vlb,vub);disp(投资方案)X=xdisp(最优年收益)Z=-z/100附录九：M文件%yueshu.mfunction g,ceq=yueshu(x) %定义非线性的约束条件 g(1)=x(1)-400; g(2)=1.1*x(1)+x(2)-440; g(3)=1.21*x(1)+1.1*x(2

19、)+x(3)-484; g(4)=1.331*x(1)+1.21*x(2)+1.1*x(3)+x(4)-532.4; ceq=0;endM文件%fun8.mfunction y=fun8(x)y=-(sqrt(x(1)+sqrt(x(2)+sqrt(x(3)+sqrt(x(4);end代码：%zs13_6.mx0=1;1;1;1;vlb=0;0;0;0;vub=;A=;b=;Aeq=;beq=; x,z=fmincon(fun8,x0,A,b,Aeq,beq,vlb,vub,yueshu);X=xZ=-z附录十：%zs13_7.mc=200*11 15 25 20 10 12 3;a=(-1)

20、*0 1 1 1 1 1 0;0 0 1 1 1 1 1;1 0 0 1 1 1 1;1 1 0 0 1 1 1;1 1 1 0 0 1 1;1 1 1 1 0 0 1;1 1 1 1 1 0 0;b=(-1)*18 15 12 16 19 14 12;aeq=;beq=;vlb=0;0;0;0;0;0;0;0;vub=;x,z=linprog(c,a,b,aeq,beq,vlb,vub);disp(周一到周七休息的人数分别为)X=xdisp(需要最少的费用)z附录十一：%zs13_8.mc=(-1)*6 8;a=-1 2;1 1;-4 -4;b=0 60 -36;aeq=;beq=;vlb=0;0;vub=;x,z=linprog(c,a,b,aeq,beq,vlb,vub);disp(甲乙产品分别生产)X=xdisp(需要最少的费用)Z=-z