用二分法求下面方程在(-10,10)之间的根: 2x3-4x2+3x-6=0
- 为了把方程式简单化,可以把方程改写成 F=((2x-4)x+3)x-6;
- 二分法的思路为:先指定一个区间[x1,x2],如果函数在此区间是单调变化的,可以根据f(x1)和f(x2)是否同符号来确定方程f(x)=0在[x1,x2]区间是否只有一个实根。
- 若f(x1)和f(x2)之间不同号,则有一个确定的f(x0)=0,使得f(x)在区间范围内,有且只有一个实根。
- 若f(x1)和f(x2)之间同号,则无实根,则重新取x1,x2的值来继续判断。有实根后将[x1,x2]区间一分为二,来一半一半的判断在哪一半区间有值。
- 实现的大致步骤:
- 1:先输入一个x1,x2的值,通过函数来求出f1,f2的值,判断f1和f2是否异号,如果同号就继续输入x1,x2的值,不同号就跳出进入下一个程序。
- 2:现在的x1和x2已经异号了,然后取x0,使得x0=(x1+x2)/2,用函数求出f0的函数值。现在用目前的f0和f1(或者f2比较,哪个都可以,我就拿x1比较了)如果f0和f1异号,则说明x1和x0之间有根,则把x0赋值给x2,f0赋值给f2,(目的是把找出来新的f0赋值给右端点,因为跟只在x1和x0之间,之前的那个x2已经没有意义了,所以我把新的f0给f2,新的x0给x2.这样又是有左右端点了,x0又是x1和x2之间的值,就可以继续运算) 如果x1和x0之间没有根,则x0赋值给x1,f0赋值给f1;
- 最后判断f(x0)的绝对值是否小于某一个指定的值(例如10-5)
- 最后输出x0的值就是所求的近似根。
- 代码如下:
#include<stdio.h>
#include<math.h>
int main()
{
double x1,x2,f1,f2,x0,f0;
do
{
printf("Input x1,x2;");
scanf("%lf,%lf",&x1,&x2);
f1=x1*((2*x1-4)*x1+3)-6;
f2=x2*((2*x2-4)*x2+3)-6;
}while(f1*f2>0);
do
{
x0=(x1+x2)/2;
f0=x0*((2*x0-4)*x0+3)-6;
if(f0*f1<0){
x2=x0;
f2=f0;
}
else{
x1=x0;
f1=f0;
}
}while(fabs(f0)>=1e-5);
printf("x0=%5.2lf",x0);
return 0;
}