机械优化设计习题及答案1

1、. 机械优化设计习题及参考答案 简述优化设计问题数学模型的表达形式。1-1.在明确设计优化问题的数学模型是实际优化设计问题的数学抽象。答：优化设计问题就可以表示成一般数学目标函数之后，变量、约束条件、T? 形式。求设计变量向量使 xxx?Lxmin?f(x)n12 且满足约束条件 )0l(k?1,2,Lh(x)?k )m(j?1,2,g(x)?0Lj 何谓函数的梯度？梯度对优化设计有何意义？2-1.,xx点处的方向导数的表达式可以改写成下面x)答：二元函数f(在210?1cos?f?f?f?f?f? 的形式：?2?1?coscos? ?xoxocosxo2xo2xx1?x?1?x2?d?f?T

2、ff? ?1x? ，令?0)?f(x? f?xo2x1?x? 2x? ）在x点处的梯度。x则称它为函数f（，x012 ）梯度方向是函数值变化最快方向，梯度模是函数变化率的最大值。（1从而推得梯度方向为等值面的法线方向。d垂直，（2）梯度与切线方向负梯度方向为函数变化率最大方向，也就是最速上升方向。梯度)0(x?f 方向为函数变化率最小方向，即最速下降方向。-)0(x?f 22 在=2x+x处函数变化率最-2x+x），求二元函数2-2.f（xxT00x?,2211210 大的方向和数值。p解：由于函数变化率最大的方向就是梯度的方向，这里用单位向量表示，函数变化率最大和数值时梯度的模。求f（x1，

3、x2）在)(f?x0. . 点处的梯度方向和数值，计算如下：x0f?2?2?4x1? 1?x? ?x0?f?f?10?1x2x2? 2x?22ff? =5?(x0)?f?21?x?x?2?2?1)x0?f(?5?p? 1)0?f(x5?5? ?T0X处的最速下降在点=1,0 2-3.试求目标函数22x?3x4?fxx,xx?221112 方向，并求沿着该方向移动一个单位长度后新点的目标函数值。 解：求目标函数的偏导数 f?f? x24x?4?6x?x,? 2121xx?21T0=X 1,0则函数在处的最速下降方向是f? xx?4?66?x? 210?f(X)?P?1?f?x24x?4?1?x1

4、210x? ?2x?x?1210?x2 这个方向上的单位向量是：TT?3,4P?6,2 ?e?P132246)?(? 新点是3?1?13 01?e?XX?2?13? 新点的目标函数值94 113?Xf()?2 13. . 2-4.何谓凸集、凸函数、凸规划？（要求配图） 答：一个点集（或区域），如果连接其中任意两点x1、x2的线段都全部包含在该集合内，就称该点集为凸集，否则为非凸集。 及凸集域内为凸集定义域内的函数，若对任何的函数f(x）?10? ，存在如下不等式：x1、x2的任意两点?xxx?f1?xf1? 2121 x）是定义在图集上的一个凸函数。（称f 对于约束优化问题. . f(x)、g

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

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

7、56 4.72125 x 36 6 6 4.73656 y 10.909297 -0.756802 -0.987572 -0.987572 y 2-0.756802 -0.987572 -0.999708 -0.999961 y 3-0.279415 -0.279415 -0.279415 -0.999708 x p4.55457 4.73656 4.72125 4.71236 y p-0.987572 -0.999708 -0.999961 -1 迭代次数K= 4 ，极小点为 4.71236 ，最小值为 -1 y?yy?yc?c， 131221?c?cc?132 x?xx?xx?x13321

8、21c 1)?(x?x?x3p1 2c3y?y p2 收敛的条件： ?y2. . 4-1.简述无约束优化方法中梯度法、共轭梯度法、鲍威尔法的主要区别。 答：梯度法是以负梯度方向作为搜索方向，使函数值下降最快，相邻两个迭代点上的函数相互垂直即是相邻两个搜索方向相互垂直。这就是说在梯度法中，迭代点向函数极小点靠近的过程，走的是曲折的路线。这一次的搜索方向与前一次的搜索过程互相垂直，形成“之”字形的锯齿现象。从直观上可以看到，在远离极小点的位置，每次迭代可使函数值有较多的下降。可是在接近极小点的位置，由于锯齿现象使每次迭代行进的距离缩短，因而收敛速度减慢。这种情况似乎与“最速下降”的名称矛盾，其实不

9、然，这是因为梯度是函数的局部性质。从局部上看，在一点附近函数的下降是最快的，但从整体上看则走了许多弯路，因此函数的下降并不算快。 共轭梯度法是共轭方向法中的一种，因为在该方法中每一个共轭的量都是依赖于迭代点处的负梯度而构造出来的，所以称作共轭梯度法。该方法的第一个搜索方向取作负梯度方向，这就是最速下降法。其余各步的搜索方向是将负梯度偏转一个角度，也就是对负梯度进行修正。所以共轭梯度法实质上是对最速下降法进行的一种改进，故它又被称作旋转梯度法。 鲍威尔法是直接利用函数值来构造共轭方向的一种共轭方向法，这种方法是在研究其有正1TTx?cGx(x)?bxfG的极小化问题时形成的。定矩阵其基本思想是在

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

11、?)（k=0，1,2，） k 由于最速下降法是以负梯度方向作为搜索方向，所以最速下降法有称为梯度法ka?f(x)应取一维搜 为了使目标函数值沿搜索方向其步长因子-能获得最大的下降值，k 索的最佳步长。即有k?1kkkk?)min(x?fmin?f(x)?()?fx?afx?afx? ?k? 根据一元函数极值的必要条件和多元复合函数求导公式得；k?1Tkk?1Tk?f(x0?dd?0()?fx) 或写成?由此可知，在最速下降法中，相邻两个迭代点上的函数梯度相互垂直。而搜索方向就是负梯度方向，因此相邻的两个搜索方向相互垂直。这就是说在最速下降法中，迭代点向函数极小 点靠近的过程。 0T数函求法顿值

12、初给4-3.定始牛用使，=-7,11x . . 的极小值点和极小值。 22)?22)x?(x)f(x,x?(x?21211 解： 梯度函数、海赛矩阵分别为 2(x?2)?2(x?2x)?112 （2分） ?),(xx?f?21?4(x?2x)?2111?4?4? 421?22 （4分）?f?,?(x,x?)f?11218?4? 44?T 0=-7,11假设初始值x76?0 分） （1则 ,)f(x?116?2?1?0210 分） （2?x?x?f?f(x)?1?0?1 分） （1则 ,)(?fx?0?1 x满足极值的必要条件，海赛矩阵是正定的，所以是极小点1?*1 （2分） 。 1f)?(x?

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

14、扩张。 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所确定的可行域内选择一个初始点然后决定可行搜索方向，d1x，即完成一个

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

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