生产策略问题-数学实验

1、一、实验题目：生产策略问题二、实验内容：问题重述现代化生产过程中，生产部门面临的突出问题之一，便是如何选取合理的生产率。生产率过高，导致产品大量积压，使流动资金不能及时回笼；生产率过低，产品不能满足市场需要，使生产部门失去获利的机会。可见，生产部门在生产过程中必须时刻注意市场需求的变化，以便适时调整生产率，获取最大收益。某生产厂家年初要制定生产策略，已预知其产品在年初的需求量为a=6万单位，并以b=1万单位/月速度递增。若生产产品过剩，则需付单位产品单位时间（月）的库存保管费C2=0.2元；若产品短缺，则单位产品单位时间的短期损失费C3=0.4元。假定生产率每调整一次带有固定的调整费C1=1万

2、元，试问工厂如何制定当年的生产策略，使工厂的总损失最小?三、数学模型：生产率过高，导致产品大量积压，使流动资金不能及时回笼；生产率过低，产品不能满足市场需要，使生产部门失去获利的机会。可见，为使工厂的总损失最少，生产部门在生产过程中必须时刻注意市场需求的变化，从而制定出使工厂总损失最小的生产策略。文章把此求工厂总损失最小生产策略问题化为最短路问题的多阶段决策问题。设每个顶点代表各月，且以每个顶点为转折点进行生产策略调整，求出每个阶段的最小损耗。最后，使用Matlab软件求出最短的路径，此路径即为使工厂-损失最小的生产策略每月社会需求量见下表:月份123456789101112需求（万兀）678

3、91011121314151617四、模型假设与符号说明：1、市场的需求量严格按照年初的需求量为a=6万单位，并以b=1万单位/月速度递增；2、单位产品单位时间的库存保管费、短期损失费以及生产率每调整一次带有固定的调整费均不变；3、工厂可以严格按照生产率生产产品符号说明顶点XX121月至12月初；顶点X1312月末；弧XXia从i月至ia1月不调整生产策略，12ia2,11i1；SXiXia从i月至ia1月库存保管费和短期损失费的最小值以及第ia月的调整费用之和，12ia2,11i1;SXiX13从i月至12月库存保管费和短期损失费的最小值，11i1;s工厂一年的总损失；X不调整前每月生产X万

