Next Article in Journal
Dynamic Analysis of a Stochastic Rumor Propagation Model with Regime Switching
Next Article in Special Issue
Solution of the Problem P = L
Previous Article in Journal
Steffensen Type Inequalites for Convex Functions on Borel σ-Algebra
Previous Article in Special Issue
Roots of Characteristic Polynomial Sequences in Iterative Block Cyclic Reductions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Approximation Properties of Chebyshev Polynomials in the Legendre Norm

1
Department of Mathematics, Shanghai University, Shanghai 200444, China
2
School of Computer Science and Technology, Shandong Technology and Business University, Yantai 264000, China
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(24), 3271; https://doi.org/10.3390/math9243271
Submission received: 22 November 2021 / Revised: 10 December 2021 / Accepted: 11 December 2021 / Published: 16 December 2021
(This article belongs to the Special Issue Polynomial Sequences and Their Applications)

Abstract

:
In this paper, we present some important approximation properties of Chebyshev polynomials in the Legendre norm. We mainly discuss the Chebyshev interpolation operator at the Chebyshev–Gauss–Lobatto points. The cases of single domain and multidomain for both one dimension and multi-dimensions are considered, respectively. The approximation results in Legendre norm rather than in the Chebyshev weighted norm are given, which play a fundamental role in numerical analysis of the Legendre–Chebyshev spectral method. These results are also useful in Clenshaw–Curtis quadrature which is based on sampling the integrand at Chebyshev points.

1. Introduction

Orthogonal polynomials are useful in many areas of numerical analysis and are powerful for function approximation, numerical integration and numerical solution of differential and integral equations [1,2]. The core idea of spectral methods is that any nice enough function can be expanded in a series of orthogonal polynomials so that orthogonal polynomials play a fundamental role in spectral methods [3,4,5,6]. Particularly, Chebyshev polynomials and Legendre polynomials are frequently used in spectral methods and are two important sequences in numerical analysis.
The related approximation results of typical Chebyshev and Legendre spectral approximation are discussed in many literatures [3,4,7,8,9,10]. These results of Chebyshev spectral approximation are usually in the weighted norm forms. The Legendre–Chebyshev spectral method is a popular numerical method, which enjoys advantages of better stability of the Legendre method and easy implementation of the Chebyshev method. Therefore, it is necessary to develop the approximation properties of Chebyshev polynomials in the Legendre norm. In [11,12], the approximation result of the Chebyshev interpolation operator without the Chebyshev weighted norm was first given. Some other valuable results related to Chebyshev polynomials can be referred to [2,13,14,15,16,17,18,19] and references therein.
In addition, Chebyshev polynomials have a special connection with Clenshaw–Curtis quadrature, which uses Chebyshev points instead of optimal nodes. Clenshaw–Curtis quadrature can be implemented in O ( N log N ) operations using the fast Fourier transform (FFT) and is used in numerical integration and numerical analysis [20,21,22,23,24,25]. As we known, Gauss quadrature is a beautiful and powerful idea. Zeros of orthogonal polynomials are chosen as the nodes of Gauss-type quadratures and used to generate computational grids for spectral methods. Yet, the Clenshaw–Curtis formula has essentially the same performance for most integrands and can be implemented effortlessly by the FFT [26]. Thus, the Clenshaw–Curtis and Gauss formulas are employed in the numerical solution of Ordinary differential equations and Partial differential equations by spectral methods [5,26,27,28]. And, Chebyshev polynomials also have an important connection with the mock–Chebyshev subset interpolation exploited to cutdown the Runge phenomenon [29,30], which takes advantages of the optimality of the interpolation processes on Chebyshev–Lobatto nodes.
The purpose of this paper is to present some essential approximation results related to Chebyshev polynomials in the Legendre norm. The first fundamental result of orthogonal polynomials is the Weierstrass Theorem, which is an important element of the classical polynomial approximation theory [31,32]. In numerical analysis of the Legendre–Chebyshev spectral method, we need to consider the stability and approximation properties of the Chebyshev interpolation operator in the L 2 -norm rather than in the Chebyshev weighted norm [13]. In the paper, we consider the Chebyshev interpolation operator at the Chebyshev–Gauss–Lobatto (CGL) points. The cases of single domain and multidomain for both one dimension and multi-dimensions are discussed. Some approximation results in the Legendre norm rather than in the Chebyshev weighted norm are given. These results serve as preparations for polynomial-based spectral methods.
The rest of the paper is organized as follows. In Section 2, Chebyshev polynomials are described, and some related notations are introduced. In Section 3, some approximation properties of Chebyshev interpolation operators in one dimension are given. The cases of single domain and multidomain are discussed, respectively. In Section 4, some approximation properties in multi-dimensions are given. The conclusion is given in Section 5.

