数学建模实验4-线性规划模型求解.docVIP

湖南第一师范学院数学系实验报告 姓名： 学号： 专业： 数学与应用数学 班级： 课程名称： 线性规划与数学建模 实验名称： 线性规划模型的Matlab求解 实验类型： 基础实验 实验室名称： 数学建模实验室 实验地点： 实A302 实验时间： 2015年6月7日 指导教师： 成绩评定： 一、实验目的与要求： 掌握线性规划模型求解方法。 熟练掌握函数linprog等的用法。 会利用matlab优化工具箱求解简单的优化问题。 二、实验环境（实验器材、环境要求）： 计算机 Matlab软件 三、实验内容（实验原理、任务等）： 1、求解下列线性规划问题： 2、求解下列二次规划问题： 3、（连续投资问题）某部门有现金10万元，在今后五年内考虑给下面项目投资，已知： 项目A，从第一年到第四年每年年初需要投资，并于次年末回收本利115%; 项目B，第三年初需要投资，到第五年末回收本利125%，但最大投资额不超过4万元； 项目C，第二年初需要投资，到第五年末回收本利140%，但最大投资额不超过3万元； 项目D，五年内每年初可购买公债，于当年末归还，并加利息6%； 问：如何投资到第五年末时总收益最大？ 四、实验具体步骤： 1. min 6x1+3x2+4x3 f=[6 3 4] A=[]; b=[]; Aeq=[1 1 1]; beq=[120]; lb=[30 15 20]; [xopt,fopt]=linprog(f,A,b,Aeq,beq,lb) ans = 43 f = 6 3 4 Optimization terminated. xopt = 30.0000 70.0000 20.0000 fopt = 470.0000 2. H=[2 -2 3;-2 4 0;3 0 0]; C=[0 4 0]; A=[2 3 1]; b=[8]; Aeq=[1 -1 2]; beq=[8]; lb=zeros(3,1); [xopt,fopt]=quadprog(H,C,A,b,Aeq,beq,lb) Optimization terminated. xopt = 0 0 4 fopt = 0 3. 一 二 三 四 五 A X1A X2A X3A X4A B X3B C X2C D X1D X2D X3D X4D X5D Max z = 1.15*X4A+1.25*X3B+1.40*X2C+1.06*X5D S.t X1A+X1D=100000 X2A+X2C+X2D=1.06*X1D X3A+X3B+X3D=1.06*X2D+1.15*X1A X4A+X4D=1.06*X3D+1.15*X2A X5D=1.15*X3A+1.06*X4D X3B=40000 X2C=30000 X1A,X2A,…,X5D=0 编制M文件： max 1.15x4a+1.40x2c+1.25x3b+1.06x5d f=[0 0 0 -1.15 -1.25 -1.40 0 0 0 0 -1.06]; A=[0 0 0 0 1 0 0 0 0 0 0;0 0 0 0 0 1 0 0 0 0 0]; b=[40000 30000]; Aeq=[1 0 0 0 0 0 1 0 0 0 0;0 1 0 0 0 1 -1.06 1 0 0 0;-1.15 0 1 0 1 0 0 -1.06 1 0 0;0 -1.15 0 1 0 0 0 0 -1.06 1 0;0 0 -1.15 0 0 0 0 0 0 -1.06 1]; beq=[100000 0 0 0 0]; lb=zeros(11,1); [xopt,fopt]=linprog(f,A,b,Aeq,beq,lb) ans = 120 Optimization terminated. xopt = 1.0e+004 * 5.8912 1.3553 1.6367 2.7651 4.0000 3.0000 4.1088 0.0000 1.1381 0.0000 1.8822 fopt = -1.4375e+005 五、实验心得（质疑、建议）： Matlab不仅功能强大，其接近数学表达式的自然化语言，也使我们易于学习和掌握。本节课学习的是linprog等函数的用法。在学习的过程中，我发现在使用linprog命令时，系统默认它的参数至少为3个，但是当我