Gibbs phenomenon

Gibbs phenomenon

In mathematics, the Gibbs phenomenon (also known as ringing artifacts), named after the American physicist J. Willard Gibbs, is the peculiar manner in which the Fourier series of a piecewise continuously differentiable periodic function "f" behaves at a jump discontinuity: the "n"th partial sum of the Fourier series has large oscillations near the jump, which might increase the maximum of the partial sum above that of the function itself. The overshoot does not die out as the frequency increases, but approaches a finite limit.cite book
author=H S Carslaw
title=Introduction to the theory of Fourier's series and integrals
page=Chapter IX
year= 1930
edition=Third Edition
publisher=Dover Publications Inc.
location=New York
url=http://books.google.com/books?id=JNVAAAAAIAAJ&printsec=frontcover&dq=intitle:Introduction+intitle:to+intitle:the+intitle:theory+intitle:of+intitle:Fourier%27s+intitle:series+intitle:and+intitle:integrals+inauthor:carslaw&lr=&as_brr=0#PPA264,M1
]

Description

The three pictures on the right demonstrate the phenomenon for a square wave whose Fourier expansion is

: sin(x)+frac{1}{3}sin(3x)+frac{1}{5}sin(5x)+dotsb

More precisely, this is the function "f" which equals pi/4 between 2npi and (2n+1)pi and -pi/4 between (2n+1)pi and (2n+2)pi for every integer "n"; thus this square wave has a jump discontinuity of height pi/2 at every integer multiple of pi.

As can be seen, as the number of terms rises, the error of the approximation is reduced in width and energy, but converges to a fixed height. A calculation for the square wave (see Zygmund, chap. 8.5., or the computations at the end of this article) gives an explicit formula for the limit of the height of the error. It turns out that the Fourier series exceeds the height pi/4 of the square wave by

:frac{1}{2}int_0^pi frac{sin t}{t}, dt - frac{pi}{4} = frac{pi}{2}cdot (0.089490dots)

or about 9 percent. More generally, at any jump point of a piecewise continuously differentiable function with a jump of "a", the "n"th partial Fourier series will (for "n" very large) overshoot this jump by approximately a cdot (0.089490dots) at one end and undershoot it by the same amount at the other end; thus the "jump" in the partial Fourier series will be about 18% larger than the jump in the original function. At the location of the discontinuity itself, the partial Fourier series will converge to the midpoint of the jump (regardless of what the actual value of the original function is at this point). The quantity:int_0^pi frac{sin t}{t} dt = (1.851937052dots) = frac{pi}{2} + pi cdot (0.089490dots)is sometimes known as the "Wilbraham-Gibbs constant".

The Gibbs phenomenon was first noticed and analyzed by the obscure Henry Wilbraham. He published a paper on it in 1848 that went unnoticed by the mathematical world. It was not until Albert Michelson observed the phenomenon via a mechanical graphing machine that interest arose. Michelson developed a device in 1898 that could compute and re-synthesize the Fourier series. When the Fourier coefficients for a square wave were input to the machine, the graph would oscillate at the discontinuities. This would continue to occur even as the number of Fourier coefficients increased.

Because it was a physical device subject to manufacturing flaws, Michelson was convinced that the overshoot was caused by errors in the machine. In 1898 J. Willard Gibbs published a paper on Fourier series in which he discussed the example of what today would be called a "sawtooth wave", and described the graph obtained as a limit of the graphs of the partial sums of the Fourier series. Interestingly in this paper he failed to notice the phenomenon that bears his name, and the limit he described was incorrect. In 1899 he published a correction to his paper in which he describes the phenomenon and points out the important distinction between the limit of the graphs and the graph of the function that is the limit of the partial sums of the Fourier series. Maxime Bôcher gave a detailed mathematical analysis of the phenomenon in 1906 and named it the Gibbs phenomenon.

Informally, it reflects the difficulty inherent in approximating a discontinuous function by a "finite" series of continuous sine and cosine waves. It is important to put emphasis on the word "finite" because even though every partial sum of the Fourier series overshoots the function it is approximating, the limit of the partial sums do not exhibit the same overshoot. In fact, for a piecewise "C"1 function the Fourier series will converge to the function at every point except at jump discontinuities. At jump the discontinuities the limit will converge to the average of the values of the function on either side of the the jump. This is a consequences of the Dirichlet-Jordan theorem cite book
author=M. Pinsky
title=Introduction to Fourier Analysis and Wavelets
page=p 27
year= 2002
publisher=Brooks/Cole
location=United states of America
] . The Gibbs phenomenon is also closely related to the principle that the decay of the Fourier coefficients of a function at infinity is controlled by the smoothness of that function; very smooth functions will have very rapidly decaying Fourier coefficients (resulting in the rapid convergence of the Fourier series), whereas discontinuous functions will have very slowly decaying Fourier coefficients (causing the Fourier series to converge very slowly). Note for instance that the Fourier coefficients 1, 1/3, 1/5, dots of the discontinuous square wave described above decay only as fast as the harmonic series, which is not absolutely convergent; indeed, the above Fourier series turns out to be only conditionally convergent for almost every value of x. This provides a partial explanation of the Gibbs phenomenon, since Fourier series with absolutely convergent Fourier coefficients would be uniformly convergent by the Weierstrass M-test and would thus be unable to exhibit the above oscillatory behavior. By the same token, itis impossible for a discontinuous function to have absolutely convergent Fourier coefficients, since the function would thus be the uniform limit of continuous functions and therefore be continuous, a contradiction. See more about absolute convergence of Fourier series.

