Mathematics · Book 4 · Bachelor Year 2

University Mathematics — Year 2

University Mathematics — Year 2 · Bachelor Year 2

14Fourier Series

Can every periodic signal be rebuilt from pure sines and cosines? Fourier’s audacious “yes” created a century of analysis. This chapter proves the two pillars within MP* reach: Dirichlet’s theorem (pointwise reconstruction for piecewise C1C^1 functions, via the Dirichlet kernel) and Parseval’s identity (the energy of a signal is the sum of the energies of its harmonics), and harvests the classical numerical series — Basel’s 1/n2=π2/6\sum 1/n^2 = \pi^2/6 first among them.

Throughout, functions are 2π2\pi-periodic, piecewise continuous, complex-valued; C\mathcal{C} denotes the continuous ones.

14.1 Fourier coefficients

Definition 14.1

The Fourier coefficients of ff are

cn(f)=12πππf(t)eint ⁣dt(nZ),c_n(f) = \frac{1}{2\pi}\int_{-\pi}^{\pi} f(t)\,\eu^{-\iu n t}\,\dd t \qquad (n \in \Z),

and the real-form coefficients an=cn+cna_n = c_n + c_{-n}, bn=i(cncn)b_n = \iu(c_n - c_{-n}), so that the Fourier partial sums are

SN(f)(t)=n=NNcn(f)eint=a02+n=1N(ancosnt+bnsinnt).S_N(f)(t) = \sum_{n=-N}^{N} c_n(f)\,\eu^{\iu nt} = \frac{a_0}{2} + \sum_{n=1}^{N}\bigl(a_n\cos nt + b_n \sin nt\bigr).

On C\mathcal{C}, define the Hermitian inner product f,g=12πππfg\langle f, g\rangle = \frac{1}{2\pi}\int_{-\pi}^{\pi}\conj f\,g: the exponentials en(t)=einte_n(t) = \eu^{\iu nt} are orthonormal (em,en=δmn\langle e_m, e_n\rangle = \delta_{mn}, direct computation), and cn(f)=en,fc_n(f) = \langle e_n, f\rangle: Fourier analysis is Hermitian geometry (Chapter 13) in infinite dimension.

Proposition 14.2 (Bessel’s inequality)

SN(f)S_N(f) is the orthogonal projection of ff onto the space TN\mathcal{T}_N of trigonometric polynomials of degree N\leq N, and

n=NNcn(f)2f22=12πππf2:\sum_{n=-N}^{N} \abs{c_n(f)}^2 \leq \norm f_2^2 = \frac{1}{2\pi}\int_{-\pi}^{\pi} \abs f^2 :

the series cn2\sum \abs{c_n}^2 converges, and cn(f)0c_n(f) \to 0 as n\abs n \to \infty (Riemann–Lebesgue for coefficients).

Proof. fSN(f)f - S_N(f) is orthogonal to each eke_k, kN\abs k \leq N (ek,fSNf=ckck=0\langle e_k, f - S_N f\rangle = c_k - c_k = 0): SNfS_Nf is the orthogonal projection onto TN=Vect(eN,,eN)\mathcal{T}_N = \operatorname{Vect}(e_{-N}, \dots, e_N) (Year 1 volume’s projection theorem, verbatim in the Hermitian setting). Pythagoras: f22=SNf22+fSNf22SNf22=nNcn2\norm f_2^2 = \norm{S_Nf}_2^2 + \norm{f - S_Nf}_2^2 \geq \norm{S_Nf}_2^2 = \sum_{\abs n \leq N} \abs{c_n}^2; let NN \to \infty.

Example 14.3 (Best approximation, measured)

How well do low-degree trigonometric polynomials approximate the sawtooth f(t)=tf(t) = t (on (π,π)\intoo{-\pi}{\pi}) in the quadratic mean? By Proposition 14.2 the best degree-NN approximation is SN(f)S_N(f), with squared error

fSNf22=f22nNcn2.\norm{f - S_Nf}_2^2 = \norm f_2^2 - \sum_{\abs n\leq N}\abs{c_n}^2 .

Here f22=12πππt2 ⁣dt=π23\norm f_2^2 = \frac{1}{2\pi}\int_{-\pi}^\pi t^2\dd t = \frac{\pi^2}{3}, and from bn=2(1)n+1nb_n = \frac{2(-1)^{n+1}}{n} (Example 14.12): cn2+cn2=bn22=2n2\abs{c_n}^2 + \abs{c_{-n}}^2 = \frac{b_n^2}{2} = \frac{2}{n^2}. So

fSNf22=π23n=1N2n2: numerically 1.29, 0.79, 0.57, 0.44\norm{f - S_Nf}_2^2 = \frac{\pi^2}{3} - \sum_{n=1}^{N}\frac{2}{n^2} \qquad\text{: numerically } 1.29,\ 0.79,\ 0.57,\ 0.44

for N=1,2,3,4N = 1, 2, 3, 4 — decreasing, but slowly: the tail n>N2n22N\sum_{n>N}\frac2{n^2} \sim \frac2N is governed by the slow 1n\frac1n decay of the coefficients, itself the signature of the jump (Exercise 14.6 read backwards). Closing insight: Parseval turns approximation quality into a tail of a numerical series — and predicts, before any picture, that jumps make Fourier series converge reluctantly.

Method 14.4 (Computing Fourier coefficients efficiently)

Before integrating anything:

  1. Parity: even ff has bn=0b_n = 0, odd ff has an=0a_n = 0, and the surviving integrals reduce to 2π0π\frac2\pi \int_0^\pi — half the work, twice the reliability.
  2. Trigonometric polynomials are done already: linearize products (cos3\cos^3, sin2cos\sin^2\cos, …) and read the coefficients off (Exercise 14.9); orthonormality makes any further integration redundant.
  3. Complex exponentials for exponentials: for factors eat\eu^{at} or damped oscillations, compute cnc_n directly — one integral of e(ain)t\eu^{(a - \iu n)t} beats two integrations by parts (Exercise 14.10).
  4. Differentiate a known expansion: if ff' has known coefficients and ff is continuous, cn(f)=cn(f)inc_n(f) = \frac{c_n(f')}{\iu n} (n0n \neq 0) recovers all but c0c_0, which is the mean — often the fastest route, and legitimate exactly under the hypotheses of Theorem 14.10 (1).

14.2 Dirichlet’s theorem

Lemma 14.5 (Dirichlet kernel)

SN(f)(x)=12πππf(x+u)DN(u) ⁣duS_N(f)(x) = \frac{1}{2\pi}\int_{-\pi}^{\pi} f(x + u)\,D_N(u)\,\dd u, where

DN(u)=n=NNeinu=sin((N+12)u)sinu2(u2πZ),12πππDN=1.D_N(u) = \sum_{n=-N}^{N} \eu^{\iu nu} = \frac{\sin\bigl((N + \frac12)u\bigr)}{\sin\frac u2} \quad (u \notin 2\pi\Z), \qquad \frac{1}{2\pi}\int_{-\pi}^{\pi} D_N = 1 .

Proof. Insert the definition of cnc_n into SNS_N and swap sum and integral (legitimate: the sum is finite):

SN(f)(x)=n=NN(12πππf(t)eint ⁣dt)einx=12πππf(t)n=NNein(xt) ⁣dt;S_N(f)(x) = \sum_{n=-N}^{N}\Bigl(\frac{1}{2\pi}\int_{-\pi}^{\pi} f(t)\,\eu^{-\iu nt}\dd t\Bigr)\eu^{\iu nx} = \frac{1}{2\pi}\int_{-\pi}^{\pi} f(t) \sum_{n=-N}^{N}\eu^{\iu n(x-t)}\,\dd t ;

substitute u=txu = t - x and slide the integration segment back to [π,π]\intcc{-\pi}{\pi} by 2π2\pi-periodicity of the integrand; the symmetric index range makes neinu=DN(u)\sum_n\eu^{-\iu nu} = D_N(u). The closed form: geometric sum with ratio eiu\eu^{\iu u},

DN(u)=eiNuei(2N+1)u1eiu1=ei(N+12)uei(N+12)ueiu/2eiu/2,D_N(u) = \eu^{-\iu Nu}\,\frac{\eu^{\iu(2N+1)u} - 1}{\eu^{\iu u} - 1} = \frac{\eu^{\iu(N + \frac12)u} - \eu^{-\iu(N+\frac12)u}} {\eu^{\iu u/2} - \eu^{-\iu u/2}} ,

which is the sine quotient. Its mean is 11: only n=0n = 0 contributes.

Theorem 14.6 (Riemann–Lebesgue lemma)

For gg piecewise continuous on a segment, abg(t)sin(λt+φ) ⁣dt0\int_a^b g(t)\sin(\lambda t + \varphi)\,\dd t \to 0 as λ+\lambda \to +\infty.

Proof. Approximate gg uniformly by step functions (Theorem 10.16 is not needed — the elementary step-function approximation of piecewise continuous functions suffices) and integrate each step explicitly: every piece contributes O(1λ)O\bigl(\frac1\lambda\bigr), and the approximation error contributes ε(ba)\varepsilon(b - a). This argument was carried out in full as the final exercise of the Year 1 volume’s integration chapter; for C1C^1 pieces one can instead integrate by parts and bound by Cλ\frac C\lambda.

Example 14.7 (How fast do coefficients die?)

Riemann–Lebesgue says the coefficients tend to 00; their rate is a smoothness meter. Three specimens from this chapter and its exercises:

square wave: bn=4πn (n odd),t:an=4πn2 (n odd),sint:a2k=4π(4k21).\text{square wave: } b_n = \frac{4}{\pi n}\ (n\ \text{odd}), \qquad \abs t : a_n = \frac{-4}{\pi n^2}\ (n\ \text{odd}), \qquad \abs{\sin t} : a_{2k} = \frac{-4}{\pi(4k^2-1)} .

