从 Sympy 中的级数展开中获取 n 阶系数

Question

我有一个使用 numpy 的 poly1d 定义的多项式：

import numpy as np

p = np.poly1d([-1.37970733e+07,  8.02055768e+07, -1.98442743e+08,  2.70981281e+08,

       -2.20611918e+08,  1.07019785e+08, -2.85873914e+07,  3.24251782e+06])

然后我对这个多项式关于我称为 的点进行级数展开radMax = 0.80868373。以下是我进行系列扩展的方法：

from sympy import series, Symbol, Poly

x = Symbol('x')

polyExpand = p(x).series(x,radMax,3)/(p(radMax)/0.0523772)

这给了我想要的展开式，我想要二阶项的系数。如果运行上面的代码，您将看到该系数为“0.2843162879664”。我想提取这个系数。

Answer 1

如果您想按项处理系数：from sympy import sympifyprint(polyExpand.as_coefficients_dict()[sympify(f"(x - {radMax}) ** 2")])如果您想通过索引获取系数：coeffs = list(polyExpand.as_coefficients_dict().values())print(coeffs[2]) # order n = 2这些将正确输出-0.284316287966400。你忘记了负号。