二分法函数仅对部分函数返回正确值的问题求助

核心问题：二分法初始条件不满足 #

你的二分法函数没有检查初始区间两端的函数值是否异号——这是二分法能正常工作的前提（根据中间值定理，只有当f(a)*f(b) < 0时，才能保证区间[a,b]内至少存在一个根）。

针对x_sin_exp函数，计算区间[-3,3]两端的函数值：

• f(-3) ≈ 0.5004（正数）
• f(3) ≈ 0.4996（正数）
  两端同号，二分法无法定位到根的位置，导致迭代过程中不断将左边界a向右移动，直到触发最大迭代次数，最终返回右边界3.0。

修复步骤 #

1. 给二分法函数添加初始条件校验

在函数开头加入检查逻辑，当初始区间不满足异号条件时，抛出异常或提示，避免无效迭代：

import numpy as np

def bisection(f: callable, a: float, b: float, tol=0.001, maxiter=100)-> tuple[float, int]:
    f_a = f(a)
    f_b = f(b)
    
    # 检查初始区间是否满足二分法前提
    if np.sign(f_a) == np.sign(f_b):
        raise ValueError("初始区间两端函数值同号，无法使用二分法")
    
    r = (a + b) / 2.0
    for nit in range(maxiter):
        f_r = f(r)
        if np.abs(f_r) < tol:
            return (r, nit+1)  # 迭代次数从1开始计数更合理
        elif np.sign(f_a) != np.sign(f_r):
            b = r
            f_b = f_r
        else:
            a = r
            f_a = f_r
        r = (a + b) / 2.0
    
    # 达到最大迭代次数，返回当前结果并提示
    print(f"警告：已达到最大迭代次数{maxiter}，可能未收敛")
    return (r, maxiter)

优化点说明：

• 提前计算f(a)和f(b)，避免循环内重复计算函数值，提升效率
• 使用np.sign替代直接相乘，避免极端值导致的数值溢出
• 调整迭代次数计数逻辑，从1开始更符合直觉
• 添加最大迭代次数警告

2. 为x_sin_exp选择正确的初始区间

先分析x_sin_exp(x)=0的根位置：

• 当x=1时，f(1)≈-0.177（负数）
• 当x=2时，f(2)≈0.447（正数）
  区间[1,2]满足两端异号，调用二分法即可得到正确结果：

result = bisection(x_sin_exp, 1, 2)
print(result)  # 输出示例：(1.29970703125, 10)

验证测试 #

• 测试x_2函数：bisection(x_2, -3, 3)仍能正确返回约(-2.236, 12)
• 测试x_sin_exp函数：使用[1,2]区间可收敛到正确根

内容的提问来源于stack exchange，提问作者Aledelmon