Python 从sympy'导出所有系数；s系列

所以我一直在试着从辛相中的级数展开式推导出所有的系数。考虑

import numpy as np
from sympy import *
import matplotlib.pyplot as plt
from numpy.linalg import matrix_power

L     = symbols("L")
alpha = 1 -1.3*L + 0.4*L**2
roots = solve(alpha,L)
#Sanity checks
factorz_1 = factor(L**2 - 3.25*L + 2.5)
factorz_2 = factor(0.4*L**2 - 1.3*L + 1)
##GET COEFFICIENTS
coeff_1 = 1/roots[0]
coeff_2 = 1/roots[1]

##GET COEFFICIENTS OF SERIES
series_1 = series(1/(1-coeff_1*L))
series_2 = series(1/(1-coeff_2*L))
Lag_inv  = expand(series(1/(1-coeff_1*L),n=7)*series(1/(1-coeff_2*L),n=7))
Lag_poly = Poly(Lag_inv,L)
coeffs   = Lag_poly.coeffs()

我从定义alpha开始，通过solve函数得到根。忽略健全性检查。“获取系数”部分是获取数字，稍后我将根据这些数字构建序列。之后，我使用series函数定义幂级数。我想要的是提取其中给定的系数，直到我选择的顺序。这里的问题是Lag_inv是一个包含O（L^7）表达式的序列，我假设它不允许Poly函数工作。这有什么关系吗？得到我表达式的系数，就是

---

Y

您可以使用

Lag\u inv.removeO（）

。