机械优化设计题目答案

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

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

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

4、法，都分别具有针对无约束条件和有约束条件的具体方法。可以按照对函数倒数计算的要求，把数值方法分为需要计算函数的二阶导数、一阶导数和零阶导数(即只要计算函数值而不需计算其导数)的方法。2-1.何谓函数的梯度？梯度对优化设计有何意义?答：二元函数f(x1,x2)在x0点处的方向导数的表达式可以改写成下面的形式f苕|xo=fcose-互cX1xocx2cos22=xox1次2xo开令f(x0)=应=ffT并称它为函数f(x1,x2)在x0点处的梯度。开Ix1Fx2xo；2(x0)TJ即函数f(x1，x2)在x0点处沿某一方向d的xo假设d=cos911为d方向上的单位向量，则有3|lcosu2Fd方

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

6、下：f(x0)-f2f2=5(),;X1；X2:1匕1f(X0)_|1一p1f(X°)=一52-3.试求目标函数f（X1,X2）=3xj4x1X2+X2在点X0=1,0T处的最速下降方向，并求沿着该方向移动一个单位长度后新点的目标函数值。解：求目标函数的偏导数f=6X1-4x.Xi>2=-4Xi2X2则函数在X0=1,0T处的最速下降方向是这个方向上的单位向量是:-6,4T4,2T2-4.何谓凸集、凸函数、凸规划？（要求配图）一个点集（或区域），如果连接其中任意两点称f（X）是定义在图集上的一个凸函数。P一(242新点是X1=X0e=新点的目标函数值i94f(X)=-21313

7、Xi、X2的线段都全部包含在该集合内，就称该点集为凸集，否则为非凸集。函数f（x）为凸集定义域内的函数，若对任何的0<0t<1及凸集域内的任意两点x1、x2,存在如下不等式:f1：X1-1：工）X2<：.fX1|：1-1X2对于约束优化问题minf(x)s.t.gj(x)<0(j=1,2,川，m)若f(x)、gj(x)j=1,2,.,m都是凸函数，则称此问题为凸规划3-1.简述一维搜索区间消去法原理。(要配图)a,b)内答：搜索区间(a,b)确定之后，采用区间逐步缩短搜索区间，从而找到极小点的数值近似解。假设搜索区间(任取两点a1,1) f(a1)«f2) f

8、(a1)»f3) f(a1)=fb1,a1«b1,并计算函数值f(a1),f(b1)(b1)由于函数为单谷，所以极小点必在区间(b1),同理，极小点应在区间(a1,b)内(b1),这是极小点应在(a1,b1)内将有下列三种可能情形;a,bl)内fCbl)blftbl)hlb3-2.简述黄金分割法中0.618的来由，搜索过程及程序框图。黄金分割法适用于匕,b区间上的任何单谷函数求极小值问题。对函数除要求“单谷”外不作其他要求，甚至可以不连续。即在搜索区间b,b内适当插入值近似解。黄金分割法要求插入点因此，这种方法的适应面相当广。黄金分割法也是建立在区间消去法原理基础上的试探方

9、法,两点0(1、a2,并计算其函数值。0(1、a2将区间分成三段。应用函数的单谷性质，通过函数值大小的比较，删去其中一段,使搜索区间得以缩短。然后再在保留下来的区间上作同样的处置，如此迭代下去，使搜索区间无限缩小，从而得到极小点的数必久2的位置相对于区间A,b两端点具有对称性，即工1=b"(ba)a2=a+儿(ba)其中，九为待定常数。除对称要求外，黄金分割法还要求在保留下来的区间内再插入一点所形成的区间新三段，与原来区间的三段具有相同的比例分布。设原区间a,b和度为1,如图a所示，保留下来的区间1a,0(2长度为,区间缩短率为九为了保持相同的比例分布，新插入点仪3应在入(1九)位置

10、上，豆1在原区间的1一九位置应相当于在彳留区间的人2位置。故有取方程正数解，得5-10.6182若保留下来的区间为"，b,根据插入点的对称性，也能推得同样的九值。所谓“黄金分割”是指将一线段分成两段的方法，使整段长与较长段的长度比值等于较长段与较短段长度的比值，即1: '-.:(1-)同样算得九0.618°可见黄金分割法能使相邻两次搜索区间都具有相同的缩短率0.618,所以黄金分割法又被称作0.618法4给定口、尻s口二aa*St.结能)图b黄金分割法的搜索过程是：(1) 给出初始搜索区间b,b鹿收敛精度名，将九赋以0.618。(2) 按坐标点计算公式0tl=b九(

11、ba)、豆2=a+九(b-a)计算s和口2，并计算其对应的函数值f(%),f(a2)o(3) 根据区间消去法原理缩短搜索区间。为了能用原来的坐标点计算公式，需进行区间名称的代换，并在保留区间中计算一个新的试验点及其函数值。(4) 检查区间是否缩短到足够小和函数值收敛到足句近，如果条件不满足则返回到步骤(2)。(5) 如果条件满足，则取最后两试验点的平均值作为极小点的数值近似解。(6)黄金分割法的程序框图如图b所示。23-3.对函数f(a)=a+2a,当给定搜索区间-5<a<5时，写出用黄金分割法求极小点a的前三次搜索过程。(要列表)解；此时的a=-5,b=5首先插入两点a1和a2。

12、可得a1=b-1(b-a)=-1.18,a2=a+,(b-a)=1.18再计算相应插入点的函数值，得y1=f(a1)=-0.9676,y2=f(a2)=3.7524因为y2>y1,所以消去区间a2,b,则新的搜索区间a,b的端点a=-5不变，而端点b=a2=1.18第一次迭代；此时插入点a1=b-Mba)=-2.639,a2=-1.181。相应插入点的函数值y1=f(a1)=1.686,y2=f(a2)=-0.967,由于y1>y2,故消去区间a,a1,新的搜索区间为-2.639,1.18,如此继续迭代下去列出前三次迭代结果黄金分割法的搜索过程迭代序号aa1a2bY1比较Y20-5

