python - 用 SymPy 求解方程的两个边

Question

我有以下公式来计算电子电路中的边界电感。使用 SymPy，当我有 D 时，我可以使用 Lb.evalf(subs={...}) 计算 Lb 值，f 和 R。

Lb = ((1-D)**2*D*R)/(2*f)

如果供应Lb 值？例如，我知道等式中的所有值，除了 D。

Answer 1

您可以使用 sympy 求解变量。例如:

import sympy
from sympy.abc import d,f,r,l
equation = sympy.Eq( ((1-d)**2*d*r)/(2*f), l )

print sympy.solve(equation,"f")
print sympy.solve(equation,"r")
print sympy.solve(equation,"d")

这给出了分别根据 f 和 r 解析求解的方程。

[d*r*(d - 1)**2/(2*l)]
[2*f*l/(d*(d - 1)**2)]

请注意，求解 d 的结果是 mess of solutions .这是意料之中的，因为方程是 d 的三次方。

[-(-1/2 - sqrt(3)*I/2)*(-f*l/r + sqrt((-2*f*l/r + 2/27)**2/4 - 1/729) + 1/27)**(1/3) + 2/3 - 1/(9*(-1/2 - sqrt(3)*I/2)*(-f*l/r + sqrt((-2*f*l/r + 2/27)**2/4 - 1/729) + 1/27)**(1/3)), -(-1/2 + sqrt(3)*I/2)*(-f*l/r + sqrt((-2*f*l/r + 2/27)**2/4 - 1/729) + 1/27)**(1/3) + 2/3 - 1/(9*(-1/2 + sqrt(3)*I/2)*(-f*l/r + sqrt((-2*f*l/r + 2/27)**2/4 - 1/729) + 1/27)**(1/3)), -(-f*l/r + sqrt((-2*f*l/r + 2/27)**2/4 - 1/729) + 1/27)**(1/3) + 2/3 - 1/(9*(-f*l/r + sqrt((-2*f*l/r + 2/27)**2/4 - 1/729) + 1/27)**(1/3))]