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
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
-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
and are orthogonal with respect to the weight function
The three-term recurrence relation for the Chebyshev polynomials is as follows [
6]:
which satisfies
where
,
. As we known, there have been many useful properties of Chebyshev polynomials [
4,
6,
28].
Denote and be the inner product and the norm of the space , respectively. We will drop the subscript Q whenever . Let be the space of polynomials with the degree at most N on an interval I. And let be the classical Sobolev space with norm .
Define the Chebyshev interpolation operator at the CGL points by
satisfying
where
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
, we give the following lemma.
Lemma 1 ([
11,
15]).
If , thenIn addition, if and , then 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
where the nodes
depend on
N. Since the weights
are defined uniquely by the property that
is equal to the integral of the degree
polynomial interpolation through the data points. Then we have
For the Clenshaw-Curtis numerical integration in [
26], the unique best approximation to
u on
of degree
with respect to the
-norm.
The following lemma shows that we simply use the
-norm estimation result (
3) to get the desired error estimate.
Lemma 2. If and , then 3.2. Case of Multidomain in One Dimension
For
, we denote
and set
Let be the space of polynomials with the degree at most on the interval . Denote
Define the following space
Set the relationship between
and
I as follows:
Define the operator
such that
where
is the CGL interpolation operator defined as (
1).
Lemma 3. If (), then Denote We arrive at the following approximation result.
Theorem 1. If then Proof. Applying Lemma 1 and Lemma 3, we get
Thus, the theorem is proved. □
4. Approximation Properties of Chebyshev Interpolation Operator in Multi-Dimensions
Set If , we use instead of
Define the following space
Denote
. With each function
v in
, we associate the
d-function
defined by
Define the operator
by
where
is the CGL interpolation operator
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 .
Theorem 2. If and then Proof. Applying (
3) in Lemma 1 and noting
we have
Repeating the above discussion
d times for
and
and by the following imbedding relationship
the desired result is obtained. □
Theorem 3. If and then Proof. By (
3) in Lemma 1 and Theorem 2, we get
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 .
For simplicity, we make the same subdivision in each direction of space. Similar to the case of multidomain in one dimension, for
, denote
and set
Let be the space of polynomials with the degree at most on the interval .
We introduce the space
Define the Chebyshev-Gauss–Lobatto interpolation operator
by
where
is the CGL interpolation operator defined as (
5).
By the assumption of the same subdivision in each direction of space, we set and give the following approximation results.
Theorem 4. Assume that . If and then Proof. By Theorem 1 and Lemma 1, we get
Thus, the proof is completed. □
Theorem 5. Assume that . If and then Proof. By Theorem 1 and Lemma 1, we have
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
-error used in the following experiments is defined as
where
,
,
, and
.
Example 1. We consider the following continuous function in
:
which is approximated by the multidomain Chebyshev-Gauss–Lobatto interpolation
. We make the same subdivision in
x and
y directions as follows:
Figure 1 displays the shape of
with low frequency
and high frequency
, respectively.
Table 1 gives the discrete
-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 : where . The functions are approximated by the multidomain Chebyshev-Gauss–Lobatto interpolation .
Suppose that the parameters ϵ and are piecewise constants: and . The functions are discontinuous at . The domain is decomposed as follows: Figure 2 displays the shape of and . It is clear that is discontinuous and is weak discontinuous at . Table 2 gives the discrete -errors of the Chebyshev-Gauss–Lobatto interpolation for functions . 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
- 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]
- Xiang, S.; Chen, X.; Wang, H. Error bounds for approximation in Chebyshev points. Numer. Math. 2010, 116, 463–491. [Google Scholar] [CrossRef]
- 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]
- Boyd, J.P. Chebyshev and Fourier Spectral Methods, 2nd ed.; Dover Publications Inc.: Mineola, NY, USA, 2001; p. xvi+668. [Google Scholar]
- 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]
- 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]
- Maday, Y.; Quarteroni, A. Legendre and Chebyshev spectral approximations of Burgers’ equation. Numer. Math. 1981, 37, 321–332. [Google Scholar] [CrossRef]
- 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]
- Shen, J. Efficient Chebyshev-Legendre Galerkin methods for elliptic problems. Proc. ICOSAHOM95 Houst. J. Math. 1996, 95, 233–239. [Google Scholar]
- Wang, H.; Xiang, S. On the convergence rates of Legendre approximation. Math. Comp. 2012, 81, 861–877. [Google Scholar] [CrossRef] [Green Version]
- Ma, H. Chebyshev-Legendre spectral viscosity method for nonlinear conservation laws. SIAM J. Numer. Anal. 1998, 35, 869–892. [Google Scholar] [CrossRef] [Green Version]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- Xiang, S. On error bounds for orthogonal polynomial expansions and Gauss-type quadrature. SIAM J. Numer. Anal. 2012, 50, 1240–1263. [Google Scholar] [CrossRef]
- 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]
- 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]
- 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]
- Johnson, L.W.; Riess, R.D. Numerical Analysis, 2nd ed.; Addison-Wesley Publishing Co.: Reading, MA, USA, 1982; p. xii+563. [Google Scholar]
- Mason, J.C.; Handscomb, D.C. Chebyshev Polynomials; Chapman & Hall/CRC: Boca Raton, FL, USA, 2003; p. viii+341. [Google Scholar]
- Neumaier, A. Introduction to Numerical Analysis; Cambridge University Press: Cambridge, UK, 2001; p. viii+356. [Google Scholar] [CrossRef] [Green Version]
- 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]
- Trefethen, L.N. Is Gauss quadrature better than Clenshaw-Curtis? SIAM Rev. 2008, 50, 67–87. [Google Scholar] [CrossRef]
- 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]
- 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]
- 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]
- 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]
- Powell, M.J.D. Approximation Theory and Methods; Cambridge University Press: Cambridge, UK; New York, NY, USA, 1981; p. x+339. [Google Scholar]
- Trefethen, L.N. Approximation Theory and Approximation Practice; Society for Industrial and Applied Mathematics (SIAM): Philadelphia, PA, USA, 2013; p. viii+305. [Google Scholar]
| Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).