机械优化上机测试.docxVIP

fminsearch 函数 （1）min f(X)=0.1935x 1 x2 2 x32（4＋6x4） X∈R4 2 x32（4＋6x4） X∈R4 在 MA TLAB 的 Command 窗口输入下列代码，即可求解： x0=[1;1;1;1]; [x,fval]=fminsearch(@examplefsearch,x0) x =1.4592 0.7102 0.0000 1.2960 fval = 5.0671e-012 （2） min f(X)= x 3＋cos x 2＋log x3） / e 1 x1 X∈R3 在 MA TLAB 的 Command 窗口输入下列代码，即可求解： x0=[1;1;1]; [x,fval]=fminsearch(@ex,x0) x = -5.3618 5.8769 0.0006 fval = -3.4230e+004 （3）min f(X)=2x 3＋ 4x 3 －10x 3 X ∈R3 1x2 1x2＋x3 1 在 MA TLAB 的 Command 窗口输入下列代码，即可求解： x0=[1;1;1]; [x,fval]=fminsearch(@exs,x0) x = 1.0e+026 * -0.0701 -0.0000 -2.0289 fval =-3.4489e+078 函数 linprog 用 linprog 函数求下列数学模型的最优解。 min f(X)= －7 x1－ 5 x 2 s.t. 3x1＋x2≤ 90 4x1＋6 x2≤ 200 x1≤ 20 x2≤ 30 x1≥ 0，x2≥ 0 在 MA TLAB 的 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) Optimization terminated. x = 20.0000 20.0000 fval =-240.0000 quadprog 函数 2＋x22－10 x1－16 x 2 ＋ 89 min f(X)= x 1 s.t. x1－x2＋1≥ 0 10－x1≥ 0 x2－1≥ 0 10－x2≥ 0 x2－x1=0 在 MA TLAB 的 Command 窗口输入下列代码，即可求解： 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.2929 0.2929 fval =-7.4437+89 函数 fmincon 第一题：min f(X)=25/ （x1 x2 3） s.t. 30/（x1 x2 3）－ 50≤ 0 10－4x1 x2≥ 0 2≤ x1≤ 4 0.5≤ x2≤ 1 （1）编制非线性约束的 .m 文件，并保存为myconfun.m ，代码如下： function [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 ，代码如下： function f=myobjfun(x) f=25/(x(1)*x(2)^3); (3) 在 MA TLAB 的 Command 窗口输入下列代码，即可求解： x0=[1;1]; A=[]; b=[]; Aeq=[]; beq=[]; lb=[2;0.5]; ub=[4;1]; [x,fval,exitflag,output]=fmincon(@myobjfun,x0,A,b,Aeq,beq,lb,ub,@myconfun) x = 2.5000 1.0000 fval =10 exitflag =1 output = iterations: 2 funcCount: 11 stepsize: 1 algorithm: medium-scale: SQP, Quasi-Newton, line-search firstorderopt: 1.8771e-007 cgiterations: [] message: [1x144 char] 第二题 ：min f ( x ) x x x 1 2 3 s.t. 0 x1 2 x 2 x 72 2 3 （1）先编制目标函数的 m 文件，并保存为 myobjfun.m ，代码如下： function f=myobjfun(x)