问题描述:
- 求下面方程的一个根: f ( x ) = x 3 − 5 x 2 + 10 x − 80 = 0 f(x) = x^3-5x^2+10x-80 = 0 f(x)=x3−5x2+10x−80=0,若求出的根是啊,则要求 ∣ f ( a ) ∣ ≤ 1 / 1 0 6 |f(a)| \leq 1/10^6 ∣f(a)∣≤1/106
问题解析:
-
对 f ( x ) f(x) f(x)求导,得 f ‘ ( x ) = 3 x 2 − 10 x + 10 f^`(x) = 3x^2-10x+10 f‘(x)=3x2−10x+10。由一元二次方程求根公式知方程 f ‘ ( x ) = 0 f^`(x)=0 f‘(x)=0无解,因此 f ‘ ( x ) f^`(x) f‘(x)恒大于0.故 f ( x ) f(x) f(x)是单调递增的。易知 f ( x ) < 0 f(x) <0 f(x)<0且 f ( 100 ) > 0 f(100)>0 f(100)>0,所以区间[0,100]内必然有且只有一个根。由于 f ( x ) f(x) f(x)在[0,100]内是单调的,所以可以用二分的办法在区间[0,100]中寻找根。
-
C语言示例
#include<stdio.h>
#include<math.h>
double EPS = 1e-6;
double f(double x) { return x*x*x - 5*x*x + 10*x - 80; }
int main(){
double root, x1 = 0, x2 = 100, y;
int triedTimes = 1;
root = x1+(x2-x1)/2;
y = f(root);
while( fabs(y) > EPS) {
if( y>0 )
x2 = root;
else
x1 = root;
root = x1+(x2 - x1)/2;
y = f(root);
triedTimes++;
}
printf("%.8f\n", root);
printf("%d\n", triedTimes);
return 0;
}
-
运行示例
注:文中问题及代码参考 MOOC——《程序设计与算法》(北京大学 郭炜) -
算法系列目录:
- 一、递归形式的问题
- 二、程序或算法的时间复杂度
- 三、二分算法