最优化方法课程实验报告.docx

下载可编辑 下载可编辑 . .专业.整理. 项目一一维搜索算法 项目一 一维搜索算法 （ 一） [ 实验目的 ] 编写加步探索法 、对分法、Newton 法的程序 。 [ 实验准备 ] 1． 掌握一维收搜索中搜索区间的加步探索法的思想及迭代步骤 ； 掌握对分法的思想及迭代步骤 ； 掌握 Newton 法的思想及迭代步骤 。 [ 实验内容及步骤 ] 编程解决以下问题 ： 1． 用加步探索法确定一维最优化问题 min (t ) 的搜索区间 ，要求选取 t 0 t 0 0, h0 1, t 3 2 ． 2t 1 加步探索法算法的计算步骤 (1) 选取初始点 t 0 ： [ 0， )(或t 0 [ 0， t max ]) ， 计算 0 (t 0 ) ． 给出初始步长 h0 0 加步系数 1 ， 令 k (2) 比较目标函数值 ． 令 tk 1 则转 (4)。 0 。 tk hk ， 计算 k 1 (tk 1 ) ， 若 k 1 k ， 转(3)， 否 (3) 加大探索步长 ． 令 hk 1 hk ， 同时， 令 t tk , tk hk , t k 1 , k k 1， 转(2)。 (4) 反向探索 ． 若 k 0 ， 转换探索方向 ，令 hk t tk 1 ， 转（ 2）。 否则， 停止迭 代， 令 a min{ t，t k 1 }， b max{t，tk 1} 。 加步探索法算法的计算框图 程序清单 加步探索法算法程序见附录 1 实验结果 运行结果为 ： 用对分法求解  min  (t )  t( t  2) ， 已知初始单谷区间 [ a, b] [ 3, 5] ，要求按精度 0.3 ， 分别计算 ． 对分法迭代的计算步骤 ： 确定初始搜索区间 [ a, b] ， 要求 (a ) 0， ( b) 0 。 计算 [ a, b] 的中点 c 1 (a 2 b) ． 若 (c) 0 ， 则 a c ， 转(4) ； 若 (c) 0 ， 则 t * c ， 转 (5)； 若 (c) 0 ， 则 b c ，转 (4)． 若 | a b | ， 则t * (a 2 b) ， 转(5)； 否则转 (2)． 打印 t* ， 结束 对分法的计算框图 程序清单 对分法程序见附录 2 实验结果 运行结果为 ： 用 Newton 法求解  min  (t)  3t 2t 1 3 ， 已知初始单谷区间 [ a, b] [0, 1] ，要求精度 0.01 ． Newton 法的计算步骤 确定初始搜索区间 [a, b] ， 要求 ( a ) 0 ， ( b ) 0 选定 计算 t0 t t 0  (t0 ) / ( t0 ) 若 | t t0 | ，则 t0 t ， 转（ 3 ）； 否则转 （ 5）． 打印 t， (t ) ， 结束 ． Newton 法的计算框图 程序清单 Newton 法程序见附录 3 实验结果 运行结果为 ： 下载可编辑 下载可编辑 项目二 一维搜索算法 （ 二） [ 实验目的 ] 编写黄金分割法 、抛物线插值法的程序 。 [ 实验准备 ] 掌握黄金分割法的思想及迭代步骤 ； 掌握抛物线插值法的思想及迭代步骤 。 [ 实验内容及步骤 ] 编程解决以下问题 ： 用黄金分割法求解 min (t ) t(t 2) ， 已知初始单谷区间 [a, b] [ 3, 5] ，要求精度 0.001． 黄金分割法迭代步骤 ： 确定 (t) 的初始搜索区间 [ a，b] ． 计算 t 2 a 计算 t1 a 若 | t t | 0.382(b a) 0.618(b a) *， 则打印 t t1 *  t2 ， 结束 ； 否则转 (5)． 1 2 2 判别是否满足 1 ： 若满足 ， 则置 a t2，t 2 t1， 2 1 然后转 (3) ；否则 ， 置 b  t1，t1  t2， 1  2，t2  (b a)， 2  (t2 ) 然后转 (4) ． 黄金分割法的计算框图 ： .专业.整理. . .专业.整理. 程序清单 黄金分割法程序见附录 4 实验结果 运行结果为 ： 用抛物线插值法求解  min  f (x)  8x 3  2 x2  7 x 3 ， 已知初始单谷区间 [ a，b] [ 0，2]， 0.001． 抛物线插值法的计算步骤 ： （1 ） (t ) (t0 ) ， 所以相对 t0 来说 t 是好点 ， 故划掉区间 [t 0 , t2 ] ， 保留 [t1, t0 ] 为新区 间， 故置 t2 t0， (t2 ) (t0 ) ， t0 t， (t ) (t0 ) ， t1 保持