Fourier Analysis

From the discrete Fourier transform to singular integrals. The page starts from two things you already know, characters of finite abelian groups and uniform convergence from undergraduate analysis, then builds Fourier series, the transform on ℝⁿ, distributions, and the real-variable methods of harmonic analysis: maximal functions, the Hilbert transform, Littlewood–Paley theory and BMO. Every module has numbered examples and a long list of exercises, with the important ones starred.

Follows L. Grafakos, Classical Fourier Analysis (3rd ed., Springer GTM 249, 2014), Chapters 1–6, and Modern Fourier Analysis (3rd ed., GTM 250, 2014), Chapters 1 and 3; section and result numbers refer to those books. The first two modules are shared and appear identically on other pages: Abelian groups and the Fourier transform from Representation Theory, and Uniform convergence from Real Analysis, which supplies the undergraduate analysis the classical theory is built on. ★ marks the exercises worth doing first. The examples and exercises are written for this page. For how this subject connects to the others, see the roadmap.

InteractiveEtingof et al. §3 (the abelian case); Grafakos, Classical Fourier Analysis, §3.1

Abelian groups and the Fourier transform

For an abelian group every irreducible representation is 1-dimensional, and breaking the space of functions on the group into those pieces is exactly the Fourier transform. The discrete Fourier transform, Fourier series and the Fourier integral are the same theorem for three different groups.

The idea

Characters of an abelian group. A character of a finite abelian group GG is a homomorphism χ:G→ℂ×\chi:G\to\mathbb C^\times; by Exercise 2.2 of the representation theory page these are all its irreducible representations. They form a group G^\hat G under pointwise multiplication, the dual group, and G^≅G\hat G\cong G (Exercise 4.4). For G=ℤ/nG=\mathbb Z/n, χk(x)=e2πikx/n\chi_k(x)=e^{2\pi ikx/n}, k=0,…,n−1k=0,\dots,n-1.
Fourier transform and inversion. For f:G→ℂf:G\to\mathbb C put f^(χ)=1|G|∑xf(x)χ(x)\hat f(\chi)=\frac1{|G|}\sum_xf(x)\overline{\chi(x)}. Since the characters are an orthonormal basis of functions on GG for ⟨f,g⟩=1|G|∑fg‾\langle f,g\rangle=\frac1{|G|}\sum f\bar g (Exercise 4.1),
f=∑χ∈G^f^(χ)χ,1|G|∑x|f(x)|2=∑χ|f^(χ)|2(Plancherel).f=\sum_{\chi\in\hat G}\hat f(\chi)\,\chi,\qquad \frac1{|G|}\sum_x|f(x)|^2=\sum_\chi|\hat f(\chi)|^2\quad\text{(Plancherel)}.
In the language of representations: functions on GG are the regular representation, and f^(χ)χ\hat f(\chi)\chi is the component of ff in the copy of χ\chi.
Convolution becomes multiplication. With (f*g)(x)=∑yf(y)g(x−y)(f*g)(x)=\sum_yf(y)g(x-y), one has f*g^=|G|f^g^\widehat{f*g}=|G|\,\hat f\,\hat g (Exercise 4.2). Translation-invariant linear operators (convolutions, circulant matrices) are intertwiners of the regular representation, so by Schur's lemma each acts on each character by a scalar: the Fourier transform diagonalizes all of them at once.
Three groups, one theorem. For the circle 𝕋=ℝ/2πℤ\mathbb T=\mathbb R/2\pi\mathbb Z the continuous characters are eikθe^{ik\theta}, k∈ℤk\in\mathbb Z, and the transform is the sequence of Fourier coefficients (Grafakos §3.1). For ℝ\mathbb R they are e2πixξe^{2\pi ix\xi}, ξ∈ℝ\xi\in\mathbb R, and the transform is the Fourier integral. Sums become integrals against Haar measure; the algebra is the same. For non-abelian compact groups the irreducibles are no longer 1-dimensional, and the Peter–Weyl theorem replaces characters by matrix entries of irreducible representations.

Explore: a signal and its Fourier transform on ℤ/n\mathbb Z/n

Signal f
ff (drag the bars) , and f*bf*b for a box bb of the chosen width
|f^(χk)||\hat f(\chi_k)| and |f*b^(χk)|=n|f^||b^||\widehat{f*b}(\chi_k)|=n|\hat f||\hat b| , for k=0,…,n−1k=0,\dots,n-1

Key terms in plain words

The ideas above without the symbols. Switch to Quiz me to test yourself: you see an explanation and pick the term it describes.

Dual group
The collection of all "pure frequencies" of an abelian group (homomorphisms to the unit circle), which is itself a group under pointwise multiplication.
Fourier transform
Writing a function on a group as a combination of pure frequencies, and recording how much of each frequency it contains.
Plancherel identity
The total energy of a signal equals the total energy of its frequency components (with the right normalisation): nothing is lost or gained by changing to frequencies.
Convolution
Sliding one function along another and adding up the products: a weighted moving average. In frequencies it becomes plain multiplication.
Circulant matrix
A square matrix in which each row is the previous row shifted one place to the right, wrapping around. All of them share the same eigenvectors, the pure frequencies.
Peter–Weyl theorem
The version of Fourier analysis for compact groups that need not be abelian: functions are expanded using the matrix entries of irreducible representations instead of single frequencies.

Examples

Example 4.1 (ℤ/2\mathbb Z/2).

The characters are χ0=(1,1)\chi_0=(1,1) and χ1=(1,−1)\chi_1=(1,-1), so f^(χ0)=f(0)+f(1)2\hat f(\chi_0)=\frac{f(0)+f(1)}2 and f^(χ1)=f(0)−f(1)2\hat f(\chi_1)=\frac{f(0)-f(1)}2: the average and half the difference: the splitting of a representation of ℤ/2\mathbb Z/2 into the parts where the generator acts by +1+1 and by −1-1.

Example 4.2 (the DFT matrix).

For ℤ/n\mathbb Z/n the matrix Fkx=1nω−kxF_{kx}=\frac1n\omega^{-kx}, ω=e2πi/n\omega=e^{2\pi i/n}, sends ff to f^\hat f. Orthogonality of characters says nF\sqrt nF is unitary. The fast Fourier transform computes it in O(nlog⁡n)O(n\log n) steps instead of O(n2)O(n^2), by splitting ℤ/2m\mathbb Z/2m into even and odd residues.

Example 4.3 (a spike and a constant).

The spike δ0\delta_0 has δ^0(χ)=1|G|\hat\delta_0(\chi)=\frac1{|G|} for every χ\chi: it contains every frequency equally. The constant 1 has 1^=δχ0\hat1=\delta_{\chi_0}. Concentrated on one side means spread out on the other (Exercise 4.7). Try "spike" in the figure.

Example 4.4 (moving averages).

Convolving with a box b=1w𝟏{0,…,w−1}b=\frac1w\mathbf 1_{\{0,\dots,w-1\}} replaces ff by its average over a window. Its transform nb^(χk)=1w∑j<wω−jkn\hat b(\chi_k)=\frac1w\sum_{j<w}\omega^{-jk} is close to 1 for small kk and small for most large kk, and vanishes exactly at the k≠0k\ne0 with n|kwn\mid kw: averaging keeps low frequencies and damps high ones.

Example 4.5 (Boolean functions).

For G=(ℤ/2)nG=(\mathbb Z/2)^n, written multiplicatively as {±1}n\{\pm1\}^n, the characters are χS(x)=∏i∈Sxi\chi_S(x)=\prod_{i\in S}x_i for S⊆{1,…,n}S\subseteq\{1,\dots,n\}. The transform is the Walsh–Hadamard transform, the basic tool of the analysis of Boolean functions. For example Maj(x1,x2,x3)=12(x1+x2+x3)−12x1x2x3\mathrm{Maj}(x_1,x_2,x_3)=\frac12(x_1+x_2+x_3)-\frac12x_1x_2x_3.

Example 4.6 (Fourier series).

On the circle, f^(k)=12π∫02πf(θ)e−ikθdθ\hat f(k)=\frac1{2\pi}\int_0^{2\pi}f(\theta)e^{-ik\theta}\,d\theta, and Plancherel reads 12π∫|f|2=∑k|f^(k)|2\frac1{2\pi}\int|f|^2=\sum_k|\hat f(k)|^2. The finite theory has no convergence questions; on the circle, whether and in what sense ∑f^(k)eikθ\sum\hat f(k)e^{ik\theta} converges back to ff is the subject of Grafakos Chapters 3 and 4.

Check your understanding

Why does the discrete Fourier transform diagonalize every circulant matrix?

  1. Circulant matrices commute with translation, i.e. they intertwine the regular representation of ℤ/n\mathbb Z/n, and each character spans its own 1-dimensional piece, so each is an eigenvector.
  2. Because circulant matrices are symmetric.Most are not, e.g. the shift itself.
  3. Because the DFT matrix is unitary.Unitarity alone doesn't make it diagonalize anything in particular.
The pieces of the regular representation of an abelian group are distinct 1-dimensional representations, so an intertwiner must map each into itself, acting by a scalar (Exercise 4.5).

How many characters does ℤ/2×ℤ/4\mathbb Z/2\times\mathbb Z/4 have?

  1. 6That is 2+42+4; characters of a product are products of characters.
  2. 8
  3. 4An abelian group has as many characters as elements.
Every element is its own conjugacy class, so there are |G|=8|G|=8 irreducibles, all 1-dimensional: χ(a,b)=(−1)jaikb\chi(a,b)=(-1)^{ja}i^{kb}.

Exercises

Exercise 4.1Characters are an orthonormal basis

For a finite abelian group, show directly that distinct characters are orthonormal for ⟨f,g⟩=1|G|∑fg‾\langle f,g\rangle=\frac1{|G|}\sum f\bar g, that there are |G||G| of them, and hence that Fourier inversion holds.

Hint

If χ≠ψ\chi\ne\psi, then χψ‾\chi\bar\psi is a nontrivial character η\eta; pick aa with η(a)≠1\eta(a)\ne1 and compare ∑xη(x)\sum_x\eta(x) with ∑xη(ax)\sum_x\eta(ax). For the count, use ∑di2=|G|\sum d_i^2=|G| with all di=1d_i=1.

Solution

⟨χ,χ⟩=1\langle\chi,\chi\rangle=1 since |χ|=1|\chi|=1. If η=χψ‾\eta=\chi\bar\psi is nontrivial, S=∑xη(x)S=\sum_x\eta(x) satisfies S=∑xη(ax)=η(a)SS=\sum_x\eta(ax)=\eta(a)S with η(a)≠1\eta(a)\ne1, so S=0S=0. By Schur's lemma (Exercise 2.2 on the representation theory page) all irreducibles are characters, and ∑di2=|G|\sum d_i^2=|G| gives |G||G| of them. |G||G| orthonormal functions in a |G||G|-dimensional space are a basis, and the coefficient of χ\chi in ff is ⟨f,χ⟩=f^(χ)\langle f,\chi\rangle=\hat f(\chi).

Exercise 4.2The convolution theorem

With the conventions above, prove f*g^=|G|f^g^\widehat{f*g}=|G|\,\hat f\,\hat g.

Hint

Substitute x=y+zx=y+z and use χ(y+z)=χ(y)χ(z)\chi(y+z)=\chi(y)\chi(z).

Solution

f*g^(χ)=1|G|∑x,yf(y)g(x−y)χ(x)=1|G|∑yf(y)χ(y)∑zg(z)χ(z)=1|G|⋅|G|f^(χ)⋅|G|g^(χ)\widehat{f*g}(\chi)=\frac1{|G|}\sum_{x,y}f(y)g(x-y)\overline{\chi(x)}=\frac1{|G|}\sum_yf(y)\overline{\chi(y)}\sum_zg(z)\overline{\chi(z)}=\frac1{|G|}\cdot|G|\hat f(\chi)\cdot|G|\hat g(\chi).

Exercise 4.3Plancherel

Prove ∑x|f(x)|2=|G|∑χ|f^(χ)|2\sum_x|f(x)|^2=|G|\sum_\chi|\hat f(\chi)|^2, and the polarized form for ∑xf(x)g(x)\sum_xf(x)\overline{g(x)}.

Hint

Expand f=∑f^(χ)χf=\sum\hat f(\chi)\chi and use orthonormality.

Solution

∑x|f(x)|2=∑χ,ψf^(χ)f^(ψ)∑xχ(x)ψ(x)=|G|∑χ|f^(χ)|2\sum_x|f(x)|^2=\sum_{\chi,\psi}\hat f(\chi)\overline{\hat f(\psi)}\sum_x\chi(x)\overline{\psi(x)}=|G|\sum_\chi|\hat f(\chi)|^2, since the inner sum is |G|δχψ|G|\delta_{\chi\psi}. The same computation with gg in the second slot gives ∑fg‾=|G|∑f^g^\sum f\bar g=|G|\sum\hat f\,\overline{\hat g}.

Exercise 4.4The dual group

Show that G^\hat G is a group under pointwise multiplication, that the dual group of ℤ/n\mathbb Z/n is isomorphic to ℤ/n\mathbb Z/n, and that G×H^≅G^×H^\widehat{G\times H}\cong\hat G\times\hat H. Deduce G^≅G\hat G\cong G for every finite abelian group (using the structure theorem). Note that the isomorphism is not canonical, while G≅G^^G\cong\hat{\hat G} is.

Hint

A character of ℤ/n\mathbb Z/n is determined by χ(1)\chi(1), an nnth root of unity.

Exercise 4.5Eigenvalues of circulant matrices

Let Cxy=c(x−y)C_{xy}=c(x-y) on ℤ/n\mathbb Z/n. Show Cχk=nc^(χk)χkC\chi_k=n\hat c(\chi_k)\,\chi_k.

Hint

(Cχ)(x)=∑yc(x−y)χ(y)(C\chi)(x)=\sum_yc(x-y)\chi(y); substitute z=x−yz=x-y.

Solution

(Cχ)(x)=∑zc(z)χ(x−z)=χ(x)∑zc(z)χ(z)=nc^(χ)χ(x)(C\chi)(x)=\sum_zc(z)\chi(x-z)=\chi(x)\sum_zc(z)\overline{\chi(z)}=n\hat c(\chi)\,\chi(x), using χ(−z)=χ(z)\chi(-z)=\overline{\chi(z)}. So CC is c*⋅c*\cdot, and Exercise 4.2 is the same statement.

Exercise 4.6Majority on three bits

Verify the expansion of Maj\mathrm{Maj} in Example 4.5 by computing f^(S)=18∑xf(x)χS(x)\hat f(S)=\frac18\sum_xf(x)\chi_S(x) for every S⊆{1,2,3}S\subseteq\{1,2,3\}.

Hint

By symmetry the coefficients depend only on |S||S|; even-size sets give 0 because Maj(−x)=−Maj(x)\mathrm{Maj}(-x)=-\mathrm{Maj}(x).

Exercise 4.7An uncertainty principle

Show that for f≠0f\ne0 on a finite abelian group, |suppf|⋅|suppf^|≥|G||\mathrm{supp}\,f|\cdot|\mathrm{supp}\,\hat f|\ge|G|. When is it an equality?

Hint

‖f‖∞≤∑χ|f^(χ)|≤|suppf^|1/2(∑|f^|2)1/2\lVert f\rVert_\infty\le\sum_\chi|\hat f(\chi)|\le|\mathrm{supp}\,\hat f|^{1/2}\big(\sum|\hat f|^2\big)^{1/2}, then Plancherel and ∑|f|2≤|suppf|‖f‖∞2\sum|f|^2\le|\mathrm{supp}\,f|\,\lVert f\rVert_\infty^2. Equality holds for indicator functions of subgroups (and their translates and modulations).

Exercise 4.8Characters of the circle

Show that the continuous homomorphisms ℝ/2πℤ→ℂ×\mathbb R/2\pi\mathbb Z\to\mathbb C^\times are θ↦eikθ\theta\mapsto e^{ik\theta}, k∈ℤk\in\mathbb Z.

Hint

A continuous homomorphism ϕ:ℝ→ℂ×\phi:\mathbb R\to\mathbb C^\times is t↦ectt\mapsto e^{ct} (smooth it by integrating: ∫0εϕ(t+s)ds=ϕ(t)∫0εϕ(s)ds\int_0^\varepsilon\phi(t+s)ds=\phi(t)\int_0^\varepsilon\phi(s)ds). Being 2π2\pi-periodic forces c=ikc=ik.

InteractiveRudin, Ch. 7: Sequences and Series of Functions

Uniform convergence

Pointwise limits of continuous functions can be discontinuous. Uniform convergence is the fix, and it is just convergence in the metric d(f,g)=sup⁡|f−g|d(f,g)=\sup|f-g|: the story of continuity, applied to a space whose points are functions.

The idea

