In mathematics, the question of whether the Fourier series of a periodic function converges to a given function is researched by a field known as classical harmonic analysis, a branch of pure mathematics. Convergence is not necessarily given in the general case, and certain criteria must be met for convergence to occur.

Determination of convergence requires the comprehension of pointwise convergence, uniform convergence, absolute convergence, L^p spaces, summability methods and the Cesàro mean.

Preliminaries

Consider f an integrable function on the interval [0, 2π]. For such an f the Fourier coefficients ${\displaystyle {\widehat {f}}(n)}$ are defined by the formula

${\displaystyle {\widehat {f}}(n)={\frac {1}{2\pi }}\int _{0}^{2\pi }f(t)e^{-int}\,\mathrm {d} t,\quad n\in \mathbb {Z} .}$

It is common to describe the connection between f and its Fourier series by

${\displaystyle f\sim \sum _{n}{\widehat {f}}(n)e^{int}.}$

The notation ~ here means that the sum represents the function in some sense. To investigate this more carefully, the partial sums must be defined:

${\displaystyle S_{N}(f;t)=\sum _{n=-N}^{N}{\widehat {f}}(n)e^{int}.}$

The question of whether a Fourier series converges is: Do the functions ${\displaystyle S_{N}(f)}$ (which are functions of the variable t we omitted in the notation) converge to f and in which sense? Are there conditions on f ensuring this or that type of convergence?

Before continuing, the Dirichlet kernel must be introduced. Taking the formula for ${\displaystyle {\widehat {f}}(n)}$, inserting it into the formula for ${\displaystyle S_{N}}$ and doing some algebra gives that

${\displaystyle S_{N}(f)=f*D_{N}}$

where ∗ stands for the periodic convolution and ${\displaystyle D_{N}}$ is the Dirichlet kernel, which has an explicit formula,

${\displaystyle D_{n}(t)={\frac {\sin((n+{\frac {1}{2}})t)}{\sin(t/2)}}.}$

The Dirichlet kernel is not a positive kernel, and in fact, its norm diverges, namely

${\displaystyle \int |D_{n}(t)|\,\mathrm {d} t\to \infty }$

a fact that plays a crucial role in the discussion. The norm of D_n in L¹(T) coincides with the norm of the convolution operator with D_n, acting on the space C(T) of periodic continuous functions, or with the norm of the linear functional f → (S_nf)(0) on C(T). Hence, this family of linear functionals on C(T) is unbounded, when n → ∞.

Magnitude of Fourier coefficients

In applications, it is often useful to know the size of the Fourier coefficient.

If ${\displaystyle f}$ is an absolutely continuous function,

${\displaystyle \left|{\widehat {f}}(n)\right|\leq {K \over |n|}}$

for ${\displaystyle K}$ a constant that only depends on ${\displaystyle f}$.

If ${\displaystyle f}$ is a bounded variation function,

${\displaystyle \left|{\widehat {f}}(n)\right|\leq {{\rm {var}}(f) \over 2\pi |n|}.}$

If ${\displaystyle f\in C^{p}}$

${\displaystyle \left|{\widehat {f}}(n)\right|\leq {\|f^{(p)}\|_{L_{1}} \over |n|^{p}}.}$

If ${\displaystyle f\in C^{p}}$ and ${\displaystyle f^{(p)}}$ has modulus of continuity ${\displaystyle \omega _{p}}$,

${\displaystyle \left|{\widehat {f}}(n)\right|\leq {\omega (2\pi /n) \over |n|^{p}}}$

and therefore, if ${\displaystyle f}$ is in the α-Hölder class

${\displaystyle \left|{\widehat {f}}(n)\right|\leq {K \over |n|^{\alpha }}.}$

Pointwise convergence

Superposition of sinusoidal wave basis functions (bottom) to form a sawtooth wave (top); the basis functions have wavelengths λ/k (k=integer) shorter than the wavelength λ of the sawtooth itself (except for k=1). All basis functions have nodes at the nodes of the sawtooth, but all but the fundamental have additional nodes. The oscillation about the sawtooth is called the Gibbs phenomenon

There are many known sufficient conditions for the Fourier series of a function to converge at a given point x, for example if the function is differentiable at x. Even a jump discontinuity does not pose a problem: if the function has left and right derivatives at x, then the Fourier series converges to the average of the left and right limits (but see Gibbs phenomenon).

