第四讲 基于MATLAB优化工具箱优化计算.ppt

本文观看结束！！！ 控制参数options 序号 功能 默认值及其含义 说明 1 输出形式 0，无中间结果输出 Options 1 1,按照表格输出结果 Options 1 -1,隐藏警告信息 2 解x的精度 1e-4 Options 2 设置x解的终止条件 3 函数f的精度 1e-4 Options 3 设置函数f的终止条件 4 约束g的精度 1e-6 Options 4 设置约束g的终止条件 5 选择主要算法 0 Options 5 选择主要优化算法 6 搜索方向算法 0 fmin 函数为无约束优化搜索方向提供3种算法： Options 6 0,拟牛顿法BFGS公式 Options 6 1,拟牛顿法DFP公式 Options 6 2,梯度法 7 步长一维搜索 0 fmin 函数为无约束优化的步长一维搜索提供2种算法： Options 7 0,二次和三次混合插值法 Options 7 1,三次多项式插值法 控制参数options 序号 功能 默认值及其含义 说明 8 函数值输出 Options 8 输出最终迭代函数值 9 梯度检验 0，不检验 Options 9 比较梯度 10 函数计算次数 Options 10 输出函数计算次数 11 梯度计算次数 Options 11 输出函数梯度计算次数 12 约束计算次数 Options 12 输出约束计算次数 13 等式约束个数 0，等式约束为0 Options 13 输入等式约束个数 14 最大迭代次数 100n （n为变量维数） Options 14 输入最大迭代次数 15 目标个数 0 Options 15 输入目标个数 16 差分步长 最小值 1e-8 Options 16 步长的下限或变量的最小梯度值 17 差分步长 最大值 0.1 Options 17 步长的上限或变量的最大梯度值 18 步长 Options 18 步长参数，第1次迭代时置1 例 某二级斜齿圆柱齿轮减速器，高速级输入功率P1 6.2kW，转速n1 1450r/min；总传动比i 31.5，齿轮宽度系数ψa 0.4，齿轮材料和热处理：大齿轮45号钢正火187~207HB，小齿轮45号钢调质228~255HB，工作寿命10年以上。要求按照总中心距a∑最小来确定齿轮传动方案。 解: 1 建立优化设计的数学模型 ①设计变量： 将涉及总中心距a∑齿轮传动方案的6个独立参数作为设计变量 X [mn1,mn2,z1,z3,i1,β]T [x1,x2,x3,x4,x5,x6]T ②目标函数： 减速器总中心距a∑最小为目标函数 二、例题 ③约束条件:含性能约束和边界约束 高速级齿轮接触强度条件 低速级齿轮接触强度条件 高速级大齿轮弯曲强度条件 低速级大齿轮弯曲强度条件 大齿轮与轴不干涉 低速级齿轮副模数的下限和上限 高速级小齿轮齿数的下限和上限 高速级小齿轮齿数的下限和上限 低速级小齿轮齿数的下限和上限 高速级传动比的下限和上限 齿轮副螺旋角的下限和上限 性 能 约 束 边 界 约 束 二、例题 2 编制优化设计的M文件 [x,fval,exitflag,output, grad,hessian] fmincon @fun,x0,A,b,Aeq,beq,Lb,Ub,’Nlc’,options,P1,P2… %两级斜齿轮减速器总中心距目标函数 function f jsqyh_f x ; hd pi/180; a1 x 1 *x 3 * 1+x 5 ; a2 x 2 *x 4 * 1+31.5/x 5 ; cb 2*cos x 6 *hd ; f a1+a2 /cb; %两级斜齿轮减速器优化设计的非线性不等式约束函数 function[g,ceq] jsqyh_g x ; hd pi/180; g 1 cos x 6 *hd ^3-3.079e-6*x 1 ^3*x 3 ^3*x 5 ; g 2 x 5 ^2*cos x 6 *hd ^3-1.701e-4*x 2 ^3*x 4 ^3; g 3 cos x 6 *hd ^2-9.939e-5* 1+x 5 *x 1 ^3*x 3 ^2; g 4 x 5 ^2.*cos x 6 *hd ^2-1.076e-4* 31.5+x 5 *x 2 ^3*x 4 ^2; g 5 x 5 * 2* x 1 +50 *cos x 6 *hd +x 1 *x 2 *x 3 -x 2 *x 4 * 31.5+x 5 ; ceq []; [ ] 文件名为： Jsqyh_g.m 文件名为： Jsqyh_f.m 二、例题 2 编制优化设计的M文件 [x,fval,exitflag,output, grad,hessian] fmincon @f