python - 在 Numpy 中转换多项式变量

Question

我正在尝试学习如何使用 Numpy。假设我有一个多项式的根。我用coeff = np.polynomial.polynomial.polyfromroots(roots)以数组形式获取多项式的系数。然后我用 打印 np.poly1d(coeff)打印出多项式。令该多项式为

x^2 +3x + 2

现在如何转换变量使得

x is now 2/x

这就是等式变为

(2/x)^2 + 6/x + 2

在 scilab 中，我可以使用 horner 函数来完成此操作。在numpy中可以吗？

Answer 1

在 SymPy 中，这很简单:

from sympy.abc import x

f = x**2 + 3*x + 2
g = f.subs({x:2/x})

结果:

print(g)
#2 + 6/x + 4/x**2