Python scipy：**或pow()的不支持操作数类型：'list'和'list'。

Question

我需要将函数适配到数据数组中，并获取该函数方程的最优系数。我使用scipy库中的curve_fit方法，该方法基于最小二乘法。

import numpy as np 
from scipy.optimize import curve_fit

#This is my function from which i need to get optimal coefficients 'a' and 'b'
def func(x, a, b):  
return a*x**(b*x)

#the arrays of input data                               
x = [1,2,3,4,5]
y =[6,7,8,9,10]

#default (guess) coefficients
p0 = [1, 1] 

popt, pcov = curve_fit(func, x, y, p0)
print popt

执行代码后出现以下错误：

TypeError: 不支持的操作数类型，使用 ** 或 pow()：'list' 和 'list'

但是，当我使用另一个更简单的函数而没有幂运算时，它可以正常工作。

def func(x, a, b):  
return a*x + b

它必须试图将一个数的幂次应用于整个输入数据数组

怎么办？帮帮忙...

- Vladimir

---

如何将一个数组的元素按照另一个数组的幂次进行计算？ - Ignacio Vazquez-Abrams

Answer 1

你看起来是想进行逐元素的幂运算？

比如对于每个i执行 a*x[i]**(b*x[i])？

在这种情况下，你需要使用 np.power 函数：

def func(x,a,b):
    return a*np.power(x,b*x)

然后它就工作了。

(顺便提一下，将 x 和 y 从列表转换为numpy数组可能是值得的: np.array(x))。

- mathematical.coffee

---

我尝试了这个： def f(x): return 1 / (x ** 5 * (np.power(np.exp, 1/x))) ，但对我没有用。 - Sergo