Fourier Series

A Fourier series is an expansion of a periodic function

in terms of an infinite sum of sines and cosines. Fourier series make use of theorthogonality relationships of the sine and cosine functions. The computation and study of Fourier series is known as harmonic analysis and is extremely useful as a way to break up an arbitrary periodic function into a set of simple terms that can be plugged in, solved individually, and then recombined to obtain the solution to the original problem or an approximation to it to whatever accuracy is desired or practical. Examples of successive approximations to common functions using Fourier series are illustrated above.

In particular, since the superposition principle holds for solutions of a linear homogeneous ordinary differential equation, if such an equation can be solved in the case of a single sinusoid, the solution for an arbitrary function is immediately available by expressing the original function as a Fourier series and then plugging in the solution for each sinusoidal component. In some special cases where the Fourier series can be summed in closed form, this technique can even yield analytic solutions.

Any set of functions that form a complete orthogonal system have a corresponding generalized Fourier series analogous to the Fourier series. For example, using orthogonality of the roots of a Bessel function of the first kind gives a so-called Fourier-Bessel series.

The computation of the (usual) Fourier series is based on the integral identities

			(1)
			(2)
			(3)
			(4)
			(5)

for

, where

is the Kronecker delta.

Using the method for a generalized Fourier series, the usual Fourier series involving sines and cosines is obtained by taking

and

. Since these functions form a complete orthogonal system over

, the Fourier series of a function

is given by

(6)

where

			(7)
			(8)
			(9)

and

, 2, 3, .... Note that the coefficient of the constant term

has been written in a special form compared to the general form for a generalized Fourier series in order to preserve symmetry with the definitions of

and

The Fourier cosine coefficient

and sine coefficient

are implemented in the Wolfram Language as FourierCosCoefficient[expr, t, n] andFourierSinCoefficient[expr, t, n], respectively.

A Fourier series converges to the function

(equal to the original function at points of continuity or to the average of the two limits at points of discontinuity)

$f^_={1/2[lim_(x->x_0^-)f(x)+lim_(x->x_0^+)f(x)] for -pi<x_0<pi; 1/2[lim_(x->-pi^+)f(x)+lim_(x->pi_-)f(x)] for x_0=-pi,pi$

(10)

if the function satisfies so-called Dirichlet boundary conditions. Dini's test gives a condition for the convergence of Fourier series.

As a result, near points of discontinuity, a "ringing" known as the Gibbs phenomenon, illustrated above, can occur.

For a function

periodic on an interval

instead of

, a simple change of variables can be used to transform the interval of integration from

. Let

			(11)
			(12)

Solving for

gives

, and plugging this in gives

(13)

Therefore,

			(14)
			(15)
			(16)

Similarly, the function is instead defined on the interval

, the above equations simply become

			(17)
			(18)
			(19)

In fact, for

periodic with period

, any interval

can be used, with the choice being one of convenience or personal preference (Arfken 1985, p. 769).

The coefficients for Fourier series expansions of a few common functions are given in Beyer (1987, pp. 411-412) and Byerly (1959, p. 51). One of the most common functions usually analyzed by this technique is the square wave. The Fourier series for a few common functions are summarized in the table below.

function		Fourier series
Fourier series--sawtooth wave
Fourier series--square wave
Fourier series--triangle wave

If a function is even so that

, then

is odd. (This follows since

is odd and an even function times an odd function is an odd function.) Therefore,

for all

. Similarly, if a function is odd so that

, then

is odd. (This follows since

is even and an even function times an odd function is an odd function.) Therefore,

for all

The notion of a Fourier series can also be extended to complex coefficients. Consider a real-valued function

. Write

(20)

Now examine

			(21)
			(22)
		$sum_(n=-infty)^(infty)A_nint_(-pi)^pi{cos[(n-m)x]+isin[(n-m)x]}dx$	(23)
			(24)
			(25)

(26)

The coefficients can be expressed in terms of those in the Fourier series

			(27)
		${1/(2pi)int_(-pi)^pif(x)[cos(nx)+isin(\|n\|x)]dx n<0; 1/(2pi)int_(-pi)^pif(x)dx n=0; 1/(2pi)int_(-pi)^pif(x)[cos(nx)-isin(nx)]dx n>0$	(28)
		${1/2(a_n+ib_n) for n<0; 1/2a_0 for n=0; 1/2(a_n-ib_n) for n>0.$	(29)