2. Preliminaries and Notations

In this section, we give a brief description of Chebyshev polynomials and define the Chebyshev interpolation operators. Some notations are also given, which will be used in the following sections.
We consider orthogonal polynomials—Chebyshev polynomials, which are proportional to Jacobi polynomials J n 1 2 , 1 2 and are orthogonal with respect to the weight function
ω ( x ) = ( 1 x 2 ) 1 2 , x [ 1 , 1 ] .
The three-term recurrence relation for the Chebyshev polynomials is as follows [6]:
T 0 ( x ) = 1 , T 1 ( x ) = x , T n + 1 ( x ) = 2 x T n ( x ) T n 1 ( x ) , n 1 ,
which satisfies
1 1 T i ( x ) T j ( x ) ( 1 x 2 ) 1 2 d x = c i 2 δ i j ,
where c 0 = 2 , c i = 1 ( i 1 ) . As we known, there have been many useful properties of Chebyshev polynomials [4,6,28].
Denote ( · , · ) Q and · Q be the inner product and the norm of the space L 2 ( Q ) , respectively. We will drop the subscript Q whenever Q = I = ( 1 , 1 ) . Let P N ( I ) be the space of polynomials with the degree at most N on an interval I. And let H σ ( Q ) ( σ > 0 ) be the classical Sobolev space with norm · H σ ( Q ) .
Define the Chebyshev interpolation operator at the CGL points by I N C : C ( I ¯ ) P N ( I ) satisfying
I N C φ ( x j ) = φ ( x j ) , 0 j N ,
where x j = cos π j N .

3. Approximation Properties of Chebyshev Interpolation Operator in One Dimension

In this section, some approximation properties of the Chebyshev interpolation operator in one dimension are derived. The cases of both single domain and multidomain are considered respectively.

3.1. Case of Single Domain in One Dimension

Similar to the approximation results presented in [11] for the Chebyshev interpolation operator I N C , we give the following lemma.
Lemma 1
([11,15]). If u H 1 ( I ) , then
N I N C u u + x I N C u C x u .
In addition, if u H σ ( I ) and σ > 1 / 2 , then
I N C u u H l ( I ) C N l σ u H σ ( I ) , 0 l 1 .
We note that the norm in the approximation results (2) and (3) is already without the Chebyshev weighted function and is in Legendre norm. The lemma is important in numerical analysis of Legendre–Chebyshev spectral method.
Next, the applications of the result of interpolation (3) to connect with the Clenshaw-Curtis quadrature are presented as follows. Given
I ( u ) = 1 1 u ( x ) d x , I N ( u ) = k = 0 N ω k u x k ,
where the nodes x k depend on N. Since the weights ω k are defined uniquely by the property that I N is equal to the integral of the degree N polynomial interpolation through the data points. Then we have
I N ( u ) = k = 0 N ω k u x k = 1 1 I N C u ( x ) d x .
For the Clenshaw-Curtis numerical integration in [26], the unique best approximation to u on [ 1 , 1 ] of degree N with respect to the L -norm.
The following lemma shows that we simply use the L 2 -norm estimation result (3) to get the desired error estimate.
Lemma 2.
If u H σ ( I ) and σ > 1 / 2 , then
| I N ( u ) I ( u ) | C N σ u H σ ( I ) .

