数学建模实习

1、数学建模实习、数学规划模型1 问题得提出某企业将铝加工成A,B 两种铝型材，每5 吨铝原料就能在甲设备上用12小时加工成3吨A型材，每吨A获利2400元，或者在乙设备上用8小时加工成4吨 B 型材，每吨B 获利 1600 元。现在加工厂每天最多能得到250 吨铝原料，每天工人得总工作时间不能超过为480 小时，并且甲种设备每天至多能加工100吨A,乙设备得加工能力没有限制。（ 1）请为该企业制定一个生产计划，使每天获利最大。（ 2） 若用 1000元可买到1 吨铝原料，就是否应该做这项投资？若投资，每天最多购买多少吨铝原料？（ 3）如果可以聘用临时工人以增加劳动时间，付给工人得工资最多就是每小

2、时几元？（ 4）如果每吨A 型材得获利增加到3000 元，应否改变生产计划？2 问题分析与假设2、 1 问题分析我们为该企业制定得生产计划要使得每天获利最大，也就就是也就就是要确定分别用多少吨铝原料分配给甲、乙设备使得总利润最大，因此分配给甲、乙设备铝材料得吨数就就是我们这次线性规划得决策变量，由此就确定了获利得目标函数。同时目标函数又要满足一些约束条件，如每天最多能得到250 吨铝原料、每天工人得总工作时间不能超过为480小时、甲种设备每天至多能加工100吨 A，由此可以建立求解利润最大化得数学模型。对于问题2,3 上，仅仅改变相关参数，就可以得到最优解所得到最大值得变化，也就就是所谓得影子

3、价格，通过与影子价格比较，制定相应得策略。对于问题 4 可以从两个角度进行分析，一种直接改变参数，观察最优解就是否变化，另外一种就是对其进行灵敏度分析，观察其系数就是否落在取值范围内。2、 2 问题假设1假设加工A,B 型材得铝原料数就是满足铝原料供应得非负实数；2、假设就是在完全市场经济得情形下，进行问题分析得3假设A， B 型材每吨得获利就是与产量，所用时间就是相互独立得，即两两之间就是没有关系得。这三条假设就是进行线性规划，影子价格分析得基础3模型建立在模型建立之前，我们先给出如下记号Xi :分配给甲设备Xi个5吨铝材料X2：分配给乙设备X2个5吨铝材料W:每天得生产获利现在我们建立数学

4、模型各给甲、乙分配5吨原材料得情况下，原材料得生产 能力、消耗时间、获利之间得关系如下表：甲设备消耗5吨铝材料/ x1甲设备消耗5吨铝材料/ x2所用时间128厂量3吨从产品4吨B产品获利3 2400=72004 1600=6400由于5xi吨产品给甲设备生产，能够得到 3xi吨A产品，能够获利7200Xi元,并且5x2吨产品给乙设备生产，能够得到 4X2吨B产品，能够获利6400X2元，则建立目标函数 maXW 7200X1 6400x2并且由于1）加工厂每天最多能得到 250吨铝原料，2）每天工人得总工作 时间不能超过为480小时，并且3）甲种设备每天至多能加工100吨A,乙设备 得加工能

5、力没有限制，同时应满足非负约束，则其应满足约束条件：5 5xi 5X2 250J 12x1 8x2 4803xi 100L x1 0, X2 04模型求解4、1问题1求解根据建立得模型，我们可以知道maxW 7200x1 6400x2 ,我们利用Matlab 来进行求解，将问题极小化以便处理，即M W ,则当ox2时maxW得最优解时，也就就是min M得最有解，进而将问题转化为:min M7200X, 6400x2L 5X15x225012x18x24803x1100-X10, X20解得 x1 20, x2 30minM 3.36?105那么 maxW 3.36?105此时最优解就是Xi

6、20, X2 30,即分配给甲100吨，分配给乙150吨，此时获得最大利润3.36?105万元。4、2问题得进一步求解4. 2、1若用1000元可买到1吨铝原料，即加工厂每天最多能得到 251吨 铝原料，则可以原问题转化为min M7200X1 6400X25x1 5x2 25112X1 8x2 4803X1 100j X1 0, x2 0此时最优解X1 19, X2 36并且 M=-3.3696 ?105 ,则 maxW=3.3696 ?105,利润相比增力口了 960 1000。那么此时说明不应该做这项投资。其实这960元就就是铝材料得影子价格，在完全市 场经济得条件下，由于该资源得价格高

