C语言基础--二分法

概念

二分法是一种常用的数值计算方法，用于求解单调函数的根或者最值。其原理是不断缩小求解区间的范围，直到找到满足精度要求的解。

注意！

二分法只能求解单调函数的根，因此在使用二分法求解方程时，需要先确定方程在求解区间内是否单调。同时，二分法的收敛速度比较慢，因此需要根据实际情况选择合适的精度要求和求解区间。

格式

1-库函数-...stdio....math...

2-定义函数-f（x）表达式

3-定义二分法-二分法基础格式

4-主内容-按所需要求填

举个例子

 

#include <stdio.h>
#include <math.h>

/*格式1-这两行是预处理指令，用于包含两个头文件：stdio.h 和 math.h。stdio.h 是C标准输入输出库，用于输入输出函数（如printf和scanf）。math.h 是数学函数库，提供了各种数学运算，例如绝对值函数fabs*/

double f(double x)
{
    return x * x - 2;

/*格式2-这是一个名为f的函数定义，它接受一个double类型的参数x，返回x的平方减去2。这个函数用于计算方程x^2 - 2 = 0的函数值*/
}

double bisection(double a, double b)
{
    double c;

/*以下皆为格式3-这是一个名为bisection的函数定义，它采用二分法求解方程的根。函数接受两个double类型的参数a和b，表示搜索根的初始区间。*/
 

  while (fabs(b - a) > 1e-8)

/*这是一个while循环，条件是区间[a, b]的长度大于1e-8（即0.00000001）。这个值是解的精度，表示我们希望找到的根与实际根的误差不超过这个值*/

{
        c = (a + b) / 2;

/*这行代码计算区间的中点c，这是二分搜索法的核心步骤。*/
        if (f(c) == 0) {
            return c;
        } else if (f(a) * f(c) < 0) {
            b = c;
        } else {
            a = c;
        }

/*这个if...else if...else结构用于检查中点c是不是精确的根（f(c) == 0）。如果不是，它将检查中点的函数值是否与区间一个端点的函数值符号不同（f(a) * f(c) < 0）。如果符号不同，这意味着根在区间[a, c]内，我们将b更新为c。如果符号相同，根在区间[c, b]内，我们更新a为c。*/

        c = (a + b) / 2;//最后更新c的值
    }
    
    return (a + b) / 2;

/*当while循环结束时（即区间的长度小于或等于指定精度），我们返回a和b的中点作为近似解。*/
}

int main()
{
    double a, b;
    scanf("%lf %lf", &a, &b);

/*格式4-这是程序的main函数，它首先定义了两个double类型的变量a和b。然后，scanf从标准输入读取两个double值，并将它们分别存储在a和b中。*/
 

   double x = bisection(a, b);

/*这行代码调用bisection函数，传入用户输入的a和b作为参数，并将返回的近似解存储在x变量中。*/
    printf("The solution of the equation is: %lf\n", x);
    return 0;
}

/*最后打印出x的解，main函数返回0，表示程序正常结束*/

当然你也可以给打印的x解设定保留几位小数点例如将"%lf"改为"%.2lf"表示保留两位小数。

完整代码

​
#include <stdio.h>

#include <math.h>



double f(double x)

{

    return x * x - 2; // 求解 x^2 - 2 = 0 的近似解

}



double bisection(double a, double b)

{

    double c;

    while (fabs(b - a) > 1e-8) {

        c = (a + b) / 2;

        if (f(c) == 0) {

            return c;

        } else if (f(a) * f(c) < 0) {

            b = c;

        } else {

            a = c;

        }

        c = (a + b) / 2;

    }

    

    return (a + b) / 2;

}



int main()

{

    double a, b;

    scanf("%lf %lf", &a, &b);

    double x = bisection(a, b);

    printf("The solution of the equation is: %lf\n", x);

    return 0;

}

​