In practice, the difficulties associated with the Gibbs phenomenon can be ameliorated by using a smoother method of Fourier series summation, such as Fejér summation or Riesz summation, or by using sigma-approximation. Using a wavelet transform with Haar basis functions, the Gibbs phenomenon does not occur.

Formal mathematical description of the phenomenon

Let f: {Bbb R} o {Bbb R} be a piecewise continuously differentiable function which is periodic with some period L > 0. Suppose that at some point x_0, the left limit f(x_0^-) and right limit f(x_0^+) of the function f differ by a non-zero gap a:

: f(x_0^+) - f(x_0^-) = a eq 0.

For each positive integer "N" ≥ 1, let "S""N" "f" be the "N"th partial Fourier series

: S_N f(x) := sum_{-N leq n leq N} hat f(n) e^{2pi i n x/L}= frac{1}{2} a_0 + sum_{n=1}^N left( a_n cosleft(frac{2pi nx}{L} ight) + b_n sinleft(frac{2pi nx}{L} ight) ight)

where the Fourier coefficients hat f(n), a_n, b_n are given by the usual formulae

: hat f(n) := frac{1}{L} int_0^L f(x) e^{-2pi i n x/L}, dx

: a_n := frac{2}{L} int_0^L f(x) cosleft(frac{2pi nx}{L} ight), dx

: b_n := frac{2}{L} int_0^L f(x) sinleft(frac{2pi nx}{L} ight), dx.

Then we have

: lim_{N o infty} S_N fleft(x_0 + frac{L}{2N} ight) = f(x_0^+) + acdot (0.089490dots)

and

: lim_{N o infty} S_N fleft(x_0 - frac{L}{2N} ight) = f(x_0^-) - acdot (0.089490dots)

but

: lim_{N o infty} S_N f(x_0) = frac{f(x_0^-) + f(x_0^+)}{2}.More generally, if x_N is any sequence of real numbers which converges to x_0 as N o infty, and if the gap "a" is positive then: limsup_{N o infty} S_N f(x_N) leq f(x_0^+) + acdot (0.089490dots)and: liminf_{N o infty} S_N f(x_N) geq f(x_0^-) - acdot (0.089490dots)If instead the gap "a" is negative, one needs to interchange limit superior with limit inferior, and also interchangethe ≤ and ≥ signs, in the above two inequalities.

The square wave example

We now illustrate the above Gibbs phenomenon in the case of the square wave described earlier. In this case the period "L" is 2pi, the discontinuity x_0 is at zero, and the jump "a" is equal to pi/2.For simplicity let us just deal with the case when "N" is even (the case of odd "N" is very similar). Thenwe have

:S_N f(x) = sin(x) + frac{1}{3} sin(3x) + cdots + frac{1}{N-1} sin((N-1)x).

Substituting x=0, we obtain