7、于影子价格，则此时企业应当卖掉该资源， 而不就是扩大生产。4、2、2同样对于问题3,当工时增加1小时，企业能够获得多获得大得利 益，同样这也就是工人工资得影子价格，与上述方法类似计算得到工人得影子价 格就是200元，即付给工人得工资最多不能超过 200元，否则还不如不要这个工 人所带来得收益更大。4、2、3对于问题4,当A产品得获利变为每吨为3000元时，改变相关参数，与上述3个问题类似，我么可以得到最优解仍就是 X, 20,X2 30,并且我们求得maxW=3.72 ?105元。所以不需要改变生产计划。4、2、4接下来进行参数得灵敏度分析，使用 LINGO11软件进行灵敏度分析数学建模实习较

8、好，所以在这使用了这个软件进行分析，软件运行截图如下：Ranges in which the basis is unchanged:VariableX1 0000X20000CurrentAllowableAllowableCoefficientIncreaseDecrease7200、 0002400、 0000800、6400、 000800、 00001600、Objective Coefficient Ranges2 问题得分析与假设Righthand Side RangesRowCurrentAllowableAllowableRHSIncreaseDecreaseSTUFFTIME

9、CAPACITY250、 0000480、 0000100、 000050、 0000053、 33333INFINITY33、 3333380、 0000040、 00000得 到x1 得 取 值 范 围 7200 800,7200 2400 ， 得 到x2 得 取 值 范 围6400 1600,6400 800 ，在这样得取值范围内，最优解不改变。显然问题4 求解得 9000 就落在取值范围，从这角度也验证了此时最优解不改变。铝原料最多可增加50吨，劳动时间最多可增加53、3h,在这些取值范围内，进行影子价格得讨论才就是有意义得。5 模型评价该模型就是在在完全市场经济得情形下，假设A， B

10、 型材每吨得获利就是与产量， 所用时间相互独立得情况下，建立线性规划模型，前者就是影子价格分析得前提，后者就是线性得模型。对于问题1,2,3 ，笔者都就是利用matlab 分析得， matlab 进行求最优解与影子价格得求解还就是很方便得，但就是不适合对于参数得敏感性得分析，这里笔者采用了lingo 软件， 得到内容丰富得输出，虽然输出得内容很多，这里笔者仅仅使用了一部分笔者关心得数据，事实上lingo就是解决线性规划问题很好得软件，这对于未来解决实际问题也就是十分有用得。该模型就是个简单得线性规划问题，运用了单纯形法、对偶问题、影子价格等问题，事实上该模型不仅适应与生产计划得制定，同时也可以

11、应用于市场销售、库存管理、运输问题等领域，此外还有合理下料、配料问题、物料管理等方面，分析方法、建模方法都就是类似得。由于该模型仅仅适应线性规划，实际上， 我们还可以引进正负偏差量d ,d等， 将线性规划变为目标规划，这就能够适应更多得情形，而这些引进都就是当市场条件改变得时候才发生改变。7 程序附录一8 参考文献1 运筹学甘应爱、田丰等清华大学出版社背景 20052 数学建模案例分析白其峥海洋出版社2000 年 北京3 数学实验初步肖海军 科学出版社2007 北京二、鱼得捕捞问题1 问题提出作为鱼塘得管理者来说，在相同得初始鱼苗量得情况下，养殖费基本就是固定得，但就是经济效益会随着捕捞策略得

12、改变而改变，那么采取何种捕捞策略，就会使得在固定得投入下，产生更大得效益呢？捕捞问题就就是在这样得背景下产生得。1、1在鱼塘中投放n。尾鱼苗，随着时间得增长，尾数将减少而每尾得重量将增加。设尾数n(t)得(相对)减少率为常数；由于喂养引起得每尾鱼重量得增加率与鱼得表面积成正比，由于消耗引起得每尾鱼重量得减少率与质量本身成正比。分别建立尾数与每尾鱼重得微分方程，并求解。1、2用控制网眼得办法不捕小鱼，到时刻T才开始捕捞，捕捞能力用尾数得.相对减少量|n/n|表示，记作E,即单位时间捕获量就是En(t)。问如何选择T与E,使从T开始得捕获量最大。数学建模实习2、 1 分析鱼得捕捞策略会受到很多因素

