机械优化设计题目答案

1、1- 1.简述优化设计问题数学模型的表达形式。答：优化问题的数学模型是实际优化设计问题的数学抽象。在明确设计变量、约束条件、目标函数之后，优化设计问题就可以表示成一般数学形式。求设计变量向量x=LxxLxA使f(x)Tmin且满足约束条件12nh(x)=0(k=1,2,Ll)kg.(x)<0(j=1,2,Lm)j利用可行域概念，可将数学模型的表达进一步简练。设同时满足gj(x)<0(j=1,2,Lm)和h(x)=0(k=1,2,L1)的设计点集合为R,即R为优化问题的可行域，则优化问题的数学模型可简练地写成k求x使minf(x)符号“G”表示“从属于”。xeR在实际优化问题中，对目

2、标函数一般有两种要求形式：目标函数极小化/(x)Tmin或目标函数极大化f(x)Tmax。由于求f(x)的极大化与求-f(x)的极小化等价，所以今后优化问题的数学表达一律采用目标函数极小化形式。1-2简述优化设计问题的基本解法。(不要抄书，要归纳)答：求解优化问题可以用解析解法，也可以用数值的近似解法。解析解法就是把所研究的对象用数学方程(数学模型)描述出来，然后再用数学解析方法(如微分、变分方法等)求出有化解。但是，在很多情况下，优化设计的数学描述比较复杂，因而不便于甚至不可能用解析方法求解；另外，有时对象本身的机理无法用数学方程描述，而只能通过大量试验数据用插值或拟合方法构造一个近似函数式

3、，再来求其优化解，并通过试验来验证；或直接以数学原理为指导，从任取一点出发通过少量试验(探索性的计算)，并根据试验计算结果的比较，逐步改进而求得优化解。这种方法是属于近似的、迭代性质的数值解法。数值解法不仅可用于求复杂函数的优化解，也可以用于处理没有数学解析表达式的优化问题。因此，它是实际问题中常用的方法，很受重视。其中具体方法较多，并且目前还在发展。但是，应当指出，对于复杂问题，由于不能把所有参数都完全考虑并表达出来，只能是一个近似的最后的数学描述。由于它本来就是一种近似，那么，采用近似性质的数值方法对它们进行解算，也就谈不到对问题的精确性有什么影响了。不管是解析解法，还是数值解法，都分别具

4、有针对无约束条件和有约束条件的具体方法。可以按照对函数倒数计算的要求，把数值方法分为需要计算函数的二阶导数、一阶导数和零阶导数(即只要计算函数值而不需计算其导数)的方法。2- 1.何谓函数的梯度？梯度对优化设计有何意义？答：二元函数f(xl,x2)在X0点处的方向导数的表达式可以改写成下面的形式ddxofcos91+6x1xo6x2cos02=r6f6f:cos01xo_6xl6x2_xocos02令Vf(x0)=呼=6x26f6f6x16x2并称它为函数f(X1，x2)在x0点处的梯度。xo假设d=方向导数6dcos01cos02xo为D方向上的单位向量，则有笙6dxo=Vf(x0)Td即函

5、数f(x1,x2)在xO点处沿某一方向d的等于函数在该点处的梯度Vf(x0)与d方向单位向量的内积。梯度方向是函数值变化最快的方向，而梯度的模就是函数变化率的最大值。梯度与切线方向d垂直，从而推得梯度方向为等值面的法线方向。梯度Vf(x0)方向为函数变化率最大方向，也就是最速上升方向。负梯度-Vf(x0)方向为函数变化率最小方向，即最速下降方向。2-2.求二元函数/(x1,x2)=2x：+x2-2x1+x2在x0=0,0丁处函数变化率最大的方向和数值。解；由于函数变化率最大的方向就是梯度的方向，这里用单位向量p表示，函数变化率最大和数值时梯度的模“Vf(x0)ll求f(x1，x2)在xO点处的

