运筹学课设某企业和用户签订了设备交货合同

1、.课程设计报告课程设计名称运筹课程设计专业电子商务班级130511 班学 生 姓名雷涵博第四组指 导 教师亚君2016 年6月24 日z.课程设计任务书课程设计题目：第二十八题起止日期：设计地点：教室、电子商务中心z.设计任务及日程安排：1、设计任务1.1 通过课程设计可以增强学生解决实际运筹学问题的能力1.2 通过课程设计可以使学生巩固、拓展和深化所学的基础理论、专业理论和知识。1.3 通过课程设计可以使学生初步掌握用运筹学法解决实际问题的过程和技巧，树立理论联系实际的工作作风。1.4 通过课程设计可以使学生初步建立正确的设计思想和法，进一步提高运算、计算机应用技能和综合分析、解决问题的能力

2、。2 、设计进度安排本课程设计时间分为两：第一（ 2016 年 6 月 13 日-2016 年 6 月 17 日）：建模阶段。此阶段各小组根据给出的题目完成模型的建立。主要环节包括：1.16 月 13 日上午：发指导书；按组布置设计题目；说明进度安排。1.2 6 月 13 日下午至 6 月 15 日：各小组审题，查阅资料，进行建模前的必要准备 (包括求解程序的编写与查找 )。1.3 6 月 16 日至 6 月 17 日：各个小组进行建模， 并根据题目及设计要求拟定设计提纲，指导教师审阅；同时阅读，理解求解程序，为上机求解做好准备。第二（ 2016 年 6 月 20 日 -6 月 24 日）：上

3、机求解，结果分析及答辩。主要环节包括1.16 月 20 日至 6 月 21 日：上机调试程序1.26 月 22 日：完成计算机求解与结果分析。1.36 月 23 日：撰写设计报告。1.46 月 24 日：设计答辩及成绩评定。z.运筹学课程设计报告组别：第四组题号：28题设计人员 ： 黄灵洁黄晓娜雷涵博设计时间 ： 2016 年 6 月 13日至 2016 年 6 月24 日1.设计进度计划第一（ 2016 年 6 月 13 日-2016 年 6 月 17 日）：建模阶段。此阶段各小组根据给出的题目完成模型的建立。主要环节包括：1.16 月 13 日上午：发指导书；按组布置设计题目；说明进度安排

4、。1.2 6 月 13 日下午至 6 月 15 日：各小组审题，查阅资料，进行建模前的必要准备 (包括求解程序的编写与查找 )。1.3 6 月 16 日至 6 月 17 日：各个小组进行建模， 并根据题目及设计要求拟定设计提纲，指导教师审阅；同时阅读，理解求解程序，为上机求解做好准备。第二（ 2016 年 6 月 20 日 -6 月 24 日）：上机求解，结果分析及答辩。主要环节包括1.16 月 20 日至 6 月 21 日：上机调试程序1.26 月 22 日：完成计算机求解与结果分析。1.36 月 23 日：撰写设计报告。1.46 月 24 日：设计答辩及成绩评定。2.设计题目二十八、某企业

5、和用户签订了设备交货合同，已知该企业各季度的生产能力、每台设备的生产成本和每季度末的交货量（见下表），若生产出的设备当季度不交货，每台设备每季度需支付保管维护费 0.2 万元，试问在遵守合同的条件下，企业应如安排生产计划，才能使年消耗费用最低？并按要求分别完成下列分析：（ 1）2 季度每台设备的生产成本在围变化时最有生产计划不变？（ 2）每台设备每季度需支付保管维护费在围变化时最优生产计划不变？（ 3） 1 季度生产能力在围变化变化时最优基不变？（ 4） 4 季度交货量在围变化时最优基不变？z.季度工厂生产能力交货量每台设备生产成本（万元）1251512.02352011.03302511.5

6、4202012.53.建模3.1 题目分析 ,变量设定这个优化问题的目标是使年消耗费用最低，因此此问题需要做的是有关生产决策的灵敏度分析问题，其受到三个约束条件的限制：第 1、 2、 3 、 4季度工厂生产能力；各个季度工厂交货量与要求的交货量相等；生产的产品个数应该为非负整数。所以变量设定如下：Z ：年消耗费用决策变量 Xij ： 表示第 i 季度生产第 j 季度交货的设备数量3.2 建模分析目标函数分析：第1季度生产的消耗费用+ 第 2 季度的消耗费用+ 第 3 季度的消耗费用+第 3 季度的消耗费用 + 第4 季度的消耗费用由生产能力所限列出下列式子：第一季度生产能力所限：X11+X 1

