Lemma数学ロードマップ

08 微分方程式 · 第 6 章

フーリエ級数

目安 8〜11 時間定理など 13演習 8 問実験 1

この章の目標

  • フーリエ係数を計算でき、ベッセルの不等式とリーマン–ルベーグの補題を証明できる
  • ディリクレ核を用いて、区分的に滑らかな関数のフーリエ級数の各点収束を証明できる
  • フェイエールの定理とパーセヴァルの等式を証明し、三角関数系の完全性を説明できる
  • 滑らかさと係数の減衰の関係を理解し、一様収束・項別微分の条件を述べられる
  • バーゼル問題・等周不等式・ワイエルシュトラスの近似定理にフーリエ級数を応用できる

前提:第5章、微分積分学 第5章、第6章

フーリエは 1807 年頃、棒の中の熱の伝わり方を研究し、「任意の」関数が三角関数の和 ∑(ancos⁡nx+bnsin⁡nx)\sum(a_n\cos nx + b_n\sin nx) で表せると主張した。熱方程式を変数分離すると sin⁡nx\sin nx という固有関数が現れる(第5章、第7章)ので、初期温度分布をこれらの和に展開できれば、熱方程式が解けるからである。この主張がどのような意味で正しいかを明らかにする試みから、関数の概念、リーマン積分、一様収束、集合論、ルベーグ積分が生まれた。本章では、フーリエ級数が「どの関数について」「どの意味で」元の関数を表すかを調べる。

以下、関数は R\mathbb{R} 上の周期 2π2\pi の複素数値関数で、[−π,π][-\pi, \pi] 上リーマン積分可能なものとする(このような関数は有界である)。内積とノルムを

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

とおく(∥f∥2=0\lVert f \rVert_2 = 0 でも f=0f = 0 とは限らないので、厳密には半ノルムである)。en(x)=einxe_n(x) = e^{inx} とすると ⟨em,en⟩\langle e_m, e_n \rangle は m=nm = n なら 1、そうでなければ 0 である。

6.1 三角級数とフーリエ係数

定義 6.1(フーリエ級数, Fourier series)n∈Zn \in \mathbb{Z} に対し

f^(n)=⟨f,en⟩=12π∫−ππf(x)e−inx dx\hat{f}(n) = \langle f, e_n \rangle = \frac{1}{2\pi}\int_{-\pi}^\pi f(x)e^{-inx}\,dx

を ff のフーリエ係数、形式的な級数 ∑n∈Zf^(n)einx\sum_{n \in \mathbb{Z}}\hat{f}(n)e^{inx} をフーリエ級数といい、f∼∑f^(n)einxf \sim \sum\hat{f}(n)e^{inx} と書く。部分和を SNf(x)=∑∣n∣≤Nf^(n)einxS_Nf(x) = \sum_{\lvert n \rvert \leq N}\hat{f}(n)e^{inx} とおく。

実数値関数では、an=1π∫−ππf(x)cos⁡nx dxa_n = \frac{1}{\pi}\int_{-\pi}^\pi f(x)\cos nx\ dx, bn=1π∫−ππf(x)sin⁡nx dxb_n = \frac{1}{\pi}\int_{-\pi}^\pi f(x)\sin nx\ dx とおくと f^(±n)=(an∓ibn)/2\hat{f}(\pm n) = (a_n \mp ib_n)/2 で、SNf=a02+∑n=1N(ancos⁡nx+bnsin⁡nx)S_Nf = \frac{a_0}{2} + \sum_{n=1}^N(a_n\cos nx + b_n\sin nx) である。ff が奇関数なら an=0a_n = 0、偶関数なら bn=0b_n = 0。

例 6.2 周期 2π2\pi に拡張した次の関数を考える(部分積分で計算できる)。

  • のこぎり波 f(x)=xf(x) = x(−π<x<π-\pi < x < \pi):bn=2π∫0πxsin⁡nx dx=2(−1)n+1nb_n = \frac{2}{\pi}\int_0^\pi x\sin nx\ dx = \frac{2(-1)^{n+1}}{n} で、x∼2∑n≥1(−1)n+1nsin⁡nxx \sim 2\sum_{n \geq 1}\frac{(-1)^{n+1}}{n}\sin nx。
  • 方形波 f(x)=sgn⁡xf(x) = \operatorname{sgn}x:bn=2π∫0πsin⁡nx dx=2(1−(−1)n)πnb_n = \frac{2}{\pi}\int_0^\pi\sin nx\ dx = \frac{2(1 - (-1)^n)}{\pi n} で、sgn⁡x∼4π∑k≥0sin⁡(2k+1)x2k+1\operatorname{sgn}x \sim \frac{4}{\pi}\sum_{k \geq 0}\frac{\sin(2k + 1)x}{2k + 1}。
  • f(x)=x2f(x) = x^2(∣x∣≤π\lvert x \rvert \leq \pi):a0=2π23a_0 = \frac{2\pi^2}{3}、an=2π∫0πx2cos⁡nx dx=4(−1)nn2a_n = \frac{2}{\pi}\int_0^\pi x^2\cos nx\ dx = \frac{4(-1)^n}{n^2} で、x2∼π23+4∑n≥1(−1)nn2cos⁡nxx^2 \sim \frac{\pi^2}{3} + 4\sum_{n \geq 1}\frac{(-1)^n}{n^2}\cos nx。

注意 6.3(区間 [0,ℓ][0, \ell] 上の正弦級数・余弦級数)[0,ℓ][0, \ell] 上の関数 ff を奇関数として [−ℓ,ℓ][-\ell, \ell] に拡張し、周期 2ℓ2\ell に延長して変数を πx/ℓ\pi x/\ell に取り替えると、正弦級数 f(x)∼∑n≥1bnsin⁡nπxℓf(x) \sim \sum_{n \geq 1}b_n\sin\frac{n\pi x}{\ell}, bn=2ℓ∫0ℓf(x)sin⁡nπxℓdxb_n = \frac{2}{\ell}\int_0^\ell f(x)\sin\frac{n\pi x}{\ell}dx が得られる。偶関数として拡張すれば余弦級数 a02+∑ancos⁡nπxℓ\frac{a_0}{2} + \sum a_n\cos\frac{n\pi x}{\ell}, an=2ℓ∫0ℓf(x)cos⁡nπxℓdxa_n = \frac{2}{\ell}\int_0^\ell f(x)\cos\frac{n\pi x}{\ell}dx が得られる。これらは −y′′=λy-y'' = \lambda y のディリクレ条件・ノイマン条件の固有関数展開(例 5.8)であり、第7章で使う。奇関数拡張が連続になるのは f(0)=f(ℓ)=0f(0) = f(\ell) = 0 のときである。