The Dirichlet–Dini Criterion states that: if ƒ is 2π–periodic, locally integrable and satisfies

${\displaystyle \int _{0}^{\pi }\left|{\frac {f(x_{0}+t)+f(x_{0}-t)}{2}}-\ell \right|{\frac {\mathrm {d} t}{t}}<\infty ,}$

then (S_nf)(x₀) converges to ℓ. This implies that for any function f of any Hölder class α > 0, the Fourier series converges everywhere to f(x).

It is also known that for any periodic function of bounded variation, the Fourier series converges everywhere. See also Dini test. In general, the most common criteria for pointwise convergence of a periodic function f are as follows:

• If f satisfies a Holder condition, then its Fourier series converges uniformly.
• If f is of bounded variation, then its Fourier series converges everywhere.
• If f is continuous and its Fourier coefficients are absolutely summable, then the Fourier series converges uniformly.

There exist continuous functions whose Fourier series converges pointwise but not uniformly; see Antoni Zygmund, Trigonometric Series, vol. 1, Chapter 8, Theorem 1.13, p. 300.

However, the Fourier series of a continuous function need not converge pointwise. Perhaps the easiest proof uses the non-boundedness of Dirichlet's kernel in L¹(T) and the Banach–Steinhaus uniform boundedness principle. As typical for existence arguments invoking the Baire category theorem, this proof is nonconstructive. It shows that the family of continuous functions whose Fourier series converges at a given x is of first Baire category, in the Banach space of continuous functions on the circle.

So in some sense pointwise convergence is atypical, and for most continuous functions the Fourier series does not converge at a given point. However Carleson's theorem shows that for a given continuous function the Fourier series converges almost everywhere.

It is also possible to give explicit examples of a continuous function whose Fourier series diverges at 0: for instance, the even and 2π-periodic function f defined for all x in [0,π] by^[1]

${\displaystyle f(x)=\sum _{n=1}^{\infty }{\frac {1}{n^{2}}}\sin \left[\left(2^{n^{3}}+1\right){\frac {x}{2}}\right].}$

Uniform convergence

Suppose ${\displaystyle f\in C^{p}}$, and ${\displaystyle f^{(p)}}$ has modulus of continuity ${\displaystyle \omega }$; then the partial sums of the Fourier series converge to the function with speed^[2]

${\displaystyle |f(x)-(S_{N}f)(x)|\leq K{\ln N \over N^{p}}\omega (2\pi /N)}$

for a constant ${\displaystyle K}$ that does not depend upon ${\displaystyle f}$, nor ${\displaystyle p}$, nor ${\displaystyle N}$.

This theorem, first proved by D Jackson, tells, for example, that if ${\displaystyle f}$ satisfies the ${\displaystyle \alpha }$-Hölder condition, then

${\displaystyle |f(x)-(S_{N}f)(x)|\leq K{\ln N \over N^{\alpha }}.}$

If ${\displaystyle f}$ is ${\displaystyle 2\pi }$ periodic and absolutely continuous on ${\displaystyle [0,2\pi ]}$, then the Fourier series of ${\displaystyle f}$ converges uniformly, but not necessarily absolutely, to ${\displaystyle f}$.^[3]

Absolute convergence

A function ƒ has an absolutely converging Fourier series if

${\displaystyle \|f\|_{A}:=\sum _{n=-\infty }^{\infty }|{\widehat {f}}(n)|<\infty .}$

Obviously, if this condition holds then ${\displaystyle (S_{N}f)(t)}$ converges absolutely for every t and on the other hand, it is enough that ${\displaystyle (S_{N}f)(t)}$ converges absolutely for even one t, then this condition holds. In other words, for absolute convergence there is no issue of where the sum converges absolutely — if it converges absolutely at one point then it does so everywhere.

The family of all functions with absolutely converging Fourier series is a Banach algebra (the operation of multiplication in the algebra is a simple multiplication of functions). It is called the Wiener algebra, after Norbert Wiener, who proved that if ƒ has absolutely converging Fourier series and is never zero, then 1/ƒ has absolutely converging Fourier series. The original proof of Wiener's theorem was difficult; a simplification using the theory of Banach algebras was given by Israel Gelfand. Finally, a short elementary proof was given by Donald J. Newman in 1975.

