如何在SymPy积分中指定特定哑变量？报错问题求解

问题场景 #

运行自定义变量的SymPy积分代码时，始终报错ValueError: specify dummy variables for (...)，代码如下：

import numpy as np
import seaborn
import math
import pandas as pd
import matplotlib as plt
import scipy.integrate as integrate
import scipy.special as special
from sympy import init_session

init_session()
%matplotlib inline

dP, dT = symbols('dP dT', real = True, positive = True)
dr = symbols('dr', real = True, positive = True)
T_9, T = symbols('T_9 T', real = True, positive = True)
rho, mu, m_u, k, P, gamma = symbols('rho mu m_u k P gamma', real = True, positive = True)
r, R_solar, R = symbols('r R_solar R', real = True, positive = True)
G, M_solar, M = symbols('G M_solar M', real = True, positive = True)
a, c, kappa = symbols('a c kappa', real = True, positive = True)
F = symbols('F', real = True, positive = True)

cc = (T / (rho * k * T / (mu * m_u)) * (1 - (1 / gamma)) * ((-rho * G * M / r**2)))

r_int = integrate(((cc), (0, r, R)))

错误原因 #

核心问题是积分调用格式错误：

• SymPy的integrate函数要求明确指定「积分变量（哑变量）」，语法为integrate(被积函数, (积分变量, 下限, 上限))
• 你写的integrate(((cc), (0, r, R)))既没指定积分变量，还错误嵌套了括号，导致SymPy无法识别积分对象。

解决步骤 #

1. 先化简被积函数（可选但推荐）

cc表达式存在可约分项，先化简能大幅简化计算：

from sympy import simplify
cc_simplified = simplify(cc)
# 化简结果：-(1 - 1/gamma)*G*M*mu*m_u/(k*r**2)

2. 正确调用积分函数

直接明确指定积分变量为r，修正积分限格式：

r_int = integrate(cc_simplified, (r, 0, R))
# 积分结果：(G*M*mu*m_u*(1 - 1/gamma))/(k*R)

3. 用哑变量避免变量冲突（可选）

如果担心当前r变量与积分变量产生冲突，可以单独定义哑变量：

r_dummy = symbols('r_dummy', real=True, positive=True)
cc_dummy = cc.subs(r, r_dummy)
r_int = integrate(cc_dummy, (r_dummy, 0, R))

修正后的完整代码 #

import numpy as np
import seaborn
import math
import pandas as pd
import matplotlib.pyplot as plt  # 修正原代码的导入错误
import scipy.integrate as integrate
import scipy.special as special
from sympy import init_session, symbols, simplify, display

init_session()
%matplotlib inline

dP, dT = symbols('dP dT', real = True, positive = True)
dr = symbols('dr', real = True, positive = True)
T_9, T = symbols('T_9 T', real = True, positive = True)
rho, mu, m_u, k, P, gamma = symbols('rho mu m_u k P gamma', real = True, positive = True)
r, R_solar, R = symbols('r R_solar R', real = True, positive = True)
G, M_solar, M = symbols('G M_solar M', real = True, positive = True)
a, c, kappa = symbols('a c kappa', real = True, positive = True)
F = symbols('F', real = True, positive = True)

cc = (T / (rho * k * T / (mu * m_u)) * (1 - (1 / gamma)) * ((-rho * G * M / r**2)))
cc_simplified = simplify(cc)

# 正确的积分调用
r_int = integrate(cc_simplified, (r, 0, R))
display(r_int)  # 显示格式化的数学结果

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