Is Gauss quadrature better than Clenshaw-Curtis?

被引:441
作者
Trefethen, Lloyd N. [1 ]
机构
[1] Oxford Comp Lab, Oxford OX1 3QD, England
关键词
Gauss quadrature; Newton-Cotes; Clenshaw-Curtis; Chebyshev expansion; rational approxiination; FFT; spectral methods;
D O I
10.1137/060659831
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
We compare the convergence behavior of Gauss quadrature with that of its younger brother, Clenshaw-Curtis. Seven-line MATLAB codes are presented that implement both methods, and experiments show that the supposed factor-of-2 advantage of Gauss quadrature is rarely realized. Theorems are given to explain this effect. First, following O'Hara and Smith in the 1960s, the phenomenon is explained as a consequence of aliasing of coefficients in Chebyshev expansions. Then another explanation is offered based on the interpretation of a quadrature formula as a, rational approximation of log((z + 1)/(z - 1)) in the complex plane. Gauss quadrature corresponds to Pade approximation at z = infinity. Clenshaw-Curtis quadrature corresponds to an approximation whose order of accuracy at z = infinity is only half as high, but which is nevertheless equally accurate near [-1, 1].
引用
收藏
页码:67 / 87
页数:21
相关论文
共 61 条
[1]  
Acton F.S., 1990, Numerical Methods That Work
[2]   Hybrid Gauss-trapezoidal quadrature rules [J].
Alpert, BK .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 1999, 20 (05) :1551-1584
[3]  
[Anonymous], 1971, P STEKLOV I MATH+
[4]  
[Anonymous], 1964, NUMER MATH, DOI DOI 10.1007/BF01386080
[5]  
[Anonymous], 1982, NUMERICAL ANAL
[6]  
[Anonymous], 1925, B AM MATH SOC, DOI DOI 10.1090/S0002-9904-1925-03998-1
[7]  
[Anonymous], NUMERICAL COMPUTATIO
[8]  
[Anonymous], 1975, METHODS NUMERICAL IN
[9]  
[Anonymous], SPECTRAL METHODS FUN
[10]  
Baker Jr G. A., 1996, PADE APPROXIMANTS