3.2. Case of Multidomain in One Dimension

For 1 k K , we denote 1 = a 0 < a 1 < < a K = 1 , and set
I k = ( a k 1 , a k ] , I = k = 0 K I k , h k = a k a k 1 , v k v | I k .
Let P N k ( I k ) be the space of polynomials with the degree at most N k on the interval I k . Denote N = ( N 1 , , N K ) .
Define the following space
P N ( I ) = { u : u | I k P N k ( I k ) , 1 k K } .
Set the relationship between I k and I as follows:
v ( x ) = v ^ ( x ^ ) , x = 1 2 ( h k x ^ + a k 1 + a k ) , x I k , x ^ I .
Define the operator I N C : C ( I ¯ ) P N such that
( I N C u ) ^ k = I N k C u ^ k ( x ^ ) , 1 k K ,
where I N k C : C ( I ¯ ) P N k ( I ) is the CGL interpolation operator defined as (1).
Lemma 3.
If v H σ ( I k ) ( σ 0 ), then
v ^ H σ ( I ) C h k σ 1 2 v H σ ( I k ) ,
v H σ ( I k ) C h k 1 2 σ v ^ H σ ( I ) .
Denote = max 1 k K h k N k . We arrive at the following approximation result.
Theorem 1.
If u H σ ( I ) ( σ 1 ) , then
I N C u u H l ( I ) C σ l u H σ ( I ) , 0 l 1 .
Proof. 
Applying Lemma 1 and Lemma 3, we get
I N C u u H l ( I ) 2 = 1 k K ( I N C u ) k u k H l ( I k ) 2 C 1 k K h k 1 2 l I N k C u ^ k u ^ k H l ( I ) 2 C 1 k K h k 1 2 l N k 2 ( l σ ) u ^ k H σ ( I ) 2 C 1 k K h k 1 2 l N k 2 ( l σ ) h k 2 σ 1 u k H σ ( I k ) 2 C 2 ( σ l ) u H σ ( I ) 2 .
Thus, the theorem is proved. □

4. Approximation Properties of Chebyshev Interpolation Operator in Multi-Dimensions

Set I i = ( 1 , 1 ) ( i = 1 , , d ) , Ω d = I 1 × I 2 × × I d . If d = 2 , we use Ω = I x × I y instead of Ω 2 = I 1 × I 2 .
Define the following space
P N ( Ω d ) = P N ( I 1 ) P N ( I 2 ) P N ( I d ) .
Denote x = ( x 1 , , x d ) Ω d . With each function v in C ( Ω ¯ d ) , we associate the d-function v j defined by
v j ( x j ) ( x 1 , , x j 1 , x j + 1 , , x d ) = v ( x 1 , , x d ) , 1 j d .
Define the operator I N C by
I N C = I N , 1 C I N , d C ,
where I N , i C is the CGL interpolation operator I N , i C : C ( I i ¯ ) P N ( I i ) defined as (1).

4.1. Case of Single Domain in Multi-Dimensions

