Noname manuscript No. (will be inserted by the editor) How much faster does the best polynomial approximation converge than Legendre projection? Haiyong Wang Received: date / Accepted: date Abstract We compare the convergence behavior of best polynomial approx- imations and Legendre and Chebyshev projections and derive optimal rates of convergence of Legendre projections for analytic and differentiable func- tions in the maximum norm. For analytic functions, we show that the best polynomial approximation of degree n is better than the Legendre projection of the same degree by a factor of n1/2. For differentiable functions such as piecewise analytic functions and functions of fractional smoothness, however, we show that the best approximation is better than the Legendre projection by only some constant factors. Our results provide some new insights into the approximation power of Legendre projections. Keywords Legendre projection best polynomial approximation Chebyshev projection optimal rate of convergence· analytic functions differentiable· functions · · · Mathematics Subject Classification (2010) 41A25 41A10 · 1 Introduction The Legendre polynomials are one of the most important sequences of orthog- onal polynomials which have been extensively used in many branches of sci- entific computing such as Gauss-type quadrature, special functions, p-version of the finite element method and spectral methods for differential and integral This work was supported by National Natural Science Foundation of China under grant number 11671160. Haiyong Wang arXiv:2001.01985v3 [math.NA] 31 Dec 2020 School of Mathematics and Statistics, Huazhong University of Science and Technology, Wuhan 430074, P. R. China. E-mail: [email protected] Hubei Key Laboratory of Engineering Modeling and Scientific Computing, Huazhong Uni- versity of Science and Technology, Wuhan 430074, P. R. China. 2 Haiyong Wang equations (see, e.g., [6,9,11,13,14,21,24,25,27]). Among these applications, Legendre polynomials are particularly appealing owing to their superior prop- erties: (i) they have excellent error properties in the approximation of a glob- ally smooth function; (ii) quadrature rules based on their zeros or extrema are optimal in the sense of maximizing the exactness of integrated polynomials; (iii) they are orthogonal with respect to the uniform weight function ω(x)=1 which makes them preferable in Galerkin methods for PDEs. Let n 0 be an integer and let P (x) denote the Legendre polynomial of ≥ n degree n which is normalized by Pn(1) = 1. The sequence of Legendre poly- nomials Pn(x) forms a system of polynomials orthogonal over Ω = [ 1, 1] and { } − 1 2 P (x)P (x)dx = δ , (1) n m 2n +1 mn Z−1 where δmn is the Kronecker delta [20, p. 14]. Given a real-valued function f(x) which belongs to a Lipschitz class of order larger than 1/2 on Ω, then it has the following uniformly convergent Legendre series expansion [26] ∞ 1 1 f(x)= akPk(x), ak = k + f(x)Pk(x)dx. (2) 2 −1 Xk=0 Z Let (f) denote the truncated Legendre expansion of degree n, i.e., Pn n (f)= a P (x). (3) Pn k k Xk=0 which is also known as the Legendre projection. It is well known that this polynomial is the best polynomial approximation to f(x) in the L2 norm with respect to the Legendre weight ω(x) = 1. The computation of the first n +1 n Legendre coefficients ak k=0 has received much attention over the past decade and fast algorithms have{ } been developed in [2,15] that require only O(n log n) operations and in [29] that require O(n log2 n) operations. Besides Legendre polynomials, another widely used sequence of orthogo- nal polynomials is the Chebyshev polynomials, i.e., Tk(x) = cos(k arccos(x)). Suppose that f(x) is Dini-Lipschitz continuous on Ω, then it has the following uniformly convergent Chebyshev series [18, Theorem 5.7] ∞ 1 ′ 2 f(x)Tk(x) f(x)= ckTk(x), ck = dx, (4) 2 π −1 √1 x kX=0 Z − where the prime indicates that the first term of the sum is halved. Let n(f) denote the truncated Chebyshev expansion of degree n, i.e., T n (f)= a T (x), (5) Tn k k Xk=0 How much faster does the best approximation converge than Legendre 3 which is also known as the Chebyshev projection. It is well known that n(f) is the best polynomial approximation to f(x) in the L2 norm with respectT to the Chebyshev weight ω(x) = (1 x2)−1/2 and the first n + 1 Chebyshev n − coefficients ck k=0 can be evaluated efficiently by making use of the FFT in only O(n log{n)} operations (see, e.g., [18, Section 5.2.2]). Let n(f) denote the best approximation polynomial of degree n to f on Ω = [ 1B, 1] in the maximum norm, i.e., − f n(f) ∞ = min f p ∞, k −B k p∈Πn k − k where Πn denotes the space of polynomials of degree at most n. If f is contin- uous on Ω, it is well known that n(f) exists and is unique. From the point of view of polynomial approximationB in the maximum norm, it is clear that (f) is more accurate than (f) and (f). However, explicit expressions Bn Pn Tn for n(f) are generally impossible to obtain since the dependence of n(f) on f isB nonlinear and Remez-type algorithms, which are realized by iterativB e pro- cedures, have been developed for computing n(f) (see, e.g., [30, Chapter 10]). Although algorithms are available, they areB still time-consuming when n is in the thousands or higher. Obviously, this leads us to face an inevitable dilemma of whether the increase in accuracy is sufficient to justify the extra cost of com- puting n(f). WithB these three approaches, a natural question is: How much better is the accuracy of n(f) than n(f) and n(f) in the maximum norm? For the case of (f) whereB f C(Ω),T it has beenP shown in [22, Theorem 2.2] that the Tn ∈ maximum error of n(f) is inferior to that of n(f) by at most a logarithmic factor, i.e., T B 4 f (f) log n +4 f (f) . (6) k − Tn k∞ ≤ π2 k −Bn k∞ For the case of n(f), it has been widely reported that the maximum error of P 1/2 n(f) is inferior to that of n(f) by at most a factor of n (see, e.g., [19,31, 33,35]).P We summarize hereB existing results from two perspectives: – For f C(Ω), it is well known that ∈ f (f) (1 + Λ ) f (f) , (7) k −Pn k∞ ≤ n k −Bn k∞ where Λn = supf6=0 n(f) ∞/ f ∞ is the Lebesgue constant of n(f). Furthermore, Qu andkP Wongk in [19,k k Equation (1.10)] showed that P 1 3/2 n +1 (1,0) 2 1/2 Λn = Pn (x) dx = n + O(1), 2 −1 √π Z (1,0) where Pn (x) is the Jacobi polynomial of degree n with α = 1 and β =0 and α, β are the parameters in Jacobi polynomials. Hence we can conclude that the rate of convergence of n(f) is slower than that of n(f) by a factor of n1/2. P B 4 Haiyong Wang – Under the assumption that f,f ′,...,f (m−1) are absolutely continuous, f (m) is of bounded variation and f (m) < where m 1 is an in- k kT ∞ ≥ teger and T denotes some weighted semi-norm. It has been shown in k·k −m−1/2 [31,33] that the Legendre coefficients of f satisfy ak = O(k ). As a direct consequence we obtain | | ∞ f (f) a = O(n−m+1/2), (8) k −Pn k∞ ≤ | k| k=Xn+1 where we have used the inequality Pk(x) 1 (see, e.g., [25, p. 94]). | | ≤ −m Notice that the rate of convergence of n(f) for such functions is O(n ) as n [28, Chapter 7]. Again, weB see that the rate of convergence of (f)→ is ∞ slower than that of (f) by a factor of n1/2. Pn Bn Is the rate of convergence of n(f) really slower than n(f) by a factor of n1/2? Let us consider a motivatingP example f(x) = x B, which is absolutely continuous on Ω and its first-order derivative is of bounded| | variation. More- over, it has been shown in [33, Equation (2.11)] that the Legendre coefficients of f satisfy the following sharp bound −1 4 1 −3/2 ak k = O(k ), (9) | |≤ π(2k 3) − 2 − where k 2 is even andp a = 0 when k is odd. We now consider the rate ≥ k of convergence of n(f), n(f) and n(f). For n(f) and n(f), it is well know that their ratesB of convergenceT P are O(n−1)B as n T (see, e.g., [30, → ∞ Chapter 7]). For n(f), however, from (7) and (8) we can deduce that the predicted rate ofP convergence of (f) is only O(n−1/2). Unexpectedly, we Pn observed in [33, Figure 3] that the rate of convergence of n(f) is actually −1 P O(n ) as n , which is the same as that of n(f) and n(f). This unex- pected observation→ ∞ suggests that existing resultsB on the rateT of convergence of n(f) may be suboptimal. P In this paper, we aim to investigate the optimal rate of convergence of (f) in the maximum norm.