:S_N f(0) = 0 = frac{-frac{pi}{4} + frac{pi}{4{2} = frac{f(0^-) + f(0^+)}{2}

as claimed above. Next, we compute

:S_N f(frac{2pi}{2N}) = sinleft(frac{pi}{N} ight) + frac{1}{3} sinleft(frac{3pi}{N} ight)+ cdots + frac{1}{N-1} sinleft( frac{(N-1)pi}{N} ight).

If we introduce the normalized sinc function, operatorname{sinc}(x),, we can rewrite this as

:S_N fleft(frac{2pi}{2N} ight) = frac{pi}{2} left [ frac{2}{N} operatorname{sinc}left(frac{1}{N} ight) + frac{2}{N} operatorname{sinc}left(frac{3}{N} ight)+ cdots + frac{2}{N} operatorname{sinc}left( frac{(N-1)}{N} ight) ight] .

But the expression in square brackets is a numerical integration approximation to the integral int_0^1 operatorname{sinc}(x) dx (more precisely, it is a midpoint rule approximation with spacing 2/N). Since the sinc function is continuous, this approximation converges to the actual integral as N o infty. Thus we have

:

which was what was claimed in the previous section. A similar computation shows

:lim_{N o infty} S_N fleft(-frac{2pi}{2N} ight) = -frac{pi}{2} int_0^1 operatorname{sinc}(x) dx = -frac{pi}{4} -frac{pi}{2} cdot (0.089490dots).

Consequences

The Gibbs phenomenon causes a compression artifact called pre-echo in audio compression algorithms that use Fourier-related transforms – particularly the modified discrete cosine transform (MDCT) – such as MP3, MPEG-4 AAC, and Vorbis.

Notes and references

ee also

* Compare with Runge's phenomenon for polynomial approximations
* Sigma approximation
* Henry Wilbraham

Publications

* Gibbs, J. W., "Fourier Series". Nature 59, 200 (1898) and 606 (1899).
* Antoni Zygmund, "Trigonometrical series", Dover publications, 1955.
* Wilbraham, H. "On a certain periodic function", Cambridge and Dublin Math. J., 3 (1848), pp. 198-201.
* Paul J. Nahin, "Dr. Euler's Fabulous Formula," Princeton University Press, 2006. Ch. 4, Sect. 4.

External links and references

* Weisstein, Eric W., " [http://mathworld.wolfram.com/GibbsPhenomenon.html Gibbs Phenomenon] ". From MathWorld--A Wolfram Web Resource.
* Prandoni, Paolo, " [http://lcavwww.epfl.ch/~prandoni/dsp/gibbs/gibbs.html Gibbs Phenomenon] ".
* Radaelli-Sanchez, Ricardo, and Richard Baraniuk, " [http://cnx.rice.edu/content/m10092/latest/ Gibbs Phenomenon] ". The Connexions Project. (Creative Commons Attribution License)


Wikimedia Foundation. 2010.

Игры ⚽ Нужно сделать НИР?

Look at other dictionaries:

  • Gibbs — may refer to:People*Cecil Armstrong Gibbs, composer *Cory Gibbs, soccer player *Frederic A. Gibbs, neurologist *George Gibbs (mineralogist), (1776 1833) *George Gibbs (geologist), (1815 1873) *Herschelle Gibbs, South African cricketer *Humphrey… …   Wikipedia

  • Gibbs-Effekt — Gibbssches Phänomen bei einer Rechteckskurve Als Gibbssches Phänomen oder „Ringing“ bezeichnet man in der Mathematik das typische Verhalten von Fourierreihen in der Umgebung von Sprungsstellen. Entwickelt man eine Fourierreihe aus einer… …   Deutsch Wikipedia

  • Gibbs-Phänomen — Gibbssches Phänomen bei einer Rechteckskurve Als Gibbssches Phänomen oder „Ringing“ bezeichnet man in der Mathematik das typische Verhalten von Fourierreihen in der Umgebung von Sprungsstellen. Entwickelt man eine Fourierreihe aus einer… …   Deutsch Wikipedia

  • Josiah Willard Gibbs — Infobox Scientist box width = 300px name = J. Willard Gibbs image size = 300px caption = Josiah Willard Gibbs birth date = birth date|1839|2|11|mf=y birth place = New Haven, Connecticut, USA death date = death date and… …   Wikipedia

  • Pinsky phenomenon — The Pinsky phenomenon is a result in Fourier analysis, a branch of mathematics [Cite journal volume = 27 issue = 3 pages = 565 605 last = Taylor first = M. E. title = THE GIBBS PHENOMENON, THE PINSKY PHENOMENON, AND VARIANTS FOR EIGENFUNCTION… …   Wikipedia

  • Josiah Willard Gibbs — J. Willard Gibbs Josiah Willard Gibbs Nacimiento 11 de febrero de 1839 New Haven, Connecticut …   Wikipedia Español

  • Runge's phenomenon — red|the Runge function. The blue curve is blue|a 5th order interpolating polynomial (using six equally spaced interpolating points). The green curve is green|a 9th order interpolating polynomial (using ten equally spaced interpolating points). At …   Wikipedia

  • Fredia Gibbs — Fredia “The Cheetah” Gibbs (born July 8, 1964) a.k.a. “The Most Dangerous Woman in the World” is a sports icon and Muay Thai kickboxing legend. Gibbs is a 3 Time World Champion in Kickboxing and the first African American female to hold the world …   Wikipedia

  • Scientific phenomena named after people — This is a list of scientific phenomena and concepts named after people (eponymous phenomena). For other lists of eponyms, see eponym. NOTOC A* Abderhalden ninhydrin reaction Emil Abderhalden * Abney effect, Abney s law of additivity William de… …   Wikipedia

  • Henry Wilbraham — Infobox Person name=Henry Wilbraham caption=Mathematician birth date=birth date|1825|7|25|mf=y birth place= dead=dead death date=death date and age|1883|2|13|1825|7|25|mf=y death place=Henry Wilbraham (July 25, 1825 – February 13, 1883) was an… …   Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”