6、梯度方向和数值，计算如下:Vf(x0)=6x1f6x24x1-2=一22x2+1x01目(x0)p_lbf(xo)广_2=、5丄2-3.试求目标函数，x/=3x1-4Xx2在点Xo=l,0T处的最速下降方向，并求沿着该方向移动一个单位长度后新点的目标函数值。解：求目标函数的偏导数=6x一4x,=-4x+2xdx12dx1212则函数在Xo=1,Ot处的最速下降方向是这个方向上的单位向量是：P=-Vf(X0)=-dx2x1=1x12=0一6x+4x124x-2x12x1=1x12=0-64P-6,4t-3,2t'=P、石忑=云新点是1-亠13新点的目标函数值94f(X1=B2d3就称该点

7、集为凸集，否则为非凸集。函数f(x)为凸集定义域内的函数，若对任何的0-a-1及凸集域内的任意两点Xl、x2，存在如下不等式:fax+(1-a)x1<af(x)+(1-a)x1212对于约束优化问题minf(x)s.t.g.(x)<0(j=1,2,m)j若f(x)、g(x)j=l,2,.，m都是凸函数，则称此问题为凸规划。3- 1.简述一维搜索区间消去法原理。(要配图)答：搜索区间(a,b)确定之后，采用区间逐步缩短搜索区间，从而找到极小点的数值近似解。假设搜索区间a,b)内任取两点al，bl，albl,并计算函数值f(al),f(bl)。将有下列三种可能情形；1) f(al)f(

8、bl)由于函数为单谷，所以极小点必在区间(a,bl)内X2黄金分割法适用于a,b区间上的任何单谷函数求极小值问题。对函数除要求“单谷”外不作其他要求，甚至可以不连续。因此，这种方法的适应面相当广。黄金分割法也是建立在区间消去法原理基础上的试探方法，即在搜索区间a,b内适当插入两点a；、a2，并计算其函数值。a；、a2将区间分成三段。应用函数的单谷性质，通过函数值大小的比较，删去其中一段，使搜索区间得以缩短。然后再在保留下来的区间上作同样的处置，如此迭代下去，使搜索区间无限缩小，从而得到极小点的数值近似解。黄金分割法要求插入点a1、a2的位置相对于区间a，b两端点具有对称性，即a=b九(ba)1

9、a2=a+X(b-a)其中，九为待定常数。ab.1圮1予"除对称要求外黄分布。设原区间“新插入点a3应在貢金分割法还要求在保留下来的区间内再插a,b长度为l,如图a所示，保留下来的区间a，a2长度为，区间缩短率为九。为了保持相同的比例分"l1-X)位置上，a1取方程正数解，得入一点所形成的区间新三段，与原来区间的三段具有相同的比'可a,a2长度为，区间缩短率为"。为了保持相同的比例分在原区间的1九位置应相当于在保留区间的九2位置。故有1X=九2九2+X1=0<51«0.618若保留下来的区间为L,b,根据插入点的对称性，也能推得同样的九值。

10、所谓“黄金分割”是指将一线段分成两段的方法,使整段长与较长段的长度比值等于较长段与较短段长度的比鸟，：即1尢）按坐标点计算公式H=b一九（b一a）、a2=a+九（b一a）计算a1和a2，并计算其对应的函数值f），f）。12（3）根据区间消去法原理缩短搜索区间。为了能用原来的坐标点计算公式，需进行区间名称的代换，并在保留区间中计算一个新的试验点及其函数值。（4）检查区间是否缩短到足够小和函数值收敛到足够近，如果条件不满足则返回到步骤（2）。（5）如果条件满足，则取最后两试验点的平均值作为极小点的数值近似解。（6）黄金分割法的程序框图如图b所示。3-3.对函数f（a）=a2+2a_,当给定搜索区间

