C语言用二分法实现函数在区间范围内的根

用二分法求下面方程在(-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;
}

  • 9
    点赞
  • 65
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
这里提供一种用C语言二分法求解方程根的代码思路: 1.定义函数f(x),用来表示要求根的方程,根据题目中的式子,即为f(x) = 2*x*x*x - 4*x*x + 3*x - 6。 2.定义二分法函数,传入参数为左端点left、右端点right、所需精度precision,返回值为方程的一个根root。 3.在二分法函数内,利用while循环进行二分查找: (1)定义变量mid为left和right的中点,计算f(mid)的值; (2)如果f(mid)的绝对值小于等于所需精度precision,即可认为找到了一个根,将其赋值给root并返回; (3)如果f(mid)小于0,则根在[mid, right]区间内,更新left为mid; (4)如果f(mid)大于0,则根在[left, mid]区间内,更新right为mid。 4.在主函数内,调用二分法函数求出[-10,10]范围内的一个根,并输出结果。 具体代码如下: #include <stdio.h> double f(double x); //定义方程的函数 double find_root(double left, double right, double precision); //定义二分法函数 int main() { double root = find_root(-10, 10, 0.0001); //在[-10,10]范围内查找一个根,精度为0.0001 printf("The root of the equation is: %lf\n", root); return 0; } double f(double x) { return 2*x*x*x - 4*x*x + 3*x - 6; } double find_root(double left, double right, double precision) { double mid, root; while (1) //不断进行二分查找 { mid = (left + right) / 2; if (f(mid) <= precision && f(mid) >= -precision) //找到一个根 { root = mid; break; } else if (f(mid) < 0) //根在[mid,right]区间内 { left = mid; } else //根在[left,mid]区间内 { right = mid; } } return root; }
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值