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

《机械优化设计》复习题解答一、填空题1、用最速下降法求f(X)=100(x2-x12)2+(1-x1)2的最优解时，设X（0）＝[-0.5,0.5]T，第一步迭代的搜索方向为[-47,-50]T。2、机械优化设计采用数学规划法，其核心一是寻找搜索方向，二是计算最优步长。3、当优化问题是凸规划的情况下，任何局部最优解就是全域最优解。4、应用进退法来确定搜索区间时，最后得到的三点，即为搜索区间的始点、中间点和终点，它们的函数值形成高－低－高趋势。5、包含n个设计变量的优化问题，称为n维优化问题。6、函数CXBHXXTT21的梯度为B。7、设G为n×n对称正定矩阵，若n维空间中有两个非零向量d0，d1，满足(d0)TGd1=0，则d0、d1之间存在共轭关系。8、设计变量、目标函数、约束条件是优化设计问题数学模型的基本要素。9、对于无约束二元函数),(21xxf，若在),(x20100xx点处取得极小值，其必要条件是∇f(x10,x20)=0，充分条件是∇2f(x10,x20)=0正定。10、K-T条件可以叙述为在极值点处目标函数的梯度为起作用的各约束函数梯度的非负线性组合。11、用黄金分割法求一元函数3610)(2xxxf的极小点，初始搜索区间]10,10[],[ba，经第一次区间消去后得到的新区间为[-2.3610]。12、优化设计问题的数学模型的基本要素有设计变量、目标函数、约束条件。13、牛顿法的搜索方向dk=，其计算量大，且要求初始点在极小点附近位置。14、将函数f(X)=x12+x22-x1x2-10x1-4x2+60表示成CXBHXXTT21的形式12[x1x2][2−1−12][x1x2]+[−10−4][x1x2]+60。15、存在矩阵H，向量d1，向量d2，当满足d1THd2=0，向量d1和向量d2是关于H共轭。16、采用外点法求解约束优化问题时，将约束优化问题转化为外点形式时引入的惩罚因子r数列，具有单调递增特点。17、采用数学规划法求解多元函数极值点时，根据迭代公式需要进行一维搜索，即求最优步长。1kkHg二、选择题1、下面C方法需要求海赛矩阵。A、最速下降法B、共轭梯度法C、牛顿型法D、DFP法2、对于约束问题2212221122132min44g10g30g0fXxxxXxxXxXx根据目标函数等值线和约束曲线，判断1[1,1]TX为，251[,]22TX为。DA．内点；内点B.外点；外点C.内点；外点D.外点；内点3、内点惩罚函数法可用于求解B优化问题。A无约束优化问题B只含有不等式约束的优化问题C只含有等式的优化问题D含有不等式和等式约束的优化问题4、对于一维搜索，搜索区间为[a，b]，中间插入两个点a1、b1，a1b1，计算出f(a1)f(b1)，则缩短后的搜索区间为D。A[a1，b1]B[b1，b]C[a1，b]D[a，b1]5、D不是优化设计问题数学模型的基本要素。A设计变量B约束条件C目标函数D最佳步长6、变尺度法的迭代公式为xk+1=xk-αkHk▽f(xk)，下列不属于Hk必须满足的条件的是C。A.Hk之间有简单的迭代形式B.拟牛顿条件C.与海塞矩阵正交D.对称正定7、函数)(Xf在某点的梯度方向为函数在该点的A。A、最速上升方向B、上升方向C、最速下降方向D、下降方向8、下面四种无约束优化方法中，D在构成搜索方向时没有使用到目标函数的一阶或二阶导数。A梯度法B牛顿法C变尺度法D坐标轮换法9、设)(Xf为定义在凸集R上且具有连续二阶导数的函数，则)(Xf在R上为凸函数的充分必要条件是海塞矩阵G(X)在R上处处B。A正定B半正定C负定D半负定10、下列关于最常用的一维搜索试探方法——黄金分割法的叙述，错误的是D，假设要求在区间[a，b]插入两点α1、α2，且α1α2。A、其缩短率为0.618B、α1=b-λ（b-a）C、α1=a+λ（b-a）D、在该方法中缩短搜索区间采用的是外推法。11、与梯度成锐角的方向为函数值A方向，与负梯度成锐角的方向为函数值B方向，与梯度成直角的方向为函数值C方向。A、上升B、下降C、不变D、为零12、二维目标函数的无约束极小点就是B。A、等值线族的一个共同中心B、梯度为0的点C、全局最优解D、海塞矩阵正定的点13、最速下降法相邻两搜索方向dk和dk+1必为B向量。A相切B正交C成锐角D共轭14、下列关于内点惩罚函数法的叙述，错误的是A。A可用来求解含不等式约束和等式约束的最优化问题。B惩罚因子是不断递减的正值C初始点应选择一个离约束边界较远的点。D初始点必须在可行域内三、问答题（看讲义）1、试述两种一维搜索方法的原理，它们之间有何区别？2、惩罚函数法求解约束优化问题的基本原理是什么？3、试述数值解法求最佳步长因子的基本思路。4、试述求解无约束优化问题的最速下降法与牛顿型方法的优缺点。5、写出用数学规划法求解优化设计问题的数值迭代公式，并说明公式中各变量的意义，并说明迭代公式的意义。6、什么是共轭方向？满足什么关系？共轭与正交是什么关系？四、解答题1、试用梯度法求目标函数f(X)=1.5x12+0.5x22-x1x2-2x1的最优解，设初始点x(0)=[-2，4]T,选代精度ε=0.02（迭代一步）。解：首先计算目标函数的梯度函数∇f=[3∗x1−x2−2x2−x1],计算当前迭代点的梯度向量值∇f(X(0))=[−3∗2−4−24+2]=[−126]梯度法的搜索方向为S（k）=−∇f,因此在迭代点x(0)的搜索方向为[12，－6]T在此方向上新的迭代点为：X（k+1）=X（k）+αS（k）=X（0）+αS（0）=[−24]+α[12−6]=[−2+12α4−6α]把新的迭代点带入目标函数，目标函数将成为一个关于单变量α的函数F(α)f(X(k+1))=f([−2+12α4−6α])=1.5(−2+12α)2+0.5(4−6α)2−(−2+12α)(4−6α)−2(−2+12α)＝F(α)令dF(α)dα=−180+612α＝0，可以求出当前搜索方向上的最优步长α＝517≈0.2941新的迭代点为X（0）+αS（0）＝[1.52922.2354]当前梯度向量的长度‖∇f‖=√12x12+6x6=13.4164ε,因此继续进行迭代。第一迭代步完成。2、试用牛顿法求f(X)=(x1-2)2+(x1-2x2)2的最优解，设初始点x(0)=[2,1]T。解1：（注：题目出题不当，初始点已经是最优点，解2是修改题目后解法。）牛顿法的搜索方向为S(k)=−∇2(f)−1∇(f),因此首先求出当前迭代点x(0)的梯度向量、海色矩阵及其逆矩阵∇(f)＝[4∗x1−4∗x2−48∗x2−4∗x1]∇(f（x(0)）)=[00]∇2(f)＝[4−4−48]∇2(f)−1=14[2111]S(k)=−∇2(f)−1∇(f)=[00]不用搜索，当前点就是最优点。解2：上述解法不是典型的牛顿方法，原因在于题目的初始点选择不当。以下修改求解题目的初始点，以体现牛顿方法的典型步骤。以非最优点x(0)=[1,2]T作为初始点，重新采用牛顿法计算牛顿法的搜索方向为S(k)=−∇2(f)−1∇(f),因此首先求出当前迭代点x(0)的梯度向量、以及海色矩阵及其逆矩阵梯度函数：∇(f)＝[4∗x1−4∗x2−48∗x2−4∗x1]初始点梯度向量：∇(f（x(0)）)=[−812]海色矩阵：∇2(f)＝[4−4−48]海色矩阵逆矩阵：∇2(f)−1=14[2111]当前步的搜索方向为：S(k)=−∇2(f)−1∇(f)＝−14[2111][−812]＝[−11]新的迭代点位于当前的搜索方向上：X（k+1）=X（k）+αS（k）=X（0）+αS（0）＝[12]＋α[−11]＝[1−α2＋α]把新的迭代点带入目标函数，目标函数将成为一个关于单变量α的函数F(α)f(X(k+1))=f([1−α2＋α])=(α+1)2+(3α+3)2＝F(α)令dF(α)dα=20α+20＝0，可以求出当前搜索方向上的最优步长α＝−1新的迭代点为X(1)=X（0）+αS（0）＝[12]–[−11]＝[21]当前梯度向量的长度‖∇f‖=√12x12+8x8=14.4222ε,因此继续进行迭代。第二迭代步：∇(f)＝[4∗x1−4∗x2−48∗x2−4∗x1]∇(f（x(1)）)=[00]‖∇f‖=0ε因此不用继续计算，第一步迭代已经到达最优点。这正是牛顿法的二次收敛性。对正定二次函数，牛顿法一步即可求出最优点。3、设有函数f(X)=x12+2x22-2x1x2-4x1，试利用极值条件求其极值点和极值。解：首先利用极值必要条件∇(f)=[00]找出可能的极值点：令∇(f)＝[2∗x1−2∗x2−44∗x2−2∗x1]＝[00]求得[x1x2]=[42],是可能的极值点。再利用充分条件∇2(f)正定（或负定）确认极值点。∇2(f)＝[2−2−24]|2|＝20|2−2−24|＝8−4＝40因此∇2(f)正定,X∗=[x1x2]=[42]是极小点，极值为f(X*)=-84、求目标函数f(X)=x12+x1x2+2x22+4x1+6x2+10的极值和极值点。解法同上5、试证明函数f(X)=2x12+5x22+x32+2x3x2+2x3x1-6x2+3在点[1，1，-2]T处具有极小值。解：必要条件：∇(f)＝[4∗x1+2∗x310∗x2+2∗x3−62∗x1+2∗x2+2∗x3]将点[1，1，-2]T带入上式，可得∇(f)＝[000]充分条件∇2(f)＝[4020102222]|4|＝40|40010|＝400|4020102222|＝80−40−16＝240∇2(f)正定。因此函数在点[1，1，-2]T处具有极小值6、给定约束优化问题minf(X)=(x1-3)2+(x2-2)2s.t.g1(X)=－x12－x22＋5≥0g2(X)=－x1－2x2＋4≥0g3(X)=x1≥0g4(X)=x2≥0验证在点TX]2[，１Kuhn-Tucker条件成立。解：首先，找出在点TX]2[，１起作用约束：g1(X)＝0g2(X)＝0g3(X)＝2g4(X)＝1因此起作用约束为g1(X)、g2(X)。然后，计算目标函数、起作用约束函数的梯度，检查目标函数梯度是否可以表示为起作用约束函数梯度的非负线性组合。∇(f)＝[2∗x1−62∗x2−4]＝[−2−2]∇(g1)＝[−2∗x1−2∗x2]＝[−4−2],∇(g2)＝[−1−2]求解线性组合系数∇(f)=λ1∇(g1)+λ2∇(g2)[−2−2]=λ1[−4−2]+λ2[−1−2]得到λ1＝13,λ2=23,均大于0因此在点TX]2[，１Kuhn-Tucker条件成立7、设非线性规划问题01)(0)(0)(..)2()(min2221322112221xxXgxXgxXgtsxxXf用K-T条件验证TX0,1*为其约束最优点。解法同上8、已知目标函数为f(X)=x1+x2，受约束于：g1(X)=-x12+x2≥0g2(X)=x1≥0写出内点罚函数。解：内点罚函数的一般公式为其中：r(1)r(2)r(3)…r(k)…0是一个递减的正值数列r(k)＝Cr(k-1)，0＜C＜1因此罚函数为：∅(X,r(k))=x1+x2+r(k)(1−x12+x2+1x1)9、已知目标函数为f(X)=(x1-1)2+(x2+2)2受约束于：g1(X)=-x2-x1-1≥0g2(X)=2-x1-x2≥0g3(X)=x1≥0g4(X)=x2≥0试写出内点罚函数。解法同上10、如图，有一块边长为6m的正方形铝板，四角截去相等的边长为x的方块并折转，造一个无盖的箱子，问如何截法（x取何值）才能获得最大容器的箱子。试写出这一优化问题的数学模型以及用MATLAB软件求解的程序。11、某厂生产一个容积为8000cm3的平底无盖的圆柱形容器，要求设计此容器消耗原材料最少，试写出这一优