如何使用broyden1算法计算指数值

Question

我必须在等式中找到温度：u = 1.2*e^(-.025*T)u = viscosity in N*s/m^2T = degree C使用scipy.optimize.broyden1()当 u=.001 时，您将如何使用它来找到 T？

Answer 1

要使用 broyden1 函数，只需给它提供您想要找到根的函数以及一个初始值。如果您只对求解 u=.001 感兴趣，您可以定义一个以 T 作为输入的函数。如果您可能想更改 u 的值，partial 可能是一个方便的工具。

import scipy.optimize

import numpy as np

from functools import partial

def findtemp(T, u):

    return u - 1.2 * np.exp(-.025*T)

sol = scipy.optimize.broyden1(partial(findtemp, u=.001), 0)

print(sol)

给我一个接近理论值的值 283.4828690696808。