运筹学课程设计

工业大学课程设计报告课程设计名称 运筹课程设计专 业 班 级 学生姓名 指导教师 2011年7月8日课程设计任务书课程设计题目第二题起止日期：2011、6、27〜2011、7、8|设计地点：教室、电子商务中心~设计任务及日程安排：1、 设计任务1.1通过课程设计可以增强学生解决实际运筹学问题的能力1.2通过课程设计可以使学生巩固、拓展和深化所学的基础理论、专业理论和知识。1.3通过课程设计可以使学生初步掌握用运筹学方法解决实际问题的过程和技巧，树立理论联系实际的工作作风。1.4通过课程设计可以使学生初步建立正确的设计思想和方法，进一步提高运算、计算机应用技能和综合分析、解决问题的能力。2、 设计进度安排本课程设计时间分为两周：第一周（2011年6月27日----2011年7月1日）：建模阶段。此阶段各小组根据给出的题目完成模型的建立。主要环节包括：6月27日上午：发指导书；按组布置设计题目；说明进度安排。6月27日下午至29日：各小组审题，查阅资料，进行建模前的必要准备包括求解程序的编写与查找）。6月30日至7月1日：各个小组进行建模，并根据题目及设计要求拟定设计提纲，指导教师审阅；同时阅读，理解求解程序，为上机求解做好准备。第二周（2011年7月4日---7月8日）：上机求解，结果分析及答辩。主要环节包括7月4日至7月5日：上机调试程序7月6日：完成计算机求解与结果分析。7月7日：撰写设计报告。7月8日：设计答辩及成绩评定。（答辩时间：7月8日上午8：30〜12：00，下午2：00〜5：00）组别：第九组设计人员：设计时间：2011年6月27日---2011年7月8日1、 设计进度：本课程设计时间分为两周：第一周（2011年6月27日——2011年6月29日）：建模阶段。此阶段各小组根据给出的题目完成模型的建立。主要环节包括：（1） 6月27日上午：发指导书；按组布置设计题目；说明进度安排。（2） 6月27日下午至29日：各小组审题，查阅资料，进行建模前的必要准备包括求解程序的编写与查找）。（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年内每年年初均可投资，且金额不限，投资期限两年，年投资回报率10%（不计复利）；方案C：5年内每年年初均可投资，且金额不限，投资期限三年，年投资回报率12%（不计复利）；方案D：只在第一年初有一次投资机会，最大投资金额为50万元，投资期限四年，年投资回报率20%（不计复利）；方案E：在第二年和第四年初有一次投资机会，最大投资额为30万元，投资期限一年，年投资回报率30%；方案F：在第四年年初有一次投资机会，金额不限，投资期限两年，年投资回报率25%。假设当年的投资金额及其收益均可用于下一年投资，问公司应如何投资才能使第五年末收回的资金最多？并按要求分别完成下列分析：（1） 方案C的年投资回报率在何范围内变化时最优投资方案不变？（2） 方案E的最大资金金额在何范围内变化时最优投资方案不变？（3） 最初投资额为200万元时的最优投资方案。3、建模过程3.1定义变量：第一年：方案A投资额为X］］，方案B投资额为x12，方案C投资额为x13，方案D投资额为x］4；第二年：方案A投资额为x2］，方案B投资额为x22，方案C投资额为x23，方案E投资额为x25；第三年：方案A投资额为x3］，方案B投资额为x32，方案C投资额为x33；第四年：方案A投资额为x4］，方案B投资额为x42，方案E投资额为x45；第五年：方案A投资额为x5］，方案F投资额为x46；3.2约束条件：第一年，方案A、B、C、D均可以投资，投资总额为最大100万元，方案A的投资期限为一年，方案B为二年，方案C为三年，方案D为四年，且只能在第一年年初投资。所以有x11+x12+x13+x14=100方案D的最大投资额不能超过50万元，所以有x14W50；第二年，方案A、B、C、E均可以投资，方案E的投资期限为一年，投资总额为方案A在第一年的收益，所以有x21+x22+x23+x25-1.07x11=0方案E的最大投资额不的超过30万元，所以有x25W30；第三年，方案A、B、C均可以投资，且投资总额为方案A在第二年的收益，方案B在前两年的收益和方案E在第二年的收益的总和，所以有x,+xo+xo-1.07x,-1.2xo-1.3xc=0；31 32 33 21 12 25

