我对简化有问题。请参阅下面的代码。
当我用某物替换sympy下的表达式时,它就变得“牢不可破”。例如,在最后一行代码中,我将表达式乘以sqrt并简化了结果。分子中的f0和分母中的f0没有得到简化,即使sqrt(f0**2)被声明为f0。我做错了什么,如何在没有这种效果的情况下替换表达式?
代码(pretty version with rendered equations):

>>> from __future__ import division
>>> from sympy import *
>>> from sympy.abc import pi

>>> init_printing(use_unicode=True)
>>> L0, Csum0, f0 = symbols("L0 C_{{\\Sigma}0} f0", nonnegative=True)
>>> equation = sqrt(L0)
>>> equation
  ____
╲╱ L₀

>>> substitute = solve(Eq(1/(2*pi*sqrt(L0)*sqrt(Csum0)), f0), L0, dict=1)[0]
>>> substitute
⎧              1           ⎫
⎪L₀: ──────────────────────⎪
⎨                      2  2⎬
⎪    4⋅C_{{\Sigma}0}⋅f₀ ⋅π ⎪
⎩                          ⎭

>>> equation = equation.subs(substitute)
>>> equation
     ______________________
    ╱          1
   ╱  ────────────────────
  ╱                   2  2
╲╱    C_{{\Sigma}0}⋅f₀ ⋅π
───────────────────────────
             2

>>> simplify(equation*f0)
        ______________________
       ╱          1
f₀⋅   ╱  ────────────────────
     ╱                   2  2
   ╲╱    C_{{\Sigma}0}⋅f₀ ⋅π
──────────────────────────────
               2

最佳答案

当f0为0时,取消f0/sqrt(f0**2)是不合法的。要确保允许取消,请将f0声明为正,而不仅仅是非负。
另外,从sympy.abc导入pi会使pi成为一个没有特定含义的通用符号;特别是,它不知道是正的。SymPy已经内置了pi(数学常数),并且知道它是一个正数。因此,从sympy.abc中删除导入可以改进简化。

from sympy import *
L0, Csum0, f0 = symbols("L0 C_{{\\Sigma}0} f0", positive=True)
equation = sqrt(L0)
substitute = solve(Eq(1/(2*pi*sqrt(L0)*sqrt(Csum0)), f0), L0, dict=1)[0]
equation = equation.subs(substitute)
simplify(equation*f0)

回报
1/(2*pi*sqrt(C_{{\Sigma}0}))

关于python - Sympy不会用符号除以其平方的平方根来简化表达式,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/49956507/

10-11 22:05