11、.-5<«<5时，写出用黄金分割法求极小点a*的前三次搜索过程。（要列表）解；此时的a=-5，b=5。首先插入两点al和a2。可得al=b尢（方一a）=-1.18,a2=a+尢（方一a）=1.18再计算相应插入点的函数值，得y1=f（a1）=-0.9676，y2=f（a2）=3.7524因为y2y1,所以消去区间a2,b，则新的搜索区间a,b的端点a=-5不变，而端点b=a2=1.18第一次迭代；此时插入点a1=b-尢（ba）=-2.639，a2=T.181。相应插入点的函数值y1=f（a1）=1.686，y2=f（a2）=-0.967，由于y1>y2，故消去区间a

12、，a1,新的搜索区间为-2.639,1.18,如此继续迭代下去列出前三次迭代结果黄金分割法的搜索过程迭代序号aa1a2bY1比较Y20-5-1181185-09676<375241-52639-11811181686-096722639-118-0279118-09676<-0483-2.639-1.737-1.181-0279-0.457-0.4823-4.使用二次插值法求f（x）=sin（x）在区间2,6的极小点，写出计算步骤和迭代公式，给定初始点X1=2，x?=4,X3=6，s=10-4»解1234x244554574554571x44554574736564721

13、252x666473656y0909297-0756802-0987572-0987572Jiy-0756802-0987572-0999708-09999612y-0279415-0279415-0279415-09997083x455457473656472125471236P|y1-0987572-0999708-0999961-1p迭代次数K=，极小点为一4.71236，最小值为一1yy=J+xx31I,=(x+xcc=T1xx23c=2xx212)C3收敛的条件：儿二爲y24-1简述无约束优化方法中梯度法、共轭梯度法、鲍威尔法的主要区别。答：梯度法是以负梯度方向作为搜索方向，使函数值

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

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

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

17、kk.=fxaVf(x)=minfxaVf(x)=minp(a)Ik根据一元函数极值的必要条件和多元复合函数求导公式得；k+11TkVf(x厂Vf(x')=0或写成fdk+八Tk、丿在最速下降法中，相邻两个迭代点上的函数梯度相互垂直。而搜索方向就是负梯度方向，因此相邻的两个搜索这就是说在最速下降法中，迭代点向函数极小点靠近的过程。由此可知，方向相互垂直。4-3.给定初始值xo二-7,llh使用牛顿法求函数f(x,x)=(x2)2+(x2x)2的极小值点和极小值。12112解：梯度函数、海赛矩阵分别为2(x2)+2(x2x)1124(x2x)12Vf(X1,X2)=2分)213V2f(X

18、1,X2)=121414144分)假设初始值x0=-7,11T则Vf(xo)=-76,(1分)2分)116X1=x0-V2f-1Vf(x0)=则Vf(x1)=1分)xi满足极值的必要条件，海赛矩阵是正定的，所以是极小点x11，f(x*)=-1。2分)4- 4以二元函数f(X1，X2)_为例说明单形替换法的基本原理。答：如图所示在平面上取不在同一直线上的三个点xl,x2,x3,以它们为顶点组成一单纯形。计算各顶点函数值，设f(xl)>f(x2)>f(x3)，这说明x3点最好，xl点最差。为了寻找极小点，一般来说。应向最差点的反对称方向进行搜索，即通过xl并穿过x2x3的中点x4的方向

19、上进行搜索。在此方向上取点x5使x5=x4+a(x4-x1)x5称作xl点相对于x4点的反射点，计算反射点的函数值f(X5),可能出现以下几种情形；1) f(x5)<f(x3)即反射点比最好点好要好，说明搜索方向正确，可以往前迈一步，也就是扩张。2) f(x3)<f(x5)<f(x2)即反射点比最好点差，比次差点好，说明反射可行，一反射点代替最差点构成新单纯形3) f(x2)<f(x5)f(xl),即反射点比次差点差，比最差点好，说明x5走的太远，应缩回一些，即收缩。答：直接解法通常适用于仅含不等式约束的问题，它的基本思路是在m个不等式约束条件所确定的可行域内选择一个初

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