现代设计方法与应用课件作者房亚东第六章节.pptVIP

* 主讲：房亚东 Email:fangyadong@ * 现代设计方法及其应用 西安工业大学机电学院 The Institute of Mechanical and Electrical Engineer, Xian Technological University 尚辅网 / 第六章 一维搜索方法 一维搜索的基本思想 一维搜索方法的Matlab程序及实例 3 1 2 确定单峰区间的方法——进退法 4 5 二次插值法 黄金分割法 ■ 内容 ：本章对一维搜索的基本思想和主要处理步骤进行论述，将介绍一维搜索的典型方法。首先介绍如何通过进退法实现单峰区间的确定，然后讨论了黄金分割法和二次插值法求极值点的算法描述及程序框图，最后通过三个实例说明了一维搜索方法的Matlab实现。 ■ 目的 ：了解一维搜索的基本思想和步骤，掌握进退法、黄金分割法和二次插值法的处理方法，能够通过Matlab环境实现一维搜索方法的处理。 优化问题涉及从当前点X(k)沿着S(k)方向，寻找最优步长因子α(k)，使本次搜索达到最小值，即f(X(k+1)) f(X(k))。上述问题可描述成以下优化问题： 6.1 一维搜索的基本思想 O X1 f(X) X2 f(X(k)) X(k) X(k+1)= X(k)+α(k) S(k) S(k) S(k)所在平面 minf(X(k)+α(k) S(k)) 图 二维问题转化一维问题时的几何关系 6.1 一维搜索的基本思想（续） 图 单峰函数与单峰区间 O x f(x) (a) O x (b) f(x) a b x* a x* b 一维搜索分为两个步骤：首先确定搜索区间，然后在基础上搜索最小点。 6.1 一维搜索的基本思想（续） 图 区间消去法示意图 O x f(x) (a) a b x1 x2 f1 f2 O x f(x) (c) a b x1 x2 f1 f2 O x f(x) (b) a b x1 x2 f1 f2 （1）如果f1 f2如图6-3(a)所示，则将单峰区间缩短为[x1 b]； （2）如果f1= f2如图6-3(b)所示，则将单峰区间缩短为[x1 x2]； （3）如果f1 f2如图6-3(c)所示，则将单峰区间缩短为[a x2]。 消去法的过程为：在区间中任取两点，计算和比较它们函数值的大小，消去大函数值一边的区间，剩下的区间中一定包含极小点。 6.2 确定单峰区间的方法——进退法 O x f(x) (a) x0 O x f(x) x0 x0+h x0+3h x0+h x0+3h x0-h (b) 前进 后退 图 进退示意图 进退法是一种通过比较函数值大小来确定单峰区间的方法，由单峰函数的性质可知，极小点左侧函数值严格下降，右侧严格上升。 6.2 确定单峰区间的方法——进退法（续） 图 进退法程序框图 黄金分割法又称为0.618法，它通过不断缩短区间长度来寻求一维函数f(x)的极小点。其基本思想是在已确定的单峰区间中，适当插入两点，利用区间消去法，寻求极小点。选点一般遵循以下原则： （1）等比收缩原则，区间每一次的缩短率不变； （2）对称取点原则，所插入两点在区间中位置对称。 6.3 黄金分割法 O x f(x) a x2 b x1 x3 图 黄金分割法示意图 6.3 黄金分割法（续） 图 黄金分割法程序框图 二次插值法又叫做近似抛物线法，属于利用拉格朗日多项式逼近的近似法。其基本思想是利用目标函数在不同三点的函数值，构成一个与原函数f(x)相近似的二次多项式 p(x)，以函数p(x)的极值xp*作为目标函数f(x)的近似极值点。经过多次迭代，逐步缩短区间，直至满足计算精度，最终逼近函数f(x)的极值点 6.4 二次插值法 6.4 二次插值法（续） 6.4 二次插值法（续） 二次函数插值法算法描述如下： STEP 1 给定单峰区间[x1 x3]和收敛精度ε，并设定x2=0.5( x1+x3)，进行算法初始化操作，使得f(x1) f(x2) f(x3)； STEP 2 构造过x1、x2、x3构造二次曲线p(x)，并按公式6-7计算二次曲线极小值xp*； STEP 3判断二次函数插入法是否满足精度要求： 若? x2- xp*?ε，则将x2和xp*代入f(x)中，并输出函数值较小的点的坐标及对应函数值，并停止求解； 若不满足精度要求，则转向STEP 4; STEP 4 由图6-7所示的x2和xp*的相对位置和函数值f2和f(xp*)大小6种情况，缩短单峰区间： 6.4 二次插值