如何解决从 sympy 的系列中导出所有系数
所以我一直在尝试从 sympy 中的级数展开中推导出所有系数。考虑
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 函数工作。有没有办法解决?得到我的表达式的系数!
是
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。