如何解决使用Scipy的Python代码中的除零错误？

1. 核心错误分析 #

你的代码存在两个关键问题：

• 导数调用错误：A(x).derivative是无效用法——A是普通Python数值函数，返回的是浮点数，没有derivative方法，这本身会触发AttributeError；你遇到的ZeroDivisionError是代码执行到x=0时的除法操作导致的，但导数调用错误是先存在的问题。
• x=0处的奇点：被积函数里的2/x在积分下限x=0时直接触发除以0错误，不过这个奇点是可去的，通过化简表达式就能消除。

2. 手动化简与修正代码 #

先手动推导A(z)的各阶导数并化简被积函数：

• A(z) = -a z²，一阶导数A’(z) = -2a z，二阶导数A''(z) = -2a
• 将导数代入原被积函数并化简，可完全消去分母的x，得到无奇点的表达式：

修正后的代码：

import numpy as np
from scipy.integrate import quad

a = 0.15

def integrand(x, z, B):
    # 化简后的被积函数，消除了x分母，避免除以0
    term_const = 18 * a - 2 * B**2
    term_x2 = 12 * a**2 - 2 * B**4
    inside_sqrt = 2 * (term_const + term_x2 * x**2)
    return np.sqrt(inside_sqrt)

def Phi(z, B):
    result, _ = quad(integrand, 0, z, args=(z, B))
    return result

phi_values = [Phi(z, 0) for z in range(101)]
print(phi_values[:5])  # 测试输出前5个值

3. 用SymPy自动推导化简（可选） #

如果你的函数更复杂，手动推导麻烦，可以用SymPy自动推导并化简表达式，再转成数值函数：

import numpy as np
from scipy.integrate import quad
import sympy as sp

# 定义符号变量
x, a_sym, B_sym = sp.symbols('x a B')
A_sym = -a_sym * x**2

# 计算各阶导数
A_prime = sp.diff(A_sym, x)
A_double_prime = sp.diff(A_sym, x, 2)

# 构建原被积函数表达式并化简
inside_expr = 3*x*A_double_prime - 3*x*A_prime**2 + 6*A_prime + 2*B_sym**4*x**3 + 2*B_sym**2*x
integrand_sym = sp.sqrt(-(2/x)*inside_expr)
simplified_integrand = sp.simplify(integrand_sym)

# 转换成numpy可调用的数值函数
integrand_num = sp.lambdify((x, a_sym, B_sym), simplified_integrand, 'numpy')

a = 0.15

def integrand(x, z, B):
    return integrand_num(x, a, B)

def Phi(z, B):
    result, _ = quad(integrand, 0, z, args=(z, B))
    return result

phi_values = [Phi(z, 0) for z in range(101)]
print(phi_values[:5])

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