基于MATLAB优化工具箱的优化计算分析报告.ppt

* * Optimization terminated successfully: Search direction less than 2*options. xopt = 1.0e+003 * 0.1000 1.2269 fopt = 100.0000 281.5296 3.5957 (4) M文件运算结果 9.6.1 函数fgoalattain 三、例题 (5) 优化结果处理 dd1 100mm 1227mm Ld 小带轮基准直径 带传动中心距 带的根数 1250mm 圆整 4 圆整 %优化结果数据处理后部分参数计算 Dd1=100;Dd2=Dd1*i;z=4;Ld=1250; v=pi*Dd1*n1/6e4; a1=Ld/4-pi*Dd1*(i+1)/8; a2=Dd1^2*(i-1)^2/8; a=a1+sqrt(a1^2-a2); alpha=180-180*Dd1*(i-1)/pi/a; disp disp ***************计算结果***************** fprintf(1, 小带轮基准直径 Dd1=%3.0fmm\n,Dd1); fprintf(1, 大带轮基准直径 Dd2=%3.0fmm\n,Dd2); fprintf(1, V带基准长度 Ld=%3.0fmm\n,Ld); fprintf(1, 传动中心距 a=%3.2fmm\n,a); fprintf(1, 小带轮包角 alpha=%3.2f度\n,alpha); fprintf(1, V带根数 z=%3.0fmm\n,z); (4) M文件运算结果 9.6.1 函数fgoalattain 三、例题 (5) 优化结果处理 ***************计算结果***************** 小带轮基准直径 Dd1=100mm 大带轮基准直径 Dd2=300mm V带基准长度 Ld=1250mm 传动中心距 a=293.82mm 小带轮包角 alpha=141.00度 V带根数 z= 4mm (4) M文件运算结果 9.6.1 函数fgoalattain 三、例题 (5) 优化结果处理 (6) 最终方案 9.6.2 函数fminimax min max {f1,f2,…,f3} s.t. AX≤b （线性不等式约束） AeqX=beq （线性等式约束） C(X)≤0 （非线性不等式约束条件） Ceq(X)=0 （非线性等式约束） Lb ≤X ≤Ub （边界约束条件） 一、多目标优化问题数学模型 各分目标函数 二、优化函数使用格式 [x,fval,exitflag,output, grad,hessian]= fminimax(@fun,x0, A,b,Aeq,beq,Lb,Ub,’Nlc’,options,P1,P2…) 设置优化选项参数 目标函数在最优解的海色矩阵 返回目标函数在最优解的梯度 优化算法信息的一个数据结构 返回算法的终止标志 返回目标函数的最优值 返回目标函数的最优解 附加参数 非线性约束条件的函数名 设计变量的下界和上界 线性等式约束的常数向量 线性等式约束的系数矩阵 线性不等式约束的常数向量 线性不等式约束的系数矩阵 无定义时以空矩阵 符号“[ ]”代替 9.6.2 函数fminimax 初始点 目标函数文件名 三、例题 已知直径为1单位长度的圆柱梁，要求将它制成矩形截面梁，满足重量最轻和强度最大的条件，试确定矩形截面尺寸。 解:(1)建立优化设计的数学模型 ①设计变量： 矩形截面的宽和高 X=[x1,x2]T ②目标函数： 重量→截面积： 弯曲强度→ 矩形截面矩量: 9.6.2 函数fmini