Hi all, This Blog is an English archive of my PhD experience in Imperial College London, mainly logging my research and working process, as well as some visual records.

Wednesday 1 August 2007

Fourier series

The Fourier series is a mathematical tool used for analyzing periodic functions by decomposing such a function into a weighted sum of much simpler sinusoidal component functions sometimes referred to as normal Fourier modes, or simply modes for short. The weights, or coefficients, of the modes, are a one-to-one mapping of the original function. Generalizations

include generalized Fourier series and other expansions over orthonormal bases.

Definition

General form

Given a complex-valued function f of real argument t, f: RC, where f(t) is piecewise smooth and continuous, periodic with period T, and square-integrable over the interval from t1 to t2 of length T, that is,

 \int_{t_1}^{t_2} |f(t)|^2\, dt<+\infty

where

  • T = t2t1 is the period,
  • t1 and t2 are integration bounds.

The Fourier series expansion of f is

  •  f(t) = \frac{1}{2} a_0 + \sum_{n=1}^{\infty}[a_n \cos(\omega_n t) + b_n \sin(\omega_n t)]

where, for any non-negative integer n,

  •  \omega_n = n\frac{2\pi}{T} is the nth harmonic (in radians) of the function f,
  • a_n = \frac{2}{T}\int_{t_1}^{t_2} f(t) \cos(\omega_n t)\, dt are the even Fourier coefficients of f, and
  • b_n = \frac{2}{T}\int_{t_1}^{t_2} f(t) \sin(\omega_n t)\, dt are the odd Fourier coefficients of f.

Equivalently, in complex exponential form,

  • f(t) = \sum_{n=-\infty}^{+\infty} c_n e^{i \omega_n t}

where:

For a formal justification, see Modern derivation of the Fourier coefficients below.

Canonical form

In the special case where the period T = 2π, we have

\omega_n = n \,

In this case, the Fourier series expansion reduces to a particularly simple form:

f(t) = \frac{1}{2} a_0 +\sum_{n=1}^{\infty}[a_n \cos(nt) + b_n \sin(nt)]

where

  • a_n = \frac{1}{\pi}\int_{-\pi}^{\pi} f(t) \cos(nt)\, dt
  • b_n = \frac{1}{\pi}\int_{-\pi}^{\pi} f(t) \sin(nt)\, dt

for any non-negative integer n.

or, equivalently:

f(t) = \sum_{n=-\infty}^{+\infty} c_n e^{i nt}

where

  • c_n = \frac{1}{2 \pi}\int_{-\pi}^{\pi} f(t) e^{-i nt}\, dt = \frac{1}{2}(a_n-ib_n).

Choice of the form

The form for period T can be easily derived from the canonical one with the change of variable defined by x=\frac{2\pi}{T}t. Therefore, both formulations are equivalent. However, the form for period T is used in most practical cases because it is directly applicable. For the theory, the canonical form is preferred because it is more elegant and easier to interpret mathematically, as will later be seen.

Examples

Simple Fourier series

Let f be periodic of period , with f(x) = x for x from −π to π. Note that this function is a periodic version of the identity function.

Plot of a periodic identity function - a sawtooth wave.

Plot of a periodic identity function - a sawtooth wave.
Animated plot of the first five successive partial Fourier series

Animated plot of the first five successive partial Fourier series

We will compute the Fourier coefficients for this function.

\begin{align} a_n &{}= \frac{1}{\pi}\int_{-\pi}^{\pi}f(x)\cos(nx)\,dx \\ &{}= \frac{1}{\pi}\int_{-\pi}^{\pi}x \cos(nx)\,dx \\ &{}= 0. \end{align}


\begin{align} b_n &{}= \frac{1}{\pi}\int_{-\pi}^{\pi}f(x)\sin(nx)\,dx \\ &{}= \frac{1}{\pi}\int_{-\pi}^{\pi} x \sin(nx)\, dx \\ &{}= \frac{2}{\pi}\int_{0}^{\pi} x\sin(nx)\, dx \\ &{}= \frac{2}{\pi} \left(\left[-\frac{x\cos(nx)}{n}\right]_0^{\pi} + \left[\frac{\sin(nx)}{n^2}\right]_0^{\pi}\right) \\ &{}= 2\frac{(-1)^{n+1}}{n}.\end{align}

Notice that an are 0 because the x\mapsto x\cos(nx) are even functions. Hence the Fourier series for this function is:

f(x)=\frac{a_0}{2} + \sum_{n=1}^{\infty}\left[a_n\cos\left(nx\right)+b_n\sin\left(nx\right)\right]
=2\sum_{n=1}^{\infty}\frac{(-1)^{n+1}}{n} \sin(nx), \quad \forall x\in [-\pi,\pi].

One application of this Fourier series is to compute the value of the Riemann zeta function at s = 2; by Parseval's theorem, we have:

0}\left[2\frac{(-1)^n}{n}\right]^2" src="http://upload.wikimedia.org/math/e/7/7/e777681f5b32083c015c6ba54578e77f.png">

which yields: 0}\frac{1}{n^2}=\frac{\pi^2}{6}" src="http://upload.wikimedia.org/math/3/5/9/359981037d4391ebe85b915baae0b473.png">.

No comments: