运筹课程设计3

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.36月30日至7月1日：各个小组进行建模，并根据题目及设计要求拟定设计提纲，指导教师审阅；同时阅读，理解求解程序，为上机求解做好准备。第二周（2011年7月4日-7月8日）：上机求解，结果分析及答辩。主要环节包括1.17月4日至7月5日：上机

3、调试程序1.27月6日：完成计算机求解与结果分析。1.37月7日：撰写设计报告。 1.47月8日：设计答辩及成绩评定。（答辩时间：7月8日上午8：3012：00，下午2：005：00）组别：第九组设计人员：吴维 熊众威 于菲菲 设计时间：2011年6月27日-2011年7月8日1、 设计进度：本课程设计时间分为两周：第一周（2011年6月27日-2011年6月29日）：建模阶段。此阶段各小组根据给出的题目完成模型的建立。主要环节包括：（1) 6月27日上午：发指导书；按组布置设计题目；说明进度安排。（2) 6月27日下午至29日：各小组审题，查阅资料，进行建模前的必要准备(包括求解程序的编写与

4、查找)。（3) 6月30日至7月1日：各个小组进行建模，并根据题目及设计要求拟定设计提纲，指导教师审阅；同时阅读，理解求解程序，为上机求解做好准备。第二周（2011年7月4日-7月6日）：上机求解，结果分析及答辩。主要环节包括：（1) 7月4日至5日：上机调试程序（2) 7月6日：完成计算机求解与结果分析。（3) 7月7日：撰写设计报告。（4) 7月8日：设计答辩及成绩评定。2、设计题目某投资公司有100万元资金用于投资，投资方案有六种，现要做一个5年期的投资计划，具体可选择的投资方案如下：方案A：5年内每年年初均可投资，且金额不限，投资期限一年，年投资回报率7%；方案B：5年内每年年初均可投

5、资，且金额不限，投资期限两年，年投资回报率10%（不计复利）；方案C：5年内每年年初均可投资，且金额不限，投资期限三年，年投资回报率12%（不计复利）；方案D：只在第一年初有一次投资机会，最大投资金额为50万元，投资期限四年，年投资回报率20%（不计复利）；方案E：在第二年和第四年初有一次投资机会，最大投资额为30万元，投资期限一年，年投资回报率30%；方案F：在第四年年初有一次投资机会，金额不限，投资期限两年，年投资回报率25%。假设当年的投资金额及其收益均可用于下一年投资，问公司应如何投资才能使第五年末收回的资金最多？并按要求分别完成下列分析：（1） 方案C的年投资回报率在何范围内变化时最

6、优投资方案不变？（2） 方案E的最大资金金额在何范围内变化时最优投资方案不变？（3） 最初投资额为200万元时的最优投资方案。3、建模过程定义变量：第一年：方案A投资额为x11，方案B投资额为x12，方案C投资额为x13，方案D投资额为x14；第二年：方案A投资额为x21，方案B投资额为x22，方案C投资额为x23，方案E投资额为x25；第三年：方案A投资额为x31，方案B投资额为x32，方案C投资额为x33；第四年：方案A投资额为x41，方案B投资额为x42，方案E投资额为x45；第五年：方案A投资额为x51，方案F投资额为x46；约束条件：第一年，方案A、B、C、D均可以投资，投资总额为最

7、大100万元，方案A的投资期限为一年，方案B为二年，方案C为三年，方案D为四年，且只能在第一年年初投资。所以有x11+x12+x13+x14=100 方案D 的最大投资额不能超过50万元。所以有x1450第二年，方案A、B、C、E均可以投资，方案E的投资期限为一年，投资总额为方案A在第一年的收益，所以有x21+x22+x23+x25-1.07x11=0方案E的最大投资额不的超过30万元，所以有x2530第三年，方案A、B、C均可以投资，且投资总额为方案A在第二年的收益，方案B在前两年的收益和方案E在第二年的收益的总和，所以有x31+x32+x33-1.07x21-1.2x12-1.3x25=0

