---
title: Chebyshev Interpolation
url: https://www.emergentmind.com/topics/chebyshev-interpolation
type: topic
---

# Chebyshev Interpolation

Chebyshev interpolation is a polynomial interpolation scheme based on evaluation at specially chosen, non-uniform Chebyshev nodes, resulting in near-minimax uniform approximation and stable numerical properties. The method utilizes the structure of Chebyshev polynomials and their extremal distribution of nodes to suppress the Runge phenomenon and achieve exponential convergence for analytic functions. Chebyshev interpolation generalizes efficiently to higher dimensions and serves as the foundation for numerous modern algorithms in spectral methods, uncertainty quantification, numerical analysis, and computational finance. This article details the mathematical structure, algorithmic implementation, theoretical guarantees, practical extensions, and modern research related to Chebyshev interpolation.

## 1. Mathematical Foundations and Node Selection

Chebyshev interpolation exploits the properties of Chebyshev polynomials of the first kind, defined for $x \in [-1, 1]$ as $T_n(x) = \cos(n\,\arccos x)$. These polynomials have explicit extremal and optimality structures on $[-1,1]$:
- The zeros of $T_{n+1}$ give the Chebyshev–Gauss nodes, while the extrema provide the Chebyshev–Gauss–Lobatto nodes.
- For degree $N$, Chebyshev–Gauss nodes are $x_k = \cos\left(\frac{2k+1}{2N+2}\pi\right)$, $k=0,1,\dots,N$ [2305.04845].
- Chebyshev nodes cluster quadratically near the endpoints, a property central for controlling large oscillations in high-degree interpolation (Runge phenomenon) [2404.00414].

The optimality of Chebyshev or scaled Chebyshev grids follows from classical minimax theory; among all node distributions interpolating in $C_M^{s+1}[-1,1]$, scaled Chebyshev nodes minimize the max-norm of the error factor $\max_{x \in[-1,1]}|\prod_i(x-c_i)|$ [1305.6104].

## 2. Polynomial Construction and Algorithms

Given function samples $f_k$ at Chebyshev nodes $x_k$, the interpolating polynomial $P_N(x)$ can be expressed in several mathematically equivalent but numerically distinct forms:
- **Lagrange form**: $P_N(x) = \sum_{k=0}^N f_k \ell_k(x)$, where $\ell_k(x)$ is the fundamental Lagrange basis [1505.04648].
- **Chebyshev (cosine) expansion**: $P_N(x) = \sum_{k=0}^N a_k T_k(x)$, with coefficients computed via a discrete cosine transform (DCT) [2404.00414].
- **Barycentric form**: The numerically preferred form,
  $$
  P_N(x) = \frac{\sum_{k=0}^N w_k f_k / (x - x_k)}{\sum_{k=0}^N w_k / (x - x_k)}
  $$
  with barycentric weights $w_k = (-1)^k \sin\left(\frac{2k+1}{2N+2}\pi\right)$ for Gauss nodes [2305.04845, 1805.00898]. Evaluation is $O(N)$, insensitive to rounding error, and vectorizable in modern implementations.

Key features:
- At node $x_j$, $P_N(x_j) = f_j$ (Lagrange property).
- If $x$ is very close to a node, direct evaluation with $P_N(x_j)$ avoids cancellation.

FFT-accelerated algorithms and DCTs enable $O(N \log N)$ coefficient computation for large-scale problems [1001.1185, 2404.00414].

## 3. Error Analysis and Convergence

The error of Chebyshev interpolation admits both classical and modern bounds:
- **Smoothness-based (Sobolev) bound**: If $f \in C^r[-1,1]$, $E_N := \max_{|x| \le 1}|f(x) - P_N(x)| \le C N^{-r} \max_{|x| \le 1}|f^{(r)}(x)|$ [2305.04845].
- **Analyticity-based (Bernstein ellipse) bound**: For $f$ analytic in the open Bernstein ellipse of parameter $\rho > 1$,
  $$
  E_N \le 2 M \frac{\rho^{-N}}{\rho-1}
  $$
  with $M = \max_{z \in \text{ellipse}}|f(z)|$, yielding geometric (spectral) convergence [1805.00898, 1611.08706, 2404.00414, 1505.04648].

The pointwise interpolation error is near-minimax: for the Chebyshev grid, the interpolant $P_N$ is within a fixed constant of the best uniform approximation polynomial of degree $N$ [1305.6104].

Lebesgue constants, which govern the maximal error amplification from nodal noise, grow logarithmically for Chebyshev grids: $\Lambda_N = \frac{2}{\pi}\log N + O(1)$, as opposed to exponential growth for equispaced nodes [2404.00414, 2511.17911].

## 4. Multidimensional Tensorization and Generalizations

Chebyshev interpolation extends naturally to multivariate settings via tensor-product constructions:
- For $d$ dimensions, nodes are defined as $x_{j_i}^{(i)} = \cos\left(\frac{2j_i+1}{2(N_i+1)}\pi\right)$, $j_i=0,\dots,N_i$ for each dimension [2305.04845].
- The full tensor grid comprises $\prod_{i=1}^d (N_i+1)$ nodes.
- The interpolant is expressed as a sum over tensor-product basis functions: $P(x_1,\ldots,x_d) = \sum_{j_1,\ldots,j_d} c_{j_1,\ldots,j_d} \prod_{i=1}^d T_{j_i}(x_i)$ [1505.04648, 1611.08706].

