生产与库存的动态规划模型

课程设计(论文)题目名称生产与库存的动态规划模型课程名称数学模型学生姓名黄初学号0940802016系、专业理学系信息与计算科学指导教师杜超雄2011年12月18日邵阳学院课程设计（论文）任务书年级专业09信息与计算科学学生姓名黄初学号0940802016题目名称生产与库存的动态规划模型设计时间2011.12.1-12.16课程名称数学模型课程编号080801201设计地点实验室2210一、课程设计（论文）目的1.利用我们所学的数学模型的知识对工厂的生产与库存的问题进行分析，模型的假设，模型的说明，及模型的建立和求解，2.对模型进行评价，对模型的优点及缺点分析，对优点加以利用到实际中，对缺点改进，使模型更加完美。二、已知技术参数和条件1.电子计算机。2．LINDO软件或者LINGO软件。三、任务和要求任务：1.熟悉数学模型书中所学的经典模型。2.理解数学知识在实际生活的广泛应用。要求：1.巩固和加深对数学模型基本知识的理解，提高综合运用课程知识的能力。2.培养学生自学参考书籍，查阅手册、图表和文献资料的能力。3.通过实际课程设计，初步掌握简单软件的分析方法和设计方法。4.了解与课程有关的数学模型，能正确解释和分析实验结果。5.题目具有足够的工作量。注：1．此表由指导教师填写，经系、教研室审批，指导教师、学生签字后生效；2．此表1式3份，学生、指导教师、教研室各1份。四、参考资料和现有基础条件（包括实验室、主要仪器设备等）[1]谢金星等，《数学模型》第三版，高等教育出版社，2003[2]胡运权等，《运筹学基础及应用》第五版，高等教育出版社，2008[3]柳振航等，《数学建模》第一版，中国人民大学出版社2004五、进度安排2011年11月25日-30日：收集和课程设计有关的资料，熟悉课题任务何要求2011年12月1日-4日：总体方案设计2011年12月5日-8日：数学模型设计2011年12月9日-12日：模型设计与程序调试2011年12月13日-16日：整理书写设计说明书六、教研室审批意见教研室主任（签字）：年月日七|、主管教学主任意见主管主任（签字）：年月日八、备注指导教师（签字）：学生（签字）：邵阳学院课程设计（论文）评阅表学生姓名黄初学号0940802016系理学系专业班级信息与计算科学题目名称生产与库存的动态规划模型课程名称数学模型一、学生自我总结通过这次课程设计，让我更好地理解并掌握所学课程的理论知识，同时也锻炼了自己的实践能力，分析问题，解决问题的能力，在课程设计的过程中，老师精心的指导，让我们很快就明确了题目，然后我们查阅资料，确定设计方案，充分锻炼了我们实践能力。在设计的过程中遇到问题，可以说得是困难重重，这毕竟第一次做的，难免会遇到过各种各样的问题，同时在设计的过程中发现了自己的不足之处，对以前所学过的知识理解得不够深刻，掌握得不够牢固，但是通过这次课程设计之后，一定把以前所学过的知识重新温故。这次课程设计终于顺利完成了，在设计中遇到了很多编程问题，最后在杜超雄老师的辛勤指导下，终于游逆而解。同时，对给过我帮助的所有同学和各位指导老师再次表示忠心的感谢！学生签名：黄初2011年12月18日二、指导教师评定评分项目资料查阅编写规范基本技能设计能力科学素养工作量综合成绩权重101225301310单项成绩指导教师评语：指导教师（签名）：年月日注：1、本表是学生课程设计（论文）成绩评定的依据，装订在设计说明书（或论文）的“任务书”页后面；2、表中的“评分项目”及“权重”根据各系的考核细则和评分标准确定。生产与库存的动态规划模型摘要本文讨论了关于生产与存储的问题，这是一个多阶段决策的生产问题，就此可建立一个动态规划的数学模型．利用运筹学和计算机的数学软件等相关知识，应用动态规划方法解决了这一问题，达到生产、需求与库存之间的平衡，以及在资源限制条件下的最优化的生产方案．并建立混合整数规划模型用LINDON数学软件进行检验.问题的提出生产与库存最有问题。设某工厂调查了解市场情况，估计在今后四个时期市场对产品的去求见表1表1时期1234需求量2324假定不论在任何时期，生产每批草坪的固定成本费为3(万元)，若不生产，则为零。每单位生产的固定成本费为1(万元)。同时任何一个时期生产能力所允许的最大生产批量为不超过6个单位。有设每时期的每个单位产品库存费为0.5（万元），同时规定在第一期期初几第四期期末均无产品库存。试问，该厂如何安排各个时期的生产与库存，才能使所花的总成本费用最低？符号说明生产过程划分为四个阶段，阶段变量.4,3,2,1k即：1、状态变量ks表示第k阶段末的库存量，由已知得040ss2、决策变量kx表示第k阶段的生产量，kd表示第k阶段的需求量.3、状态转移方程：kkkkdxss1，4、阶段指标函数),(kkkxsv表示第k阶段的总成本，它由两部分构成一部分是第k阶段的生产成本)(kkxc，另一部分是第k阶段的存贮费)(kksh.最优指标函数)(kksf问题重述已知时段k某产品的需求量为kd(k=1,2,……K)，任一时段若生产该产品，需付出生产准备费0c，且生产每单位产品的生产成本为n，若满足本时段需求后有剩余，每时段每单位产品需付出存贮费0h.设每时段最大生产能力为mX，最大存贮量为mI，且第1时段初有库存量0s，试制订产品的生产计划，即每时段的产量，使K个时段的总费用最小.为了通过具体的计算说明解决这问题的方法，现设4k，,21d,32d,23d,43d30c千元，n=1千元/单位，5.00h千元/单位.时期.01s，6mX单位，mI没有给出，视为存贮量不受限制.模型的建立建立模型Ⅰ在提出生产与存贮问题时，忽略生产准备费用，首先考虑到生产、需求与库存之间存在着的平衡关系，这是一个一般的线性规划问题，可假设生产量为1x，2x，3x，4x，由于存贮费用取决于库存量，则记第一、二、三时期末的库存量为1s，2s，3s，由此可以用生产成本与存贮费之和（记作Z）作为问题为目标函数，在已知的第一期期初及第四期期末均无产品库存040ss，得到一个简单的线性规模型：41415.0kkkksxzMin..ts0...,.....6.....42323141413432321211ssxxxxsxssxssxsx此模型可用单纯形法求解，或用数学软件Maple求解,也可将上模型输入LINDON求解，就可得到最优解（略）.注意：这是在忽略生产准备费用时的最优解.建立模型Ⅱ以上用混合整数规划求解过多阶段生产计划，实际上，这是一类典型的动态优化问题，与用变分法建立连续动态优化模型不同的是，多阶段生产计划属于离散动态优化问题，动态规划模型是解决这类问题的有效方法.本文先讨论确定需求下的最优生产计划，并将它转化为典型的动态优化模型——最短路问题，然后研究随机需求下如何求解最优生产计划.由上述数据、假设，可建立一个动态规划的数学模型.由题可知：6................6,....3,2,1,........30......,.........0)(kkkkkkxxxxxckkkssh5.0)(所以：)()(),(kkkkkkkshxcxsv基本方程为：6,min,0)()4,3,2,1,..()(),(min)(00110kkkkkkkkxkkdssfksfxsvsfkk而模型Ⅱ的求解动态规划的寻优方向一般有用逆序算法（反向递归）或顺序算法（正向递归）进行求解.当问题的第一阶段初和第三阶段末的状态方程均已知时，即040ss，可采用两种方法求解.下面用顺序算法求解：为了简化这个多阶段生产计划问题，可以将它从前向后地分解为一个个单时段问题.（1）首先看第一个时期，为使4个时期的总费用最小，对于第一时期期初的存贮量00s，则可由状态转移方程：kkkkdxss1，考虑到1s，在最大生产能力为6mX与第一时期的需求量21d出发，则可能存在的1s的5种情况：当1k时，有.)()(min)(1111111shxcsfx这时状态集合为：.4,3,2,1,0,26,9min0|,6;min0|1111142111为整数且为整数且ssssddssskk下面就各状态分别计算：505.0213)0()2(min)0(11211hcfx，所以21x5.615.0313)1()3(min)1(11311hcfx，所以31x825.0413)2()4(min)2(11411hcfx，所以41x，同理可得：5.9)3(1f，所以51x，11)4(1f，所以61x（2）当2k时，由)()()(min)()()(min)(2221222201122220222222xdsfshxcsfshxcsfxx其中由：6,min222ds，而状态集合是：.3,2,1,0,36,6min0|,6;min0|2222243222为整数且为整数且ssssddssskk下面就各状态分别计算：5.9565.65845.90min)0()0()3()1()0()2()2()0()1()3()0()0(min)3()0()(min)0(122122122122212223022fhcfhcfhcfhcxfhxcfx所以02x，5.1155.75.65.685.55.95.4115.0min)0()0()4()1()1()3()2()1()2()3()1()1()4()1()0(min)4()1()(min)1(122122122122122212224022fhcfhcfhcfhcfhcxfhxcfx所以02x，同理可得：14)5()2()(min)2(212225022xfhxcfx，所以52x5.15)6()3()(min)3(212226022xfhxcfx，所以62x注意：在计算)2(2f和)3(2f时，需要用到)5(1f和)6(1f，由于每个时期的最大生产批量为6单位，故)5(1f和)6(1f没有意义的，就取)6()5(11ff，其余类推.（3）当3k时，由：33323333033()()(min)(33xdsfshxcsfx,其中6,2min33s，而状态集合为：4,3,2,1,0,6,min0|334333为整数且sddsss下面就各状态分别计算：14)2()0()(min)0(323332033xfhxcfx，所以03x；16)3()1()(min)1(323333033xfhxcfx，所以03x或3；5.17)4()2()(min)2(323334033xfhxcfx，所以43x19)5()3()(min)3(323335033xfhxcfx，所以53x5.20)6()4()(min)4(323336033xfhxcfx，所以63x（4）当4k时，因为要求第4时期期末的库存量为0，即为04s，故有：5.201471665.1751945.200min)0()4()1()3()2()2()3()1()4()0(min)4()0()(min)0(3434343434434444044