二分法（Bisection Method）是一种求解方程根的数值方法，适用于连续函数在给定区间内存在单根的情况。其核心思想是通过不断缩小区间范围逼近方程的根。以下是详细步骤：

算法原理

1. 前提条件

   • 函数 ( f(x) ) 在闭区间 ([a, b]) 上连续。
   • ( f(a) ) 与 ( f(b) ) 异号（即 ( f(a) \cdot f(b) < 0 )），由介值定理可知区间内至少有一个根。

2. 迭代过程

   • 步骤1：计算中点 ( c = \frac{a + b}{2} )。
   • 步骤2：计算函数值 ( f(c) )。
   • 步骤3：判断根的位置：
     • 若 ( f(c) = 0 )，则 ( c ) 为根，停止迭代。
     • 若 ( f(a) \cdot f(c) < 0 )，则根在 ([a, c]) 内，更新 ( b = c )。
     • 若 ( f(b) \cdot f(c) < 0 )，则根在 ([c, b]) 内，更新 ( a = c )。
   • 步骤4：重复上述步骤，直到区间长度小于预设精度 ( \epsilon )。

3. 终止条件

   • 区间长度 ( |b - a| < \epsilon )（(\epsilon) 为允许误差，如 ( 10^{-6} )）。
   • 此时近似根为最终区间的中点 ( c )。

示例：求解 ( f(x) = x^3 - x - 1 = 0 )

• 初始区间：( [1, 2] )（因 ( f(1) = -1 < 0 )，( f(2) = 5 > 0 )）。
• 精度要求：( \epsilon = 0.001 )。

结果：近似根为 ( x \approx 1.3247 )（误差 ( < 0.001 )）。

Python 代码实现

def bisection(f, a, b, tol=1e-6, max_iter=100):
    if f(a) * f(b) >= 0:
        raise ValueError("区间两端点函数值需异号！")
    for i in range(max_iter):
        c = (a + b) / 2
        if abs(b - a) < tol:  # 达到精度要求
            return c
        if f(c) == 0:         # 恰好找到根
            return c
        if f(a) * f(c) < 0:   # 根在左半区间
            b = c
        else:                 # 根在右半区间
            a = c
    return (a + b) / 2         # 返回最终中点

# 示例函数
def f(x):
    return x**3 - x - 1

# 调用二分法
root = bisection(f, 1, 2)
print(f"方程的近似根为: {root:.6f}")

优缺点

• 优点：简单可靠，收敛性有保证（线性收敛）。
• 缺点：
  • 仅适用于单根且函数值异号的区间。
  • 收敛速度较慢（相比牛顿法）。
  • 无法直接处理重根或复数根。

应用场景

• 工程计算中需要快速估计根的初值。
• 函数连续但导数复杂或不可求的情况。
• 对精度要求不高的实时计算。