Python数值积分报错：function与Symbol不支持/运算问题求助

报错原因 #

• 你混用了SymPy符号变量和Python普通函数：r是SymPy定义的符号，但b、f是lambda函数（函数对象），当你写b/r时，相当于把函数和符号变量做除法，这两种类型不兼容，直接触发TypeError。
• 另外，scipy.integrate.quad要求被积函数是数值函数（接收数值输入，返回数值结果），但你传入的Int是SymPy的符号表达式，quad无法直接处理。

修正思路 #

放弃SymPy的符号变量，全程用数值运算适配quad的要求——因为你明确提到积分解析不可解，数值积分才是正确方向。

修正后的代码 #

import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import quad

# 定义常数
s = 1
Mu = 0.5
lamb = 0.6

# 定义b函数：接收数值r，返回数值结果
def b(r):
    return ((r - s) * lamb + (s**2)*r) / (s*r)

# 定义f函数：接收数值r，返回数值结果
def f(r):
    return (1 - lamb/(s*r)) ** ((s**2 + lamb)/lamb)

# 定义被积函数：接收数值r，返回数值结果
def integrand(r):
    denominator = r * np.sqrt( (1 - b(r)/r) * ( (r**2)/(Mu**2) * (1/f(r)) - 1 ) )
    return 1 / denominator

# 测试单点积分：下限4，上限无穷大
result, error = quad(integrand, 4, np.inf)
print("积分结果：", result)
print("误差估计：", error)

代码解释 #

• 把lambda函数改成普通def函数，更直观，也方便调试。
• 用np.sqrt()替代SymPy的smp.Rational(1/2)次方，数值运算场景下numpy的sqrt更高效。
• 被积函数integrand完全是数值运算，接收数值r、返回数值，符合quad的要求。
• quad返回两个值：积分结果和误差估计，需要同时接收。

后续绘图思路 #

如果要绘制函数F(x) = ∫(x到∞) integrand(r) dr，可以参考以下代码：

# 生成x的取值数组（避开被积函数无意义的区间，比如分母为0的情况）
x_values = np.linspace(0.1, 10, 50)
# 计算每个x对应的积分结果
y_values = [quad(integrand, x, np.inf)[0] for x in x_values]

# 绘图
plt.figure(figsize=(8, 5))
plt.plot(x_values, y_values, label=r'$F(x) = \int_{x}^{\infty} \text{integrand}(r) dr$')
plt.xlabel('x')
plt.ylabel('F(x)')
plt.title('积分定义的函数曲线')
plt.legend()
plt.grid(True)
plt.show()

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