8、第四年，方案A、B、E、F均可以投资，方案F的投资期限为两年，且投资总额为方案A在第三年的收益，方案B在第二和第三年的收益，方案C在前三年的收益的总和，所以有x41+x42+x45+x46-1.07x31-1.2x22-1.36x13=0方案E的最大投资额不得超过30万元，所以有x4530第五年，只有方案A可以投资，且投资总额为方案A在第四年的收益，方案B在第三和第四年的收益，方案D在前四年的收益，方案E在第三和第四年的收益，方案C在第二、三、四年得收益的总和，所以有x51-1.07x41-1.2x32-1.8x14-1.3x45-1.36x23=0xi,j0 i=1,2,3,4,5 j=1,

9、2,3,4,5,6maxz=1.07x51+1.2x42+1.36x33+1.5x46x11+x12+x13+x14=100 x1450x21+x22+x23+x25-1.07x11=0x2530x31+x32+x33-1.07x21-1.2x12-1.3x25=0x41+x42+x45+x46-1.07x31-1.2x22-1.36x13=0x4530x51-1.07x41-1.2x32-1.8x14-1.3x45-1.36x23=0xi,j0 i=1,2,3,4,5 j=1,2,3,4,5,64、程序功能介绍4.1、求解程序功能简介4.1.1、概述：该软件由本小组（熊众威、吴维、

10、于菲菲）专门为运筹学课程设计开发！该软件是所有较简单线性规划问题设计的，运用通用的单纯形法的解题思路，经过反复迭代而得最终解。4.1.2、功能简介：开发所用语言为Java，程序包括两个类，即：LP类和LP_Test类。其中LP_Test类中有main()方法为程序入口，负责基础数据的接收，如： * param minmax 表示：求函数的最大值或最小值； * param m 表示：约束条件的个数； * param n表示：变量个数； * param m1表示：<=的约束条件个数； * param m2表示：=的约束条件个数； * param m3表示：>=的约束条件个数； * pa

11、ram a表示：约束条件的系数矩阵； * param x表示：目标函数的价值系数；接收方式为从键盘接收。在接收到这些数据后，在LP_Test类中实例化LP类，再通过得到类的对象lp来调用LP类的的入口方法solve()方法，由此进入LP类，进行求解。程序在接收到了以上参数后，首先会反馈完整的系数矩阵到控制台上；然后进行反复的迭代，并将迭代过程全部显示到控制台上，最终的一个迭代表为最终表，表的最左边标明了所有的基变量，最上面一行为Z和检验数值，求解完成，如下图（运行截图1）：运行截图1最后程序给出最优解和最优值，程序结束。4.1.3、功能总结：这是一个经典的用单纯形法解线性规划问题的程序，运行程

12、序时只需按照程序的友好提示就可以很容易地得到模型结果，简单易用！开始4.2、程序功能介绍变量个数16、约束条件个数8、方程组系数矩阵A、约束条件常数项b值、选择目标函数类型、目标函数系数化成标准形式：加入松弛变量根据人工变量构造辅助LP问题max w，构造单纯形表进行迭代是否所有检验数j0 是是否存在检验数所在列中有pj>0无解 否 否 是进行比较选择，确定主元和换入、换出变量 是换基迭代 是人工变量所在行原始变量系数不全为零人工变量所在行原始变量系数全为零输出最优表和最优解、最优值W是否=0人工变量是否为基变量 否 是 是以非零系数其中之一为主元进行换基迭代，把人工变量变为非基变量删去

13、人工变量诸列，用Z代替w，用单纯形法求解删去相应行是 否是是否存在检验数所在列中有pj>0是否所有检验数j0 否进行比较选择，确定主元和换入、换出变量 否换基迭代输出最优表和最优解、最优值5、结果分析5.1计算机求解前的手工数据准备将模型代码化：目标函数类型：1。约束条件个数：8。变量个数：16。<=的约束条件个数：3。=的约束条件个数：5。>=的约束条件个数：0。系数矩阵：0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,500,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,300,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,30

