二分法与黄金分割法求函数方程最小值
程序设计语言:C++
输入:线性搜索模型(目标函数系数,搜索区间,误差限等)
输出:最优解及对应目标函数值
实验数据
区间[0.3,1],误差ε=1e-4
结果
二分法最后结果为:导数在0.637009处取零点,函数在区间【0.3,1】取得最小值11.252754.
(二分法结果图)
黄金分割法最后结果为:输入初始值0.3,导数在0.636939处取零点,函数在区间【0.46,0.94】取得最小值11.252756.
(黄金分割法结果图)
两种方法所求无论是极小值点还是最小值点结果在所限制误差内基本一致,表明结果正确。
代码:
二分法:
#include <stdio.h>
#include <math.h>
double fun(double x)
{
return -2*(x-0.3)/ (((x-0.3)*(x-0.3) + 0.01)*((x-0.3)*(x-0.3) + 0.01) )- 2*(x-0.9)/ (((x-0.9)*(x-0.9) + 0.04)*((x-0.9)*(x-0.9) + 0.04));
}
double solve(double l, double r)
{
double t = fun((l + r) / 2.0);
while(fabs(t) >= 1e-4)
{
l = t < 0 ? (l + r) / 2.0 : l;
r = t > 0 ? (l + r) / 2.0