Error bounds for the multivariate case generalize the 1D result. If $f$ is analytic in a poly-ellipse with radii $\rho_i$, then
$$
\|f - I_{N_1,\ldots,N_d} f\|_{L^\infty} \le C\, \max_i \rho_i^{-N_i}
$$
and sharper mixed-derivative bounds (see [1611.08706] for advanced technical refinements).

High-dimensional tensor grids face the curse of dimensionality; modern approaches use adaptive, sparse, or low-rank tensor decompositions (e.g., functional Tucker expansion, as in [2007.16126]), and alternative node sets such as Lissajous–Chebyshev or Padua points to balance complexity and accuracy [1511.04564, 2001.09144].

## 5. Stability, Filtering, and Gibbs Phenomenon

Chebyshev interpolation at classical nodes avoids the Runge phenomenon and is stable in both the uniform and weighted norms. However, classical Lagrange interpolation at Chebyshev nodes still incurs Lebesgue constants $\sim \log N$ and suffers from Gibbs oscillations near jump discontinuities.

Filtered polynomial interpolation, notably through de la Vallée–Poussin (VP) operators, combines local Lagrange interpolation with spectral filtering:
- The VP interpolant $V_n^m$ at Chebyshev nodes applies filter coefficients $\mu_{n,j}^m$ to damp high-degree modes while preserving interpolation at the nodes [2008.00240, 2101.04551].
- VP interpolation uniformly bounds the Lebesgue constant—no logarithmic growth—provided Jacobi weight exponents satisfy explicit inequalities [2008.00240].
- Filtered Chebyshev interpolants exhibit strong suppression of the Gibbs phenomenon, preserving spectral convergence in smooth pieces and rapidly damping overshoots at discontinuities.

In the context of weighted spaces (Jacobi weights), necessary and sufficient conditions characterize when the VP operator’s Lebesgue constant is uniformly bounded, ensuring near-best approximation order, i.e., $\|f - V_n^m[f]\|_\infty \le C E_{n-m}(f)$, where $E_n(f)$ is the best approximation error by degree-$n$ polynomials [2008.00240].

## 6. Applications and Adaptive Strategies

Chebyshev interpolation underpins a diverse array of applications:
- **Numerical PDEs and spectral methods**: Used for constructing differentiation matrices, pseudospectral solvers, and projection operators in both regular and singularly perturbed settings [1708.02611].
- **Sparse and compressed reconstruction**: Sparse Chebyshev interpolation recovers functions with few nonzero Chebyshev coefficients with sample complexity $O(r^2 \log^{2d-2} r)$ in the multivariate setting, dramatically reducing the number of required samples relative to tensor grids [2001.09144, 1912.05719].
- **Signal processing**: Chebyshev interpolation's minimax and cosine-series structure enables stable signal reconstruction from nonequispaced samples, outperforming equispaced and even some Fourier methods in presence of endpoint singularities or Gibbs artifacts [2404.00414, 1001.1185].
- **Derivative computation (“Greeks”) in finance**: Once built, the derivative of a Chebyshev interpolant is available in closed form, allowing ultra-efficient and stable sensitivity calculations for option pricing and risk metrics [1805.00898, 1710.01797].
- **Algorithm architectures**: Digital and hardware-efficient architectures, such as systolic and word-serial DCT engines, implement Chebyshev interpolation efficiently for real-time systems, often with major reductions in sampling rate and power consumption compared to equispaced approaches [1001.1185].

Adaptive partition and overlapping partition of unity methods resolve non-analyticities and local singularities by locally splitting the interval and gluing together local Chebyshev interpolants with smooth weights, preserving spectral convergence globally and reducing node counts by orders of magnitude compared to a global interpolant [1708.02611].

## 7. Modern Developments, Generalizations, and Stability

Recent research addresses both theoretical and practical aspects:
- **Stable polynomial interpolation on nonequispaced nodes**: Grünwald-type operators and symmetric-wave interpolation establish strong convergence and stability even for certain perturbed or equidistant node sets, provided filtering or node warping is applied [2511.16927, 2511.17911].
- **Generalized Chebyshev bases**: Extensions to generalized Chebyshev polynomials associated with root systems and Weyl groups facilitate sparse interpolation for high-dimensional function spaces with algebraic or symmetry constraints [2001.09144, 1511.04564].
- **Error bounds in high dimensions**: Recent sharp bounds (Glau & Mahlstedt [1611.08706]) provide tight, dimension-adapted error estimates for tensor-product Chebyshev interpolation.
- **Practical guidelines**: For high-accuracy or high-dimensional tasks, careful selection of node distribution (e.g., scaled Chebyshev, adaptive partitioning) and the use of filtering (VP or similar) optimizes both approximation order and computational stability [1305.6104, 2008.00240, 2101.04551].

Filtered and adaptively constructed Chebyshev interpolants, as well as low-rank tensorized surrogates, remain the methods of choice whenever uniform accuracy, high stability, and computational efficiency are critical.

---

**References**: All results and constructions detailed above are grounded in the cited arXiv works [2305.04845], [1305.6104], [2008.00240], [2101.04551], [1611.08706], [2001.09144], [2511.17911], [2511.16927], [1511.04564], [1805.00898], [1505.04648], [1710.01797], [1001.1185], [2404.00414], [1708.02611], [2007.16126], [1810.04282], [1912.05719], [1605.02776].

Source: https://www.emergentmind.com/topics/chebyshev-interpolation