scipy.special.

roots_hermitenorm#

scipy.special.roots_hermitenorm(n, mu=False)[源代码]#

高斯-埃尔米特(统计学家)求积。

计算高斯-埃尔米特求积的采样点和权重。采样点是 n 次埃尔米特多项式 \(He_n(x)\) 的根。这些采样点和权重可以正确积分区间 \([-\infty, \infty]\) 上次数为 \(2n - 1\) 或更小的多项式,权重函数为 \(w(x) = e^{-x^2/2}\)。有关更多详细信息,请参见 [AS] 中的 22.2.15。

参数:
nint

求积阶数

mubool,可选

如果为 True,则返回权重的总和,可选。

返回:
xndarray

采样点

wndarray

权重

mufloat

权重之和

注释

对于小于等于 150 的小 n,使用 Golub-Welsch 算法的修改版本。节点从特征值问题计算得出,并通过牛顿迭代的一步进行改进。权重从众所周知的解析公式计算得出。

对于大于 150 的 n,使用最优渐近算法,该算法以数值稳定的方式计算节点和权重。该算法具有线性运行时,使得计算非常大的 n(几千甚至更多)成为可能。

参考文献

[AS]

Milton Abramowitz 和 Irene A. Stegun 编辑。《数学函数手册,包含公式、图形和数学表》。纽约:Dover,1972 年。