《食油生产问题》word版.doc

食油生产问题目录1问题重述2问题分析3基本假设与符号说明4建模思想与过程5算法设计与计算结果参考文献附录问题重述食油厂通过精炼两种硬质原料油和软质原料油，得到一种食油，以下简称产品油。硬质原料油来自两个产地：产地1和产地2，而软质原料油来自另外三个产地：产地3，产地4和产地5。据预测，这5种原料油的价格从一至六月分别为表1所示，产品油售价200元/吨。硬质油和软质油需要由不同生产线来精炼。硬质油生产线最大处理能力为200吨/月，软质油生产线最大处理能力为250吨/月。五种原料油都备有贮罐，每个贮罐容量均为1000吨，每吨原料每月存贮费用为5元。而各种精制油以及产品无油罐可存贮。精炼的加工费用可略去不计，产品销售没有任何问题。产品油的硬度有一定的技术要求，它取决于各种原料油的硬度以及混合比例。产品油的硬度与各种成分的硬度以及所占比例成线性关系。根据技术要求，产品油的硬度必须不小于3.0而不大于6.0。硬质1、硬质2、软质3、软质4、软质5等各种原料油的硬度为8.8，6.1，2.0，4.2，5.0，其中硬度单位是无量纲的，并且这里假定精制过程不会影响硬度。假设在一月初，每种原料油有500吨存贮而要求在六月底仍保持同样贮备。（1） 根据表1预测的原料油价格，编制逐月各种原料油采购量、耗用量及库存量计划，使半年【注】内的利润最大。（2） 考虑原料油价格上涨对利润的影响。据市场预测分析，若二月份硬质原料油价格比表1中的数字上涨%，则软质油在二月份的价格将比表1中的数字上涨2%。相应地，三月份，硬质原料油将上涨2%，软质原料油将上涨4%，依此类推至六月份。试分析从1到20的各情况下，利润将如何变化？（3） 附加以下3个条件后，再求解上面的问题：1 每一个月所用的原料油不多于三种；2 若在某一个月用一种原料油，那么这种油不能少于20吨；3 若在一个月中用硬质油1或硬质油2，则这个月就必须用软质油5。注：书中所述为“本年”，由题意可知，应为“半年”。 表1 原料油的价格（元/吨）硬质1 硬质2 软质3 软质4 软质5一月 110 120 130 110 115二月 130 130 110 90 115三月 110 140 130 100 95四月 120 110 120 120 125五月 100 120 150 110 105六月 90 110 140 80 135问题分析问题一这是一个线性规划问题，要求根据表1预测的原料油价格，编制逐月各种原料油采购量、耗用量及库存量计划，使本年内的利润最大；可将利润作为目标函数，各种要求作为约束条件。由题意，利润=销售额-原料费-存贮费。约束条件有以下几类： （1）一月初和六月底的各种原料油存贮量相等。 （2）每月硬、软两类原料油不能超过生产线最大处理能力。 （3）各个月产品油的硬度不小于3.0且不大于6.0。 （4）六个月内任何时刻任何一种油的量不能超过1000吨。问题二在问题一的基础上，分析利润对原料油价格增量的灵敏度。问题三在问题一的基础上，增加以下约束条件即可（5）每月所用的原料油不多于三种。（6）若在某个月用一种原料油，则这种油不能少于20吨。（7）若在一个月中用硬质油1或硬质油2，则这个月就必须用软质油5。基本假设与符号说明基本假设1. 在每个月的月初购买各种原料油。2. 精炼过程中没有原油的损失，即精炼前后油的质量守恒。3. 产品油是五种原油的混合。4. 精炼的加工费用略去不计。5. 产品销售没有任何问题，即生产出来后可立即销售出去（这样，产品油就不会有存贮，符合题干中“各种精制油以及产品无油罐可存贮”）。6. 在每个月内（不考虑月初、月末这两个时间点），原料油的库存量随时间流逝均匀减少，即精炼所得的产品油的量随时间流逝均匀增加（均为一次函数）。7. 产品油的硬度与各种成分的硬度以及所占比例成线性关系。8. 精制过程不会影响硬度。9. 六月底每种原料油的存贮等于一月初的。符号说明表2符号说明M半年内的总利润P产品油的售价（200元/吨）H硬质油生产线最大处理能力（200吨/月）S软质油生产线最大处理能力（250吨/月）C原料油贮罐容量（1000吨）p每吨原料油每月贮存费用（5元/月）hj第j种原料油的硬度（8.8，6.1，2.0，4.2，5.0）xij第i个月月初第j种原料油的采购量yij第i个月第j种原料油的耗用量zij第i个月月末第j种原料油的库存量（zij=zi-1，j+xij-yij）A价格矩阵（65，数据见由表1，见注）B价格增量矩阵（65，数据由问题二得出，见注）zoij耗用量的0-1变量注：A=110,120,130,110,115, 130,130,110,90,115, 110,140,130,100,95, 120,110,120,120,125, 100,120,150,110,105, 90,100,140,80,135;B=0.00,0.00,0.00,0.00,0.00, 0.01,0.01,0.02,0.02,0.02, 0.02,0.02,0.04,0.04,0.04, 0.04,0.04,0.08,0.08,0.08, 0.16,0.16,0.32,0.32,0.32, 0.32,0.32,0.64,0.64,0.64;建模思想与过程问题一显然，本问题为最优化模型，易知，本题的目标函数为利润的函数，建立合理的假设后，利润与三部分有关，即销售额，原料费和存贮费，且，利润=销售额-原料费-存贮费，即： （1） 题目中要求比较多，约束条件也就比较多，有，（1） 六月底每种原料油的存贮等于一月初的。 （2）（2） 每个月硬质油的耗用量小于等于其每月生产线最大处理能力（200吨），每个月软质油的耗用量小于等于其每月生产线最大处理能力（250吨）。 （3） （4）（3） 各个月产品油的硬度不小于3.0且不大于6.0。 （5）（4） 六个月内，任何一种油的量不曾超过1000吨。 (6) (7)问题二 在问题一的基础上，分析利润对原料油价格增量的灵敏度。构建一个价格增量矩阵B，则Aij + (Aij Bij )即为第j种原料油在第i个月的价格。至此，可将问题一与问题二的模型合二为一，目标函数为:（8）约束条件同问题一中约束条件。对于问题一，只需将（8）式中的取为0即可。问题三 首先定义耗用量的0-1变量： 1 (yij0)zoij= 0 (yij=0) （9）在问题一的基础上，增加以下约束条件即可 （5）每月所用的原料油不多于三种。 （10）（6） 若在某个月用一种原料油，则这种油不能少于20吨。 （11）（7）若在一个月中用硬质油1或硬质油2，则这个月就必须用软质油5。 （12） （13）算法设计与计算结果计算结果问题一表3 采购量硬质1 硬质2 软质3 软质4 软质5一月 0 0 0 0 0二月 0 0 0 750 0三月 0 0 0 0 250四月 0 90.74074 0 0 0五月 409.2593 0 0 0 0六月 409.2593 290.7407 0 500 0表4 耗用量硬质1 硬质2 软质3 软质4 软质5一月 159.2593 40.74074 0 250 0二月 85.18519 114.8148 0 0 250三月 159.2593 40.74074 0 250 0四月 96.29630 103.7037 0 250 0五月 159.2593 40.74074 0 250 0六月 159.2593 40.74074 0 250 0表5 库存量硬质1 硬质2 软质3 软质4 软质5一月 340.7407459.25926 500 250 0二月255.55551 344.44446 500 1000 250三月 96.29621303.70372 500 750 500四月 -0.00009 290.74076 500 500 500五月249.99991250.00002 500 250 500六月499.99991499.99998 500 500 500 解得半年内的总利润为M=222231.5。问题二M1176913.72163920.03155733.34147866.75140500.06134386.77129673.38125160.09122346.710119600.011116860.012114270.013112480.014110790.015109900.016109010.017108120.018107230.019106340.020105450.0问题三所编程序调试未成功。（详情见附录）参考文献1徐玖平，胡知能，王緌.运筹学M.北京：科学出版社，2004.2肖华勇.实用数学建模与软件应用M.西安：西北工业大学，2008.附录问题一代码：model:sets:R/1.6/;C/1.5/; JZ(R,C):A,x,y;endsetsdata:A=110,120,130,110,115, 130,130,110,90,115, 110,140,130,100,95, 120,110,120,120,125, 100,120,150,110,105, 90,100,140,80,135;enddatamax=200*(sum(R(i):sum(C(j):x(i,j)-sum(JZ:A*x)-75000-5*(sum(C(j):12*x(1,j)-11*y(1,j)+10*x(2,j)-9*y(2,j)+8*x(3,j)-7*y(3,j)+6*x(4,j)-5*y(4,j)+4*x(5,j)-3*y(5,j)+2*x(6,j)-y(6,j)/2;for(C(j):sum(R(i):x(i,j)-y(i,j)=0);for(R(i):y(i,1)+y(i,2)=200);for(R(i):y(i,3)+y(i,4)+y(i,5)=0);for(R(i):(8.8*y(i,1)+6.1*y(i,2)+2.0*y(i,3)+4.2*y(i,4)+5.0*y(i,5)-6*(y(i,1)+y(i,2)+y(i,3)+y(i,4)+y(i,5)=0);for(C(j):x(1,j)=500);for(C(j):x(1,j)-y(1,j)+x(2,j)=500);for(C(j):x(1,j)-y(1,j)+x(2,j)-y(2,j)+x(3,j)=500);for(C(j):x(1,j)-y(1,j)+x(2,j)-y(2,j)+x(3,j)-y(3,j)+x(4,j)=500);for(C(j):x(1,j)-y(1,j)+x(2,j)-y(2,j)+x(3,j)-y(3,j)+x(4,j)-y(4,j)+x(5,j)=500);for(C(j):x(1,j)-y(1,j)+x(2,j)-y(2,j)+x(3,j)-y(3,j)+x(4,j)-y(4,j)+x(5,j)-y(5,j)+x(6,j)=0);for(C(j):500+x(1,j)-y(1,j)+x(2,j)-y(2,j)=0);for(C(j):500+x(1,j)-y(1,j)+x(2,j)-y(2,j)+x(3,j)-y(3,j)=0);for(C(j):500+x(1,j)-y(1,j)+x(2,j)-y(2,j)+x(3,j)-y(3,j)+x(4,j)-y(4,j)=0);for(C(j):500+x(1,j)-y(1,j)+x(2,j)-y(2,j)+x(3,j)-y(3,j)+x(4,j)-y(4,j)+x(5,j)-y(5,j)=0);for(C(j):500+x(1,j)-y(1,j)+x(2,j)-y(2,j)+x(3,j)-y(3,j)+x(4,j)-y(4,j)+x(5,j)-y(5,j)=0+x(6,j)-y(5,j);end问题二代码：model:sets:R/1.6/;C/1.5/; JZ(R,C):A,x,y,B;endsetsdata:A=110,120,130,110,115, 130,130,110,90,115, 110,140,130,100,95, 120,110,120,120,125, 100,120,150,110,105, 90,100,140,80,135;B=0.00,0.00,0.00,0.00,0.00, 0.01,0.01,0.02,0.02,0.02, 0.02,0.02,0.04,0.04,0.04, 0.04,0.04,0.08,0.08,0.08, 0.16,0.16,0.32,0.32,0.32, 0.32,0.32,0.64,0.64,0.64;enddatamax=200*(sum(R(i):sum(C(j):x(i,j)-sum(JZ:(A+A*B*1)*x)-75000-5*(sum(C(j):12*x(1,j)-11*y(1,j)+10*x(2,j)-9*y(2,j)+8*x(3,j)-7*y(3,j)+6*x(4,j)-5*y(4,j)+4*x(5,j)-3*y(5,j)+2*x(6,j)-y(6,j)/2;for(C(j):sum(R(i):x(i,j)-y(i,j)=0);for(R(i):y(i,1)+y(i,2)=200);for(R(i):y(i,3)+y(i,4)+y(i,5)=0);for(R(i):(8.8*y(i,1)+6.1*y(i,2)+2.0*y(i,3)+4.2*y(i,4)+5.0*y(i,5)-6*(y(i,1)+y(i,2)+y(i,3)+y(i,4)+y(i,5)=0);for(C(j):x(1,j)=500);for(C(j):x(1,j)-y(1,j)+x(2,j)=500);for(C(j):x(1,j)-y(1,j)+x(2,j)-y(2,j)+x(3,j)=500);for(C(j):x(1,j)-y(1,j)+x(2,j)-y(2,j)+x(3,j)-y(3,j)+x(4,j)=500);for(C(j):x(1,j)-y(1,j)+x(2,j)-y(2,j)+x(3,j)-y(3,j)+x(4,j)-y(4,j)+x(5,j)=500);for(C(j):x(1,j)-y(1,j)+x(2,j)-y(2,j)+x(3,j)-y(3,j)+x(4,j)-y(4,j)+x(5,j)-y(5,j)+x(6,j)=0);for(C(j):500+x(1,j)-y(1,j)+x(2,j)-y(2,j)=0);for(C(j):500+x(1,j)-y(1,j)+x(2,j)-y(2,j)+x(3,j)-y(3,j)=0);for(C(j):500+x(1,j)-y(1,j)+x(2,j)-y(2,j)+x(3,j)-y(3,j)+x(4,j)-y(4,j)=0);for(C(j):500+x(1,j)-y(1,j)+x(2,j)-y(2,j)+x(3,j)-y(3,j)+x(4,j)-y(4,j)+x(5,j)-y(5,j)=0);for(C(j):500+x(1,j)-y(1,j)+x(2,j)-y(2,j)+x(3,j)-y(3,j)+x(4,j)-y(4,j)+x(5,j)-y(5,j)=0+x(6,j)-y(5,j);end问题三代码（调试失败）：model:sets:R/1.6/;C/1.5/; JZ(R,C):A,x,y,B;endsetsdata:A=110,120,130,110,115, 130,130,110,90,115, 110,140,130,100,95, 120,110,120,120,125, 100,120,150,110,105, 90,100,140,80,135;B=0.00,0.00,0.00,0.00,0.00, 0.01,0.01,0.02,0.02,0.02, 0.02,0.02,0.04,0.04,0.04, 0.04,0.04,0.08,0.08,0.08, 0.16,0.16,0.32,0.32,0.32, 0.32,0.32,0.64,0.64,0.64;enddatamax=200*(sum(R(i):sum(C(j):x(i,j)-sum(JZ:(A+A*B*0)*x)-75000-5*(sum(C(j):12*x(1,j)-11*y(1,j)+10*x(2,j)-9*y(2,j)+8*x(3,j)-7*y(3,j)+6*x(4,j)-5*y(4,j)+4*x(5,j)-3*y(5,j)+2*x(6,j)-y(6,j)/2;for(C(j):sum(R(i):x(i,j)-y(i,j)=0);for(R(i):y(i,1)+y(i,2)=200;);for(R(i):y(i,3)+y(i,4)+y(i,5)=0);for(R(i):(8.8*y(i,1)+6.1*y(i,2)+2.0*y(i,3)+4.2*y(i,4)+5.0*y(i,5)-6*(y(i,1)+y(i,2)+y(i,3)+y(i,4)+y(i,5)=0);for(C(j):x(1,j)=500);for