用matlab化工具箱解线性规划.doc

佛山科学技术学院 上 机 报 告 课程名称 数学应用软件 上机项目 用matlab优化工具箱解线性规划 专业班级 一. 上机目的 1．了解二次线性规划和非线性规划的基本理论知识。 2．学习用matlab求解线性规划和非线性规划的解法问题。 二. 上机内容 1、用quadprog求解二次规划问题min f(x)： 2、求解优化问题： min S.T. 注：取初值为（10,10,10）。 3、求表面积为常数150 m2的体积最大的长方体体积及各边长。 注：取初值为（4,5,6）。 三. 上机方法与步骤 第1题：要用quadprog求解二次规划问题min f(x)，运用matlab命令[x,f]=quadprog(H,c,A,b,Aeq,beq,vlb,vub)编程求解。 H=[1 -1;-1 2]; c=[-2;-6]; A=[1 1;-1 2;2 1]; b=[2;2;3]; Aeq=[]; beq=[]; vlb=[0;0]; vub=[]; [x,f]=quadprog(H,c,A,b,Aeq,beq,vlb,vub) 第2题：求解优化问题minf(x),对于非线性规划问题，运用matlab命令[x,fval,exitflag,output]=fmincon(fun2,x0,A,b,Aeq,beq,vlb,vub,mycon2)求解。 function f=fun2(x); f=-x(1)*x(2)*x(3); function [g,ceq]=mycon2(x) g=[];1、 ceq=[]; x0=[10;10;10]; A=[1 2 2;-1 -2 -2]; b=[72;0]; Aeq=[]; beq=[]; vlb=[]; vub=[]; [x,fval,exitflag,output]=fmincon(fun2,x0,A,b,Aeq,beq,vlb,vub,mycon2) 第3题：求解优化问题maxf(x),对于非线性规划问题，运用matlab命令[x,fval]=fmincon(fun3,x0,A,b,Aeq,beq,vlb,vub,mycon3)求解。 function f=fun3(x); f=-x(1)*x(2)*x(3); function [g,ceq]=mycon3(x) ceq=[2*x(1)*x(2)+2*x(1)*x(3)+2*x(2)*x(3)-150]; g=[]; x0=[4;5;6]; A=[]; b=[]; Aeq=[]; beq=[]; vlb=[0;0;0]; vub=[]; [x,fval]=fmincon(fun3,x0,A,b,Aeq,beq,vlb,vub,mycon3) fval=-fval 上机结果 第1题： 第2题： 第3题： loan management (post) review: 1, approval criteria have been implemented; 2, the contract is correct; 3, in line with the contract, payment; 4, no other major cases against loan security;、Reviewed by: date of record form 9: qilu bank borrowing by single years of maximum amount a maximum no-loan borrowers: borrowers: ... Borrower name number, line of credit expiry date year month day to year month day maximum number of years in the loan contract are the top line opened, the borrowed and requested your bank opened credit lines above, is hereby requested. Borrower signature: date agreed to open credit lines. Customer