python - 在其他多项式上评估 numpy 多项式

Question

numpy.lib.polynomial.polyval 允许您使用另一个多项式评估多项式:

numpy.polyval(poly1d([1, 2, 3]), 2)
Out[832]: 11

numpy.polyval(poly1d([1, 1]), poly1d([1, 1, 1]))
Out[820]: poly1d([ 1.,  1.,  2.])

例如。如何使用 numpy.polynomial.polynomial.polyval 执行相同操作？

numpy.polynomial.polynomial.polyval(2, [3, 2, 1])
Out[833]: 11.0

numpy.polynomial.polynomial.polyval(Polynomial([1, 1, 1]), Polynomial([1, 1]))
Out[834]: Polynomial([ 1.,  1.], [-1.,  1.], [-1.,  1.])

Answer 1

最简单的方法是使用 Polynomial 类。

In [1]: from numpy.polynomial import Polynomial as P

In [2]: p1 = P([1,1])

In [3]: p2 = P([1,1,1])

In [4]: p2(p1)
Out[4]: Polynomial([ 3.,  3.,  1.], [-1.,  1.], [-1.,  1.])

In [5]: p1(p2)
Out[5]: Polynomial([ 2.,  1.,  1.], [-1.,  1.], [-1.,  1.])

如果您坚持使用 polyval，则需要系数和用作 x 的多项式。

In [12]: import numpy.polynomial.polynomial as poly

In [13]: poly.polyval(p1, [1,1,1])
Out[13]: Polynomial([ 3.,  3.,  1.], [-1.,  1.], [-1.,  1.])

In [14]: poly.polyval(p2, [1,1])
Out[14]: Polynomial([ 2.,  1.,  1.], [-1.,  1.], [-1.,  1.])