scipy.special.

roots_hermite#

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

高斯埃尔米特(物理学家)正交。

计算高斯埃尔米特正交的样本点和权重。样本点是 n 次埃尔米特多项式的根,\(H_n(x)\)。这些样本点和权重可以正确积分区间 \([-\infty, \infty]\) 上度数小于或等于 \(2n - 1\) 的多项式,其中权重函数为 \(w(x) = e^{-x^2}\)。详情请参阅 [AS] 中的 22.2.14。

参数:
nint

正交阶

mubool,可选

如果为 True,返回权重之和(可选)。

返回:
xndarray

样本点

wndarray

权重

mufloat

权重之和

备注

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

对于大于 150 的 n,将应用以数值稳定的方式计算节点和权重的最佳渐近算法。该算法具有线性运行时,使得对很大的 n(几千个或更多)进行计算变得可行。

参考文献

[townsend.trogdon.olver-2014]

A. Townsend;T. Trogdon;S. Olver(2014 年)实数线上的高斯积分节点和权重的快速计算arXiv:1410.5286

[townsend.trogdon.olver-2015]

A. Townsend,T. Trogdon 和 S. Olver(2015 年)实数线上的高斯积分节点和权重的快速计算。IMA 数值分析杂志 DOI:10.1093/imanum/drv002

[AS]

Milton Abramowitz 和 Irene A. Stegun 编辑, 带有公式、图表和数学表格的数学函数手册。纽约:多佛,1972 年。