6.2 ベッセルの不等式とリーマン–ルベーグの補題

命題 6.4(最良近似)TN\mathcal{T}_N を ene_n(∣n∣≤N\lvert n \rvert \leq N)の一次結合(NN 次以下の三角多項式)全体とする。

∥f−SNf∥22=∥f∥22−∑∣n∣≤N∣f^(n)∣2\lVert f - S_Nf \rVert_2^2 = \lVert f \rVert_2^2 - \sum_{\lvert n \rvert \leq N}\lvert \hat{f}(n) \rvert^2

であり、任意の P∈TNP \in \mathcal{T}_N に対し ∥f−SNf∥2≤∥f−P∥2\lVert f - S_Nf \rVert_2 \leq \lVert f - P \rVert_2 が成り立つ。

証明. ∣n∣≤N\lvert n \rvert \leq N なら ⟨f−SNf,en⟩=f^(n)−f^(n)=0\langle f - S_Nf, e_n \rangle = \hat{f}(n) - \hat{f}(n) = 0 なので、f−SNff - S_Nf は TN\mathcal{T}_N と直交する。ピタゴラスの定理より ∥f∥22=∥f−SNf∥22+∥SNf∥22\lVert f \rVert_2^2 = \lVert f - S_Nf \rVert_2^2 + \lVert S_Nf \rVert_2^2 で、∥SNf∥22=∑∣n∣≤N∣f^(n)∣2\lVert S_Nf \rVert_2^2 = \sum_{\lvert n \rvert \leq N}\lvert \hat{f}(n) \rvert^2。また SNf−P∈TNS_Nf - P \in \mathcal{T}_N なので ∥f−P∥22=∥f−SNf∥22+∥SNf−P∥22\lVert f - P \rVert_2^2 = \lVert f - S_Nf \rVert_2^2 + \lVert S_Nf - P \rVert_2^2。□\square

SNfS_Nf は ff の TN\mathcal{T}_N への直交射影であり、L2L^2 の意味での最良近似である。

