The resolution of the Gibbs phenomenon for "spliced" functions in one and two dimensions

Anne Gelb, D. Gottlieb

Research output: Contribution to journalArticle

25 Citations (Scopus)

Abstract

In this paper we study approximation methods for analytic functions that have been "spliced" into nonintersecting subdomains. We assume that we are given the first 2N + 1 Fourier coefficients for the functions in each subdomain. The objective is to approximate the "spliced" function in each subdomain and then to "glue" the approximations together in order to recover the original function in the full domain. The Fourier partial sum approximation in each subdomain yields poor results, as the convergence is slow and spurious oscillations occur at the boundaries of each subdomain. Thus once we "glue" the subdomain approximations back together, the approximation for the function in the full domain will exhibit oscillations throughout the entire domain. Recently methods have been developed that successfully eliminate the Gibbs phenomenon for analytic but nonperiodic functions in one dimension. These methods are based on the knowledge of the first 2N + 1 Fourier coefficients and use either the Gegenbauer polynomials (Gottlieb et al.) or the Bernoulli polynomials (Abarbanel, Gottlieb, Cai et al., and Eckhoff). We propose a way to accurately reconstruct a "spliced" function in a full domain by extending the current methods to eliminate the Gibbs phenomenon in each nonintersecting subdomain and then "gluing" the approximations back together. We solve this problem in both one and two dimensions. In the one-dimensional case we provide two alternative options, the Bernoulli method and the Gegenbauer method, as well as a new hybrid method, the Gegenbauer-Bernoulli method. In the two-dimensional case we prove, for the very first time, exponential convergence of the Gegenbauer method, and then we apply it to solve the "spliced" function problem.

Original languageEnglish (US)
Pages (from-to)35-58
Number of pages24
JournalComputers and Mathematics with Applications
Volume33
Issue number11
StatePublished - Jun 1997
Externally publishedYes

Fingerprint

Gibbs Phenomenon
One Dimension
Two Dimensions
Approximation
Fourier coefficients
Glues
Bernoulli
Eliminate
Oscillation
Gegenbauer Polynomials
Polynomials
Bernoulli Polynomials
Gluing
Exponential Convergence
Partial Sums
Hybrid Method
Approximation Methods
Analytic function
Entire
Alternatives

Keywords

  • Bernoulli polynomials
  • Exponential accuracy
  • Fourier series
  • Gegenbauer polynomials
  • Gibbs phenomenon

ASJC Scopus subject areas

  • Applied Mathematics
  • Computational Mathematics
  • Modeling and Simulation

Cite this

The resolution of the Gibbs phenomenon for "spliced" functions in one and two dimensions. / Gelb, Anne; Gottlieb, D.

In: Computers and Mathematics with Applications, Vol. 33, No. 11, 06.1997, p. 35-58.

Research output: Contribution to journalArticle

@article{c3311b1343974d20bd7828876b6eb745,
title = "The resolution of the Gibbs phenomenon for {"}spliced{"} functions in one and two dimensions",
abstract = "In this paper we study approximation methods for analytic functions that have been {"}spliced{"} into nonintersecting subdomains. We assume that we are given the first 2N + 1 Fourier coefficients for the functions in each subdomain. The objective is to approximate the {"}spliced{"} function in each subdomain and then to {"}glue{"} the approximations together in order to recover the original function in the full domain. The Fourier partial sum approximation in each subdomain yields poor results, as the convergence is slow and spurious oscillations occur at the boundaries of each subdomain. Thus once we {"}glue{"} the subdomain approximations back together, the approximation for the function in the full domain will exhibit oscillations throughout the entire domain. Recently methods have been developed that successfully eliminate the Gibbs phenomenon for analytic but nonperiodic functions in one dimension. These methods are based on the knowledge of the first 2N + 1 Fourier coefficients and use either the Gegenbauer polynomials (Gottlieb et al.) or the Bernoulli polynomials (Abarbanel, Gottlieb, Cai et al., and Eckhoff). We propose a way to accurately reconstruct a {"}spliced{"} function in a full domain by extending the current methods to eliminate the Gibbs phenomenon in each nonintersecting subdomain and then {"}gluing{"} the approximations back together. We solve this problem in both one and two dimensions. In the one-dimensional case we provide two alternative options, the Bernoulli method and the Gegenbauer method, as well as a new hybrid method, the Gegenbauer-Bernoulli method. In the two-dimensional case we prove, for the very first time, exponential convergence of the Gegenbauer method, and then we apply it to solve the {"}spliced{"} function problem.",
keywords = "Bernoulli polynomials, Exponential accuracy, Fourier series, Gegenbauer polynomials, Gibbs phenomenon",
author = "Anne Gelb and D. Gottlieb",
year = "1997",
month = "6",
language = "English (US)",
volume = "33",
pages = "35--58",
journal = "Computers and Mathematics with Applications",
issn = "0898-1221",
publisher = "Elsevier Limited",
number = "11",

}

