數(shù)學應用軟件作業(yè)5-用MATLAB求解非線性規(guī)劃問題

佛山科學技術學院上機報告課程名稱數(shù)學應用軟件上機工程用MATLAB求解非線性規(guī)劃問題專業(yè)班級姓名學號一.上機目的1．了解非線性規(guī)劃的根本理論知識。2．比照Matlab求解線性規(guī)劃，學習用Matlab求解非線性規(guī)劃的問題。二.上機內容1、用quadprog求解二次規(guī)劃問題minf(x)：2、求解優(yōu)化問題：minS.T.注：取初值為〔10,10,10〕。3、求外表積為常數(shù)150m2注：取初值為〔4,5,6〕。上機方法與步驟1、可用兩種方法解題：方法一：Matlab程序:H=[1-1;-12];c=[-2;-6];A=[11;-12;21];b=[2;2;3];Aeq=[];beq=[];vlb=[0;0];vub=[];[x,z]=quadprog(H,c,A,b,Aeq,beq,vlb,vub)方法二：Matlab程序如下：先建立fun.m文件，程序為：functionf=fun(x);f=1/2*x(1)^2+x(2)^2-x(1)*x(2)-2^x(1)-6*x(2);再建立chushi.m文件，程序為：x0=[1;1];A=[11;-12;21];b=[2;2;3];Aeq=[];beq=[];vlb=[0;0];vub=[];[x,fval]=fmincon('fun',x0,A,b,Aeq,beq,vlb,vub)2、Matlab程序:先建立fun1.m文件，程序為：functionf=fun1(x);f=-x(1)*x(2)*x(3);再建立chushi1.m文件，程序為：x0=[10;10;10];A=[122;-1-2-2];b=[72;0];Aeq=[];beq=[];vlb=[];vub=[];[x,fval]=fmincon('fun1',x0,A,b,Aeq,beq,vlb,vub)3、假設長方形的長、寬、高分別為、、，那么長方形的體積為f,那么有Matlab程序:先建立fun2.m文件，程序為：functionf=fun2(x);f=-x(1)*x(2)*x(3);然后建立mycon.m文件，程序為：function[g,ceq]=mycon(x)g=[];ceq=[2*x(1)*x(2)+2*x(1)*x(3)+2*x(2)*x(3)-150];再建立chushi2.m文件，程序為：x0=[4;5;6];A=[];b=[];Aeq=[];beq=[];vlb=[0;0;0];vub=[];[x,fval]=fmincon('fun2',x0,A,b,Aeq,beq,vlb,vub,'mycon')f=-fval上機結果結果：〔1〕方法一結果：x=0.66671.3333z=-8.2222〔2〕方法二結果：x=0.66671.3333fval=-8.47632、結果：x=24.000012.000012.0000fval=-3.4