主题
Search

勒让德-高斯求积

DOWNLOAD Mathematica Notebook下载 Wolfram 笔记本

勒让德-高斯求积是一种数值积分方法，也称为“高斯”求积或勒让德求积。在区间 [-1,1] 上，权重函数为 W(x)=1 的 高斯求积。求积阶数 n 的 求积节点 由 勒让德多项式 P_n(x) 的根给出，这些根关于 0 对称分布。权重为

w_i=-(A_(n+1)gamma_n)/(A_nP_n^'(x_i)P_(n+1)(x_i))
(1)
=(A_n)/(A_(n-1))(gamma_(n-1))/(P_(n-1)(x_i)P_n^'(x_i)),
(2)

其中 A_n 是 系数 x^n 在 P_n(x) 中。对于 勒让德多项式，

A_n=((2n)!)/(2^n(n!)^2)
(3)

（Hildebrand 1956, 第 323 页），因此

(A_(n+1))/(A_n)=([2(n+1)]!)/(2^(n+1)[(n+1)!]^2)(2^n(n!)^2)/((2n)!)
(4)
=(2n+1)/(n+1).
(5)

此外，

gamma_n=int_(-1)^1[P_n(x)]^2dx
(6)
=2/(2n+1)
(7)

（Hildebrand 1956, 第 324 页），因此

w_i=-2/((n+1)P_(n+1)(x_i)P_n^'(x_i))
(8)
=2/(nP_(n-1)(x_i)P_n^'(x_i)).
(9)

使用 递推关系

(1-x^2)P_n^'(x)=-nxP_n(x)+nP_(n-1)(x)
(10)
=(n+1)xP_n(x)-(n+1)P_(n+1)(x)
(11)

（更正了 Hildebrand 1956, 第 324 页）得到

w_i=2/((1-x_i^2)[P_n^'(x_i)]^2)
(12)
=(2(1-x_i^2))/((n+1)^2[P_(n+1)(x_i)]^2)
(13)

（Hildebrand 1956, 第 324 页）。

权重 w_i 满足

sum_(i=1)^nw_i=2,
(14)

这由恒等式得出

sum_(nu=1)^n(1-x_nu^2)/((n+1)^2[P_(n+1)(x_nu)]^2)=1.
(15)

误差项为

E=(2^(2n+1)(n!)^4)/((2n+1)[(2n)!]^3)f^((2n))(xi).
(16)

Beyer (1987) 给出了高达 n=16 的 求积节点 和权重的表格，Chandrasekhar (1960) 给出了高达 n=8（对于 n 偶数）的表格。

nx_iw_i
2+/-0.577351.000000
300.888889
+/-0.7745970.555556
4+/-0.3399810.652145
+/-0.8611360.347855
500.568889
+/-0.5384690.478629
+/-0.906180.236927

精确的 求积节点 在下表中给出。

nx_iw_i
2+/-1/3sqrt(3)1
308/9
+/-1/5sqrt(15)5/9
4+/-1/(35)sqrt(525-70sqrt(30))1/(36)(18+sqrt(30))
+/-1/(35)sqrt(525+70sqrt(30))1/(36)(18-sqrt(30))
50(128)/(225)
+/-1/(21)sqrt(245-14sqrt(70))1/(900)(322+13sqrt(70))
+/-1/(21)sqrt(245+14sqrt(70))1/(900)(322-13sqrt(70))

阶数为 n 的求积的求积节点是 勒让德多项式 P_n(x) 的根，这意味着它们是次数为 1, 2, 2, 4, 4, 6, 6, 8, 8, 10, 10, 12, ... 的代数数，对于 n>1，其等于 2|_n/2_| (OEIS A052928)。

类似地，阶数为 n 的求积的权重可以表示为次数为 1, 1, 1, 2, 2, 3, 3, 4, 4, 5, 5, ... 的多项式的根，对于 n>1，其等于 |_n/2_| (OEIS A008619)。 根确定权重的多项式三角形是

x-2
(17)
x-1
(18)
9x-5
(19)
216x^2-216x+49
(20)
45000x^2-32200x+5103
(21)
2025000x^3-2025000x^2+629325x-58564
(22)
142943535000x^3-113071253400x^2+27510743799x-1976763932
(23)
1707698764800000x^4-1707698764800000x^3+606530263046400x^2-88878097916608x+4373849390625
(24)

（OEIS A112734）。

使用 探索

参考文献

Abbott, P. "Tricks of the Trade: Legendre-Gauss Quadrature." Mathematica J. 9, 689-691, 2005.Beyer, W. H. CRC Standard Mathematical Tables, 28th ed. Boca Raton, FL: CRC Press, pp. 462-463, 1987.Chandrasekhar, S. Radiative Transfer. New York: Dover, pp. 56-62, 1960.Hildebrand, F. B. Introduction to Numerical Analysis. New York: McGraw-Hill, pp. 323-325, 1956.Sloane, N. J. A. 序列 A008619, A052928, 和 A112734，出自 "The On-Line Encyclopedia of Integer Sequences."

在 中被引用

勒让德-高斯求积

请这样引用

Weisstein, Eric W. "勒让德-高斯求积。" 来自 Web 资源。 https://mathworld.net.cn/Legendre-GaussQuadrature.html

主题分类