TY - JOUR

T1 - The resolution of the Gibbs phenomenon for "spliced" functions in one and two dimensions

AU - Gelb, Anne

AU - Gottlieb, D.

PY - 1997/6

Y1 - 1997/6

N2 - In this paper we study approximation methods for analytic functions that have been "spliced" into nonintersecting subdomains. We assume that we are given the first 2N + 1 Fourier coefficients for the functions in each subdomain. The objective is to approximate the "spliced" function in each subdomain and then to "glue" the approximations together in order to recover the original function in the full domain. The Fourier partial sum approximation in each subdomain yields poor results, as the convergence is slow and spurious oscillations occur at the boundaries of each subdomain. Thus once we "glue" the subdomain approximations back together, the approximation for the function in the full domain will exhibit oscillations throughout the entire domain. Recently methods have been developed that successfully eliminate the Gibbs phenomenon for analytic but nonperiodic functions in one dimension. These methods are based on the knowledge of the first 2N + 1 Fourier coefficients and use either the Gegenbauer polynomials (Gottlieb et al.) or the Bernoulli polynomials (Abarbanel, Gottlieb, Cai et al., and Eckhoff). We propose a way to accurately reconstruct a "spliced" function in a full domain by extending the current methods to eliminate the Gibbs phenomenon in each nonintersecting subdomain and then "gluing" the approximations back together. We solve this problem in both one and two dimensions. In the one-dimensional case we provide two alternative options, the Bernoulli method and the Gegenbauer method, as well as a new hybrid method, the Gegenbauer-Bernoulli method. In the two-dimensional case we prove, for the very first time, exponential convergence of the Gegenbauer method, and then we apply it to solve the "spliced" function problem.

AB - In this paper we study approximation methods for analytic functions that have been "spliced" into nonintersecting subdomains. We assume that we are given the first 2N + 1 Fourier coefficients for the functions in each subdomain. The objective is to approximate the "spliced" function in each subdomain and then to "glue" the approximations together in order to recover the original function in the full domain. The Fourier partial sum approximation in each subdomain yields poor results, as the convergence is slow and spurious oscillations occur at the boundaries of each subdomain. Thus once we "glue" the subdomain approximations back together, the approximation for the function in the full domain will exhibit oscillations throughout the entire domain. Recently methods have been developed that successfully eliminate the Gibbs phenomenon for analytic but nonperiodic functions in one dimension. These methods are based on the knowledge of the first 2N + 1 Fourier coefficients and use either the Gegenbauer polynomials (Gottlieb et al.) or the Bernoulli polynomials (Abarbanel, Gottlieb, Cai et al., and Eckhoff). We propose a way to accurately reconstruct a "spliced" function in a full domain by extending the current methods to eliminate the Gibbs phenomenon in each nonintersecting subdomain and then "gluing" the approximations back together. We solve this problem in both one and two dimensions. In the one-dimensional case we provide two alternative options, the Bernoulli method and the Gegenbauer method, as well as a new hybrid method, the Gegenbauer-Bernoulli method. In the two-dimensional case we prove, for the very first time, exponential convergence of the Gegenbauer method, and then we apply it to solve the "spliced" function problem.

KW - Bernoulli polynomials

KW - Exponential accuracy

KW - Fourier series

KW - Gegenbauer polynomials

KW - Gibbs phenomenon

UR - http://www.scopus.com/inward/record.url?scp=0031167838&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0031167838&partnerID=8YFLogxK

M3 - Article

AN - SCOPUS:0031167838

VL - 33

SP - 35

EP - 58

JO - Computers and Mathematics with Applications

JF - Computers and Mathematics with Applications

SN - 0898-1221

IS - 11

ER -