《机械优化设计》习题及答案

1、机械优化设计习题及参考答案1-1.简述优化设计问题数学模型的表达形式。答：优化问题的数学模型是实际优化设计问题的数学抽象。在明确设计变量、约束条件、目标函数之后，优化设计问题就可以表示成一般数学 形式。求设计变量向量x x1 x2 L xn ,使f(x) min且满足约束条件hk(x) 0 (k 1,2,Ll)gj(x) 0 (j 1,2,L m)2-1.何谓函数的梯度梯度对优化设计有何意义答：二元函数f( x1,x 2)在汽点处的方向导数的表达式可以改写成下面的形式:cos 1d xo x1 xof令 f(x0) Hx122cos 2 x2 xof _L tx1 x2 xo'f f

2、cos 1x1 x2 xo cos 2则称它为函数f (Xi, x2)在Xo点处的梯度(1)梯度方向是函数值变化最快方向，梯度模是函数变化率的最大值。(2)梯度与切线方向d垂直，从而推得梯度方向为等值面的法线方向。梯度f(x0)方向为函数变化率最大方向，也就是最速上升方向。负梯度-f(x0)方向为函数变化率最小方向，即最速下降方向2-2.求二元函数f (Xi, X2) =2xi2+X22-2xi+X2在Xo 0,0处函数变化率最 大的方向和数值。解：由于函数变化率最大的方向就是梯度的方向，这里用单位向量p表示，函数变化率最大和数值时梯度的模 | f (x0)| o求f (x1, x2)在x0点

3、处的梯度方向和数值，计算如下：f x074x1 2f2x2 1 x0 1x2f(x0)2,51,56x1x14x2, x24x1 2x22n _f(x0)1p | f(x0)|疾2-3.试求目标函数f x1,x2 3x12 4x1x2 x2在点X0=1,0 T处的最速下降方向，并求沿着该方向移动一个单位长度后新点的目标函数值。解：求目标函数的偏导数则函数在X0=1,0 T处的最速下降方向是P f(X0)这个方向上的单位向量是：P e IP新点是fx16x14x26f4x1 2x2 x1 14x2 0x9 x1 12 x2 06,4T 3,2T(6)242"3X1 X0 e新点的目标函

4、数值f(X1) 94 2.13132-4.何谓凸集、凸函数、凸规划(要求配图)答：一个点集(或区域)，如果连接其中任意两点x1、 包含在该集合，就称该点集为凸集，否则为非凸集。x2的线段都全部函数f(X )为凸集定义域的函数，若对任何的01及凸集域的任意两点x1、x2,存在如下不等式：f X1 1x2f x11x2称f (x)是定义在图集上的一个凸函数。对于约束优化问题若f(x)、gj(x)j=1,2,m都是凸函数，则称此问题为凸规划3-1.简述一维搜索区间消去法原理。(要配图)答：搜索区间(a, b)确定之后，采用区间逐步缩短搜索区间，从而找 到极小点的数值近似解。假设搜索区间(a, b)任

