Question

我在使用Sympy来解决方程时遇到了麻烦。当我运行代码时，例如：

print(correction(10))

我希望它打印一个f的数字。相反，它给了我错误：执行中止。

def correction(r):

    from sympy import cosh, log, exp, symbols, solve
    f = symbols('f')

    def equation():
        return cosh(((r - 1.0)/(r + 1.0))*(log(2.0)/f)) - 0.5*exp(log(2.0)/f)

    correction = solve(equation(),f)
    return correction

有什么问题？

Answer 1

你的等式是高度非线性的，我的猜测是无法找到封闭形式的解决方案。这就是sympy.solve失败的原因。您唯一的选择是以数字方式求解方程。 Sympy为此目的提供nsolve函数，这在数值求解器中是典型的，需要对解决方案进行估计。

import sympy as sp
r, f = sp.symbols('r, f')
expr = sp.cosh(((r - 1)/(r + 1))*(sp.log(2)/f)) - sp.Rational(1,2)*sp.exp(sp.log(2)/f)

sol = sp.nsolve(expr.subs({r:10}), f, 0.5)
print(sol)

0.699259455239414

solve（）执行中止

1 个答案: