如何找到擬合指數函數的x？

Question

fp.append(np.polyfit(train_x, train_y, 2)) 
f.append(np.poly1d(fp)) 
print(np.poly1d(fp)) 
threshold = fsolve(f, 50) 

上述代碼成功找到y = 50的x值。但是當我試圖對擬合的指數函數做同樣的事情時，我無法理解如何做到這一點。

def f_exp(x, a, b, c): 
    y = a * np.exp(-b * x) + c 
    return y 
popt, pcov = curve_fit(f_exp, train_x, train_y) 
print(popt) 
threshold = fsolve(f_exp, 50) fails with :TypeError: f_exp() missing 3 required positional arguments: 'a', 'b', and 'c' 

如果我添加* POPT然後我得到

threshold = fsolve(f_exp(*popt), 50) fails with: TypeError: f_exp() missing 1 required positional argument: 'c' 

我認爲我需要補充的x值，但它是我試圖找到...無論如何價值，增加了一些價值代替X，導致另一個錯誤：

threshold = fsolve(f_exp(1, *popt), 50) fails with: TypeError: 'numpy.float64' object is not callable 

Answer 1

我想你需要與優化參數傳遞f_exp功能fsolve（即，具有一個，從獲得設置爲值b和c ARGS）。爲了做到這一點，你可以使用functools.partial功能：

popt, pcov = curve_fit(f_exp, train_x, train_y) 
print(popt) 
import functools 
# preparing arguments 
kwargs = dict(zip(['a', 'b', 'c'], popt)) 
optimized_f_exp = functools.partial(f_exp, **kwargs) 
threshold = fsolve(optimized_f_exp, 50) 

我們做了這裏的部分固定a，b和你原來的功能c ARGS基本上做出了新的功能optimized_f_exp是popt（什麼也正是爲什麼稱爲部分）。