机械优化设计——鲍威尔法.docVIP

机械优化设计 鲍 威 尔 法 班级：0841001 成员：张波2010213217 张建2010213214 潘阳瑞2010213227 2013年6月 鲍威尔（Powell）法是直接利用函数值来构造共轭方向的一种方法。 基本思想：在不用导数的前提下，在迭代中逐次构造G 的共轭方向。 一．基本算法：（二维情况描述鲍威尔的基本算法） 1）任选一初始点x0，再选两个线性无关的向量，如坐标轴单位向量e1=[1,0]T和e2=[0,1]T作为初始搜索方向。 2）从x0出发，顺次沿、作一维搜索，得 、 点，两点连线得一新 方向 用 代替e1，形成两个线性无关向量,e1，作为下一轮迭代的搜索方向。再从出发，沿作一维搜索得点，作为下一轮迭代的初始点。 3）从出发，顺次沿、作一维搜索，得到点、 ，两点连线得一新方向: 。 沿作一维搜索得点，即是二维问题的极小点。 把二维情况的基本算法扩展到n维，则鲍威尔基本算法的要点是： 在每一轮迭代中总有一个始点（第一轮的始点是任选的初始点）和n个线性独立的搜索方向。从始点出发顺次沿n个方向作一维搜索得一终点，由始点和终点决定了一个新的搜索方向。 用这个方向替换原来n个方向中的一个，于是形成新的搜索方向组。替换的原则是去掉原方向组的第一个方向而将新方向排在原方向的最后。此外规定，从这一轮的搜索终点出发沿新的搜索方向作一维搜索而得到的极小点，作为下一轮迭代的始点。这样就形成算法的循环。 图1.二维情况下的鲍威尔法 二．改进算法 在鲍威尔基本算法中，每一轮迭代都用连结始点和终点所产生出的搜索方向去替换原向量组中的第一个向量，而不管它的“好坏”，这是产生向量组线性相关的原因所在。 在改进的算法中首先判断原向量组是否需要替换。如果需要替换，还要进一步判断原向量组中哪个向量最坏，然后再用新产生的向量替换这个最坏的向量，以保证逐次生成共轭方向。 为此，要解决两个关键问题： （1）是否较好？是否应该进入新的方向组？即方向组是否进行更新？ （2）如果应该更新方向组，不一定替换方向，而是有选择地替换某一方向。 改进算法的具体步骤： （1）给定初始点（记作），选取初始方向组，它由n个线性无关的向量，,...,所组成，置。 （2）从出发，顺次沿，，...，作一维搜索得，，...，。接着以为起点，沿方向 移动一个的距离，得到 、、分别称为一轮迭代的始点、终点和反射点。始点、终点和反射点所对应的函数值分别为： 同时计算各中间点的函数值，并记为：（i=1,2,...,n) 因此有， 计算n个函数值之差，，...，。 记作： （i=1,2,...,n) 其中最大者记作： (3)根据是否满足判断条件和，来确定是否要对原方向组进行替换。 若不满足判别条件，则下轮迭代仍使用原方向组，并以、中函数值小者作为下轮迭代的始点。 若满足上述判别条件，则下轮迭代应对原方向组进行替换，将补充到原方向组的最后位置，而除掉。即新方向组为作为下轮迭代的搜索方向。下轮迭代的始点取为沿方向进行一维搜索的极小点。 判断是否满足收敛准则。若满足则取为极小点，否则应置，返回2，继续进行下一轮迭代。 这样重复迭代的结果，后面加进去的向量都彼此对G共轭，经 n轮迭代即可得到一个由n个共轭方向所组成的方向组。对于二次函数，最多n次就可找到极小点，而对一般函数，往往要超过n次才能找到极小点（这里“n”表示设计空间的维数）。 图2.鲍威尔法的程序框图 题目：用改进的鲍威尔法求目标函数的最优解。已知初始点[1,1]T，迭代精度。 解：初始点，初始搜索方向，。初始点处的函数值。 第一轮迭代： 1）沿方向进行一维搜索，得 ： 为一维搜索最佳步长，应满足极值必要条件： 所以算出一维搜索最佳步长 得 从而算出点处的函数值及沿走步后函数值的增量 再沿方向进行一维搜索，得 为一维搜索最佳步长，应满足极值必要条件： 所以算出一维搜索最佳步长 得 从而算出第一轮终点处的函数值及沿走步后函数值的增量 取沿，