If ${\displaystyle f}$ belongs to a α-Hölder class for α > 1/2 then

${\displaystyle \|f\|_{A}\leq c_{\alpha }\|f\|_{{\rm {Lip}}_{\alpha }},\qquad \|f\|_{K}:=\sum _{n=-\infty }^{+\infty }|n||{\widehat {f}}(n)|^{2}\leq c_{\alpha }\|f\|_{{\rm {Lip}}_{\alpha }}^{2}}$

for ${\displaystyle \|f\|_{{\rm {Lip}}_{\alpha }}}$ the constant in the Hölder condition, ${\displaystyle c_{\alpha }}$ a constant only dependent on ${\displaystyle \alpha }$; ${\displaystyle \|f\|_{K}}$ is the norm of the Krein algebra. Notice that the 1/2 here is essential—there are 1/2-Hölder functions, which do not belong to the Wiener algebra. Besides, this theorem cannot improve the best known bound on the size of the Fourier coefficient of a α-Hölder function—that is only ${\displaystyle O(1/n^{\alpha })}$ and then not summable.

If ƒ is of bounded variation and belongs to a α-Hölder class for some α > 0, it belongs to the Wiener algebra.

Norm convergence

The simplest case is that of L², which is a direct transcription of general Hilbert space results. According to the Riesz–Fischer theorem, if ƒ is square-integrable then

${\displaystyle \lim _{N\rightarrow \infty }\int _{0}^{2\pi }\left|f(x)-S_{N}(f)(x)\right|^{2}\,\mathrm {d} x=0}$

i.e.,  ${\displaystyle S_{N}f}$ converges to ƒ in the norm of L². It is easy to see that the converse is also true: if the limit above is zero, ƒ must be in L². So this is an if and only if condition.

If 2 in the exponents above is replaced with some p, the question becomes much harder. It turns out that the convergence still holds if 1 < p < ∞. In other words, for ƒ in L^p,  ${\displaystyle S_{N}(f)}$ converges to ƒ in the L^p norm. The original proof uses properties of holomorphic functions and Hardy spaces, and another proof, due to Salomon Bochner relies upon the Riesz–Thorin interpolation theorem. For p = 1 and infinity, the result is not true. The construction of an example of divergence in L¹ was first done by Andrey Kolmogorov (see below). For infinity, the result is a corollary of the uniform boundedness principle.

If the partial summation operator S_N is replaced by a suitable summability kernel (for example the Fejér sum obtained by convolution with the Fejér kernel), basic functional analytic techniques can be applied to show that norm convergence holds for 1 ≤ p < ∞.

Convergence almost everywhere

The problem whether the Fourier series of any continuous function converges almost everywhere was posed by Nikolai Lusin in the 1920s. It was resolved positively in 1966 by Lennart Carleson. His result, now known as Carleson's theorem, tells the Fourier expansion of any function in L² converges almost everywhere. Later on, Richard Hunt generalized this to L^p for any p > 1.

Contrariwise, Andrey Kolmogorov, as a student at the age of 19, in his very first scientific work, constructed an example of a function in L¹ whose Fourier series diverges almost everywhere (later improved to diverge everywhere).

Jean-Pierre Kahane and Yitzhak Katznelson proved that for any given set E of measure zero, there exists a continuous function ƒ such that the Fourier series of ƒ fails to converge on any point of E.

Summability

