主題
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 求積

使用 探索

參考文獻

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.

在 上引用

高斯求積

引用為

Weisstein, Eric W. “高斯求積。” 來自 Web 資源。 https://mathworld.tw/GaussianQuadrature.html

主題分類