scipy.special.

laguerre

scipy.special.laguerre(n, monic=False)

拉盖尔多项式。

定义为以下方程的解：

\[x\frac{d^2}{dx^2}L_n + (1 - x)\frac{d}{dx}L_n + nL_n = 0;\]

\(L_n\) 是一个 \(n\) 次多项式。

参数:

nint

多项式的次数。

monicbool，可选

如果为 True，则缩放前导系数为 1。默认值为 False。

返回:

Lorthopoly1d

拉盖尔多项式。

另请参阅

genlaguerre

广义（伴随）拉盖尔多项式。

注释

多项式 \(L_n\) 在 \([0, \infty)\) 上正交，权重函数为 \(e^{-x}\)。

参考文献

[AS]

Milton Abramowitz 和 Irene A. Stegun 编辑。《数学函数手册，包含公式、图表和数学表》。纽约：多佛出版社，1972 年。

示例

在您的浏览器中尝试！

拉盖尔多项式 \(L_n\) 是广义拉盖尔多项式 \(L_n^{(\alpha)}\) 的特殊情况 \(\alpha = 0\)。让我们在区间 \([-1, 1]\) 上验证一下

>>> import numpy as np
>>> from scipy.special import genlaguerre
>>> from scipy.special import laguerre
>>> x = np.arange(-1.0, 1.0, 0.01)
>>> np.allclose(genlaguerre(3, 0)(x), laguerre(3)(x))
True

多项式 \(L_n\) 还满足以下递推关系

\[(n + 1)L_{n+1}(x) = (2n +1 -x)L_n(x) - nL_{n-1}(x)\]

这可以很容易地在 \([0, 1]\) 上针对 \(n = 3\) 进行检查

>>> x = np.arange(0.0, 1.0, 0.01)
>>> np.allclose(4 * laguerre(4)(x),
...             (7 - x) * laguerre(3)(x) - 3 * laguerre(2)(x))
True

这是前几个拉盖尔多项式 \(L_n\) 的图

>>> import matplotlib.pyplot as plt
>>> x = np.arange(-1.0, 5.0, 0.01)
>>> fig, ax = plt.subplots()
>>> ax.set_ylim(-5.0, 5.0)
>>> ax.set_title(r'Laguerre polynomials $L_n$')
>>> for n in np.arange(0, 5):
...     ax.plot(x, laguerre(n)(x), label=rf'$L_{n}$')
>>> plt.legend(loc='best')
>>> plt.show()

返回在新标签页中打开