5、取两点al, bl , aibi,并计算函数值f (ai), f (bi)。将有下列三种可能情形；1) f (a) «f (bi)由于函数为单谷，所以极小点必在区间(a, bi)2) f (a)f (bi),同理，极小点应在区间(a, b)3) f (a) =f (bi),这是极小点应在(ai, bi)ffti)f门鱼1£(bn-IUvvjI-U立 Jbl b 耳 5 bl bgbl(b a)3-2.简述黄金分割法搜索过程及程序框图b (b a)3-3.对函数f( ) 2 2,当给定搜索区间 55时，写出用黄金分割法求极小点 的前三次搜索过程。(要列表)黄金分割法的搜索过程

6、序号aaia?bYi比较Y20-55<1-5>2<3>3-4.使用二次插值法求f(x)=sin( x)在区间2,6的极小点，写出计算 步骤和迭代公式，给定初始点 Xi=2, X2=4, X3=6,e=104。解：i234Xi24X24X3666yiy2y3xpyp-i迭代次数K= 4,极小点为,最小值为 二1一c V3yic '2yc C2CiCi , c2 , C3 X3XiX2XiX2X31 /Ci、Xp (Xi X3 ) 2C3收敛的条件：F2广4-1.简述无约束优化方法中梯度法、共钝梯度法、鲍威尔法的主要区别c答：梯度法是以负梯度方向作为搜索方向，使函数

7、值下降最快，相邻两个迭代点上的函数相互垂直即是相邻两个搜索方向相互垂直。这就是说在梯度法中，迭代点向函数极小点靠近的过程，走的是曲折的路线。这一次的搜索方向与前一次的搜索过程互相垂直，形成“之”字形 的锯齿现象。从直观上可以看到，在远离极小点的位置，每次迭代可使函数值有较多的下降。可是在接近极小点的位置，由于锯齿现象使每次迭代行进的距离缩短，因而收敛速度减慢。这 种情况似乎与“最速下降”的名称矛盾，其实不然，这是因为梯度是函数的局部性质。从局部 上看，在一点附近函数的下降是最快的，但从整体上看则走了许多弯路，因此函数的下降并不 算快。共轲梯度法是共轲方向法中的一种，因为在该方法中每一个共轲的量

8、都是依赖于迭代点处的负梯度而构造出来的，所以称作共轲梯度法。该方法的第一个搜索方向取作负梯度方向，这 就是最速下降法。其余各步的搜索方向是将负梯度偏转一个角度，也就是对负梯度进行修正。 所以共辗梯度法实质上是对最速下降法进行的一种改进，故它又被称作旋转梯度法。鲍威尔法是直接利用函数值来构造共轲方向的一种共轲方向法，这种方法是在研究其有正1 T T定矩阵G的一次函数f(x) -xTGx bTx C的极小化问题时形成的。 其基本思想是在不用导数的前提下，在迭代中逐次构造 G的共轲方向。在该算法中，每一轮迭代都用连结始点和终点所产生出的搜索方向去替换原向量组中的第一个向量，而不管它的“好坏” ，这是

9、产生向量组线性相关的原因所在。因此在改进的算法中首先判断原向量组是否需要替换。如果需要替换，还要进一步判断原向量组中哪个向量最坏，然后再用新产生的向量替换这个最坏的向量，以保证逐次生成共轭方向。4-2. 如何确定无约束优化问题最速下降法的搜索方向答：优化设计是追求目标函数值最小，因此搜所方向 d 取该点的负梯度方向 - f ( x) 。使函数值在该点附近的围下降最快。按此规律不断走步，形成以下迭代的算法k1 kkx xf (x ) (k=0, 1,2 ,)k由于最速下降法是以负梯度方向作为搜索方向，所以最速下降法有称为梯度法k为了使目标函数值沿搜索方向 - f ( x ) 能获得最大的下降值，

10、 其步长因子a 应取一维搜k索的最佳步长。即有k1 kkk kf x f x a f (x ) min f x a f(x ) min ( ) k根据一元函数极值的必要条件和多元复合函数求导公式得；k1T kk1Tkf (x ) f ( x ) 0 或写成 d d 0由此可知，在最速下降法中，相邻两个迭代点上的函数梯度相互垂直。而搜索方向就是负梯度方向，因此相邻的两个搜索方向相互垂直。这就是说在最速下降法中，迭代点向函数极小点靠近的过程。4-3. 给 定 初 始 值x0=-7,11T ， 使 用 牛 顿 法 求 函 数f (x1 ,x2) (x1 2)2 (x1 2x2)2 的极小值点和极小值

11、。解： 梯度函数、海赛矩阵分别为f (Xi , X2 )2(xi 2) 2(xi 2x2)4(Xi 2x2)2 f (Xi,X2)112f 12411 14 4(2分)(4分)假设初始值x0=-7,11则 f(x0)76 ,(1 分)116“cc 1C2X1 x0 2ff(x0) 1(2 分)则 f(x1)0，(1 分)x1满足极值的必要条件，海赛矩阵是正定的，所以是极小点*11 r */八、x X1, f(x )1。2 分14-4.以二元函数f(xX2)为例说明单形替换法的基本原理。答：如图所示在平面上取不在同一直线上的三个点x1, x2, x3,以它们为顶点组成一单纯形。计算各顶点函数值，

12、设 f (x1) >f (x2) >f (x3),这说明x3点最好，x1点最差。为了寻找极小点，一般来说。应向最差点的反对称方向进行搜索，即通过x1并穿过x2x3的中点x4的方向上进行搜索。在此方向上取点x5使 x5=x4+(x4-x1 )x5称彳x1点相对于X4点的反射点，计算反射点的函数值f (X5),可能出现以下几种情形；1) f (x5) <f (x3)即反射点比最好点好要好，说明搜索方向正确，可以往前迈一步， 也就是扩。2) f (x3) <f (x5) <f (x2)即反射点比最好点差，比次差点好，说明反射可行，一反射点代替最差点构成新单纯形3) f

