C语言编程 二分法求方程的根

Python014

C语言编程 二分法求方程的根,第1张

这段代码是求解方程f(x)=0在区间[-10,10]上的根的数值解。

方法的思想就是:一直选取区间中间的数值,如果发现中间的函数值与一侧函数值,异号,那么说明解在这个更小的区间中,采用eps=1e-5作为区间的极限大小,通过迭代的方法求解这个方程的数值解。

所以了解了上述思想,那么else

if(f(a)*f(c)<0)

b=c

说明的是

f(a)和f(c)异号,那么使用b=(a+b)/2缩小迭代区间,继续迭代;同理else

a=c说明f(a)和f(c)同号,那么使用a(a+b)/2缩小迭代区间,继续迭代!

楼主你好~

提供以下代码。

Fedora10下gcc4.3.2编译通过,手动测试通过。。

#include <stdio.h>

#include <math.h>

#define DEFAULT_UPPER(10)

#define DEFAULT_LOWER(-10)

#define DEFAULT_E(0.00000001)

#define _MID(x,y)((x+y)/2)

#define _VALUE(x)(2*x*x*x-4*x*x+3*x-6)

double _e

int getRoot(double lower, double upper, double *result)

main()

{

double root

printf("Enter a deviation:")

scanf("%lf",&_e)

if(_e == 0.0)

_e = DEFAULT_E

if(getRoot(DEFAULT_LOWER, DEFAULT_UPPER, &root))

printf("Root:%2.8lf\n", root)

else

printf("Root:No Solution.\n")

}

int getRoot(double lower, double upper, double *result)

{

*result = _MID(lower,upper)

if(upper - lower <= _e)

return 1

if(_VALUE(lower)*_VALUE(*result) <= 0)

return getRoot(lower, *result, result)

else if(_VALUE(*result)*_VALUE(upper) <= 0)

return getRoot(*result, upper, result)

else

return 0

}

很简单的递归。

不过对于楼主的这道题,还要说一点,提示中的二分法求根,只能适用于在给定函数在给定区间中呈单调连续情况下才行得通哦~我也是按照给定的提示来写的上述代码。

如果连续函数在给定区间不单调,很有可能中值*下界值和中值*上界值都大于0,那么会跳出认为没有根,而事实上很有可能这个中值点靠近函数极点。

而真正用二分法求给定区间的思路是:

首先为函数求导,算出导函数的零点,然后再判断零点性质,最后将函数区间分为单调递增和单调递减间隔的形式,对每一段进行二分法求根。

请追问~