使用具有sency函数的fSolutions

Question

我遇到了以下问题与scipy.fsolve，但我不知道该怎么办：

U = 0.00043
ThC =1.19
Dist = 7   
IncT = 0.2  
pcw = 1180000          
k = 1.19                    
B = U * pcw / (2 * k)

fugato = fsolve((((Ql/(2*math.pi* k))*math.exp(B * x)*special.kv(0, B * x))-IncT),0.01)

print fugato

我得到的错误TypeError：'numpy.float64'对象是不可调用在解决。

如何解决此问题？

Answer 1

fsolve的参数必须是函数。

我假设你想解x？若有，请书面说明：

fugato = fsolve(lambda x: Ql/(2*math.pi* k)*math.exp(B * x)*special.kv(0, B * x)-IncT,
                0.01)

作品

为了解释这里发生了什么，构造lambda x:2*x是一个函数定义。这类似于写作：

def f(x):
    return 2*x

lambda构造通常用于定义只需要一次的函数。注册回调或表示数学表达式时经常会出现这种情况。例如，如果你想集成f（x）=2*x，你可以写：

from scipy.integrate import quad
integral = quad(lambda x: 2*x, 0., 3.)

类似地，如果要求解2*x=1，可以编写：

from scipy.optimize import fsolve
fsolve(lambda x: 2*x-1, 0.)