数学建模实习

一、数学规划模型1问题的提出某企业将铝加工成A,B两种铝型材，每5吨铝原料就能在甲设备上用12小时加工成3吨A型材，每吨A获利2400元，或者在乙设备上用8小时加工成4吨B型材，每吨B获利1600元。现在加工厂每天最多能得到250吨铝原料，每天工人的总工作时间不能超过为480小时，并且甲种设备每天至多能加工100吨A，乙设备的加工能力没有限制。（1）请为该企业制定一个生产计划，使每天获利最大。（2）若用1000元可买到1吨铝原料，是否应该做这项投资？若投资，每天最多购买多少吨铝原料？（3）如果可以聘用临时工人以增加劳动时间，付给工人的工资最多是每小时几元？（4）如果每吨A型材的获利增加到3000元，应否改变生产计划？2问题分析与假设2.1问题分析我们为该企业制定的生产计划要使得每天获利最大，也就是也就是要确定分别用多少吨铝原料分配给甲、乙设备使得总利润最大，因此分配给甲、乙设备铝材料的吨数就是我们这次线性规划的决策变量，由此就确定了获利的目标函数。同时目标函数又要满足一些约束条件，如每天最多能得到250吨铝原料、每天工人的总工作时间不能超过为480小时、甲种设备每天至多能加工100吨A，由此可以建立求解利润最大化的数学模型。对于问题2,3上，仅仅改变相关参数，就可以的到最优解所得到最大值的变化，也就是所谓的影子价格，通过与影子价格比较，制定相应的策略。对于问题4可以从两个角度进行分析，一种直接改变参数，观察最优解是否变化，另外一种是对其进行灵敏度分析，观察其系数是否落在取值范围内。2.2问题假设1．假设加工A,B型材的铝原料数是满足铝原料供应的非负实数；2.假设是在完全市场经济的情形下，进行问题分析的3．假设A，B型材每吨的获利是与产量，所用时间是相互独立的，即两两之间是没有关系的。这三条假设是进行线性规划，影子价格分析的基础3模型建立在模型建立之前，我们先给出如下记号1x：分配给甲设备1x个5吨铝材料2x：分配给乙设备2x个5吨铝材料W：每天的生产获利现在我们建立数学模型各给甲、乙分配5吨原材料的情况下，原材料的生产能力、消耗时间、获利之间的关系如下表：甲设备消耗5吨铝材料/1x甲设备消耗5吨铝材料/2x所用时间128产量3吨A产品4吨B产品获利32400=720041600=6400由于51x吨产品给甲设备生产，能够得到31x吨A产品，能够获利72001x元，并且52x吨产品给乙设备生产，能够得到42x吨B产品，能够获利64002x元，则建立目标函数12max72006400Wxx并且由于1）加工厂每天最多能得到250吨铝原料，2）每天工人的总工作时间不能超过为480小时，并且3）甲种设备每天至多能加工100吨A，乙设备的加工能力没有限制，同时应满足非负约束，则其应满足约束条件：12121125525012848031000,0xxxxxxx4模型求解4.1问题1求解根据建立的模型，我们可以知道12max72006400Wxx，我们利用Matlab来进行求解，将问题极小化以便处理，即MW，则当12,xx时maxW的最优解时，也就是minM的最有解，进而将问题转化为：12min72006400Mxx12121125525012848031000,0xxxxxxx解得1220,30xx53.36mn10iM那么53.36mx10aW此时最优解是1220,30xx，即分配给甲100吨，分配给乙150吨，此时获得最大利润53.3610万元。4.2问题的进一步求解4．2.1若用1000元可买到1吨铝原料，即加工厂每天最多能得到251吨铝原料，则可以原问题转化为12min72006400Mxx12121125525112848031000,0xxxxxxx此时最优解1219,36xx并且5M=-3.369610，则5maxW=3.369610，利润相比增加了9601000。那么此时说明不应该做这项投资。其实这960元就是铝材料的影子价格，在完全市场经济的条件下，由于该资源的价格高于影子价格，则此时企业应当卖掉该资源，而不是扩大生产。4.2.2同样对于问题3，当工时增加1小时，企业能够获得多获得大的利益，同样这也是工人工资的影子价格，和上述方法类似计算得到工人的影子价格是200元，即付给工人的工资最多不能超过200元，否则还不如不要这个工人所带来的收益更大。4.2.3对于问题4，当A产品的获利变为每吨为3000元时，改变相关参数，与上述3个问题类似，我么可以得到最优解仍是1220,30xx，并且我们求得5maxW=3.7210元。所以不需要改变生产计划。4.2.4接下来进行参数的灵敏度分析，使用LINGO11软件进行灵敏度分析较好，所以在这使用了这个软件进行分析，软件运行截图如下：Rangesinwhichthebasisisunchanged:ObjectiveCoefficientRangesCurrentAllowableAllowableVariableCoefficientIncreaseDecreaseX17200.0002400.0000800.0000X26400.000800.00001600.0000RighthandSideRangesRowCurrentAllowableAllowableRHSIncreaseDecreaseSTUFF250.000050.0000033.33333TIME480.000053.3333380.00000CAPACITY100.0000INFINITY40.00000得到1x的取值范围[7200800,72002400]，得到2x的取值范围[64001600,6400800]，在这样的取值范围内，最优解不改变。显然问题4求解的9000就落在取值范围，从这角度也验证了此时最优解不改变。铝原料最多可增加50吨，劳动时间最多可增加53.3h，在这些取值范围内，进行影子价格的讨论才是有意义的。5模型评价该模型是在在完全市场经济的情形下，假设A，B型材每吨的获利是与产量，所用时间相互独立的情况下，建立线性规划模型，前者是影子价格分析的前提，后者是线性的模型。对于问题1,2,3，笔者都是利用matlab分析的，matlab进行求最优解和影子价格的求解还是很方便的，但是不适合对于参数的敏感性的分析，这里笔者采用了lingo软件，得到内容丰富的输出，虽然输出的内容很多，这里笔者仅仅使用了一部分笔者关心的数据，事实上lingo是解决线性规划问题很好的软件，这对于未来解决实际问题也是十分有用的。6模型推广与优化该模型是个简单的线性规划问题，运用了单纯形法、对偶问题、影子价格等问题，事实上该模型不仅适应与生产计划的制定，同时也可以应用于市场销售、库存管理、运输问题等领域，此外还有合理下料、配料问题、物料管理等方面，分析方法、建模方法都是类似的。由于该模型仅仅适应线性规划，实际上，我们还可以引进正负偏差量,dd等，将线性规划变为目标规划，这就能够适应更多的情形，而这些引进都是当市场条件改变的时候才发生改变。7程序附录一8参考文献[1]运筹学甘应爱、田丰等清华大学出版社背景2005[2]数学建模案例分析白其峥海洋出版社2000年北京[3]数学实验初步肖海军科学出版社2007北京二、鱼的捕捞问题1问题提出作为鱼塘的管理者来说，在相同的初始鱼苗量的情况下，养殖费基本是固定的，但是经济效益会随着捕捞策略的改变而改变，那么采取何种捕捞策略，就会使得在固定的投入下，产生更大的效益呢？捕捞问题就是在这样的背景下产生的。1.1在鱼塘中投放0n尾鱼苗，随着时间的增长，尾数将减少而每尾的重量将增加。设尾数()nt的(相对)减少率为常数;由于喂养引起的每尾鱼重量的增加率与鱼的表面积成正比，由于消耗引起的每尾鱼重量的减少率与质量本身成正比。分别建立尾数和每尾鱼重的微分方程，并求解。1.2用控制网眼的办法不捕小鱼，到时刻T才开始捕捞，捕捞能力用尾数的相对减少量表示，记作E，即单位时间捕获量是E()nt。问如何选择T和E，使从T开始的捕获量最大。2问题的分析与假设2.1分析鱼的捕捞策略会受到很多因素的影响，如气候、水温、天敌、中间斗争等因素都在影响捕捞策略的制定，但是这些因素都不是主要的，主要因素还是捕捞能力及捕捞时机这两个。我们在建模的过程中，没有必要考虑所有的因素，只要抓住主要的、关键的因素做出合理的假设，我们这个模型就是在抓住主要因素，忽略次要因素的基础上建立起来的。对于问题1.1，在不考虑捕捞的情况下，鱼的尾数()nt的(相对)减少率为常数，初始0n尾鱼苗，那么可以依据此条件建立起微分方程，通过求解可以得到鱼的尾数()nt关于t的函数表达式。由于喂养引起的每尾鱼重量的增加率与鱼的表面积成正比，并且消耗引起的每尾鱼重量的减少率与质量本身成正比，这样我们就可以建立鱼重量的变化率和质量以及表面积的关系，同时我们又知道质量与表面积也存在一定的关系，这是由于质量与体积成正比（我们在这里假设鱼的密度是一个常数），并且我们在假设鱼是椭球体的情况下，鱼的体积与表面积存在某种固定关系，通过这些关系，我们可以建立鱼质量的变化率和鱼质量间的函数关系，通过取极限，我们就可以得到关于鱼重的微分方程。通过求解我们就可以.|/|nn的到鱼的重量()mt与t的关系。对于问题1.2，在有捕捞也有自然死亡的情况下，尾数的相对减少量不仅与鱼的自然死亡有关，还与捕捞能力有关，我们先假设捕捞能力E为一常数，从而可以建立起微分方程，求得通过求解可以得到鱼的尾数()nt关于t的函数表达式。同时对于捕捞量W，要使得从T开始的捕获量最大，那么就要()()TWEntmtdt达到最大值，那么如何求得T，将是后面建模将要解决的问题。2.2假设在建立模型之前，我们要进行合理的假设：2.2.1假设鱼的尾数()nt是关于t的连续可微函数2.2.2假设鱼的重量()mt是关于t的连续可微函数2.2.3假设鱼的密度是一个常数2.2.4假设鱼是椭球体的，三个方向上的半径分别是,,xyz2.2.5在不考虑捕捞的情况下，鱼的尾数()nt的(相对)减少率为常数a2.2.5喂养引起的每尾鱼重量的增加率与鱼的表面积成正比,比例为b2.2.6消耗引起的每尾鱼重量的减少率与质量本身成正比，比例为c2.2.7鱼的尾数的减少仅自然死亡和捕捞有关，排除其他因素，如天敌、气候、环境污染、种内斗争等因素3模型的建立3.1在模型建立之前，我们先在这里给出相应的记号：()nt鱼的尾数()mt鱼的质量a鱼的尾数()nt的(相对)减少率为常数0m鱼的初始质量b喂养引起的每尾鱼重量的增加率与鱼的表面积成正比的比例系数c消耗引起的每尾鱼重量的减少率与质量本身成正比的比例系数,,xyz鱼是椭球体的，三个方向上的半径分别是,,xyz鱼的密度是一个常数V鱼的体积S鱼的表面积W总的捕捞量E捕捞能力T捕捞时机.|/|nn3.2无捕捞情况下，建立尾数()nt的微分方程模型在鱼塘中投放0n尾鱼苗，鱼的尾数()nt的(相对)减少率为常数为a则()()/()nttntntat并且0(0)nn两边令0t有()dnantdt0(0)nn很容易解得该方程的解为0()atntne该模型符合鱼尾数是逐渐减少情形3.3建立鱼重()mt与t的时间关系首先我们知道由于喂养引起的每尾鱼重量的增加率与鱼的表面积成正比，比例系数为b。由于消耗引起的每尾鱼重量的减少率与质量本身成正比，比例系数为c，则()()()mttmtbScmtt两边同时令0t，得到()dmbScmtdt,0(0)mm现在我们来建立Sm与之间的关系由于mV43Vxyz由于目前对于椭球表面积并没有确定的公式，在,,xyz相差不大的情况下，这里仅给出一个近似公式2/34()Sxyz则3/266mSV2223336(6)()mSm则22336()()dmbmcmtdt，我们再