13、得影响，如气候、水温、天敌、中间斗争等因素都在影响捕捞策略得制定，但就是这些因素都不就是主要得，主要因素还就是捕捞能力及捕捞时机这两个。我们在建模得过程中，没有必要考虑所有得因素，只要抓住主要得、关键得因素做出合理得假设，我们这个模型就就是在抓住主要因素，忽略次要因素得基础上建立起来得。对于问题1、 1， 在不考虑捕捞得情况下，鱼得尾数n(t) 得 (相对 )减少率为常数，初始no尾鱼苗，那么可以依据此条件建立起微分方程，通过求解可以得到鱼得尾数n关于t得函数表达式。由于喂养引起得每尾鱼重量得增加率与鱼得表面积成正比，并且消耗引起得每尾鱼重量得减少率与质量本身成正比，这样我们就可以建立鱼重量得

14、变化率与质量以及表面积得关系，同时我们又知道质量与表面积也存在一定得关系，这就是由于质量与体积成正比(我们在这里假设鱼得密度 就是一个常数)，并且我们在假设鱼就是椭球体得情况下，鱼得体积与表面积存在某种固定关系，通过这些关系，我们可以建立鱼质量得变化率与鱼质量间得函数关系，通过取极限，我们就可以得到关于鱼重得微分方程。通过求解我们就可以得到鱼得重量m与t得关系。.对于问题1、 2，在有捕捞也有自然死亡得情况下，尾数得相对减少量| n/n |不仅与鱼得自然死亡有关，还与捕捞能力有关，我们先假设捕捞能力E 为一常数，从而可以建立起微分方程，求得通过求解可以得到鱼得尾数n关于t得函数表达式。同时对于

15、捕捞量W ,要使得从T开始得捕获量最大，那么就要 W t En(t)?mdt达到最大值，那么如何求得 T,将就是后面建模将要解决得 问题。2、 2 假设在建立模型之前，我们要进行合理得假设：2、2、1假设鱼彳马尾数n(t)就是关于t得连续可微函数2、2、2假设鱼彳#重量m就是关于t得连续可微函数2、 2、 3 假设鱼得密度就是一个常数2、 2、 4 假设鱼就是椭球体得，三个方向上得半径分别就是x, y,z2、2、5在不考虑捕捞得情况下，鱼得尾数 n(t)得(相对)减少率为常数a2、2、5喂养引起得每尾鱼重量得增加率与鱼得表面积成正比 ，比例为b2、2、6消耗引起得每尾鱼重量得减少率与质量本身成

16、正比，比例为 c2、2、7鱼得尾数得减少仅自然死亡与捕捞有关， 排除其她因素，如天敌、气候、环境污染、种内斗争等因素3模型得建立3、1在模型建立之前，我们先在这里给出相应得记号：n(t)鱼得尾数m(t)鱼得质量 a鱼得尾数n(t)得(相对)减少率为常数 mo鱼得初始质量b喂养引起得每尾鱼重量得增加率与鱼得表面积成正比得比例系数 c消耗引起得每尾鱼重量得减少率与质量本身成正比得比例系数x,y,z鱼就是椭球体得,三个方向上得半径分别就是 x,y,z鱼得密度就是一个常数V鱼得体积S鱼得表面积W总得捕捞量E捕捞能力T捕捞时机3、2无捕捞情况下，建立尾数n(t)得微分方程模型在鱼塘中投放no尾鱼苗,鱼得