Pointwise vs uniform (Rudin 7.1, 7.7). fn→ff_n\to f pointwise if fn(x)→f(x)f_n(x)\to f(x) for each xx (the NN may depend on xx). Uniformly if for every ε\varepsilon one NN works for all xx: equivalently supx⁡|fn(x)−f(x)|→0\sup_x|f_n(x)-f(x)|\to0 (Rudin 7.9).
Uniform limits keep continuity (Rudin 7.12). If continuous fn→ff_n\to f uniformly, ff is continuous. Uniform limits also commute with integrals (Rudin 7.16), but not with derivatives without more hypotheses (Rudin 7.17).

In the language of the story: C([a,b])C([a,b]) with the sup metric is a complete metric space (Rudin 7.15), and the Weierstrass theorem (Rudin 7.26) says polynomials are dense in it. Pointwise convergence, by contrast, is convergence in a topology that comes from no metric.

Key terms in plain words

The ideas above without the symbols. Switch to Quiz me to test yourself: you see an explanation and pick the term it describes.

Pointwise convergence
At each single input the values settle down eventually, but some inputs may take far longer than others.
Uniform convergence
The whole graph settles down together: from some step on, every input is within ε\varepsilon of the limit at once.
Sup norm
The largest gap between two functions anywhere on the domain. Convergence in it is exactly the "whole graph settles down together" kind.
Uniform limit theorem
If continuous functions converge with the whole graph settling down together, the limit is continuous too. Limits taken one point at a time can break continuity, as xnx^n on [0,1][0,1] shows.
Weierstrass approximation theorem
Any continuous function on a closed interval can be matched as closely as you like, everywhere at once, by a polynomial.
Complete metric space
A space with no holes: every sequence whose terms bunch up together has a limit inside the space.

Explore

Sequence fn
The limit f in ink, with an ε-tube around it. fn is in blue where it stays in the tube and red where it leaves.

Worked example

Prove Rudin 7.12: a uniform limit of continuous functions is continuous.

  1. Fix pp and ε\varepsilon. Which nn do you use?

    One with sup⁡|fn−f|<ε/3\sup|f_n-f|<\varepsilon/3. This is where uniformity matters: the same nn is good at pp and at every nearby xx.

  2. Now use continuity of that fnf_n at pp.

    Pick δ\delta with |fn(x)−fn(p)|<ε/3|f_n(x)-f_n(p)|<\varepsilon/3 when |x−p|<δ|x-p|<\delta. Then |f(x)−f(p)|≤|f(x)−fn(x)|+|fn(x)−fn(p)|+|fn(p)−f(p)|<ε|f(x)-f(p)|\le|f(x)-f_n(x)|+|f_n(x)-f_n(p)|+|f_n(p)-f(p)|<\varepsilon. This "ε/3\varepsilon/3 argument" is the template for many proofs in analysis.

Examples

Example 7.1 (a discontinuous limit).

xnx^n on [0,1][0,1] converges pointwise to the function that is 0 on [0,1)[0,1) and 1 at 1. The limit is discontinuous, so the convergence cannot be uniform: supx<1⁡xn=1\sup_{x<1}x^n=1 for every nn.

Example 7.2 (uniform on bounded sets only).

x/n→0x/n\to0 uniformly on every bounded interval but not on ℝ\mathbb R, where sup⁡|x/n|=∞\sup|x/n|=\infty. Many convergence statements are uniform only on compact sets.

Example 7.3 (limits and integrals disagree).

fn(x)=2n2xe−n2x2f_n(x)=2n^2xe^{-n^2x^2} tends to 0 at every x∈[0,1]x\in[0,1], but ∫01fn=1−e−n2→1\int_0^1f_n=1-e^{-n^2}\to1. The bump escapes towards 0, which uniform convergence would forbid.

Example 7.4 (the M-test).

∑n2−ncos⁡(3nx)\sum_n2^{-n}\cos(3^nx) converges uniformly by comparison with ∑2−n\sum2^{-n}, so it is continuous; by Hardy's theorem it is differentiable nowhere. Uniform limits keep continuity but not differentiability.

Check your understanding

fn(x)=xnf_n(x)=x^n on [0,1][0,1] converges pointwise to ff with f(1)=1f(1)=1 and f=0f=0 elsewhere. Is the convergence uniform?

  1. Yes, since each fnf_n is continuous.Continuity of the fnf_n doesn't force uniform convergence.
  2. No: a uniform limit of continuous functions would be continuous.
Rudin 7.12 in reverse: ff is discontinuous at 1, so the convergence can't be uniform. Directly, sup[0,1)⁡xn=1\sup_{[0,1)}x^n=1 for every nn. On [0,0.9][0,0.9] it is uniform, with sup⁡=0.9n\sup=0.9^n.

fn(x)=nxe−nxf_n(x)=nxe^{-nx} on [0,∞)[0,\infty) tends to 0 pointwise. Then supx⁡fn(x)\sup_x f_n(x)

  1. tends to 0.Look at x=1/nx=1/n.
  2. equals 1/e1/e for every nn.
The maximum is at x=1/nx=1/n with value e−1e^{-1}. The bump slides toward 0 and never shrinks: pointwise but not uniform convergence, even though the limit is continuous.

Practice

Exercise 7.1The Weierstrass M-test (Rudin 7.10)

Show that ∑n≥1sin⁡(nx)n2\sum_{n\ge1}\dfrac{\sin(nx)}{n^2} converges uniformly on ℝ\R to a continuous function.

Solution

|sin⁡(nx)/n2|≤1/n2|\sin(nx)/n^2|\le1/n^2 and ∑1/n2<∞\sum1/n^2<\infty. The partial sums are uniformly Cauchy: supx⁡|∑k=mn…|≤∑k=mn1/k2→0\sup_x|\sum_{k=m}^n\dots|\le\sum_{k=m}^n1/k^2\to0. So they converge uniformly, and the limit is continuous by Rudin 7.12.

Exercise 7.2Limits and integrals: when can you swap them?

For fn(x)=n2xe−nxf_n(x)=n^2xe^{-nx} on [0,1][0,1], compute lim⁡∫01fn\lim\int_0^1f_n and ∫01lim⁡fn\int_0^1\lim f_n.

Solution

fn→0f_n\to0 pointwise, so the second is 0. But ∫01n2xe−nxdx=∫0nue−udu→1\int_0^1n^2xe^{-nx}dx=\int_0^nue^{-u}du\to1. The convergence is not uniform (sup⁡fn=n/e\sup f_n=n/e), so Rudin 7.16 doesn't apply, and the swap fails.

From Rudin, Principles of Mathematical Analysis (3rd ed., McGraw–Hill, 1976)

Numbers match the book. Each question is restated in our own words; the hints and solutions are written for this page.

Rudin 7.1Uniform limits of bounded functions

Prove that every uniformly convergent sequence of bounded functions is uniformly bounded.

Hint

Use one NN from the uniform Cauchy condition, then handle the finitely many earlier functions separately.

Solution

Let fn→ff_n\to f uniformly on EE. Choose NN with |fn(x)−fN(x)|<1|f_n(x)-f_N(x)|<1 for all n≥Nn\ge N and all xx (the uniform Cauchy condition, Rudin 7.8). fNf_N is bounded by some BB, so |fn|<B+1|f_n|<B+1 for n≥Nn\ge N. The finitely many f1,…,fN−1f_1,\dots,f_{N-1} are bounded by B1,…,BN−1B_1,\dots,B_{N-1}. Then max⁡(B+1,B1,…,BN−1)\max(B+1,B_1,\dots,B_{N-1}) bounds all fnf_n.

Rudin 7.2Sums and products

If fn→ff_n\to f and gn→gg_n\to g uniformly on EE, prove fn+gn→f+gf_n+g_n\to f+g uniformly. If in addition the fnf_n and gng_n are bounded, prove fngn→fgf_ng_n\to fg uniformly.

Hint

For the product, add and subtract fngf_ng and use Rudin 7.1.

Solution

Sums: sup⁡|(fn+gn)−(f+g)|≤sup⁡|fn−f|+sup⁡|gn−g|→0\sup|(f_n+g_n)-(f+g)|\le\sup|f_n-f|+\sup|g_n-g|\to0.

Products: by Rudin 7.1 there is BB with |fn|,|gn|≤B|f_n|,|g_n|\le B, and then |f|,|g|≤B|f|,|g|\le B too. Then |fngn−fg|≤|fn||gn−g|+|g||fn−f|≤Bsup⁡|gn−g|+Bsup⁡|fn−f|→0|f_ng_n-fg|\le|f_n||g_n-g|+|g||f_n-f|\le B\sup|g_n-g|+B\sup|f_n-f|\to0.

Without boundedness it fails: on ℝ\R, fn=gn=x+1/nf_n=g_n=x+1/n converge uniformly to xx, but fngn−x2=2x/n+1/n2f_ng_n-x^2=2x/n+1/n^2 is unbounded for each nn.

Rudin 7.9Uniform convergence and moving points

Let fnf_n be continuous and fn→ff_n\to f uniformly on EE. Prove fn(xn)→f(x)f_n(x_n)\to f(x) whenever xn→xx_n\to x in EE. Is the converse true?

Hint

Split |fn(xn)−f(x)||f_n(x_n)-f(x)| through f(xn)f(x_n), using that ff is continuous (Rudin 7.12).

Solution

|fn(xn)−f(x)|≤|fn(xn)−f(xn)|+|f(xn)−f(x)|≤sup⁡|fn−f|+|f(xn)−f(x)||f_n(x_n)-f(x)|\le|f_n(x_n)-f(x_n)|+|f(x_n)-f(x)|\le\sup|f_n-f|+|f(x_n)-f(x)|. The first term tends to 0 by uniformity, the second by continuity of ff.

The converse fails. Take fn(x)=x/nf_n(x)=x/n on E=ℝE=\R. If xn→xx_n\to x, then fn(xn)=xn/n→0=f(x)f_n(x_n)=x_n/n\to0=f(x), since xnx_n is bounded. But supℝ⁡|x/n|=∞\sup_{\R}|x/n|=\infty, so the convergence is not uniform. (On a compact EE the converse is true: this condition is called continuous convergence.)

InteractiveGrafakos, Classical, §§1.2.2–1.2.4

Convolution and approximate identities

Convolving with a narrow bump of total mass 1 replaces a function by local averages. As the bump narrows, the averages return the function: uniformly if it is uniformly continuous, in LpL^p if it is in LpL^p. Almost every approximation theorem in Fourier analysis, Fejér's included, is a case of this one.

The idea

Convolution (Grafakos §1.2.2). (f*g)(x)=∫f(y)g(x−y)dy(f*g)(x)=\int f(y)g(x-y)\,dy, on ℝn\mathbb R^n or on the torus 𝕋n=[0,1]n\mathbb T^n=[0,1]^n with periodic functions. It is commutative and associative, commutes with translations, and ∂(f*g)=f*∂g\partial(f*g)=f*\partial g when gg is smooth: convolution makes things as smooth as the smoother factor.
Basic inequalities (Grafakos §1.2.3). Minkowski's integral inequality gives ‖f*g‖p≤‖f‖1‖g‖p\lVert f*g\rVert_p\le\lVert f\rVert_1\lVert g\rVert_p for 1≤p≤∞1\le p\le\infty; Hölder gives ‖f*g‖∞≤‖f‖p‖g‖p′\lVert f*g\rVert_\infty\le\lVert f\rVert_p\lVert g\rVert_{p'}. Young's inequality interpolates: ‖f*g‖r≤‖f‖p‖g‖q\lVert f*g\rVert_r\le\lVert f\rVert_p\lVert g\rVert_q when 1+1r=1p+1q1+\frac1r=\frac1p+\frac1q.
Approximate identities (Grafakos §1.2.4). A family KεK_\varepsilon with ∫Kε=1\int K_\varepsilon=1, supε⁡‖Kε‖1<∞\sup_\varepsilon\lVert K_\varepsilon\rVert_1<\infty, and ∫|x|≥δ|Kε|→0\int_{|x|\ge\delta}|K_\varepsilon|\to0 as ε→0\varepsilon\to0 for each δ>0\delta>0. Then Kε*f→fK_\varepsilon*f\to f uniformly if ff is bounded and uniformly continuous, in LpL^p if f∈Lpf\in L^p with 1≤p<∞1\le p<\infty, and at each point of continuity of a bounded ff. The proof is the ε/3\varepsilon/3 argument from the uniform convergence module: split the integral into |y|<δ|y|<\delta and |y|≥δ|y|\ge\delta.
The standard kernels. The box 12ε𝟏[−ε,ε]\frac1{2\varepsilon}\mathbf 1_{[-\varepsilon,\varepsilon]}; the Gaussian (heat) kernel ε−ne−π|x/ε|2\varepsilon^{-n}e^{-\pi|x/\varepsilon|^2}; the Poisson kernel Pε(x)=1πεx2+ε2P_\varepsilon(x)=\frac1\pi\frac{\varepsilon}{x^2+\varepsilon^2} on ℝ\mathbb R; and on the torus, the Fejér kernel of module 4. The Dirichlet kernel is not on the list, and that is the whole difficulty of Fourier series.

Explore: smoothing a function

Function f
Kernel K
ff and Kε*fK_\varepsilon*f
The kernel KεK_\varepsilon , total area 1

Key terms in plain words

The ideas above without the symbols. Switch to Quiz me to test yourself: you see an explanation and pick the term it describes.

Convolution
Replacing the value at each point by a weighted average of nearby values, with the weights given by a second function slid along the first.
Approximate identity
A family of averaging weights, each of total mass one, that concentrate more and more tightly around the origin, so that averaging with them gets closer and closer to doing nothing.
Young's inequality
The rule that tells you how integrable a smoothed function is, given how integrable the function and the smoothing weight are.
Mollifier
A smooth bump with compact support and total mass one, used to turn rough functions into infinitely differentiable ones by averaging.
Poisson kernel
The averaging weight that produces the harmonic function in the upper half-plane with given boundary values; its width is the height above the boundary.

Examples

Example 3.1 (moving averages).

With the box kernel, (Kε*f)(x)=12ε∫x−εx+εf(K_\varepsilon*f)(x)=\frac1{2\varepsilon}\int_{x-\varepsilon}^{x+\varepsilon}f. By the fundamental theorem of calculus it tends to f(x)f(x) at every point of continuity, and it is continuous even when ff is not.

Example 3.2 (the heat equation).

u(x,t)=(Gt*f)(x)u(x,t)=(G_t*f)(x) with Gt(x)=(4πt)−1/2e−x2/4tG_t(x)=(4\pi t)^{-1/2}e^{-x^2/4t} solves ut=uxxu_t=u_{xx} with u(⋅,0)=fu(\cdot,0)=f, the initial value attained in the sense of approximate identities as t→0t\to0.

Example 3.3 (the Poisson kernel).

1πyx2+y2=1πIm−1x+iy\frac1\pi\frac y{x^2+y^2}=\frac1\pi\,\mathrm{Im}\frac{-1}{x+iy} is harmonic in (x,y)(x,y), has integral 1 in xx, and is an approximate identity as y→0y\to0. So Py*fP_y*f is the harmonic extension of ff to the upper half-plane. Its slow decay, like 1/x21/x^2, shows in the figure as long tails.

Example 3.4 (smoothing).

If ϕ∈Cc∞\phi\in C_c^\infty with ∫ϕ=1\int\phi=1 and ϕε(x)=ε−nϕ(x/ε)\phi_\varepsilon(x)=\varepsilon^{-n}\phi(x/\varepsilon), then ϕε*f∈C∞\phi_\varepsilon*f\in C^\infty for every f∈Lloc1f\in L^1_{\mathrm{loc}}, and ϕε*f→f\phi_\varepsilon*f\to f in LpL^p. Hence Cc∞C_c^\infty is dense in LpL^p for p<∞p<\infty.

Example 3.5 (Weierstrass by convolution).

The Landau kernels Ln(x)=cn(1−x2)nL_n(x)=c_n(1-x^2)^n on [−1,1][-1,1], normalized to integral 1, form an approximate identity. For continuous ff on [0,1][0,1] vanishing at the ends, Ln*fL_n*f is a polynomial on [0,1][0,1] and converges uniformly to ff: the proof of the Weierstrass theorem behind the uniform convergence module.

Example 3.6 (jumps).

A continuous function cannot be uniformly close to a step function across its jump. So Kε*𝟏[0,∞)→𝟏[0,∞)K_\varepsilon*\mathbf 1_{[0,\infty)}\to\mathbf 1_{[0,\infty)} uniformly away from 0 but not on any neighbourhood of 0. At 0 itself a symmetric kernel gives exactly 12\frac12, the average of the two sides (Exercise 3.5).

Check your understanding

