机械优化设计实例及matlab优化工具.pptVIP

* * * * * * * * 1.4.2 非函数文件 无函数头的M文件，由若干命令和注释构成。相当于主程序 如： %Filename is a sine.m x = 0:0.1:2*pi; y = sin(x); plot(x,y) 1.4.3 M-文件的操作 1.4.3 M-文件的操作 1.4.3 M-文件的操作 1.4.3 M-文件的操作 第2部分 优化计算工具 2.1 线性规划优化函数 2.2 无约束非线性优化函数 2.3 约束优化函数 MATLAB解决的线性规划问题的标准形式为： 2.1 线性规划优化函数 [函数] linprog [格式] x = linprog(f,A,b,Aeq,beq) x = linprog(f,A,b,Aeq,beq,lb,ub,x0) f——是优化参数x的系数矩阵； A——线性不等式约束系数阵 b——线性不等式约束常数向量 Aeq——线性等式约束系数阵 Beq——线性等式约束常数向量 2.1 线性规划及其优化函数 [应用举例] 求使函数 取最小值的x值， 且满足约束条件： 2.1线性规划及其优化函数 [代码] f = [-5; -4; -6]; A = [1 -1 1;3 2 4;3 2 0]; b = [20; 42; 30]; lb = zeros(3,1); [x,fval] = linprog(f,A,b,[],[],lb) [结果] x = 0.0000 15.0000 3.0000 fval = -78.0000 2.2 无约束非线性优化函数 [函数]fminunc [格式] x = fminunc(fun,x0) x = fminunc(fun,x0,options) [x,fval] = fminunc(…) 2.2 无约束非线性优化函数 [应用举例] 求 的最小值 [代码]%首先编写目标的.m文件 function f=myfun(x) f=3*x(1)^2+2*x(1)*x(2)+x(2)^2 %然后调用函数 fminunc x0=[1,1]; [x,fval]=fminunc(myfun,x0) 2.2 无约束非线性优化函数 [结果] x = 1.0e-008 * -0.7512 0.2479 fval = 1.3818e-016 2.3 约束优化函数 [函数] fmincon [格式] x = fmincon(fun,x0,A,b,Aeq,beq,lb,ub,nonlcon) x = fmincon(fun,x0,A,b,Aeq,beq,lb,ub,nonlcon,options) [x,fval] = fmincon(…) 2.3 约束优化函数—盖板优化实例： 目标函数： 约束： 盖板优化实例 盖板优化实例 盖板优化实例 运行结果： x = 0.6332 25.3264 fval = 101.3056 盖板优化实例 前面空心轴的问题： clear all x0=[23,19,4]; options=optimset(largescale,off,display,iter,tolx,1e-6); [x,fval,exitflag,output]=fmincon(myfun1,x0,[],[],[],[],[],[],confun1,options) function f=myfun1(x) f=6.12*(x(1)^2-x(2)^2)*x(3)*10e-6 x = 33.7505 12.8830 3.0000 fval = 0.1787 exitflag = 4