13、(x2) <f (x5) <f(x1),即反射点比次差点差，比最差点好，说明 x5走的太远，应 缩回一些，即收缩。4) f(x5)>f(x1),反射点比最差点还差，说明收缩应该多一些。将新点收缩在 x1x4之间5) f(x)>f(x1),说明x1x4方向上所有点都比最差点还要差，不能沿此方向进行搜索。5-1.简述约束优化方法的分类。(简述约束优化问题的直接解法、间接解法的原理、特点及主要方法。)答：直接解法通常适用于仅含不等式约束的问题，它的基本思路是在 m个不等式约束条件0所确定的可行域选择一个初始点 X ,然后决定可行搜索方向 d,且以适当的步长 沿d方向进1行搜索

14、，得到一个使目标函数值下降的可行的新点x,即完成一个迭代。再以新点为起点，重复上述搜索过程，满足收敛条件后，迭代终止。所谓可行搜索方向是指，当设计点沿该方向作微量移动时，目标函数值将下降，且不会越出可行域。产生可行搜索方向的方法将由直接解法 中的各种算法决定。直接解法的原理简单，方法实用。其特点是：1)由于整个求解过程在可行域进行，因此迭代计算不论何时终点，都可以获得一个比初始点好的设计点。2)若目标函数为凸函数，可行域为凸集，则可保证获得全域最优解。否则，因存在多个局部最优解，当选择的初始点不相同时，可能搜索到不同的局部最优解。为此，常在可行域选择几个差别较大的初始点分别进行计算，以便从求得

15、多个局部最优解中选择最好的最优解。 3)要求可行域为有界的非空集，即在有界可行域存在满足全部约束条件的点，且目标函数有定义。直接解法有：随机方向法、复合形法、可行方向法、广义简约梯度法等。间接解法有不同的求解策略， 其中一种解法的基本思路是将约束优化问题中的约束函数进行特殊的加权处理后，和目标函数结合起来，构成一个新的目标函数，即将原约束优化问题转化成一个或一系列的无约束优化问题。再对新的目标函数进行无约束优化计算，从而间接地搜索到原约束问题的最优解。间接解法是目前在机械优化设计中得到广泛应用的一种有效方法。其特点是： 1 )由于无约束优化方法的研究日趋成熟， 已经研究出不少有效的无约束最优化方法和程序， 使得间接解法有了可靠的基础。目前，这类算法的计算效率和数值稳定性也都有了较大提高。 2 )可以有效地处理具有等式约束的约束优化问题。3)间接算法存在的主要问题是，选取加权因子比较困难，加权因子选取不当，不但影响收敛速度和计算精度，甚至会导致计算失败。间接解法有惩罚函数法和增广乘子法。5-2. 用点法求下列问题的最优解：2