Python SymPy库积分运算两类异常问题的解决方法咨询

问题1：运行时报错expr not of form a*x**b: %s #

• 错误原因：代码中使用了Python原生的浮点数指数3/2，SymPy进行符号积分时无法识别浮点数的精确代数结构，触发规则匹配错误。
• 解决方法：将浮点数指数替换为SymPy的精确有理数类型Rational
• 修复后代码：

from sympy import *
x = Symbol('x')
f = 1 / (25 + x**2)**Rational(3,2)
print(integrate(f, x))

• 运行后输出正确结果：x/(25*sqrt(x**2 + 25))

问题2：积分结果为复数对数形式，和预期的反正切形式不符 #

• 错误原因：定义符号a时未指定域约束，SymPy默认a属于复数域，因此返回复数域下的通用积分结果，该结果和实数域下的反正切结果等价但表现形式不同。
• 两种可选解决方法：
  1. 定义符号时指定a为实数类型，SymPy会直接返回实数域下的反正切结果
     修复代码：

  from sympy import *
  x, a = symbols('x a', real=True)
  f = 1 / (a**2 + x**2)
  print(integrate(f, x))
  

  2. 对已得到的复数对数结果调用rewrite('atan')方法，手动转换为反正切形式
     示例代码：

  from sympy import *
  x, a = symbols('x a')
  f = 1 / (a**2 + x**2)
  res = integrate(f, x)
  print(res.rewrite('atan'))
  

两种方法均可得到预期结果atan(x/a)/a。

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