使用SymPy实现三角代换积分的问题求助

问题背景 #

学习三角代换积分时，尝试用SymPy复现该方法，但遇到两个具体问题：

1. 积分操作后出现额外微分项
2. 回代变量后不知如何应用边界条件

原尝试代码如下：

import sympy as sy
import sympy.diffgeom as dg
sy.init_printing(use_latex=True)
x, r, t = sy.symbols('x r t')
dx = dg.Differential(x)
dt = dg.Differential(t)
y = sy.sqrt(r ** 2 - x ** 2) * dx
y = y.subs({x: r * sy.sin(t), dx: r * sy.cos(t) * dt})
y = sy.simplify(y)
#here's where the code breaks down
y = sy.integrate(y, t)
y = y.subs(t, sy.asin(x / r))
sy.pprint(y)

问题分析与解决 #

1. 积分后出现额外微分的解决 #

问题根源是误用了sympy.diffgeom.Differential对象，SymPy的普通积分功能不需要处理微分形式，直接对表达式进行代换即可。三角代换的核心是将x = r sin(t)、dx = r cos(t) dt转换为被积表达式的替换，无需构造微分对象。

修改后的核心代码片段：

import sympy as sy
sy.init_printing(use_latex=True)
x, r, t, C = sy.symbols('x r t C')

# 原被积表达式
expr = sy.sqrt(r ** 2 - x ** 2)
# 执行三角代换：替换x并乘上dx对应的表达式
substituted_expr = expr.subs(x, r * sy.sin(t)) * r * sy.cos(t)
# 化简代换后的表达式
simplified_expr = sy.simplify(substituted_expr)

# 对t积分，得到含常数项C的原函数
integrated_result = sy.integrate(simplified_expr, t) + C

2. 回代变量与应用边界条件 #

回代变量t = asin(x/r)后，可通过两种方式应用边界条件：

方式一：计算定积分

直接代入积分上下限，用上限函数值减去下限函数值：

# 回代变量t
original_function = integrated_result.subs(t, sy.asin(x / r))
simplified_original = sy.simplify(original_function)

# 示例：计算x从0到r的定积分
definite_result = simplified_original.subs(x, r) - simplified_original.subs(x, 0)
sy.pprint(sy.simplify(definite_result))

方式二：根据初始条件确定常数C

若已知某点的函数值，代入方程求解C，得到特解：

# 示例：初始条件x=0时，函数值为0
equation = sy.Eq(simplified_original.subs(x, 0), 0)
C_solution = sy.solve(equation, C)[0]

# 代入常数C得到特解
specific_solution = simplified_original.subs(C, C_solution)
sy.pprint(specific_solution)

完整可运行代码 #

import sympy as sy
sy.init_printing(use_latex=True)
x, r, t, C = sy.symbols('x r t C')

# 1. 三角代换与积分
expr = sy.sqrt(r ** 2 - x ** 2)
substituted_expr = expr.subs(x, r * sy.sin(t)) * r * sy.cos(t)
simplified_expr = sy.simplify(substituted_expr)
integrated_result = sy.integrate(simplified_expr, t) + C

# 2. 回代变量
original_function = integrated_result.subs(t, sy.asin(x / r))
simplified_original = sy.simplify(original_function)
print("不定积分结果：")
sy.pprint(simplified_original)

# 3. 计算定积分（示例：0到r）
definite_result = simplified_original.subs(x, r) - simplified_original.subs(x, 0)
print("\n定积分（0到r）结果：")
sy.pprint(sy.simplify(definite_result))

# 4. 应用初始条件求特解（示例：x=0时函数值为0）
equation = sy.Eq(simplified_original.subs(x, 0), 0)
C_solution = sy.solve(equation, C)[0]
specific_solution = simplified_original.subs(C, C_solution)
print("\n满足初始条件的特解：")
sy.pprint(specific_solution)

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