库存问题数学模1.doc

库存问题数学模型摘要本文主要针对某商店在鱼竿经营过程中，各方面因素对利润和成本的影响进行了综合分析。在鱼杆销售过程中，商店的利润由多方面原因组成：市场需求量的变化，订货费用，进货成本，库存费用以及厂家给的优惠条件等等。通过建立合理的模型，对库存问题建立合理的订货方案。 对于问题一，给定一组一年中各个月鱼杆的需求量值，由于在哪个月份订货，一年的订货次数，每次的订货量、库存都是不确定的，而且不同月份一个批量的订货费不同，每支鱼杆每月还需一定的贮存费，所以需要设出所有变量，通过建立数学函数表达式得到数学模型，最后在中实现，找到合理的订货方案。 对于问题二，在问题一的基础上，对其进行优化，由题意可得，增加了一个约束条件即如果鱼杆的订货数量超过250支，厂家将给予优惠，每支鱼杆的购置费降至120元，这是就需要通过设置0、1变量，根据订货数量的多少来确定每支鱼杆的购置费用，然后利用对目标函数进行优化，求出订货方案。最后将此订货方案与问题一中的订货方案进行比较，若在此约束条件下，成本降低了，则说明可以采取此订货方案，反之，则不采用此方案。 关键词： 成本最小 软件 0、1变量题目重述对鱼杆的需求量在12月份最小，4月份达到最大。一家鱼具商店预测，12月份的需求量为50支，随后每个月增加10支，到4月份达到90支。除了从2月到4月的高峰需求月份以外，普通月份一个批量的订货费为2500元，而高峰月份的订货费增加到3000元。每支鱼杆的购置费用大约是150元，全年不变，而每支鱼杆每月的存贮费用为10元。该鱼具商店正在制定下一年度（1月到12月）鱼杆订货计划。1. 商店的经理认为，鱼具属季节性商品，因此不允许缺货。试为该鱼具商店制订一份下一年度的订货方案。2. 其他条件不变。如果鱼杆的订货数量超过250支，厂家将给予优惠，每支鱼杆的购置费降至120元。那么，是否利用此项优惠。如果利用，全年的订货方案将如何改变？一、问题分析 问题1： 这个库存问题的目标是使成本最低,成本包括：订货费、每个月库存鱼杆的存贮费等，同时，我们还需要考虑的是订货次数、每个月鱼杆的需求量；约束条件为：商店对鱼杆数量的限制即总存货量不小于总需求量，这样才能保证全年不缺货；非高峰期的鱼杆需求量小于高峰期的鱼杆需求量； 第一个月的存货不小于第一个月的需求量。于是建立规划模型，可以利用LINGO进行求解。问题2：根据问题一的思想方法，增加一个约束条件即如果鱼杆的订货数量超过250支，厂家将给予优惠，每支鱼杆的购置费降至120元。再次利用进行求解。二、模型假设1、把各月的鱼杆的需求量假设出来；2、为了保证成本最小，第十二月不能有库存；3、每个月的订购量得最大值不受限制；4、本月订购量 + 上月库存量 = 本月需求量 + 本月库存量；三、符号说明符号说明每个月鱼杆的需求量，每个月鱼杆的订货量，每个月鱼杆的存货量，0、1变量，当b等于1时，这个月订购；当b等于0时，这个月不订购，0、1变量，当等于1时，鱼杆的订货数量超过250支，享受鱼杆购置费降至120元的优惠；当等于0时，鱼杆的订货数量小于250支，不享受购置费降价的优惠 四、模型建立 问题1的模型：Min= （1） 0 1 变量约束：为1或0(表示第 月份的订货，表示第 月份的不订货)； （2）第个月的存货量等于第个月的订货量加上第个月的存货量减去第个月的需求量, ； （3）第一个月的存货量等于第一个月的订货量同第一个月的需求量的差 ； （4）存货量不为0，保证全年不缺货 ； （5）为了保证花费最小，要是这12个月的订货量与需求量相等 ； （6）问题2 的模型：Min=（7）约束条件同问题1中的（2）（3）（4）（5）（6）；另外增加0 1 变量约束条件：为1或0(表示第 月份的订货数量大于250，表示第 月份的订货数量小于250)； （8）五、模型求解 将表达式输入到进行求解，模型在实际当中决策变量只能取整数，对得到的计算结果进行必要的取整，得到的问题1的鱼杆商店下一年订货方案如下所示月份123456789101112订货费（元）250030002500需求量607080906864626058565250订货量210002220 01800015800在此方案中最少花费133540元问题2的订货方案如下所示月份123456789101112订货费（元）250030002500需求量607080906864626058565250订货量251002510002680000在此方案中最少花费113330元六、检验评价本文考虑了鱼杆的需求量值、鱼杆的存贮费、订货次数、订货数量，第二个模型是依据第一个模型进行建立的。从理论上来说比较合理，但是考虑到在现实可能会有偏差，不过，此方法也有一定的意义，可以帮助商家赢得更大的收益。七、参考文献【1】姜启源,谢金星 ,叶俊.数学建模（第三版），高等教育出版社；【2】数学建模培训资料 济南大学理学院数学建模指导中心；【3】概率论与数理统计（第四版） 主编 盛骤 高等教育出版社。 八、附录问题一的程序如下：sets:month/1.12/:x,b,y,z; endsetsmin=SUM(month(i)|i#gt#4#and#i#lt#13:150*y(i)+2500*b(i)+150*y(1)+2500+ sum(month(i)|i#gt#1#and#i#lt#5:150*y(i)+3000*b(i)+sum(month:10*s);z(1)=y(1)-x(1); for(month(i)|i#gt#1:z(i)=y(i)-x(i)+z(i-1); y(1)-x(1)=0;for(month(i)|i#gt#1:y(i)-x(i)+z(i-1)=0);sum(month:x)=sum(month:y);for(month:y=b);for(month:y=0;for(month(i)|i#gt#1:y(i)-x(i)+z(i-1)=0);sum(month:x)=sum(month:y);for(month:y=b);for(month:y=c);for(month:y-2501000*c);data:a=60,70,80,90,68,64,62,60,58,56,52,50 ;enddataend问题二的运行结果如下：Global optimal solution found. Objective value: 113330.0 Extended solver steps: 1 Total solver iterations: 5312 Variable Value Reduced Cost Y( 1) 251.0000 190.0000 Y( 2) 0.000000 0.000000 Y( 3) 0.000000 0.000000 Y( 4) 251.0000 160.0000 Y( 5) 0.000000 0.000000 Y( 6) 0.000000 0.000000 Y( 7) 0.000000 0.000000 Y( 8) 268.0000 120.0000 Y( 9) 0.000000 0.000000 Y( 10) 0.000000 0.000000 Y( 11) 0.000000 0.000000 X( 12) 0.000000 0.000000 B( 1) 1.000000 0.000000 B( 2) 0.000000 0.000000 B( 3) 0.000000 0.000000 B( 4) 1.000000 3000.000 B( 5) 0.000000 0.000000 B( 6) 0.000000 0.000000 B( 7) 0.000000 0.000000 B( 8) 1.000000 2500.000 B( 9) 0.000000 0.000000 B( 10) 0.000000 0.000000 B( 11) 0.000000 0.000000 B( 12) 0.000000 0.000000 Z( 1) 191.0000 0.000000 Z( 2) 121.0000 0.000000 Z( 3) 41.00000 0.000000 Z( 4) 202.0000 0.000000 Z( 5) 134.0000 0.000000 Z( 6) 70.00000 0.000000 Z( 7) 8.000000 0.000000 Z( 8) 216.0000 0.000000 Z( 9) 158.0000 0.000000 Z( 10) 102.0000 0.000000 Z( 11) 50.00000 0.000000 Z( 12) 0.000000 0.000000 X( 1) 60.00000 0.000000 X( 2) 70.00000 0.000000 X( 3) 80.00000 0.000000 X( 4) 90.00000 0.000000 X( 5) 68.00000 0.000000 X( 6) 64.00000 0.000000 X( 7) 62.00000 0.000000 X( 8) 60.00000 0.000000 X( 9) 58.00000 0.000000 X( 10) 56.00000 0.000000 X( 11) 52.00000 0.000000 X( 12) 50.00000 0.000000 C( 1) 1.000000 -7530.000 C( 2) 0.000000 0.000000 C( 3) 0.000000 0.000000C( 4) 1.000000 -7530.000 C( 5) 0.000000 0.000000 C( 6) 0.000000 0.000000 C( 7) 0.000000 0.000000 C( 8) 1.000000 -8040.000 C( 9) 0.000000 0.000000 C( 10) 0.000000 0.000000 C( 11) 0.000000 0.000000 C( 12) 0.000000 0.000000 Row Slack or Surplus Dual Price 1 113330.0 -1.000000 2 0.000000 -70.00000 3 0.000000 -60.00000 4 0.000000 -50.00000 5 0.000000 -40.00000 6 0.000000 -30.00000 7 0.000000 -20.00000 8 0.000000 -10.00000 9 0.000000 0.000000 10 0.000000 10.00000 11 0.000000 20.00000 12 0.000000 30.00000 13 0.000000 -10.00000 14 191.0000 0.000000 15 121.0000 0.000000 16 41.00000 0.000000 17 202.0000 0.000000 18 134.0000 0.000000 19 70.00000 0.000000 20 8.000000 0.000000 21 216.0000 0.000000 22 158.0000 0.000000 23 102.0000 0.000000 24 50.00000 0.000000 25 0.000000 -50.00000 26 0.000000 0.000000 27 250.0000 0.000000 28 0.000000 -213.2132 29 0.000000 -203.2032 30 250.0000 0.000000 31 0.000000 -182.6827 32 0.000000 -172.6727 33 0.000000 -162.6627 34 267.0000 0.000000 35 0.000000 -142.6426 36 0.000000 -132.6326 37 0.000000 -122.6226 38 0.000000 -112.6126 39 749.0000 0.000000 40 0.000000 3.213213 41 0.000000 3.203203 42 74