如何获得为函数的每个项提供指数数的多项式?

时间:2017-12-06 11:57:30

标签: python numpy scipy polynomials

我试图写一个多项式。然而,我发现scipy方法给出了多项式格式,如x ^ 4 + 2x ^ 3 - x ^ 2 + 3,系数为[3,0,-1,2,1],指数为4 ,3和2。

有没有办法通过给定的数字列表来操纵指数?有点像,给予:

nj = [0.0, 1.0, 2.0, 2.0, 2.0, 3.0, 3.0, 3.0, 4.0, 4.0, 4.0, 4.0, 4.0, 4.0, 5.0, 5.0, 5.0, 6.0, 6.0, 6.0, 7.0, 8.0]

例如,对于给定的唯一系数,我可以计算多项式的项,如coef*x^nj = 0

谢谢!

1 个答案:

答案 0 :(得分:0)

如果您有NumPy系数数组coef和NumPy指数数组nj,则x处的多项式值可以计算为

coef.dot(x**nj)

示例:

coef = np.array([3.1, 0, -2, 5, 1.1])
nj = np.array([0, 1, 2, 2, 5])   
x = 1.23
print(coef.dot(x**nj))   # 10.7355362527