14、1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,100-1.07,0,0,0,1,1,1,1,0,0,0,0,0,0,0,0,00,-1.2,0,0,-1.07,0,0,-1.3,1,1,1,0,0,0,0,0,00,0,-1.36,0,0,-1.2,0,0,-1.07,0,0,1,1,1,1,0,00,0,0,-1.8,0,0,-1.36,0,0,-1.2,0,-1.07,0,-1.3,0,1,0目标函数系数：0,0,0,0,0,0,0,0,0,0,1.36,0,1.2,0,1.5,1.075.2计算机求解结果截图如下（运行截图2到运行截图9）：第一步：进入程序运行截图2第二

15、步：输入预定义数据运行截图3第三步：输入系数矩阵运行截图4第四步：输入目标函数系数运行截图5第五步：反馈系数矩阵运行截图6第六步：进入迭代过程，第一次迭代运行截图7中间迭代省略，最后一次迭代运行截图8第七步：得到最优解和最优值，求解完成，程序结束运行截图9最优解如下：x11=28.04，x13=71.96，x25=30,x31=39，x46=139.6，其他变量取值为零，即在第一年方案A投资28.04万元，方案C投资71.96万元，第二年方案E投资30万元，第三年方案A投资39万元，第四年方案F投资139.6万元，最优值是209.3987，即第五年年末可收回的最大资金金额是209.3987万元

16、。5.3结果分析思路5.3.1问题分析通过对问题的理解和分析，根据题意可得到一个最大投资收益的模型，以这个模型为基础可以快速的求解各个年度的投资额以及投资收益，使在第五年年末可以得到最大的收益；然后通过做灵敏度分析，判断方案C的投资回报率何方案E的最大资金金额在何范围内变化时，使得最优投资方案不变，以及在投入资金额为200万元时的最优投资方案。5.3.2、结果分析最优解如下：x11=28.04，x13=71.96，x25=30,x31=39，x46=139.6，其他变量取值为零，即在第一年方案A投资28.04万元，方案C投资71.96万元，第二年方案E投资30万元，第三年方案A投资39万元，第

