机械优化设计习题及答案.doc

1、机械 优化设计习题 及 参考答案 1-1.简述优化设计问题数学模型的表达形式。 答： 优化问题的数学模型是实际优化设计问题的数学抽象。在明确设计变量、约束条件、目标函数之后，优化设计问题就可以表示成一般数学形式。求设计变量向量 12 Tnx x x x L 使 ( ) minfx 且满足约束条件 ( ) 0 ( 1, 2 , )kh x k l L ( ) 0 ( 1, 2 , )jg x j m L 2-1.何谓函数的梯度？梯度对优化 设计有何意义？ 答： 二元函数 f(x1,x2)在 x0点处的方向导数的表达式可以改写成下面的形式 ： 2c o s 1c o s212c o s21c o

2、s1 xoxfxfxoxfxoxfxodf令xoTxfxfxfxfxf 2121)0( ， 则 称它为函数 f（ x1， x2）在 x0点处的梯度。 （ 1） 梯度方向 是函数值变化最快 方向， 梯度 模 是函数变化率的最大值。 （ 2） 梯度与切线方向 d 垂直，从而推得梯度方向为等值面的法线方向。梯度 )0(xf 方向为函数变化率最大方向，也就是最速上升方向。负梯度- )0(xf 方向为函数变化率最小方向，即最速下降方向。 2-2.求二元函数 f（ x1， x2） =2x12+x22-2x1+x2在 Tx 0,00 处函数变化率最 大的方向和数值。 解 ： 由于函数变化率最大的方向就是梯度

3、的方向，这里用单位向量 p表示，函数变化率最大和数值时梯度的模 )0(xf 。求 f（ x1， x2）在x0 点处的梯度方向和数值，计算如下： 1 20122 214210xxxxfxfxf 2221)0( xfxfxf = 5 5152512)0()0(xfxfp 2-3.试求目标函数 22212121 43, xxxxxxf 在点 X0=1,0T 处的最速下降方向，并求沿着该方向移动一个单位长度后新点的目标函数值。 解： 求目标函数的偏导数 212211 24,46 xxxfxxxf 则函数在 X0=1,0T处的最速下降方向是 462446)(012121012102121xxxxxxxx

4、xfxfXfP 这个方向上的单位向量是： 13 2,34)6( 4,6T22T PPe新点是 132133101 eXX 新点的目标函数值 1321394)( 1 Xf 2-4.何谓凸集、凸函数、凸规划？ （ 要求 配图） 答 ： 一个点集（或区域），如果连接其中任意两点 x1、 x2 的线段都全部包含在该集合内，就称该点集为凸集，否则为非凸集。 函数 f(x）为凸集定义域内的函数，若对任何的 01及 凸集 域内的任意两点 x1、 x2，存在如下不等式： 称 f（ x）是定义在图集上的一个凸函数。 对于约束优化问题 若 ( ) j =jf x g x、 （ ） 1,2,.,m 都是凸函数，则称

5、此问题为凸规划。 1 2 1 211f x x f x x 3-1.简述一维搜索区间消去法原理。（要配图） 答： 搜索区间（ a， b）确定之后，采用区间逐步缩短搜索区间，从而找到极小点的数值近似解。假设搜索区间（ a， b）内任取两点 a1， b1 ，a1 b1，并计算函数值 f（ a1）， f（ b1）。将有下列三种可能情形； 1） f（ a1） f（ b1）由于函数为单谷，所以极小点必在区间（ a， b1）内 2） f（ a1） f（ b1），同理，极小点应在区间（ a1， b）内 3） f（ a1） =f（ b1），这是极小点应在（ a1， b1）内 3-2.简述黄金分割法 搜索过程及

6、程序框图。 1 ()b b a 2 ()a b a 其中， 为待定常数。 3-3.对函数 2)( 2 f ，当给定搜索区间 55 时，写出用黄金 分割法求极小点 的前三次搜索过程。（要列表） 黄金分割法的搜索过程 序号 a a1 a2 b Y1 比较 Y2 0 -5 -1.18 1.18 5 -0.9676 -0.967 2 ？ -1.18 -0.279 1.18 -0.9676 -0.482 3-4.使用二次插值法求 f(x)=sin(x)在区间 2,6的极小点，写出计算步骤和迭代公式，给定初始点 x1=2， x2=4， x3=6， =10-4。 解： 1 2 3 4 x1 2 4 4.55

7、457 4.55457 x2 4 4.55457 4.73656 4.72125 x3 6 6 6 4.73656 y1 0.909297 -0.756802 -0.987572 -0.987572 y2 -0.756802 -0.987572 -0.999708 -0.999961 y3 -0.279415 -0.279415 -0.279415 -0.999708 xp 4.55457 4.73656 4.72125 4.71236 yp -0.987572 -0.999708 -0.999961 -1 迭代次数 K= 4 ，极小点为 4.71236 ，最小值为 -1 13131 xx y

8、yc ，12122 xx yyc ，32123 xx ccc )(21 3131 ccxxxp 收敛的条件： 22y yy p 4-1.简述无约束优化方法中梯度法、共轭梯度法、鲍威尔法的主要区别。 答：梯度法是以负梯度方向作为搜索方向，使函数值下降最快，相邻两个迭代点上的函数相互垂直即是相邻两个搜索方向相互垂直。这就是说在梯度法中，迭代点向函数极小点靠近的过程，走的是曲折的路线。这一次的搜索方向与前一次的搜索过程互相垂直，形成“之”字形的锯齿现象。从直观上可以看到，在远离极小点的位置，每次迭代可使函数值有较多的下降。可是在接近极小点的位置，由于锯齿现象使每次迭代行进的距离缩 短，因而收敛速度减

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

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

11、该点附近的范围下降最快。按此规律不断走步，形成以下迭代的算法 )(1 kxfkkxkx （ k=0， 1,2， ） 由于最速下降法是以负梯度方向作为搜索方向，所以最速下降法有称为梯度法 为了使目标函数值 沿搜索方向 - )( kxf 能获得最大的下降值，其步长因子ka应取一维搜索的最佳步长。即有 )(m i n)(m i n)(1 kxfakxfkxfkakxfkxf 根据一元函数极值的必要条件和多元复合函数求导公式得； 0)()1( kxfTkxf 或写成 01 kdTkd 由此可知，在最速下降法中，相邻两个迭代点上的函数梯度相互垂直。而搜索方向就是负梯度方向，因此相邻的两个搜索方向相互垂直

12、。这就是说在最速下降法中，迭代点向函数极小点靠近的过程。 4-3. 给 定初始值 x0=-7,11T ，使用牛顿法求函数2212121 )2()2(),( xxxxxf 的极小值点和极小值。 解： 梯度函数、海赛矩阵分别为 )2(4 )2(2)2(2),( 21 21121 xx xxxxxf （ 2 分） 4141 4121,84 44),( 12212 fxxf （ 4 分） 假设初始值 x0=-7,11T 则 ,11676)( 0 xf（ 1 分） 12)( 01201 xxx ff （ 2 分） 则 ,00)( 1 xf（ 1 分） x1 满足极值的必要条件，海赛矩阵是正定的，所以是极

13、小点 1)(,11 *1* xxx f。 （ 2 分） 4-4.以二元函数 ),( 21 xxf 为例说明单形替换法的基本原理。 答： 如图所示在平面上取不在同一直线上的三个点 x1， x2， x3，以它们为顶点组成一单纯形。 计算各顶点函数值，设 f（ x1） f（ x2） f（ x3） ， 这说明 x3 点最好， x1 点最差。 为了寻找极小点，一般来说。应向最差点的反对称方向进行搜索，即通过 x1 并穿过 x2x3的中点 x4 的方向上进行搜索。在此方向上取点 x5 使 x5=x4+ （ x4-x1） x5 称作 x1 点相对于 x4 点的反射点，计算反射点的函数值 f（ X5），可能出

14、现以下几种情形； 1） f（ x5） f(x1)，反射点比最差点还差，说明收缩应该多一些。将新点收缩在 x1x4 之间 5) f(x)f(x1)，说明 x1x4 方向上所有点都比最差点还要差，不能沿此方向进行搜索。 5-1.简述约束优化方法的分类。 （ 简述约束优化问题的直接解法、间接解法的原理、特点及主要方法。） 答 : 直接解法通常适用于仅含不等式约束的问题，它的基本思路是在 m 个不 等式约束条件所确定的可行域内选择一个初始点 0x ，然后决定可行搜索方向 d，且以适当的步长 沿 d 方向进行搜索，得到一个使目标函数值下降的可行的新点 1x ，即完成一个迭代。再以新点为起点，重复上述搜索

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

16、可行域内存在满足全部约束条件的点，且目标函数有定义。 直接解法有：随机方向法、复合形法、可 行方向法、广义简约梯度法等。 间接解法有不同的求解策略，其中一种解法的基本思路是将约束优化问题中的约束函数进行特殊的加权处理后，和目标函数结合起来，构成一个新的目标函数，即将原约束优化问题转化成一个或一系列的无约束优化问题。再对新的目标函数进行无约束优化计算，从而间接地搜索到原约束问题的最优解。 间接解法是目前在机械优化设计中得到广泛应用的一种有效方法。其特点是： 1）由于无约束优化方法的研究日趋成熟，已经研究出不少有效的无约束最优化方法和程序，使得间接解法有了可靠的基础。目前，这类算法的计算效率和数值

17、稳定 性也都有了较大提高。 2）可以有效地处理具有等式约束的约束优化问题。 3）间接算法存在的主要问题是，选取加权因子比较困难，加权因子选取不当，不但影响收敛速度和计算精度，甚至会导致计算失败。 间接解法有惩罚函数法和增广乘子法。 5-2.用内点法求下列问题的最优解： 03 12)( 21 12221 xgtsxxxxfm i n （提示：可构造惩罚函数 )(ln)(),( xgrxfrx，然后用解析法求解。） 解 构造内点惩罚函数： 2 1 )()(),( xrxxxxgrxfrx )3l n (12ln 令惩罚函数对 x的极值等于零： 0)3/()(2 22 xrxxdxd 得： 48366121rxx 舍去负根后，得 4 83662 rx 当 xxr 3130 2 该问题的最优解为，时， 。