Gibbs phenomenon Guide, Meaning , Facts, Information and Description
The Gibbs phenomenon (also known as ringing artifacts) is the peculiar manner in which the Fourier series of a piecewise continuously differentiable periodic function f behaves at a jump discontinuity: the nth 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.
The three pictures on the right demonstrate this for a square wave whose Fourier expansion is
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 of the square wave by
The Gibbs phenomenon was first observed by Albert Michelson via a mechanical graphing machine. Michelson developed a device in 1898 that could compute and re-synthesize the Fourier series. When a square wave was input into the machine, the graph would move to and from around the discontinuities. This would occur, and continue to occur, as the number of Fourier coefficients approached infinity.
The phenomenon was first explained mathematically by J. Willard Gibbs in 1899. Informally, it reflects the difficulty inherent in approximating a discontinuous function by a series of continuous sine and cosine waves. This 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 (and thus very rapidly convergent Fourier series), whereas discontinuous functions will have very slowly decaying Fourier coefficients (and thus very badly convergent Fourier series). Note for instance that the Fourier coefficients 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, it is 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. If one uses a wavelet transform instead of the Fourier transform then the Gibbs phenomenon no longer occurs.
| Table of contents |
|
2 The square wave example 3 See also 4 Publications 5 External links and references |
Let be a piecewise continuously differentiable function which is periodic with some period . Suppose that at some point , the left limit and right limit of the function differ by a non-zero gap :
Formal mathematical description of the phenomenon
For each positive integer , let be the th partial Fourier series
where the Fourier coefficients are given by the usual formulae
Then we have
and
but
More generally, if is any sequence of real numbers which converges to as , and if the gap a is positive then
and
If instead the gap a is negative, one needs to interchange limit superior with limit inferior, and also interchange
the ≤ and ≥ signs, in the above two inequalities.
We now illustrate the above Gibbs phenomenon in the case of the square wave described earlier. In this case the period L is , the discontinuity is at zero, and the jump a is equal to .
For simplicity let us just deal with the case when N is even (the case of odd N is very similar). Then
we have
The square wave example
Substituting , we obtain
as claimed above. Next, we compute
If we introduce the sinc function , we can rewrite this as
See also
- List of mathematical topics
- Square wave
- Scientific phenomena named after people
- Yale University
- Compare with Runge's phenomenon for polynomial approximations
Publications
- Gibbs, J. W., "Fourier Series". Nature 59, 200 and 606, 1899.
- Antoni Zygmund, Trigonometrical series, Dover publications, 1955.
External links and references
- Braennlund, Johan, "Why are sine waves fundamental".
- Weisstein, Eric W., "Gibbs Phenomenon". From MathWorld--A Wolfram Web Resource.
- Prandoni, Paolo, "Gibbs Phenomenon".
- Radaelli-Sanchez, Ricardo, and Richard Baraniuk, "Gibbs Phenomenon". The Connexions Project. (Creative Commons Attribution License)
- Pavel, "Gibbs phenomenon". math.mit.edu. (Java applet)
This is an Article on Gibbs phenomenon. Page Contains Information, Facts Details or Explanation Guide About Gibbs phenomenon