For f∈L1(ℝ)f\in L^1(\mathbb R) and g∈L2(ℝ)g\in L^2(\mathbb R), which bound holds?

  1. ‖f*g‖2≤‖f‖1‖g‖2\lVert f*g\rVert_2\le\lVert f\rVert_1\lVert g\rVert_2
  2. ‖f*g‖1≤‖f‖1‖g‖2\lVert f*g\rVert_1\le\lVert f\rVert_1\lVert g\rVert_2f*gf*g need not be integrable.
  3. ‖f*g‖∞≤‖f‖1‖g‖2\lVert f*g\rVert_\infty\le\lVert f\rVert_1\lVert g\rVert_2That needs g∈L∞g\in L^\infty, or f∈L2f\in L^2.
Minkowski's integral inequality: f*gf*g is an average of translates of gg with weights |f||f|, and each translate has the same L2L^2 norm.

Is there an f∈L1(ℝ)f\in L^1(\mathbb R) with f*g=gf*g=g for every g∈L1g\in L^1?

  1. Yes, the delta function.δ\delta is not a function (module 7).
  2. No: taking Fourier transforms, f^=1\hat f=1, contradicting the Riemann–Lebesgue lemma.
That is why approximate identities are families: the identity itself is not in L1L^1 (Exercise 3.8).

Exercises

Exercise 3.1Convolving with an integrable function

Prove ‖f*g‖p≤‖f‖1‖g‖p\lVert f*g\rVert_p\le\lVert f\rVert_1\lVert g\rVert_p for 1≤p≤∞1\le p\le\infty.

Hint

f*g=∫f(y)τygdyf*g=\int f(y)\,\tau_yg\,dy with τyg(x)=g(x−y)\tau_yg(x)=g(x-y), an LpL^p-valued integral; ‖τyg‖p=‖g‖p\lVert\tau_yg\rVert_p=\lVert g\rVert_p.

Solution

By Minkowski's integral inequality, ‖∫f(y)g(⋅−y)dy‖p≤∫|f(y)|‖g(⋅−y)‖pdy=‖f‖1‖g‖p\big\lVert\int f(y)g(\cdot-y)\,dy\big\rVert_p\le\int|f(y)|\,\lVert g(\cdot-y)\rVert_p\,dy=\lVert f\rVert_1\lVert g\rVert_p. For p=∞p=\infty it is direct: |f*g(x)|≤∫|f(y)|‖g‖∞dy|f*g(x)|\le\int|f(y)|\,\lVert g\rVert_\infty dy.

Exercise 3.2Uniform convergence

Let KεK_\varepsilon be an approximate identity and ff bounded and uniformly continuous. Show Kε*f→fK_\varepsilon*f\to f uniformly.

Hint

Kε*f(x)−f(x)=∫Kε(y)(f(x−y)−f(x))dyK_\varepsilon*f(x)-f(x)=\int K_\varepsilon(y)(f(x-y)-f(x))\,dy. Split at |y|=δ|y|=\delta.

Solution

Let A=sup⁡‖Kε‖1A=\sup\lVert K_\varepsilon\rVert_1. Given η>0\eta>0, choose δ\delta with |f(x−y)−f(x)|<η|f(x-y)-f(x)|<\eta for |y|<δ|y|<\delta and all xx. Then |Kε*f(x)−f(x)|≤Aη+2‖f‖∞∫|y|≥δ|Kε||K_\varepsilon*f(x)-f(x)|\le A\eta+2\lVert f\rVert_\infty\int_{|y|\ge\delta}|K_\varepsilon|, uniformly in xx, and the last integral tends to 0.

Exercise 3.3LpL^p convergence

Show Kε*f→fK_\varepsilon*f\to f in LpL^p for f∈Lpf\in L^p, 1≤p<∞1\le p<\infty.

Hint

For f∈Ccf\in C_c, use Exercise 3.2 plus a uniform bound on the tails. For general ff, approximate by g∈Ccg\in C_c and use Exercise 3.1: ‖Kε*(f−g)‖p≤A‖f−g‖p\lVert K_\varepsilon*(f-g)\rVert_p\le A\lVert f-g\rVert_p.

Exercise 3.4The Poisson kernel

Show ∫Pε=1\int P_\varepsilon=1 and that PεP_\varepsilon is an approximate identity.

Hint

∫εdxx2+ε2=arctan⁡(x/ε)|−∞∞=π\int\frac{\varepsilon\,dx}{x^2+\varepsilon^2}=\arctan(x/\varepsilon)\big|_{-\infty}^\infty=\pi. Tails: ∫|x|≥δPε=1−2πarctan⁡(δ/ε)→0\int_{|x|\ge\delta}P_\varepsilon=1-\frac2\pi\arctan(\delta/\varepsilon)\to0.

Exercise 3.5Symmetric kernels at a jump

If KεK_\varepsilon is an even approximate identity and ff has one-sided limits f(0±)f(0^\pm) at 0, show Kε*f(0)→12(f(0+)+f(0−))K_\varepsilon*f(0)\to\frac12(f(0^+)+f(0^-)).

Hint

Evenness gives ∫0∞Kε=∫−∞0Kε=12\int_0^\infty K_\varepsilon=\int_{-\infty}^0K_\varepsilon=\frac12. Treat the two half-lines separately, as in Exercise 3.2.