17、四年方案F投资139.6万元，最优值是209.3987，即第五年年末可收回的最大资金金额是209.3987万元。5.3.3、灵敏度分析方案C的投资回报率的变化属于目标函数中x33的系数的变化，其变化范围通过lindo分析结果可知，要使最优投资方案保持不变，其投资回报率的变化范围是：(0,20.17%。方案E的最大资金金额的变化属于约束条件右端项b的变化，其变化范围通过lindo分析结果可知，要使最优基保持不变其最大资金金额的变化范围如下：0,107。最初投资金额是200万元时，应属于约束条件右端项b的变化，其变化范围通过lindo分析结果可知，要使最优基保持不变其最大资金金额的变化范围如下：0

18、,100。当最初投资金额是200万元时，其最优投资方案已经改变，结果截图如下（运行截图10）：运行截图10最优解是x11=28.04,x13=171.96,x25=30,x31=39,x46=275.6，即在第一年方案A投资28.04万元，方案C投资171.96万元，第二年方案E投资30万元，第四年方案F投资275.6万元，最优值是413.3987万元，即第五年年末可收回的最大资金金额是413.3987万元。6、综合题目的创新的分析内容我针对这个投资问题做了进一步的分析，设计了几个问题，如下：（1） 假设方案C的投资回报率是40%时，最优投资方案如何变化？方案C的投资回报率的变化属于目标函数中

19、x33的系数的变化，其变化范围通过lindo分析结果可知，当方案C的投资回报率是40%时，其最优方案已经改变，模型如下：maxz=1.07x51+1.2x42+2.2x33+1.5x46x11+x12+x13+x14=100x1450x21+x22+x23+x25-1.07x11=0x2530x31+x32+x33-1.07x21-1.2x12-1.3x25=0x41+x42+x45+x46-1.07x31-1.2x22-2.2x13=0x4530x51-1.07x41-1.2x32-1.8x14-1.3x45-2.2x23=0xi,j0 i=1,2,3,4,5 j=1,2,3,4,5,6结果

20、如下截图如下（运行截图11）：运行截图11最优解如下：x13=100,x46=220，其他变量取值为零，即在第一年方案C投资100万元，第四年方案F投资220万元，最优值是330，即最终在第五年年末可收回的最大资金金额是330万元。（2） 假设方案B的最大投资金额是30万元时，最优投资方案如何变化？假设方案B的最大投资金额是30万元属于灵敏度分析中增加一个约束条件，模型如下：maxz=1.07x51+1.2x42+1.36x33+1.5x46x11+x12+x13+x14=100x1450x21+x22+x23+x25-1.07x11=0x2530x31+x32+x33-1.07x21-1.2

21、x12-1.3x25=0x41+x42+x45+x46-1.07x31-1.2x22-1.36x13=0x4530x51-1.07x41-1.2x32-1.8x14-1.3x45-1.36x23=0x1230x2230x3230x4230xi,j0 i=1,2,3,4,5 j=1,2,3,4,5,6结果截图如下（运行截图12）：运行截图12最优解如下：x11=28.04，x13=71.96，x25=30,x31=39，x46=139.6，其他变量取值为零，即在第一年方案A投资28.04万元，方案C投资71.96万元，第二年方案E投资30万元，第三年方案A投资39万元，第四年方案F投资139.6

22、万元，最优值是209.3987，即第五年年末可收回的最大资金金额是209.3987万元，最优方案不变。（3） 假设方案D的投资期限为三年时，最优投资方案如何变化？假设方案D的投资期限为三年属于灵敏度分析中A的变化，模型如下：maxz=1.07x51+1.2x42+1.36x33+1.5x46x11+x12+x13+x14=100 x1450x21+x22+x23+x25-1.07x11=0x2530x31+x32+x33-1.07x21-1.2x12-1.3x25=0x41+x42+x45+x46-1.07x31-1.2x22-1.36x13-1.6x14=0x4530x51-1.0

23、7x41-1.2x32-1.3x45-1.36x23=0xi,j0 i=1,2,3,4,5 j=1,2,3,4,5,6结果如下：运行截图13最优解如下：x11=28.04,x13=21.96,x14=50,x25=30,x31=39,x46=151.6，其他变量取值为零，即在第一年方案A投资28.04万元，第二年方案C投资21.96万元，方案D投资50万元，第三年方案E投资30万元，第四年方案F投资151.6万元。最优值是227.3987，即第五年年末可收回的资金为227.3987万元。（4） 增加一个方案G，可以在第一年年初投资，投资期限是两年，投资回报率是50%，最优投资方案如何变化？增加

24、一个方案G属于灵敏度分析中增加一个新变量，模型如下：maxz=1.07x51+1.2x42+1.36x33+1.5x46x11+x12+x13+x14+x17=100x1450x21+x22+x23+x25-1.07x11=0 x2530x31+x32+x33-1.07x21-1.2x12-1.3x25-2x17=0x41+x42+x45+x46-1.07x31-1.2x22-1.36x13=0 x4530x51-1.07x41-1.2x32-1.8x14-1.3x45-1.36x23=0xi,j0 i=1,2,3,4,5 j=1,2,3,4,5,6结果截图如下（运行截图14）：运行截图14最优解如下：x11=100,x31=200,x46=214,其他变量取值为零，即在第一年方案A投资100万元，第三年方案A投资200万元，第四年方案F投资214万元。最优值是321，即第五年年末可收回