13、-1.181.185-0.9676<3.75241-5-2.639-1.1811.181.686>-0.9672-2.639-1.18-0.2791.18-0.9676<-0.483-2.639-1.737-1.181-0.279-0.457>-0.4823-4.使用二次插值法求f(x)=sin(x)在区间2,6的极小点，写出计算步骤和迭代公式，给定初始点Xi=2,x?=4,X3=6,£=104解:1234Xi244.554574.55457X244.554574.736564.72125X36664.73656y0.909297-0.756802-0.987

14、572-0.987572V2-0.756802-0.987572-0.999708-0.999961y3-0.279415-0.279415-0.279415-0.999708Xp4.554574.736564.721254.71236yp-0.987572-0.999708-0.999961-1迭代次数K=4,极小点为4.71236,最小值为-1y3-y1y2-y1C2-GCi=,C2=，C3=X3-X1X2XX2-X31/C1.Xp=-(X1X3-)2C3收敛的条件:G的共钝方向。-Vf(x)。使函数值在该点附近的范围下降a应取一维搜索的最佳步长。即有kfixk1k1kfx-aIf(x)必

15、要条k和多元复合E根据一元函数极值k1=minfx-a'f(x的数求导公k得；k)=min.(二)4-1.简述无约束优化方法中梯度法、共辗梯度法、鲍威尔法的主要区别。答：梯度法是以负梯度方向作为搜索方向，使函数值下降最快，相邻两个迭代点上的函数相互垂直即是相邻两个搜索方向相互垂直。这就是说在梯度法中，迭代点向函数极小点靠近的过程，走的是曲折的路线。这一次的搜索方向与前一次的搜索过程互相垂直，形成“之”字形的锯齿现象。从直观上可以看到，在远离极小点的位置，每次迭代可使函数值有较多的下降。可是在接近极小点的位置，由于锯齿现象使每次迭代行进的距离缩短，因而收敛速度减慢。这种情况似乎与“最速下

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

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

18、以负梯度方向作为搜索方向，所以最速下降法有称为梯度法k为了使目标函数值沿搜索方向-环(X)能获得最大的下降值，其步长因子kf(xkTVf(xk)=0或写成idk411dk=0由此可知，在最速下降法中，相邻两个迭代点为勺函数梯假相互垂直。而搜索方向就是负梯度方向，因此相邻的两个搜索方向相互垂直。这就是说在最速下降法中，迭代点向函数极小点靠近的过程。4-3.给定初始值x0=-7,11:使用牛顿法求函数f(X1,X2)=(X12)2+(X1-2X2)2的极小值点和极小值。解：梯度函数、海赛矩阵分别为kf(X1,X2)二2(X1-2)2(X1-2X2)IL4(为2X2)(2分)V2f(x1,X2)=l

19、|4卜4假设初始值x0=-7,11则if(x。)1-76,116x答：直接解法通常适用于仅含不等式约束的问题，它的基本思路是在m个不等式约束条件所确定的可行域内选择一个初0始点x,然后决定可行搜索方向d,且以适当的步长a沿d方向进行搜索，得到一个使目标函数值下降的可行的新点即完成一个迭代。再以新点为起点，重复上述搜索过程，满足收敛条件后，迭代终止。所谓可行搜索方向是指，当设计点沿该方向作微量移动时，目标函数值将下降，且不会越出可行域。产生可行搜索方向的方法将由直接解法中的各种算法决定。直接解法的原理简单，方法实用。其特点是：1)由于整个求解过程在可行域内进行，因此迭代计算不论何时终点，都可以获

20、得一个比初始点好的设计点。2)若目标函数为凸函数，可行域为凸集，则可保证获得全域最优解。否则，因存在多个局部最优解，当选择的初始点不相同时，可能搜索到不同的局部最优解。为此，常在可行域内选择几个差别较大的初始点分别进=x0一.,2fr-f则5(xb=j。X1满足极值的必曩奈件，x=x1=,f(x4-4.以二元函薮f(x1海赛矩阵是正定的，所以是极小点彳。(2分)=1,x2)为例说明单形替换法的基本原理。答：如图所示在平面上取不在同一直线上的三个点x1,x2,x3,以它们为顶点组成一单纯形。计算各顶点函数值，设f(x1)>f(x2)>f(x3),这说明x3点最女f,x1点最差。为了寻

21、找极小点，一般来说。应向最差点的反对称方向进行搜索，即通过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(x2)<f(x5)<f(x1),即反射点比次差点差，比最差点好，说明x5走的太远，应缩回一些，即收缩。4)f(x5)>f(x1),反射点比最差点还差，说明收缩应该多一些。落新点收缩在x1x4之间5)f(x)>f(x1),说明x1x4方向上所有点都比最差点还要差，不能沿此方向进行搜索。行计算，以便从求得多个局部最优解中选择最好的最优解口约束条件的点，且目标函数有定义口3)要求可行域为有界