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

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

2、x2 xox2则称它为函数f (x, x)在x0点处的梯度。(1)梯度方向是函数值变化最快方向，梯度模是函数变化率的最大值。(2)梯度与切线方向d垂直，从而推得梯度方向为等值面的法线方向。梯度f(x0)方向为函数变化率最大方向，也就是最速上升方向。负梯度-f(x0)方向为函数变化 率最小方向，即最速下降方向_o o_2-2,求二元函数 f (x, x)=2x 1 +x -2x i+x在x。0,0t处函数变化率 最大的方向和数值。解：由于函数变化率最大的方向就是梯度的方向，这里用单位向量f(xo)。求 f (x1, x2)在 x02-3.试求目标函数f x1x3x24x1x2x2 在点x0=1,

3、0t处的最速卜降方向，并求沿着该方向移动一个单位长度后新点的目标函数值。解：求目标函数的偏导数6x1 4x2,_4x1 2x2x1 x2则函数在x0=1,0 t处的最速下降方向是p f( x0)这个方向上的单位向量是新点是xi6x1 4x24x12x2 文 210x2 x2；6,43,2t(6)2420 .1*3-3.对函数f()2 2,当给定搜索区间55时，写出用黄金分割法求极小点的前三次搜索过程。(要列表)黄金分割法的搜索过程序号aaia2b丫 1比较y 20-5-1.181.185-0.9676-0.9672?o-1.18-0.279?1.18-0.9676-0.4823-4.使用二次插

4、值法求f(x)=sin( x)在区间2,6的极小点，写出计算步骤和迭代公式，给定初始点 xi=2, x2=4, x3=6,e=10”4。 i234xi244.554574.55457x244.554574.736564.72i25x36664.73656yi0.909297-0.756802-0.987572-0.987572y2-0.756802-0.987572-0.999708-0.99996iy3-0.2794i5-0.2794i5-0.2794i5-0.999708xp4.554574.736564.72i254.7i236yp-0.987572-0.999708-0.99996i-

5、ic j ,c2 j ,c3 c-cij x3 xix2 xix2 x3解:迭代次数k= _4,极小点为4.71236 ,最小值为jxp 2 (xi x3收敛的条件：4- 1.简述无约束优化方法中梯度法、共苑梯度法、鲍威尔法的主要区别答：梯度法是以负梯度方向作为搜索方向，使函数值下降最快，相邻两个迭代点上的函数相互垂直即是相邻两个搜索方向相互垂直。这就是说在梯度法中，迭代点向函数极小点靠近的过程，走的是曲折的路线。这一次的搜索方向与前一次的搜索过程互相垂直，形成“之”字形的锯齿现象。从直观上可以看到，在远离极小点的位置，每次迭代可使函数值有较多的下降。可是在接近极小点的位置，由于锯齿现象使每次

6、迭代行进的距离缩短，因而收敛速度减慢。这种情况似乎与“最速下降”的名称矛盾，其实不然，这是因为梯度是函数的局部性质。从局部上看,在一点附近函数的下降是最快的，但从整体上看则走了许多弯路，因此函数的下降弁不算快。 共轲梯度法是共轲方向法中的一种，因为在该方法中每一个共轲的量都是依赖于迭代点处的负梯度而构造出来的，所以称作共轲梯度法。该方法的第一个搜索方向取作负梯度方向，这就是最速下降法。其余各步的搜索方向是将负梯度偏转一个角度，也就是对负梯度进行修正。所以共轲梯度法实质上是对最速下降法进行的一种改进，故它又被称作旋转梯度法。 鲍威尔法是直接利用函数值来构造共轲方向的一种共轲方向法，这种方法是在研

7、究具有正1定矩阵g的二次函数f(x)xtgx btx c的极小化问题时形成的。其基本思想是在不用2导数的前提下，在迭代中逐次构造 g的共轲方向。在该算法中，每一轮迭代都用连结始点和终点所产生出的搜索方向去替换原向量组中的第一个向量，而不管它的“好坏”，这是产生向量组线性相关的原因所在。因此在改进的算法中首先判断原向量组是否需要替换。如果需要替换，还要进一步判断原向量组中哪个向量最坏，然后再用新产生的向量替换这个最坏的向量，以保证逐次生成共轲方向。4-2.如何确定无约束优化问题最速下降法的搜索方向?答：优化设计是追求目标函数值最小，因此搜所方向d取该点的负梯度方向-f (x)。使函数值在该点附近

8、的范围下降最快。按此规律不断走步，形成以下迭代的算法k 1 kkx xf (x ) (k=0 , 1,2，)k由于最速下降法是以负梯度方向作为搜索方向，所以最速下降法有称为梯度法k为了使目标函数值沿搜索方向-f(x )能获得最大的下降值，其步长因子a应取一维搜k索的最佳步长。即有k 1kkkkf xf x a f (x ) min f x a f(x ) min ()k根据一元函数极值的必要条件和多元复合函数求导公式得；k 1 t kk 1 t kf (x ) f (x )0 或写成 dd 0由此可知，在最速下降法中，相邻两个迭代点上的函数梯度相互垂直。而搜索方向就是负梯度方向，因此相邻的两个

9、搜索方向相互垂直。这就是说在最速下降法中，迭代点向函数极小 点靠近的过 程。4-3.给定初始值x0=-7,11 t,使用牛顿法求函数f （xi,x2）（xi 2） 2 （xi2x2）2的极小值点和极小值。解：梯度函数、 海赛矩阵分别为f （xi,x2）2 心庄）2x2)2x22 ( x1 2 ) 2)112 4114 4假设初始值x0=-7,11 t76f （ x0）（ 1 分）116，x1x0f（ x0 ）（ 2 分）f （ x1）（ 1 分）x1 满足极值的必要条件，海赛矩阵是正定的，所以是极小点1 *1,f（ x ）（ 2 分）4-4. 以二元函数f （ x x2 ） 为例说明单形替换法

10、的基本原理。答：如图所示在平面上取不在同一直线上的三个点 x1, x2, x3, 以它们为顶点组成一单纯 形。计算各顶点函数值，设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 ）即反射点比最好点好要好，说明搜索方向正确，可以往前迈一步，

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

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

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