7、2+X 13+X 1425第二季度生产能力所限：X22+X 23 +X 24 35第三季度生产能力所限：X33+X 34 30第四季度生产能力所限： X4420得到前四个约束变量由交货量所限列出下列式子：第一季度交货量：X11=15第二季度交货量：X12+X 22=20第三季度交货量：X13+X 23+X 33 =25第四季度交货量：X14+X 24+X 34 +X 44 =20得到后四个约束变量第 i 季度生产第 j 季度交货的每台设备所消耗的费用Cij 应等于生产成本加上保管维护费z.用之和其值如下表i1234112.012.212.412.6211.011.211.4311.511.74

8、12.53.3 数学模型用 Ai 表示该企业第 i季度的生产能力， Bj 表示第 j季度的交货量，则可将这一问题的数学模型写成：Min Z = 12.0X 11 + 12.2X 12 + 12.4X 13 + 12.6X 14 + 11.0X 22 + 11.2X 23 + 11.4X 24 + 11.5X 33 + 11.7X 34 + 12.5X 44X 11+X 12 +X 13+X 1425X 22 +X 23 +X 24 35X 33 +X 34 30X 44 20X 11=15X 12+X 22 =20X 13+X 23 +X 33 =25X 14+X 24 +X 34 +X 44

9、 =20X 11,X12 ,X13 ,X14 ,X22 ,X23 ,X24 ,X33 ,X34 ,X44 0z.4.程序4.1 求解程序流程图开始根据题目设定变量X 和约束条件，方程组系数矩阵A、约束条件常数项b 值、选择目标函数类型、目标函数系数化成标准形式：调整目标函数为max Z ，加入松弛变量、剩余变量和人工变量，构造人造基。根据人工变量构造辅助LP 问题 max w，构造单纯形表进行迭是否存在检验数所是否在列中有 pj0是否所有检验数j 0否无解解是输出最优表和最优解、最优值进行比较选择，确定主元和换入、换出变量换基迭代W 是否=0否无 解是是人工变量所在行原是人工变量所在行原人工变

10、量是否为基变量始变量系数不全为始变量系数全为零否以非零系数其中之一为主元进行换基迭代， 把人工变量变为非删去人工变量诸列，用Z 代替 w，用单纯形法求解删去相应行基变量得出最优值Zb 发生变化 ,确定 b 的变化范围否-1-1C,b 是否在变化范围之内重新计算 CB Bb、Bb 代入最优表中重新迭代。是。最优解（基）不变，用C B B-1 b 计算出现在的最优值得出最优值并与原最优值进行比较z.结束4.2 求解程序功能介绍Java 是一种简单的，跨平台的，面向对象的，分布式的，解释的，健壮的，安全的，结构的中立的，可移植的，性能很优异的多线程的，动态的语言。我们用java 语言设计编辑了一个解

11、题程序来解答这个问题，运用 eclipse 环境运行的， 我们首先程序要求输入目标函数类型以确定使用哪种程序解决当前问题，然后函数要求输入约束条件的个数和变量个数，以用来确定所创建的数组，然后函数要求小于等于、大于等于、和等于的条件个数，同样用来创建数组和确定解题程序，然后函数要求输入题的系数矩阵，函数记录下来用于解题，当输入结束时函数会输出刚刚所录入的系数矩阵，以便让输入者确定所输入的矩阵是否正确，避免错误的产生，此时进行到程序的最后一步，函数要求输入目标函数系数，输入者正确输入后函数就会运算然后输出结果LINGO 的注意事项：1. Lingo 中不能省略乘号 * 以及结束符分号 ;2. 目

12、标函数在 Lingo 写成 max=2*x+3*y;3. Lingo 中所用的符号均必须是英文状态下的，使用中文状态下的这些符号，运行会报错的。4. Lingo 默认变量非负。5. 即有约束条件 X 0,Y ij 0 之类的，在 Lingo 中都可以省略不写。6. 有需要变量取负值的话，可以用free(); 函数，令变量取全体实数。7. Lingo 中没有这个常数。在要求不是那么精确的话可以用PI=3.1415926 （取几位小数自行估计）。通过三角函数来得到精确的值，如PI=acos(-1);8.Lingo 中没有格大或格小的概念（除集下标约束或条件判断中的#gt#,#lt# ）。在 Lin

