使用SymPy的fraction方法提取嵌套分式分子分母遇异常结果求助

问题原因 #

你定义的表达式a=1/(1+1/x)在SymPy中是未化简的原始形式，fraction()方法只会直接读取当前表达式字面意义上的分子和分母，不会主动进行通分、化简操作，所以才会返回(1, 1+1/x)，而不是你期望的(x, x+1)。

解决方法 #

在调用fraction()之前，先对表达式进行有理式化简，推荐使用ratsimp()方法（专门用于有理表达式的化简，比通用的simplify()更高效），或者用simplify()也能达到效果。

修改后的代码 #

from sympy import *
x = symbols('x')
a = 1/(1+1/x)
# 先对表达式进行有理化化简
simplified_a = ratsimp(a)
# 再提取分子分母
num, den = fraction(simplified_a)
print((num, den))

输出结果 #

(x, x + 1)

补充说明 #

• ratsimp()会自动处理有理表达式的通分、约分，把表达式转化为最简分式形式，这时候fraction()就能正确提取到你想要的分子和分母了。
• 如果用simplify()替换ratsimp()，代码同样可以正常运行，结果一致，只是simplify()是通用化简函数，对于复杂表达式可能会做更多额外处理，而ratsimp()更专注于有理式的化简，效率更高。

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