什麼是高斯-勒讓德求積法計算器?
這個工具可在參考區間 [-1, 1] 上計算 n 點高斯-勒讓德求積法的節點(橫坐標)與權重。高斯-勒讓德求積法是一種數值積分方法,它將定積分近似為函數值的加權總和:f(x) 從 -1 到 1 的積分約等於各項 \(w_i\) 乘上 \(f(x_i)\) 的總和。只要 \(n\) 個節點,就能精確積分次數達 \(2n-1\) 的任意多項式,因此遠比梯形法、辛普森法(Simpson)等等距節點方法更為準確。
使用方式
選擇階數 \(n\)(即節點數,範圍 2 至 100),並可選擇要顯示的有效位數。計算器會回傳一個含 \(n\) 列的表格,每列包含一個節點 \(x_i\) 及其對應權重 \(w_i\)。節點關於 0 對稱,且全部嚴格落在 (-1, 1) 區間內;權重皆為正值,總和恰為 2,也就是該區間的長度。若要在任意區間 [a, b] 上積分,請以 \(t_i = \frac{b-a}{2} \times x_i + \frac{a+b}{2}\) 對各節點做映射,並將每個權重乘上 \(\frac{b-a}{2}\)。
公式說明
節點即勒讓德多項式 \(P_n(x)\) 的 \(n\) 個根,可用 Bonnet 遞迴式建立:\(P_0=1\),\(P_1=x\),且 $$P_k = \frac{(2k-1)\, x\, P_{k-1} - (k-1)\, P_{k-2}}{k}.$$每個權重為 $$w_i = \frac{2}{\left(1 - x_i^{2}\right)\left(P_n^{\prime}(x_i)\right)^{2}},$$其中導數 $$P_n^{\prime}(x) = \frac{n\,(x\, P_n(x) - P_{n-1}(x))}{x^2 - 1}.$$求根採用牛頓法,以 \(x = \cos\!\left(\frac{\pi (i - 0.25)}{n + 0.5}\right)\) 作為初始值,通常只需幾次迭代即可收斂。
範例演算(n = 3)
\(P_3\) 的根為 \(x = 0\) 與 \(x = \pm\sqrt{3/5} = \pm 0.7745966692\)。\(x = 0\) 處的權重為 \(\frac{8}{9} = 0.888888889\),而 \(x = \pm 0.7745966692\) 處的權重各為 \(\frac{5}{9} = 0.555555556\)。權重總和為 $$\frac{5}{9} + \frac{8}{9} + \frac{5}{9} = 2,$$且此 3 點規則可精確積分次數達 5 的多項式。
常見問題
為什麼權重總和為 2?將常數函數 \(f(x) = 1\) 在 [-1, 1] 上積分結果為 2,而求積法必須能精確重現常數函數,因此所有權重的總和必定等於區間長度。
計算結果有多準確?本計算器以雙精度浮點數重建,可提供約 15 位有效數字。若顯示選項要求更多位數,超出部分會四捨五入至雙精度所能表示的範圍。
可精確積分的最高次數是多少?\(n\) 點規則對所有次數最高達 \(2n-1\) 的多項式都能精確積分。