13、go 中使用 (=(=5+0.000001;10. 变量为某几个不连续的数值之一，例如 x 为 0 或 3 或 7 ，可表示为 x*(x-3)*(x-7)=0;z.4.3 手工数据准备Java 程序的录入界面：Lingo 录入界面：z.5.结果分析5.1 结果分析思路季度每台设备的生产成本在围变化时最优生产计划不变？此问题为目标函数系数Ci 的变化围的处理：如果变化的系数为非基变量系数，确定非基变量系数变化围，非基变量系数变化只影响自身的检验数，因此，设 Ci 为非基变量 Xi 的系数，令它在当前最优表中的检验数 i =C BB-1Pi-Ci 0,当 Ci 发生了 Ci 变化后，要保证当前最优

14、表中相应的检验数仍大于或等于0 ，即当 Ci i 均满足时，最优生产计划不变。如果变化的系数为基变量的系数，则要确定基变量系数的变化围，基变量系数的变化影响所有非基变量的检验数和目标函数值。Ci 的变化围在基变量Ci 增量的变化围之，则其最优解、最优值均不变，即要求基变量系数满足一定的条件，即当max i /b ri |b ri 0 Ci min i/ b ri |b ri 0 ，其中 i 取非基变量检验数， 此时，最优生产计划不变。每台设备每季度需支付保管维护费在围变化时最优生产计划不变？此问题为目标函数系数Ci 的变化围的处理：如果变化的系数为非基变量系数，确定非基变量系数变化围，非基变量

15、系数变化只影响自身的检验数， 因此，设 Ci 为非基变量 Xi 的系数， 令它在当前最优表中的检验数 i=C B B-1Pi -Ci 0, 当 Ci 发生了 Ci 变化后，要保证当前最优表中相应的检验数仍大于或等于 0 ，即当 Ci i 均满足时，最优生产计划不变。z.如果变化的系数为基变量的系数，则要确定基变量系数的变化围，基变量系数的变化影响所有非基变量的检验数和目标函数值。Ci 的变化围在基变量Ci 增量的变化围之，则其最优解、最优值均不变，即要求基变量系数满足一定的条件，即当max i /b ri |b ri 0 Ci min i/ b ri |b ri 0 ，其中 i 取非基变量检验

16、数， 此时，最优生产计划不变。季度生产能力在围变化变化时最优基不变？此问题为约束条件常数项系数b 变化围的处理：根据公式max-b i/ ir | ir 0 b min-b i/ ir | ir 0 b min-b i/ ir | ir 0 ，确定 b 的变化围，如果 b 的变化在该围之最优基不变最优解变化，最优解由公式XB =B -1b 求得；如果b 的变化超出该围最优基最优解均发生变化，重新计算CB B-1b、 B-1b 代入最优表中重新迭代。5.2 求解结果截图结果z.Global optimal solution found.Objective value:917.0000Infeas

17、ibilities:0.000000Total solver iterations:4Elapsed runtime seconds:0.04Model Class:LPTotal variables:9Nonlinear variables:0Integer variables:0Total constraints:8Nonlinear constraints:0Total nonzeros:27Nonlinear nonzeros:0VariableValueReduced CostX1115.000000.000000X120.0000000.9000000X130.0000000.90

18、00000X140.0000000.9000000X2220.000000.000000X2315.000000.000000X240.0000000.000000X3310.000000.000000X3420.000000.000000X440.0000000.8000000RowSlack or SurplusDual Price1917.0000-1.000000210.000000.00000030.0000000.300000040.0000000.000000520.000000.00000060.000000-12.0000070.000000-11.3000080.00000

19、0-11.50000z.90.000000-11.700005.3 灵敏度分析结果Ranges in which the basis is unchanged:Objective Coefficient Ranges:CurrentAllowableAllowableVariableCoefficientIncreaseDecreaseX1212.20000INFINITY0.9000000X1312.40000INFINITY0.9000000X1412.60000INFINITY0.9000000X2211.000000.9000000INFINITYX2311.200000.000000

20、0.9000000X2411.40000INFINITY0.000000X3311.500000.90000000.000000X3411.700000.000000INFINITYX4412.50000INFINITY0.8000000Righthand Side Ranges:CurrentAllowableAllowableRowRHSIncreaseDecrease210.00000INFINITY10.00000335.0000010.000000.000000430.00000INFINITY0.000000520.00000INFINITY20.00000720.000000.0

21、0000010.00000825.000000.00000010.00000920.000000.00000020.000005.4 结果分析季度每台设备的生产成本在围变化时最有生产计划不变？即 X22 、 X23 、 X24 的系数变化时最优生产计划不变时的系数变化围：z.VariableCoefficientIncreaseDecreaseX2211.000000.9000000INFINITYX2311.200000.0000000.9000000X2411.40000INFINITY0.000000由 lingo 灵敏度分析表得： 2 季度每台设备的生产成本在 11.2,11.4围变化

22、时最优生产计划不变。每台设备每季度需支付保管维护费在围变化时最优生产计划不变？根据 lingo灵敏度分析表VariableCoefficientIncreaseDecreaseX1212.20000INFINITY0.9000000X1312.40000INFINITY0.9000000X1412.60000INFINITY0.9000000X2211.000000.9000000INFINITYX2311.200000.0000000.9000000X2411.40000INFINITY0.000000X3311.500000.90000000.000000X3411.700000.0000

23、00INFINITYX4412.50000INFINITY0.8000000C12 变化围 11.3,+ C13 变化围 11.5,+ C14 变化围 11.7,+ C22 变化围 - ,11.9C23 变化围 11.2,+ C24 变化围 11.4,+ C33 变化围 11.5,12.4C34 变化围 - ,11.7C44 变化围 11.7,+ 季度生产能力在围变化变化时最优基不变？即约束条件右端常数项的变化对最优解的变化，第一个约束条件右端常数项的变化对最有基的影响。RowRHSIncreaseDecrease210.00000INFINITY10.00000由分析结果得：0,+ 。季度交

24、货量在围变化时最优基不变？即约束条件右端常数项的变化对最优解的变化，第八个约束条件右端常数项的变化对最有基的影响。z.RowRHSIncreaseDecrease920.000000.00000020.00000由分析结果得：4 季度交货量在 0,20 的围变化时最优基不变。6.创新容6.1 2 季度交货量为 30 时的最优生产计划？由灵敏度表分析可得：RowRHSIncreaseDecrease720.000000.00000010.00000第五个约束条件右端常数项的变化围为10 到 20 ，超出围，输入程序重新计算得：最优案如上图所示。6.2 第 3 季度生产第4 季度交货的每台设备所消

25、耗的费用C34 在什么围变化时最优生产计划不变？目标函数X 34 的系数 C34 的变化围，由lingo 灵敏度分析可得：VariableCoefficientIncreaseDecreaseX3411.700000.000000INFINITYC34 的变化围为：- ,11.7。6.3 2 季度生产能力在围变化变化时最优基不变？RowRHSIncreaseDecrease335.0000010.000000.000000根据第二约束条件可知2季度生产能力在35,45 之最优基不变。z.7.课程设计总结学习理论的目的就是为了解决实际问题。通过两的课程设计，让我对运筹学的认识更加深刻，也知道了运

26、筹学的在实践中的重要作用。在之前的运筹学课程学习中，一直都很认真的在学习每个算法，掌握了解体的法。运筹学的学习的最终目的是运用量化和模型化法，针对有有关管理问题最优化配置的现实需求，来寻找系统的优化的解决案，系统求得一个合理应用人力，节省财力的最佳案，为决策者提供科学的决策依据，而不是平时学习中的简单的数据计算和公式运用。而在实际问题的解决中，建立一个合理实际的模型尤为困难，也是解决问题的关键。只有建立了正确的模型，解决问题才能水到渠成。在开始看到这道题目时，大家的首要任务就是讨论如建立正确的数学模型，找到正确的法来解决问题。建立了正确的数学模型之后，才能快速的编程。作为科学的数学法， 它与计算机技术结合解决实际问题。 在这次实践中我意识到，我们应该尽可能多的将先进的管理中的科学法与先进的计算机技术结合起来，切实提高我们适应当代社会信息、运用运筹学解决实际问题的能力，我想这才是我们这次课程设计所要达到的最终目的。在现实生活中，我们如果能运用科学的法和最简单的办法来处理遇到的问题，那就会大大提高效率。运筹学就是这样一个为我们提供简便的算法来分析一些生活中的问