4、单位；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,k1(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=,停止.否则kk+1,转(2).六、算法编程：附表1：functionD,R=floyd(a)n=size(a,1);D=afori=1:nforj=1:nR(i,j)=j;endendR

5、fork=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;1 011.11.4234.25.87.79.912.414;1.1 1011.11.4234.25.87.79.911.4;1.4 1.11011.11.4234.25.87.78.9;2 1.41.11011.11.4234.25.86.7;3 21.41.11011.11.4234.24.8;4 .2321.41.110

6、11.11.4233.2;5 .84.2321.41.11011.11.422;6 .75.84.2321.41.11011.11.41;7 .97.75.84.2321.41.11011.10.4;8 2.49.97.75.84.2321.41.11010.1;15 12.49.97.75.84.2321.41.1100;17 1411.48.96.74.83.2210.40.100七、实验过程及结果：7.1 计算1月的库存保管费和短期损失费的最小值0以及2月的调整费用1万，因此为最小损耗s为1（万元）。x1x2同理，可得s（11i1）皆为1（万元），s为0。xixi1x12x137.2 计

7、算1月至2月的库存保管费和短期损失费的最小值以及3月的调整费用1万最小值计算（1）6<=X<6.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,sx1x3为1.1（万元）0同理，可得s（10i1）皆为1.1（万元），s为0.1（万元）xixi2x11x13从上式我们可以看出不论在何种情况下，因Yi是一次函数，而s为xixiaYi的和加1（除1月至12月），所以s也为一次函数，所以最小损耗必在xixia端点处取值。7.3计算1月至3月的

8、库存保管费和短期损失费的最小值以及4月的调整费用1万分X>=7，6.5<=X<7，6<=X<6.5三种情况讨论；得X=7，因此s为1.4（万元）。sx1x4同理，可得s（9i1）皆为1.4（万元），s为0.4（万元）。xixi3sx10x137.4 计算1月至4月的库存保管费和短期损失费的最小值以及5月的调整费用1万分X>=7.5，7=<X<7.5,6.5=<X<7，6=<X<6.5四种情况讨论；得X=7.5，因此s为2（万元）。sx1x5同理，可得s（8i1）皆为2（万元），s为1（万元）。xixi4sx9x137.5

9、计算1月至5月的库存保管费和短期损失费的最小值以及6月的调整费用1万分X>=8，7.5<=X<8,7<=X<7.5,6.5<=X<7，6<=X<6.5五种情况讨论；得X=8，因此s为3（万元）。sx1x6同理，可得s（7i1）皆为3（万元），s为2（万元）。xixi5x8x137.6 计算1月至6月的库存保管费和短期损失费的最小值以及7月的调整费用1万分X>=8,5,8<=X<8,5,7.5<=X<8,7<=X<7,5,6.5<=X<7,6<=X<6,5六种情况讨论；得X=8

10、,因此s为4.2（万元）。sXiX7同理，可得s（6i1）皆为4.2（万元），&为3.2（万元）。XiXi6sX7X137.7 计算1月至7月的库存保管费和短期损失费的最小值以及8月的调整费用1万分X>=9,8.5<=X<9,8<=X<8.5,7.5<=X<8,7<=X<7.5,6.5<=X<7,6<=X<6.5七种情况讨论；得X=8,因此s为5.8（万元）。sX1X8同理，可得s（5i1）皆为5.8（万元），s为4.8（万元）。XiXi7sX6X137.8 计算1月至8月的库存保管费和短期损失费的最小值以及

11、9月的调整费用1万，分X>=9.5,9<=X<9.5,8.5<=X<9,8<=X<8.5,7.5<=X<8,7<=X<7.5,6.5<=X<7,6<=X<6.5八种情况讨论；得X=8.5,因此sX1X9为7.7.（万元）。同理，可得sXiXi8（4i1）皆为7.7（万元）,sX5X3为6.7（万元）。7.9 计算1月至9月的库存保管费和短期损失费的最小值以及10月的调整费用1万，分X>=10,9.5=<X<10,9=<X<9.5,8.5=<X<9,8=<X&

12、lt;8.5,7.5=<X<8,7=<X<7.5,6.5=<X<7,6=<X<6.5九种情况讨论；得X=8.5,因此s为9.9（万元）。sXiX10同理，可得s（3i1）皆为9.9（万元）,s为8.9（万元）。wXiXi92X4X137.10 计算1月至10月的库存保管费和短期损失费的最小值以及11月的调整费用1万，分X>=10.5,10=<X<10.5,9.5=<X<10,9=<X<9.5,8.5=<X<9,8=<X<8.5,7.5=<X<8,7=<X<7

13、.5,6.5=<X<7,6=<X<6.5十种情况讨论；得X=9.5,因此sXx“为12.4（万元）。同理，可得s（2i1）皆为12.4（万元）,s为11.4（万元）csXiXi10sX3X137.11 计算1月至11月的库存保管费和短期损失费的最小值以及12月的调整费用1万，分X>=11,10.5=<X<11,10=<X<10.5,9.5=<X<10,9=<X<9.5,8.5=<X<9,8=<X<8.5,7.5=<X<8,7=<X<7.5,6.5=<X<7,

14、6=<X<6.5H一种情况讨论；得X=9.5,因此s为15（万元）。s为14（万元）。sX1X12sX2X137.12 计算1月至12月的库存保管费和短期损失费的最小值，分X>=11.5,11=<X<11.5,10.5=<X<11,10=<X<10.5,9.5=<X<10,9=<X<9.5,8.5=<X<9,8=<X<8.5,7.5=<X<8,7=<X<7.5,6.5=<X<7,6=<X<6.5十二种情况讨论；得X=9.5,s=17万。sX1X1

15、3邻接矩阵zlx2X'l工5:<6>:7工9virzllzl2zl011.11.42R4.25.3T.79.912.41517_工21urij1.11,423_4.25.87.7%92.414工31.11n11.11.424.25.87.79.9LL41,41.1ri二01二L.11.4234.25.87.7工52L4i.i1011.1lg234.2186,7工615J211.4j1工1J011.11.4234,24.814.2321.4：J1011.1L4233,215.3J4.213J21.4ILI1011.11.422_内L75,84.2321.41.11011.1L41kIOJ7.715.3j4.23J21.4Ll101Ll0.4,xll12,49.97.74,232L41.11010.1i121