二分法与黄金分割法求函数方程最小值

最新推荐文章于 2024-03-24 23:31:13 发布

原创

最新推荐文章于 2024-03-24 23:31:13 发布 · 9.7k 阅读

收录于

当前文章被以下社区和专栏收录：

这篇博客介绍了使用C++编程实现二分法和黄金分割法来寻找函数方程的最小值。通过实验数据，如区间[0.3, 1]和误差ε=1e-4，展示了两种方法的计算结果。二分法得到的最小值点在0.637009附近，黄金分割法得到的最小值点在0.636939附近，两者结果在允许误差范围内基本一致，证明了算法的准确性。" 123237365,12436646,Vue3组件与动画实战,"['Vue.js', '动画', '组件开发']

二分法与黄金分割法求函数方程最小值

程序设计语言：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 : r;   

 t = fun((l + r) / 2.0); 
 }  

 re