机械优化上机测试

1、WORD恪式fminsearch 函数(1) minf (X)=0. 1935x ix22x32 ( 4+ 6x4) X e R*2x32 (4+6x0 XeR5在MATLAB的Command窗口输入下列代码，即可求解：xO二x,fval=fminsearch(examplefsearch, xO)x=l.45920.71020. 00001.2960fval=5. 0671e-012(2) minf(X)=XxlXeR33 + cosx2+ logx3) /e1在MATLAB的Command窗口输入卜列代码，即可求解：xO二x, fval=fminsearch (ex,xO)x=-5. 36

2、185. 87690. 0006fval=-3. 4230e+004(3) minf(X)二2x3+4x3-10x3Xe R21X2ix2 + x31在MATLAB的Command窗口输入下列代码，即可求解：xO二x, fval=fminsearch(exs,xO)x 二1. 0e+026*-0. 0701-0. 0000-2.0289fval=-3.4489e+078函数 1 inprog用linprog函数求下列数学模型的最优解。minf (X)二一7xi 5x2s. t. 3xi+ X2W9O4xi+ 6X2W2OOxW20专业资料整理xW30x&O, x20在MATLAB的Co

3、mmand窗口输入下列代码，即可求解 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)Optimizationterminated20. 000020. 0000fval=-240. 0000quadprog 函数' + xj lOxi 16x2+89minf (X)=xls. t. Xi x2+ 12010 XiOXg 12010 x?20X2 Xi=0在MATLAB的Command窗口输入下列代码，即可求解H 二2,0;0,2;f二-10;-16;A 二-1,

4、1;b 二-1;lb 二;1;ub二10;10;Aeq二T, 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/ ( XiX2)s.t.30/ (xix2 3) -50W010 4x1x202Wxi W40. 5WxWl(1)编制非线性腿文件，并保存d/yconfun. m,代码如下: functionc, ceq=myconfun(x)c(l)=30/ (x(l)*x(2)"3)-50;c(2)二(4*x(l

5、)*x (2) )-10;ceq二;(2) 先编制目标函数的ni文件，并保存为myobjfun.m,代码如下：functionf二myobjfim (x)f 二 25/(x(l)*x(2 厂 3);(3) 在MATLAB的Command窗口输入下列代码，即可求解：x0=l;l;A 二;b 二;Aeq二; beq=;lb=2;0. 5;ub=4;l;x, fval, exitflag, output二fmincon(inyobjfim, xO, A, b, Aeq, beq, lb, ub, myconfun) x=2. 50001. 0000fval=10exitflag=loutput二it

6、erations:2funcCount:11stepsize:1algorithm:? medium-scale:SQP, Quasi-Newton, linesearch'firstorderopt:1. 8771e-007 cgiterations: message:lxl44char第一题：minf (x) xxx="1 Q3S +<s. t. 0X12x2x7223(1)先编制目标函数的ni文件，并保存为myobjfun. m,代码如下:functionf=myobjfun(x)f二-x(l)*x(2)*x(3);在MATLAB的Command窗口输入卜列代码，

7、即可求解：x0二A=-l, -2, -2; 1, 2, 2;b二0;72;x,fval, exitflag, output二fmincon(myobjfun, xO, A, b) x=1. 0e+013*7. 6034-1. 9009-1.9009fval=-2. 7473e+040exitflag=0output二iterations:74funcCount:304stepsize:2algorithm:' medium-scale:SQP, Quasi-Newton, line-search5firstorderopt:l 1128e+027cgiterations:message

8、:lx79charmin+ x； -xyx9 一 2xy -5x9s.t(Xi l)? + n 0 a第三题:2Xj -3x2 +6 20卩(1)先编制非线性絢文件，并保存n/yconfun. m,代码如下:functionc, ceq=myconfun(x) c二(x(l)T厂2-x(2);ceq二;(2)编制目标函数ni文件，并保存u/yobjfun. m,代码如下：f unc t i onf=myob j fun(x)f=x (1厂2+x(2厂2-x(1)*x(2)-2*x(1) -5*x(2);(3)在MATLAB的Command窗口输入下列代码，即可求解xO 二1;1; A二-2,

9、3;b二-6; Aeq二; beq二; lb 二; ub=;x, fval, exitflag, output =fmincon (myobjfun, xO, A, b, Aeq, beq, lb, ub, myconfun)1. 3333-1. 1111fval=7. 3827exitflag=0output二iterations:13funcCount:203stepsize:1. 2207e-004algorithm:' medium-scale:SQP, Quasi-Newton,line-search'firstorderopt:1.2948e+017cgiterat

10、ions:message:lx79char第四题:minf (X) =Xi2+x22 IOxi16x2+89s. t. Xi x2+ 12010 x&OX2 12010 x?20X2 Xi 二 0(1)编制目标函數m文件，并保存询objfun. m,代码如下：functionf二myobjfun (x)f二x(1厂2+x(2厂2-10*x(1)-16*x (2)+89;(2)在MATLAB的Command窗口输入下列代码，即可求解：xO 二A 二-1,1;b 二-1;Aeq二;beq二; lb 二;1;ub=10;10;x, fval, exit flag, output二 fmincon(niyobjfun, xO, A, b, Aeq, beq, lb, ub) x=7. 00006. 0000fval=8exitfla