According to the one dimensional approximation results, we give some approximation results of the Chebyshev interpolation operators for the case of single domain in multi-dimensions ( d 2 ) .
Theorem 2.
If u H σ ( Ω d ) and σ > d 2 , then
I N C u u Ω d C N σ u H σ ( Ω d ) .
Proof. 
Applying (3) in Lemma 1 and noting L 2 ( Ω d ) = L 2 ( I j ; L 2 ( Ω d 1 ) ) , we have
u I N C u Ω d = u I N , 1 C I N , d C u Ω d u I N , 1 C u L 2 ( I 1 ; L 2 ( Ω d 1 ) ) + u I N , 2 C I N , d C u L 2 ( I 1 ; L 2 ( Ω d 1 ) ) + ( I I N , 1 C ) ( u I N , 2 C I N , d C u ) L 2 ( I 1 ; L 2 ( Ω d 1 ) ) C N σ u H σ ( I 1 ; L 2 ( Ω d 1 ) ) + u I N , 2 C I N , d C u L 2 ( I 1 ; L 2 ( Ω d 1 ) ) + C N σ d u I N , 2 C I N , d C u H σ d ( I 1 ; L 2 ( Ω d 1 ) ) .
Repeating the above discussion d times for u I N , 2 C I N , d C u L 2 ( I 1 ; L 2 ( Ω d 1 ) ) and u I N , 2 C I N , d C u H σ d ( I 1 ; L 2 ( Ω d 1 ) ) , and by the following imbedding relationship
H s ( Ω d ) H k ( I j ; H s k ( Ω d 1 ) ) , s k ,
the desired result is obtained. □
Theorem 3.
If u H σ ( Ω d ) and σ > d + 1 2 , then
I N C u u H 1 ( Ω d ) C N 1 σ u H σ ( Ω d ) .
Proof. 
For 1 j d , we have
u I N C u H 1 ( Ω d ) 2 j = 1 d u I N C H 1 ( I j ; L 2 ( Ω d 1 ) ) 2 .
By (3) in Lemma 1 and Theorem 2, we get
u I N C u H 1 ( I j ; L 2 ( Ω d 1 ) ) u I N , j C u H 1 ( I j ; L 2 ( Ω d 1 ) ) + I N , j C ( u I N , 1 C I N , j 1 C I N , j + 1 C I N , d C u ) H 1 ( I j ; L 2 ( Ω d 1 ) ) C N 1 σ u H σ ( I j ; L 2 ( Ω d 1 ) ) + C u I N , 1 C I N , j 1 C I N , j + 1 C I N , d C u H 1 ( I j ; L 2 ( Ω d 1 ) ) C N 1 σ u H σ ( I j ; L 2 ( Ω d 1 ) ) + C N 1 σ u H 1 ( I j ; H σ 1 ( Ω d 1 ) ) C N 1 σ u H σ ( Ω d ) .
Thus, the theorem is proved. □

4.2. Case of Multidomain in Multi-Dimensions

