Robust reprojection methods for the resolution of the Gibbs phenomenon

Anne Gelb, Jared Tanner

Research output: Contribution to journalArticlepeer-review

66 Scopus citations

Abstract

The classical Gibbs phenomenon exhibited by global Fourier projections and interpolants can be resolved in smooth regions by reprojecting in a truncated Gegenbauer series, achieving high resolution recovery of the function up to the point of discontinuity. Unfortunately, due to the poor conditioning of the Gegenbauer polynomials, the method suffers both from numerical round-off error and the Runge phenomenon. In some cases the method fails to converge. Following the work in [D. Gottlieb, C.W. Shu, Atti Conv. Lincei 147 (1998) 39-48], a more general framework for reprojection methods is introduced here. From this insight we propose an additional requirement on the reprojection basis which ameliorates the limitations of the Gegenbauer reconstruction. The new robust Gibbs complementary basis yields a reliable exponentially accurate resolution of the Gibbs phenomenon up to the discontinuities.

Original languageEnglish (US)
Pages (from-to)3-25
Number of pages23
JournalApplied and Computational Harmonic Analysis
Volume20
Issue number1
DOIs
StatePublished - Jan 2006

Keywords

  • Fourier series
  • Freud polynomials
  • Gegenbauer post-processing
  • Robust Gibbs complementary
  • Round-off error
  • Weight functions

ASJC Scopus subject areas

  • Applied Mathematics

Fingerprint

Dive into the research topics of 'Robust reprojection methods for the resolution of the Gibbs phenomenon'. Together they form a unique fingerprint.

Cite this