系 6.5(ベッセルの不等式, Bessel's inequality)∑n∈Z∣f^(n)∣2≤12π∫−ππ∣f(x)∣2dx\sum_{n \in \mathbb{Z}}\lvert \hat{f}(n) \rvert^2 \leq \frac{1}{2\pi}\int_{-\pi}^\pi\lvert f(x) \rvert^2dx。特に f^(n)→0\hat{f}(n) \to 0(∣n∣→∞\lvert n \rvert \to \infty)。

より一般に次が成り立つ。

定理 6.6(リーマン–ルベーグの補題, Riemann–Lebesgue lemma)ff が [a,b][a, b] 上リーマン積分可能なら、λ∈R\lambda \in \mathbb{R} について

lim⁡∣λ∣→∞∫abf(x)eiλx dx=0\lim_{\lvert \lambda \rvert \to \infty}\int_a^b f(x)e^{i\lambda x}\,dx = 0

証明. 区間の特性関数 1[c,d]\mathbf{1}_{[c, d]} については ∣∫cdeiλxdx∣=∣eiλd−eiλciλ∣≤2∣λ∣\left\lvert \int_c^d e^{i\lambda x}dx \right\rvert = \left\lvert \frac{e^{i\lambda d} - e^{i\lambda c}}{i\lambda} \right\rvert \leq \frac{2}{\lvert \lambda \rvert} なので、階段関数 s=∑jmj1Ijs = \sum_j m_j\mathbf{1}_{I_j} については積分は 2∑j∣mj∣/∣λ∣→02\sum_j\lvert m_j \rvert/\lvert \lambda \rvert \to 0 である。ε>0\varepsilon > 0 に対し、上ダルブー和と下ダルブー和の差が ε\varepsilon より小さい分割をとり、各小区間で ff の下限をとる階段関数を ss とすると、∫ab∣f−s∣ dx<ε\int_a^b\lvert f - s \rvert\ dx < \varepsilon。よって lim sup⁡∣∫abfeiλxdx∣≤ε+lim sup⁡∣∫abseiλxdx∣=ε\limsup\left\lvert \int_a^b fe^{i\lambda x}dx \right\rvert \leq \varepsilon + \limsup\left\lvert \int_a^b se^{i\lambda x}dx \right\rvert = \varepsilon。□\square

速く振動する関数との積の積分は、正負が打ち消し合って 0 に近づく。

6.3 ディリクレ核と各点収束

補題 6.7 SNf(x)=12π∫−ππf(x−t)DN(t) dtS_Nf(x) = \frac{1}{2\pi}\int_{-\pi}^\pi f(x - t)D_N(t)\ dt。ここでディリクレ核 (Dirichlet kernel) は

DN(t)=∑∣n∣≤Neint=sin⁡(N+12)tsin⁡t2D_N(t) = \sum_{\lvert n \rvert \leq N}e^{int} = \frac{\sin\left(N + \frac{1}{2}\right)t}{\sin\frac{t}{2}}

であり、偶関数で 12π∫−ππDN(t) dt=1\frac{1}{2\pi}\int_{-\pi}^\pi D_N(t)\ dt = 1 をみたす。

証明. SNf(x)=∑∣n∣≤N12π∫f(y)ein(x−y)dy=12π∫f(y)DN(x−y) dyS_Nf(x) = \sum_{\lvert n \rvert \leq N}\frac{1}{2\pi}\int f(y)e^{in(x - y)}dy = \frac{1}{2\pi}\int f(y)D_N(x - y)\ dy で、t=x−yt = x - y と置換し周期性を使う。等比数列の和から DN(t)=ei(N+1)t−e−iNteit−1D_N(t) = \frac{e^{i(N + 1)t} - e^{-iNt}}{e^{it} - 1} で、分子分母に e−it/2e^{-it/2} を掛ければ右の形になる。積分が 1 なのは n=0n = 0 の項だけが残るからである。□\square

定義 6.8 周期 2π2\pi の関数 ff が区分的に滑らか (piecewise smooth) であるとは、分割 −π=x0<x1<⋯<xm=π-\pi = x_0 < x_1 < \cdots < x_m = \pi があって、各開区間 (xj−1,xj)(x_{j-1}, x_j) で ff は C1C^1 級、かつ ff と f′f' が各端点で有限の片側極限をもつことをいう。

定理 6.9(ディリクレの定理)ff が区分的に滑らかなら、各点 xx で

lim⁡N→∞SNf(x)=f(x+0)+f(x−0)2\lim_{N \to \infty}S_Nf(x) = \frac{f(x + 0) + f(x - 0)}{2}

が成り立つ。特に ff が xx で連続なら SNf(x)→f(x)S_Nf(x) \to f(x)。

証明. DND_N は偶関数なので SNf(x)=12π∫0π(f(x+t)+f(x−t))DN(t) dtS_Nf(x) = \frac{1}{2\pi}\int_0^\pi\bigl(f(x + t) + f(x - t)\bigr)D_N(t)\ dt、また 12π∫0πDN=12\frac{1}{2\pi}\int_0^\pi D_N = \frac{1}{2}。よって

SNf(x)−f(x+0)+f(x−0)2=12π∫0πg(t)sin⁡(N+12)t dt,g(t)=f(x+t)−f(x+0)+f(x−t)−f(x−0)sin⁡(t/2)S_Nf(x) - \frac{f(x + 0) + f(x - 0)}{2} = \frac{1}{2\pi}\int_0^\pi g(t)\sin\left(N + \frac{1}{2}\right)t\,dt, \qquad g(t) = \frac{f(x + t) - f(x + 0) + f(x - t) - f(x - 0)}{\sin(t/2)}

ff の xx の右側の C1C^1 級の部分に平均値の定理を使うと、f′f' の片側極限が有限なので、小さい t>0t > 0 で ∣f(x+t)−f(x+0)∣≤Ct\lvert f(x + t) - f(x + 0) \rvert \leq Ct。左側も同様で、0<t≤π0 < t \leq \pi では sin⁡(t/2)≥t/π\sin(t/2) \geq t/\pi なので gg は (0,π](0, \pi] 上有界であり、区分的に連続なので [0,π][0, \pi] 上リーマン積分可能である。リーマン–ルベーグの補題(sin⁡λt=(eiλt−e−iλt)/(2i)\sin\lambda t = (e^{i\lambda t} - e^{-i\lambda t})/(2i) に適用する)により右辺は 0 に収束する。□\square

収束が xx の近くの ff の値だけで決まることに注意する(局所化原理)。

例 6.10 のこぎり波の級数で x=π/2x = \pi/2 とおくと、π2=2(1−13+15−⋯ )\frac{\pi}{2} = 2\left(1 - \frac{1}{3} + \frac{1}{5} - \cdots\right)、すなわちライプニッツの級数 ∑k≥0(−1)k2k+1=π4\sum_{k \geq 0}\frac{(-1)^k}{2k + 1} = \frac{\pi}{4} を得る。x=πx = \pi ではのこぎり波は π\pi から −π-\pi に跳び、級数の値は sin⁡nπ=0\sin n\pi = 0 より 0、すなわち跳びの中点である。

6.4 ギブス現象

不連続点の近くで部分和はどう振る舞うか。方形波 f=sgn⁡xf = \operatorname{sgn}x の部分和 S(x)=4π∑k=0M−1sin⁡(2k+1)x2k+1S(x) = \frac{4}{\pi}\sum_{k=0}^{M-1}\frac{\sin(2k + 1)x}{2k + 1} を考える。2sin⁡x∑k=0M−1cos⁡(2k+1)x=sin⁡2Mx2\sin x\sum_{k=0}^{M-1}\cos(2k + 1)x = \sin 2Mx より S′(x)=2sin⁡2Mxπsin⁡xS'(x) = \frac{2\sin 2Mx}{\pi\sin x} なので、SS は xM=π2Mx_M = \frac{\pi}{2M} で最初の極大をとる。uk=(2k+1)π2Mu_k = (2k + 1)\frac{\pi}{2M} とおくと

S(xM)=2π∑k=0M−1sin⁡ukuk⋅πM⟶2π∫0πsin⁡uu du≈1.1790(M→∞)S(x_M) = \frac{2}{\pi}\sum_{k=0}^{M-1}\frac{\sin u_k}{u_k}\cdot\frac{\pi}{M} \longrightarrow \frac{2}{\pi}\int_0^\pi\frac{\sin u}{u}\,du \approx 1.1790 \qquad (M \to \infty)

(区間 [0,π][0, \pi] を MM 等分した中点でのリーマン和)。部分和は跳びの近くで、跳びの大きさ 2 のおよそ 9% だけ行き過ぎ、項を増やしてもこの行き過ぎは消えない。行き過ぎる位置 xMx_M が不連続点に近づくだけである。これをギブス現象 (Gibbs phenomenon) という。特に、フーリエ級数は不連続点の近くで一様収束しない。

実験フーリエ級数の部分和項を増やすと方形波に近づく一方、不連続点の近くに残る「ギブス現象」を観察します。

この実験は JavaScript を有効にすると動きます。

6.5 フェイエールの定理

ディリクレ核は正負に大きく振動し、12π∫∣DN∣\frac{1}{2\pi}\int\lvert D_N \rvert は N→∞N \to \infty で発散する(問題 6.4)。これが部分和の収束を難しくしている。そこで部分和の算術平均(チェザロ平均)

σNf=1N+1∑k=0NSkf=f∗FN,FN(t)=1N+1∑k=0NDk(t)\sigma_Nf = \frac{1}{N + 1}\sum_{k=0}^N S_kf = f \ast F_N, \qquad F_N(t) = \frac{1}{N + 1}\sum_{k=0}^N D_k(t)

を考える(f∗FN(x)=12π∫−ππf(x−t)FN(t) dtf \ast F_N(x) = \frac{1}{2\pi}\int_{-\pi}^\pi f(x - t)F_N(t)\ dt)。FNF_N をフェイエール核 (Fejér kernel) という。

補題 6.11 FN(t)=1N+1(sin⁡(N+1)t2sin⁡t2)2F_N(t) = \dfrac{1}{N + 1}\left(\dfrac{\sin\frac{(N + 1)t}{2}}{\sin\frac{t}{2}}\right)^2。したがって (i) FN≥0F_N \geq 0、(ii) 12π∫−ππFN=1\frac{1}{2\pi}\int_{-\pi}^\pi F_N = 1、(iii) 0<δ≤∣t∣≤π0 < \delta \leq \lvert t \rvert \leq \pi で FN(t)≤1(N+1)sin⁡2(δ/2)→0F_N(t) \leq \frac{1}{(N + 1)\sin^2(\delta/2)} \to 0。

証明. 2sin⁡t2sin⁡(k+12)t=cos⁡kt−cos⁡(k+1)t2\sin\frac{t}{2}\sin\left(k + \frac{1}{2}\right)t = \cos kt - \cos(k + 1)t を k=0,…,Nk = 0, \dots, N で足すと 1−cos⁡(N+1)t=2sin⁡2(N+1)t21 - \cos(N + 1)t = 2\sin^2\frac{(N + 1)t}{2}。両辺を 2sin⁡2t22\sin^2\frac{t}{2} で割れば ∑kDk(t)\sum_k D_k(t) の式を得る。(ii) は各 DkD_k の積分が 1 だから。(iii) は sin⁡2(t/2)≥sin⁡2(δ/2)\sin^2(t/2) \geq \sin^2(\delta/2) から。□\square

(i)〜(iii) をみたす核の列を近似単位元 (approximate identity, good kernel) という。質量 1 の非負の山が原点に集中していく。

定理 6.12(フェイエールの定理, Fejér's theorem)ff が連続なら、σNf\sigma_Nf は ff に一様収束する。

証明. (ii) より σNf(x)−f(x)=12π∫−ππ(f(x−t)−f(x))FN(t) dt\sigma_Nf(x) - f(x) = \frac{1}{2\pi}\int_{-\pi}^\pi\bigl(f(x - t) - f(x)\bigr)F_N(t)\ dt。ε>0\varepsilon > 0 に対し、一様連続性から ∣t∣<δ\lvert t \rvert < \delta なら ∣f(x−t)−f(x)∣<ε\lvert f(x - t) - f(x) \rvert < \varepsilon(xx によらない)となる δ\delta をとる。∣t∣<δ\lvert t \rvert < \delta の部分は (i)(ii) より ε\varepsilon 以下、δ≤∣t∣≤π\delta \leq \lvert t \rvert \leq \pi の部分は (iii) より 2sup⁡∣f∣/((N+1)sin⁡2(δ/2))2\sup\lvert f \rvert/\bigl((N + 1)\sin^2(\delta/2)\bigr) 以下である。よって lim sup⁡Nsup⁡x∣σNf(x)−f(x)∣≤ε\limsup_N\sup_x\lvert \sigma_Nf(x) - f(x) \rvert \leq \varepsilon。□\square

系 6.13

  1. (一意性)連続関数 ff のフーリエ係数がすべて 0 なら f=0f = 0。
  2. 三角多項式全体は、周期 2π2\pi の連続関数全体の中で一様ノルムについて稠密である。
  3. 連続関数 ff について、SNf(x)S_Nf(x) がある点 xx で収束すれば、その極限は f(x)f(x) である。

証明. 1:σNf=0\sigma_Nf = 0 で、f=lim⁡σNf=0f = \lim\sigma_Nf = 0。2:σNf\sigma_Nf は三角多項式である。3:収束する数列の算術平均は同じ極限に収束するので、lim⁡SNf(x)=lim⁡σNf(x)=f(x)\lim S_Nf(x) = \lim\sigma_Nf(x) = f(x)。□\square

系 6.14(ワイエルシュトラスの近似定理)[a,b][a, b] 上の連続関数は、多項式で一様に近似できる。

証明. 1 次変換で [a,b]=[0,π][a, b] = [0, \pi] としてよい。ff を偶関数として [−π,π][-\pi, \pi] に拡張し周期 2π2\pi に延長すると連続関数 gg を得る。系 6.13 の 2 より sup⁡∣g−P∣<ε/2\sup\lvert g - P \rvert < \varepsilon/2 となる三角多項式 PP がある。einxe^{inx} のテイラー級数は [−π,π][-\pi, \pi] 上で一様収束するので、PP のテイラー展開を十分先で打ち切った多項式 QQ は sup⁡[−π,π]∣P−Q∣<ε/2\sup_{[-\pi, \pi]}\lvert P - Q \rvert < \varepsilon/2 をみたす。実部をとれば、[0,π][0, \pi] 上で ∣f−Re⁡Q∣<ε\lvert f - \operatorname{Re}Q \rvert < \varepsilon。□\square

6.6 L2L^2 理論:完全性とパーセヴァルの等式

定理 6.15(パーセヴァルの等式, Parseval's identity)リーマン積分可能な ff について ∥f−SNf∥2→0\lVert f - S_Nf \rVert_2 \to 0 であり、

12π∫−ππ∣f(x)∣2dx=∑n∈Z∣f^(n)∣2\frac{1}{2\pi}\int_{-\pi}^\pi\lvert f(x) \rvert^2dx = \sum_{n \in \mathbb{Z}}\lvert \hat{f}(n) \rvert^2

が成り立つ。より一般に ⟨f,g⟩=∑nf^(n)g^(n)‾\langle f, g \rangle = \sum_n\hat{f}(n)\overline{\hat{g}(n)}。

証明. まず ff が連続のとき:フェイエールの定理より、ε>0\varepsilon > 0 に対し sup⁡∣f−σMf∣<ε\sup\lvert f - \sigma_Mf \rvert < \varepsilon となる MM がある。N≥MN \geq M なら σMf∈TN\sigma_Mf \in \mathcal{T}_N なので、命題 6.4 より ∥f−SNf∥2≤∥f−σMf∥2≤sup⁡∣f−σMf∣<ε\lVert f - S_Nf \rVert_2 \leq \lVert f - \sigma_Mf \rVert_2 \leq \sup\lvert f - \sigma_Mf \rvert < \varepsilon。

一般の ff(∣f∣≤B\lvert f \rvert \leq B)では、定理 6.6 の証明の階段関数 ss(∣s∣≤B\lvert s \rvert \leq B, ∫∣f−s∣<ε\int\lvert f - s \rvert < \varepsilon)の跳びを短い線分でつなぎ、両端の値をそろえると、周期的な連続関数 hh で ∣h∣≤B\lvert h \rvert \leq B, ∫−ππ∣f−h∣<2ε\int_{-\pi}^\pi\lvert f - h \rvert < 2\varepsilon となるものが得られる。すると ∥f−h∥22≤2B2π∫∣f−h∣<2Bεπ\lVert f - h \rVert_2^2 \leq \frac{2B}{2\pi}\int\lvert f - h \rvert < \frac{2B\varepsilon}{\pi}。命題 6.4 より ∥f−SNf∥2≤∥f−SNh∥2≤∥f−h∥2+∥h−SNh∥2\lVert f - S_Nf \rVert_2 \leq \lVert f - S_Nh \rVert_2 \leq \lVert f - h \rVert_2 + \lVert h - S_Nh \rVert_2 で、第 2 項は前半より 0 に収束する。よって ∥f−SNf∥2→0\lVert f - S_Nf \rVert_2 \to 0 で、命題 6.4 の等式からパーセヴァルの等式を得る。一般形は ∣⋅∣2\lvert \cdot \rvert^2 の等式を f+gf + g, f+igf + ig などに適用すれば(偏極恒等式)得られる。□\square

{en}\lbrace e_n \rbrace は「完全」である:どの関数も L2L^2 の意味で展開でき、係数の 2 乗和がノルムの 2 乗に等しい。これは有限次元の正規直交基底の性質そのものである。

例 6.16(バーゼル問題)f(x)=xf(x) = x に適用すると、例 6.2 より ∣f^(±n)∣2=bn2/4=1/n2\lvert \hat{f}(\pm n) \rvert^2 = b_n^2/4 = 1/n^2 なので

12π∫−ππx2 dx=π23=2∑n=1∞1n2,よって∑n=1∞1n2=π26\frac{1}{2\pi}\int_{-\pi}^\pi x^2\,dx = \frac{\pi^2}{3} = 2\sum_{n=1}^\infty\frac{1}{n^2}, \qquad \text{よって} \quad \sum_{n=1}^\infty\frac{1}{n^2} = \frac{\pi^2}{6}

注意 6.17 逆に、∑∣cn∣2<∞\sum\lvert c_n \rvert^2 < \infty となる任意の数列 (cn)(c_n) に対し、フーリエ係数が cnc_n である関数が存在するか。リーマン積分可能な関数の範囲では答えは否定的だが、ルベーグ積分を用いた L2(−π,π)L^2(-\pi, \pi) の中では肯定的であり(リース–フィッシャーの定理)、L2L^2 は数列空間 ℓ2(Z)\ell^2(\mathbb{Z}) と同型なヒルベルト空間になる(測度と積分 第4章、関数解析 第2章)。

6.7 一様収束と項別微分

定理 6.18 ff が連続かつ区分的に滑らかなら、∑n∣f^(n)∣<∞\sum_n\lvert \hat{f}(n) \rvert < \infty であり、SNfS_Nf は ff に一様収束する。

証明. f′f'(有限個の点を除いて定義される)はリーマン積分可能である。各小区間で部分積分すると、ff の連続性と周期性から境界項は打ち消し合い、f′^(n)=inf^(n)\widehat{f'}(n) = in\hat{f}(n) を得る。コーシー–シュワルツの不等式とベッセルの不等式(f′f' に適用)より

∑n≠0∣f^(n)∣=∑n≠0∣f′^(n)∣∣n∣≤(∑n≠01n2)1/2(∑n∣f′^(n)∣2)1/2<∞\sum_{n \neq 0}\lvert \hat{f}(n) \rvert = \sum_{n \neq 0}\frac{\lvert \widehat{f'}(n) \rvert}{\lvert n \rvert} \leq \left(\sum_{n \neq 0}\frac{1}{n^2}\right)^{1/2}\left(\sum_n\lvert \widehat{f'}(n) \rvert^2\right)^{1/2} < \infty

よってワイエルシュトラスの M 判定法によりフーリエ級数は一様収束し、その和 gg は連続である。一様収束なので項別に積分でき、g^(n)=f^(n)\hat{g}(n) = \hat{f}(n)。系 6.13 の 1 より g=fg = f。□\square

例 6.19 f(x)=x2f(x) = x^2(∣x∣≤π\lvert x \rvert \leq \pi)の周期拡張は連続で区分的に滑らかなので、例 6.2 の級数はすべての xx で ff に収束する。x=πx = \pi とおくと π2=π23+4∑1n2\pi^2 = \frac{\pi^2}{3} + 4\sum\frac{1}{n^2} となり、再び ∑1n2=π26\sum\frac{1}{n^2} = \frac{\pi^2}{6}。x=0x = 0 とおけば ∑n≥1(−1)n+1n2=π212\sum_{n \geq 1}\frac{(-1)^{n+1}}{n^2} = \frac{\pi^2}{12}。

滑らかさとフーリエ係数の減衰は、次のように対応する。

命題 6.20 (1) ff が CkC^k 級(周期 2π2\pi)なら、f(k)^(n)=(in)kf^(n)\widehat{f^{(k)}}(n) = (in)^k\hat{f}(n) であり、∣n∣k∣f^(n)∣→0\lvert n \rvert^k\lvert \hat{f}(n) \rvert \to 0 かつ ∑∣n∣2k∣f^(n)∣2<∞\sum\lvert n \rvert^{2k}\lvert \hat{f}(n) \rvert^2 < \infty。(2) 連続関数 ff について ∑∣n∣k∣f^(n)∣<∞\sum\lvert n \rvert^k\lvert \hat{f}(n) \rvert < \infty なら、ff は CkC^k 級で、フーリエ級数を kk 回まで項別微分できる。

証明. (1) 部分積分を kk 回行い、f(k)f^{(k)} にリーマン–ルベーグの補題とベッセルの不等式を適用する。(2) j≤kj \leq k 回項別微分した級数は M 判定法で一様収束するので、和は CkC^k 級で項別微分できる(微分積分学 第6章)。和は ff とフーリエ係数が一致する連続関数なので、系 6.13 の 1 より ff に等しい。□\square

注意

項別微分はいつでも許されるわけではない。のこぎり波の級数 2∑(−1)n+1nsin⁡nx2\sum\frac{(-1)^{n+1}}{n}\sin nx を項別微分した 2∑(−1)n+1cos⁡nx2\sum(-1)^{n+1}\cos nx は、一般項が 0 に収束しないのでどの点でも収束しない。周期拡張したのこぎり波は x=±πx = \pm\pi で跳ぶので、定理 6.18 の仮定をみたさないのである。一方、項別積分はリーマン積分可能な関数についてつねに許される。実際、コーシー–シュワルツの不等式と定理 6.15 より、∣x∣≤π\lvert x \rvert \leq \pi で一様に ∣∫0x(f−SNf) dt∣≤∫−ππ∣f−SNf∣ dt≤2π∥f−SNf∥2→0\left\lvert \int_0^x (f - S_Nf)\ dt \right\rvert \leq \int_{-\pi}^\pi\lvert f - S_Nf \rvert\ dt \leq 2\pi\lVert f - S_Nf \rVert_2 \to 0 である。

6.8 応用:等周不等式

長さ LL の閉曲線で囲める面積 AA は最大でいくらか。答えは円のときで A=L2/(4π)A = L^2/(4\pi) である。

定理 6.21(等周不等式, isoperimetric inequality)C1C^1 級の単純閉曲線(速さが 0 にならないもの)の長さを LL、囲む領域の面積を AA とすると 4πA≤L24\pi A \leq L^2 であり、等号は円のときに限る。

証明.(フルヴィッツ)弧長に比例するパラメータで、反時計回りに z(s)=x(s)+iy(s)z(s) = x(s) + iy(s)(0≤s≤2π0 \leq s \leq 2\pi)と表すと、∣z′(s)∣=L/(2π)\lvert z'(s) \rvert = L/(2\pi) である。zz は周期 2π2\pi の C1C^1 級関数なので、z=∑cneinsz = \sum c_ne^{ins}、z′∼∑incneinsz' \sim \sum inc_ne^{ins}。パーセヴァルの等式より

L24π2=12π∫02π∣z′∣2ds=∑nn2∣cn∣2\frac{L^2}{4\pi^2} = \frac{1}{2\pi}\int_0^{2\pi}\lvert z' \rvert^2ds = \sum_n n^2\lvert c_n \rvert^2

グリーンの定理より A=12∫02π(xy′−yx′) ds=12Im⁡∫02πz‾z′ dsA = \frac{1}{2}\int_0^{2\pi}(xy' - yx')\ ds = \frac{1}{2}\operatorname{Im}\int_0^{2\pi}\overline{z}z'\ ds であり、パーセヴァルの等式の一般形より 12π∫02πz′z‾ ds=∑nincncn‾\frac{1}{2\pi}\int_0^{2\pi}z'\overline{z}\ ds = \sum_n inc_n\overline{c_n} なので A=π∑nn∣cn∣2A = \pi\sum_n n\lvert c_n \rvert^2。よって

L2−4πA=4π2∑n(n2−n)∣cn∣2≥0L^2 - 4\pi A = 4\pi^2\sum_n(n^2 - n)\lvert c_n \rvert^2 \geq 0

(整数 nn について n2−n≥0n^2 - n \geq 0)。等号は n≠0,1n \neq 0, 1 で cn=0c_n = 0、すなわち z=c0+c1eisz = c_0 + c_1e^{is}(円)のときに限る。□\square

6.9 発散する例

フェイエールの定理と定理 6.18 は、連続関数のフーリエ級数が各点で収束することまでは主張していない。実際、次が知られている(主張のみ)。

  • 連続関数でフーリエ級数がある点で発散するものが存在する(デュ・ボア–レイモン, 1873 年)。SNf(0)=12π∫fDNS_Nf(0) = \frac{1}{2\pi}\int fD_N は連続関数の空間上の有界線形汎関数で、そのノルムは 12π∫∣DN∣\frac{1}{2\pi}\int\lvert D_N \rvert(ルベーグ定数)に等しく、N→∞N \to \infty で対数的に発散する(問題 6.4)。関数解析の一様有界性原理(バナッハ–シュタインハウスの定理)から、sup⁡N∣SNf(0)∣=∞\sup_N\lvert S_Nf(0) \rvert = \infty となる連続関数 ff の存在が従う(関数解析 第3章)。
  • ルベーグ可積分な関数で、フーリエ級数がいたるところ発散するものが存在する(コルモゴロフ)。
  • 一方、L2L^2 に属する関数のフーリエ級数はほとんどいたるところ収束する(カールソン, 1966 年)。これは非常に深い定理である。

各点収束は微妙だが、チェザロ平均による収束(フェイエール)と L2L^2 収束(パーセヴァル)は非常に扱いやすい。現代の解析学がフーリエ級数を L2L^2 の枠組みで扱う理由の一つである。

まとめ

  • フーリエ係数は正規直交系 einxe^{inx} に関する内積であり、部分和 SNfS_Nf は NN 次三角多項式の中での L2L^2 最良近似である。ベッセルの不等式が成り立ち、係数は 0 に収束する(リーマン–ルベーグ)。
  • 部分和はディリクレ核との畳み込みで、区分的に滑らかな関数では各点で左右の極限の平均に収束する。不連続点の近くでは約 9% の行き過ぎが残る(ギブス現象)。
  • チェザロ平均はフェイエール核(近似単位元)との畳み込みで、連続関数に一様収束する。ここから一意性と三角多項式の稠密性、ワイエルシュトラスの近似定理が従う。
  • パーセヴァルの等式 12π∫∣f∣2=∑∣f^(n)∣2\frac{1}{2\pi}\int\lvert f \rvert^2 = \sum\lvert \hat{f}(n) \rvert^2 が成り立ち、三角関数系は完全である。バーゼル問題 ∑1/n2=π2/6\sum 1/n^2 = \pi^2/6 はその帰結である。
  • 連続で区分的に滑らかな関数のフーリエ級数は絶対かつ一様に収束する。滑らかなほど係数は速く減衰する。
  • パーセヴァルの等式から等周不等式 4πA≤L24\pi A \leq L^2 が得られる。
  • 連続関数のフーリエ級数が発散することもある。

演習問題

問題 6.1 ★ f(x)=∣sin⁡x∣f(x) = \lvert \sin x \rvert のフーリエ級数を求め、∑k=1∞14k2−1=12\sum_{k=1}^\infty\frac{1}{4k^2 - 1} = \frac{1}{2} を導け。

解答

偶関数なので bn=0b_n = 0。an=2π∫0πsin⁡xcos⁡nx dx=1π∫0π(sin⁡(1+n)x+sin⁡(1−n)x)dxa_n = \frac{2}{\pi}\int_0^\pi\sin x\cos nx\ dx = \frac{1}{\pi}\int_0^\pi\bigl(\sin(1 + n)x + \sin(1 - n)x\bigr)dx。nn が偶数なら 1±n1 \pm n は奇数で an=1π(21+n+21−n)=4π(1−n2)a_n = \frac{1}{\pi}\left(\frac{2}{1 + n} + \frac{2}{1 - n}\right) = \frac{4}{\pi(1 - n^2)}、nn が奇数なら an=0a_n = 0(n=1n = 1 も直接確かめられる)。よって ∣sin⁡x∣=2π−4π∑k≥1cos⁡2kx4k2−1\lvert \sin x \rvert = \frac{2}{\pi} - \frac{4}{\pi}\sum_{k \geq 1}\frac{\cos 2kx}{4k^2 - 1}(連続で区分的に滑らかなので定理 6.18 により等号)。x=0x = 0 とおけば主張を得る。

問題 6.2 ★★ f(x)=exf(x) = e^x(−π<x<π-\pi < x < \pi)を周期 2π2\pi に拡張した関数のフーリエ係数を求め、ディリクレの定理を x=πx = \pi で用いて ∑n∈Z11+n2=πcoth⁡π\sum_{n \in \mathbb{Z}}\frac{1}{1 + n^2} = \pi\coth\pi を示せ。

解答

f^(n)=12π∫−ππe(1−in)xdx=e(1−in)π−e−(1−in)π2π(1−in)=(−1)nsinh⁡ππ(1−in)\hat{f}(n) = \frac{1}{2\pi}\int_{-\pi}^\pi e^{(1 - in)x}dx = \frac{e^{(1 - in)\pi} - e^{-(1 - in)\pi}}{2\pi(1 - in)} = \frac{(-1)^n\sinh\pi}{\pi(1 - in)}。x=πx = \pi で左極限は eπe^\pi、右極限は e−πe^{-\pi} なので、定理 6.9 より SNf(π)→cosh⁡πS_Nf(\pi) \to \cosh\pi。一方 SNf(π)=∑∣n∣≤Nf^(n)(−1)n=sinh⁡ππ∑∣n∣≤N1+in1+n2=sinh⁡ππ∑∣n∣≤N11+n2S_Nf(\pi) = \sum_{\lvert n \rvert \leq N}\hat{f}(n)(-1)^n = \frac{\sinh\pi}{\pi}\sum_{\lvert n \rvert \leq N}\frac{1 + in}{1 + n^2} = \frac{\sinh\pi}{\pi}\sum_{\lvert n \rvert \leq N}\frac{1}{1 + n^2}(虚部は nn と −n-n で打ち消し合う)。よって ∑11+n2=πcosh⁡π/sinh⁡π\sum\frac{1}{1 + n^2} = \pi\cosh\pi/\sinh\pi。

問題 6.3 ★★ f(x)=x2f(x) = x^2 にパーセヴァルの等式を適用して ∑n=1∞1n4=π490\sum_{n=1}^\infty\frac{1}{n^4} = \frac{\pi^4}{90} を示せ。

解答

f^(0)=π23\hat{f}(0) = \frac{\pi^2}{3}、n≠0n \neq 0 で f^(n)=an2=2(−1)nn2\hat{f}(n) = \frac{a_n}{2} = \frac{2(-1)^n}{n^2}。12π∫−ππx4dx=π45\frac{1}{2\pi}\int_{-\pi}^\pi x^4dx = \frac{\pi^4}{5} なので π45=π49+2∑n≥14n4\frac{\pi^4}{5} = \frac{\pi^4}{9} + 2\sum_{n \geq 1}\frac{4}{n^4}、よって ∑1n4=18(π45−π49)=π490\sum\frac{1}{n^4} = \frac{1}{8}\left(\frac{\pi^4}{5} - \frac{\pi^4}{9}\right) = \frac{\pi^4}{90}。

問題 6.4 ★★(ルベーグ定数)12π∫−ππ∣DN(t)∣ dt≥4π2∑k=1N1k\frac{1}{2\pi}\int_{-\pi}^\pi\lvert D_N(t) \rvert\ dt \geq \frac{4}{\pi^2}\sum_{k=1}^N\frac{1}{k} を示せ。特に左辺は N→∞N \to \infty で発散する。

解答

∣sin⁡(t/2)∣≤∣t∣/2\lvert \sin(t/2) \rvert \leq \lvert t \rvert/2 より ∣DN(t)∣≥2∣sin⁡(N+12)t∣∣t∣\lvert D_N(t) \rvert \geq \frac{2\lvert \sin(N + \frac{1}{2})t \rvert}{\lvert t \rvert}。偶関数なので

12π∫−ππ∣DN∣≥2π∫0π∣sin⁡(N+12)t∣tdt=2π∫0(N+12)π∣sin⁡u∣udu≥2π∑k=1N1kπ∫(k−1)πkπ∣sin⁡u∣ du=4π2∑k=1N1k\frac{1}{2\pi}\int_{-\pi}^\pi\lvert D_N \rvert \geq \frac{2}{\pi}\int_0^\pi\frac{\lvert \sin(N + \frac{1}{2})t \rvert}{t}dt = \frac{2}{\pi}\int_0^{(N + \frac{1}{2})\pi}\frac{\lvert \sin u \rvert}{u}du \geq \frac{2}{\pi}\sum_{k=1}^N\frac{1}{k\pi}\int_{(k-1)\pi}^{k\pi}\lvert \sin u \rvert\,du = \frac{4}{\pi^2}\sum_{k=1}^N\frac{1}{k}

問題 6.5 ★★(ヴィルティンガーの不等式)ff が周期 2π2\pi の C1C^1 級関数で ∫02πf dx=0\int_0^{2\pi}f\ dx = 0 なら ∫02π∣f∣2dx≤∫02π∣f′∣2dx\int_0^{2\pi}\lvert f \rvert^2dx \leq \int_0^{2\pi}\lvert f' \rvert^2dx であり、等号は f=acos⁡x+bsin⁡xf = a\cos x + b\sin x のときに限ることを示せ。

解答

f^(0)=0\hat{f}(0) = 0 と f′^(n)=inf^(n)\widehat{f'}(n) = in\hat{f}(n) およびパーセヴァルの等式より 12π∫∣f′∣2=∑n≠0n2∣f^(n)∣2≥∑n≠0∣f^(n)∣2=12π∫∣f∣2\frac{1}{2\pi}\int\lvert f' \rvert^2 = \sum_{n \neq 0}n^2\lvert \hat{f}(n) \rvert^2 \geq \sum_{n \neq 0}\lvert \hat{f}(n) \rvert^2 = \frac{1}{2\pi}\int\lvert f \rvert^2。等号は ∣n∣≥2\lvert n \rvert \geq 2 で f^(n)=0\hat{f}(n) = 0、すなわち f=c1eix+c−1e−ixf = c_1e^{ix} + c_{-1}e^{-ix} のときに限る(定理 6.18 より級数は ff に収束する)。これは第5章の最小固有値の評価の周期境界条件版である。

問題 6.6 ★★ リーマン–ルベーグの補題を用いて lim⁡N→∞∫0(N+12)πsin⁡uudu=π2\lim_{N \to \infty}\int_0^{(N + \frac{1}{2})\pi}\frac{\sin u}{u}du = \frac{\pi}{2} を示せ。広義積分 ∫0∞sin⁡uudu\int_0^\infty\frac{\sin u}{u}du が収束することを認めて、その値を求めよ。

解答

12π∫−ππDN=1\frac{1}{2\pi}\int_{-\pi}^\pi D_N = 1 と偶関数性から ∫0πsin⁡(N+12)tsin⁡(t/2)dt=π\int_0^\pi\frac{\sin(N + \frac{1}{2})t}{\sin(t/2)}dt = \pi。h(t)=1sin⁡(t/2)−2th(t) = \frac{1}{\sin(t/2)} - \frac{2}{t} は (0,π](0, \pi] で連続で、t→0t \to 0 で h(t)=t−2sin⁡(t/2)tsin⁡(t/2)→0h(t) = \frac{t - 2\sin(t/2)}{t\sin(t/2)} \to 0 なので [0,π][0, \pi] 上の連続関数に拡張できる。リーマン–ルベーグの補題より ∫0πh(t)sin⁡(N+12)t dt→0\int_0^\pi h(t)\sin(N + \frac{1}{2})t\ dt \to 0、よって ∫0π2sin⁡(N+12)ttdt→π\int_0^\pi\frac{2\sin(N + \frac{1}{2})t}{t}dt \to \pi。u=(N+12)tu = (N + \frac{1}{2})t と置換すれば主張を得る。広義積分が収束すれば、その値は部分列の極限 π2\frac{\pi}{2} に等しい。

問題 6.7 ★★★(ワイルの一様分布定理)α\alpha を無理数とする。周期 2π2\pi の任意の連続関数 ff について lim⁡N→∞1N∑k=1Nf(2πkα)=12π∫02πf(x) dx\lim_{N \to \infty}\frac{1}{N}\sum_{k=1}^N f(2\pi k\alpha) = \frac{1}{2\pi}\int_0^{2\pi}f(x)\ dx を示せ。

解答

f=emf = e_m(m≠0m \neq 0)なら、e2πimα≠1e^{2\pi im\alpha} \neq 1(α\alpha は無理数)なので

∣1N∑k=1Ne2πimkα∣=1N∣e2πimα(e2πimNα−1)e2πimα−1∣≤2N∣e2πimα−1∣→0\left\lvert \frac{1}{N}\sum_{k=1}^N e^{2\pi imk\alpha} \right\rvert = \frac{1}{N}\left\lvert \frac{e^{2\pi im\alpha}(e^{2\pi imN\alpha} - 1)}{e^{2\pi im\alpha} - 1} \right\rvert \leq \frac{2}{N\lvert e^{2\pi im\alpha} - 1 \rvert} \to 0

で、右辺の積分も 0。m=0m = 0 なら両辺とも 1。線形性から三角多項式について成り立つ。一般の連続関数 ff と ε>0\varepsilon > 0 に対し、系 6.13 の 2 より sup⁡∣f−P∣<ε\sup\lvert f - P \rvert < \varepsilon となる三角多項式 PP をとると、左辺の平均と右辺の積分はそれぞれ PP の場合と ε\varepsilon 以内の差しかないので、lim sup⁡\limsup で 2ε2\varepsilon 以内となり主張が従う。区間の定義関数を連続関数で上下から近似すれば、kαk\alpha の小数部分が [0,1)[0, 1) に一様に分布することがわかる。

問題 6.8 ★ [0,π][0, \pi] 上の f(x)=x(π−x)f(x) = x(\pi - x) の正弦級数を求め、∑k=0∞(−1)k(2k+1)3=π332\sum_{k=0}^\infty\frac{(-1)^k}{(2k + 1)^3} = \frac{\pi^3}{32} を示せ。

解答

部分積分を 2 回行うと ∫0πx(π−x)sin⁡nx dx=2n2∫0πsin⁡nx dx=2(1−(−1)n)n3\int_0^\pi x(\pi - x)\sin nx\ dx = \frac{2}{n^2}\int_0^\pi\sin nx\ dx = \frac{2(1 - (-1)^n)}{n^3} なので bn=8πn3b_n = \frac{8}{\pi n^3}(nn 奇数)、0(nn 偶数)。奇関数拡張は連続で区分的に滑らかなので、定理 6.18 より x(π−x)=8π∑k≥0sin⁡(2k+1)x(2k+1)3x(\pi - x) = \frac{8}{\pi}\sum_{k \geq 0}\frac{\sin(2k + 1)x}{(2k + 1)^3}(0≤x≤π0 \leq x \leq \pi)。x=π/2x = \pi/2 とおくと π24=8π∑k(−1)k(2k+1)3\frac{\pi^2}{4} = \frac{8}{\pi}\sum_k\frac{(-1)^k}{(2k + 1)^3}。

この章を読み終えたら

「読了」にすると学習記録とロードマップに反映されます。演習の自己採点もお忘れなく。

この章の誤りを報告GitHub で見る