Does the sequence 0,1,0,1,0,1,... (the partial sums of Grandi's series) converge to ½? This does not seem like a very unreasonable generalization of the notion of convergence. Hence we say that any sequence ${\displaystyle (a_{n})_{n=1}^{\infty }}$ is Cesàro summable to some a if

${\displaystyle \lim _{n\to \infty }{\frac {1}{n}}\sum _{k=1}^{n}s_{k}=a.}$

Where with ${\displaystyle s_{k}}$ we denote the kth partial sum:

${\displaystyle s_{k}=a_{1}+\cdots +a_{k}=\sum _{n=1}^{k}a_{n}}$

It is not difficult to see that if a sequence converges to some a then it is also Cesàro summable to it.

To discuss summability of Fourier series, we must replace ${\displaystyle S_{N}}$ with an appropriate notion. Hence we define

${\displaystyle K_{N}(f;t)={\frac {1}{N}}\sum _{n=0}^{N-1}S_{n}(f;t),\quad N\geq 1,}$

and ask: does ${\displaystyle K_{N}(f)}$ converge to f? ${\displaystyle K_{N}}$ is no longer associated with Dirichlet's kernel, but with Fejér's kernel, namely

${\displaystyle K_{N}(f)=f*F_{N}\,}$

where ${\displaystyle F_{N}}$ is Fejér's kernel,

${\displaystyle F_{N}={\frac {1}{N}}\sum _{n=0}^{N-1}D_{n}.}$

The main difference is that Fejér's kernel is a positive kernel. Fejér's theorem states that the above sequence of partial sums converge uniformly to ƒ. This implies much better convergence properties

• If ƒ is continuous at t then the Fourier series of ƒ is summable at t to ƒ(t). If ƒ is continuous, its Fourier series is uniformly summable (i.e. ${\displaystyle K_{N}(f)}$ converges uniformly to ƒ).
• For any integrable ƒ, ${\displaystyle K_{N}(f)}$ converges to ƒ in the ${\displaystyle L^{1}}$ norm.
• There is no Gibbs phenomenon.

Results about summability can also imply results about regular convergence. For example, we learn that if ƒ is continuous at t, then the Fourier series of ƒ cannot converge to a value different from ƒ(t). It may either converge to ƒ(t) or diverge. This is because, if ${\displaystyle S_{N}(f;t)}$ converges to some value x, it is also summable to it, so from the first summability property above, x = ƒ(t).

Order of growth

The order of growth of Dirichlet's kernel is logarithmic, i.e.

${\displaystyle \int |D_{N}(t)|\,\mathrm {d} t={\frac {4}{\pi ^{2}}}\log N+O(1).}$

See Big O notation for the notation O(1). The actual value ${\displaystyle 4/\pi ^{2}}$ is both difficult to calculate (see Zygmund 8.3) and of almost no use. The fact that for some constant c we have

${\displaystyle \int |D_{N}(t)|\,\mathrm {d} t>c\log N+O(1)}$

is quite clear when one examines the graph of Dirichlet's kernel. The integral over the n-th peak is bigger than c/n and therefore the estimate for the harmonic sum gives the logarithmic estimate.

This estimate entails quantitative versions of some of the previous results. For any continuous function f and any t one has

${\displaystyle \lim _{N\to \infty }{\frac {S_{N}(f;t)}{\log N}}=0.}$

However, for any order of growth ω(n) smaller than log, this no longer holds and it is possible to find a continuous function f such that for some t,

${\displaystyle \varlimsup _{N\to \infty }{\frac {S_{N}(f;t)}{\omega (N)}}=\infty .}$

The equivalent problem for divergence everywhere is open. Sergei Konyagin managed to construct an integrable function such that for every t one has

${\displaystyle \varlimsup _{N\to \infty }{\frac {S_{N}(f;t)}{\sqrt {\log N}}}=\infty .}$

It is not known whether this example is best possible. The only bound from the other direction known is log n.

Multiple dimensions

Upon examining the equivalent problem in more than one dimension, it is necessary to specify the precise order of summation one uses. For example, in two dimensions, one may define

${\displaystyle S_{N}(f;t_{1},t_{2})=\sum _{|n_{1}|\leq N,|n_{2}|\leq N}{\widehat {f}}(n_{1},n_{2})e^{i(n_{1}t_{1}+n_{2}t_{2})}}$

which are known as "square partial sums". Replacing the sum above with

${\displaystyle \sum _{n_{1}^{2}+n_{2}^{2}\leq N^{2}}}$

lead to "circular partial sums". The difference between these two definitions is quite notable. For example, the norm of the corresponding Dirichlet kernel for square partial sums is of the order of ${\displaystyle \log ^{2}N}$ while for circular partial sums it is of the order of ${\displaystyle {\sqrt {N}}}$.

Many of the results true for one dimension are wrong or unknown in multiple dimensions. In particular, the equivalent of Carleson's theorem is still open for circular partial sums. Almost everywhere convergence of "square partial sums" (as well as more general polygonal partial sums) in multiple dimensions was established around 1970 by Charles Fefferman.

Notes

References