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
积分阶数
- mu布尔值,可选
如果为 True,则返回权重之和,可选。
- 返回值:
- xndarray
样本点
- wndarray
权重
- mu浮点数
权重之和
注意
对于 150 以下的小 n,将使用 Golub-Welsch 算法的修改版。从特征值问题计算节点,并通过牛顿迭代的一步来改进。权重由众所周知的分析公式计算。
对于大于 150 的 n,使用了经过数次稳定计算节点和权重的理想渐进算法。该算法拥有线性的运行时间,使非常大的 n(几千或更多)的计算变得可行。
参考资料
[AS]Milton Abramowitz 和 Irene A. Stegun 编辑。带公式、图形和数学表的数学函数手册。纽约:多佛,1972 年。