04第四章_无约束问题直接法01.ppt

04第四章_无约束问题直接法01

 (2)  这表示 点走得太远,应缩回一些.若以 表示压缩因子,则有      (4.1)  常取为0.5. (3) 这时应压缩更多一些,将新点压缩至 至 之间,有     (4.2)  注意,将式(4.1)中的 代之以 ,即可得式(4.2). ⑦若 ,则转⑤;否则,置 转③. ⑧若 ,则输出 ,停止计算;否则,置 ,转②. 19 22 18 注* 22点模移动失败 重置18点为探测初始点 注意:算法中的模式为 。当由3产生时,模式 ;但当由6产生时,模式才为 ,这是加速 既为 模式。 例4.2 用步长加速法求解 已知 ,收缩系数 要求一直迭代到步长向量满足 为止。 。 解 令 ,计算 开始 Ⅰ型探测 ,√ 记 ,接着计算 ,× ,√ ,× 得 及 。因为 于是进行第一次模式移动 , ,并计算 。开始 Ⅱ型探测 ,× ,√ 记 ,接着计算 ,× ,√ 得 及 。因为 第二次模式移动 ,于是进行 ,并计算 。开始 Ⅱ型探测 × √ × 得 及 。因为 可以进行第三次模式移动 ,于是又 ,并计算 。又开始 Ⅱ型探测 √ 记 ,接着计算 √ 得 及 。因为 上次模式移动作废。 ,所以 重令 ,则 ,并又开始 Ⅰ型探测 × × × 探测失败,需缩小步长,新的步长向量 计算终止。 × , 以上过程的路径如下图所示。 * 第四章 无约束问题 ——直接法 第二节 步长加速法(轴向、模式搜索法) 第三节 方向加速法(Powell) 第一节 单纯形替换法 第一节 单纯形替换法 1.单纯形替换法: Spendley(斯彭德莱) 、Hext(赫克斯特) 和Himsworth(希姆斯沃思 )于1962年提出; Nelder(尼尔德 ) 和Mead(米德)1965年改进 问 题: 3.算法思想 4. 单纯形概念 (2)单纯形的定义 现以求二元函数的极小点为例,说明单纯形法形成原理. 设二元函数 在 平面上取不在同一条直线上的三个点 , 和 ,并以它们为顶点构成一单纯形——三角形.算出各顶点的函数值 , , ,比较其大小,现假定比较后有 5.如何构造单纯形? 这说明点 最差,点 最好,点 次差.为了寻找极小点,一般来说应向最差点的反对称方向进行搜索.以 记为 的中点(如图4.1所示),在 的延长线上取点 ,使 称为 关于 的反射点. 算出 的函数值 ,可能出现以下几种情形: (1) 这说明搜索方向正确,可进一步扩大效果,继续沿 向前搜索,也就是向前扩张.这时取 其中 为扩张因子,一般取 . 如果 ,说明扩张有利,就可以点 代替点 构成新的单纯形 .如果 ,说明扩张不利,舍去 ,仍以 代替 构成新的单纯形   图4.1 返 回 (4) 若 方向上所有点的函数值 都大于 ,则不能沿此方向搜索.这时,可以以 为中心进行缩边,若使顶点 和 向 移近一半距离(如图4.2所示),得新单纯形 .以此单纯形为基础再进行寻优. 图 4.2 返 回 以上说明,不管哪种情况,我们都可以得到一个新的单纯形,其中至少有一顶点的函数值比原单纯形为小.如此继续,直至满足收敛终止准则. 在

文档评论(0)

1亿VIP精品文档

相关文档