Asymptotic Coefficients and Errors for Chebyshev Polynomial Approximations with Weak Endpoint Singularities: Effects of Different Bases
Abstract
When solving differential equations by a spectral method, it is often convenient to shift from Chebyshev polynomials Tn(x) with coefficients an to modified basis functions that incorporate the boundary conditions. For homogeneous Dirichlet boundary conditions, u( 1)=0, popular choices include the ``Chebyshev difference basis", n(x) Tn+2(x) - Tn(x) with coefficients here denoted bn and the ``quadratic-factor basis functions" n(x) (1-x2) Tn(x) with coefficients cn. If u(x) is weakly singular at the boundaries, then an will decrease proportionally to O(A(n)/n) for some positive constant , where the A(n) is a logarithm or a constant. We prove that the Chebyshev difference coefficients bn decrease more slowly by a factor of 1/n while the quadratic-factor coefficients cn decrease more slowly still as O(A(n)/n-2). The error for the unconstrained Chebyshev series, truncated at degree n=N, is O(|A(N)|/N) in the interior, but is worse by one power of N in narrow boundary layers near each of the endpoints. Despite having nearly identical error norms, the error in the Chebyshev basis is concentrated in boundary layers near both endpoints, whereas the error in the quadratic-factor and difference basis sets is nearly uniform oscillations over the entire interval in x. Meanwhile, for Chebyshev polynomials and the quadratic-factor basis, the value of the derivatives at the endpoints is O(N2), but only O(N) for the difference basis.
Turn this paper into a lesson
ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.