In this subsection, we give some approximation properties of the CGL interpolation operator for the case of multidomain in multi-dimensions ( d 2 ) .
For simplicity, we make the same subdivision in each direction of space. Similar to the case of multidomain in one dimension, for 1 k K , denote 1 = a 0 < a 1 < < a K = 1 , and set
I k i = ( a k 1 , a k ] , I i = k = 0 K I k i , i = 1 , , d .
Let P N k ( I k i ) be the space of polynomials with the degree at most N k on the interval I k i .
We introduce the space P N ( Ω d ) = P N ( I 1 ) P N ( I d ) . Define the Chebyshev-Gauss–Lobatto interpolation operator I N C by
I N C = I N , 1 C I N , 2 C I N , d C ,
where I N , i C : C ( I i ¯ ) P N ( I i ) is the CGL interpolation operator defined as (5).
By the assumption of the same subdivision in each direction of space, we set = max 1 k K h k N k and give the following approximation results.
Theorem 4.
Assume that d = 2 . If u H σ ( Ω ) and σ > 1 , then
I N C u u Ω C σ u H σ ( Ω ) .
Proof. 
By Theorem 1 and Lemma 1, we get
u I N C u = u I N , 1 C I N , 2 C u u I N , 1 C u L 2 ( I x ; L 2 ( I y ) ) + u I N , 2 C u L 2 ( I x ; L 2 ( I y ) ) + ( I I N , 1 C ) ( u I N , 2 C u ) L 2 ( I x ; L 2 ( I y ) ) C σ u H σ ( I x ; L 2 ( I y ) ) + u I N , 2 C u L 2 ( I x ; L 2 ( I y ) ) + C u I N , 2 C u H 1 ( I x ; L 2 ( I y ) ) C σ u H σ ( I x ; L 2 ( I y ) ) + C σ u L 2 ( I x ; H σ ( I y ) ) + C σ u H 1 ( I x ; H σ 1 ( I y ) ) C σ u H σ ( Ω ) .
Thus, the proof is completed. □
Theorem 5.
Assume that d = 2 . If u H σ ( Ω ) and σ > d + 1 2 = 3 2 , then
u I N C u H 1 ( Ω ) C σ 1 u H σ ( Ω ) .
Proof. 
By Theorem 1 and Lemma 1, we have
u I N C u H 1 ( I x ; L 2 ( I y ) ) u I N , 1 C u H 1 ( I x ; L 2 ( I y ) ) + I N , 1 C ( u I N , 2 C u ) H 1 ( I x ; L 2 ( I y ) ) C σ 1 u H σ ( I x ; L 2 ( I y ) ) + ( I I N , 1 C ) ( u I N , 2 C u ) H 1 ( I x ; L 2 ( I y ) ) + u I N , 2 C u H 1 ( I x ; L 2 ( I y ) ) C σ 1 u H σ ( I x ; L 2 ( I y ) ) + C u I N , 2 C u H 1 ( I x ; L 2 ( I y ) ) C σ 1 u H σ ( I x ; L 2 ( I y ) ) + C σ 1 u H 1 ( I x ; H σ 1 ( I y ) ) C σ 1 u H σ ( Ω ) .
Therefore, the desired result is obtained. □

5. Numerical Experiments

In this section, we give some numerical experiments to confirm the theoretical results. The cases of continuous and discontinuous functions are considered, respectively.
The discrete L 2 -error used in the following experiments is defined as
Err ( u ) = i = 0 n 1 j = 0 n 1 | u ( x i , y j ) I N C u ( x i , y j ) | 2 Δ x Δ y 1 2 ,
where x i = i Δ x , y j = j Δ y , Δ x = Δ y = 1 n , and n = 100 .
Example 1.
We consider the following continuous function in Ω ¯ = [ 0 , 1 ] 2 :
u ( x , y ) = k y w cos k x π x sin k y π y , w = ( k x 2 + k y 2 ) 1 2 ,
which is approximated by the multidomain Chebyshev-Gauss–Lobatto interpolation I N C u . We make the same subdivision in x and y directions as follows:
Ω ¯ = [ 0 , 0.5 ] [ 0.5 , 1 ] × [ 0 , 0.5 ] [ 0.5 , 1 ] .
Figure 1 displays the shape of u ( x , y ) with low frequency k x = k y = 1 and high frequency k x = k y = 5 , respectively. Table 1 gives the discrete L 2 -errors of the Chebyshev-Gauss–Lobatto interpolation for function u. The results show the spectral accuracy of the multidomain interpolation.
Example 2.
We consider the following discontinuous functions in Ω ¯ = [ 0 , 1 ] 2 :
u 1 ( x , y ) = k y ϵ w cos k x π x sin k y π y ,
u 2 ( x , y ) = k y ϵ w sin k x π x cos k y π y ,
where w = k x 2 + k y 2 ϵ 1 2 . The functions are approximated by the multidomain Chebyshev-Gauss–Lobatto interpolation I N C u .
Suppose that the parameters ϵ and k x are piecewise constants:
ϵ = 1 , 0 x 0.5 , 0 y 1 , 4 , 0.5 x 1 , 0 y 1 , k x = 4 , 0 x 0.5 , 0 y 1 , 16 , 0.5 x 1 , 0 y 1 ,
and k y = 8 . The functions are discontinuous at x = 0.5 . The domain is decomposed as follows:
Ω ¯ = [ 0 , 0.5 ] [ 0.5 , 1 ] × [ 0 , 1 ] .
Figure 2 displays the shape of u 1 ( x , y ) and u 2 ( x , y ) . It is clear that u 1 ( x , y ) is discontinuous and u 2 ( x , y ) is weak discontinuous at x = 0.5 . Table 2 gives the discrete L 2 -errors of the Chebyshev-Gauss–Lobatto interpolation for functions u 1 , u 2 . The results show the spectral accuracy of the multidomain interpolation for the discontinuous functions.

