《运筹课程设计1.doc》由会员分享，可在线阅读，更多相关《运筹课程设计1.doc（11页珍藏版）》请在金锄头文库上搜索。
1、长 春 工 业 大 学课 程 设 计 报 告课程设计名称 运筹课程设计 专 业 信息管理与信息系统 班 级 080505 班 学 生 姓 名 王盛之 指 导 教 师 王亚君 2011年7月8日课 程 设 计 任 务 书课程设计题目：第 一 题起止日期：2011、6、272011、7、8设计地点：教室、电子商务中心设计任务及日程安排：1、设计任务1.1通过课程设计可以增强学生解决实际运筹学问题的能力1.2通过课程设计可以使学生巩固、拓展和深化所学的基础理论、专业理论和知识。1.3通过课程设计可以使学生初步掌握用运筹学方法解决实际问题的过程和技巧，树立理论联系实际的工作作风。1.4通过课程设计可以
2、使学生初步建立正确的设计思想和方法，进一步提高运算、计算机应用技能和综合分析、解决问题的能力。2、设计进度安排本课程设计时间分为两周：第一周（2011年6月27日-2011年7月1日）：建模阶段。此阶段各小组根据给出的题目完成模型的建立。主要环节包括 ：1.16月27日上午：发指导书；按组布置设计题目；说明进度安排。1.26月27日下午至29日：各小组审题，查阅资料，进行建模前的必要准备(包括求解程序的编写与查找)。1.3 6月30日至7月1日：各个小组进行建模，并根据题目及设计要求拟定设计提纲，指导教师审阅；同时阅读，理解求解程序，为上机求解做好准备。第二周（2011年7月4日-7月8日）：
3、上机求解，结果分析及答辩。主要环节包括1.1 7月4日至7月5日：上机调试程序1.2 7月6日：完成计算机求解与结果分析。1.3 7月7日：撰写设计报告。 1.4 7月8日：设计答辩及成绩评定。（答辩时间：7月8日上午8：3012：00，下午2：005：00）组 别：第八组设计人员：王盛之、王有琳、王泽旭 设计时间：2007年6月27日-2007年7月8日1.设计进度：第一周（2011年6月27日-2011年7月1日）：建模阶段。此阶段各小组根据给出的题目完成模型的建立。主要环节包括：1.1 6月27日上午：发指导书；按组布置设计题目；说明进度安排。1.2 6月27日下午至29日：各小组审题，
4、查阅资料，进行建模前的必要准备(包括求解程序的编写与查找)。1.3 6月30日至7月1日：各个小组进行建模，并根据题目及设计要求拟定设计提纲，指导教师审阅；同时阅读，理解求解程序，为上机求解做好准备。第二周（2011年7月4日-7月8日）：上机求解，结果分析及答辩。主要环节包括：1.1 7月4日至7月5日：上机调试程序1.2 7月6日：完成计算机求解与结果分析。1.3 7月7日：撰写设计报告。 1.4 7月8日：设计答辩及成绩评定。2.设计题目 (第一题)某工厂生产三种原料，用B1、B2、B3来表示.该工厂用这三种原料采用不同的配料方法还可以生产出三种不同的成品,其数量分别为E1、E2、E3.
5、这三种不同的配料方法用A1.A2.A3表示.假定在每一种配料方法A1中,要求B1原料所占比例不低于1/3,要求B2原料所占比例不高于1/4,对B3原料没有要求;在第二种配料方法A2中,要求B1原料所占比例不低于1/4,要求B2原料所占比例不高于1/3,对B3原料仍没有要求；不符合前二种配料方法要求的即为第三种配料方法。A1配料方法所产大品的价格为8（千克/吨），A2配料方法所产大品的价格为5（千克/吨），A3配料方法所产大品的价格为4（千克/吨）。原料B1的价格为7（千克/吨），原料B2的价格为3（千克/吨），原料B3的价格为4（千克/吨）。该工厂每月可以生产这三种原料的数量分别为B1=200
6、吨，B2=200吨，B3=120吨。上述情况如表1所示。用所生产的B1，B2，B3三种原料再造成品，所获得的总产值一定比只生产B1、B2、B3三种原料所获得的总产值高，那么应该采用哪种配料方法，在所采用的各配料方法中，又应当使用多少原料，才可以使获得的总产值最大？并按要求分别完成下列分析：（1）A2成品的价格在何范围内变化时最优生产方案不变？（2）B2原料的价格在何范围内变化时最优生产方案不变？（3）B2原料的数量在何范围内变化时最优基不变？（4）当成品数量E2为100时的最优生产方案。表1配料方法原料B1 B2 B3成 品数 量成品价格（千元/吨）A11/3 1/4 不定E18A21/4 1
7、/3 不定E25A3不满足以上条件的配料E34原料数量（吨）200 200 120原料价格（千元/吨）7 3 43.建模过程：设Xij表示采用Ai种配料方针时，使用Bj种原料的数量；Ei表示采用第Ai种配料方法时得到的数量；i表示尚剩余Bj种原料的数量。即 X11+X12+X13 =E1 三种原料用第一种配料方法制成第一种成品的数量， X21+X22+X23 =E2 三种原料用第二种配料方法制成第二种成品的数量， X31+X32+X33 =E3 三种原料用第三种配料方法制成第三种成品的数量。而此时的总产值为Z=8E1+5E2+4E3+71+32+43=X11+5X12+4X13-2X21+2X
8、22+X23-3X31+X32+2480此问题的LP模型如下MaxZ=X11+5X12+4X13-2X21+2X22+X23-3X31+X32+2480 X11+X21+X31=200 B1原料的总库存量为200 X12+X22+X32=200 B2原料的总库存量为200 X13+X23+X33=120 B3原料的总库存量为120-2/3X11+1/3X12+1/3X13=0 原料X11的数量要多于A1方法中总原料的1/3-1/4X11+3/4X12-1/4X13=0 原料X12的数量要少于A1方法中总原料的1/4-3/4X21+1/4X22+1/4X23=0 原料X21的数量要多于A2方法中
9、总原料的1/4-1/3X21+2/3X22-1/3X23=0 (i=1,2,3,. j=1,2,3,)为了使我们所编的程序在上机时能够计算出结果，根据程序要求，加上松弛变量， 我们需输入如下的数据信息：MaxZ=X11+5X12+4X13-2X21+2X22+X23-3X31+X32+2480 X11+X21+X31+X1=200 X12+X22+X32+X2=200 X13+X23+X33+X3=120 -2/3X11+1/3X12+1/3X13+X4=0 -1/4X11+3/4X12-1/4X13+X5=0 -3/4X21+1/4X22+1/4X23+X6=0 -1/3X21+2/3X22
10、-1/3X23+X7=0 Xij=0 (i=1,2,3,. j=1,2,3,) Xi=0 (i=1,2,3,4,5,6,7)4.程序功能介绍4.1.程序及语言介绍Java是一种简单的，跨平台的，面向对象的，分布式的，解释的，健壮的安全的，结构的中立的，可移植的，性能很优异的多线程的，动态的语言。能运行于不同的平台，对程序提供了安全管理器，防止程序的非法访问。类的封装性、继承性等有关对象的特性，使程序代码只需一次编译，然后通过上述特性反复利用。Java提供了众多的一般对象的类，通过继承即可使用父类的方法。LINGO是一种专门用于求解数学规划问题的软件包。LINGO可以用于求解非线性规划，也可以用
11、于一些线性和非线性方程组的求解等，功能十分强大，是求解优化模型的最佳选择。其特色在于内置建模语言，提供十几个内部函数，可以允许决策变量是整数，方便灵活，而且执行速度非常快。Lingo 是使建立和求解线性、非线性和整数最佳化模型更快更简单更有效率的综合工具, 提供强大的语言和快速的求解引擎来阐述和求解最佳化模型开始4.2.程序流程图变量个数9、约束条件个数7、方程组系数矩阵A、约束条件常数项b值、选择目标函数类型、目标函数系数化成标准形式：调整目标函数为max，加入松弛变量，将数据输入相应的数组中，形成初始单纯形表基变量是否有非0人工变量变量是否所有检验数j0 是 否是否有非基变量检验数为0 否
12、 是否存在检验数所在列中有pj0 是无可行解是 否 否 否唯一最优解多个最优解 无界解 是输出最优表和最优解、最优值进行比较选择，确定主元和换入、换出变量 换基迭代结束4.3.数据录入输入方程组的系数矩阵A（7行16列）： 1 0 0 1 0 0 1 0 0 1 0 0 0 0 0 0 0 1 0 0 1 0 0 1 0 0 1 0 0 0 0 0 0 0 1 0 0 1 0 0 1 0 0 1 0 0 0 0 -2 1 1 0 0 0 0 0 0 0 0 0 1 0 0 0 -1 3 -1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 -3 1 1 0 0 0 0 0 0
13、0 0 1 0 0 0 0 -1 2 -1 0 0 0 0 0 0 0 0 0 1 输入初始基变量的数字代码num矩阵： 10 11 12 13 14 15 16输入方程组右边的值矩阵B： 200 200 120 0 0 0 0输入目标函数各个变量的系数所构成的系数阵C： 1 5 4 -2 2 1 -3 1 0 0 0 0 0 0 0 04.4.程序运行 用JAVA语言编程工具运行所编程序，在程序的各种指示下求得最优解、最优值，用LINGO对程序进行灵敏度分析分析。5. 结果分析5. 1.计算机运行结果 通过运行JAVA语言程序，经过5次迭代，在控制台得到了最优解。显示结果如下：X1=200.000 X2=106.666 X3
