我有一個問題，我需要解決兩個方程組：B1=2.51/(Re(f^(0.5))) f=(1/(-2*log10((epsilon/D/3.7) B1)))^2 我嘗試使用 sympy、numpy 和 fsolve，但使用 math.log10 時出現問題，它回傳一個浮點數：
from sympy import*
import math
def f_cole(Epsilon, D, Re):
#Fattore di attrito Equazione di Colebrook per flusso turbolento
B1,f=symbols('B1,f')
eq1 =Eq(( 2.51 / (Re*(f**0.5))-B1),0)
eq2=Eq(((1/(-2* math.log10((Epsilon/D/ 3.7) B1)))**2 -f),0)
solveset((eq1,eq2), (B1,f))
return(f)
那回傳：
TypeError: can't convert expression to float
uj5u.com熱心網友回復：
您不能將模塊中的函式math與sympy運算式混合。math來自expect 的函式作為輸入浮動，而不是運算式。
相反，使用sympy.log.
import sympy
x = sympy.Symbol('x', real=True)
print('log(x): ', sympy.log(x))
print('log(exp(x)): ', sympy.simplify(sympy.log(sympy.exp(x))))
print('log10(x): ', sympy.log(x, 10))
print('log10(10**x): ', sympy.simplify(sympy.log(10**x, 10)))
輸出：
log(x): log(x)
log(exp(x)): x
log10(x): log(x)/log(10)
log10(10**x): x
轉載請註明出處，本文鏈接：https://www.uj5u.com/qita/449149.html