生产的策略规划

组员：韩河江生产策略问题一问题现代化生产过程中，生产部门面临的突出问题之一，便是如何选取合理的生产率。生产率过高，导致产品大量积压，使流动资金不能及时回笼；生产率过低，产品不能满足市场需要，使生产部门失去获利的机会。可见，生产部门在生产过程中必须时刻注意市场需求的变化，以便适时调整生产率，获取最大收益。某生产厂家年初要制定生产策略，已预知其产品在年初的需求量为a=6万单位，并以b=1万单位/月速度递增。若生产产品过剩，则需付单位产品单位时间（月）的库存保管费C2=0.2元；若产品短缺，则单位产品单位时间的短期损失费C3=0.4元。假定生产率每调整一次带有固定的调整费C1=1万元，试问工厂如何制定当年的生产策略，使工厂的总损失最小？二分析生产率过高，导致产品大量积压，使流动资金不能及时回笼；生产率过低，产品不能满足市场需要，使生产部门失去获利的机会。可见，为使工厂的总损失最少，生产部门在生产过程中必须时刻注意市场需求的变化，从而制定出使工厂总损失最小的生产策略。我们可把此求工厂总损失最小生产策略问题化为最短路问题的多阶段决策问题。设每个顶点代表各月，且以每个顶点为转折点进行生产策略调整，求出每个阶段的最小损耗，最后，使用Matlab软件求出最短的路径，此路径即为使工厂损失最小的生产策略。三假设3.1市场的需求量严格按照年初的需求量为a=6万单位，并以b=1万单位/月速度递增。3.2单位产品单位时间的库存保管费、短期损失费以及生产率每调整一次带有固定的调整费均不变。3.3工厂可以严格按照生产率生产产品。四分析与建模把此求工厂总损失最小生产策略问题化为最短路问题的多阶段决策问题，计算各阶段的最小损耗，及为它们之间的权值。符号说明符号说明顶点1月至12月初；顶点12月末；弧从月至月不调整生产策略，；从月至月库存保管费和短期损失费的最小值以及第月的调整费用之和，;从月至12月库存保管费和短期损失费的最小值,；工厂一年的总损失；X不调整前每月生产X万单位；Yii月库存保管费和短期损失费；每月社会需求量见下表：月份123456789101112需求量（万元）678910111213141516175.1计算1月的库存保管费和短期损失费的最小值0以及2月的调整费用1万，因此为最小损耗为1（万元）。同理，可得（）皆为1（万元），为0。5.2计算1月至2月的库存保管费和短期损失费的最小值以及3月的调整费用1万最小值计算（1）6=X=6.5Y1=（X-6)*0.2Y2=(2X-13)*0.2S=(0.6X-3.8)+1当X=6.5，因此为1.1（万元）。同理，可得（）皆为1.1（万元），为0.1（万元）从上式我们可以看出不论在何种情况下，因Yi是一次函数，而为Yi的和加1（除1月至12月），所以也为一次函数，所以最小损耗必在端点处取值。5.3计算1月至3月的库存保管费和短期损失费的最小值以及4月的调整费用1万分X=7，6.5=X7，6=X=7.5，7=X7.5,6.5=X7，6=X=8，7.5=X8,7=X7.5,6.5=X7，6=X=8.5,8=X8.5,7.5=X8,7=X7.5,6.5=X7，6=X=9,8.5=X9,8=X8.5,7.5=X8,7=X7.5,6.5=X7，6=X=9.5,9=X9.5,8.5=X9,8=X8.5,7.5=X8,7=X7.5,6.5=X7，6=X=10,9.5=X10,9=X9.5,8.5=X9,8=X8.5,7.5=X8,7=X7.5,6.5=X7，6=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=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=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=Xl(u)+w(u,v),则更新l(v),f(v),即l(v)?l(u)+w(u,v),f(v)?u;（3）重复步骤2),直到所有顶点都在S中为止.MATLAB程序（Dijkstra算法）见附表1：MATLAB求解程序见附表2：六结论调整三次，四月初七月初十月初各调整一次，s=1.4*4-1=4.6万元。13月，产量为7万单位每月；46月，产量为10万单位每月，79月，产量为13万单位每月；1012月，产量为16万单位每月。附表1；functionmin,path=dijkstra(w,start,terminal)n=size(w,1);label(start)=0;f(start)=start;fori=1:nifi=startlabel(i)=inf;end,ends(1)=start;u=start;whilelength(s)(label(u)+w(u,v)label(v)=(label(u)+w(u,v);f(v)=u;end,end,endv1=0;k=inf;fori=1:nins=0;forj=1:length(s)ifi=s(j)ins=1;end,endifins=0v=i;ifklabel(v)k=label(v);v1=v;end,end,ends(length(s)+1)=v1;u=v1;endmin=label(terminal);path(1)=terminal;i=1;whilepath(i)=startpath(i+1)=f(path(i);i=i+1;endpath(i)=start;L=length(path);path=path(L:-1:1);附表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.100dis,path=dijkstra(w,1,13)