Fourier Series — L2 Theory and Convergence

⚡ Key Takeaways

The Fourier series of an L2L^2 function converges to it in the L2L^2 norm — always, unconditionally — because {einx/2π}n∈Z\{e^{inx}/\sqrt{2\pi}\}_{n\in\mathbb{Z}} is a complete orthonormal basis of L2([−π,π])L^2([-\pi,\pi]). This is the Riesz–Fischer theorem in disguise. Pointwise/uniform convergence is a separate, harder, and conditional story (Carleson's theorem). Foundational to harmonic analysis, signal processing, and spectral methods; the L2L^2 theory crystallized 1900–1910 (Fischer, Riesz, Parseval), pointwise a.e. convergence resolved by Carleson in 1966.

💡 The Core Insight

L2([−π,π])L^2([-\pi,\pi]) is a separable Hilbert space, and the trigonometric system {en(x)=12πeinx}\{e_n(x) = \tfrac{1}{\sqrt{2\pi}}e^{inx}\} is an orthonormal basis of it. Once you accept this, the entire convergence theory is the abstract Hilbert-space theory of orthonormal bases: partial sums are orthogonal projections onto finite-dimensional subspaces, Bessel's inequality becomes Parseval's equality, and "convergence" means norm convergence — which holds for free. No regularity of ff is needed. All the genuine difficulty in Fourier analysis lives outside the L2L^2 norm.

Formal Statement

DefinitionFourier Coefficients and Partial Sums

For f∈L1([−π,π])f \in L^1([-\pi,\pi]) (in particular f∈L2f \in L^2, since the interval has finite measure), the nn-th Fourier coefficient is

1
f^(n)=12π∫−ππf(x) e−inx dx,n∈Z.\quad \hat f(n) = \frac{1}{2\pi}\int_{-\pi}^{\pi} f(x)\,e^{-inx}\,dx, \qquad n \in \mathbb{Z}.

Derivation hint: This is the inner product ⟨f,en⟩\langle f, e_n\rangle with en=12πeinxe_n = \tfrac{1}{\sqrt{2\pi}}e^{inx}, divided by 2π\sqrt{2\pi}; it is the coefficient that minimizes mean-square error onto ene_n. The NN-th symmetric partial sum is

2
SNf(x)=∑n=−NNf^(n) einx.\quad S_N f(x) = \sum_{n=-N}^{N} \hat f(n)\,e^{inx}.

Derivation hint: Orthogonal projection of ff onto span{einx:∣n∣≤N}\mathrm{span}\{e^{inx}: |n|\le N\} expanded in the basis {einx}\{e^{inx}\}. Here L2([−π,π])L^2([-\pi,\pi]) carries the inner product and norm

3
⟨f,g⟩=12π∫−ππf(x)g(x)‾ dx,∥f∥2=⟨f,f⟩1/2.\quad \langle f,g\rangle = \frac{1}{2\pi}\int_{-\pi}^{\pi} f(x)\overline{g(x)}\,dx, \qquad \|f\|_2 = \langle f,f\rangle^{1/2}.

Derivation hint: The normalization 12π\tfrac{1}{2\pi} makes {einx}\{e^{inx}\} orthonormal: ⟨eimx,einx⟩=δmn\langle e^{imx}, e^{inx}\rangle = \delta_{mn}, by direct integration.

TheoremRiesz–Fischer / L2L^2 Convergence

The system {einx}n∈Z\{e^{inx}\}_{n\in\mathbb{Z}} is a complete orthonormal basis of L2([−π,π])L^2([-\pi,\pi]) (with inner product [3]). Consequently, for every f∈L2([−π,π])f \in L^2([-\pi,\pi]):

4
lim⁡N→∞∥f−SNf∥2=0,\quad \lim_{N\to\infty} \|f - S_N f\|_2 = 0,

Derivation hint: Completeness ⇔\Leftrightarrow the only function orthogonal to all einxe^{inx} is 00 a.e.; combine with Bessel to force the partial-sum tail to zero. and Parseval's identity holds:

5
∥f∥22=∑n=−∞∞∣f^(n)∣2=12π∫−ππ∣f(x)∣2 dx.\quad \|f\|_2^2 = \sum_{n=-\infty}^{\infty} |\hat f(n)|^2 = \frac{1}{2\pi}\int_{-\pi}^{\pi}|f(x)|^2\,dx.

Derivation hint: Take N→∞N\to\infty in ∥f∥22=∥SNf∥22+∥f−SNf∥22\|f\|_2^2 = \|S_Nf\|_2^2 + \|f-S_Nf\|_2^2 (Pythagoras), using [4]. Moreover, the map f↦(f^(n))n∈Zf \mapsto (\hat f(n))_{n\in\mathbb{Z}} is a unitary isomorphism L2([−π,π])≅ℓ2(Z)L^2([-\pi,\pi]) \cong \ell^2(\mathbb{Z}). In particular, for any sequence (cn)∈ℓ2(Z)(c_n)\in\ell^2(\mathbb{Z}) there exists f∈L2f\in L^2 with f^(n)=cn\hat f(n) = c_n.

Why It Holds / Derivation

ProofCompleteness ⟹ L² Convergence and Parseval

Step 1 — Best approximation. Fix NN and let VN=span{einx:∣n∣≤N}V_N = \mathrm{span}\{e^{inx}: |n|\le N\}. Among all g∈VNg\in V_N, the partial sum SNfS_N f minimizes ∥f−g∥2\|f-g\|_2. Indeed for g=∑∣n∣≤Ncneinxg = \sum_{|n|\le N} c_n e^{inx},

∥f−g∥22=∥f∥22−2 Re∑∣n∣≤Ncn‾f^(n)+∑∣n∣≤N∣cn∣2,\|f-g\|_2^2 = \|f\|_2^2 - 2\,\mathrm{Re}\sum_{|n|\le N}\overline{c_n}\hat f(n) + \sum_{|n|\le N}|c_n|^2,

which (completing the square per coordinate) is minimized exactly at cn=f^(n)c_n = \hat f(n). Hence SNfS_Nf is the orthogonal projection of ff onto VNV_N and f−SNf⊥VNf - S_N f \perp V_N.

Step 2 — Pythagoras and Bessel. Since f−SNf⊥SNff - S_Nf \perp S_Nf,

∥f∥22=∥SNf∥22+∥f−SNf∥22=∑∣n∣≤N∣f^(n)∣2+∥f−SNf∥22.\|f\|_2^2 = \|S_N f\|_2^2 + \|f-S_Nf\|_2^2 = \sum_{|n|\le N}|\hat f(n)|^2 + \|f-S_Nf\|_2^2.

Dropping the nonnegative remainder gives Bessel's inequality ∑∣n∣≤N∣f^(n)∣2≤∥f∥22\sum_{|n|\le N}|\hat f(n)|^2 \le \|f\|_2^2, so ∑n∣f^(n)∣2<∞\sum_n |\hat f(n)|^2 < \infty.

Step 3 — Completeness. We must show ∥f−SNf∥2→0\|f-S_Nf\|_2\to 0. Equivalently, the only obstruction is whether the einxe^{inx} span a dense subspace. By the Stone–Weierstrass theorem, trigonometric polynomials are uniformly dense in the continuous 2π2\pi-periodic functions C(T)C(\mathbb{T}); uniform density on a finite-measure interval implies L2L^2-density. And C(T)C(\mathbb{T}) is dense in L2(T)L^2(\mathbb{T}). So trigonometric polynomials are dense in L2L^2.

Step 4 — Density ⟹ convergence. Given ε>0\varepsilon>0, choose a trig polynomial PP of degree MM with ∥f−P∥2<ε\|f-P\|_2<\varepsilon. For N≥MN\ge M, P∈VNP\in V_N, and since SNfS_N f is the best VNV_N-approximation (Step 1),

∥f−SNf∥2≤∥f−P∥2<ε.\|f - S_N f\|_2 \le \|f-P\|_2 < \varepsilon.

Thus ∥f−SNf∥2→0\|f-S_Nf\|_2\to 0, establishing [4].

Step 5 — Parseval. Let N→∞N\to\infty in the Pythagoras identity of Step 2; the remainder vanishes by [4], yielding [5].

Step 6 — Unitarity / surjectivity (Riesz–Fischer). Given (cn)∈ℓ2(c_n)\in\ell^2, the partial sums gN=∑∣n∣≤Ncneinxg_N = \sum_{|n|\le N} c_n e^{inx} form a Cauchy sequence in L2L^2 because ∥gN−gM∥22=∑M<∣n∣≤N∣cn∣2→0\|g_N - g_M\|_2^2 = \sum_{M<|n|\le N}|c_n|^2 \to 0. By completeness of L2L^2 (Riesz–Fischer's metric core), gN→gg_N\to g in L2L^2, and g^(n)=cn\hat g(n) = c_n by continuity of the inner product. ■\qquad\blacksquare

∎
⚠️ Two distinct "completenesses"

The proof uses two completeness facts with the same name. (a) L2L^2 is a complete metric space (Cauchy sequences converge) — this gives surjectivity in Step 6. (b) The trigonometric system is complete as an orthonormal family (no nonzero vector orthogonal to all of it) — this gives [4]. Do not conflate them; both are essential and they are logically independent statements.

mathematician · Frigyes Riesz | 1880–1956

Riesz, working in Hungary, and Ernst Fischer in Austria independently proved in 1907 the theorem now bearing both names: L2L^2 is complete and isomorphic to ℓ2\ell^2 via Fourier coefficients. This was a foundational moment — it justified Hilbert's emerging abstract theory of square-summable sequences by anchoring it in concrete function spaces, and it made "convergence in mean" the natural mode for Fourier series. Riesz went on to lay much of the groundwork for functional analysis (the Riesz representation theorem, LpL^p duality, compact operators).

Building Intuition

🧠 Intuition

Think of L2L^2 as an infinite-dimensional Euclidean space and {einx}\{e^{inx}\} as its coordinate axes. The Fourier coefficient f^(n)\hat f(n) is literally the component of the "vector" ff along the nn-th axis. Parseval [5] is the Pythagorean theorem: the squared length of ff equals the sum of squares of its components. L2L^2 convergence [4] just says that if you keep more axes, the projection SNfS_N f gets closer to ff — exactly as in Rk\mathbb{R}^k, projecting onto more coordinate axes recovers the vector. The reason this is "easy" while pointwise convergence is "hard": the L2L^2 norm only sees average discrepancy. A reconstructed signal can overshoot wildly at a single point (Gibbs phenomenon) yet have tiny mean-square error, because the overshoot occupies vanishing measure. The norm forgives pointwise sins; pointwise convergence does not.

Figure
Orthogonal projection: SNfS_N f is the closest point in VNV_N to ff; the error is perpendicular
V_N = span{e^{inx}, |n|≤N} f S_N f f − S_N f ⟂ V_N 0 ‖f‖² = ‖S_N f‖² + ‖f − S_N f‖² (Pythagoras → Parseval as N→∞)

Worked Examples

ExampleSawtooth: coefficients, Parseval, and a series sum

Given. f(x)=xf(x) = x on [−π,π][-\pi,\pi], extended 2π2\pi-periodically. Compute f^(n)\hat f(n), verify L2L^2 membership, and use Parseval to evaluate ∑n≥11/n2\sum_{n\ge1} 1/n^2.

Step 1 — Coefficients. f^(0)=12π∫−ππx dx=0\hat f(0) = \frac{1}{2\pi}\int_{-\pi}^\pi x\,dx = 0 (odd integrand). For n≠0n\ne 0, integrate by parts:

f^(n)=12π∫−ππxe−inx dx=12π[xe−inx−in]−ππ+12πin∫−ππe−inxdx.\hat f(n) = \frac{1}{2\pi}\int_{-\pi}^{\pi} x e^{-inx}\,dx = \frac{1}{2\pi}\left[\frac{x e^{-inx}}{-in}\right]_{-\pi}^{\pi} + \frac{1}{2\pi in}\int_{-\pi}^\pi e^{-inx}dx.

The integral term vanishes; the boundary term gives f^(n)=12π⋅2π(−1)n−in⋅(sign bookkeeping)=(−1)n in\hat f(n) = \frac{1}{2\pi}\cdot\frac{2\pi(-1)^{n}}{-in}\cdot(\text{sign bookkeeping}) = \frac{(-1)^{n}\, i}{n}. So ∣f^(n)∣=1/∣n∣|\hat f(n)| = 1/|n|.

Step 2 — L2L^2 norm. ∥f∥22=12π∫−ππx2 dx=12π⋅2π33=π23.\|f\|_2^2 = \frac{1}{2\pi}\int_{-\pi}^\pi x^2\,dx = \frac{1}{2\pi}\cdot\frac{2\pi^3}{3} = \frac{\pi^2}{3}. Finite, so f∈L2f\in L^2. ✓

Step 3 — Parseval [5].

π23=∑n≠01n2=2∑n=1∞1n2.\frac{\pi^2}{3} = \sum_{n\ne 0}\frac{1}{n^2} = 2\sum_{n=1}^\infty \frac{1}{n^2}.

Answer.

  ∑n=1∞1n2=π26  \boxed{\;\sum_{n=1}^{\infty}\frac{1}{n^2} = \frac{\pi^2}{6}\;}

The Fourier series converges to xx in L2L^2, even though pointwise it converges to 00 at x=±πx=\pm\pi (the average of the jump) — a vivid reminder that L2L^2 convergence ignores point values.

ExampleA pure ℓ2\ell^2 coefficient sequence with no classical function feel

Given. Define cn=1nc_n = \frac{1}{n} for n≥1n\ge 1 and cn=0c_n = 0 otherwise. Is there an L2L^2 function with these as Fourier coefficients? Estimate ∥f∥2\|f\|_2.

Step 1 — ℓ2\ell^2 check. ∑n∣cn∣2=∑n=1∞1/n2=π2/6<∞\sum_n |c_n|^2 = \sum_{n=1}^\infty 1/n^2 = \pi^2/6 < \infty, so (cn)∈ℓ2(Z)(c_n)\in\ell^2(\mathbb{Z}). ✓

Step 2 — Riesz–Fischer surjectivity. By Theorem [4]/[5] (Step 6 of the proof), there exists f∈L2([−π,π])f\in L^2([-\pi,\pi]) with f^(n)=cn\hat f(n)=c_n, namely the L2L^2-limit of gN=∑n=1N1neinxg_N = \sum_{n=1}^N \frac{1}{n}e^{inx}.

Step 3 — Norm via Parseval.

∥f∥22=∑n∣cn∣2=π26≈1.6449.\|f\|_2^2 = \sum_{n}|c_n|^2 = \frac{\pi^2}{6} \approx 1.6449.

Answer.

  ∥f∥2=π6≈1.2825  \boxed{\;\|f\|_2 = \frac{\pi}{\sqrt 6}\approx 1.2825\;}

Note: this ff is not given by an elementary closed form, yet it is a perfectly legitimate L2L^2 function — the existence is guaranteed abstractly by completeness, with no formula required. (In fact ∑1neinx\sum \frac1n e^{inx} relates to −log⁡(1−eix)-\log(1-e^{ix}), but the point stands: ℓ2\ell^2 coefficients ⟺ L2L^2 function, unconditionally.)

Common Misconceptions

⚠️ L2L^2 convergence does NOT imply pointwise convergence

∥f−SNf∥2→0\|f-S_Nf\|_2\to0 guarantees convergence in mean square only. It implies a subsequence SNkf→fS_{N_k}f\to f pointwise a.e., but not the full sequence. There exist continuous functions whose Fourier series diverge at a point (du Bois-Reymond, 1873) and integrable functions whose series diverge everywhere (Kolmogorov, 1926) — yet all such ff in L2L^2 still satisfy [4]. Never write "SNf(x)→f(x)S_Nf(x)\to f(x)" as a consequence of the Riesz–Fischer theorem.

CounterexampleGibbs phenomenon — L2L^2 fine, pointwise overshoot persists

Take the square wave f=sgn(x)f = \mathrm{sgn}(x) on [−π,π][-\pi,\pi]. Its partial sums SNfS_N f converge to ff in L2L^2 (since f∈L2f\in L^2, automatic by [4]). But near the jump at x=0x=0, the partial sums overshoot the limit by a fixed factor ≈8.95%\approx 8.95\% of the jump height, for every NN — the overshoot does not shrink, it merely migrates toward the discontinuity. The maximum of SNfS_Nf converges to 2π∫0πsin⁡tt dt≈1.0895\frac{2}{\pi}\int_0^\pi \frac{\sin t}{t}\,dt \approx 1.0895 times the half-jump. The L2L^2 norm is undisturbed because the overshoot region has width O(1/N)→0O(1/N)\to 0: it contributes O(1/N)→0O(1/N)\to 0 to ∥f−SNf∥22\|f-S_Nf\|_2^2. Mean-square convergence and uniform convergence are genuinely different.

⚠️ Bessel vs. Parseval

Bessel's inequality ∑∣f^(n)∣2≤∥f∥22\sum|\hat f(n)|^2 \le \|f\|_2^2 holds for any orthonormal system. Equality (Parseval, [5]) holds iff the system is complete. For an incomplete system (e.g. drop ei0xe^{i0x} from the basis), Bessel is strict for functions with a nonzero mean. The leap from inequality to equality is the content of completeness — it is not automatic.

Connections

connection

The L2L^2 theory sits at a crossroads. It is the prototypical instance of Orthonormal Bases in Hilbert Spaces: every separable Hilbert space looks like ℓ2\ell^2, and Fourier series is the canonical realization. Via the Plancherel Theorem it generalizes from the circle to R\mathbb{R}, where the Fourier transform is the unitary map L2(R)→L2(R)L^2(\mathbb{R})\to L^2(\mathbb{R}). Abstractly it is a special case of the Spectral Theorem — the einxe^{inx} are simultaneous eigenfunctions of translation (and of −d2/dx2-d^2/dx^2), and Fourier series diagonalizes these operators. The convergence dichotomy connects to Lp Convergence of Fourier Series (M. Riesz: SNf→fS_Nf\to f in LpL^p for 1<p<∞1<p<\infty, false at p=1,∞p=1,\infty) and to Carleson's Theorem on a.e. convergence. The coefficient decay rate ties to Sobolev Spaces: f∈Hsf\in H^s iff ∑(1+n2)s∣f^(n)∣2<∞\sum (1+n^2)^s|\hat f(n)|^2<\infty.

Real-World Applications

🔬 Where this appears

Historical Development

⏳ History
- 1807 / 1822 — Joseph Fourier asserts (in his heat-conduction work, published fully in Théorie analytique de la chaleur) that arbitrary functions admit trigonometric expansions, provoking controversy among Lagrange, Laplace, and Poisson over convergence.
- 1829 — Dirichlet gives the first rigorous pointwise convergence theorem for piecewise-monotone functions with finitely many discontinuities.
- 1873 — Paul du Bois-Reymond constructs a continuous function whose Fourier series diverges at a point, shattering naive optimism about pointwise convergence.
- 1899 — J. Willard Gibbs (and earlier Wilbraham, 1848) explains the persistent overshoot near jumps — the Gibbs phenomenon.
- 1900s — Hilbert and Schmidt develop the abstract theory of ℓ2\ell^2 and integral equations, setting the stage for Hilbert space.
- 1907 — Frigyes Riesz and Ernst Fischer independently prove the Riesz–Fischer theorem: L2L^2 is complete and isometric to ℓ2\ell^2 via Fourier coefficients, establishing unconditional mean-square convergence.
- 1926 — Andrey Kolmogorov constructs an L1L^1 function whose Fourier series diverges everywhere — emphatically separating L1L^1 from L2L^2 behavior.
- 1966 — Lennart Carleson proves that the Fourier series of any L2L^2 function converges almost everywhere, resolving a conjecture open since Luzin (1915).
- 1968 — Richard Hunt extends Carleson's theorem to LpL^p for all p>1p>1; Kahane–Katznelson (1966) show a.e. convergence can fail on any prescribed measure-zero set, pinning down the sharp result.

Subtleties & Limits

open-problem

The L2L^2 convergence theory is complete and closed — there are no open problems in [4]/[5] itself; they are theorems. The frontier lies in adjacent refinements. (a) Sharp constants and rates: the exact dependence of ∥f−SNf∥2\|f-S_Nf\|_2 on smoothness classes interacts with multiplier theory and remains subtle for endpoint and weighted spaces. (b) Convergence in higher dimensions: for Td\mathbb{T}^d, d≥2d\ge2, spherical partial sums ∑∣n∣≤R\sum_{|n|\le R} exhibit the notorious Bochner–Riesz problem — for which LpL^p and which summation order does convergence hold? This is open in full generality for d≥3d\ge3 and connects to the Kakeya and restriction conjectures, among the deepest unsolved problems in harmonic analysis. (c) Carleson-type results in higher dimensions and on general groups remain only partially understood. So while the Hilbert-space core is settled, replacing the L2L^2 norm by LpL^p or by pointwise/multidimensional notions reopens some of the hardest questions in modern analysis.

See Also

Hilbert Spaces · Orthonormal Bases in Hilbert Spaces · Plancherel Theorem · Carleson's Theorem · Lp Convergence of Fourier Series · Spectral Theorem · Sobolev Spaces · Gibbs Phenomenon · Riesz–Fischer Theorem · Bochner–Riesz Summability