A jump in ff (square wave, sawtooth) leaves coefficients of order 1n\frac1n: no normal convergence, Gibbs overshoot at the jumps. Continuity with a corner — a jump in ff' only — improves the order to 1n2\frac{1}{n^2}: normal convergence, uniform reconstruction. In general kk derivatives buy cn=O(nk)c_n = O(n^{-k}) (Exercise 14.6), and conversely a spectrum decaying faster than every power forces ff to be CC^\infty (differentiate term by term, now legitimately). Closing insight: regularity of the signal and decay of the spectrum are the same information — an engineer reads one off the slope of the other without ever plotting the function.

Theorem 14.8 (Dirichlet)

Let ff be 2π2\pi-periodic and piecewise C1C^1. Then for every xx,

SN(f)(x)Nf(x+)+f(x)2S_N(f)(x) \xrightarrow[N \to \infty]{} \frac{f(x^+) + f(x^-)}{2}

(the mean of one-sided limits) — in particular SN(f)(x)f(x)S_N(f)(x) \to f(x) at every point of continuity.

Proof. By the kernel lemma and its unit mean, splitting the integral into u>0u > 0 and u<0u < 0 halves (each of mean 12\frac12):

SN(f)(x)f(x+)+f(x)2=12π0π(f(x+u)f(x+))DN(u) ⁣du+12ππ0(f(x+u)f(x))DN(u) ⁣du.\begin{align*} S_N(f)(x) - \frac{f(x^+) + f(x^-)}{2} &= \frac{1}{2\pi}\int_{0}^{\pi} \bigl(f(x+u) - f(x^+)\bigr)D_N(u)\,\dd u\\ &\quad+ \frac{1}{2\pi}\int_{-\pi}^{0}\bigl(f(x+u) - f(x^-)\bigr)D_N(u)\,\dd u . \end{align*}

Treat the first (the second is symmetric). Write

(f(x+u)f(x+))DN(u)=f(x+u)f(x+)sinu2g(u)sin((N+12)u).\bigl(f(x + u) - f(x^+)\bigr)\,D_N(u) = \underbrace{\frac{f(x+u) - f(x^+)}{\sin\frac u2}}_{g(u)}\, \sin\Bigl(\Bigl(N + \frac12\Bigr)u\Bigr) .

The function gg is piecewise continuous on (0,π]\intoc{0}{\pi} and has a finite limit at 0+0^+: writing

g(u)=f(x+u)f(x+)uusinu2,g(u) = \frac{f(x+u) - f(x^+)}{u}\cdot\frac{u}{\sin\frac u2} ,

the first factor tends to f(x+)f'(x^+) (one-sided differentiability, from piecewise C1C^1) and the second to 22 (the standard limit sinvv1\frac{\sin v}{v} \to 1 at v=u2v = \frac u2): g(0+)=2f(x+)g(0^+) = 2f'(x^+) exists. So gg extends piecewise continuously to [0,π]\intcc{0}{\pi}, and Riemann–Lebesgue (Theorem 14.6) sends the integral to 00. This is the whole point of the hypothesis: without one-sided derivatives, the factor 1sin(u/2)\frac{1}{\sin(u/2)} blows up at 00 faster than Riemann–Lebesgue can compensate, and pointwise convergence can genuinely fail for merely continuous ff — the gap that Fejér’s theorem (weekend problem) closes by averaging.

Example 14.9 (Dirichlet at a jump)

For the sawtooth f(t)=tf(t) = t on (π,π)\intoo{-\pi}{\pi} (Example 14.12 below), the periodic extension jumps at t=πt = \pi from f(π)=πf(\pi^-) = \pi to f(π+)=πf(\pi^+) = -\pi. Dirichlet promises the value π+(π)2=0\frac{\pi + (-\pi)}{2} = 0 there, and indeed every term of 2(1)n+1nsinnt\sum \frac{2(-1)^{n+1}}{n}\sin nt vanishes at t=πt = \pi: the series politely converges to the midpoint, ignoring both one-sided values. Moving the evaluation point to t=π2t = \frac\pi2 instead (a continuity point) turns the same series into Leibniz’s π4\frac\pi4. One series, two behaviors — exactly the two clauses of the theorem.

