Sympy 可以通过收集多个项来简化有理表达式吗？

Question

给定一个有理表达式 E 如下所示，我希望使用 Sympy 将其简化为看起来像 F 的东西（在 [=44 的第二个块中定义） =] 下面的代码):

import sympy as sp

a, b, c, d, n, t, A, B, C = sp.symbols('a, b, c, d, n, t, A, B, C', real = True)

E = n/(c-b) * ( B - (c-b)/(c-a)*A - (b-a)/(c-a)*B ) * (c-t)/(c-b) + n/(c-b) * ( (d-c)/(d-b)*B + (c-b)/(d-b)*C - B ) * (t-b)/(c-b)

print(sp.pretty( E ))
print(sp.pretty( E.simplify() ))

这会打印

           ⎛     B⋅(-c + d)   C⋅(-b + c)⎞             ⎛  A⋅(-b + c)   B⋅(-a + b)    ⎞
n⋅(-b + t)⋅⎜-B + ────────── + ──────────⎟   n⋅(c - t)⋅⎜- ────────── - ────────── + B⎟
           ⎝       -b + d       -b + d  ⎠             ⎝    -a + c       -a + c      ⎠
───────────────────────────────────────── + ─────────────────────────────────────────
                        2                                           2                
                (-b + c)                                    (-b + c)
                
                
-n⋅((a - c)⋅(b - t)⋅(-B⋅(b - d) + B⋅(c - d) + C⋅(b - c)) + (b - d)⋅(c - t)⋅(A⋅(b - c) + B⋅(a - b) - B⋅(a - c))) 
────────────────────────────────────────────────────────────────────────────────────────────────────────────────
                                                           2                                                    
                                            (a - c)⋅(b - c) ⋅(b - d)

但是，可以手动进一步简化表达式，我将其标记为 F:

F = n/(c-a) * (B - A) * (c-t)/(c-b) + n/(d-b) * (C - B) * (t-b)/(c-b)

print(sp.pretty( F ))
print((F-E).simplify())

这输出

n⋅(-A + B)⋅(c - t)   n⋅(-B + C)⋅(-b + t)
────────────────── + ───────────────────
(-a + c)⋅(-b + c)     (-b + c)⋅(-b + d) 


0

我研究了各种选项，包括 factor()、collect() 和 apart()，但其中 none 似乎产生了与 F。关于如何进行的任何指示？

此外，我想知道是否可以通过某种方式将 Sympy 的漂亮打印功能调整为

保持分子和分母中变量的原始顺序（例如 B - A 而不是 -A + B）。目前大多数情况下顺序是颠倒的，前导减号看起来很丑。
将复合分数显示为简单分数的乘积（例如 a/b c/d 而不是 ac/bd），但在某些情况下，where/how“拆分”此类复合分数当然可能不明确分数。

Answer 1

这里的情况是 Add 有两个术语。可以使用 factor 分别简化每个术语，但每个术语的取消因素不同，因此调用 factor 总体上 Add 无法找到可能的取消。

考虑到这一点，我们需要小心地独立处理 Add 的条款，我们可以通过访问 .args:

In [122]: E.func(*(factor(term) for term in E.args))
Out[122]: 
n⋅(A - B)⋅(-c + t)   n⋅(B - C)⋅(-b + t)
────────────────── - ──────────────────
 (a - c)⋅(b - c)      (b - c)⋅(b - d)

变量的顺序实际上是由打印机在显示表达式时确定的，不一定与参数的内部顺序相同，也不一定与创建表达式时使用的顺序相同。尽管

调用 signsimp 可以规范化表达式中的减号

In [123]: signsimp(_)
Out[123]: 
  n⋅(A - B)⋅(c - t)   n⋅(B - C)⋅(b - t)
- ───────────────── + ─────────────────
   (a - c)⋅(b - c)     (b - c)⋅(b - d)

Sympy 可以通过收集多个项来简化有理表达式吗？

Can Sympy simplify a rational expression by collecting multiple terms?

python

sympy

simplify