6. Conclusions

In the paper, we have given some important approximation results of Chebyshev interpolation operators in Legendre norm. The Chebyshev interpolation operators at the Chebyshev–Gauss–Lobatto points is discussed mainly. Moreover, we considered the cases of single domain and multidomain for both one dimension and multi-dimensions, respectively. The approximation results in the Legendre norm are derived. These results play an important role in numerical integration and numerical analysis of the Legendre–Chebyshev spectral method and the Clenshaw–Curtis quadrature.

Author Contributions

Conceptualization, H.M. and C.N.; Data curation, H.M. and H.L.; Formal analysis, H.M. and C.N.; Writing–original draft, C.N., H.L., H.M. and H.W.; Writing–review and editing, C.N., H.M., H.L. and H.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors sincerely thank the anonymous referees for their very helpful comments and suggestions, which are very helpful for improving the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dahlquist, G.; Björck, A.K. Numerical Methods in Scientific Computing; Society for Industrial and Applied Mathematics (SIAM): Philadelphia, PA, USA, 2008; Volume I, p. xxviii+717. [Google Scholar] [CrossRef]
  2. Xiang, S.; Chen, X.; Wang, H. Error bounds for approximation in Chebyshev points. Numer. Math. 2010, 116, 463–491. [Google Scholar] [CrossRef]
  3. Bernardi, C.; Maday, Y. Spectral methods. In Handbook of Numerical Analysis; North-Holland: Amsterdam, The Netherlands, 1997; Volume V, pp. 209–485. [Google Scholar] [CrossRef]
  4. Boyd, J.P. Chebyshev and Fourier Spectral Methods, 2nd ed.; Dover Publications Inc.: Mineola, NY, USA, 2001; p. xvi+668. [Google Scholar]
  5. Canuto, C.; Hussaini, M.Y.; Quarteroni, A.; Zang, T.A. Spectral Methods; Evolution to Complex Geometries and Applications to Fluid Dynamics; Scientific Computation; Springer: Berlin, Germany, 2007; p. xxx+596. [Google Scholar]
  6. Shen, J.; Tang, T.; Wang, L.L. Spectral Methods; Springer Series in Computational Mathematics; Algorithms, Analysis and Applications; Springer: Berlin/Heidelberg, Germany, 2011; Volume 41, p. xvi+470. [Google Scholar] [CrossRef]
  7. Maday, Y.; Quarteroni, A. Legendre and Chebyshev spectral approximations of Burgers’ equation. Numer. Math. 1981, 37, 321–332. [Google Scholar] [CrossRef]
  8. Don, W.S.; Gottlieb, D. The Chebyshev-Legendre method: Implementing Legendre methods on Chebyshev points. SIAM J. Numer. Anal. 1994, 31, 1519–1534. [Google Scholar] [CrossRef]
  9. Shen, J. Efficient Chebyshev-Legendre Galerkin methods for elliptic problems. Proc. ICOSAHOM95 Houst. J. Math. 1996, 95, 233–239. [Google Scholar]
  10. Wang, H.; Xiang, S. On the convergence rates of Legendre approximation. Math. Comp. 2012, 81, 861–877. [Google Scholar] [CrossRef] [Green Version]
  11. Ma, H. Chebyshev-Legendre spectral viscosity method for nonlinear conservation laws. SIAM J. Numer. Anal. 1998, 35, 869–892. [Google Scholar] [CrossRef] [Green Version]
  12. Ma, H. Chebyshev-Legendre super spectral viscosity method for nonlinear conservation laws. SIAM J. Numer. Anal. 1998, 35, 893–908. [Google Scholar] [CrossRef] [Green Version]
  13. Wu, H.; Ma, H.; Li, H. Optimal error estimates of the Chebyshev-Legendre spectral method for solving the generalized Burgers equation. SIAM J. Numer. Anal. 2003, 41, 659–672. [Google Scholar] [CrossRef]
  14. Li, H.; Wu, H.; Ma, H. The Legendre Galerkin-Chebyshev collocation method for Burgers-like equations. IMA J. Numer. Anal. 2003, 23, 109–124. [Google Scholar] [CrossRef]
  15. Wu, H.; Ma, H.; Li, H. Chebyshev-Legendre spectral method for solving the two-dimensional vorticity equations with homogeneous Dirichlet conditions. Numer. Methods Partial. Differ. Equ. 2009, 25, 740–755. [Google Scholar] [CrossRef]
  16. Zhao, T.; Wu, Y.; Ma, H. Error analysis of Chebyshev-Legendre pseudo-spectral method for a class of nonclassical parabolic equation. J. Sci. Comput. 2012, 52, 588–602. [Google Scholar] [CrossRef]
  17. Qin, Y.; Li, J.; Ma, H. The Legendre Galerkin Chebyshev collocation least squares for the elliptic problem. Numer. Methods Partial. Differ. Equ. 2016, 32, 1689–1703. [Google Scholar] [CrossRef]
  18. Xiang, S. On error bounds for orthogonal polynomial expansions and Gauss-type quadrature. SIAM J. Numer. Anal. 2012, 50, 1240–1263. [Google Scholar] [CrossRef]
  19. Liu, W.; Wang, L.L.; Li, H. Optimal error estimates for Chebyshev approximations of functions with limited regularity in fractional Sobolev-type spaces. Math. Comp. 2019, 88, 2857–2895. [Google Scholar] [CrossRef]
  20. Davis, P.J.; Rabinowitz, P. Methods of Numerical Integration; Computer Science and Applied Mathematics; Academic Press [Harcourt Brace Jovanovich, Publishers]: New York, NY, USA; London, UK, 1975; p. xii+459. [Google Scholar]
  21. Sloan, I.H.; Smith, W.E. Product-integration with the Clenshaw-Curtis and related points. Convergence properties. Numer. Math. 1978, 30, 415–428. [Google Scholar] [CrossRef]
  22. Johnson, L.W.; Riess, R.D. Numerical Analysis, 2nd ed.; Addison-Wesley Publishing Co.: Reading, MA, USA, 1982; p. xii+563. [Google Scholar]
  23. Mason, J.C.; Handscomb, D.C. Chebyshev Polynomials; Chapman & Hall/CRC: Boca Raton, FL, USA, 2003; p. viii+341. [Google Scholar]
  24. Neumaier, A. Introduction to Numerical Analysis; Cambridge University Press: Cambridge, UK, 2001; p. viii+356. [Google Scholar] [CrossRef] [Green Version]
  25. Kythe, P.K.; Schäferkotter, M.R. Handbook of Computational Methods for Integration; With 1 CD-ROM (Windows, Macintosh and UNIX); Chapman & Hall/CRC: Boca Raton, FL, USA, 2005; p. xxii+598. [Google Scholar]
  26. Trefethen, L.N. Is Gauss quadrature better than Clenshaw-Curtis? SIAM Rev. 2008, 50, 67–87. [Google Scholar] [CrossRef]
  27. Trefethen, L.N. Spectral Methods in MATLAB. In Software, Environments, and Tools; Society for Industrial and Applied Mathematics (SIAM): Philadelphia, PA, USA, 2000; Volume 10, p. xviii+165. [Google Scholar] [CrossRef]
  28. Canuto, C.; Hussaini, M.Y.; Quarteroni, A.; Zang, T.A. Spectral Methods; Fundamentals in Single Domains; Scientific Computation, Springer: Berlin, Germany, 2006; p. xxii+563. [Google Scholar]
  29. De Marchi, S.; Dell’Accio, F.; Mazza, M. On the constrained mock-Chebyshev least-squares. J. Comput. Appl. Math. 2015, 280, 94–109. [Google Scholar] [CrossRef]
  30. Dell’Accio, F.; Di Tommaso, F.; Nudo, F. Generalizations of the constrained mock-Chebyshev least squares in two variables: Tensor product vs total degree polynomial interpolation. Appl. Math. Lett. 2022, 125, 107732. [Google Scholar] [CrossRef]
  31. Powell, M.J.D. Approximation Theory and Methods; Cambridge University Press: Cambridge, UK; New York, NY, USA, 1981; p. x+339. [Google Scholar]
  32. Trefethen, L.N. Approximation Theory and Approximation Practice; Society for Industrial and Applied Mathematics (SIAM): Philadelphia, PA, USA, 2013; p. viii+305. [Google Scholar]
