资源描述:
《机械优化上机测试.doc》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、fminsearch函数(1)minf(X)=0.1935x1x22x32(4+6x4)X∈R4在MATLAB的Command窗口输入下列代码,即可求解:x0=[1;1;1;1];[x,fval]=fminsearch(@examplefsearch,x0)x=1.45920.71020.00001.2960fval=5.0671e-012(2)minf(X)=x13+cosx2+logx3)/ex1X∈R3在MATLAB的Command窗口输入下列代码,即可求解:x0=[1;1;1];[x,fval]=fminsearch(@ex,x0)
2、x=-5.36185.87690.0006fval=-3.4230e+004(3)minf(X)=2x13+4x1x23-10x1x2+x33X∈R3在MATLAB的Command窗口输入下列代码,即可求解:x0=[1;1;1];[x,fval]=fminsearch(@exs,x0)x=1.0e+026*-0.0701-0.0000-2.0289fval=-3.4489e+078函数linprog用linprog函数求下列数学模型的最优解。minf(X)=-7x1-5x2s.t.3x1+x2≤904x1+6x2≤200x1≤20x2≤30
3、x1≥0,x2≥0在MATLAB的Command窗口输入下列代码,即可求解:f=[-7;-5];A=[3,1;4,6;1,0;0,1];b=[90;200;20;30];lb=[0;0];[x,fval]=linprog(f,A,b,[],[],lb)Optimizationterminated.x=20.000020.0000fval=-240.0000quadprog函数minf(X)=x12+x22-10x1-16x2+89s.t.x1-x2+1≥010-x1≥0x2-1≥010-x2≥0x2-x1=0在MATLAB的Command窗
4、口输入下列代码,即可求解:H=[2,0;0,2];f=[-10;-16];A=[-1,1];b=[-1];lb=[;1];ub=[10;10];Aeq=[-1,1];beq=[0];[x,fval]=quadprog(H,f,A,b,Aeq,beq,lb,ub)x=0.29290.2929fval=-7.4437+89函数fmincon第一题:minf(X)=25/(x1x23)s.t.30/(x1x23)-50≤010-4x1x2≥02≤x1≤40.5≤x2≤1(1)编制非线性约束的.m文件,并保存为myconfun.m,代码如下:fun
5、ction[c,ceq]=myconfun(x)c(1)=30/(x(1)*x(2)^3)-50;c(2)=(4*x(1)*x(2))-10;ceq=[];(2)先编制目标函数的m文件,并保存为myobjfun.m,代码如下:functionf=myobjfun(x)f=25/(x(1)*x(2)^3);(3)在MATLAB的Command窗口输入下列代码,即可求解:x0=[1;1];A=[];b=[];Aeq=[];beq=[];lb=[2;0.5];ub=[4;1];[x,fval,exitflag,output]=fmincon(@m
6、yobjfun,x0,A,b,Aeq,beq,lb,ub,@myconfun)x=2.50001.0000fval=10exitflag=1output=iterations:2funcCount:11stepsize:1algorithm:'medium-scale:SQP,Quasi-Newton,line-search'firstorderopt:1.8771e-007cgiterations:[]message:[1x144char]第二题:mins.t.(1)先编制目标函数的m文件,并保存为myobjfun.m,代码如下:func
7、tionf=myobjfun(x)f=-x(1)*x(2)*x(3);(2)在MATLAB的Command窗口输入下列代码,即可求解:x0=[2;1;1];A=[-1,-2,-2;1,2,2];b=[0;72];[x,fval,exitflag,output]=fmincon(@myobjfun,x0,A,b)x=1.0e+013*7.6034-1.9009-1.9009fval=-2.7473e+040exitflag=0output=iterations:74funcCount:304stepsize:2algorithm:'mediu
8、m-scale:SQP,Quasi-Newton,line-search'firstorderopt:1.1128e+027cgiterations:[]message:[1x7