Research & Teaching Faculty

Towards the resolution of the Gibbs phenomena

TitleTowards the resolution of the Gibbs phenomena
Publication TypeJournal Article
Year of Publication2003
AuthorsShizgal, BD, Jung, JH
JournalJournal of Computational and Applied Mathematics
Volume161
Pagination41-65
Date PublishedDec
Type of ArticleArticle
ISBN Number0377-0427
Keywordsequation, Gibbs phenomenon, spectral methods
Abstract

It is well known that the expansion of an analytic nonperiodic function on a finite interval in a Fourier series leads to spurious oscillations at the interval boundaries. This result is known as the Gibbs phenomenon. The present paper introduces a new method for the resolution of the Gibbs phenomenon which follows on the reconstruction method of Gottlieb and coworkers (SIAM Rev. 39 (1997) 644) based on Gegenbauer polynomials orthogonal with respect to weight function (1-x(2))(lambda-1/2). We refer to their approach as the direct method and to the new methodology as the inverse method. Both methods use the finite set of Fourier coefficients of some given function as input data in the re-expansion of the function in Gegenbauer polynomials or in other orthogonal basis sets. The finite partial sum of the new expansion provides a spectrally accurate approximation to the function. In the direct method, this requires that certain conditions are met concerning the parameter lambda in the weight function, the number of Fourier coefficients, N and the number of Gegenbauer polynomials, m. We show that the new inverse method can give exact results for polynomials independent of lambda and with m=N. The paper presents several numerical examples applied to a single domain or to subdomains of the main domain so as to illustrate the superiority of the inverse method in comparison with the direct method. (C) 2003 Elsevier B.V. All rights reserved.

URL<Go to ISI>://000186545000003