主题
Search

高斯求积


旨在通过选择最优横坐标 x_i 来获得积分的最佳数值估计,在这些横坐标处评估函数 f(x)高斯求积基本定理 指出,m 点高斯求积 公式 的最优 横坐标 正好是对于相同区间和 权函数 的正交 多项式 的根。 高斯求积是最优的,因为它精确地拟合了所有高达 2m-1 次的多项式。 从 Radau 求积Laguerre-Gauss 求积 获得稍逊的最优拟合。

W(x)区间x_i 是...的根
1(-1,1)P_n(x)
e^(-t)(0,infty)L_n(x)
e^(-t^2)(-infty,infty)H_n(x)
(1-t^2)^(-1/2)(-1,1)T_n(x)
(1-t^2)^(1/2)(-1,1)U_n(x)
x^(1/2)(0,1)x^(-1/2)P_(2n+1)(sqrt(x))
x^(-1/2)(0,1)P_(2n)(sqrt(x))

为了确定对应于高斯横坐标 x_i的权重,通过令以下公式计算 f(x)拉格朗日插值多项式

 pi(x)=product_(j=1)^m(x-x_j)
(1)

(其中 Chandrasekhar 1967 使用 F 代替 pi),因此

 pi^'(x_j)=[(dpi)/(dx)]_(x=x_j)=product_(i=1; i!=j)^m(x_j-x_i).
(2)

然后通过 m 个点拟合拉格朗日插值多项式得到

 phi(x)=sum_(j=1)^m(pi(x))/((x-x_j)pi^'(x_j))f(x_j)
(3)

对于任意点 x。 因此,我们正在寻找一组点 x_j 和权重 w_j,使得对于 权函数 W(x)

int_a^bphi(x)W(x)dx=int_a^bsum_(j=1)^(m)(pi(x)W(x))/((x-x_j)pi^'(x_j))dxf(x_j)
(4)
=sum_(j=1)^(m)w_jf(x_j),
(5)

其中 权重

 w_j=1/(pi^'(x_j))int_a^b(pi(x)W(x))/(x-x_j)dx.
(6)

权重 w_j 有时也称为 Christoffel 数 (Chandrasekhar 1967)。 对于 正交多项式 phi_j(x),其中 j=1, ..., n,

 phi_j(x)=A_jpi(x)
(7)

(Hildebrand 1956, p. 322),其中 A_nx^nphi_n(x) 中的 系数,那么

w_j=1/(phi_n^'(x_j))int_a^bW(x)(phi(x))/(x-x_j)dx
(8)
=-(A_(n+1)gamma_n)/(A_nphi_n^'(x_j)phi_(n+1)(x)),
(9)

其中

 gamma_m=int[phi_m(x)]^2W(x)dx.
(10)

使用关系式

 phi_(n+1)(x_i)=-(A_(n+1)A_(n-1))/(A_n^2)(gamma_n)/(gamma_(n-1))phi_(n-1)(x_i)
(11)

(Hildebrand 1956, p. 323) 给出

 w_j=(A_n)/(A_(n-1))(gamma_(n-1))/(phi_n^'(x_j)phi_(n-1)(x_j)).
(12)

(注意 Press et al. 1992 省略了因子 A_n/A_(n-1)。) 在高斯求积中,权重均为 正数。 误差由下式给出

E_n=(f^((2n))(xi))/((2n)!)int_a^bW(x)[pi(x)]^2dx
(13)
=(gamma_n)/(A_n^2)(f^((2n))(xi))/((2n)!),
(14)

其中 a<xi<b (Hildebrand 1956, pp. 320-321)。

其他有趣的恒等式是

 sum_(k=0)^n([phi_k(x)]^2)/(gamma_k)=(A_n)/(A_(n+1)gamma_n)[phi_(n+1)^'(x)phi_n(x)-phi_n^'(x)phi_(n+1)(x)]
(15)

sum_(k=0)^(n)([phi_k(x_j)]^2)/(gamma_k)=-(A_nphi_n^'(x_j)phi_(n+1)(x_j))/(A_(n+1)gamma_n)
(16)
=1/(w_j)
(17)

(Hildebrand 1956, p. 323)。

在 Szegö (1975) 的 符号 中,令 x_(1n)<...<x_(nn)[a,b] 中的有序点集,令 lambda_(1n), ..., lambda_(nn) 是一组 实数。 如果 f(x)闭区间 [a,b] 上的任意函数,则将高斯求积写为

 Q_n(f)=sum_(nu=1)^nlambda_(nun)f(x_(nun)).
(18)

这里 x_(nun)横坐标lambda_(nun)Cotes 数


另请参阅

Chebyshev 求积, Chebyshev-Gauss 求积, Chebyshev-Radau 求积, 高斯求积基本定理, Hermite-Gauss 求积, Jacobi-Gauss 求积, Laguerre-Gauss 求积, Legendre-Gauss 求积, Lobatto 求积, Radau 求积

使用 Wolfram|Alpha 探索

参考文献

Abramowitz, M. and Stegun, I. A. (编). Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, 第 9 版. New York: Dover, pp. 887-888, 1972.Acton, F. S. Numerical Methods That Work, 第 2 版. Washington, DC: Math. Assoc. Amer., p. 103, 1990.Arfken, G. "Appendix 2: Gaussian Quadrature." Mathematical Methods for Physicists, 第 3 版. Orlando, FL: Academic Press, pp. 968-974, 1985.Beyer, W. H. CRC Standard Mathematical Tables, 第 28 版. Boca Raton, FL: CRC Press, p. 461, 1987.Chandrasekhar, S. An Introduction to the Study of Stellar Structure. New York: Dover, 1967.Gauss, C. F. "Methodus nova integralium valores per approximationem inveniendi." Commentationes Societatis regiae scientarium Gottingensis recentiores 3, 39-76, 1814. Reprinted in Werke, 卷 3. New York: George Olms, p. 163, 1981.Golub, G. H. and Welsh, J. H. "Calculation of Gauss Quadrature Rules." Math. Comput. 23, 221-230, 1969.Hildebrand, F. B. Introduction to Numerical Analysis. New York: McGraw-Hill, pp. 319-323, 1956.Press, W. H.; Flannery, B. P.; Teukolsky, S. A.; and Vetterling, W. T. "Gaussian Quadratures and Orthogonal Polynomials." §4.5 in Numerical Recipes in FORTRAN: The Art of Scientific Computing, 第 2 版. Cambridge, England: Cambridge University Press, pp. 140-155, 1992.Stroud, A. H. and Secrest, D. Gaussian Quadrature Formulas. Englewood Cliffs, NJ: Prentice-Hall, 1966.Szegö, G. Orthogonal Polynomials, 第 4 版. Providence, RI: Amer. Math. Soc., pp. 37-48 and 340-349, 1975.Whittaker, E. T. and Robinson, G. "Gauss's Formula of Numerical Integration." §80 in The Calculus of Observations: A Treatise on Numerical Mathematics, 第 4 版. New York: Dover, pp. 152-163, 1967.

在 Wolfram|Alpha 上引用

高斯求积

引用为

Weisstein, Eric W. “高斯求积。” 来自 MathWorld--Wolfram Web 资源。 https://mathworld.net.cn/GaussianQuadrature.html

主题分类