solve（）执行中止

Question

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

print(correction(10))

我希望它能打印出一个f的数字，但它却给了我一个错误：executionaborted。在

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