17、尾数n(t)得(相对)减少率为常数为a则并且n(0)no很容易解得该方程得解为n(t)n°eatn(t V n/n(t) aVt两边令Vt 0有dnan(t) n(0)n0dt该模型符合鱼尾数就是逐渐减少情形3、3建立鱼重m(t)与t得时间关系首先我们知道由于喂养引起得每尾鱼重量得增加率与鱼得表面积成正比，比例系数为b。由于消耗引起得每尾鱼重量得减少率与质量本身成正比，比例系数 为c ,则m(t Vt) m(t)/八bS cm(t)Vt两边同时令Vt 0 ,得到dm bS cm(t) ,m(0) m0dt现在我们来建立S m之间的关系由于m V74V xyz 3由于目前对于椭球表面积

18、并没有确定得公式，在x,y,z相差不大得情况下，这里仅给出一个近似公式S 4 (xyz)2/3则 S3/26c 61mS (61m« (6)3m则 dm b(6匚)3m3 cm(t),我们再令 bi b();3dtdm 2得到 Dm3 cm(t), m(0) m0 dt3、4建立捕捞量W与捕捞时机T之间得关系首先由于有捕捞得情况下，鱼得尾数之间得微分方程关系就需要考虑捕捞量，并且由于单位时间捕获量就是 En(t),则建立微分方程有：dn (a E)n(t) (t T) dt并且由于T时刻后才有捕捞，则n(T) no?eaT很容易求得该微分方程得解为n(t) noeET?e(aE)t(

19、t T)又对于捕捞量W ,要使得从T开始得捕获量最大，那么就要W En(t)?m(t)dt T达到最大值则通过代入有捕捞量时得n(t),并代入质量m(t),就可以得到W得具体表达式。4模型得求解4、1求解无捕捞情况下，建立尾数n(t)得微分方程在无捕捞得情形下，之前已经求得n(t) n°e at ,现在我们要对这个参数a进行 合理得讨论，通过Matlab作图进行分析在不同得参数下尾数关于t得大致走势，这里应用no 5000 a 0.05,0.5,1,5这4个进行分析，做出图形。(程序附录2)信 & jj* jpa g： gMh jjjjiM jje 日虹密工gi gea 总,

20、触她卦i 加: gjul jtfiM * "寓 £© 如 ln& g RkIi 置g 也显然根据实际情况，a就是一个接近与0较小得数，否则养殖者将无法盈利， 当a大于1得情况在现实生活中也有可能出现，比如出现大规模得病害，或者天 敌得入侵，这都就是我们所不需要讨论得极端情况， 我们只就是针对一般情况作 出讨论，所以a就是一个较小得数。4、2求解鱼重m(t)与t得微分方程由之前建立好得微分方程，有dm .2.一Am3 cm(t)m(0) m0dt该方程利用变量分离得办法，很容易求出其解，笔者在这里利用Matlab来进行求解。(程序见附录2)程序求解为(b +

21、 1/exp(c*(C3 + t)/3)A3/cA3 ,代入初值化简可以得到 m(t)如下：m(t)但(师 b1)e即,由于m随着t得增大而增大，e晟关于t就是递减函 cc数，则要求m(t)就是递增函数，则要求床 b1) 0,这也就是说明满足该式 c就是鱼得固有属性。4、3求解捕捞量W与捕捞时机T由于捕捞量W T En(t)?m(t)dt其中t T ,有捕捞得情况下，n(t)noeET ?e (a E)t并且m(t)bict 3-)e 3 dt c则 WEnoeET ?e (a E)t?b1 (3 m0TcWEnoeET ? e (a E)t?b1 (/ b1)e 3t3dtccEnoeET(

22、 o "aE)t?但(晒 b1)eM3dt 00ccccET (a E)tbl/3 bl .3t,3ET T (a E)tbi, 3bi3t 13En°ee ?一 (3m0)e 3 dtEn°e? e?(3m0-)e3 dt0cc0c ' c由于就是实际问题e(aE)t?b1 (3/m； b1)e3dt存在，并且与T无关，仅与E0c ' c有关我们不妨设其值为M(E),为一常数贝 UW En0eET?Mbst 3-)e 3 3dt cT dt0 g c?e我们利用这个积分公式Inec3G.ec3T eln(G Q)可以将Te(a E)t?b1 (而

23、0 b1)e13出转变为一个关于T, E得函数，里面 0c ' c将不含积分，由于没有具体数据，这里给出具体得函数也就是没有多大意义得， 但就是很容易通过计算就是知道这个函数得具体得表达式得，我们不妨设其为f(T,E)ET -则 W En0e ?M(E) f(T,E)这就就是W关于E, T得二元函数，通过对二元函数求极值得方法，我们就 可以确定最佳得E,T,使得W达到最大，我们刚开始就是假设 E为一常数，那就 是为了更方便得建立方程，后来我们可以瞧到实际上E就是一个独立得参数，建 立起了 W关于E, T得函数，由于问题得规模已经超出了笔者得能力，这里仅 仅就是提出一个思路，具体得求解有

24、待未来知识储备加深后进行进一步得突破。5模型评价该模型将鱼得体积简化成为一个椭球体，以此为基础建立起了m(t)与t之间得关系，这就是该模型得出发点，并且假设鱼得质量就是分布均匀得， 并想从网上 找到鱼得密度等相关参数，发现没有前人对其具体测量过，这需要我们在对具体 养殖得时候要测出该鱼种类得密度等参数，将余简化为椭球体方便我们建立起微 分方程，事实上这样得假设也就是十分合理得。但就是我们这个模型并没有给出合适得捕捞时机 T与捕捞能力E , 一方面就是 由于条件得限制，相关参数，具体数据缺乏，另一方面也就是笔者能力有限，虽 知道无穷积分得存在性，但就是却无法通过已学知识给出具体得表达式， 这也就 是本次建模不足得地方，未来将在这个地方予以突破。数学建模实习6 模型得推广与改进该模型虽然就是建立对鱼得捕捞得基础上得，但