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

佛山科學(xué)技術(shù)學(xué)院上機(jī)報(bào)告課程名稱數(shù)學(xué)應(yīng)用軟件上機(jī)工程用MATLAB求解非線性規(guī)劃問題專業(yè)班級(jí)姓名學(xué)號(hào)一.上機(jī)目的1．了解非線性規(guī)劃的根本理論知識(shí)。2．比照Matlab求解線性規(guī)劃，學(xué)習(xí)用Matlab求解非線性規(guī)劃的問題。二.上機(jī)內(nèi)容1、用quadprog求解二次規(guī)劃問題minf(x)：2、求解優(yōu)化問題：minS.T.注：取初值為〔10,10,10〕。3、求外表積為常數(shù)150m2注：取初值為〔4,5,6〕。上機(jī)方法與步驟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文件，程序?yàn)椋篺unctionf=fun(x);f=1/2*x(1)^2+x(2)^2-x(1)*x(2)-2^x(1)-6*x(2);再建立chushi.m文件，程序?yàn)椋簒0=[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文件，程序?yàn)椋篺unctionf=fun1(x);f=-x(1)*x(2)*x(3);再建立chushi1.m文件，程序?yàn)椋簒0=[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、假設(shè)長(zhǎng)方形的長(zhǎng)、寬、高分別為、、，那么長(zhǎng)方形的體積為f,那么有Matlab程序:先建立fun2.m文件，程序?yàn)椋篺unctionf=fun2(x);f=-x(1)*x(2)*x(3);然后建立mycon.m文件，程序?yàn)椋篺unction[g,ceq]=mycon(x)g=[];ceq=[2*x(1)*x(2)+2*x(1)*x(3)+2*x(2)*x(3)-150];再建立chushi2.m文件，程序?yàn)椋簒0=[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上機(jī)結(jié)果結(jié)果：〔1〕方法一結(jié)果：x=0.66671.3333z=-8.2222〔2〕方法二結(jié)果：x=0.66671.3333fval=-8.47632、結(jié)果：x=24.000012.000012.0000fval=-3.4