生产策略问题数学实验

一、实验题目：生产策略问题二、实验内容：问题重述现代化生产过程中，生产部门面临的突出问题之一，便是如何选取合理的生产率。生产率过高，导致产品大量积压，使流动资金不能及时回笼；生产率过低，产品不能满足市场需要，使生产部门失去获利的机会。可见，生产部门在生产过程中必须时刻注意市场需求的变化，以便适时调整生产率，获取最大收益。某生产厂家年初要制定生产策略，已预知其产品在年初的需求量为a=6万单位，并以b=1万单位/月速度递增。若生产产品过剩，则需付单位产品单位时间（月）的库存保管费C2=0.2元；若产品短缺，则单位产品单位时间的短期损失费C3=0.4元。假定生产率每调整一次带有固定的调整费C1=1万元，试问工厂如何制定当年的生产策略，使工厂的总损失最小?三、数学模型：生产率过高，导致产品大量积压，使流动资金不能及时回笼；生产率过低，产品不能满足市场需要，使生产部门失去获利的机会。可见，为使工厂的总损失最少，生产部门在生产过程中必须时刻注意市场需求的变化，从而制定出使工厂总损失最小的生产策略。文章把此求工厂总损失最小生产策略问题化为最短路问题的多阶段决策问题。设每个顶点代表各月，且以每个顶点为转折点进行生产策略调整，求出每个阶段的最小损耗。最后，使用Matlab软件求出最短的路径，此路径即为使工厂损失最小的生产策略。每月社会需求量见下表：月份123456789101112需求（万元）67891011121314151617四、模型假设与符号说明：1、市场的需求量严格按照年初的需求量为a=6万单位，并以b=1万单位/月速度递增；2、单位产品单位时间的库存保管费、短期损失费以及生产率每调整一次带有固定的调整费均不变；3、工厂可以严格按照生产率生产产品。符号说明顶点xx1211月至12月初；顶点x1312月末；弧xxaii从i月至1ai月不调整生产策略，111,212iai；sxxaii从i月至1ai月库存保管费和短期损失费的最小值以及第ai月的调整费用之和，111,212iai;sxxi13从i月至12月库存保管费和短期损失费的最小值,111i；s工厂一年的总损失；X不调整前每月生产X万单位；Yii月库存保管费和短期损失费；五、算法步骤：Floyd算法：求任意两点间的最短路．D(i,j)：i到j的距离．R(i,j)：i到j之间的插入点．输入:带权邻接矩阵w(i,j)赋初值：对所有i,j,d(i,j)w(i,j),r(i,j)j,k1(2)更新d(i,j),r(i,j)对所有i,j，若d(i,k)+d(k,j)d(i,j)，则d(i,j)d(i,k)+d(k,j),r(i,j)k(3)若k=，停止．否则kk+1，转（２）．六、算法编程：附表1：function[D,R]=floyd(a)n=size(a,1);D=afori=1:nforj=1:nR(i,j)=j;endendRfork=1:nfori=1:nforj=1:nifD(i,k)+D(k,j)D(i,j)D(i,j)=D(i,k)+D(k,j);R(i,j)=R(i,k);endendendkDRend附表2：w=[011.11.4234.25.87.79.912.41517;1011.11.4234.25.87.79.912.414;1.11011.11.4234.25.87.79.911.4;1.41.11011.11.4234.25.87.78.9;21.41.11011.11.4234.25.86.7;321.41.11011.11.4234.24.8;4.2321.41.11011.11.4233.2;5.84.2321.41.11011.11.422;7.75.84.2321.41.11011.11.41;9.97.75.84.2321.41.11011.10.4;12.49.97.75.84.2321.41.11010.1;1512.49.97.75.84.2321.41.1100;171411.48.96.74.83.2210.40.100]七、实验过程及结果：7.1计算1月的库存保管费和短期损失费的最小值0以及2月的调整费用1万，因此为最小损耗sxx21为1（万元）。同理，可得sxxii1（111i）皆为1（万元），sxx1312为0。7.2计算1月至2月的库存保管费和短期损失费的最小值以及3月的调整费用1万最小值计算（1）6=X6.5Y1=（X-6)*0.2Y2=(13-2X)*0.4S=(4-0.6*X)+1(2)X=6.5Y1=（X-6)*0.2Y2=(2X-13)*0.2S=(0.6X-3.8)+1当X=6.5，sxx31为1.1（万元）。同理，可得sxxii2（110i）皆为1.1（万元），sxx1311为0.1（万元）从上式我们可以看出不论在何种情况下，因Yi是一次函数，而sxxaii为Yi的和加1（除1月至12月），所以sxxaii也为一次函数，所以最小损耗必在端点处取值。7.3计算1月至3月的库存保管费和短期损失费的最小值以及4月的调整费用1万分X=7，6.5=X7，6=X6.5三种情况讨论；得X=7，因此sxx41为1.4（万元）。同理，可得sxxii3（19i）皆为1.4（万元），sxx1310为0.4（万元）。7.4计算1月至4月的库存保管费和短期损失费的最小值以及5月的调整费用1万分X=7.5，7=X7.5,6.5=X7，6=X6.5四种情况讨论；得X=7.5，因此sxx51为2（万元）。同理，可得sxxii4（18i）皆为2（万元），sxx139为1（万元）。7.5计算1月至5月的库存保管费和短期损失费的最小值以及6月的调整费用1万分X=8，7.5=X8,7=X7.5,6.5=X7，6=X6.5五种情况讨论；得X=8，因此sxx61为3（万元）。同理，可得sxxii5（17i）皆为3（万元），sxx138为2（万元）。7.6计算1月至6月的库存保管费和短期损失费的最小值以及7月的调整费用1万分X=8.5,8=X8.5,7.5=X8,7=X7.5,6.5=X7，6=X6.5六种情况讨论；得X=8，因此sxx71为4.2（万元）。同理，可得sxxii6（16i）皆为4.2（万元），sxx137为3.2（万元）。7.7计算1月至7月的库存保管费和短期损失费的最小值以及8月的调整费用1万分X=9,8.5=X9,8=X8.5,7.5=X8,7=X7.5,6.5=X7，6=X6.5七种情况讨论；得X=8，因此sxx81为5.8（万元）。同理，可得sxxii7（15i）皆为5.8（万元），sxx136为4.8（万元）。7.8计算1月至8月的库存保管费和短期损失费的最小值以及9月的调整费用1万，分X=9.5,9=X9.5,8.5=X9,8=X8.5,7.5=X8,7=X7.5,6.5=X7，6=X6.5八种情况讨论；得X=8.5，因此sxx91为7.7.（万元）。同理，可得sxxii8（14i）皆为7.7（万元），sxx135为6.7（万元）。7.9计算1月至9月的库存保管费和短期损失费的最小值以及10月的调整费用1万，分X=10,9.5=X10,9=X9.5,8.5=X9,8=X8.5,7.5=X8,7=X7.5,6.5=X7，6=X6.5九种情况讨论；得X=8.5，因此sxx101为9.9（万元）。同理，可得sxxii9（13i）皆为9.9（万元），sxx134为8.9（万元）。7.10计算1月至10月的库存保管费和短期损失费的最小值以及11月的调整费用1万,分X=10.5,10=X10.5，9.5=X10,9=X9.5,8.5=X9,8=X8.5,7.5=X8,7=X7.5,6.5=X7，6=X6.5十种情况讨论；得X=9.5，因此sxx111为12.4（万元）。同理，可得sxxii10（12i）皆为12.4（万元），sxx133为11.4（万元）。7.11计算1月至11月的库存保管费和短期损失费的最小值以及12月的调整费用1万，分X=11,10.5=X11，10=X10.5，9.5=X10,9=X9.5,8.5=X9,8=X8.5,7.5=X8,7=X7.5,6.5=X7，6=X6.5十一种情况讨论；得X=9.5，因此sxx121为15（万元）。sxx132为14（万元）。7.12计算1月至12月的库存保管费和短期损失费的最小值，分X=11.5,11=X11.5，10.5=X11，10=X10.5，9.5=X10,9=X9.5,8.5=X9,8=X8.5,7.5=X8,7=X7.5,6.5=X7，6=X6.5十二种情况讨论；得X=9.5，sxx131=17万。总权值表:调整三次，四月初七月初十月初各调整一次，s=1.4*4-1=4.6万元。1—3月，产量为7万单位每月；4—6月，产量为10万单位每月，7—9月，产量为13万单位每月；10—12月，产量为16万单位每月。七、实验结论：把此求最少损耗的问题转化为最短路径的多阶段问题非常形象，让人容易理解。在计算出最低损耗的同时也表示出了最短损耗的路径，可以清楚的得出工厂生产的策略。