机械优化设计-一维搜索方法

以以 f x x 2 x 为例为例 黄金分割 include define e 0 001 收敛精度 main float a b x1 x2 y1 y2 float minx miny x1 3 x2 5 x1 x2 初始化 a x1 b x2 a1 a2 初始化 do x1 a 0 382 b a x2 a 0 618 b a y1 x1 x1 2 x1 y2 x2 x2 2 x2 if y1 y2 a x1 else b x2 while b a e 两个数值小于收敛精度 e 时终止 minx x1 x2 2 取两数值平均值 miny minx minx 2 minx printf The minx is f and the miny is f n minx miny 极小值 二次插值法 include include int main void double a1 a2 a3 ap y1 y2 y3 yp c1 c2 m double j 3 int i h 1 取正向搜索 h 值为负则为反向搜索 void finding double a 3 一维搜索确定区间 finding j 返回所得的区间值 a1 a2 a3 a1 j 0 以数组 j 3 为过渡 衔接搜索区间的值并赋值给主函数的 a1 a2 a3 a2 j 1 实际上是不会这么用的 此处是方便衔接一维搜索函数与主 函数 以后可类似 a3 j 2 m 0 001 输入 m 值 以确定精确度 double f double x 方程函数说明 y1 f a1 初始化 y2 f a2 y3 f a3 for i 1 1 1 i c1 y3 y1 a3 a1 c2 y2 y1 a2 a1 c1 a2 a3 ap 0 5 a1 a3 c1 c2 yp f ap if fabs y2 yp y2 0 if y2 yp a1 a2 y1 y2 a2 ap y2 yp else a3 ap y3 yp else if y2 yp a3 a2 y3 y2 a2 ap y2 yp else a1 ap y1 yp double x y 用 x y 储存所得极值点 if y2y 0 将逆方向寻找变向 h h a 2 a 0 y 2 y 0 do a 0 a 1 a 1 a 2 y 0 y 1 y 1 y 2 a 2 a 1 h y 2 f a 2 h 2 h while y 2 1 i a 2 a 1 h y 2 f a 2 if y 2 y 1 break 寻找成功 打印输出 h 2 h a 0 a 1 y 0 y 1 a 1 a 2 y 1 y 2 return 外推法 include float f float t float z z t t 2 t return z void main float h0 1 h h0 a1 0 a2 h y1 y2 a3 y3 y1 f a1 y2 f a2 if y2 y1 h h a3 a1 y3 y1 a1 a2 y1 y2 a2 a3 y2 y3 a3 a2 h y3 f a3 while y3 y2 h 2 h a1 a2 y1 y2 a2