机械优化设计经典实例.pptVIP

4). 数组点正除（右除） 使两数组的对应元素正除 a./b 结果为: （a、b维数必须相同） 1.3.3 数组运算 第三十页，共五十五页。 1.4 源文件（M－文件） 分为两类： 函数文件和非函数文件 都用扩展名.M 1.4.1 函数文件（相当于子程序） 格式 function [输出表]＝函数名(输入表) 函数体 第三十一页，共五十五页。 1.4.2 非函数文件 无函数头的M文件，由若干命令和注释构成。相当于主程序 如： % 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(…) 第