Exercise 3.6Convolution of LpL^p and Lp′L^{p'} is continuous

Show that for f∈Lpf\in L^p, g∈Lp′g\in L^{p'}, 1≤p<∞1\le p<\infty, f*gf*g is uniformly continuous.

Hint

|f*g(x+h)−f*g(x)|≤‖τ−hf−f‖p‖g‖p′|f*g(x+h)-f*g(x)|\le\lVert\tau_{-h}f-f\rVert_p\lVert g\rVert_{p'}, and translation is continuous in LpL^p for p<∞p<\infty (approximate by CcC_c).

Exercise 3.7Weierstrass via Landau kernels

Carry out Example 3.5: show LnL_n is an approximate identity on [−1,1][-1,1], and that Ln*fL_n*f restricted to [0,1][0,1] is a polynomial when ff is continuous on [0,1][0,1], zero at the ends, and extended by 0.

Hint

cn≤nc_n\le\sqrt n (bound (1−x2)n≥1−nx2(1-x^2)^n\ge1-nx^2 on |x|≤1/n|x|\le1/\sqrt n). For δ≤|x|≤1\delta\le|x|\le1, Ln≤n(1−δ2)n→0L_n\le\sqrt n(1-\delta^2)^n\to0. For the polynomial: ∫01f(y)(1−(x−y)2)ndy\int_0^1f(y)(1-(x-y)^2)^ndy is a polynomial in xx. This is Rudin's proof of Theorem 7.26.

Solution

Since (1−x2)n≥1−nx2(1-x^2)^n\ge1-nx^2, ∫−11(1−x2)n≥2∫01/n(1−nx2)=43n>1n\int_{-1}^1(1-x^2)^n\ge2\int_0^{1/\sqrt n}(1-nx^2)=\frac4{3\sqrt n}>\frac1{\sqrt n}, so cn<nc_n<\sqrt n. The kernels are positive with integral 1, and supδ≤|x|≤1⁡Ln≤n(1−δ2)n→0\sup_{\delta\le|x|\le1}L_n\le\sqrt n(1-\delta^2)^n\to0, so the tails vanish. For x∈[0,1]x\in[0,1] and y∈[0,1]y\in[0,1], x−y∈[−1,1]x-y\in[-1,1], so Ln*f(x)=cn∫01f(y)(1−(x−y)2)ndyL_n*f(x)=c_n\int_0^1f(y)(1-(x-y)^2)^ndy, and expanding (1−(x−y)2)n(1-(x-y)^2)^n shows it is a polynomial in xx. Exercise 3.2 gives uniform convergence (the extension of ff by 0 is uniformly continuous because f(0)=f(1)=0f(0)=f(1)=0).

Exercise 3.8No identity element

Show there is no f∈L1(ℝ)f\in L^1(\mathbb R) with f*g=gf*g=g for all g∈L1g\in L^1.

Hint

Take gg a Gaussian and Fourier transform: f^g^=g^\hat f\hat g=\hat g forces f^≡1\hat f\equiv1, but f^(ξ)→0\hat f(\xi)\to0 (Exercise 6.8).

InteractiveGrafakos, Classical, §§3.1.3, 3.2.1, 3.4, 3.5.4

Fourier series: Dirichlet, Fejér and Gibbs

The partial sums of a Fourier series are convolution with the Dirichlet kernel, which is not an approximate identity: its L1L^1 norm grows like log⁡N\log N. Averaging the partial sums gives the Fejér kernel, which is one. That single difference explains Fejér's theorem, the Gibbs overshoot, and why pointwise convergence of Fourier series is hard.

The idea

Coefficients and partial sums (Grafakos §§3.1.2, 3.2.1). For f∈L1(𝕋)f\in L^1(\mathbb T), 𝕋=[0,1]\mathbb T=[0,1]: f^(m)=∫01f(x)e−2πimxdx\hat f(m)=\int_0^1f(x)e^{-2\pi imx}dx, and SNf(x)=∑|m|≤Nf^(m)e2πimx=(DN*f)(x)S_Nf(x)=\sum_{|m|\le N}\hat f(m)e^{2\pi imx}=(D_N*f)(x) with the Dirichlet kernel
DN(x)=∑|m|≤Ne2πimx=sin⁡((2N+1)πx)sin⁡(πx).D_N(x)=\sum_{|m|\le N}e^{2\pi imx}=\frac{\sin((2N+1)\pi x)}{\sin(\pi x)}.
(Grafakos uses [0,1][0,1] and e2πimxe^{2\pi imx}; the first module used eikθe^{ik\theta} on [0,2π][0,2\pi]. Only the bookkeeping differs.)
Fejér's theorem (Grafakos §3.1.3). The Cesàro means σNf=1N+1∑j=0NSjf=FN*f\sigma_Nf=\frac1{N+1}\sum_{j=0}^NS_jf=F_N*f use the Fejér kernel
FN(x)=1N+1(sin⁡((N+1)πx)sin⁡(πx))2≥0,F_N(x)=\frac1{N+1}\Big(\frac{\sin((N+1)\pi x)}{\sin(\pi x)}\Big)^2\ge0,
an approximate identity. So σNf→f\sigma_Nf\to f uniformly for continuous ff, and in LpL^p for f∈Lpf\in L^p, p<∞p<\infty. Consequences: trigonometric polynomials are dense, and a continuous function is determined by its Fourier coefficients.
Why partial sums are harder (Grafakos §§3.4.3–3.4.4). ‖DN‖L1∼4π2log⁡N\lVert D_N\rVert_{L^1}\sim\frac4{\pi^2}\log N (the Lebesgue constants), so the uniform boundedness principle produces continuous functions whose Fourier series diverge at a point. Positive results need more regularity: if ff is Hölder continuous at x0x_0, or merely satisfies Dini's condition there, then SNf(x0)→f(x0)S_Nf(x_0)\to f(x_0). Carleson's theorem (for L2L^2; Hunt for LpL^p, p>1p>1) gives almost-everywhere convergence.
The Gibbs phenomenon (Grafakos §3.5.4). At a jump of size JJ, SNfS_Nf overshoots by about 0.0895J0.0895J on each side, however large NN is: limN⁡max⁡SN(square wave)=2π∫0πsin⁡ttdt≈1.17898\lim_N\max S_N(\text{square wave})=\frac2\pi\int_0^\pi\frac{\sin t}t\,dt\approx1.17898. The overshoot moves closer to the jump but does not shrink. Fejér means never overshoot, because FN≥0F_N\ge0.

Explore: partial sums against Fejér means

Function f on [0, 1]
Show kernel
ff , partial sum SNfS_Nf , Fejér mean σNf\sigma_Nf
The chosen kernel on [−12,12][-\frac12,\frac12]

Key terms in plain words

The ideas above without the symbols. Switch to Quiz me to test yourself: you see an explanation and pick the term it describes.

Dirichlet kernel
The sum of all pure frequencies up to a cutoff: averaging against it gives the truncated Fourier series. It oscillates and takes negative values, and its total absolute area grows without bound.
Fejér kernel
The average of the first few truncation kernels: a non-negative bump of area one that narrows as more terms are used.
Cesàro mean
The average of the first NN partial sums of a series, which can converge even when the partial sums themselves don't.
Lebesgue constants
The total absolute areas of the truncation kernels, which grow like the logarithm of the number of terms; their growth is why some continuous functions have divergent Fourier series.
Gibbs phenomenon
The overshoot of about nine percent of the jump that truncated Fourier series make next to a jump discontinuity, which narrows but never shrinks.
Dini's test
A mild condition on how a function behaves near one point, which is enough for its Fourier series to converge to the right value at that point.

Examples

Example 4.1 (the square wave).

f=1f=1 on (0,12)(0,\frac12), −1-1 on (12,1)(\frac12,1) has f^(m)=2πim\hat f(m)=\frac{2}{\pi im} for odd mm and 0 for even m≠0m\ne0, so f(x)∼4π∑k≥0sin⁡(2π(2k+1)x)2k+1f(x)\sim\frac4\pi\sum_{k\ge0}\frac{\sin(2\pi(2k+1)x)}{2k+1}. At x=14x=\frac14 this is Leibniz's series: 1=4π(1−13+15−⋯)1=\frac4\pi(1-\frac13+\frac15-\cdots).

Example 4.2 (the Dirichlet kernel changes sign).

DN(0)=2N+1D_N(0)=2N+1, ∫𝕋DN=1\int_{\mathbb T}D_N=1, but DND_N oscillates with zeros at x=k2N+1x=\frac k{2N+1}. Its positive and negative lobes both have total area of order log⁡N\log N, cancelling except for the net 1.

Example 4.3 (the Fejér kernel).

FN≥0F_N\ge0, ∫FN=1\int F_N=1, and FN(x)≤1(N+1)sin2⁡(πδ)F_N(x)\le\frac1{(N+1)\sin^2(\pi\delta)} for δ≤|x|≤12\delta\le|x|\le\frac12: the three conditions of an approximate identity, visible in the right panel.

Example 4.4 (Gibbs, numerically).

For the square wave, max⁡S9f≈1.1823\max S_9f\approx1.1823 and max⁡S59f≈1.1791\max S_{59}f\approx1.1791, against the limit 1.178981.17898: about 9%9\% of the jump of 2 on each side, for every large NN. The sawtooth shows the same 9%9\% at its single jump; the continuous triangle wave shows none.

Example 4.5 (Weierstrass on the circle).

Fejér's theorem says every continuous periodic function is a uniform limit of trigonometric polynomials, the periodic version of the Weierstrass theorem in the uniform convergence module.

Example 4.6 (convergence at a smooth point).

The square wave is smooth on (0,12)(0,\frac12), so Dini's test gives SNf(x)→1S_Nf(x)\to1 there, and at the jump SNf(0)=0S_Nf(0)=0 for all NN, the average of the two sides. Convergence holds everywhere, but not uniformly.

Check your understanding

Why is the Fejér kernel an approximate identity while the Dirichlet kernel is not?

  1. FN≥0F_N\ge0 has ‖FN‖1=1\lVert F_N\rVert_1=1, while ‖DN‖1→∞\lVert D_N\rVert_1\to\infty.
  2. ∫DN≠1\int D_N\ne1.Both integrate to 1.
  3. DND_N is not even.Both are even.
An approximate identity needs sup⁡‖K‖1<∞\sup\lVert K\rVert_1<\infty; the Lebesgue constants ‖DN‖1\lVert D_N\rVert_1 grow like log⁡N\log N.

As N→∞N\to\infty, the Gibbs overshoot of SNfS_Nf at a jump:

  1. tends to 0.Increase N in the figure.
  2. tends to about 9%9\% of the jump, while moving closer to it.
  3. grows like log⁡N\log N.That's ‖DN‖1\lVert D_N\rVert_1, not the overshoot.
The overshoot converges to (1π∫0πsin⁡ttdt−12)J≈0.0895J\big(\frac1\pi\int_0^\pi\frac{\sin t}t\,dt-\frac12\big)J\approx0.0895J.

Exercises

Exercise 4.1The Dirichlet kernel

Prove the closed form of DND_N.

Hint

Geometric series: ∑m=−NNzm=zN+1−z−Nz−1\sum_{m=-N}^Nz^m=\frac{z^{N+1}-z^{-N}}{z-1} with z=e2πixz=e^{2\pi ix}; multiply top and bottom by z−1/2z^{-1/2}.

Solution

With z=e2πixz=e^{2\pi ix}, DN=zN+1−z−Nz−1=zN+1/2−z−N−1/2z1/2−z−1/2=2isin⁡((2N+1)πx)2isin⁡(πx)D_N=\frac{z^{N+1}-z^{-N}}{z-1}=\frac{z^{N+1/2}-z^{-N-1/2}}{z^{1/2}-z^{-1/2}}=\frac{2i\sin((2N+1)\pi x)}{2i\sin(\pi x)}.

Exercise 4.2The Fejér kernel

Prove the closed form of FNF_N, and that FNF_N is an approximate identity on 𝕋\mathbb T.

Hint

sin⁡(πx)Dj(x)=sin⁡((2j+1)πx)\sin(\pi x)\,D_j(x)=\sin((2j+1)\pi x), and 2sin⁡(πx)sin⁡((2j+1)πx)=cos⁡(2jπx)−cos⁡(2(j+1)πx)2\sin(\pi x)\sin((2j+1)\pi x)=\cos(2j\pi x)-\cos(2(j+1)\pi x) telescopes.

Solution

2sin2⁡(πx)∑j=0NDj=∑j(cos⁡(2jπx)−cos⁡(2(j+1)πx))=1−cos⁡(2(N+1)πx)=2sin2⁡((N+1)πx)2\sin^2(\pi x)\sum_{j=0}^ND_j=\sum_j\big(\cos(2j\pi x)-\cos(2(j+1)\pi x)\big)=1-\cos(2(N+1)\pi x)=2\sin^2((N+1)\pi x). Divide by 2(N+1)sin2⁡(πx)2(N+1)\sin^2(\pi x). Then FN≥0F_N\ge0, ∫FN=1N+1∑∫Dj=1\int F_N=\frac1{N+1}\sum\int D_j=1, and on δ≤|x|≤12\delta\le|x|\le\frac12, FN≤1(N+1)sin2⁡(πδ)→0F_N\le\frac1{(N+1)\sin^2(\pi\delta)}\to0.

Exercise 4.3Fejér's theorem and uniqueness

Deduce that σNf→f\sigma_Nf\to f uniformly for continuous ff on 𝕋\mathbb T, and that if all f^(m)=0\hat f(m)=0 for a continuous ff, then f=0f=0.

Hint

Exercises 4.2 and 3.2 (on the torus, continuous means uniformly continuous). If all coefficients vanish, every σNf\sigma_Nf is 0.

Exercise 4.4Lebesgue constants

Show ‖DN‖L1(𝕋)≥clog⁡N\lVert D_N\rVert_{L^1(\mathbb T)}\ge c\log N for some c>0c>0.

Hint

Use |sin⁡πx|≤π|x||\sin\pi x|\le\pi|x| and integrate |sin⁡((2N+1)πx)|/(π|x|)|\sin((2N+1)\pi x)|/(\pi|x|) over the intervals [k2N+1,k+12N+1]\big[\frac k{2N+1},\frac{k+1}{2N+1}\big], k=0,…,N−1k=0,\dots,N-1. Each contributes about ck+1\frac c{k+1}.

Solution

‖DN‖1≥2∫01/2|sin⁡((2N+1)πx)|πxdx\lVert D_N\rVert_1\ge2\int_0^{1/2}\frac{|\sin((2N+1)\pi x)|}{\pi x}dx. On Ik=[k2N+1,k+12N+1]I_k=\big[\frac k{2N+1},\frac{k+1}{2N+1}\big], x≤k+12N+1x\le\frac{k+1}{2N+1}, and ∫Ik|sin⁡((2N+1)πx)|dx=2π(2N+1)\int_{I_k}|\sin((2N+1)\pi x)|dx=\frac2{\pi(2N+1)}. So the integral over IkI_k is at least 2π2(k+1)\frac{2}{\pi^2(k+1)}, and summing over k<Nk<N gives 4π2∑k≤N1k≥4π2log⁡N\frac4{\pi^2}\sum_{k\le N}\frac1k\ge\frac4{\pi^2}\log N.

Exercise 4.5Coefficients of the square wave

Verify the coefficients in Example 4.1, and derive Leibniz's series.

Hint

f^(m)=∫01/2e−2πimxdx−∫1/21e−2πimxdx=1−(−1)mπim\hat f(m)=\int_0^{1/2}e^{-2\pi imx}dx-\int_{1/2}^1e^{-2\pi imx}dx=\frac{1-(-1)^m}{\pi im}. Use Dini's test at x=14x=\frac14.

Exercise 4.6The Gibbs constant

Show that the first maximum of S2K−1fS_{2K-1}f for the square wave is at x=14Kx=\frac1{4K}, with value tending to 2π∫0πsin⁡ttdt\frac2\pi\int_0^\pi\frac{\sin t}t\,dt.

Hint

ddxS2K−1f=8∑k<Kcos⁡(2π(2k+1)x)=4sin⁡(4πKx)sin⁡(2πx)\frac d{dx}S_{2K-1}f=8\sum_{k<K}\cos(2\pi(2k+1)x)=\frac{4\sin(4\pi Kx)}{\sin(2\pi x)}, first zero at 14K\frac1{4K}. The value there is a Riemann sum for 2π∫0πsin⁡ttdt\frac2\pi\int_0^\pi\frac{\sin t}t\,dt.

Exercise 4.7Dini's test

Show that if ∫−1/21/2|f(x0−t)−f(x0)||t|dt<∞\int_{-1/2}^{1/2}\frac{|f(x_0-t)-f(x_0)|}{|t|}dt<\infty, then SNf(x0)→f(x0)S_Nf(x_0)\to f(x_0).

Hint

SNf(x0)−f(x0)=∫g(t)sin⁡((2N+1)πt)dtS_Nf(x_0)-f(x_0)=\int g(t)\sin((2N+1)\pi t)\,dt with g(t)=f(x0−t)−f(x0)sin⁡πt∈L1g(t)=\frac{f(x_0-t)-f(x_0)}{\sin\pi t}\in L^1. Apply the Riemann–Lebesgue lemma (module 5).

InteractiveGrafakos, Classical, §§3.2.2, 3.3

Decay, smoothness and L²

How fast Fourier coefficients decay measures how smooth the function is: a jump gives 1/m1/m, a corner 1/m21/m^2, and an analytic function decays exponentially. In L2L^2 there is no loss at all: the coefficients are the coordinates of ff in an orthonormal basis, exactly as for the finite groups of the first module.

The idea

Riemann–Lebesgue lemma (Grafakos §3.3.1). For f∈L1(𝕋)f\in L^1(\mathbb T), f^(m)→0\hat f(m)\to0 as |m|→∞|m|\to\infty. Nothing better holds in general: for any sequence am→0a_m\to0 there is an f∈L1f\in L^1 whose coefficients do not decay faster than ama_m.
Smoothness gives decay (Grafakos §3.3.2). Integrating by parts, f′^(m)=2πimf^(m)\widehat{f'}(m)=2\pi im\hat f(m), so if f∈Ck(𝕋)f\in C^k(\mathbb T) then |f^(m)|≤‖f(k)‖1(2π|m|)k|\hat f(m)|\le\frac{\lVert f^{(k)}\rVert_1}{(2\pi|m|)^k}, indeed o(|m|−k)o(|m|^{-k}). Conversely, if ∑|m|k|f^(m)|<∞\sum|m|^k|\hat f(m)|<\infty, the series can be differentiated kk times term by term and f∈Ckf\in C^k. A function of bounded variation has f^(m)=O(1/|m|)\hat f(m)=O(1/|m|).
Plancherel on the circle (Grafakos §3.2.2). {e2πimx}m∈ℤ\{e^{2\pi imx}\}_{m\in\mathbb Z} is an orthonormal basis of L2(𝕋)L^2(\mathbb T): orthonormality is a direct computation, completeness follows from Fejér's theorem. Hence ‖f‖22=∑m|f^(m)|2\lVert f\rVert_2^2=\sum_m|\hat f(m)|^2, SNf→fS_Nf\to f in L2L^2, and f↦f^f\mapsto\hat f is a unitary isomorphism L2(𝕋)→ℓ2(ℤ)L^2(\mathbb T)\to\ell^2(\mathbb Z). This is the infinite version of the finite Plancherel identity.

Explore: decay rates on log–log axes

|f^(m)||\hat f(m)| for m=1,…,256m=1,\dots,256 on log–log axes , with reference lines m−1m^{-1} to m−4m^{-4} . A straight line of slope −k-k means decay like m−km^{-k}.
Function

Key terms in plain words

The ideas above without the symbols. Switch to Quiz me to test yourself: you see an explanation and pick the term it describes.

Riemann–Lebesgue lemma
The statement that high-frequency content of any integrable function fades out, though possibly very slowly.
Bounded variation
Describes a function whose graph has finite total up-and-down movement, so its frequency content decays at least like one over the frequency.
Orthonormal basis
A list of mutually perpendicular unit vectors such that every vector is the (possibly infinite) sum of its components along them.
Parseval's identity
The energy of a periodic function equals the sum of the squared sizes of its frequency coefficients.
Absolutely convergent Fourier series
A series whose coefficient sizes have a finite sum, which then converges uniformly to a continuous function.

Examples

Example 5.1 (a jump).

The square wave has |f^(m)|=2π|m||\hat f(m)|=\frac2{\pi|m|} for odd mm: slope −1-1, the rate for functions of bounded variation with a jump.

Example 5.2 (a corner).

The triangle wave (the integral of a square wave) has |f^(m)|=cm2|\hat f(m)|=\frac{c}{m^2} for odd mm: one integration, one more power, by f′^(m)=2πimf^(m)\widehat{f'}(m)=2\pi im\hat f(m).

Example 5.3 (analytic).

ecos⁡2πxe^{\cos2\pi x} has f^(m)=Im(1)\hat f(m)=I_m(1), a modified Bessel function, which decays like 12mm!\frac1{2^m m!}: faster than any power. Analytic periodic functions have exponentially decaying coefficients, which is the trapezoid-rule phenomenon of the UTokyo analysis module.

Example 5.4 (Basel by Parseval).

The sawtooth f(x)=x−12f(x)=x-\frac12 on [0,1)[0,1) has f^(m)=i2πm\hat f(m)=\frac{i}{2\pi m} for m≠0m\ne0. Parseval: ∫01(x−12)2dx=112=∑m≠014π2m2\int_0^1(x-\frac12)^2dx=\frac1{12}=\sum_{m\ne0}\frac1{4\pi^2m^2}, so ∑m≥11m2=π26\sum_{m\ge1}\frac1{m^2}=\frac{\pi^2}6.

Example 5.5 (a cusp).

|x−12|1/2|x-\frac12|^{1/2} is continuous and monotone on each half, with a cusp at 12\frac12 where the derivative blows up. Its coefficients decay like |m|−3/2|m|^{-3/2}, between the jump and the corner: fractional smoothness gives fractional decay.

Example 5.6 (from decay back to smoothness).

If |f^(m)|≤C|m|−4|\hat f(m)|\le C|m|^{-4}, then ∑|m|2|f^(m)|<∞\sum|m|^2|\hat f(m)|<\infty, so f∈C2f\in C^2. The figure's x2(1−x)2x^2(1-x)^2 has exactly this decay: it is C2C^2, and the jump in its third derivative costs the fourth power.

Check your understanding

A periodic function has |f^(m)|≈m−2|\hat f(m)|\approx m^{-2}. What can you say?

  1. The series converges absolutely and uniformly, so ff is continuous; a corner (a jump in f′f') is typical of this rate.
  2. ff must have a jump.Jumps give 1/m1/m.
  3. ff is C1C^1.∑|m|m−2\sum|m|\,m^{-2} diverges; the triangle wave is not C1C^1.
∑m−2<∞\sum m^{-2}<\infty gives uniform convergence (the Weierstrass M-test from the uniform convergence module).

For f∈L2(𝕋)f\in L^2(\mathbb T), SNf→fS_Nf\to f:

  1. in L2L^2, always.
  2. uniformly.The square wave is in L2L^2.
  3. only if ff is continuous.Orthonormal expansions converge in norm for every vector.
‖f−SNf‖22=∑|m|>N|f^(m)|2→0\lVert f-S_Nf\rVert_2^2=\sum_{|m|>N}|\hat f(m)|^2\to0. Almost-everywhere convergence is also true but is Carleson's deep theorem.

Exercises

Exercise 5.1The Riemann–Lebesgue lemma

Prove f^(m)→0\hat f(m)\to0 for f∈L1(𝕋)f\in L^1(\mathbb T).

Hint

It holds for trigonometric polynomials (eventually 0). Approximate ff in L1L^1 by one, using Fejér, and use |f^(m)−g^(m)|≤‖f−g‖1|\hat f(m)-\hat g(m)|\le\lVert f-g\rVert_1.

Solution

Given ε\varepsilon, Fejér's theorem in L1L^1 gives a trigonometric polynomial gg of degree MM with ‖f−g‖1<ε\lVert f-g\rVert_1<\varepsilon. For |m|>M|m|>M, g^(m)=0\hat g(m)=0, so |f^(m)|=|f−g^(m)|≤‖f−g‖1<ε|\hat f(m)|=|\widehat{f-g}(m)|\le\lVert f-g\rVert_1<\varepsilon.

Exercise 5.2Differentiation multiplies by 2πim2\pi im

Show f′^(m)=2πimf^(m)\widehat{f'}(m)=2\pi im\hat f(m) for f∈C1(𝕋)f\in C^1(\mathbb T), and deduce f^(m)=o(|m|−k)\hat f(m)=o(|m|^{-k}) for f∈Ckf\in C^k.

Hint

Integrate by parts; periodicity kills the boundary term. Then apply Riemann–Lebesgue to f(k)f^{(k)}.

Solution

f′^(m)=[f(x)e−2πimx]01+2πim∫01f(x)e−2πimxdx=2πimf^(m)\widehat{f'}(m)=\big[f(x)e^{-2\pi imx}\big]_0^1+2\pi im\int_0^1f(x)e^{-2\pi imx}dx=2\pi im\hat f(m). Iterating, f^(m)=f(k)^(m)(2πim)k\hat f(m)=\frac{\widehat{f^{(k)}}(m)}{(2\pi im)^k} for m≠0m\ne0, and f(k)^(m)→0\widehat{f^{(k)}}(m)\to0.

Exercise 5.3Completeness of the exponentials

Show that {e2πimx}\{e^{2\pi imx}\} is an orthonormal basis of L2(𝕋)L^2(\mathbb T), and deduce Parseval's identity.

Hint

Orthonormal: integrate e2πi(m−n)xe^{2\pi i(m-n)x}. Complete: continuous functions are dense in L2L^2, and Fejér approximates each one uniformly (hence in L2L^2) by trigonometric polynomials.

Exercise 5.4Absolutely convergent series

Show that if ∑|f^(m)|<∞\sum|\hat f(m)|<\infty for f∈L1f\in L^1, then SNfS_Nf converges uniformly to a continuous function equal to ff almost everywhere.

Hint

The M-test gives a continuous limit gg with g^=f^\hat g=\hat f. Uniqueness of coefficients in L1L^1 (via Fejér in L1L^1) gives f=gf=g a.e.

Exercise 5.5C1C^1 functions have absolutely convergent series

Show ∑m|f^(m)|<∞\sum_m|\hat f(m)|<\infty for f∈C1(𝕋)f\in C^1(\mathbb T). (So Fourier series of C1C^1 functions converge uniformly.)

Hint

∑m≠0|f^(m)|=∑|f′^(m)|2π|m|≤(∑14π2m2)1/2(∑|f′^(m)|2)1/2\sum_{m\ne0}|\hat f(m)|=\sum\frac{|\widehat{f'}(m)|}{2\pi|m|}\le\big(\sum\frac1{4\pi^2m^2}\big)^{1/2}\big(\sum|\widehat{f'}(m)|^2\big)^{1/2}, and Parseval for f′f'.

Exercise 5.6Bounded variation

Show that if ff has bounded variation VV on [0,1][0,1], then |f^(m)|≤V2π|m||\hat f(m)|\le\frac{V}{2\pi|m|}.

Hint

Integrate by parts in the Riemann–Stieltjes sense (the Riemann–Stieltjes module on the Real Analysis page): f^(m)=12πim∫01e−2πimxdf(x)\hat f(m)=\frac1{2\pi im}\int_0^1e^{-2\pi imx}\,df(x), plus a boundary term that vanishes for periodic ff.

Exercise 5.7Basel and beyond

Verify Example 5.4, then apply Parseval to the triangle wave to compute ∑k≥01(2k+1)4=π496\sum_{k\ge0}\frac1{(2k+1)^4}=\frac{\pi^4}{96}.

Hint

For the sawtooth, integrate by parts. For the triangle wave t(x)=|x−12|t(x)=|x-\frac12| on [0,1][0,1] you get |t^(m)|=1π2m2|\hat t(m)|=\frac1{\pi^2m^2} for odd mm; compare ∫01(t−14)2=148\int_0^1(t-\frac14)^2=\frac1{48} with 2∑odd m>01π4m42\sum_{\text{odd }m>0}\frac1{\pi^4m^4}.

InteractiveGrafakos, Classical, §§2.2, 3.2.3

The Fourier transform on ℝⁿ

On ℝn\mathbb R^n the frequencies form a continuum, and sums become integrals. The Schwartz functions, smooth and rapidly decaying, are the setting in which every formula holds without fuss: the transform swaps differentiation with multiplication, translation with modulation, and the Gaussian is its own transform.

The idea

Schwartz functions and the transform (Grafakos §§2.2.1–2.2.2). 𝒮(ℝn)\mathcal S(\mathbb R^n) is the space of smooth ff with xα∂βfx^\alpha\partial^\beta f bounded for all multi-indices. For f∈L1f\in L^1, f^(ξ)=∫f(x)e−2πix⋅ξdx\hat f(\xi)=\int f(x)e^{-2\pi ix\cdot\xi}dx, also written ℱ(f)\mathcal F(f). The transform maps 𝒮\mathcal S to itself, with ∂jf^=2πiξjf^\widehat{\partial_jf}=2\pi i\xi_j\hat f, xjf^=−12πi∂jf^\widehat{x_jf}=-\frac1{2\pi i}\partial_j\hat f, f(⋅−a)^=e−2πia⋅ξf^\widehat{f(\cdot-a)}=e^{-2\pi ia\cdot\xi}\hat f, and f(λ⋅)^=λ−nf^(⋅/λ)\widehat{f(\lambda\cdot)}=\lambda^{-n}\hat f(\cdot/\lambda).
Inversion and Plancherel (Grafakos §§2.2.3–2.2.4). e−π|x|2e^{-\pi|x|^2} is its own transform. Using Gaussians as an approximate identity one proves f(x)=∫f^(ξ)e2πix⋅ξdξf(x)=\int\hat f(\xi)e^{2\pi ix\cdot\xi}d\xi for f∈𝒮f\in\mathcal S, and ‖f^‖2=‖f‖2\lVert\hat f\rVert_2=\lVert f\rVert_2. The transform extends to a unitary map of L2(ℝn)L^2(\mathbb R^n), and f*g^=f^g^\widehat{f*g}=\hat f\hat g.
Poisson summation (Grafakos §3.2.3). For nice ff on ℝn\mathbb R^n, ∑k∈ℤnf(x+k)=∑m∈ℤnf^(m)e2πim⋅x\sum_{k\in\mathbb Z^n}f(x+k)=\sum_{m\in\mathbb Z^n}\hat f(m)e^{2\pi im\cdot x}: periodizing ff and sampling f^\hat f are the same operation seen from the two sides. At x=0x=0, ∑f(k)=∑f^(m)\sum f(k)=\sum\hat f(m). This links the transform on ℝ\mathbb R, the series on 𝕋\mathbb T, and theta functions.
Uncertainty. For f∈𝒮(ℝ)f\in\mathcal S(\mathbb R), ‖xf‖2‖ξf^‖2≥14π‖f‖22\lVert xf\rVert_2\lVert\xi\hat f\rVert_2\ge\frac1{4\pi}\lVert f\rVert_2^2, with equality exactly for Gaussians (Exercise 6.6). A function and its transform cannot both be concentrated.

Explore: Gaussians, widths and Poisson summation

Ref\mathrm{Re}\,f and |f||f| ; ticks at the integers, where ∑f(k)\sum f(k) samples it
f^(ξ)=σe−πσ2(ξ−a)2\hat f(\xi)=\sigma e^{-\pi\sigma^2(\xi-a)^2} ; ticks at the integers, where ∑f^(m)\sum\hat f(m) samples it

Key terms in plain words

The ideas above without the symbols. Switch to Quiz me to test yourself: you see an explanation and pick the term it describes.

Schwartz function
A smooth function that, together with all its derivatives, decays faster than any power of 1/|x|1/|x|: the class where every manipulation of the Fourier transform is safe.
Modulation
Multiplying a function by a pure oscillation e2πiaxe^{2\pi iax}, which shifts its frequency content by aa.
Fourier inversion
Recovering a function by superposing its frequency components with the conjugate exponentials.
Poisson summation formula
The sum of a function's values at the integers equals the sum of its transform's values at the integers.
Uncertainty principle
A function and its frequency content can't both be sharply concentrated: the product of their spreads has a positive lower bound.

Examples

Example 6.1 (the Gaussian).

e−πx2^=e−πξ2\widehat{e^{-\pi x^2}}=e^{-\pi\xi^2}, and by dilation e−πx2/σ2^=σe−πσ2ξ2\widehat{e^{-\pi x^2/\sigma^2}}=\sigma e^{-\pi\sigma^2\xi^2}: width σ\sigma on one side, 1/σ1/\sigma on the other.

Example 6.2 (box and sinc).

ℱ(𝟏[−1/2,1/2])(ξ)=sin⁡πξπξ\mathcal F\big(\mathbf 1_{[-1/2,1/2]}\big)(\xi)=\frac{\sin\pi\xi}{\pi\xi}. It is in L2L^2 but not L1L^1, so its inverse transform exists only in the L2L^2 sense, and the box's jumps reappear as the slow 1/ξ1/\xi decay.

Example 6.3 (the two-sided exponential).

e−2π|x|^(ξ)=1π11+ξ2\widehat{e^{-2\pi|x|}}(\xi)=\frac1\pi\frac1{1+\xi^2}: a corner at 0 gives ξ−2\xi^{-2} decay, and the transform of the Poisson kernel is e−2π|ξ|e^{-2\pi|\xi|} by inversion.

Example 6.4 (the heat equation).

Transforming ut=Δuu_t=\Delta u in xx gives u^t=−4π2|ξ|2u^\hat u_t=-4\pi^2|\xi|^2\hat u, so u^(ξ,t)=e−4π2t|ξ|2f^(ξ)\hat u(\xi,t)=e^{-4\pi^2t|\xi|^2}\hat f(\xi): high frequencies die fastest, and u=Gt*fu=G_t*f with the Gaussian of Example 3.2.

Example 6.5 (theta functions).

Poisson summation for e−πtx2e^{-\pi tx^2} gives ∑ke−πtk2=t−1/2∑me−πm2/t\sum_ke^{-\pi tk^2}=t^{-1/2}\sum_me^{-\pi m^2/t}, i.e. θ(1/t)=tθ(t)\theta(1/t)=\sqrt t\,\theta(t). This modular transformation is the root of the functional equation of the zeta function, and of the Modular Forms page.

Example 6.6 (sampling).

If f^\hat f vanishes outside [−12,12][-\frac12,\frac12], Poisson summation and inversion give f(x)=∑kf(k)sin⁡π(x−k)π(x−k)f(x)=\sum_kf(k)\frac{\sin\pi(x-k)}{\pi(x-k)} (Shannon): a band-limited function is determined by its integer samples.

Check your understanding

ff is a Gaussian of width σ\sigma. Halving σ\sigma makes f^\hat f:

  1. twice as wide and half as tall.
  2. half as wide.Width σ\sigma on one side means width 1/σ1/\sigma on the other.
  3. unchanged.Try it in the figure.
f(λ⋅)^=λ−1f^(⋅/λ)\widehat{f(\lambda\cdot)}=\lambda^{-1}\hat f(\cdot/\lambda) with λ=2\lambda=2.

What does multiplying ff by e2πiaxe^{2\pi iax} do to f^\hat f?

  1. Shifts it by aa.
  2. Multiplies it by e2πiaξe^{2\pi ia\xi}.That's what translating ff does (with a minus sign).
∫f(x)e2πiaxe−2πixξdx=f^(ξ−a)\int f(x)e^{2\pi iax}e^{-2\pi ix\xi}dx=\hat f(\xi-a). Translation and modulation are swapped by the transform.

Exercises

Exercise 6.1The Gaussian is its own transform

Prove e−πx2^=e−πξ2\widehat{e^{-\pi x^2}}=e^{-\pi\xi^2}.

Hint

Both sides solve y′=−2πξyy'=-2\pi\xi y with y(0)=1y(0)=1: differentiate under the integral and integrate by parts. Or complete the square and shift the contour.

Solution

Let g(ξ)=∫e−πx2e−2πixξdxg(\xi)=\int e^{-\pi x^2}e^{-2\pi ix\xi}dx. Then g′(ξ)=∫(−2πix)e−πx2e−2πixξdx=i∫(ddxe−πx2)e−2πixξdx=i⋅(2πiξ)g(ξ)=−2πξg(ξ)g'(\xi)=\int(-2\pi ix)e^{-\pi x^2}e^{-2\pi ix\xi}dx=i\int\big(\tfrac d{dx}e^{-\pi x^2}\big)e^{-2\pi ix\xi}dx=i\cdot(2\pi i\xi)g(\xi)=-2\pi\xi g(\xi), integrating by parts. With g(0)=∫e−πx2=1g(0)=\int e^{-\pi x^2}=1, g(ξ)=e−πξ2g(\xi)=e^{-\pi\xi^2}.

Exercise 6.2The symmetry rules

Prove the rules for translation, modulation, dilation and differentiation stated above.

Hint

Change variables for the first three; integrate by parts for the last.

Exercise 6.3Fourier inversion

Prove inversion for f∈𝒮(ℝ)f\in\mathcal S(\mathbb R).

Hint

Insert e−πε2ξ2e^{-\pi\varepsilon^2\xi^2}: ∫f^(ξ)e−πε2ξ2e2πixξdξ=(f*Gε)(x)\int\hat f(\xi)e^{-\pi\varepsilon^2\xi^2}e^{2\pi ix\xi}d\xi=(f*G_\varepsilon)(x) with GεG_\varepsilon the Gaussian approximate identity, by Fubini and Exercise 6.1. Let ε→0\varepsilon\to0 on both sides.

Solution

By Fubini, ∫f^(ξ)e−πε2ξ2e2πixξdξ=∫f(y)(∫e−πε2ξ2e2πi(x−y)ξdξ)dy\int\hat f(\xi)e^{-\pi\varepsilon^2\xi^2}e^{2\pi ix\xi}d\xi=\int f(y)\Big(\int e^{-\pi\varepsilon^2\xi^2}e^{2\pi i(x-y)\xi}d\xi\Big)dy, and by Exercise 6.1 with a dilation the inner integral is ε−1e−π(x−y)2/ε2\varepsilon^{-1}e^{-\pi(x-y)^2/\varepsilon^2}, a Gaussian approximate identity. So the right side tends to f(x)f(x) (Exercise 3.2), and the left side tends to ∫f^e2πixξ\int\hat fe^{2\pi ix\xi} by dominated convergence, since f^∈𝒮⊆L1\hat f\in\mathcal S\subseteq L^1.

Exercise 6.4Plancherel

Show ∫fg‾=∫f^g^‾\int f\bar g=\int\hat f\bar{\hat g} for f,g∈𝒮f,g\in\mathcal S, and explain how this extends ℱ\mathcal F to a unitary map of L2L^2.

Hint

∫f^h=∫fh^\int\hat f\,h=\int f\,\hat h (Fubini), with h=g^‾h=\bar{\hat g}, whose transform is g‾\bar g by inversion. Then extend by density of 𝒮\mathcal S in L2L^2.

Exercise 6.5Poisson summation

Prove Poisson summation for f∈𝒮(ℝ)f\in\mathcal S(\mathbb R).

Hint

F(x)=∑kf(x+k)F(x)=\sum_kf(x+k) is smooth and 1-periodic. Compute its Fourier coefficients: F^(m)=∫01∑kf(x+k)e−2πimxdx=f^(m)\hat F(m)=\int_0^1\sum_kf(x+k)e^{-2\pi imx}dx=\hat f(m). Its Fourier series converges to it by Exercise 5.5.

Solution

The series for FF and its derivatives converge uniformly because f∈𝒮f\in\mathcal S, so F∈C∞(𝕋)F\in C^\infty(\mathbb T). Unfolding, F^(m)=∑k∫kk+1f(y)e−2πimydy=f^(m)\hat F(m)=\sum_k\int_k^{k+1}f(y)e^{-2\pi imy}dy=\hat f(m). Since F∈C1F\in C^1, its Fourier series converges to it everywhere (Exercise 5.5): F(x)=∑mf^(m)e2πimxF(x)=\sum_m\hat f(m)e^{2\pi imx}.

Exercise 6.6Heisenberg's inequality

Prove ‖f‖22≤4π‖xf‖2‖ξf^‖2\lVert f\rVert_2^2\le4\pi\lVert xf\rVert_2\lVert\xi\hat f\rVert_2 for f∈𝒮(ℝ)f\in\mathcal S(\mathbb R).

Hint

‖f‖22=∫|f|2=−∫xddx|f|2=−2Re∫xf‾f′\lVert f\rVert_2^2=\int|f|^2=-\int x\frac d{dx}|f|^2=-2\mathrm{Re}\int x\bar ff'. Use Cauchy–Schwarz, then Plancherel: ‖f′‖2=2π‖ξf^‖2\lVert f'\rVert_2=2\pi\lVert\xi\hat f\rVert_2.

Exercise 6.7The two-sided exponential

Verify Example 6.3, and deduce Py^(ξ)=e−2πy|ξ|\widehat{P_y}(\xi)=e^{-2\pi y|\xi|} for the Poisson kernel.

Hint

2∫0∞e−2πxcos⁡(2πxξ)dx=2⋅2π(2π)2+(2πξ)22\int_0^\infty e^{-2\pi x}\cos(2\pi x\xi)dx=\frac{2\cdot2\pi}{(2\pi)^2+(2\pi\xi)^2}. Then invert and dilate.

Exercise 6.8Riemann–Lebesgue on ℝ

Show that for f∈L1(ℝ)f\in L^1(\mathbb R), f^\hat f is uniformly continuous, bounded by ‖f‖1\lVert f\rVert_1, and tends to 0 at infinity.

Hint

Continuity: dominated convergence. Decay: true for Gaussians or for 𝒮\mathcal S, and 𝒮\mathcal S is dense in L1L^1.

ReadingGrafakos, Classical, §§2.3–2.4

Tempered distributions

The constant function 1 has no Fourier transform as an integral, and δ\delta is not a function, yet each should be the transform of the other. Distributions make this precise by defining everything through pairing with test functions: a distribution is whatever you can integrate a Schwartz function against.

The idea

Tempered distributions (Grafakos §2.3.3). 𝒮′(ℝn)\mathcal S'(\mathbb R^n) is the space of continuous linear functionals on 𝒮\mathcal S: uu with |⟨u,ϕ⟩|≤C∑|α|,|β|≤Nsup⁡|xα∂βϕ||\langle u,\phi\rangle|\le C\sum_{|\alpha|,|\beta|\le N}\sup|x^\alpha\partial^\beta\phi|. Every ff with ∫|f|(1+|x|)−N<∞\int|f|(1+|x|)^{-N}<\infty defines one by ⟨f,ϕ⟩=∫fϕ\langle f,\phi\rangle=\int f\phi. So do δ\delta (⟨δ,ϕ⟩=ϕ(0)\langle\delta,\phi\rangle=\phi(0)) and its derivatives.
Operations by duality (Grafakos §2.3). Every operation is defined so that it agrees with the usual one on functions: ⟨∂u,ϕ⟩=−⟨u,∂ϕ⟩\langle\partial u,\phi\rangle=-\langle u,\partial\phi\rangle (integration by parts), and ⟨u^,ϕ⟩=⟨u,ϕ^⟩\langle\hat u,\phi\rangle=\langle u,\hat\phi\rangle (from ∫f^ϕ=∫fϕ^\int\hat f\phi=\int f\hat\phi). So every tempered distribution has derivatives of all orders and a Fourier transform, and all the rules of module 6 continue to hold.
The basic transforms. δ^=1\hat\delta=1, 1^=δ\hat1=\delta, e2πiax^=δa\widehat{e^{2\pi iax}}=\delta_a, ∂αδ^=(2πiξ)α\widehat{\partial^\alpha\delta}=(2\pi i\xi)^\alpha. In one dimension, the principal value ⟨p.v.1x,ϕ⟩=limε→0⁡∫|x|>εϕ(x)xdx\langle\mathrm{p.v.}\tfrac1x,\phi\rangle=\lim_{\varepsilon\to0}\int_{|x|>\varepsilon}\frac{\phi(x)}x\,dx is a tempered distribution with ℱ(p.v.1x)=−πisgnξ\mathcal F\big(\mathrm{p.v.}\tfrac1x\big)=-\pi i\,\mathrm{sgn}\,\xi (Exercise 7.4). That is the multiplier of the Hilbert transform in module 9.
Support and homogeneity (Grafakos §2.4). A distribution supported at one point is a finite combination of derivatives of δ\delta there. |x|−a|x|^{-a} on ℝn\mathbb R^n is locally integrable for a<na<n, and its transform is cn,a|ξ|a−nc_{n,a}|\xi|^{a-n}: homogeneity of degree −a-a becomes degree a−na-n. With a=n−2a=n-2 this gives the fundamental solution of the Laplacian, −Δ1(n−2)ωn−1|x|n−2=δ-\Delta\frac{1}{(n-2)\omega_{n-1}|x|^{n-2}}=\delta.

Key terms in plain words

The ideas above without the symbols. Switch to Quiz me to test yourself: you see an explanation and pick the term it describes.

Test function
A very well-behaved function, smooth and rapidly decaying, used as a probe: a generalized function is known by the numbers it produces when integrated against these.
Tempered distribution
A rule assigning a number to each well-behaved probe function, linearly and continuously, with at most polynomial growth allowed at infinity.
Dirac delta
The generalized function that returns the value of the probe at the origin: a unit mass concentrated at a single point.
Distributional derivative
The derivative defined by moving it onto the probe function with a minus sign, as integration by parts would; it exists for every generalized function.
Principal value
A way of integrating 1/x1/x across 0 by cutting out a symmetric gap and letting the gap shrink, so the infinite parts on the two sides cancel.
Fundamental solution
A generalized function whose image under a differential operator is the Dirac delta; convolving with it solves the equation for any right-hand side.

Examples

Example 7.1 (derivatives of kinks and jumps).

ddx|x|=sgnx\frac d{dx}|x|=\mathrm{sgn}\,x and ddxsgnx=2δ\frac d{dx}\mathrm{sgn}\,x=2\delta as distributions (Exercise 7.2). A jump of size JJ differentiates to JδJ\delta, which is why jumps cost 1/m1/m decay in module 5.

Example 7.2 (1 and δ\delta).

⟨1^,ϕ⟩=∫ϕ^=ϕ(0)\langle\hat1,\phi\rangle=\int\hat\phi=\phi(0) by inversion, so 1^=δ\hat1=\delta. Read backwards: the Dirichlet kernels DND_N are partial sums of ∑me2πimx\sum_me^{2\pi imx}, which "equals" the comb ∑kδk\sum_k\delta_k (Example 7.6).

Example 7.3 (principal value).

1/x1/x is not locally integrable, but ⟨p.v.1x,ϕ⟩=∫|x|<1ϕ(x)−ϕ(0)xdx+∫|x|≥1ϕ(x)xdx\langle\mathrm{p.v.}\tfrac1x,\phi\rangle=\int_{|x|<1}\frac{\phi(x)-\phi(0)}x\,dx+\int_{|x|\ge1}\frac{\phi(x)}x\,dx: the first integrand is bounded by sup⁡|ϕ′|\sup|\phi'|, the second by sup⁡|xϕ|/x2\sup|x\phi|/x^2.

Example 7.4 (the sign function).

ℱ(sgn)=1πip.v.1ξ\mathcal F\big(\mathrm{sgn}\big)=\frac1{\pi i}\,\mathrm{p.v.}\frac1\xi, dual to the transform of p.v.1x\mathrm{p.v.}\frac1x. Check with Example 7.1: ℱ(sgn′)=2πiξ⋅1πiξ=2=2δ^\mathcal F\big(\mathrm{sgn}'\big)=2\pi i\xi\cdot\frac1{\pi i\xi}=2=\widehat{2\delta}.

Example 7.5 (the Laplacian in ℝ3\mathbb R^3).

−Δ14π|x|=δ-\Delta\frac1{4\pi|x|}=\delta. So u=14π|x|*fu=\frac1{4\pi|x|}*f solves −Δu=f-\Delta u=f: the Newtonian potential. On the Fourier side, −Δ-\Delta is multiplication by 4π2|ξ|24\pi^2|\xi|^2, and 1/(4π|x|)^=14π2|ξ|2\widehat{1/(4\pi|x|)}=\frac1{4\pi^2|\xi|^2}.

Example 7.6 (the Dirac comb).

Poisson summation (Exercise 6.5) says exactly that the Dirac comb Ш=∑k∈ℤδk\text{Ш}=\sum_{k\in\mathbb Z}\delta_k satisfies Ш^=Ш\hat{\text{Ш}}=\text{Ш}.

Check your understanding

The distributional derivative of the indicator 𝟏[0,1]\mathbf 1_{[0,1]} is:

  1. δ0−δ1\delta_0-\delta_1
  2. 0, since it is constant almost everywhere.That is the pointwise derivative; the distributional one sees the jumps.
  3. δ0+δ1\delta_0+\delta_1The jump at 1 is downwards.
−∫01ϕ′=ϕ(0)−ϕ(1)-\int_0^1\phi'=\phi(0)-\phi(1).

Is exe^{x} a tempered distribution on ℝ\mathbb R?

  1. Yes, every locally integrable function is.Tempered means at most polynomial growth, roughly.
  2. No: ∫exϕ\int e^x\phi diverges for some Schwartz ϕ\phi, e.g. ϕ=e−1+x2\phi=e^{-\sqrt{1+x^2}}.
It is a distribution in the larger space 𝒟′\mathcal D', paired with compactly supported test functions, but it has no Fourier transform in 𝒮′\mathcal S'.

Exercises

Exercise 7.1The principal value is tempered

Show that the limit defining p.v.1x\mathrm{p.v.}\frac1x exists for ϕ∈𝒮\phi\in\mathcal S, equals the formula in Example 7.3, and satisfies |⟨p.v.1x,ϕ⟩|≤2sup⁡|ϕ′|+2sup⁡|xϕ||\langle\mathrm{p.v.}\frac1x,\phi\rangle|\le2\sup|\phi'|+2\sup|x\phi|.

Hint

∫ε<|x|<1ϕ(0)xdx=0\int_{\varepsilon<|x|<1}\frac{\phi(0)}x\,dx=0 by oddness, so you may subtract ϕ(0)\phi(0) inside the inner integral.

Solution

∫|x|>εϕx=∫ε<|x|<1ϕ(x)−ϕ(0)x+∫|x|≥1ϕx\int_{|x|>\varepsilon}\frac\phi x=\int_{\varepsilon<|x|<1}\frac{\phi(x)-\phi(0)}x+\int_{|x|\ge1}\frac\phi x, using oddness. By the mean value theorem the first integrand is bounded by sup⁡|ϕ′|\sup|\phi'|, so it converges as ε→0\varepsilon\to0 to an integral over |x|<1|x|<1 of size at most 2sup⁡|ϕ′|2\sup|\phi'|. The second is at most ∫|x|≥1sup⁡|xϕ|x2=2sup⁡|xϕ|\int_{|x|\ge1}\frac{\sup|x\phi|}{x^2}=2\sup|x\phi|.

Exercise 7.2Derivatives of |x||x| and sgn\mathrm{sgn}

Prove the two identities of Example 7.1.

Hint

−∫|x|ϕ′=−∫0∞xϕ′+∫−∞0xϕ′-\int|x|\phi'=-\int_0^\infty x\phi'+\int_{-\infty}^0x\phi'; integrate by parts on each half-line.

Solution

−∫0∞xϕ′=∫0∞ϕ-\int_0^\infty x\phi'=\int_0^\infty\phi and ∫−∞0xϕ′=−∫−∞0ϕ\int_{-\infty}^0x\phi'=-\int_{-\infty}^0\phi, so ⟨|x|′,ϕ⟩=∫sgnϕ\langle|x|',\phi\rangle=\int\mathrm{sgn}\,\phi. Next, −∫sgnϕ′=−∫0∞ϕ′+∫−∞0ϕ′=ϕ(0)+ϕ(0)-\int\mathrm{sgn}\,\phi'=-\int_0^\infty\phi'+\int_{-\infty}^0\phi'=\phi(0)+\phi(0).

Exercise 7.3δ^=1\hat\delta=1 and 1^=δ\hat1=\delta

Prove both from the definitions.

Hint

⟨δ^,ϕ⟩=ϕ^(0)=∫ϕ\langle\hat\delta,\phi\rangle=\hat\phi(0)=\int\phi; ⟨1^,ϕ⟩=∫ϕ^=ϕ(0)\langle\hat1,\phi\rangle=\int\hat\phi=\phi(0) by inversion at 0.

Exercise 7.4The transform of p.v.1x\mathrm{p.v.}\frac1x

Show ℱ(p.v.1x)=−πisgnξ\mathcal F\big(\mathrm{p.v.}\frac1x\big)=-\pi i\,\mathrm{sgn}\,\xi.

Hint

p.v.1x\mathrm{p.v.}\frac1x is the limit in 𝒮′\mathcal S' of xx2+ε2\frac{x}{x^2+\varepsilon^2}, whose transform is −πisgn(ξ)e−2πε|ξ|-\pi i\,\mathrm{sgn}(\xi)e^{-2\pi\varepsilon|\xi|} (use Exercise 6.7 and xx2+ε2=Re1x−iε\frac x{x^2+\varepsilon^2}=\mathrm{Re}\frac1{x-i\varepsilon}). Or check that x⋅p.v.1x=1x\cdot\mathrm{p.v.}\frac1x=1, transform, and fix the constant by oddness.

Exercise 7.5The Dirac comb is its own transform

Show that Ш=∑kδk\text{Ш}=\sum_k\delta_k is tempered and Ш^=Ш\hat{\text{Ш}}=\text{Ш}.

Hint

|∑kϕ(k)|≤Csup⁡(1+x2)|ϕ||\sum_k\phi(k)|\le C\sup(1+x^2)|\phi|. Then ⟨\Sha^,ϕ⟩=∑kϕ^(k)=∑kϕ(k)\langle\hat\Sha,\phi\rangle=\sum_k\hat\phi(k)=\sum_k\phi(k) is Poisson summation.

Exercise 7.6δ\delta is not a function

Show there is no g∈Lloc1g\in L^1_{\mathrm{loc}} with ∫gϕ=ϕ(0)\int g\phi=\phi(0) for all ϕ∈Cc∞\phi\in C_c^\infty.

Hint

Take ϕε(x)=ϕ(x/ε)\phi_\varepsilon(x)=\phi(x/\varepsilon) with ϕ(0)=1\phi(0)=1, 0≤ϕ≤10\le\phi\le1. Then 1=∫gϕε≤∫|x|<ε|g|→01=\int g\phi_\varepsilon\le\int_{|x|<\varepsilon}|g|\to0.

Exercise 7.7The Newtonian potential

Prove −Δ14π|x|=δ-\Delta\frac1{4\pi|x|}=\delta in 𝒮′(ℝ3)\mathcal S'(\mathbb R^3).

Hint

1|x|\frac1{|x|} is harmonic away from 0. Apply Green's identity to ∫|x|>ε1|x|Δϕ\int_{|x|>\varepsilon}\frac1{|x|}\Delta\phi; the boundary terms on the sphere of radius ε\varepsilon tend to −4πϕ(0)-4\pi\phi(0) and 0.

Solution

For r=|x|>0r=|x|>0, Δ1r=0\Delta\frac1r=0. Green's second identity on {ε<r<R}\{\varepsilon<r<R\}, with RR beyond the support of ϕ\phi, gives ∫r>ε1rΔϕ=∫r=ε(1r∂νϕ−ϕ∂ν1r)dS\int_{r>\varepsilon}\frac1r\Delta\phi=\int_{r=\varepsilon}\big(\frac1r\partial_\nu\phi-\phi\,\partial_\nu\frac1r\big)dS, where ν\nu points towards the origin. The first boundary term is O(ε−1⋅ε2)→0O(\varepsilon^{-1}\cdot\varepsilon^2)\to0. The second is −∫r=εϕ⋅1ε2dS→−4πϕ(0)-\int_{r=\varepsilon}\phi\cdot\frac1{\varepsilon^2}dS\to-4\pi\phi(0). So ⟨Δ14πr,ϕ⟩=−ϕ(0)\langle\Delta\frac1{4\pi r},\phi\rangle=-\phi(0).

InteractiveGrafakos, Classical, §2.1

The maximal function and Lebesgue differentiation

The Hardy–Littlewood maximal function records, at each point, the largest average of |f||f| over balls centred there. It is not integrable, but it is "almost" integrable (weak type (1,1)), and that is exactly enough to prove that averages of an integrable function converge back to it almost everywhere.

The idea

The maximal function (Grafakos §2.1.1). Mf(x)=supr>0⁡1|B(x,r)|∫B(x,r)|f(y)|dyMf(x)=\sup_{r>0}\frac1{|B(x,r)|}\int_{B(x,r)}|f(y)|\,dy. It is sublinear, Mf≥|f|Mf\ge|f| almost everywhere, and ‖Mf‖∞≤‖f‖∞\lVert Mf\rVert_\infty\le\lVert f\rVert_\infty.
The maximal theorem (Grafakos Thm. 2.1.6). Weak type (1,1): |{Mf>λ}|≤3nλ‖f‖1|\{Mf>\lambda\}|\le\frac{3^n}\lambda\lVert f\rVert_1. Combined with the trivial L∞L^\infty bound by Marcinkiewicz interpolation (Grafakos §1.3.1), ‖Mf‖p≤Cn,p‖f‖p\lVert Mf\rVert_p\le C_{n,p}\lVert f\rVert_p for 1<p≤∞1<p\le\infty. The proof of the weak bound is a covering argument: Vitali's lemma picks disjoint balls covering a third of the measure (Exercise 8.1).
Lebesgue differentiation (Grafakos Cor. 2.1.16). For f∈Lloc1f\in L^1_{\mathrm{loc}}, 1|B(x,r)|∫B(x,r)|f−f(x)|→0\frac1{|B(x,r)|}\int_{B(x,r)}|f-f(x)|\to0 for almost every xx. The pattern of proof, "true for a dense class, controlled by a maximal function, hence true a.e.", recurs throughout harmonic analysis. In one dimension it says an indefinite integral F(x)=∫axfF(x)=\int_a^xf has F′=fF'=f almost everywhere: the Lebesgue version of the fundamental theorem of calculus.
Controlling approximate identities. If |Kε||K_\varepsilon| is bounded by a radially decreasing integrable function Φε\Phi_\varepsilon, then supε⁡|Kε*f|≤‖Φ‖1Mf\sup_\varepsilon|K_\varepsilon*f|\le\lVert\Phi\rVert_1\,Mf (Grafakos Thm. 2.1.10). So Gaussian and Poisson means of L1L^1 functions converge almost everywhere, not just in norm.

Explore: Mf and its level sets

Function f
|f||f| , MfMf , the level λ\lambda , and the set {Mf>λ}\{Mf>\lambda\} shaded on the axis.

Key terms in plain words

The ideas above without the symbols. Switch to Quiz me to test yourself: you see an explanation and pick the term it describes.

Hardy–Littlewood maximal function
At each point, the largest average of the function's size over all balls centred at that point.
Weak type (1,1)
The set where an operator's output exceeds a level has measure at most a constant times the input's total integral divided by the level: integrable "up to a logarithm".
Vitali covering lemma
From any finite collection of balls you can choose some that don't overlap, so that tripling their radii covers all the original balls.
Lebesgue point
A point where the averages of the function over shrinking balls converge to its value there; almost every point is one.
Marcinkiewicz interpolation
A way to get a genuine bound on the spaces in between from two weaker bounds at the ends, such as one in the integrable space and one in the bounded space.

Examples

Example 8.1 (one box).

For f=𝟏[−1,1]f=\mathbf 1_{[-1,1]}, Mf=1Mf=1 on [−1,1][-1,1] and Mf(x)=1|x|+1Mf(x)=\frac1{|x|+1} for |x|>1|x|>1 (the best ball just reaches the far end of the box). So Mf∉L1Mf\notin L^1, although |{Mf>λ}|=max⁡(2,2λ−2)≤2λ|\{Mf>\lambda\}|=\max\big(2,\frac2\lambda-2\big)\le\frac2\lambda for 0<λ<10<\lambda<1, well within the bound 3‖f‖1/λ=6/λ3\lVert f\rVert_1/\lambda=6/\lambda.

Example 8.2 (never integrable).

For any f≠0f\ne0 in L1(ℝn)L^1(\mathbb R^n), Mf(x)≥c|x|nMf(x)\ge\frac c{|x|^n} for large |x||x|: a ball centred at xx of radius 2|x|2|x| captures a fixed part of the mass of ff. So the L1L^1 bound genuinely fails, and weak type is the best possible.

Example 8.3 (a tall spike).

f=N𝟏[0,1/N]f=N\mathbf 1_{[0,1/N]} has ‖f‖1=1\lVert f\rVert_1=1 for all NN, while Mf(x)≈12|x|Mf(x)\approx\frac1{2|x|} away from the spike. The level sets do not depend on NN: the weak bound sees only ‖f‖1\lVert f\rVert_1.

Example 8.4 (a singularity).

f(x)=|x|−1/2𝟏|x|<1f(x)=|x|^{-1/2}\mathbf 1_{|x|<1} is in LpL^p for p<2p<2, and MfMf has the same |x|−1/2|x|^{-1/2} singularity at 0, so it is in exactly the same LpL^p spaces near 0, as the LpL^p bound predicts.

Example 8.5 (the fundamental theorem of calculus).

For f∈L1([a,b])f\in L^1([a,b]), F(x+h)−F(x)h=1h∫xx+hf→f(x)\frac{F(x+h)-F(x)}h=\frac1h\int_x^{x+h}f\to f(x) almost everywhere, by Lebesgue differentiation with one-sided intervals. Continuity of ff, as in the Riemann version, is not needed.

Example 8.6 (heat and Poisson means).

The Gaussian and Poisson kernels are radially decreasing with integral 1, so supt⁡|Gt*f|≤Mf\sup_t|G_t*f|\le Mf and supy⁡|Py*f|≤Mf\sup_y|P_y*f|\le Mf. Hence Py*f→fP_y*f\to f almost everywhere for f∈L1f\in L^1: Fatou's theorem on boundary values of harmonic functions.

Check your understanding

For nonzero f∈L1(ℝ)f\in L^1(\mathbb R), is Mf∈L1(ℝ)Mf\in L^1(\mathbb R)?

  1. Yes, by the maximal theorem.The theorem gives weak type (1,1), not an L1L^1 bound.
  2. Never: Mf≥c/|x|Mf\ge c/|x| for large |x||x|.
Example 8.2. The failure is at infinity; near any point MfMf can be integrable or not depending on ff (Example 8.4).

In the proof of Lebesgue differentiation, what is the maximal function used for?

  1. To show the set where averages misbehave is small, after approximating ff by a continuous function, for which the result is clear.
  2. To compute the limit of the averages.It bounds them; the limit comes from the continuous approximation.
Exercise 8.3: "dense class plus maximal inequality" is the template.

Exercises

Exercise 8.1Vitali's covering lemma

Show that from finitely many balls B1,…,BkB_1,\dots,B_k in ℝn\mathbb R^n one can choose pairwise disjoint ones Bi1,…,BilB_{i_1},\dots,B_{i_l} with ⋃jBj⊆⋃m3Bim\bigcup_jB_j\subseteq\bigcup_m3B_{i_m}, so |⋃Bj|≤3n∑m|Bim||\bigcup B_j|\le3^n\sum_m|B_{i_m}|.

Hint

Greedy: take the largest ball, discard every ball meeting it, repeat. A discarded ball meets a chosen one at least as large, so it lies in its triple.

Solution

Order the balls by decreasing radius. Choose Bi1B_{i_1} as the largest; remove all balls meeting it; choose the largest remaining; and so on. The chosen balls are disjoint. If BjB_j was removed because it meets a chosen BB with radius r≥rad(Bj)r\ge\mathrm{rad}(B_j), then every point of BjB_j is within r+2rad(Bj)≤3rr+2\,\mathrm{rad}(B_j)\le3r of the centre of BB, so Bj⊆3BB_j\subseteq3B. Finally |3B|=3n|B||3B|=3^n|B|.

Exercise 8.2Weak type (1,1)

Prove |{Mf>λ}|≤3nλ‖f‖1|\{Mf>\lambda\}|\le\frac{3^n}\lambda\lVert f\rVert_1.

Hint

Each xx in E={Mf>λ}E=\{Mf>\lambda\} is the centre of a ball BxB_x with ∫Bx|f|>λ|Bx|\int_{B_x}|f|>\lambda|B_x|. Cover a compact K⊆EK\subseteq E by finitely many, apply Vitali, and use inner regularity.

Solution

For compact K⊆EK\subseteq E, finitely many BxB_x cover KK. Vitali gives disjoint ones BmB_m with |K|≤3n∑|Bm|<3nλ∑∫Bm|f|≤3nλ‖f‖1|K|\le3^n\sum|B_m|<\frac{3^n}\lambda\sum\int_{B_m}|f|\le\frac{3^n}\lambda\lVert f\rVert_1. Take the supremum over KK (Lebesgue measure is inner regular, and EE is open since MfMf is lower semicontinuous).

Exercise 8.3Lebesgue differentiation

Deduce the differentiation theorem from Exercise 8.2.

Hint

Let Ωf(x)=lim supr→0⁡1|B|∫B|f−f(x)|\Omega f(x)=\limsup_{r\to0}\frac1{|B|}\int_B|f-f(x)|. For continuous gg, Ωg=0\Omega g=0. Show Ωf≤M(f−g)+|f−g|\Omega f\le M(f-g)+|f-g|, so |{Ωf>2λ}|≤3n+1λ‖f−g‖1|\{\Omega f>2\lambda\}|\le\frac{3^n+1}\lambda\lVert f-g\rVert_1, which is as small as you like.

Exercise 8.4Compute Example 8.1

Verify the formula for M𝟏[−1,1]M\mathbf 1_{[-1,1]} and the measure of its level sets.

Hint

For x>1x>1 the average over [x−r,x+r][x-r,x+r] is r−x+12r\frac{r-x+1}{2r} for x−1≤r≤x+1x-1\le r\le x+1, which increases in rr, then 1r\frac1r for larger rr. Maximum at r=x+1r=x+1.

Exercise 8.5Not integrable at infinity

Prove Example 8.2.

Hint

Pick RR with ∫|y|<R|f|=a>0\int_{|y|<R}|f|=a>0. For |x|>R|x|>R, the ball B(x,2|x|)B(x,2|x|) contains B(0,R)B(0,R), so Mf(x)≥a|B(0,1)|(2|x|)nMf(x)\ge\frac a{|B(0,1)|(2|x|)^n}.

Exercise 8.6From weak (1,1) to LpL^p

Prove ‖Mf‖p≤Cp‖f‖p\lVert Mf\rVert_p\le C_p\lVert f\rVert_p for 1<p<∞1<p<\infty directly, without quoting Marcinkiewicz.

Hint

Split f=f𝟏|f|>λ/2+f𝟏|f|≤λ/2f=f\mathbf 1_{|f|>\lambda/2}+f\mathbf 1_{|f|\le\lambda/2}. The second part has M≤λ/2M\le\lambda/2, so {Mf>λ}⊆{M(f𝟏|f|>λ/2)>λ/2}\{Mf>\lambda\}\subseteq\{M(f\mathbf 1_{|f|>\lambda/2})>\lambda/2\}. Apply the weak bound and ‖Mf‖pp=p∫0∞λp−1|{Mf>λ}|dλ\lVert Mf\rVert_p^p=p\int_0^\infty\lambda^{p-1}|\{Mf>\lambda\}|d\lambda, then Fubini.

Exercise 8.7Radial kernels are controlled by MM

Show that if Φ≥0\Phi\ge0 is radial, decreasing and integrable on ℝn\mathbb R^n, then |Φ*f|≤‖Φ‖1Mf|\Phi*f|\le\lVert\Phi\rVert_1Mf.

Hint

Write Φ\Phi as a superposition of normalized indicators of balls (the layer-cake formula). Each indicator convolved with |f||f| is at most MfMf times its integral.

InteractiveGrafakos, Classical, §§4.1, 5.1, 5.3

The Hilbert transform and singular integrals

Convolution with 1πx\frac1{\pi x} is the simplest operator that is not an average: its kernel is not integrable, and it is defined only as a principal value. It is the multiplier −isgnξ-i\,\mathrm{sgn}\,\xi, a 90° phase shift of every frequency. Whether it is bounded on LpL^p is equivalent to whether Fourier series converge in LpL^p.

The idea

Definition (Grafakos §5.1). Hf(x)=1πp.v.∫f(y)x−ydyHf(x)=\frac1\pi\,\mathrm{p.v.}\int\frac{f(y)}{x-y}dy, that is, Hf=1πp.v.1x*fHf=\frac1\pi\mathrm{p.v.}\frac1x*f. By module 7, Hf^(ξ)=−isgn(ξ)f^(ξ)\widehat{Hf}(\xi)=-i\,\mathrm{sgn}(\xi)\hat f(\xi). So HH is an isometry of L2L^2 with H2=−IH^2=-I, and Hcos⁡=sin⁡H\cos=\sin, Hsin⁡=−cos⁡H\sin=-\cos.
Riesz's theorem (Grafakos §§4.1.2, 5.1). HH is bounded on LpL^p for 1<p<∞1<p<\infty, but not on L1L^1 or L∞L^\infty: H𝟏[a,b]=1πlog⁡|x−ax−b|H\mathbf 1_{[a,b]}=\frac1\pi\log\big|\frac{x-a}{x-b}\big| is unbounded and decays only like 1/x1/x. At p=1p=1 it is weak type (1,1), proved by the Calderón–Zygmund decomposition.
Complex analysis. f+iHff+iHf is the boundary value of the holomorphic function F=(Py+iQy)*fF=(P_y+iQ_y)*f on the upper half-plane, where Qy(x)=1πxx2+y2Q_y(x)=\frac1\pi\frac x{x^2+y^2} is the conjugate Poisson kernel. On the circle the same operator is the conjugate function f~=∑m−isgn(m)f^(m)e2πimx\tilde f=\sum_m-i\,\mathrm{sgn}(m)\hat f(m)e^{2\pi imx}, and the partial sums SNS_N can be written using it. So SNf→fS_Nf\to f in LpL^p for all f∈Lpf\in L^p exactly when the conjugate function is bounded on LpL^p: true for 1<p<∞1<p<\infty, false for p=1,∞p=1,\infty (Grafakos §4.1.1).
Calderón–Zygmund decomposition (Grafakos §5.3). Given f∈L1f\in L^1 and α>0\alpha>0, write f=g+bf=g+b with |g|≤2nα|g|\le2^n\alpha and b=∑bjb=\sum b_j supported on disjoint dyadic cubes QjQ_j, each with ∫bj=0\int b_j=0 and ∑|Qj|≤1α‖f‖1\sum|Q_j|\le\frac1\alpha\lVert f\rVert_1. Bounded on gg by L2L^2, and on each bjb_j because of the cancellation, an operator like HH becomes weak type (1,1). The Riesz transforms RjR_j, with multipliers −iξj/|ξ|-i\xi_j/|\xi|, are the nn-dimensional versions.

Explore: the Hilbert transform of a few functions

Function f
ff and HfHf , computed by a principal-value quadrature (and compared with the exact formula where there is one).

Key terms in plain words

The ideas above without the symbols. Switch to Quiz me to test yourself: you see an explanation and pick the term it describes.

Hilbert transform
The operator that shifts the phase of every positive frequency by a quarter turn one way and every negative frequency a quarter turn the other way, turning cosines into sines.
Singular integral
An operator given by integrating against a kernel that blows up like one over the distance to the diagonal, too singular to be an average, but tamed by cancellation.
Conjugate function
On the circle, the function whose sum with ii times the original is the boundary value of an analytic function; it is the periodic Hilbert transform.
Calderón–Zygmund decomposition
Splitting a function at a height into a bounded part and a bad part made of pieces with zero average, living on disjoint cubes of small total size.
Fourier multiplier
An operator that multiplies each frequency component by its own number; the list of those numbers, as a function of frequency, is the symbol.
Riesz transforms
The higher-dimensional cousins of the phase-shift operator, one per coordinate, with symbols −iξj/|ξ|-i\xi_j/|\xi|.

Examples

Example 9.1 (a phase shift).

cos⁡(2πx)=12(e2πix+e−2πix)\cos(2\pi x)=\frac12(e^{2\pi ix}+e^{-2\pi ix}), and multiplying the two frequencies by −i-i and +i+i gives sin⁡(2πx)\sin(2\pi x). (These are tempered distributions, not LpL^p functions, but the multiplier formula applies.)

Example 9.2 (a box).

H𝟏[a,b](x)=1πlog⁡|x−ax−b|H\mathbf 1_{[a,b]}(x)=\frac1\pi\log\big|\frac{x-a}{x-b}\big|: logarithmic spikes of opposite signs at the two ends, and a b−aπx\frac{b-a}{\pi x} tail. So HH does not map L∞L^\infty to L∞L^\infty or L1L^1 to L1L^1, although H𝟏[a,b]∈LpH\mathbf 1_{[a,b]}\in L^p for every 1<p<∞1<p<\infty.

Example 9.3 (from complex analysis).

F(z)=iz+iF(z)=\frac i{z+i} is holomorphic in the upper half-plane, and on the real line F(x)=11+x2+ix1+x2F(x)=\frac1{1+x^2}+i\frac x{1+x^2}. So H11+x2=x1+x2H\frac1{1+x^2}=\frac x{1+x^2}: the Poisson kernel's conjugate is the conjugate Poisson kernel (at height 1).

Example 9.4 (L2L^2 norms).

For the box of length ℓ\ell, ‖𝟏‖22=ℓ\lVert\mathbf 1\rVert_2^2=\ell, and the figure's numerical ‖Hf‖22\lVert Hf\rVert_2^2 agrees, as the multiplier |−isgnξ|=1|-i\,\mathrm{sgn}\,\xi|=1 predicts.

Example 9.5 (Riesz transforms).

∑jRj2=−I\sum_jR_j^2=-I, since ∑j(−iξj/|ξ|)2=−1\sum_j(-i\xi_j/|\xi|)^2=-1. And ∂j∂ku=−RjRkΔu\partial_j\partial_ku=-R_jR_k\Delta u, so LpL^p bounds for Riesz transforms give LpL^p bounds for all second derivatives from the Laplacian alone, the starting point of elliptic regularity.

Example 9.6 (Fourier series in LpL^p).

On the circle let PP keep the frequencies m≥0m\ge0; then P=12(I+iH~)+12E0P=\frac12(I+i\tilde H)+\frac12E_0, where H~\tilde H is the conjugate function and E0f=f^(0)E_0f=\hat f(0). With Maf=e2πiaxfM_af=e^{2\pi iax}f, which shifts frequencies by aa and preserves LpL^p norms, SN=M−NPMN−MN+1PM−N−1S_N=M_{-N}PM_N-M_{N+1}PM_{-N-1}. So a bounded conjugate function gives supN⁡‖SN‖Lp→Lp<∞\sup_N\lVert S_N\rVert_{L^p\to L^p}<\infty, which is LpL^p convergence of every Fourier series.

Check your understanding

On which spaces is the Hilbert transform bounded?

  1. LpL^p for 1<p<∞1<p<\infty.
  2. LpL^p for 1≤p≤∞1\le p\le\infty.Example 9.2 rules out p=1p=1 and p=∞p=\infty.
  3. Only L2L^2.L2L^2 is the easy case; Riesz proved the rest.
At the endpoints it is weak type (1,1), and it maps L∞L^\infty into BMO (module 10).

What is H(Hf)H(Hf)?

  1. −f-f
  2. ff(−isgnξ)2=−1(-i\,\mathrm{sgn}\,\xi)^2=-1.
  3. 0HH is an isometry of L2L^2.
Two quarter-turn phase shifts make a half turn.

Exercises

Exercise 9.1The Hilbert transform of a box

Compute H𝟏[a,b]H\mathbf 1_{[a,b]} and show it is not bounded and not integrable.

Hint

For x∉[a,b]x\notin[a,b] the integral is ordinary: 1π∫abdyx−y\frac1\pi\int_a^b\frac{dy}{x-y}. For x∈(a,b)x\in(a,b), cut out (x−ε,x+ε)(x-\varepsilon,x+\varepsilon) and note the two logarithms of ε\varepsilon cancel.

Solution

1π∫abdyx−y=1π[−log⁡|x−y|]ab=1πlog⁡|x−a||x−b|\frac1\pi\int_a^b\frac{dy}{x-y}=\frac1\pi\big[-\log|x-y|\big]_a^b=\frac1\pi\log\frac{|x-a|}{|x-b|}, and inside (a,b)(a,b) the principal value gives the same formula. It tends to ±∞\pm\infty at aa and bb, and for large |x||x| it equals 1πlog⁡(1+b−ax−b)≈b−aπx\frac1\pi\log\big(1+\frac{b-a}{x-b}\big)\approx\frac{b-a}{\pi x}, which is not integrable.

Exercise 9.2The multiplier

Using Exercise 7.4, show Hf^=−isgn(ξ)f^\widehat{Hf}=-i\,\mathrm{sgn}(\xi)\hat f, and deduce that HH is an isometry of L2L^2 with H2=−IH^2=-I.

Hint

u*f^=u^f^\widehat{u*f}=\hat u\hat f for u∈𝒮′u\in\mathcal S' and f∈𝒮f\in\mathcal S; 1π⋅(−πisgnξ)\frac1\pi\cdot(-\pi i\,\mathrm{sgn}\,\xi). Then Plancherel.

Exercise 9.3Phase shifts

Show Hcos⁡(2πkx)=sin⁡(2πkx)H\cos(2\pi kx)=\sin(2\pi kx) for k>0k>0, and check Example 9.3 by evaluating the principal value integral with residues.

Hint

For Example 9.3, close the contour in the lower or upper half-plane and add half the residue at the real pole (the principal value contributes πi\pi i times the residue there).

Exercise 9.4Symmetries of HH

Show HH commutes with translations and positive dilations and anticommutes with the reflection f(x)↦f(−x)f(x)\mapsto f(-x). Conversely, show that a bounded operator on L2L^2 with these three properties is cHcH for a constant cc.

Hint

Commuting with translations makes it a multiplier m(ξ)m(\xi) (the Fourier transform diagonalizes translations, exactly as in the finite case). Dilations force mm to be constant on (0,∞)(0,\infty) and on (−∞,0)(-\infty,0); reflection forces m(−ξ)=−m(ξ)m(-\xi)=-m(\xi).

Exercise 9.5Riesz's theorem for p=4p=4

Prove ‖Hf‖4≤C‖f‖4\lVert Hf\rVert_4\le C\lVert f\rVert_4 for real f∈𝒮f\in\mathcal S, using the identity (Hf)2=f2+2H(fHf)(Hf)^2=f^2+2H(fHf).

Hint

The identity holds because (f+iHf)2(f+iHf)^2 is again a boundary value of a holomorphic function. Then ‖Hf‖42=‖(Hf)2‖2≤‖f‖42+2‖fHf‖2≤‖f‖42+2‖f‖4‖Hf‖4\lVert Hf\rVert_4^2=\lVert(Hf)^2\rVert_2\le\lVert f\rVert_4^2+2\lVert fHf\rVert_2\le\lVert f\rVert_4^2+2\lVert f\rVert_4\lVert Hf\rVert_4, a quadratic inequality in ‖Hf‖4/‖f‖4\lVert Hf\rVert_4/\lVert f\rVert_4. Iterating gives all p=2kp=2^k, and interpolation and duality give the rest.

Exercise 9.6The Calderón–Zygmund decomposition

Construct the decomposition in one dimension: given f∈L1(ℝ)f\in L^1(\mathbb R) and α>0\alpha>0, select the maximal dyadic intervals QQ with 1|Q|∫Q|f|>α\frac1{|Q|}\int_Q|f|>\alpha, and show α<1|Q|∫Q|f|≤2α\alpha<\frac1{|Q|}\int_Q|f|\le2\alpha, ∑|Q|≤1α‖f‖1\sum|Q|\le\frac1\alpha\lVert f\rVert_1, and |f|≤α|f|\le\alpha a.e. off ⋃Q\bigcup Q.

Hint

Maximality: the dyadic parent Q~\tilde Q of a chosen QQ has average ≤α\le\alpha, and |Q~|=2|Q||\tilde Q|=2|Q|. Off the union, every dyadic interval containing xx has average ≤α\le\alpha; use Lebesgue differentiation with dyadic intervals. Set g=fg=f off ⋃Q\bigcup Q and g=g= the average of ff on each QQ.

Exercise 9.7Riesz transforms and second derivatives

Verify ∑jRj2=−I\sum_jR_j^2=-I and ∂j∂ku=−RjRkΔu\partial_j\partial_ku=-R_jR_k\Delta u for u∈𝒮(ℝn)u\in\mathcal S(\mathbb R^n).

Hint

On the Fourier side: ∂j∂k\partial_j\partial_k is −4π2ξjξk-4\pi^2\xi_j\xi_k, Δ\Delta is −4π2|ξ|2-4\pi^2|\xi|^2, and RjRkR_jR_k is −ξjξk/|ξ|2-\xi_j\xi_k/|\xi|^2.

InteractiveGrafakos, Classical, §6.1; Grafakos, Modern, §§1.3, 3.1–3.2

Littlewood–Paley, Sobolev spaces and BMO

Split a function into frequency bands an octave wide. In L2L^2 the bands are orthogonal; the Littlewood–Paley theorem says that, in a square-function sense, they remain almost orthogonal in every LpL^p. Measuring the bands with weights 2js2^{js} gives Sobolev smoothness, and the right endpoint space at p=∞p=\infty is BMO, the functions of bounded mean oscillation. This is where Modern Fourier Analysis begins.

The idea

Littlewood–Paley theorem (Grafakos, Classical §6.1). Choose ψ\psi with ψ^\hat\psi supported in 12≤|ξ|≤2\frac12\le|\xi|\le2 and ∑j∈ℤψ^(2−jξ)=1\sum_{j\in\mathbb Z}\hat\psi(2^{-j}\xi)=1 for ξ≠0\xi\ne0, and put Δjf^=ψ^(2−jξ)f^\widehat{\Delta_jf}=\hat\psi(2^{-j}\xi)\hat f. Then for 1<p<∞1<p<\infty,
‖f‖p≈‖(∑j|Δjf|2)1/2‖p.\lVert f\rVert_p\approx\Big\lVert\Big(\sum_j|\Delta_jf|^2\Big)^{1/2}\Big\rVert_p.
At p=2p=2 it is Plancherel plus the finite overlap of the bands; for other pp it is a vector-valued singular integral estimate, the machinery of module 9.
Sobolev spaces (Grafakos, Modern §1.3). Hs(ℝn)={f:(1+|ξ|2)s/2f^∈L2}H^s(\mathbb R^n)=\{f:(1+|\xi|^2)^{s/2}\hat f\in L^2\}: "ss derivatives in L2L^2", for any real ss. Equivalently ∑j22js‖Δjf‖22<∞\sum_j2^{2js}\lVert\Delta_jf\rVert_2^2<\infty (plus a low-frequency part). The LpL^p versions use the square function. Sobolev embedding: Hs⊆C0H^s\subseteq C_0 when s>n/2s>n/2 (Exercise 10.1), and it fails at s=n/2s=n/2.
BMO (Grafakos, Modern §§3.1–3.2). ‖f‖BMO=supQ⁡1|Q|∫Q|f−fQ|\lVert f\rVert_{BMO}=\sup_Q\frac1{|Q|}\int_Q|f-f_Q|, with fQf_Q the average over the cube QQ. L∞⊊BMOL^\infty\subsetneq BMO: log⁡|x|\log|x| is in BMO. Singular integrals such as HH map L∞L^\infty to BMO, BMO is the dual of the Hardy space H1H^1 (Fefferman), and the John–Nirenberg inequality says BMO functions are exponentially integrable on each cube.

Explore: dyadic frequency bands and Sobolev norms

Function on the circle
Energy ∑2j−1≤|m|<2j|f^(m)|2\sum_{2^{j-1}\le|m|<2^j}|\hat f(m)|^2 of each band jj , and the same weighted by 22js2^{2js}
Partial sums ∑j≤J22js(band energy)\sum_{j\le J}2^{2js}(\text{band energy}) against JJ: bounded exactly when f∈Hsf\in H^s

Key terms in plain words

The ideas above without the symbols. Switch to Quiz me to test yourself: you see an explanation and pick the term it describes.

Littlewood–Paley decomposition
Splitting a function into pieces whose frequencies lie in bands an octave wide, so that each piece oscillates at roughly one scale.
Square function
At each point, the square root of the sum of the squared sizes of all the frequency-band pieces there; its size in any integrable power measures the size of the function.
Sobolev space
The functions with a given number of derivatives, possibly fractional, in a given integrable power; measured by weighting each frequency by a power of its size.
Sobolev embedding
Enough derivatives in an average sense, compared with the dimension, force a function to be continuous and bounded.
Bounded mean oscillation
Describes a function whose average distance from its own mean, over every cube, is bounded, even if the function itself is not, like the logarithm.
John–Nirenberg inequality
The statement that functions of bounded average oscillation deviate from their means on large sets only with exponentially small probability.

Examples

Example 10.1 (a pure frequency).

e2πimxe^{2\pi imx} lives in the single band jj with 2j−1≤|m|<2j2^{j-1}\le|m|<2^j, so Δj\Delta_j of it is itself or 0, and the square function is just |f||f|.

Example 10.2 (a jump costs half a derivative).

The square wave has |f^(m)|≈1m|\hat f(m)|\approx\frac1m, so band jj has energy about 2−j2^{-j}, and ∑j22js2−j<∞\sum_j2^{2js}2^{-j}<\infty iff s<12s<\frac12. Functions with jumps are in HsH^s exactly for s<12s<\frac12. The triangle wave reaches s<32s<\frac32.

Example 10.3 (Sobolev embedding in one dimension).

For f∈H1(ℝ)f\in H^1(\mathbb R), |f(x)−f(y)|=|∫yxf′|≤|x−y|1/2‖f′‖2|f(x)-f(y)|=\big|\int_y^xf'\big|\le|x-y|^{1/2}\lVert f'\rVert_2 by Cauchy–Schwarz: H1H^1 functions are Hölder continuous of order 12\frac12.

Example 10.4 (log⁡|x|\log|x| is in BMO).

Its oscillation over [a,b][a,b] is unchanged by rescaling x↦λxx\mapsto\lambda x (which adds a constant) so it suffices to check intervals of length 1. Those near 0 have bounded oscillation because log⁡|x|\log|x| is integrable; those far away because it is nearly constant there. On the circle, log⁡|sin⁡πx|\log|\sin\pi x| in the figure has the same singularity, |f^(m)|≈12|m||\hat f(m)|\approx\frac1{2|m|}, and is in HsH^s for s<12s<\frac12 but not bounded.

Example 10.5 (HH maps L∞L^\infty into BMO).

H𝟏[a,b]=1πlog⁡|x−ax−b|H\mathbf 1_{[a,b]}=\frac1\pi\log\big|\frac{x-a}{x-b}\big| from Example 9.2 is a difference of two translated copies of 1πlog⁡|x|\frac1\pi\log|x|: unbounded, but in BMO.

Example 10.6 (embedding fails at the endpoint).

log⁡log⁡e|x|\log\log\frac e{|x|}, cut off near |x|=1|x|=1, is in H1(ℝ2)H^1(\mathbb R^2) but unbounded at 0. With n=2n=2 and s=1=n2s=1=\frac n2, Hs⊈L∞H^s\not\subseteq L^\infty; such functions do lie in BMO.

Check your understanding

For which ss is the square wave in Hs(𝕋)H^s(\mathbb T)?

  1. s<12s<\frac12
  2. s≤1s\le1A jump can't be in H1H^1: H1H^1 functions are continuous in one dimension.
  3. s=0s=0 onlyMove the slider: the partial sums stay bounded for small positive ss.
∑mm2sm−2\sum_m m^{2s}m^{-2} converges iff 2s−2<−12s-2<-1.

Which of these is in BMO but not in L∞L^\infty?

  1. log⁡|x|\log|x|
  2. |x|−1/2|x|^{-1/2}Its oscillation over [0,r][0,r] grows like r−1/2r^{-1/2} (Exercise 10.6).
  3. sgnx\mathrm{sgn}\,xThat is bounded.
Power singularities are too strong for BMO; logarithmic ones are just right, by John–Nirenberg.

Exercises

Exercise 10.1Sobolev embedding

Show Hs(ℝn)⊆C0(ℝn)H^s(\mathbb R^n)\subseteq C_0(\mathbb R^n) for s>n/2s>n/2, with ‖f‖∞≤Cn,s‖f‖Hs\lVert f\rVert_\infty\le C_{n,s}\lVert f\rVert_{H^s}.

Hint

‖f‖∞≤‖f^‖1=∫(1+|ξ|2)−s/2⋅(1+|ξ|2)s/2|f^|\lVert f\rVert_\infty\le\lVert\hat f\rVert_1=\int(1+|\xi|^2)^{-s/2}\cdot(1+|\xi|^2)^{s/2}|\hat f|, then Cauchy–Schwarz.

Solution

By Cauchy–Schwarz, ‖f^‖1≤(∫(1+|ξ|2)−sdξ)1/2‖f‖Hs\lVert\hat f\rVert_1\le\big(\int(1+|\xi|^2)^{-s}d\xi\big)^{1/2}\lVert f\rVert_{H^s}, and the first integral is finite iff 2s>n2s>n (polar coordinates). So f^∈L1\hat f\in L^1, and ff is the inverse transform of an integrable function: continuous and vanishing at infinity by Riemann–Lebesgue (Exercise 6.8).

Exercise 10.2Littlewood–Paley in L2L^2

Show 12‖f‖22≤∑j‖Δjf‖22≤‖f‖22\frac12\lVert f\rVert_2^2\le\sum_j\lVert\Delta_jf\rVert_2^2\le\lVert f\rVert_2^2 when 0≤ψ^≤10\le\hat\psi\le1.

Hint

By Plancherel, ∑j‖Δjf‖22=∫∑jψ^(2−jξ)2|f^|2\sum_j\lVert\Delta_jf\rVert_2^2=\int\sum_j\hat\psi(2^{-j}\xi)^2|\hat f|^2. Each ξ≠0\xi\ne0 lies in at most two bands, and ∑ψ^=1\sum\hat\psi=1 with 0≤ψ^≤10\le\hat\psi\le1 gives 12≤∑ψ^2≤1\frac12\le\sum\hat\psi^2\le1.

Exercise 10.3Jumps and corners

Verify Example 10.2, and show that the triangle wave is in Hs(𝕋)H^s(\mathbb T) exactly for s<32s<\frac32. Check both with the figure.

Hint

∑mm2s−4\sum_m m^{2s-4} converges iff 2s−4<−12s-4<-1.

Exercise 10.4log⁡|x|∈BMO\log|x|\in BMO

Prove Example 10.4 on ℝ\mathbb R.

Hint

By scaling invariance, take |I|=1|I|=1. If I⊆[−2,2]I\subseteq[-2,2], compare with the constant 0: ∫I|log⁡|x||≤∫−22|log⁡|x||<∞\int_I|\log|x||\le\int_{-2}^2|\log|x||<\infty. If II is far from 0, log⁡|x|\log|x| varies by at most log⁡|c|+1|c|−1\log\frac{|c|+1}{|c|-1} on it, where cc is the centre. Use 1|I|∫I|f−fI|≤2|I|∫I|f−a|\frac1{|I|}\int_I|f-f_I|\le\frac2{|I|}\int_I|f-a| for any constant aa.

Solution

First, for any constant aa, 1|I|∫I|f−fI|≤2|I|∫I|f−a|\frac1{|I|}\int_I|f-f_I|\le\frac2{|I|}\int_I|f-a| (since |fI−a|≤1|I|∫|f−a||f_I-a|\le\frac1{|I|}\int|f-a|). Replacing II by λI\lambda I changes log⁡|x|\log|x| by the constant log⁡λ\log\lambda, which doesn't change the oscillation, so assume |I|=1|I|=1. If I⊆[−2,2]I\subseteq[-2,2], take a=0a=0: the oscillation is at most 2∫−22|log⁡|x||dx2\int_{-2}^2|\log|x||\,dx. Otherwise |x|≥1|x|\ge1 on II with |x||x| between |c|−12|c|-\frac12 and |c|+12|c|+\frac12; take a=log⁡|c|a=\log|c|, and |log⁡|x|−log⁡|c||≤log⁡3|\log|x|-\log|c||\le\log 3. Both bounds are independent of II.

Exercise 10.5BMO is a seminorm space

Show ‖⋅‖BMO\lVert\cdot\rVert_{BMO} is a seminorm vanishing exactly on constants, and ‖f‖BMO≤2‖f‖∞\lVert f\rVert_{BMO}\le2\lVert f\rVert_\infty.

Hint

If every oscillation is 0, ff equals its average on every cube, hence on nested cubes, hence is constant a.e. For the bound, use the inequality from Exercise 10.4 with a=0a=0.

Exercise 10.6Power singularities are not BMO

Show |x|−1/2𝟏|x|<1∉BMO(ℝ)|x|^{-1/2}\mathbf 1_{|x|<1}\notin BMO(\mathbb R).

Hint

On I=[0,r]I=[0,r], fI=2r−1/2f_I=2r^{-1/2}, and f>2fIf>2f_I on [0,r/16][0,r/16], so the oscillation is at least 1r∫0r/16(f−fI)≳r−1/2→∞\frac1r\int_0^{r/16}(f-f_I)\gtrsim r^{-1/2}\to\infty.

Exercise 10.7Consequences of John–Nirenberg

John–Nirenberg says |{x∈Q:|f−fQ|>λ}|≤C|Q|e−cλ/‖f‖BMO|\{x\in Q:|f-f_Q|>\lambda\}|\le C|Q|e^{-c\lambda/\lVert f\rVert_{BMO}}. Deduce that BMO functions are in LlocpL^p_{\mathrm{loc}} for all p<∞p<\infty, with (1|Q|∫Q|f−fQ|p)1/p≤Cp‖f‖BMO\big(\frac1{|Q|}\int_Q|f-f_Q|^p\big)^{1/p}\le C_p\lVert f\rVert_{BMO}.

Hint

Layer cake: ∫Q|f−fQ|p=p∫0∞λp−1|{|f−fQ|>λ}|dλ\int_Q|f-f_Q|^p=p\int_0^\infty\lambda^{p-1}|\{|f-f_Q|>\lambda\}|d\lambda, and the exponential tail makes it finite.