第四年，方案A、B、E、F均可以投资，方案F的投资期限为两年，且投资总额为方案A在第三年的收益，方案B在第二和第三年的收益，方案C在前三年的收益的总和，所以有x+x+x+x-1.07x-1.2x-1.36x=041 42 45 46 31 22 13方案E的最大投资额不得超过30万元，所以有x45W30；第五年，只有方案A可以投资，且投资总额为方案A在第四年的收益，方案B在第三和第四年的收益，方案C在第二、三、四年的收益，方案D在前四年的收益，方案E在第四年的收益的总和，所以有x-1.07x-1.2x-1.8x-1.3xc-1.36x=0；51 41 22 14 45 233.3目标函数：maxz=1.07x51+1.2x42+1.36x33+1.5x46第五年末收回的资金额为所有投资可在第五年收回的投资之和，即方案A在第五年的收益，方案B在第四和第五年的收益，方案C在第三、四、五年的收益，方案F在在第四和第五年的收益的总和，求其最大值即为最有投资方案。3.4模型：maxz=1.07x51+1.2x42+1.36x33+1.5x46x11+x12+x13+x14=100x14<50x25<30x21x21+x22+x23+x25-1.07x11=0x+x+x-1.07x-1.2x-1.3x=0TOC\o"1-5"\h\z31 32 33 21 12 25x+x+x+x-1.07x-1.2x-1.36x=041 42 45 46 31 22 13x45<30x-1.07x-1.2x-1.8x-1.3x-1.36x=051 41 32 14 45 23x尸0i={1,2,3,4,5}j={1,2,3,4,5,6}4、程序功能介绍4.1求解程序功能简介4.1.1概述:该软件由本小组（熊众威、吴维、于菲菲）专门为运筹学课程设计开发！该软件是所有较简单线性规划问题设计的，运用通用的单纯形法的解题思路，经过反复迭代而得最终解。

