主題
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.tw/Legendre-GaussQuadrature.html

主題分類