在Python中对Chebyshev多项式进行数字积分

Question

我想在Python中将Chebyshev多项式进行数字积分

这就是我正在使用的：

p = numpy.polynomial.Chebyshev.basis(5, domain = [0,1], window = [0,1])
coef = p.convert(kind=numpy.polynomial.Polynomial, domain = [0,1], window = [0,1])

我想通过使用以下积分来积分30级的切比雪夫

$I = \int_{-1}^{1} dx T_j(x)$

如何在python中完成？我找不到应该写下限和上限的方法。

Answer 1

您尝试过scipy.integrate中的功能之一吗？

例如：

>>> import numpy as np
>>> from scipy.integrate import quad
>>> p = np.polynomial.Chebyshev.basis(5, domain = [0,1], window = [0,1])
>>> quad(p,-1,1)
(0.0, 1.3873040010713506e-14)