Theorem 14.10 (Normal convergence for C1C^1; Parseval)

  1. If ff is continuous, 2π2\pi-periodic and piecewise C1C^1, then cn(f)=incn(f)c_n(f') = \iu n\,c_n(f), the Fourier series of ff converges normally on R\R, and its sum is ff.
  2. (Parseval) For every piecewise continuous 2π2\pi-periodic ff:

    12πππf2=n=cn(f)2=a024+12n1(an2+bn2).\frac{1}{2\pi}\int_{-\pi}^{\pi}\abs{f}^2 = \sum_{n=-\infty}^{\infty} \abs{c_n(f)}^2 = \frac{\abs{a_0}^2}{4} + \frac12\sum_{n\geq1} \bigl(\abs{a_n}^2 + \abs{b_n}^2\bigr).

    (Proved here for ff continuous piecewise C1C^1; admitted in general.)

Proof. (1) Integration by parts on each C1C^1 piece (boundary terms cancel by continuity and periodicity): cn(f)=incn(f)c_n(f') = \iu n c_n(f). Then, by Cauchy–Schwarz on the two square-summable families (Proposition 14.2 for ff'):

n0cn(f)=n0cn(f)n(cn(f)2)1/2(n01n2)1/2<:\sum_{n \neq 0} \abs{c_n(f)} = \sum_{n\neq0} \frac{\abs{c_n(f')}}{\abs n} \leq \Bigl(\sum \abs{c_n(f')}^2\Bigr)^{1/2} \Bigl(\sum_{n\neq0}\frac{1}{n^2}\Bigr)^{1/2} < \infty :

normal convergence of the Fourier series. Its sum is continuous and coincides with ff at every point by Dirichlet (Theorem 14.8: ff is continuous): the series converges to ff, uniformly.

(2) For such ff: SNffS_N f \to f uniformly, so fSNf2fSNf0\norm{f - S_Nf}_2 \leq \norm{f - S_Nf}_\infty \to 0, and Pythagoras (f22=nNcn2+fSNf22\norm f_2^2 = \sum_{\abs n \leq N}\abs{c_n}^2 + \norm{f - S_Nf}_2^2) passes to the limit. The real form is bookkeeping with an,bna_n, b_n.

Example 14.11 (The C1C^1 tail bound, made quantitative)

The proof of Theorem 14.10 (1) hides a usable estimate. For continuous piecewise-C1C^1 ff, the same Cauchy–Schwarz applied to the tail only gives

n>Ncn(f)=n>Ncn(f)n(n>Ncn(f)2) ⁣1/2(n>N1n2) ⁣1/2f22N,\sum_{\abs n > N}\abs{c_n(f)} = \sum_{\abs n > N}\frac{\abs{c_n(f')}}{\abs n} \leq \Bigl(\sum_{\abs n > N}\abs{c_n(f')}^2\Bigr)^{\!1/2} \Bigl(\sum_{\abs n>N}\frac{1}{n^2}\Bigr)^{\!1/2} \leq \norm{f'}_2\,\sqrt{\frac{2}{N}} ,

using Bessel for ff' and n>Nn21N\sum_{n>N}n^{-2} \leq \frac1N. So the uniform error of the partial sums obeys

fSNfn>Ncn(f)f22N.\norm{f - S_Nf}_\infty \leq \sum_{\abs n>N}\abs{c_n(f)} \leq \norm{f'}_2\,\sqrt{\frac2N} .

For f(t)=tf(t) = \abs t: f2=1\norm{f'}_2 = 1 (the derivative is ±1\pm1), so ten terms already reconstruct t\abs t uniformly within 0.20.45\sqrt{0.2} \approx 0.45, and N=104N = 10^4 within 0.0150.015. Closing insight: one derivative buys the 1N\frac{1}{\sqrt N} uniform rate; comparing with the 1n\frac1n-coefficient world of the square wave (no uniform convergence at all), the dictionary of Example 14.7 acquires numbers.

Example 14.12 (Basel and friends)

Let f(t)=tf(t) = t on (π,π)\intoo{-\pi}{\pi}, extended 2π2\pi-periodically (a sawtooth, piecewise C1C^1). Computing, an=0a_n = 0 (oddness) and

bn=1πππtsinnt ⁣dt=2(1)n+1n.b_n = \frac{1}{\pi}\int_{-\pi}^{\pi} t\sin nt\,\dd t = \frac{2(-1)^{n+1}}{n} .

Dirichlet at t=π2t = \frac\pi2 recovers Leibniz’s π4=113+15\frac\pi4 = 1 - \frac13 + \frac15 - \dots; Parseval gives

12πππt2 ⁣dt=π23=12n14n2  n11n2=π26  \frac{1}{2\pi}\int_{-\pi}^{\pi} t^2\,\dd t = \frac{\pi^2}{3} = \frac12\sum_{n\geq1}\frac{4}{n^2} \quad\Longrightarrow\quad \boxed{\;\sum_{n\geq1}\frac{1}{n^2} = \frac{\pi^2}{6}\;}

— Euler’s Basel sum, in two lines. The function f(t)=t2f(t) = t^2 similarly yields 1n4=π490\sum \frac1{n^4} = \frac{\pi^4}{90} (Exercise 14.3).

Example 14.13 (A full expansion with built-in check: sint\abs{\sin t})

The function f(t)=sintf(t) = \abs{\sin t} is continuous, even, π\pi-periodic (hence 2π2\pi-periodic), piecewise C1C^1. Evenness kills the bnb_n; a0=1πππsint ⁣dt=4πa_0 = \frac1\pi\int_{-\pi}^{\pi} \abs{\sin t}\dd t = \frac4\pi; and for n1n \geq 1, product-to-sum gives

an=2π0πsintcosnt ⁣dt=1π0π(sin(1+n)t+sin(1n)t) ⁣dt=2π1+cosnπ1n2a_n = \frac2\pi\int_0^\pi \sin t\cos nt\,\dd t = \frac{1}{\pi}\int_0^\pi\bigl(\sin(1+n)t + \sin(1-n)t\bigr)\dd t = \frac2\pi\cdot\frac{1 + \cos n\pi}{1 - n^2}

for n1n \neq 1 (and a1=0a_1 = 0 directly): zero for odd nn, and a2k=4π(4k21)a_{2k} = \frac{-4}{\pi(4k^2-1)}. By Theorem 14.10 (1) the convergence is normal, and

sint=2π4πk1cos(2kt)4k21(tR).\abs{\sin t} = \frac{2}{\pi} - \frac{4}{\pi} \sum_{k\geq1}\frac{\cos(2kt)}{4k^2 - 1} \qquad (t \in \R) .

Built-in check at t=0t = 0: the identity demands k114k21=12\sum_{k\geq1}\frac{1}{4k^2-1} = \frac12, which telescoping confirms:

k114k21=12k1(12k112k+1)=12.\sum_{k\geq1}\frac{1}{4k^2-1} = \frac12\sum_{k\geq1}\Bigl(\frac{1}{2k-1} - \frac{1}{2k+1}\Bigr) = \frac12 . \checkmark

Closing insight: the spectrum of sin\abs{\sin} lives on the even frequencies only — rectifying a sine doubles its frequency content, which is why full-wave rectifiers hum at 100100 or 120120 hertz, twice the mains frequency.

Example 14.14 (Parseval as a computing device)

Parseval turns expansions into numerical series wholesale. Apply it to f(t)=tf(t) = \abs t (Exercise 14.2: a0=πa_0 = \pi, an=4πn2a_n = \frac{-4}{\pi n^2} for odd nn, the rest zero):

12πππt2 ⁣dt=π23=a024+12n oddan2=π24+8π2n odd1n4,\frac{1}{2\pi}\int_{-\pi}^{\pi}t^2\,\dd t = \frac{\pi^2}{3} = \frac{a_0^2}{4} + \frac12\sum_{n \text{ odd}} a_n^2 = \frac{\pi^2}{4} + \frac{8}{\pi^2}\sum_{n\text{ odd}}\frac{1}{n^4} ,

whence

n odd1n4=π28(π23π24)=π496.\sum_{n\text{ odd}}\frac{1}{n^4} = \frac{\pi^2}{8}\Bigl(\frac{\pi^2}{3} - \frac{\pi^2}{4}\Bigr) = \frac{\pi^4}{96} .

Cross-check against Exercise 14.3: splitting 1n4\sum\frac1{n^4} into odd and even parts gives ζ\zeta-style bookkeeping S=Sodd+S16S = S_{\mathrm{odd}} + \frac{S}{16}, so S=1615π496=π490S = \frac{16}{15}\cdot\frac{\pi^4}{96} = \frac{\pi^4}{90} — exactly the value found there by a different function. Two expansions, one number: the consistency is Parseval’s isometry at work. Closing insight: each new Fourier expansion is a generating machine for series identities; the weekend problem’s Part II explains why the machine can never contradict itself.

The square wave (gray) and the Fourier partial sums S_1 (blue) and S_9 (red): convergence at every continuity point, but a persistent 9\% overshoot near the jumps — the Gibbs phenomenon. Uniform convergence fails exactly because the limit is discontinuous.
The square wave (gray) and the Fourier partial sums S1S_1 (blue) and S9S_9 (red): convergence at every continuity point, but a persistent 9%\sim 9\% overshoot near the jumps — the Gibbs phenomenon. Uniform convergence fails exactly because the limit is discontinuous.

Example 14.15 (Translation and modulation)

Two one-line rules generate many expansions from one. For aRa \in \R, substituting s=tas = t - a:

cn(f(a))=12πππf(ta)eint ⁣dt=einacn(f)(translation modulates the spectrum),c_n\bigl(f(\cdot - a)\bigr) = \frac{1}{2\pi}\int_{-\pi}^{\pi}f(t - a)\eu^{-\iu nt}\dd t = \eu^{-\iu na}\,c_n(f) \qquad\text{(translation modulates the spectrum)},

and directly from the definition,

cn(eiktf)=cnk(f)(modulation translates the spectrum).c_n\bigl(\eu^{\iu kt}f\bigr) = c_{n-k}(f) \qquad\text{(modulation translates the spectrum)}.

Worked instance: the sawtooth shifted by π\pi, g(t)=f(tπ)g(t) = f(t - \pi) with f(t)=tf(t) = t, has bnb_n-coefficients (1)n2(1)n+1n=2n(-1)^n\cdot\frac{2(-1)^{n+1}}{n} = -\frac2n: the expansion g2sinntng \sim -2\sum\frac{\sin nt}{n} of the sawtooth that jumps at 00 instead of at π\pi — no integral recomputed. Closing insight: time shifts only turn phases, never amplitudes (cn\abs{c_n} is shift-invariant), which is why energy (Parseval) and convergence class are properties of the signal shape, not of where the clock starts.

Remark 14.16 (Common pitfalls)

(i) Three convergences, three currencies: pointwise (Dirichlet: needs piecewise C1C^1, pays the midpoint at each jump — never the one-sided value), uniform (needs a continuous limit; impossible across a jump, Gibbs is the visible symptom), and quadratic mean (Parseval: the most robust, blind to individual points). Always name which one you are claiming. (ii) No termwise differentiation by default: differentiating the sawtooth series of Example 14.12 term by term yields 2(1)n+1cosnt\sum 2(-1)^{n+1}\cos nt, whose terms do not even tend to 00 — the transfer theorems of the function-sequences chapter need uniform convergence of the derived series, which the jump destroys. Smoothness first, differentiation second (Exercise 14.6 is the dictionary). (iii) Symmetric partial sums: Dirichlet’s theorem concerns SN=NNS_N = \sum_{-N}^{N}; rearranging or summing one side first can change divergence into convergence and back. (iv) Normalization drift: conventions differ across books (12π\frac{1}{2\pi} or 1π\frac1\pi in front, period 2π2\pi or 11); the reliable invariants are the orthonormality relations — recompute em,en\langle e_m, e_n\rangle in the convention at hand before trusting any formula.

Remark 14.17 (Where this is used)

Parseval is the germ of the L2L^2 theory of Fourier series: the Year 3 volume completes the picture (the exponentials are a Hilbert basis of L2L^2, and the map f(cn)f \mapsto (c_n) is a bijective isometry). Within this volume, the weekend problem proves Fejér’s theorem — the Cesàro means of the Fourier series converge uniformly for every continuous periodic ff — which upgrades the admitted general Parseval to a theorem, yields the trigonometric Weierstrass theorem, and pays two spectacular dividends: Weyl’s equidistribution theorem and the isoperimetric inequality. Applied mathematics reads this chapter daily: spectra of signals, harmonics of vibrating systems, and the fast Fourier transform (the finite avatar was Exercise 13.10).

Remark 14.18 (Perspectives within this volume)

Three chapters converse with this one. Backward: the Hermitian chapter supplied the geometry (orthonormal families, projections, Bessel), and the function-sequences chapter the analysis (uniform convergence, transfer theorems, approximate identities — Fejér’s kernel is to Fourier series what Bernstein’s polynomials were to Weierstrass). Sideways: the power-series chapter’s boundary theory returns through Abel summation, realized here by the Poisson kernel (Exercise 14.12) — the disk’s radius rr playing the role of the summation parameter. Forward: the differential-equations chapter decomposes periodic forcing into harmonics and feeds each to the oscillator’s frequency response; resonance happens when a Fourier mode of the input matches a natural frequency, which is why its weekend problem and this chapter’s are two halves of one story.

The Dirichlet kernel D_8 (blue) oscillates and takes negative values; the Fejér kernel F_8 (red) is nonnegative, concentrates at 0, and has mean 1: an approximate identity. Positivity is exactly what the Dirichlet kernel lacks, and what makes Fejér’s theorem in the weekend problem unconditional.
The Dirichlet kernel D8D_8 (blue) oscillates and takes negative values; the Fejér kernel F8F_8 (red) is nonnegative, concentrates at 00, and has mean 11: an approximate identity. Positivity is exactly what the Dirichlet kernel lacks, and what makes Fejér’s theorem in the weekend problem unconditional.

14.3 Exercises

Exercise 14.1

Compute the Fourier coefficients of the square wave (f=1f = -1 on (π,0)\intoo{-\pi}{0}, +1+1 on (0,π)\intoo{0}{\pi}), state Dirichlet’s conclusion at t=π2t = \frac\pi2 and at the jump t=0t = 0, and recover Leibniz’s series.

Solution

Solution of Exercise 14.1.

Oddness kills the ana_n. For n1n \geq 1:

bn=2π0πsinnt ⁣dt=2π1(1)nn={4πnn odd,0n even.b_n = \frac{2}{\pi}\int_0^{\pi} \sin nt\,\dd t = \frac{2}{\pi}\cdot\frac{1 - (-1)^n}{n} = \begin{cases} \frac{4}{\pi n} & n \text{ odd},\\ 0 & n \text{ even}. \end{cases}

So S(f)(t)=4πk0sin((2k+1)t)2k+1S(f)(t) = \frac{4}{\pi}\sum_{k\geq0} \frac{\sin\bigl((2k+1)t \bigr)}{2k+1}. Dirichlet at t=π2t = \frac\pi2 (a continuity point, value 11): sin((2k+1)π2)=(1)k\sin\bigl((2k+1)\frac\pi2\bigr) = (-1)^k, giving

1=4πk0(1)k2k+1k0(1)k2k+1=π4(Leibniz).1 = \frac4\pi \sum_{k\geq0}\frac{(-1)^k}{2k+1} \quad\Longrightarrow\quad \sum_{k\geq0}\frac{(-1)^k}{2k+1} = \frac{\pi}{4} \quad\text{(Leibniz)}.

At the jump t=0t = 0: the series sums to 0=f(0+)+f(0)20 = \frac{f(0^+) + f(0^-)}{2}, as Dirichlet prescribes (every term vanishes: consistent).

Exercise 14.2

Expand f(t)=tf(t) = \abs t (tπ\abs t \leq \pi, 2π2\pi-periodic) in Fourier series; justify normal convergence; evaluate at t=0t = 0 to obtain k01(2k+1)2=π28\sum_{k\geq0}\frac{1}{(2k+1)^2} = \frac{\pi^2}{8}, and re-derive Basel from it.

Solution

Solution of Exercise 14.2.

Evenness kills the bnb_n; a0=1πππt ⁣dt=πa_0 = \frac{1}{\pi}\int_{-\pi}^\pi\abs t\,\dd t = \pi, and for n1n \geq 1:

an=2π0πtcosnt ⁣dt=2π(1)n1n2={4πn2n odd,0n even,a_n = \frac{2}{\pi}\int_0^\pi t\cos nt\,\dd t = \frac{2}{\pi}\cdot\frac{(-1)^n - 1}{n^2} = \begin{cases} -\frac{4}{\pi n^2} & n \text{ odd},\\ 0 & n \text{ even}, \end{cases}

(one integration by parts). Hence

t=π24πk0cos((2k+1)t)(2k+1)2,\abs t = \frac{\pi}{2} - \frac{4}{\pi}\sum_{k\geq0} \frac{\cos\bigl((2k+1)t\bigr)}{(2k+1)^2} ,

with normal convergence ((2k+1)2<\sum (2k+1)^{-2} < \infty) — as Theorem 14.10 (1) predicts for this continuous piecewise-C1C^1 function. At t=0t = 0:

0=π24πk01(2k+1)2k01(2k+1)2=π28.0 = \frac\pi2 - \frac4\pi\sum_{k\geq0}\frac{1}{(2k+1)^2} \quad\Longrightarrow\quad \sum_{k\geq0}\frac{1}{(2k+1)^2} = \frac{\pi^2}{8} .

Splitting 1n2\sum \frac{1}{n^2} into odd and even parts: S=π28+S4S = \frac{\pi^2}{8} + \frac S4, so S=π26S = \frac{\pi^2}{6}: Basel again.

Exercise 14.3

Expand f(t)=t2f(t) = t^2 (tπ\abs t \leq \pi) and deduce

n1(1)n+1n2=π212,n11n4=π490(Parseval).\sum_{n\geq1}\frac{(-1)^{n+1}}{n^2} = \frac{\pi^2}{12}, \qquad \sum_{n\geq1}\frac{1}{n^4} = \frac{\pi^4}{90} \quad\text{(Parseval)}.
Solution

Solution of Exercise 14.3.

Evenness: bn=0b_n = 0; a0=1πππt2=2π23a_0 = \frac{1}{\pi}\int_{-\pi}^{\pi} t^2 = \frac{2\pi^2}{3}; two integrations by parts give an=4(1)nn2a_n = \frac{4(-1)^n}{n^2} (n1n \geq 1). Hence

t2=π23+4n1(1)nn2cosnt(tπ),t^2 = \frac{\pi^2}{3} + 4\sum_{n\geq1} \frac{(-1)^n}{n^2}\cos nt \qquad (\abs t \leq \pi),

normally convergent. At t=0t = 0: 0=π23+4(1)nn20 = \frac{\pi^2}{3} + 4\sum\frac{(-1)^n}{n^2}, i.e. n1(1)n+1n2=π212\sum_{n\geq1} \frac{(-1)^{n+1}}{n^2} = \frac{\pi^2}{12}. Parseval:

12πππt4 ⁣dt=π45=a024+12n1an2=π49+8n11n4,\frac{1}{2\pi}\int_{-\pi}^{\pi} t^4\,\dd t = \frac{\pi^4}{5} = \frac{a_0^2}{4} + \frac12\sum_{n\geq1} a_n^2 = \frac{\pi^4}{9} + 8\sum_{n\geq1}\frac{1}{n^4} ,

so 1n4=18(π45π49)=π490\sum \frac{1}{n^4} = \frac18\bigl(\frac{\pi^4}{5} - \frac{\pi^4}{9}\bigr) = \frac{\pi^4}{90}.

Exercise 14.4 ★★

Let ff be continuous 2π2\pi-periodic with cn(f)=0c_n(f) = 0 for all nn. Prove f=0f = 0 (Parseval — for which class is it proved here? justify that continuity plus piecewise C1C^1 can be dropped by admitting general Parseval, or give the density argument in outline).

Solution

Solution of Exercise 14.4.

If moreover ff is piecewise C1C^1: Parseval (proved) gives f22=cn2=0\norm f_2^2 = \sum\abs{c_n}^2 = 0, and strict positivity of the integral of the continuous f2\abs f^2 forces f=0f = 0.

For merely continuous ff, admit general Parseval: same one-line proof. (Outline of the density route: Fejér/Weierstrass-type trigonometric approximation shows trigonometric polynomials are 2\norm\cdot_2-dense among continuous periodic functions; since ff \perp all of them, f22=f,fPf2fP2\norm f_2^2 = \langle f, f - P\rangle \leq \norm f_2\norm{f - P}_2 for approximants PP, forcing f2=0\norm f_2 = 0.)

Exercise 14.5 ★★

For αZ\alpha \notin \Z, expand f(t)=cos(αt)f(t) = \cos(\alpha t) (tπ\abs t \leq \pi) and deduce the partial-fraction expansion of the cotangent:

πcot(πα)=1α+n12αα2n2.\pi\cot(\pi\alpha) = \frac{1}{\alpha} + \sum_{n\geq1} \frac{2\alpha}{\alpha^2 - n^2} .
Solution

Solution of Exercise 14.5.

Evenness: bn=0b_n = 0;

an=2π0πcos(αt)cos(nt) ⁣dt=2π(1)nαsin(πα)α2n2a_n = \frac{2}{\pi}\int_0^\pi \cos(\alpha t)\cos(nt)\,\dd t = \frac{2}{\pi}\cdot \frac{(-1)^n\,\alpha\sin(\pi\alpha)}{\alpha^2 - n^2}

(product-to-sum, then integrate; a0=2sin(πα)παa_0 = \frac{2\sin(\pi\alpha)}{\pi\alpha}). Dirichlet at t=πt = \pi (a continuity point of the periodic extension, whose one-sided values agree by evenness):

cos(πα)=sin(πα)πα+n12αsin(πα)π(α2n2)(1)ncos(nπ)=sin(πα)π(1α+n12αα2n2),\cos(\pi\alpha) = \frac{\sin(\pi\alpha)}{\pi\alpha} + \sum_{n\geq1} \frac{2\alpha\sin(\pi\alpha)}{\pi(\alpha^2 - n^2)}\,(-1)^n\cos(n\pi) = \frac{\sin(\pi\alpha)}{\pi}\Bigl(\frac{1}{\alpha} + \sum_{n\geq1}\frac{2\alpha}{\alpha^2 - n^2}\Bigr),

using (1)ncosnπ=1(-1)^n\cos n\pi = 1. Dividing by sin(πα)/π\sin(\pi\alpha)/\pi gives the cotangent expansion.

Exercise 14.6 ★★

Prove that if ff is 2π2\pi-periodic and CkC^k with f(k)f^{(k)} piecewise continuous, then cn(f)=O(nk)c_n(f) = O\bigl(\abs n^{-k}\bigr): smoothness of the signal == decay of its spectrum.

Solution

Solution of Exercise 14.6.

Iterating cn(f)=incn(f)c_n(f') = \iu n\,c_n(f) (kk times, integration by parts across the CkC^{k} pieces with matching boundary values): cn(f(k))=(in)kcn(f)c_n(f^{(k)}) = (\iu n)^k c_n(f). The coefficients of the piecewise continuous f(k)f^{(k)} are bounded (indeed 0\to 0, Bessel):

cn(f)=cn(f(k))nk=O(nk).\abs{c_n(f)} = \frac{\abs{c_n(f^{(k)})}}{\abs n^k} = O\bigl(\abs n^{-k}\bigr) .

Exercise 14.7 ★★★

(Wirtinger’s inequality) Let ff be C1C^1, 2π2\pi-periodic, with ππf=0\int_{-\pi}^{\pi} f = 0. Prove

ππf2ππf2,\int_{-\pi}^{\pi} \abs{f}^2 \leq \int_{-\pi}^{\pi} \abs{f'}^2 ,

with equality iff f(t)=acost+bsintf(t) = a\cos t + b\sin t. (Parseval on both sides; compare cn2\abs{c_n}^2 and n2cn2n^2\abs{c_n}^2.)

Solution

Solution of Exercise 14.7.

Parseval for ff and for ff' (both legitimate: ff is C1C^1, ff' piecewise continuous — indeed continuous):

12πf2=n0cn2(c0=0 by the mean-zero hypothesis),12πf2=n0n2cn2.\frac{1}{2\pi}\int \abs f^2 = \sum_{n\neq0} \abs{c_n}^2 \quad (c_0 = 0 \text{ by the mean-zero hypothesis}), \qquad \frac{1}{2\pi}\int \abs{f'}^2 = \sum_{n\neq0} n^2\abs{c_n}^2 .

Termwise, n2cn2cn2n^2\abs{c_n}^2 \geq \abs{c_n}^2 for n1\abs n \geq 1: the inequality follows. Equality forces (n21)cn2=0(n^2 - 1)\abs{c_n}^2 = 0 for all nn, i.e. cn=0c_n = 0 for n2\abs n \geq 2: f(t)=c1eit+c1eit=acost+bsintf(t) = c_1\eu^{\iu t} + c_{-1}\eu^{-\iu t} = a\cos t + b\sin t (real form); conversely such ff give equality.

Exercise 14.8 ★★★

(The Gibbs constant) For the square wave of Exercise 14.1, evaluate the partial sum at xN=π2Nx_N = \frac{\pi}{2N}: writing uk=(2k+1)π2Nu_k = \frac{(2k+1)\pi}{2N} and Δu=πN\Delta u = \frac{\pi}{N}, show that

S2N1(π2N)=4πk=0N1sinuk2k+1=2πk=0N1sinukukΔuN2π0πsinuu ⁣du1.179:S_{2N-1}\Bigl(\frac{\pi}{2N}\Bigr) = \frac{4}{\pi}\sum_{k=0}^{N-1} \frac{\sin u_k}{2k+1} = \frac{2}{\pi}\sum_{k=0}^{N-1} \frac{\sin u_k}{u_k}\,\Delta u \xrightarrow[N\to\infty]{} \frac{2}{\pi}\int_0^{\pi}\frac{\sin u}{u}\,\dd u \approx 1.179 :

a Riemann sum of 2πsinuu\frac{2}{\pi}\cdot\frac{\sin u}{u} on [0,π]\intcc{0}{\pi} at midpoints. Conclude that the overshoot beyond the jump value 11 does not vanish as NN \to \infty.

Solution

Solution of Exercise 14.8.

From Exercise 14.1, S2N1(x)=4πk=0N1sin((2k+1)x)2k+1S_{2N-1}(x) = \frac4\pi\sum_{k=0}^{N-1}\frac{\sin((2k+1)x)}{2k+1}. At xN=π2Nx_N = \frac{\pi}{2N}, with uk=(2k+1)xN=(2k+1)π2Nu_k = (2k+1)x_N = \frac{(2k+1)\pi}{2N}:

S2N1(xN)=4πk=0N1sinuk2k+1=4πk=0N1sinukukuk2k+1=2πk=0N1sinukukπN,S_{2N-1}(x_N) = \frac{4}{\pi}\sum_{k=0}^{N-1}\frac{\sin u_k}{2k+1} = \frac{4}{\pi}\sum_{k=0}^{N-1}\frac{\sin u_k}{u_k}\cdot \frac{u_k}{2k+1} = \frac{2}{\pi}\sum_{k=0}^{N-1}\frac{\sin u_k}{u_k}\cdot \frac{\pi}{N},

since uk2k+1=π2N\frac{u_k}{2k+1} = \frac{\pi}{2N}. The points uku_k are the midpoints of the NN subintervals of [0,π]\intcc{0}{\pi} of length πN\frac{\pi}{N}: the sum is a midpoint Riemann sum of the continuous usinuuu \mapsto \frac{\sin u}{u}, hence converges to

2π0πsinuu ⁣du2π×1.85191.179.\frac{2}{\pi}\int_0^\pi \frac{\sin u}{u}\,\dd u \approx \frac{2}{\pi}\times 1.8519 \approx 1.179 .

The partial sums near the jump overshoot the value 11 by 18%\approx 18\% of the half-jump forever: Gibbs’s phenomenon, quantified.

Exercise 14.9

Expand cos3t\cos^3 t and sin2tcost\sin^2 t\,\cos t in Fourier series (linearize; a trigonometric polynomial is its own Fourier series, by uniqueness of coefficients). What are cnc_n, ana_n, bnb_n for each?

Solution

Solution of Exercise 14.9.

From cos3t=4cos3t3cost\cos 3t = 4\cos^3t - 3\cos t:

cos3t=3cost+cos3t4,sin2tcost=costcos3t=costcos3t4.\cos^3 t = \frac{3\cos t + \cos 3t}{4}, \qquad \sin^2t\,\cos t = \cos t - \cos^3 t = \frac{\cos t - \cos 3t}{4} .

Each is a trigonometric polynomial, hence equal to its own Fourier series (uniqueness of coefficients: two expansions would differ by a trigonometric polynomial with all coefficients zero). For cos3t\cos^3t: a1=34a_1 = \frac34, a3=14a_3 = \frac14, all other ana_n and all bnb_n zero; c±1=38c_{\pm1} = \frac38, c±3=18c_{\pm3} = \frac18. For sin2tcost\sin^2t\cos t: a1=14a_1 = \frac14, a3=14a_3 = -\frac14; c±1=18c_{\pm1} = \frac18, c±3=18c_{\pm3} = -\frac18.

Exercise 14.10 ★★

Let a>0a > 0 and f(t)=eatf(t) = \eu^{at} on (π,π]\intoc{-\pi}{\pi}, extended 2π2\pi-periodically. Compute

cn(f)=(1)nsinh(aπ)π(ain),c_n(f) = \frac{(-1)^n\sinh(a\pi)}{\pi(a - \iu n)} ,

apply Dirichlet’s theorem at the jump t=πt = \pi, and deduce the partial-fraction expansion of the hyperbolic cotangent:

coth(πa)=1πa+n12aπ(a2+n2).\coth(\pi a) = \frac{1}{\pi a} + \sum_{n\geq1} \frac{2a}{\pi(a^2 + n^2)} .
Solution

Solution of Exercise 14.10.

Direct computation:

cn=12πππe(ain)t ⁣dt=e(ain)πe(ain)π2π(ain)=(1)nsinh(aπ)π(ain),c_n = \frac{1}{2\pi}\int_{-\pi}^{\pi}\eu^{(a - \iu n)t}\dd t = \frac{\eu^{(a-\iu n)\pi} - \eu^{-(a - \iu n)\pi}} {2\pi(a - \iu n)} = \frac{(-1)^n\sinh(a\pi)}{\pi(a - \iu n)} ,

using e±inπ=(1)n\eu^{\pm\iu n\pi} = (-1)^n. At t=πt = \pi the periodic extension jumps from eaπ\eu^{a\pi} to eaπ\eu^{-a\pi}; Dirichlet (symmetric partial sums) gives

cosh(aπ)=nZ(1)ncn=sinh(aπ)π(1a+n1(1ain+1a+in))=sinh(aπ)π(1a+n12aa2+n2),\cosh(a\pi) = \sum_{n\in\Z}(-1)^n c_n\, = \frac{\sinh(a\pi)}{\pi}\Bigl(\frac1a + \sum_{n\geq1}\Bigl(\frac{1}{a - \iu n} + \frac{1}{a + \iu n}\Bigr)\Bigr) = \frac{\sinh(a\pi)}{\pi}\Bigl(\frac1a + \sum_{n\geq1}\frac{2a}{a^2 + n^2}\Bigr) ,

the imaginary parts of the paired terms cancelling. Divide by sinh(aπ)\sinh(a\pi):

coth(πa)=1πa+n12aπ(a2+n2),\coth(\pi a) = \frac{1}{\pi a} + \sum_{n\geq1}\frac{2a}{\pi(a^2 + n^2)} ,

the hyperbolic twin of Exercise 14.5.

Exercise 14.11 ★★

(Convolution) For f,gf, g continuous and 2π2\pi-periodic define

(fg)(x)=12πππf(xt)g(t) ⁣dt.(f * g)(x) = \frac{1}{2\pi}\int_{-\pi}^{\pi} f(x - t)\,g(t)\, \dd t .

Show that fg=gff * g = g * f, that cn(fg)=cn(f)cn(g)c_n(f*g) = c_n(f)\,c_n(g) (to swap the two integrals of a continuous integrand, compare the two functions of the upper limit: both vanish at the left endpoint and have the same derivative, by continuity and differentiation under the integral sign), and that SN(f)=fDNS_N(f) = f * D_N for the Dirichlet kernel. (The weekend problem’s Fejér means are likewise convolutions, σN(f)=fFN\sigma_N(f) = f * F_N.)

Solution

Solution of Exercise 14.11.

Commutativity: substitute s=xts = x - t and use periodicity of the integrand. For cn(fg)c_n(f * g), the integrand (x,t)f(xt)g(t)einx(x, t) \mapsto f(x-t)g(t)\eu^{-\iu nx} is continuous; the two iterated integrals agree (both, as functions of the upper limit of the outer variable, vanish at the left endpoint and have the same derivative — continuity plus Theorem 9.10 justify differentiating the iterated integral). Hence

cn(fg)=12πππg(t)eint(12πππf(xt)ein(xt) ⁣dx) ⁣dt=cn(f)cn(g),c_n(f*g) = \frac{1}{2\pi}\int_{-\pi}^{\pi} g(t)\,\eu^{-\iu nt} \Bigl(\frac{1}{2\pi}\int_{-\pi}^{\pi} f(x-t)\,\eu^{-\iu n(x-t)}\dd x\Bigr)\dd t = c_n(f)\,c_n(g),

the inner integral being cn(f)c_n(f) for every tt (substitution and periodicity). Finally Lemma 14.5 says SN(f)(x)=12πf(x+u)DN(u) ⁣duS_N(f)(x) = \frac{1}{2\pi}\int f(x+u)D_N(u)\dd u; the substitution utu \mapsto -t and the evenness of DND_N turn this into (fDN)(x)(f * D_N)(x).

Exercise 14.12 ★★★

(Poisson kernel: Abel means of Fourier series) For 0r<10 \leq r < 1 set Pr(t)=nZrneintP_r(t) = \sum_{n\in\Z} r^{\abs n}\eu^{\iu nt}.

  1. Sum the two geometric series and show

    Pr(t)=1r212rcost+r2>0,12πππPr=1.P_r(t) = \frac{1 - r^2}{1 - 2r\cos t + r^2} > 0, \qquad \frac{1}{2\pi}\int_{-\pi}^{\pi}P_r = 1 .
  2. Show that for δtπ\delta \leq \abs t \leq \pi: Pr(t)1r212rcosδ+r20P_r(t) \leq \frac{1 - r^2}{1 - 2r\cos\delta + r^2} \to 0 as r1r \to 1^-, uniformly.
  3. Deduce that for every continuous 2π2\pi-periodic ff, the Abel means (fPr)(x)=nrncn(f)einx(f * P_r)(x) = \sum_n r^{\abs n}c_n(f)\,\eu^{\iu nx} converge to ff uniformly as r1r \to 1^- — the continuous-parameter sibling of Fejér’s theorem, and the Fourier incarnation of Abel summation from the power-series chapter.
Solution

Solution of Exercise 14.12.

  1. With w=reitw = r\eu^{\iu t}:

    Pr(t)=1+2w1w=1+w1w=1w21w2=1r212rcost+r2>0.P_r(t) = 1 + 2\,\Re\frac{w}{1 - w} = \Re\frac{1 + w}{1 - w} = \frac{1 - \abs w^2}{\abs{1 - w}^2} = \frac{1 - r^2}{1 - 2r\cos t + r^2} > 0 .

    Mean 11: term-by-term integration of the normally convergent series keeps only n=0n = 0.

  2. For δtπ\delta \leq \abs t \leq \pi: costcosδ\cos t \leq \cos\delta, so Pr(t)1r212rcosδ+r2P_r(t) \leq \frac{1 - r^2}{1 - 2r\cos\delta + r^2}, whose denominator tends to 22cosδ>02 - 2\cos\delta > 0 while the numerator tends to 00: uniform convergence to 00 off any neighborhood of 00.
  3. Term-by-term integration (normal convergence in tt) gives (fPr)(x)=nrncn(f)einx(f * P_r)(x) = \sum_n r^{\abs n}c_n(f)\eu^{\iu nx}. The approximate-identity argument: with mean 11 and positivity,

    (fPr)(x)f(x)12πππf(xt)f(x)Pr(t) ⁣dt,\abs{(f*P_r)(x) - f(x)} \leq \frac{1}{2\pi}\int_{-\pi}^{\pi} \abs{f(x-t) - f(x)}\,P_r(t)\,\dd t ,

    split at t=δ\abs t = \delta: at most ε\varepsilon (Heine) plus 2fsupδtπPrε2\norm f_\infty\sup_{\delta\leq\abs t\leq\pi}P_r \to \varepsilon: uniform convergence as r1r \to 1^-. This is Abel summation of the Fourier series — the Fourier twin of the power-series chapter’s boundary theory.

14.4 Problem: Fejér’s theorem and its dividends

Problem 14.1

Dirichlet’s theorem needs ff piecewise C1C^1; for merely continuous ff the partial sums SN(f)S_N(f) can misbehave. Fejér’s discovery: their Cesàro means never do. The engine is the positivity of the Fejér kernel, and the harvest is immense: uniform trigonometric approximation (Weierstrass), uniqueness of Fourier coefficients, Parseval for every function of this chapter (removing the “admitted” in Theorem 14.10), Weyl’s equidistribution theorem, and — crowning a century of geometry — the isoperimetric inequality. Throughout, ff is 2π2\pi-periodic and piecewise continuous, and

σN(f)=S0(f)+S1(f)++SN1(f)N.\sigma_N(f) = \frac{S_0(f) + S_1(f) + \dots + S_{N-1}(f)}{N} .

Part I — The Fejér kernel.

  1. Show that σN(f)(x)=12πππf(x+u)FN(u) ⁣du\sigma_N(f)(x) = \frac{1}{2\pi}\int_{-\pi}^{\pi}f(x+u)\,F_N(u)\,\dd u with FN=D0++DN1NF_N = \frac{D_0 + \dots + D_{N-1}}{N}, and that 12πππFN=1\frac{1}{2\pi}\int_{-\pi}^{\pi}F_N = 1.
  2. Prove the closed form, for u2πZu \notin 2\pi\Z:

    FN(u)=1Nsin2(Nu2)sin2(u2)    0F_N(u) = \frac{1}{N}\, \frac{\sin^2\bigl(\frac{Nu}{2}\bigr)} {\sin^2\bigl(\frac u2\bigr)} \;\geq\; 0

    (sum sin((n+12)u)\sin\bigl((n+\frac12)u\bigr) as the imaginary part of a geometric series).

  3. Show the concentration estimate: for 0<δuπ0 < \delta \leq \abs u \leq \pi,

    FN(u)1Nsin2δ2N0uniformly:F_N(u) \leq \frac{1}{N\sin^2\frac\delta2} \xrightarrow[N\to\infty]{} 0 \quad\text{uniformly} :

    (FN)(F_N) is a positive approximate identity.

  4. (Fejér’s theorem) Prove: if ff is continuous and 2π2\pi-periodic, then σN(f)f\sigma_N(f) \to f uniformly on R\R (split the integral of (f(x+u)f(x))FN(u)\bigl(f(x+u) - f(x)\bigr)F_N(u) at u=δ\abs u = \delta; use Heine and questions 1–3).
  5. For piecewise continuous ff, show the pointwise version σN(f)(x)f(x+)+f(x)2\sigma_N(f)(x) \to \frac{f(x^+) + f(x^-)}{2} at every xx, and the uniform bound σN(f)f\norm{\sigma_N(f)}_\infty \leq \norm f_\infty (positivity!).

Part II — Weierstrass, uniqueness, Parseval.

  1. (Trigonometric Weierstrass) Deduce: every continuous 2π2\pi-periodic function is a uniform limit of trigonometric polynomials.
  2. (Uniqueness) Deduce: a continuous ff with cn(f)=0c_n(f) = 0 for all nn is identically zero — two continuous periodic functions with the same Fourier coefficients coincide (Exercise 14.4, now with no admission).
  3. (Parseval, continuous case) Using the projection property of SNS_N (Proposition 14.2) and σNfTN1TN\sigma_Nf \in \mathcal T_{N-1} \subseteq \mathcal T_N, prove

    fSNf2fσNf2fσNfN0,\norm{f - S_Nf}_2 \leq \norm{f - \sigma_Nf}_2 \leq \norm{f - \sigma_Nf}_\infty \xrightarrow[N\to\infty]{} 0 ,

    and conclude Parseval’s identity for every continuous 2π2\pi-periodic ff.

  4. (Parseval, piecewise continuous case) Given ff piecewise continuous and ε>0\varepsilon > 0, construct a continuous periodic gg with fg2ε\norm{f - g}_2 \leq \varepsilon (replace ff by an affine interpolation on tiny intervals around the jumps), and deduce fSNf20\norm{f - S_Nf}_2 \to 0 (use Bessel: SNh2h2\norm{S_Nh}_2 \leq \norm h_2): Parseval holds in the full generality stated in Theorem 14.10 — the “admitted” is gone.
  5. (No Gibbs for Fejér) Contrast with Exercise 14.8: show that for the square wave ff, σN(f)1\abs{\sigma_N(f)} \leq 1 everywhere, for every NN — Cesàro averaging erases the overshoot that haunts SNS_N. Explain in one sentence which property of FNF_N is responsible.

Part III — Rates.

  1. Prove the two kernel bounds, for 0<uπ0 < \abs u \leq \pi:

    FN(u)N,FN(u)π2Nu2F_N(u) \leq N, \qquad F_N(u) \leq \frac{\pi^2}{N u^2}

    (for the first, sinNθNsinθ\abs{\sin N\theta} \leq N\abs{\sin\theta} by induction; for the second, sinu2uπ\sin\frac u2 \geq \frac{u}{\pi} on [0,π]\intcc{0}{\pi}).

  2. Deduce the first-moment estimate

    12πππuFN(u) ⁣du    ClnNN(N2)\frac{1}{2\pi}\int_{-\pi}^{\pi}\abs u\,F_N(u)\,\dd u \;\leq\; \frac{C\,\ln N}{N} \qquad (N \geq 2)

    for an explicit constant (split at u=1N\abs u = \frac1N).

  3. Conclude: if ff is LL-Lipschitz and 2π2\pi-periodic, then

    σNffCLlnNN.\norm{\sigma_N f - f}_\infty \leq \frac{C\,L\ln N}{N} .
  4. (Saturation) Compute σN(e1)\sigma_N(e_1) for e1(t)=eite_1(t) = \eu^{\iu t} and show σNe1e1=1N\norm{\sigma_Ne_1 - e_1}_\infty = \frac1N: even for the smoothest functions, Fejér converges no faster than 1N\frac1N — the exact analogue of the Bernstein saturation in the function-sequences chapter’s weekend problem.
  5. (Localization) Show: if ff (piecewise continuous) vanishes on (xδ,x+δ)\intoo{x - \delta}{x + \delta}, then σN(f)(x)0\sigma_N(f)(x) \to 0, however wild ff is elsewhere — convergence of the means at xx only sees ff near xx.

Part IV — Weyl’s equidistribution theorem. A sequence (xn)n1(x_n)_{n\geq1} in [0,1)\intco{0}{1} is equidistributed when, for every interval [a,b][0,1]\intcc{a}{b} \subseteq \intcc{0}{1},

#{nN:xn[a,b]}NNba.\frac{\#\{n \leq N : x_n \in \intcc ab\}}{N} \xrightarrow[N\to\infty]{} b - a .
  1. Show that (xn)(x_n) is equidistributed as soon as 1NnNf(xn)01f\frac1N\sum_{n\leq N}f(x_n) \to \int_0^1 f for every continuous 11-periodic ff (squeeze the indicator of [a,b]\intcc ab between two continuous piecewise-affine functions whose integrals differ by ε\varepsilon).
  2. (Weyl’s criterion, sufficiency) Suppose

    1Nn=1Ne2iπkxnN0for every kZ{0}.\frac{1}{N}\sum_{n=1}^{N}\eu^{2\iu\pi kx_n} \xrightarrow[N\to\infty]{} 0 \qquad\text{for every } k \in \Z\setminus\{0\} .

    Show 1Nf(xn)01f\frac1N\sum f(x_n) \to \int_0^1f first for trigonometric polynomials, then for all continuous 11-periodic ff by question 6 (transported to period 11): with question 16, (xn)(x_n) is equidistributed.

  3. Let α\alpha be irrational and xn={nα}x_n = \{n\alpha\} (fractional part). Bound the geometric sum

    n=1Ne2iπknα21e2iπkα(k0),\Bigl|\sum_{n=1}^{N}\eu^{2\iu\pi kn\alpha}\Bigr| \leq \frac{2}{\abs{1 - \eu^{2\iu\pi k\alpha}}} \qquad (k \neq 0),

    and conclude Weyl’s theorem: ({nα})(\{n\alpha\}) is equidistributed in [0,1)\intco{0}{1}.

  4. Deduce that ({nα})(\{n\alpha\}) is dense in [0,1]\intcc{0}{1} for irrational α\alpha, and explain in one sentence why equidistribution is strictly stronger than density.
  5. (Leading digits) Prove that the proportion of integers nNn \leq N such that 2n2^n has leading (decimal) digit 11 tends to log1020.301\log_{10}2 \approx 0.301 (leading digit 11 means {nlog102}[0,log102)\{n\log_{10}2\} \in \intco{0}{\log_{10}2}; show log102\log_{10}2 is irrational).

Part V — The isoperimetric inequality. Let Γ\Gamma be a closed simple C1C^1 curve of length LL enclosing a signed area AA, parametrized by scaled arc length: z(t)=x(t)+iy(t)z(t) = x(t) + \iu y(t), 2π2\pi-periodic, with z(t)=L2π\abs{z'(t)} = \frac{L}{2\pi} constant; the enclosed area is

A=1202π(xyyx) ⁣dt=1202πzz ⁣dtA = \frac12\int_0^{2\pi}\bigl(x\,y' - y\,x'\bigr)\dd t = \frac{1}{2}\,\Im\int_0^{2\pi}\conj{z}\,z'\,\dd t

(taken here as the definition of the signed area; the chapter on multiple integrals proves it agrees with the intuitive one, via Green’s formula).

  1. Expand z(t)=nZcneintz(t) = \sum_{n\in\Z}c_n\eu^{\iu nt} (the series of a C1C^1 function, normally convergent) and prove, by Parseval applied to zz':

    L22π=02πz2 ⁣dt=2πnZn2cn2.\frac{L^2}{2\pi} = \int_0^{2\pi}\abs{z'}^2\dd t = 2\pi\sum_{n\in\Z}n^2\abs{c_n}^2 .
  2. Prove likewise A=πnZncn2A = \pi\sum_{n\in\Z} n\,\abs{c_n}^2 (Parseval in its polarized form: 12πfg=cn(f)cn(g)\frac{1}{2\pi}\int\conj f g = \sum \conj{c_n(f)}c_n(g), applied to f=zf = z, g=zg = z').
  3. (Hurwitz) Conclude:

    L24πA=4π2nZ(n2n)cn2    0,L^2 - 4\pi A = 4\pi^2\sum_{n\in\Z} (n^2 - n)\abs{c_n}^2 \;\geq\; 0 ,

    with equality if and only if z(t)=c0+c1eitz(t) = c_0 + c_1\eu^{\iu t} — a circle. The isoperimetric inequality: among closed curves of length LL, only the circle encloses area L24π\frac{L^2}{4\pi}.

  4. Sanity checks: verify equality for the circle of radius RR and the strict inequality for the square of side aa; explain why n2n0n^2 - n \geq 0 for every integer nn, including the negative ones, and where the parametrization’s constant speed was used.
  5. Synthesis. In one sentence each: (i) the single property of FNF_N from which Parts I–III flow, and which DND_N lacks; (ii) how Cesàro summation here relates to the power-series chapter’s weekend problem (Frobenius); (iii) which dividend used only Weierstrass (question 6) and which needed full Parseval; (iv) one sentence on what the Year 3 volume adds (L2L^2 completeness: Fourier series as a Hilbert basis).
Solution

Solution of Problem 14.1.

1. Averaging Lemma 14.5 over n=0,,N1n = 0, \dots, N-1 (linearity of the integral) gives σN(f)(x)=12πf(x+u)FN(u) ⁣du\sigma_N(f)(x) = \frac{1}{2\pi}\int f(x+u)F_N(u)\dd u; each DnD_n has mean 11, so FNF_N has mean 11.

2. With eiu1=2ieiu/2sinu2\eu^{\iu u} - 1 = 2\iu\,\eu^{\iu u/2}\sin\frac u2:

n=0N1sin((n+12)u)=[eiu/2eiNu1eiu1]=1eiNu2sinu2=1cosNu2sinu2=sin2Nu2sinu2.\sum_{n=0}^{N-1}\sin\Bigl(\Bigl(n + \frac12\Bigr)u\Bigr) = \Im\Bigl[\eu^{\iu u/2}\,\frac{\eu^{\iu Nu} - 1}{\eu^{\iu u} - 1}\Bigr] = \Re\,\frac{1 - \eu^{\iu Nu}}{2\sin\frac u2} = \frac{1 - \cos Nu}{2\sin\frac u2} = \frac{\sin^2\frac{Nu}2}{\sin\frac u2} .

Dividing by Nsinu2N\sin\frac u2:

FN(u)=1Nn=0N1sin((n+12)u)sinu2=1Nsin2Nu2sin2u20.F_N(u) = \frac1N\sum_{n=0}^{N-1} \frac{\sin\bigl((n+\frac12)u\bigr)}{\sin\frac u2} = \frac{1}{N}\, \frac{\sin^2\frac{Nu}{2}}{\sin^2\frac u2} \geq 0 .

3. On δuπ\delta \leq \abs u \leq \pi: sin2u2sin2δ2\sin^2\frac u2 \geq \sin^2\frac\delta2 and sin2Nu21\sin^2\frac{Nu}2 \leq 1: FN1Nsin2(δ/2)0F_N \leq \frac{1}{N\sin^2(\delta/2)} \to 0, uniformly there.

4. By the unit mean, σNf(x)f(x)=12π(f(x+u)f(x))FN(u) ⁣du\sigma_Nf(x) - f(x) = \frac{1}{2\pi}\int\bigl(f(x+u) - f(x)\bigr)F_N(u)\dd u. Given ε\varepsilon, Heine yields δ\delta with f(x+u)f(x)ε\abs{f(x+u) - f(x)} \leq \varepsilon for uδ\abs u \leq \delta, uniformly in xx. Then, using FN0F_N \geq 0 and its unit mean,

σNf(x)f(x)ε+2f12πδuπFNε+2fNsin2δ22ε\abs{\sigma_Nf(x) - f(x)} \leq \varepsilon + 2\norm f_\infty\cdot\frac{1}{2\pi} \int_{\delta\leq\abs u\leq\pi}F_N \leq \varepsilon + \frac{2\norm f_\infty}{N\sin^2\frac\delta2} \leq 2\varepsilon

for large NN, uniformly in xx: Fejér’s theorem.

5. FNF_N is even with mean 11: each half [0,π]\intcc{0}{\pi}, [π,0]\intcc{-\pi}{0} carries mean 12\frac12. Then

σNf(x)f(x+)+f(x)2=12π0π(f(x+u)f(x+))FN ⁣du+12ππ0(f(x+u)f(x))FN ⁣du;\sigma_Nf(x) - \frac{f(x^+)+f(x^-)}{2} = \frac{1}{2\pi}\int_0^\pi\bigl(f(x+u) - f(x^+)\bigr)F_N\,\dd u + \frac{1}{2\pi}\int_{-\pi}^0\bigl(f(x+u) - f(x^-)\bigr)F_N\,\dd u ;

on each half, split at u=δ\abs u = \delta where the one-sided limit is ε\varepsilon-close, and let question 3 kill the far part: both integrals tend to 00. The bound: FN0F_N \geq 0 gives σNf(x)12πf(x+u)FNf\abs{\sigma_Nf(x)} \leq \frac{1}{2\pi}\int\abs{f(x+u)}F_N \leq \norm f_\infty.

6. Each σN(f)\sigma_N(f) is a trigonometric polynomial (an average of the Sn(f)S_n(f), n<Nn < N), and σN(f)f\sigma_N(f) \to f uniformly: the trigonometric Weierstrass theorem.

7. cn(f)=0c_n(f) = 0 for all nn makes every Sn(f)=0S_n(f) = 0, hence every σN(f)=0\sigma_N(f) = 0; by Fejér, f=limσNf=0f = \lim\sigma_Nf = 0. Applying this to a difference: continuous periodic functions are determined by their Fourier coefficients.

8. SNfS_Nf is the orthogonal projection of ff onto TN\mathcal T_N (Proposition 14.2), so it minimizes fP2\norm{f - P}_2 over PTNP \in \mathcal T_N; since σNfTN1TN\sigma_Nf \in \mathcal T_{N-1} \subseteq \mathcal T_N:

fSNf2fσNf2fσNf0\norm{f - S_Nf}_2 \leq \norm{f - \sigma_Nf}_2 \leq \norm{f - \sigma_Nf}_\infty \to 0

(the middle inequality because the mean of 2\abs\cdot^2 is at most the sup squared). Pythagoras f22=nNcn2+fSNf22\norm f_2^2 = \sum_{\abs n\leq N}\abs{c_n}^2 + \norm{f - S_Nf}_2^2 then passes to the limit: Parseval for every continuous periodic ff.

9. Let t1,,tpt_1, \dots, t_p be the jumps of ff in a period, M=fM = \norm f_\infty. For small η\eta, define g=fg = f outside the intervals (tjη,tj+η)\intoo{t_j - \eta}{t_j + \eta} and by the affine chord across each such interval: gg is continuous, periodic, gM\norm g_\infty \leq M, and

fg2212πp(2M)22ηε2\norm{f - g}_2^2 \leq \frac{1}{2\pi}\,p\cdot(2M)^2\cdot2\eta \leq \varepsilon^2

for η\eta small. Bessel makes SNS_N a contraction for 2\norm\cdot_2, so

fSNf2fg2+gSNg2+SN(gf)22ε+gSNg2,\norm{f - S_Nf}_2 \leq \norm{f - g}_2 + \norm{g - S_Ng}_2 + \norm{S_N(g - f)}_2 \leq 2\varepsilon + \norm{g - S_Ng}_2 ,

and question 8 gives lim supNfSNf22ε\limsup_N\norm{f - S_Nf}_2 \leq 2\varepsilon for every ε\varepsilon: fSNf20\norm{f - S_Nf}_2 \to 0, and Pythagoras yields Parseval for every piecewise continuous ff: the “admitted” in Theorem 14.10 is now a theorem.

10. The square wave has f=1\norm f_\infty = 1, so question 5 gives σN(f)1\abs{\sigma_N(f)} \leq 1 everywhere and for every NN — no overshoot, ever — while Exercise 14.8 shows supxS2N1(f)(x)1.179\sup_xS_{2N-1}(f)(x) \to \approx 1.179. The one responsible property: FN0F_N \geq 0, so σNf(x)\sigma_Nf(x) is a weighted average of values of ff and can never leave [minf,maxf]\intcc{\min f}{\max f}; DND_N takes negative values, so SNS_N can.

11. sinNθNsinθ\abs{\sin N\theta} \leq N\abs{\sin\theta} by induction (sin(N+1)θsinNθcosθ+cosNθsinθ(N+1)sinθ\abs{\sin(N{+}1)\theta} \leq \abs{\sin N\theta}\abs{\cos\theta} + \abs{\cos N\theta}\abs{\sin\theta} \leq (N+1)\abs{\sin\theta}): with θ=u2\theta = \frac u2,

FN(u)=sin2Nu2Nsin2u2N2sin2u2Nsin2u2=N.F_N(u) = \frac{\sin^2\frac{Nu}2}{N\sin^2\frac u2} \leq \frac{N^2\sin^2\frac u2}{N\sin^2\frac u2} = N .

Concavity of sin\sin on [0,π2]\intcc{0}{\frac\pi2} gives sinu2uπ\sin\frac u2 \geq \frac{u}{\pi} for 0uπ0 \leq u \leq \pi, so FN(u)1N(u/π)2=π2Nu2F_N(u) \leq \frac{1}{N(u/\pi)^2} = \frac{\pi^2}{Nu^2}.

12. By evenness and the split at 1N\frac1N:

12πππuFN=1π0πuFN1π(01/NuN ⁣du+1/Nππ2Nu ⁣du)=12πN+πln(πN)N.\frac{1}{2\pi}\int_{-\pi}^{\pi}\abs uF_N = \frac1\pi\int_0^\pi uF_N \leq \frac1\pi\Bigl(\int_0^{1/N}uN\,\dd u + \int_{1/N}^{\pi}\frac{\pi^2}{Nu}\,\dd u\Bigr) = \frac{1}{2\pi N} + \frac{\pi\ln(\pi N)}{N} .

For N2N \geq 2: ln(πN)(1+lnπln2)lnN2.66lnN\ln(\pi N) \leq \bigl(1 + \frac{\ln\pi}{\ln2}\bigr)\ln N \leq 2.66\ln N and 12πNlnNN\frac{1}{2\pi N} \leq \frac{\ln N}{N}, so the moment is 9lnNN\leq \frac{9\ln N}{N}.

13. For LL-Lipschitz ff:

σNf(x)f(x)12πf(x+u)f(x)FN(u) ⁣duL12πuFN9LlnNN,\abs{\sigma_Nf(x) - f(x)} \leq \frac{1}{2\pi}\int\abs{f(x+u) - f(x)}F_N(u)\dd u \leq L\cdot\frac{1}{2\pi}\int\abs uF_N \leq \frac{9L\ln N}{N},

uniformly in xx.

14. c0(e1)=0c_0(e_1) = 0 gives S0(e1)=0S_0(e_1) = 0, while Sn(e1)=e1S_n(e_1) = e_1 for n1n \geq 1: σN(e1)=N1Ne1\sigma_N(e_1) = \frac{N-1}{N}e_1, so σNe1e1=1N\norm{\sigma_Ne_1 - e_1}_\infty = \frac1N. Even for this entire, band-limited signal the rate is 1N\frac1N: Fejér saturates, exactly as Bernstein’s operator saturates at 1n\frac1n (Voronovskaya, in the function-sequences chapter’s weekend problem).

15. If ff vanishes on (xδ,x+δ)\intoo{x-\delta}{x+\delta}, then σNf(x)=12πδuπf(x+u)FN(u) ⁣du\sigma_Nf(x) = \frac{1}{2\pi}\int_{\delta \leq \abs u \leq \pi}f(x+u)F_N(u)\dd u, of absolute value at most fNsin2(δ/2)0\frac{\norm f_\infty}{N\sin^2(\delta/2)} \to 0: the Cesàro means at xx only see ff near xx.

16. Given [a,b]\intcc ab and ε\varepsilon, choose continuous 11-periodic piecewise-affine φ±\varphi^\pm with φ1[a,b]φ+\varphi^- \leq \mathbf 1_{\intcc ab} \leq \varphi^+ and 01(φ+φ)ε\int_0^1(\varphi^+ - \varphi^-) \leq \varepsilon (trapezoids with slopes over intervals of total length ε\varepsilon). Then

lim supN#{nN:xn[a,b]}NlimN1NnNφ+(xn)=01φ+ba+ε,\limsup_N\frac{\#\{n \leq N : x_n \in \intcc ab\}}{N} \leq \lim_N\frac1N\sum_{n\leq N}\varphi^+(x_n) = \int_0^1\varphi^+ \leq b - a + \varepsilon ,

and symmetrically lim infbaε\liminf \geq b - a - \varepsilon: the proportion tends to bab - a.

17. For f=e2iπkf = \eu^{2\iu\pi k\cdot} with k0k \neq 0 the hypothesis gives the limit 0=01f0 = \int_0^1f; for k=0k = 0 both sides are 11; linearity handles every trigonometric polynomial. For continuous 11-periodic ff and ε>0\varepsilon > 0, question 6 (transported by t=2πxt = 2\pi x) provides a trigonometric polynomial PP with fPε\norm{f - P}_\infty \leq \varepsilon:

1NnNf(xn)01f2ε+1NnNP(xn)01P2ε.\Bigl|\frac1N\sum_{n\leq N}f(x_n) - \int_0^1f\Bigr| \leq 2\varepsilon + \Bigl|\frac1N\sum_{n\leq N}P(x_n) - \int_0^1P\Bigr| \longrightarrow 2\varepsilon .

With question 16: (xn)(x_n) is equidistributed.

18. For k0k \neq 0 and irrational α\alpha, w=e2iπkα1w = \eu^{2\iu\pi k\alpha} \neq 1:

n=1Nwn=w(wN1)w121w,\Bigl|\sum_{n=1}^{N}w^n\Bigr| = \Bigl|\frac{w(w^N - 1)}{w - 1}\Bigr| \leq \frac{2}{\abs{1 - w}} ,

a bound independent of NN; dividing by NN gives Weyl’s criterion, and question 17 concludes: ({nα})(\{n\alpha\}) is equidistributed.

19. Every subinterval receives asymptotic proportion equal to its length, in particular infinitely many points: ({nα})(\{n\alpha\}) is dense. Equidistribution is stronger: a sequence can be dense while spending nearly all its time in one corner (density says where the sequence goes, equidistribution says how often).

20. 2n2^n has leading digit 11 iff 10m2n<210m10^m \leq 2^n < 2\cdot10^m for some mm, i.e. iff {nlog102}[0,log102)\{n\log_{10}2\} \in \intco{0}{\log_{10}2}. Irrationality: log102=pq\log_{10}2 = \frac pq would give 2q=10p=2p5p2^q = 10^p = 2^p5^p, impossible for p1p \geq 1 by unique factorization. Weyl’s theorem (question 18 with α=log102\alpha = \log_{10}2; the half-open interval is squeezed between closed ones of nearby lengths) gives the proportion log1020.301\log_{10}2 \approx 0.301: the first digits of 2n2^n follow Benford’s law.

21. zz is C1C^1, so its Fourier series converges normally with sum zz (Theorem 14.10 (1)), and cn(z)=incnc_n(z') = \iu n\,c_n. Since z=L2π\abs{z'} = \frac{L}{2\pi} is constant,

02πz2 ⁣dt=2π(L2π) ⁣2=L22π,\int_0^{2\pi}\abs{z'}^2\dd t = 2\pi\Bigl(\frac{L}{2\pi}\Bigr)^{\!2} = \frac{L^2}{2\pi} ,

and Parseval applied to the continuous zz' gives 12π02πz2=nincn2\frac{1}{2\pi}\int_0^{2\pi}\abs{z'}^2 = \sum_n\abs{\iu nc_n}^2, i.e. L22π=2πnn2cn2\frac{L^2}{2\pi} = 2\pi\sum_n n^2\abs{c_n}^2.

22. The polarized Parseval 12πfg=cn(f)cn(g)\frac{1}{2\pi}\int\conj fg = \sum\conj{c_n(f)}c_n(g) follows from Parseval applied to f+gf + g and f+igf + \iu g (polarization identity), both continuous. With f=zf = z, g=zg = z':

A=1202πzz=πncn(incn)=πnncn2.A = \frac12\,\Im\int_0^{2\pi}\conj z\,z' = \pi\,\Im\sum_n\conj{c_n}(\iu n c_n) = \pi\sum_n n\abs{c_n}^2 .

23. Combining questions 21–22:

L24πA=4π2nn2cn24π2nncn2=4π2nZ(n2n)cn20,L^2 - 4\pi A = 4\pi^2\sum_n n^2\abs{c_n}^2 - 4\pi^2\sum_n n\abs{c_n}^2 = 4\pi^2\sum_{n\in\Z}(n^2 - n)\abs{c_n}^2 \geq 0 ,

since n2n=n(n1)0n^2 - n = n(n-1) \geq 0 for every integer. Equality forces cn=0c_n = 0 for all n{0,1}n \notin \{0, 1\}: z(t)=c0+c1eitz(t) = c_0 + c_1\eu^{\iu t}, a circle of center c0c_0 and radius c1=L2π\abs{c_1} = \frac{L}{2\pi} (constant speed). Hurwitz’s proof of the isoperimetric inequality: AL24πA \leq \frac{L^2}{4\pi}, circle only.

24. Circle of radius RR: L=2πRL = 2\pi R, A=πR2A = \pi R^2: L2=4π2R2=4πAL^2 = 4\pi^2R^2 = 4\pi A: equality. Square of side aa: L2=16a2>4πa2=4πAL^2 = 16a^2 > 4\pi a^2 = 4\pi A (as 16>4π12.5716 > 4\pi \approx 12.57). For negative nn, n2n=n(n1)n^2 - n = n(n - 1) is a product of two negative integers: positive — so backward-winding modes cost area twice. Constant speed entered in question 21, converting z2\int\abs{z'}^2 into L22π\frac{L^2}{2\pi}; for a non-constant speed, Cauchy–Schwarz gives z2(z)22π=L22π\int\abs{z'}^2 \geq \frac{(\int\abs{z'})^2}{2\pi} = \frac{L^2}{2\pi}, so the inequality survives, with the circle still the only equality case.

25. (i) Everything flows from FN0F_N \geq 0 (with unit mean and concentration); DND_N has unit mean and concentration of oscillation but not positivity, and Gibbs is the price. (ii) Cesàro summability of the Fourier series implies its Abel summability with the same sum (Frobenius, proved in the power-series chapter’s weekend problem) — the Poisson-kernel route of Exercise 14.12 is exactly Abel’s method. (iii) Weyl’s theorem needed only uniform approximation (question 6); the isoperimetric inequality needed Parseval itself (questions 8, 21–22). (iv) The Year 3 volume proves completeness: the exponentials form a Hilbert basis of L2L^2, Parseval becomes an isometry of Hilbert spaces, and Fejér’s theorem becomes the statement that this isometry is computable by positive averages.