数学建模-生产与存贮问题的探讨

生产与存贮问题的探讨摘要在一定时期内,生产的成本费与库存费一直是厂家最关心的优化指标。本文根据题中的条件针对如何在成本费与库存费之和最优的情况下，使总工时最小的问题，利用了多目标动态规划的方法，建立了生产与存储的优化模型。我们知道增大生产量可以降低成本费，但如果超过市场的需求量,就会因积压增加存贮费而造成损失。相反，如果减少生产量，虽然可以降低存贮费，但又会增加生产的成本费，同样会造成损失。故可以找到一个生产计划使得生产的生产费与存贮费之和达到一个最小值,并且使他们所花的工时也最少。我们根据实际生活中生产的部件的性质可以将生产模式分成两种情况：允许有缺货的情况和不允许有缺货的情况。在模型一中,我们假设这种部件是不允许缺货的，于是目标函数为：6161)(7.03.0minkkkkkkchpakxg在模型二中,我们假设这种部件是可以缺货的，但是我们要求上个月所缺的部件必须要在本月补回来。如果中间某个月或者是某几个月出现缺货的现象,就会因为有损失费,面对这样的情况时,如果损失费比生产费少的话,对于这种方案公司还是可以考虑,根据这种情况我们可以得到目标函数为：6161)(7.03.0minkkkkkkkqphcakxg我们建立的模型一和模型二都是以动态规划为主要解题思路，在模型中我们将生产费与库存费之和赋予0.7的权重值，总耗费工时数赋予0.3的权重值，假设每件产品的单位工时费为10元，每件产品每月的存贮费为20元，每件产品每月的缺货损失费为5元，因为产品的生产量与成本费成反比，设反比系数为S，若生产量为X，则成本费为S/X元，设反比系数S为840。我们利用Lingo软件求解，在没有缺货存在的条件下得到的最小成本费为5158元，总耗费工时数最少为382小时，一到六月的逐月分配方案为：745434；在有缺货存在的条件下得到的最小成本费为4960元，总耗费工时数最少为363小时，一到六月的逐月分配方案为：634338，每月的缺货量为：021040。最后我们对模型进行了分析，并得到模型的整体评价和推广前景。关键词：允许缺货不允许缺货动态规划灵敏度分析Lingo权重一、问题重述一个生产项目，在一定时期内，增大生产量可以降低成本费，但如果超过市场的需求量,就会因积压增加存贮费而造成损失。相反，如果减少生产量，虽然可以降低存贮费，但又会增加生产的成本费，同样会造成损失.因此，如何正确地制定生产计划，使得在一定时期内,生产的成本费与库存费之和最小，这是厂家最关心的优化指标，这就是生产与存贮问题。假设某车间每月底都要供应总装车间一定数量的部件。但由于生产条件的变化，该车间每月生产单位部件所耗费的工时不同，每月的生产量除供本月需要外，剩余部分可存入仓库备用。今已知半年内，各月份的需求量及生产该部件每单位数所需工时数如下表所示：月份k123456月需求量bk853274单位用时ak111813172010设库存容量H=9，开始时库存量为2，期终库存量为0。要求制定一个半年逐月生产计划，使得既满足需求和库存容量的限制，又使得总耗费工时数最少。二、问题分析本题是典型的运筹学当中的动态规划问题。题目中要求制定一个半年逐月生产计划，使得既满足需求和库存容量的限制，又使得总耗费工时数最少。考虑到生产的成本费与库存费之和最小是厂家最关心的优化指标，而题目中没有给出相关数据，所以我们增设了有关成本费和库存费的条件来求得最优解。此外，开始时库存量为2，期终库存量为0。我们根据实际生活中部件的性质，可以分为允许缺货和不允许缺货两种情况考虑。前一种情况是：供货期间某个月或者是某几个月可以出现缺货,但是在它的下个月会补齐；第二种是：每个月都不能缺货,到最后一个月,全部销售出,没有库存。我们针对这两种情况建立允许缺货存贮模型和不允许缺货存贮模型展开不同的讨论，进而结合题目要求求出最优解。三、符号系统kx：第k个月该产品的生产量kd：第k个月该产品的需求量kv：第k个月结束时的产品库存量，则有kkkkdxvv1kc：第k个月生产产品kx时的成本费用kh：第k个月结束时有库存量kv所需的库存费用m：每个月最多能生产该产品的上限数kq：每件产品每月的缺货损失费ak：每个月每件产品的单位工时kp：第k个月的工时费四、模型建立4.1不允许缺货的存贮模型4.1.1模型假设1、每件产品每月的存贮费为常数2、每件产品单位工时的成本费为常数3、每个月都不会出现缺货的情况,即生产能力无限大,但最大不超过15件4、每个月的生产量不为零4.1.2模型建立设kd为第k个月该产品的需求量，kx为第k个月该产品的生产量，kv为第k个月结束时的产品库存量，则有kkkkdxvv1。kc为第k个月生产产品kx时的成本费用，它包括生产准备成本费K和产品成本kax（其中a是单位产品成本）两项费用，即mxmxaxKxckkkkk当当当,,2,100kh表示第k个月结束时有库存量kv所需的库存费用，kp为第k个月的工时费故第k个月的总成本费用为kkkphc，因此，不考虑缺货条件下的数学模型为)6,,2,1(,),,2,1(,0)5,,2(,0)(0,2][7.03.0min1606161kZxnkmxkdxvvvchpxakkkkjjjkkkkkkk用动态规划方法来求解，把它看作一个6阶段决策问题。令1kv为状态量，它表示第k个月开始开始时的库存量。kx为决策变量，它表示第k个月的生产量。状态转移方程为)6,,2,1(,1kdxvvkkkk最优值函数)(kkvf表示从第1个月初始库存量为2到第k个月末库存量为kv时的最小总费用。因此可写出顺序递推关系式为)6,,1(,)(min)(110kvfphcvfkkkkkxkkkk其中),min(mdvkkk。这是因为一方面每个月生产的上限为m；另一方面由于保证供应，故第k-1个月的结束时的库存量1kv必须非负，即0kkkxdv所以kkkdvx从边界条件出发，利用上面的递推关系式，对每个k，计算出)(kkvf中的kv在0至61,minkjkjdmd之间的值，最后求得的)0(6f即为所求的最小总费用。4.1.3模型求解我们假设每件产品的单位工时费为10元，每件产品每月的存贮费为20元，因为产品的生产量与成本费成反比，设反比系数为S，若生产量为X，则成本费为S/X元，设反比系数S为840，则带入模型得到：6161)(7.03.0minkkkkkkchpakx273223251623142011156654321543214321321211xxxxxxxxxxxxxxxxxxxxx用Lingo软件求得的结果如下：Localoptimalsolutionfound.Objectivevalue:3725.200Extendedsolversteps:11Totalsolveriterations:175VariableValueReducedCostX17.000000138.3000X24.000000150.6500X35.000000113.3800X44.000000115.3500X53.00000094.66662X64.00000036.24998RowSlackorSurplusDualPrice13725.200-1.00000021.0000000.00000038.0000000.00000040.0000000.00000059.0000000.00000062.0000000.00000077.0000000.00000084.0000000.00000095.0000000.00000010-0.6380800E-080.000000119.0000000.00000012-0.2648488E-070.000000根据上述结果，我们计算得到的最小成本费为5158元，总耗费工时数最少为382小时。4.2允许缺货的存贮模型4.1.1模型假设1、生产力大于需求量，允许缺货，但是缺货数量在下次生产时补足，每件产品每月的缺货损失费为kq，每个月允许缺货件数为jh，j=1,2,...,62、每件产品每月的存贮费为常数3、每件产品单位工时的成本费为常数4、每个月的生产量不为零4.1.2模型建立设kd为第k个月该产品的需求量，kx为第k个月该产品的生产量，kv为第k个月结束时的产品库存量，则有kkkkdxvv1。kc为第k个月生产产品kx时的成本费用，kh表示第k个月结束时有库存量kv所需的库存费用，kp为第k个月的工时费，kq表示每件产品每月的缺货损失费，故第k个月的总成本费用为kkkkqphc。每个月的缺货量不能超过本月的需求量，并且最后一个月应该满足其需求，即最后一个月的缺货数量为零。由于如果这个月缺货，那么在下个月必须补齐，故每个月的存储量应为本月生产量加上本月的缺货数量在减去上个月的缺货数量。而存储量要受库存限制并且要大于等于零。因此，在考虑缺货条件下的数学模型为)6,,2,1(,),,2,1(,0)5,,2(,0)(0,23.0)(7.0min1606161kZxnkmxkdxvvvakxqphcgkkkjjjkkkkkkkk用动态规划方法来求解，把它看作一个6阶段决策问题。令1kv为状态量，它表示第k个月开始开始时的库存量。kx为决策变量，它表示第k个月的生产量。状态转移方程为)6,,2,1(,1kdxvvkkkk最优值函数)(kkvf表示从第1个月初始库存量为2到第k个月末库存量为kv时的最小总费用。因此可写出顺序递推关系式为)6,,1(,)(min)(110kvfqphcvfkkkkkkxkkkk其中),min(mdvkkk。这是因为一方面每个月生产的上限为m；另一方面由于保证供应，故第k-1个月的结束时的库存量1kv必须非负，即0kkkxdv所以kkkdvx从边界条件出发，利用上面的递推关系式，对每个k，计算出)(kkvf中的kv在0至61,minkjkjdmd之间的值，最后求得的)0(6f即为所求的最小总费用。4.1.3模型求解我们假设每件产品的单位工时费为10元，每件产品每月的存贮费为20元，每件产品每月的缺货损失费为5元，因为产品的生产量与成本费成反比，设反比系数为S，若生产量为X，则成本费为S/X元，设反比系数S为840，则带入模型得到：61613.0)(7.0minkkkkkkkakxqphcg0723566,5,4,3,2,1,013016011012014015065432156645534423312211qqqqqqjqqqxqqxqqxqqxqqxqxj用Lingo软件求得的结果如下：Localoptimalsolutionfound.Objectivevalue:3580.900Extendedsolversteps:9Totalsolveriterations:228VariableValueReducedCostX16.000000116.4667X23.000000122.0666X34.000000100.1500X43.00000086.76662X53.00000094.66662X68.00000063.81250Q10.0000000.000000Q22.00000017.50000Q31.00000017.50000Q40.0000000.000000Q54.00000017.50000Q60.0000000.