4.1.2功能简介：开发所用语言为Java,程序包括两个类，即：LP类和LP_Test类。其中LP_Test类中有main()方法为程序入口，负责基础数据的接收。如：*@paramminmax表示：求函数的最大值或最小值;@paramm@paramm@paramn@paramml@paramm2@paramm3@parama@paramx表示：变量个数；表示：<=的约束条件个数；表示：=的约束条件个数；表示：>=的约束条件个数；表示：约束条件的系数矩阵;表示：目标函数的价值系数;接收方式为从键盘接收。在接收到这些数据后，在LP_Test类中实例化LP类，再通过得到类的对象lp来调用LP类的的入口方法solve()方法，由此进入LP类，进行求解。程序在接收到了以上参数后，首先会反馈完整的系数矩阵到控制台上；然后进行反复的迭代，并将迭代过程全部显示到控制台上，最终的一个迭代表为最终表，表的最左边标明了所有的基变量，最上面一行为Z和检验数值，求解完成，反复迭代过程如下：Z=209.39873831775702 0.1799579439252339 2.04K(17)50.0 0.0 0.0 0.0 -0.7222222222(1139.599158878504680.119971962616822322B.0373831775700920.9345794392523364 0.05>(1X0.0反复迭代过程如下：Z=209.39873831775702 0.1799579439252339 2.04K(17)50.0 0.0 0.0 0.0 -0.7222222222(1139.599158878504680.119971962616822322B.0373831775700920.9345794392523364 0.05>(1X0.0X(19)29.9999999999999930.00.00.00.0-1.070.01.0 1.00.00.0 0.0 1.00.01.0-0.72222222222222230.934579冬3925233640.93457943925233640.179023883696780.00.00.00.00.0 0.00.00.0 0.0 0.00.01.35999999999999990.075999999999999960.017777777777777684 1.27102803738317770.0C."7555555555555553-.-0.1.00.00.00.00.0 0.00.00.0 0.0 0.00.0 0.00.11400000000000009 0.026666666666666426 1.9065420560747661 0.10654205607476638.223 0.0 0.0 -0.7555555555555558 0.0 0.5555555555555556 -0.66666666666660.0 -0.0 -0.9345794392523364 -0.9345794392523364 -0.9345794392523364 -0.93457程序运行截图最后程序给出最优解和最优值，程序结束。4.1.3功能总结：这是一个经典的用单纯形法解线性规划问题的程序，运行程序时只需按照程序的友好提示就可以很容易地得到模型结果，简单易用！

4.2程序流程图开始变量个数16、约束条件个数8、方程组系数矩阵A、约束条件常

数项b值、选择目标函数类型、目标函数系数开始化成标准形式：加入松弛变量ojN0否是是否所有检验数是否存在检验数所在列中有pj>0无解根据人工变量构造辅助LP问题maxwojN0否是是否所有检验数是否存在检验数所在列中有pj>0无解进行比较选择，确定主元和换入、换出变量输出最优表和最优解、最优值换基迭代是是全为零为零变量是否*删去人工变量诸列，用Z代替<是删去相应行W是否=0人工变量所在行人工变量所在行输出最优表和最优解、最优值换基迭代是是全为零为零变量是否*删去人工变量诸列，用Z代替<是删去相应行W是否=0人工变量所在行人工变量所在行原始变量系数不人工变量是否为基原始变量系数全以非零系数其中之一为

主元进行换基迭代，把人

工变量变为非基变量w,用单纯形法求解进行比较选择，确定主元和换入、换出变量输出最优表和最优解、最优值换基迭代 k是否存在检验数所.在列中有pj>0是否所有检验数否<；否<；5、结果分析5.1计算机求解前的手工数据准备将原问题添加松弛变量、人工变量化成标准形式。目标函数：maxz=1.07x5i+1.2x42+1.36x33+1.5x46约束条件：x11+x12+x13+x14+y4=10051=50x21+x22+x23+x25-1.07x11+y5=0x25+y2=30x31+x32+x33-1.07x21-1.2x12-1.3x25+y6=0xA1+xA+xA+xA-1.07x^-1.2x^-1.36x1+y^=041 42 45 46 31 22 13 7x45+y3=30x51-1.07x41-1.2x22-1.8x14-1.3x45-1.36x23+y8=05.2计算机求解结果过程如下：第一步：录入数据，包括目标函数类型、约束条件个数、变量个数、W的约束条件个数、巳的约束条件个数、目标函数系数，约束条件系数。程序运行截图如下：请输入目标函数类型：2表示心;T表示min:备输入约束条件个规如输入变量个数（注意，只输入原始变量数，即不包括引入变量）:16请输沃<=的约束条件个数：蓄输入=的约束条件个数：请输笑>=的约束条件个巍据录入过程截图1

冉新.八元=的，最后制袈完>=叫）.：；：拆愉人完=的，最后渝叛完>=的》:再输人完=的,最后输入完>=的9：再输其完=的，最后:新入完X的》W•再输卷完=的，最后输久完>=的":•溥输入完=的，最后输隧完>=的》:勇-输人完=的，最后输入完>=的¥：再输入完=的，最后输;X完>=的〉:清输入目标函数系数（注:•数与数间用,广冉新.八元=的，最后制袈完>=叫）.：；：拆愉人完=的，最后渝叛完>=的》:再输人完=的,最后输入完>=的9：再输其完=的，最后:新入完X的》W•再输卷完=的，最后输久完>=的":•溥输入完=的，最后输隧完>=的》:勇-输人完=的，最后输入完>=的¥：再输入完=的，最后输;X完>=的〉:请新八承数矩阡第1行〔注：数与数同用,/间隔t升且先摘LA•完■<=的,o7o,ofi/07o/o7o/07o/o7o/07o/o7a（.5o请输入系数矩阵第w行乂注:数与数间用/间隔，并且先输密完==的'0,0,0,0,0..0,0,1,0,0,0,0,0,0,0,0,3：|请输支系数矩阵第3行（注:数与数间用/间隔，并且先输京完.的,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,30请输入系数矩阵第q行（注:数与数间用.”，,间隔，并且先输入•完。的;1,1,1,1,0,0,0,0,0.0,0,0,0,0,0,0,130请输入系数矩阵第5行（注：数与数间用>间隔t并且先输入•完一的,-1.07,0,0,0,1,1,1,1,0,0,0,0,0,0,0,0,0请输入系数矩阵第6行3注：数与数间用七,间隔，并且先输肉完一的,请输入系数矩阵第7行（注：数与数间用/间隔，并且先输入完。的,0/。『一.1.36,0,0—1.2r0,0,-■1.07,0r0,1,1.1r1r00请输入系数矩阵第8行〔注：数与数间用,，涧隔，并且先输外•完<=的;0,0,0,-1.3,0,-1.2,-1.36,0,3,0,0,-1.07,0,-1.3,0,1,0数据录入过程截图3通过录入的数据，程序反馈完整的系数矩阵如下，前8行表示按程序提示要求输入的增广矩阵，结合此题为3个“W”的约束条件系数、5个“=”的约束条件系数，共17列，前16列为约束条件系数，第17列为约束条件右端项。最后一行是目标函数系数矩阵。程序运行截图如下：辛苦了,您输入的系数矩阵为：0.00.00.01.00.00.00.00.00.00.00.00.00.00.00.00.050.00.00.00.00.00.00.00.01.00.00.00.00.00.00.00.00.030.0□.00.00.00.00.00.00.0□.00.0Q.00.00.00.01.00.00.030.01.01.01.01.00.00.00.00.00.00.00.00.00.00.00.00.0100.01.070.00.00.01.01.01.01.00.00.00.00.00.00.00.00.00.00.01.20.00.01.070.00.01.31.01.01.00.00.00.00.00.00.Q0.00.01.360.00.0-1.20.00.01.070.00.01.01.01.01.00.00.00.00.00.0-1,8Q.0-1.2-1,360.00.00.00.0-1.070.0-1.30,01.00.00.00.00.00.00.00.00.00.00.00.01.360.01.20.01.51.07反馈系数矩阵截图

第二步：二阶段法求解。第一阶段反复迭代过程如下:反复迭代过程如下:Z=-00-0.0-1.6320000000000001-0.0-0.0-1.4552000000000003-0.0-0.0-1.76800000000000021.361.361.36-0.0X(17)50.00.0 0.0 0.01.00.00.0 0.0 0.00.00.0-0.0 0.0 0.00.00.00.0X(18)30.00.0 0.0 0.00.00.00.0 0.0 1.00.00.0-0.0 0.0 0.00.00.00.0X(19)30.00.0 0.0 0.00.00.00.0 0.0 0.00.00.0-0,0 0.0 0.01.00.00.0X(20)100.01.0 1.0 1.01.00.00.0 0.0 0.00.00.0-0.0 0.0 0.00.00.00.0X(21)0.0-1.07 0.0 0.00.01.01.0 1.0 1.00.00.0-0.0 0.0 0.00.00.00.0X(ll)0.00.0 -1.2 0.00.0-1.070.0 0.0 -1.31.01.01.0 0.0 0.00.00.00.0X(23)0.00.0 0.0 -1.360.00.0-1.2 0.0 0.0-1.070.0-0,0 1.0 1.01.01.00.0X(24)0.00.0 0.0 0.0-1.80.00.0 -1.36 0.00.0-1.2-0.0 -1.07 0.0-1.30.01.0迭代运算过程截图1第二阶段反复迭代过程如下：反复迭代过程如下:Z=209.398738317757020.17995794392523392.040.114000000000000090.0266666666666664261.90654205607476610.10654205607476638X(17)50.0卜⑴X(9)X(3)fe(8>0.028.03738317757009239.0 1.371.9626168224299130.0 1.0599158878504680.00.00.00.00.0-0.75555555555555580.00.5555555555555556-0.66666666666660.93457943925233640.0-0.0-0.9345794392523364-0.9345794392523364-0.9345794392523364-0.934570.0-1.2-0.00.00.0-1.070.01.01.00.00.00.01.00.0-0.93457943925233641.0-0.72222222222222230.93457943925233640.93457943925233640.179023883696780.00.0-0.00.00.00.00.00.00.00.00.00.0X(15)139.*(冬>o.o|X(19)29.9999999999999930.1199719626168223235999999999999990.075999999999999960.0177777777777776841.27102803738317770.00.00.00.00.75555555555555580.0-0.55555555555555560.666666666666660.00.00.01.00.00.0Z=209.398738317757020.17995794392523392.040.114000000000000090.0266666666666664261.90654205607476610.10654205607476638X(17)50.0卜⑴X(9)X(3)fe(8>0.028.03738317757009239.0 1.371.9626168224299130.0 1.0599158878504680.00.00.00.00.0-0.75555555555555580.00.5555555555555556-0.66666666666660.93457943925233640.0-0.0-0.9345794392523364-0.9345794392523364-0.9345794392523364-0.934570.0-1.2-0.00.00.0-1.070.01.01.00.00.00.01.00.0-0.93457943925233641.0-0.72222222222222230.93457943925233640.93457943925233640.179023883696780.00.0-0.00.00.00.00.00.00.00.00.00.0X(15)139.*(冬>o.o|X(19)29.9999999999999930.1199719626168223235999999999999990.075999999999999960.0177777777777776841.27102803738317770.00.00.00.00.75555555555555580.0-0.55555555555555560.666666666666660.00.00.01.00.00.00.00.00.00.00.00.00.00.00.0迭代运算过程截图2第三步：求得最优解及最优值。程序运行截图如下：解题完成，最忧解与最忧值如下：xl=28.037383177570092x2=0x3=71.962616S2242991x4=0.0x5=0x6=0x7=Qx8=30.0x9=39.0xl0=0xll=0X13=Ox!4=0x25=139.59925e8785046Sxl€=0最优值为:209.39873331775702完成解题，.谢谢使用！程序运行结果截图由以上解题过程截图可知，最优解为：X]]=100,x22=107,x46=128.4,x51=128.4，即在第一年方案A投资100万元，第二年方案B投资107万元，第四年方案F投资128.4万元，第五年方案A投资128.4万元。最优值是329.988,即在第五年年末可收回最大资金金额是329.988万元。5.3结果分析思路5.3.1问题分析通过对问题的理解和分析，根据题意可得到一个最大投资收益的模型，以这个模型为基础可以快速的求解各个年度的投资额以及投资收益，使在第五年年末可以得到最大的收益；然后通过做灵敏度分析，判断方案C的投资回报率和方案E的最大资金金额在何范围内变化时，使得最优投资方案不变，以及在最初投入资金额为200万元时的最优投资方案。5.3.2结果分析最优解如下：x11=100,x22=107,x46=128.4,x51=128.4，其他变量取值为零，即在第一年方案A投资100万元，第二年方案B投资107万元，第四年方案F投资128

万元，第五年方案A投资128.4万元。最优值是329.988,即第五年年末可收回的最大资金额为329.988万元。5.3.3灵敏度分析通过lindo得到的灵敏度分析结果如下图：RANGESINUHICHTHEBASISISUNCHANGEDOBJCOEFFICIENTRANGESUARIABLECURRENTALLOWABLEALLOWABLECOEFINCREASEDECREASEX511.0700000.063333INFINITYX421.2000000.300000INFINITYX331.3600000.245000INFINITYX461.500000INFINITY0.083823X110.0000000.1140000.192555N120.0000000.114000INFINITYN130.0000000.1925550.114000N140.0000000.114000INFINITYX210.0000000.189192INFINITYX220.0000000.106542INFINITYX230.0000008.451342INFINITVX250.000000INFINITY0.179958X310.0000000.0950000.138429X320.0000000.321000INFINITYX410.0000000.355100INFINITY0.0000000.109000INFINITYRIGHTHANDSIDERANGESROMCURRENTALLOWABLEALLOWABLERHSINCREASEDECREASE2100.000000INFINITY71.962616350.000000INFINITY50.00000040.00000030.00000277.000000530.00000077.00000030.00000060.000000INFINITV39.00000070.000000INFINITY139.599152830.000000INFINITY30.00000090.000000INFINITY0.000000方案C的投资回报率的变化属于目标函数中x33的系数的变化，其变化范围通过lindo分析结果可知，要使最优投资方案保持不变，其投资回报率的变化范围是：（0,20.17%］。方案E的最大资金金额的变化属于约束条件右端项b的变化，其变化范围通过lindo分析结果可知，要使最优基保持不变其最大资金金额的变化范围如下：［0,107］。最初投资金额是200万元时，应属于约束条件右端项b的变化，其变化范围通过lindo分析结果可知，要使最优基保持不变其最大资金金额的变化范围如下：[28.04,100]。当最初投资金额是200万元时，其最优投资方案已经改变，结果如下:最优解是X]]=28.04,x13=171.96,x25=30,x31=39,x46=275.6,其他变量取值为零，即在第一年方案A投资200万元，方案C投资171.96万元，第二年方案E投资30万元，第三年方案A投资39万元，第四年方案F投资275.6万元。最优值是413.3987万元，即第五年年末可收回的最大资金金额是413.3987万元。1)413.3987VARIABLEVALUEREDUCEDCOSTX510.0000000.063333X420.0000000.300000X330.0000000.245000X46275.5991520.000000X1128.0373840.000000X120.0000000.114000X13171.9626160.000000X140.0000000.000000X210.0000000.189192X220.0000000.106542X230.0000000.365209X2530.0000000.000000X3139.0000000.000000X320.0000000.245000X410.0000000.287333X450.0000000.026667lindo运行结果截图6、综合题目创新的分析内容（1）假设方案A的投资回报率是10%时，最优投资方案如何变化？方案A的投资回报率的变化属于目标函数中x51的系数以及约束条件系数矩阵a的变化，变化后的模型如下：maxz=1.1x51+1.2x42+1.36x33+1.5x46x11+x12+x13+x14=100x14<50x21+x22+x23+x25-L1x11=0x25<30x”+x”+x”-1.1x〜-1.2x“-1.3x*=0TOC\o"1-5"\h\z31 32 33 21 12 25x八+x』c+x』u+x^-1.1x»-1.2xs-1.36x“=041 42 45 46 31 22 13x45<30x-1.1x-1.2xo-1.8x-1.3xc-1.36xo=051 41 32 14 45 23x尸0i={1,2,3,4,5}j={1,2,3,4,5,6}lindo运行求解结果如下：1)212.7136UARIABLEUALUEREDUCEDCOSTX510.0000000.000000X420.0000000.300000X330.0000000.290000X46141.8090970.000000X1127.2727280.000000X120.0000000.060000X1372.7272720.000000X140.0000000.060000X210.0000000.039545X220.0000000.054545X230.0000000.358545X2530.0000000.000000X3139.0000000.000000X320.0000000.330000X410.0000000.290000XUS0.0000000.070000最优解如下：x11=27.27,x13=72.73,x25=30,x31=39,x46=141.8，其他变量取值为零，即在第一年方案A投资27.27万元，方案C投资72.73万元，第二年方案E投资30万元，第三年方案A投资39万元，第四年方案F投资141.8万元。最优值是212.7136，即第五年年末可收回的最大资金金额为212.7136万元。（2）假设所有投入资金在第六年收回，公司应如何投资才能使第六年末收回的资金最多？投资收回期在第六年收回，属于灵敏度分析中增加一个约束条件的情况，其模型如下：maxz=1.07x61+1.2x52+1.36x43x11+x12+x13+x14=100x14<50

X21X21+X22+X23+X25-1.07x11=0X25<30x+x+xo-1.07x-1.2xo-1.3x=0X41X41+X42+X43+X45+X46-1.07X31-1.2X22-1.36X13=0X45<30x+x-1.07x-1.2x-1.8x-1.3xc-1.36x=051 52 41 32 14 45 23X61-1.07X51-1.2X42-1.36X33-1.5X46=0X尸0i={1,2,3,4,5,6}j={1,2,3,4,5,6}lindo运行求解结果如下：1)224.0567UARIABLEUALUEREDUCEDCOSTX61209.3987430..000000X520.0000000..012667X430.0000000..245000X1128.0373840..000000X120.0000000..121980X1371.9626160..000000X140.0000000..000000X210.0000000..202435X220.0000000..114000X230.0000000..390773X2530.0000000..000000X3139.0000000..000000X320.0000000..262150X330.0000000..262150X410.0000000..307447X420.0000000..321000X450.0000000..028533X46139.5991520..000000X510.0000000..067767最优解如下：x“=28.04,X|尸71.96,x“=30,x«=39,X4尸139.6,x°=209.4,其他变11 13 25 31 46 61量取值为零，即在第一年方案A投资28.04万元，方案C投资71.96万元，第二年方案E投资30万元，第三年方案A投资39万元，第四年方案F投资139.6万元，第六年方案A投资209.4万元。最优值是224.06,即第五年年末可收回的最大资金金额是224.06万元。（3）假设方案E的最大投资金额为20万元时，最优投资方案如何变化？假设方案E的投资期限为五年属于灵敏度分析中约束条件右端项b的变化，变化后的模型如下：maxz=1.07x5i+1.2x42+1.36x33+1.5x46X11+X12+X13+X14=100x14<50x21+x22+x23+x25-1.07x11=0乂25五。x+x+xo-1.07x-1.2xo-1.3x=031 32 33 21 12 25xAA+xA+xA+xA-1.07x^-1.2x^-1.36xA=041 42 45 46 31 22 13x45^20x51-1.07x41-1.2x32-1.3x45-L36x23=0xid^0i={1,2,3,4,5}j={1,2,3,4,5,6}lindo运行求解结果如下：1)207.5992VARIABLEVALUEREDUCEDCOSTX510.0000000.000000X420.0000000.300000X330.0000000.245000X46138.3994450.000000X1118.6915890.000000X120.0000000.114000X1381.3884110.000000X140.0000000.114000X210.0000000.189192X220.0000000.106542X230.0000000.451342X2520.0000000.000000X3126.0000000.000000X320.0000000.321000X410.0000000.355100X450.0000000.109000最优解如下：x11=18.7,x13=81.3,x25=20,x31=26,x46=138.4，其他变量取值为零，即在第一年方案A投资18.7万元，方案C投资81.3万元，第二年方案E投资20万元，第三年方案A投资26万元，第四年方案F投资138.4万元。最优值是207.6，即第五年年末可收回的最大资金金额为207.6万元。（4）增加一个方案G，可以在第二年年初投资，投资期限是一年，投资回报率是30%，最优投资方案如何变化？增加一个方案G属于灵敏度分析中增加一个新变量，变化后的模型如下:maxz=1.07x5i+1.2x42+1.36x33+1.5x46+1.3x57X11+X12+X13+X14+X17=100x14<50x4+x+x+x+x^-1.07x户0TOC\o"1-5"\h\z21 22 23 25 27 11X25<30x睥+Xs+x”+Xb-1.07Xf-1.2Xs-1・3x*-1.3x”=031 32 33 37 21 12 25 27xaa+xa^+xa+xa+xa-1.07xq-1.2x^-1.36xa-2.2xa-1.3xq=041 4