solve（）執行中止

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