《第五章约束优化2》由会员分享，可在线阅读，更多相关《第五章约束优化2（23页珍藏版）》请在金锄头文库上搜索。
1、5 约束非线性规划问题1. 数学模型形式：数学模型形式： min f （X） s.t. AXb （线性线性不等式约束）不等式约束） AeqX=beq （线性线性等式约束）等式约束） C(X)0 （非线性非线性不等式约束条件）不等式约束条件） Ceq(X)=0（非线性非线性等式约束）等式约束） Lb X Ub （边界约束条件）（边界约束条件）约约束束条条件件 2.使用格式：使用格式： x,fval,exitflag,output, grad,hessian= fmincon(fun,x0,A,b,Aeq,beq,Lb,Ub,Nlc,options,P1,P2)设置优化选项参数设置优化选项参数初始
2、点初始点调用目标函数的函数文件名调用目标函数的函数文件名目标函数在最优解的海色矩阵目标函数在最优解的海色矩阵返回目标函数在最优解的梯度返回目标函数在最优解的梯度优化算法信息的一个数据结构优化算法信息的一个数据结构返回算法的终止标志返回算法的终止标志返回目标函数的最优值返回目标函数的最优值返回目标函数的最优解返回目标函数的最优解附加参数附加参数非线性约束条件的函数名非线性约束条件的函数名设计变量的下界和上界设计变量的下界和上界线性等式约束的常数向量线性等式约束的常数向量线性等式约束的系数矩阵线性等式约束的系数矩阵线性不等式约束的常数向量线性不等式约束的常数向量线性不等式约束的系数矩阵线性不等式约
3、束的系数矩阵无定义时以空矩阵无定义时以空矩阵符号符号“ ”代替代替控制参数控制参数options序号功能默认值及其含义说明1输出形式0，无中间结果输出Options(1)=1,按照表格输出结果Options(1)=-1,隐藏警告信息2解x的精度1e-4Options(2)设置x解的终止条件3函数f的精度1e-4Options(3)设置函数f的终止条件4约束g的精度1e-6Options(4)设置约束g的终止条件5选择主要算法0Options(5)选择主要优化算法6搜索方向算法0fmin()函数为无约束优化搜索方向提供3种算法：Options(6)=0,拟牛顿法BFGS公式Options(6)=
4、1,拟牛顿法DFP公式Options(6)=2,梯度法7步长一维搜索0fmin()函数为无约束优化的步长一维搜索提供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)输入等式约
5、束个数14最大迭代次数100n（n为变量维数）Options(14)输入最大迭代次数15目标个数0Options(15)输入目标个数16差分步长最小值1e-8Options(16)步长的下限或变量的最小梯度值17差分步长最大值0.1Options(17)步长的上限或变量的最大梯度值18步长Options(18)步长参数，第1次迭代时置12.例题例题(1)： 求解二维约束优化问题 解解：(1)编制求解优化问题的编制求解优化问题的M文件。文件。 %建立M-函数文件myfun.m定义目标函数 function f=myfun(x) f=-x(1)*x(2)*x(3); %然后设置求解的初始点10;1
6、0;10,调用fmincon函数进行求解 x0=10;10;10; A=-1 -2 -2;1 2 2; b=0;72; optx,optf=fmincon(myfun,x0,A,b) 函数fminconfmincon(fun,x0,A,b,Aeq,beq,Lb,Ub,Nlc,options,P1,P2)mins.t.2.例题例题(2)： 求解二维约束优化问题解：解：（1）首先建立）首先建立M文件定义目标函数文件定义目标函数function f=objfun(x)f=exp(x(1)*(4*x(1)2+2*x(2)2+4*x(1)*x(2)+2*x(2)+1); (2)建立描述约束条件的建立描述
7、约束条件的M函数函数 confun.m定义该问题的非线性；定义该问题的非线性；function c,ceq=confun(x)c=1.5+x(1)*x(2)-x(1)-x(2);-x(1)*x(2)-10;%非线性不等式约束非线性不等式约束ceq=;%非线性等式约束非线性等式约束(3)编制求解优化问题的编制求解优化问题的M文件。文件。x0=-1,1; %设置求解过程中的初始点设置求解过程中的初始点x0 optx,optf=fmincon(objfun,x0,confun) %调用调用fmincon函数进函数进行求解行求解 函数fminconfmincon(fun,x0,A,b,Aeq,beq,
8、Lb,Ub,Nlc,options,P1,P2)mins.t.2.例题例题(2)： 求解二维约束优化问题解：解：（1）首先建立）首先建立M文件定义目标函数文件定义目标函数function f=objfun(x)f=exp(x(1)*(4*x(1)2+2*x(2)2+4*x(1)*x(2)+2*x(2)+1); (2)建立描述约束条件的建立描述约束条件的M函数函数 confun.m定义该问题的非线性；定义该问题的非线性；function c,ceq=confun(x)c=1.5+x(1)*x(2)-x(1)-x(2);-x(1)*x(2)-10;%非线性不等式约束非线性不等式约束ceq=;%非线
9、性等式约束非线性等式约束(3)编制求解优化问题的编制求解优化问题的M文件。文件。x0=-1,1; %设置求解过程中的初始点设置求解过程中的初始点x0lb=0,0;%设置下边界设置下边界ub=;%无上边界无上边界 optx,optf=fmincon(objfun,x0,lb,ub,confun)%调用调用fmincon函数函数进行求解进行求解 函数fminconfmincon(fun,x0,A,b,Aeq,beq,Lb,Ub,Nlc,options,P1,P2)mins.t.mins.t.单级直齿圆柱齿轮减速器的优化设计单级直齿圆柱齿轮减速器的优化设计n 在传统的设计方法中，按照相应的减速器在传
10、统的设计方法中，按照相应的减速器设计规范进行减速器的设计，对齿根弯曲强度、设计规范进行减速器的设计，对齿根弯曲强度、 齿面接触强度进行校核，用各种系数修正材料强度齿面接触强度进行校核，用各种系数修正材料强度和齿轮的载荷，并考虑齿轮精度的影响等因素；在和齿轮的载荷，并考虑齿轮精度的影响等因素；在这些步骤中，每一步都考虑安全系数，造成设计中这些步骤中，每一步都考虑安全系数，造成设计中的减速器往往的减速器往往体积过大，材料消耗增加体积过大，材料消耗增加等，导致成等，导致成本偏高。本偏高。 为解决这一问题，减速器的优化设计势为解决这一问题，减速器的优化设计势在必行。结合在必行。结合 matlabmat
11、lab 优化软件，对单级直齿圆柱优化软件，对单级直齿圆柱齿轮减速箱进行参数上的优化，为原始设计参数的齿轮减速箱进行参数上的优化，为原始设计参数的改进提供一定的修改依据。改进提供一定的修改依据。单级直齿圆柱齿轮减速器的优化设计单级直齿圆柱齿轮减速器的优化设计12示例示例13单级直齿圆柱齿轮减速器的优化设计单级直齿圆柱齿轮减速器的优化设计14（避免根切，其齿数不应小于（避免根切，其齿数不应小于17，限定，限定17x235）（传递动力，模数一般应大于（传递动力，模数一般应大于1 15mm-2mm5mm-2mm，限定，限定2 x110）（一定的载荷下，增大齿宽系数可减小齿轮直径和传动（一定的载荷下，增
12、大齿宽系数可减小齿轮直径和传动中心距，降低圆周速度；同时齿宽系数越大，载荷分中心距，降低圆周速度；同时齿宽系数越大，载荷分布越不均匀，限定布越不均匀，限定0.8x31.4）15mins.t.1718例例 某二级斜齿圆柱齿轮减速器，高速级输入功率某二级斜齿圆柱齿轮减速器，高速级输入功率P1=6.2kW，转，转速速n1=1450r/min；总传动比；总传动比i=31.5，齿轮宽度系数，齿轮宽度系数a=0.4，齿轮，齿轮材料和热处理：大齿轮材料和热处理：大齿轮45号钢正火号钢正火187207HB，小齿轮，小齿轮45号钢号钢调质调质228255HB，工作寿命，工作寿命10年以上。要求按照总中心距年以上
13、。要求按照总中心距a最最小来确定齿轮传动方案。小来确定齿轮传动方案。 解解:(1)建立优化设计的数学模型建立优化设计的数学模型 设计变量：设计变量： 将涉及总中心距将涉及总中心距a齿轮传动方案的齿轮传动方案的6个独立参数作为设计变量个独立参数作为设计变量 X=mn1,mn2,z1,z3,i1,T=x1,x2,x3,x4,x5,x6T 目标函数：目标函数： 减速器减速器总中心距总中心距a最小为目标函数最小为目标函数二、例题二、例题 约束条件约束条件:含性能约束和边界约束含性能约束和边界约束高速级齿轮接触强度条件高速级齿轮接触强度条件低速级齿轮接触强度条件低速级齿轮接触强度条件高速级大齿轮弯曲强度
14、条件高速级大齿轮弯曲强度条件低速级大齿轮弯曲强度条件低速级大齿轮弯曲强度条件大齿轮与轴不干涉大齿轮与轴不干涉低速级齿轮副模数的下限和上限低速级齿轮副模数的下限和上限高速级小齿轮齿数的下限和上限高速级小齿轮齿数的下限和上限高速级小齿轮齿数的下限和上限高速级小齿轮齿数的下限和上限低速级小齿轮齿数的下限和上限低速级小齿轮齿数的下限和上限高速级传动比的下限和上限高速级传动比的下限和上限齿轮副螺旋角的下限和上限齿轮副螺旋角的下限和上限性性 能能 约约 束束边边 界界 约约 束束X=mn1,mn2,z1,z3,i1,T=x1,x2,x3,x4,x5,x6T模模 数数二、例题二、例题(2)编制优化设计的编制
15、优化设计的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;%两级斜齿轮减速器优化设计的非线性不等式约束函数两级斜齿轮减速器优化设计的非线性不等式约束函数functiong,
16、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=;
17、文件名为：文件名为：Jsqyh_g.m文件名为：文件名为：Jsqyh_f.m二、例题二、例题(2)编制优化设计的编制优化设计的M文件文件 x,fval,exitflag,output, grad,hessian= fmincon(fun,x0,A,b,Aeq,beq,Lb,Ub,Nlc,options,P1,P2)x0=2;4;18;20;6.4;10;%设计变量的初始值设计变量的初始值lb=2;3.5;14;16;5.8;8;%设计变量的下限设计变量的下限ub=5;6;22;22;7;15;%设计变量的上限设计变量的上限x,fn=fmincon(jsqyh_f,x0,lb,ub,jsqyh_
18、g);disp *两级斜齿轮传动中心距优化设计最优解两级斜齿轮传动中心距优化设计最优解*fprintf(1, 高速级齿轮副模数高速级齿轮副模数 Mn1=%3.4fmmn,x(1)fprintf(1, 低速级齿轮副模数低速级齿轮副模数 Mn2=%3.4fmmn,x(2)fprintf(1, 高速级小齿轮齿数高速级小齿轮齿数 z1=%3.4fmmn,x(3)fprintf(1, 低速级小齿轮齿数低速级小齿轮齿数 z2=%3.4fmmn,x(4)fprintf(1, 高速级齿轮副传动比高速级齿轮副传动比 i1=%3.4fmmn,x(5)fprintf(1, 齿轮副螺旋角齿轮副螺旋角 beta=%3.
19、4fmmn,x(6)fprintf(1, 减速器总中心距减速器总中心距 a12=%3.4fmmn,fn) g=jsqyh_g(x); disp =最优点的性能约束函数值最优点的性能约束函数值=fprintf(1, 高速级齿轮副接触疲劳强度约束函数值高速级齿轮副接触疲劳强度约束函数值 g1=%3.4fmmn,g(1)fprintf(1, 低速级齿轮副接触疲劳强度约束函数值低速级齿轮副接触疲劳强度约束函数值 g2=%3.4fmmn,g(2)fprintf(1, 高速级大齿轮齿根弯曲强度约束函数值高速级大齿轮齿根弯曲强度约束函数值 g3=%3.4fmmn,g(3)fprintf(1, 低速级大齿轮齿
20、根弯曲强度约束函数值低速级大齿轮齿根弯曲强度约束函数值 g4=%3.4fmmn,g(4)fprintf(1, 大齿轮顶圆与轴不干涉几何约束函数值大齿轮顶圆与轴不干涉几何约束函数值 g5=%3.4fmmn,g(5) , , , ,文件名为文件名为eg9_6.m二、例题二、例题(3) M文件运行结果文件运行结果 *两级斜齿轮传动中心距优化设计最优解两级斜齿轮传动中心距优化设计最优解* 高速级齿轮副模数高速级齿轮副模数 Mn1=2.0461mm 低速级齿轮副模数低速级齿轮副模数 Mn2=3.6059mm 高速级小齿轮齿数高速级小齿轮齿数 z1=18.5156mm 低速级小齿轮齿数低速级小齿轮齿数 z
21、2=16.0000mm 高速级齿轮副传动比高速级齿轮副传动比 i1=5.8000mm 齿轮副螺旋角齿轮副螺旋角 beta=8.0000 减速器总中心距减速器总中心距 a12=317.4186mm =最优点的性能约束函数值最优点的性能约束函数值= 高速级齿轮副接触疲劳强度约束函数值高速级齿轮副接触疲劳强度约束函数值 g1=0.0000mm 低速级齿轮副接触疲劳强度约束函数值低速级齿轮副接触疲劳强度约束函数值 g2=0.0000mm 高速级大齿轮齿根弯曲强度约束函数值高速级大齿轮齿根弯曲强度约束函数值 g3=-1.0042mm 低速级大齿轮齿根弯曲强度约束函数值低速级大齿轮齿根弯曲强度约束函数值 g4=-15.1854mm 大齿轮顶圆与轴不干涉几何约束函数值大齿轮顶圆与轴不干涉几何约束函数值 g5=-761.8043mm2mm4mm19165.88338.292mm340mm9.838(4) 优化结果处理优化结果处理