Figure 1. The shapes of u ( x , y ) with low frequency k x = k y = 1 and high frequency k x = k y = 5 .
Figure 1. The shapes of u ( x , y ) with low frequency k x = k y = 1 and high frequency k x = k y = 5 .
Mathematics 09 03271 g001
Figure 2. The shapes of u 1 ( x , y ) and u 2 ( x , y ) .
Figure 2. The shapes of u 1 ( x , y ) and u 2 ( x , y ) .
Mathematics 09 03271 g002
Table 1. L 2 -errors of Chebyshev-Gauss–Lobatto interpolation for u ( x , y ) .
Table 1. L 2 -errors of Chebyshev-Gauss–Lobatto interpolation for u ( x , y ) .
k x = k y = 1 k x = k y = 5
N x = N y Err ( u ) Order N x = N y Err ( u ) Order
(7,7) 2.15 × 10 8 -(14,14) 1.76 × 10 8 -
(10,10) 1.28 × 10 12 27.28 (18,18) 2.93 × 10 12 34.62
(13,13) 1.61 × 10 16 34.23 (22,22) 6.31 × 10 16 42.07
Table 2. L 2 -errors of Chebyshev-Gauss–Lobatto interpolation for u 1 ( x , y ) and u 2 ( x , y ) .
Table 2. L 2 -errors of Chebyshev-Gauss–Lobatto interpolation for u 1 ( x , y ) and u 2 ( x , y ) .
N x N y Err ( u 1 ) Order Err ( u 2 ) Order
(12,26)26 8.30 × 10 8 - 4.87 × 10 8 -
(16,32)32 8.32 × 10 12 32.00 4.99 × 10 12 31.93
(20,38)38 7.52 × 10 16 54.18 7.74 × 10 16 51.04
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Niu, C.; Liao, H.; Ma, H.; Wu, H. Approximation Properties of Chebyshev Polynomials in the Legendre Norm. Mathematics 2021, 9, 3271. https://doi.org/10.3390/math9243271

AMA Style

Niu C, Liao H, Ma H, Wu H. Approximation Properties of Chebyshev Polynomials in the Legendre Norm. Mathematics. 2021; 9(24):3271. https://doi.org/10.3390/math9243271

Chicago/Turabian Style

Niu, Cuixia, Huiqing Liao, Heping Ma, and Hua Wu. 2021. "Approximation Properties of Chebyshev Polynomials in the Legendre Norm" Mathematics 9, no. 24: 3271. https://doi.org/10.3390/math9243271

APA Style

Niu, C., Liao, H., Ma, H., & Wu, H. (2021). Approximation Properties of Chebyshev Polynomials in the Legendre Norm. Mathematics, 9(24), 3271. https://doi.org/10.3390/math9243271

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop