割线法实现遇NoneType错误：无法比较None与浮点数

问题根源 #

1. solve函数未返回计算结果，默认返回None，导致f_a、f_b、f_c均为NoneType，无法与浮点数tolerance进行大小比较，直接引发循环执行错误。
2. 迭代公式错误：代码中c=(a + b)/2是二分法的中点计算，并非割线法的迭代公式，会导致算法逻辑偏离预期。
3. 循环条件未考虑函数值为负的情况，且初始值设置不合理，可能导致循环提前终止或无法正确执行。

修复步骤 #

1. 修复solve函数 #

为solve添加返回语句，同时改用更安全的变量传递方式（直接将x映射为输入值，避免字符串替换的潜在问题）：

def solve(function, n):
    return eval(function, {"x": n})

2. 修正割线法迭代公式 #

割线法的迭代公式应为：
c = b - f_b * (b - a) / (f_b - f_a)
替换原代码中的c=(a + b)/2，确保算法逻辑正确。

3. 调整循环逻辑 #

• 初始时先计算f_a和f_b，避免用假的初始值触发循环；
• 循环条件改为判断函数值的绝对值是否大于容差（函数值可能为负，直接比较大小会出错）；
• 移除不必要的f_c循环判断，割线法每次迭代只需要前两个点的函数值即可推进。

完整修复代码 #

def solve(function, n):
    return eval(function, {"x": n})

def SecantMethod(a, b, tolerance, function):
    a = float(a)
    b = float(b)
    tolerance = float(tolerance)
    iteration = 0
    
    # 先计算初始两点的函数值
    f_a = solve(function, a)
    f_b = solve(function, b)
    
    # 循环至函数值绝对值小于容差
    while abs(f_b) > tolerance:
        # 割线法迭代公式
        c = b - f_b * (b - a) / (f_b - f_a)
        f_c = solve(function, c)
        
        # 更新迭代点与对应函数值
        a, b = b, c
        f_a, f_b = f_b, f_c
        iteration += 1
    
    text = """stopped at iteration = {}
  a = {}
  b = {}
  f(a) = {}
  f(b) = {}
""" 
    print(text.format(iteration, a, b, f_a, f_b))

额外说明 #

• 使用eval存在安全风险，如果是处理不可信的函数输入，建议改用sympy等专业数学库来解析表达式；
• 可以添加迭代次数上限，避免因函数不收敛导致的死循环。

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