この章の目標
フーリエ係数を計算でき、ベッセルの不等式とリーマン–ルベーグの補題を証明できる
ディリクレ核を用いて、区分的に滑らかな関数のフーリエ級数の各点収束を証明できる
フェイエールの定理とパーセヴァルの等式を証明し、三角関数系の完全性を説明できる
滑らかさと係数の減衰の関係を理解し、一様収束・項別微分の条件を述べられる
バーゼル問題・等周不等式・ワイエルシュトラスの近似定理にフーリエ級数を応用できる
前提 :第5章 、微分積分学 第5章 、第6章
フーリエは 1807 年頃、棒の中の熱の伝わり方を研究し、「任意の」関数が三角関数の和 ∑ ( a n cos n x + b n sin n x ) \sum(a_n\cos nx + b_n\sin nx) ∑ ( a n cos n x + b n sin n x ) で表せると主張した。熱方程式を変数分離すると sin n x \sin nx sin n x という固有関数が現れる(第5章、第7章)ので、初期温度分布をこれらの和に展開できれば、熱方程式が解けるからである。この主張がどのような意味で正しいかを明らかにする試みから、関数の概念、リーマン積分、一様収束、集合論、ルベーグ積分が生まれた。本章では、フーリエ級数が「どの関数について」「どの意味で」元の関数を表すかを調べる。
以下、関数は R \mathbb{R} R 上の周期 2 π 2\pi 2 π の複素数値関数で、[ − π , π ] [-\pi, \pi] [ − π , π ] 上リーマン積分可能なものとする(このような関数は有界である)。内積とノルムを
⟨ f , g ⟩ = 1 2 π ∫ − π π f ( x ) g ( x ) ‾ d x , ∥ 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 , g ⟩ = 2 π 1 ∫ − π π f ( x ) g ( x ) d x , ∥ f ∥ 2 = ⟨ f , f ⟩ 1/2
とおく(∥ f ∥ 2 = 0 \lVert f \rVert_2 = 0 ∥ f ∥ 2 = 0 でも f = 0 f = 0 f = 0 とは限らないので、厳密には半ノルムである)。e n ( x ) = e i n x e_n(x) = e^{inx} e n ( x ) = e in x とすると ⟨ e m , e n ⟩ \langle e_m, e_n \rangle ⟨ e m , e n ⟩ は m = n m = n m = n なら 1、そうでなければ 0 である。
6.1 三角級数とフーリエ係数
定義 6.1 (フーリエ級数, Fourier series)n ∈ Z n \in \mathbb{Z} n ∈ Z に対し
f ^ ( n ) = ⟨ f , e n ⟩ = 1 2 π ∫ − π π f ( x ) e − i n x d x \hat{f}(n) = \langle f, e_n \rangle = \frac{1}{2\pi}\int_{-\pi}^\pi f(x)e^{-inx}\,dx f ^ ( n ) = ⟨ f , e n ⟩ = 2 π 1 ∫ − π π f ( x ) e − in x d x
を f f f のフーリエ係数 、形式的な級数 ∑ n ∈ Z f ^ ( n ) e i n x \sum_{n \in \mathbb{Z}}\hat{f}(n)e^{inx} ∑ n ∈ Z f ^ ( n ) e in x をフーリエ級数 といい、f ∼ ∑ f ^ ( n ) e i n x f \sim \sum\hat{f}(n)e^{inx} f ∼ ∑ f ^ ( n ) e in x と書く。部分和を S N f ( x ) = ∑ ∣ n ∣ ≤ N f ^ ( n ) e i n x S_Nf(x) = \sum_{\lvert n \rvert \leq N}\hat{f}(n)e^{inx} S N f ( x ) = ∑ ∣ n ∣ ≤ N f ^ ( n ) e in x とおく。
実数値関数では、a n = 1 π ∫ − π π f ( x ) cos n x d x a_n = \frac{1}{\pi}\int_{-\pi}^\pi f(x)\cos nx\ dx a n = π 1 ∫ − π π f ( x ) cos n x d x , b n = 1 π ∫ − π π f ( x ) sin n x d x b_n = \frac{1}{\pi}\int_{-\pi}^\pi f(x)\sin nx\ dx b n = π 1 ∫ − π π f ( x ) sin n x d x とおくと f ^ ( ± n ) = ( a n ∓ i b n ) / 2 \hat{f}(\pm n) = (a_n \mp ib_n)/2 f ^ ( ± n ) = ( a n ∓ i b n ) /2 で、S N f = a 0 2 + ∑ n = 1 N ( a n cos n x + b n sin n x ) S_Nf = \frac{a_0}{2} + \sum_{n=1}^N(a_n\cos nx + b_n\sin nx) S N f = 2 a 0 + ∑ n = 1 N ( a n cos n x + b n sin n x ) である。f f f が奇関数なら a n = 0 a_n = 0 a n = 0 、偶関数なら b n = 0 b_n = 0 b n = 0 。
例 6.2 周期 2 π 2\pi 2 π に拡張した次の関数を考える(部分積分で計算できる)。
のこぎり波 f ( x ) = x f(x) = x f ( x ) = x (− π < x < π -\pi < x < \pi − π < x < π ):b n = 2 π ∫ 0 π x sin n x d x = 2 ( − 1 ) n + 1 n b_n = \frac{2}{\pi}\int_0^\pi x\sin nx\ dx = \frac{2(-1)^{n+1}}{n} b n = π 2 ∫ 0 π x sin n x d x = n 2 ( − 1 ) n + 1 で、x ∼ 2 ∑ n ≥ 1 ( − 1 ) n + 1 n sin n x x \sim 2\sum_{n \geq 1}\frac{(-1)^{n+1}}{n}\sin nx x ∼ 2 ∑ n ≥ 1 n ( − 1 ) n + 1 sin n x 。
方形波 f ( x ) = sgn x f(x) = \operatorname{sgn}x f ( x ) = sgn x :b n = 2 π ∫ 0 π sin n x d x = 2 ( 1 − ( − 1 ) n ) π n b_n = \frac{2}{\pi}\int_0^\pi\sin nx\ dx = \frac{2(1 - (-1)^n)}{\pi n} b n = π 2 ∫ 0 π sin n x d x = π n 2 ( 1 − ( − 1 ) n ) で、sgn x ∼ 4 π ∑ k ≥ 0 sin ( 2 k + 1 ) x 2 k + 1 \operatorname{sgn}x \sim \frac{4}{\pi}\sum_{k \geq 0}\frac{\sin(2k + 1)x}{2k + 1} sgn x ∼ π 4 ∑ k ≥ 0 2 k + 1 s i n ( 2 k + 1 ) x 。
f ( x ) = x 2 f(x) = x^2 f ( x ) = x 2 (∣ x ∣ ≤ π \lvert x \rvert \leq \pi ∣ x ∣ ≤ π ):a 0 = 2 π 2 3 a_0 = \frac{2\pi^2}{3} a 0 = 3 2 π 2 、a n = 2 π ∫ 0 π x 2 cos n x d x = 4 ( − 1 ) n n 2 a_n = \frac{2}{\pi}\int_0^\pi x^2\cos nx\ dx = \frac{4(-1)^n}{n^2} a n = π 2 ∫ 0 π x 2 cos n x d x = n 2 4 ( − 1 ) n で、x 2 ∼ π 2 3 + 4 ∑ n ≥ 1 ( − 1 ) n n 2 cos n x x^2 \sim \frac{\pi^2}{3} + 4\sum_{n \geq 1}\frac{(-1)^n}{n^2}\cos nx x 2 ∼ 3 π 2 + 4 ∑ n ≥ 1 n 2 ( − 1 ) n cos n x 。
6.2 ベッセルの不等式とリーマン–ルベーグの補題
命題 6.4 (最良近似)T N \mathcal{T}_N T N を e n e_n e n (∣ n ∣ ≤ N \lvert n \rvert \leq N ∣ n ∣ ≤ N )の一次結合(N N N 次以下の三角多項式)全体とする。
∥ f − S N f ∥ 2 2 = ∥ f ∥ 2 2 − ∑ ∣ 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 ∥ f − S N f ∥ 2 2 = ∥ f ∥ 2 2 − ∣ n ∣ ≤ N ∑ ∣ f ^ ( n ) ∣ 2
であり、任意の P ∈ T N P \in \mathcal{T}_N P ∈ T N に対し ∥ f − S N f ∥ 2 ≤ ∥ f − P ∥ 2 \lVert f - S_Nf \rVert_2 \leq \lVert f - P \rVert_2 ∥ f − S N f ∥ 2 ≤ ∥ f − P ∥ 2 が成り立つ。
証明. ∣ n ∣ ≤ N \lvert n \rvert \leq N ∣ n ∣ ≤ N なら ⟨ f − S N f , e n ⟩ = f ^ ( n ) − f ^ ( n ) = 0 \langle f - S_Nf, e_n \rangle = \hat{f}(n) - \hat{f}(n) = 0 ⟨ f − S N f , e n ⟩ = f ^ ( n ) − f ^ ( n ) = 0 なので、f − S N f f - S_Nf f − S N f は T N \mathcal{T}_N T N と直交する。ピタゴラスの定理より ∥ f ∥ 2 2 = ∥ f − S N f ∥ 2 2 + ∥ S N f ∥ 2 2 \lVert f \rVert_2^2 = \lVert f - S_Nf \rVert_2^2 + \lVert S_Nf \rVert_2^2 ∥ f ∥ 2 2 = ∥ f − S N f ∥ 2 2 + ∥ S N f ∥ 2 2 で、∥ S N f ∥ 2 2 = ∑ ∣ n ∣ ≤ N ∣ f ^ ( n ) ∣ 2 \lVert S_Nf \rVert_2^2 = \sum_{\lvert n \rvert \leq N}\lvert \hat{f}(n) \rvert^2 ∥ S N f ∥ 2 2 = ∑ ∣ n ∣ ≤ N ∣ f ^ ( n ) ∣ 2 。また S N f − P ∈ T N S_Nf - P \in \mathcal{T}_N S N f − P ∈ T N なので ∥ f − P ∥ 2 2 = ∥ f − S N f ∥ 2 2 + ∥ S N f − P ∥ 2 2 \lVert f - P \rVert_2^2 = \lVert f - S_Nf \rVert_2^2 + \lVert S_Nf - P \rVert_2^2 ∥ f − P ∥ 2 2 = ∥ f − S N f ∥ 2 2 + ∥ S N f − P ∥ 2 2 。□ \square □
S N f S_Nf S N f は f f f の T N \mathcal{T}_N T N への直交射影であり、L 2 L^2 L 2 の意味での最良近似である。
系 6.5 (ベッセルの不等式, Bessel's inequality)∑ n ∈ Z ∣ f ^ ( n ) ∣ 2 ≤ 1 2 π ∫ − π π ∣ f ( x ) ∣ 2 d x \sum_{n \in \mathbb{Z}}\lvert \hat{f}(n) \rvert^2 \leq \frac{1}{2\pi}\int_{-\pi}^\pi\lvert f(x) \rvert^2dx ∑ n ∈ Z ∣ f ^ ( n ) ∣ 2 ≤ 2 π 1 ∫ − π π ∣ f ( x ) ∣ 2 d x 。特に f ^ ( n ) → 0 \hat{f}(n) \to 0 f ^ ( n ) → 0 (∣ n ∣ → ∞ \lvert n \rvert \to \infty ∣ n ∣ → ∞ )。
より一般に次が成り立つ。
定理 6.6 (リーマン–ルベーグの補題, Riemann–Lebesgue lemma)f f f が [ a , b ] [a, b] [ a , b ] 上リーマン積分可能なら、λ ∈ R \lambda \in \mathbb{R} λ ∈ R について
lim ∣ λ ∣ → ∞ ∫ a b f ( x ) e i λ x d x = 0 \lim_{\lvert \lambda \rvert \to \infty}\int_a^b f(x)e^{i\lambda x}\,dx = 0 ∣ λ ∣ → ∞ lim ∫ a b f ( x ) e iλ x d x = 0
証明. 区間の特性関数 1 [ c , d ] \mathbf{1}_{[c, d]} 1 [ c , d ] については ∣ ∫ c d e i λ x d x ∣ = ∣ e i λ d − e i λ c i λ ∣ ≤ 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} ∫ c d e iλ x d x = iλ e iλ d − e iλ c ≤ ∣ λ ∣ 2 なので、階段関数 s = ∑ j m j 1 I j s = \sum_j m_j\mathbf{1}_{I_j} s = ∑ j m j 1 I j については積分は 2 ∑ j ∣ m j ∣ / ∣ λ ∣ → 0 2\sum_j\lvert m_j \rvert/\lvert \lambda \rvert \to 0 2 ∑ j ∣ m j ∣ / ∣ λ ∣ → 0 である。ε > 0 \varepsilon > 0 ε > 0 に対し、上ダルブー和と下ダルブー和の差が ε \varepsilon ε より小さい分割をとり、各小区間で f f f の下限をとる階段関数を s s s とすると、∫ a b ∣ f − s ∣ d x < ε \int_a^b\lvert f - s \rvert\ dx < \varepsilon ∫ a b ∣ f − s ∣ d x < ε 。よって lim sup ∣ ∫ a b f e i λ x d x ∣ ≤ ε + lim sup ∣ ∫ a b s e i λ x d x ∣ = ε \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 lim sup ∫ a b f e iλ x d x ≤ ε + lim sup ∫ a b s e iλ x d x = ε 。□ \square □
速く振動する関数との積の積分は、正負が打ち消し合って 0 に近づく。
6.3 ディリクレ核と各点収束
補題 6.7 S N f ( x ) = 1 2 π ∫ − π π f ( x − t ) D N ( t ) d t S_Nf(x) = \frac{1}{2\pi}\int_{-\pi}^\pi f(x - t)D_N(t)\ dt S N f ( x ) = 2 π 1 ∫ − π π f ( x − t ) D N ( t ) d t 。ここでディリクレ核 (Dirichlet kernel) は
D N ( t ) = ∑ ∣ n ∣ ≤ N e i n t = sin ( N + 1 2 ) t sin t 2 D_N(t) = \sum_{\lvert n \rvert \leq N}e^{int} = \frac{\sin\left(N + \frac{1}{2}\right)t}{\sin\frac{t}{2}} D N ( t ) = ∣ n ∣ ≤ N ∑ e in t = sin 2 t sin ( N + 2 1 ) t
であり、偶関数で 1 2 π ∫ − π π D N ( t ) d t = 1 \frac{1}{2\pi}\int_{-\pi}^\pi D_N(t)\ dt = 1 2 π 1 ∫ − π π D N ( t ) d t = 1 をみたす。
証明. S N f ( x ) = ∑ ∣ n ∣ ≤ N 1 2 π ∫ f ( y ) e i n ( x − y ) d y = 1 2 π ∫ f ( y ) D N ( x − y ) d y S_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 S N f ( x ) = ∑ ∣ n ∣ ≤ N 2 π 1 ∫ f ( y ) e in ( x − y ) d y = 2 π 1 ∫ f ( y ) D N ( x − y ) d y で、t = x − y t = x - y t = x − y と置換し周期性を使う。等比数列の和から D N ( t ) = e i ( N + 1 ) t − e − i N t e i t − 1 D_N(t) = \frac{e^{i(N + 1)t} - e^{-iNt}}{e^{it} - 1} D N ( t ) = e i t − 1 e i ( N + 1 ) t − e − i N t で、分子分母に e − i t / 2 e^{-it/2} e − i t /2 を掛ければ右の形になる。積分が 1 なのは n = 0 n = 0 n = 0 の項だけが残るからである。□ \square □
定義 6.8 周期 2 π 2\pi 2 π の関数 f f f が区分的に滑らか (piecewise smooth) であるとは、分割 − π = x 0 < x 1 < ⋯ < x m = π -\pi = x_0 < x_1 < \cdots < x_m = \pi − π = x 0 < x 1 < ⋯ < x m = π があって、各開区間 ( x j − 1 , x j ) (x_{j-1}, x_j) ( x j − 1 , x j ) で f f f は C 1 C^1 C 1 級、かつ f f f と f ′ f' f ′ が各端点で有限の片側極限をもつことをいう。
定理 6.9 (ディリクレの定理)f f f が区分的に滑らかなら、各点 x x x で
lim N → ∞ S N f ( x ) = f ( x + 0 ) + f ( x − 0 ) 2 \lim_{N \to \infty}S_Nf(x) = \frac{f(x + 0) + f(x - 0)}{2} N → ∞ lim S N f ( x ) = 2 f ( x + 0 ) + f ( x − 0 )
が成り立つ。特に f f f が x x x で連続なら S N f ( x ) → f ( x ) S_Nf(x) \to f(x) S N f ( x ) → f ( x ) 。
証明. D N D_N D N は偶関数なので S N f ( x ) = 1 2 π ∫ 0 π ( f ( x + t ) + f ( x − t ) ) D N ( t ) d t S_Nf(x) = \frac{1}{2\pi}\int_0^\pi\bigl(f(x + t) + f(x - t)\bigr)D_N(t)\ dt S N f ( x ) = 2 π 1 ∫ 0 π ( f ( x + t ) + f ( x − t ) ) D N ( t ) d t 、また 1 2 π ∫ 0 π D N = 1 2 \frac{1}{2\pi}\int_0^\pi D_N = \frac{1}{2} 2 π 1 ∫ 0 π D N = 2 1 。よって
S N f ( x ) − f ( x + 0 ) + f ( x − 0 ) 2 = 1 2 π ∫ 0 π g ( t ) sin ( N + 1 2 ) t d t , 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)} S N f ( x ) − 2 f ( x + 0 ) + f ( x − 0 ) = 2 π 1 ∫ 0 π g ( t ) sin ( N + 2 1 ) t d t , g ( t ) = sin ( t /2 ) f ( x + t ) − f ( x + 0 ) + f ( x − t ) − f ( x − 0 )
f f f の x x x の右側の C 1 C^1 C 1 級の部分に平均値の定理を使うと、f ′ f' f ′ の片側極限が有限なので、小さい t > 0 t > 0 t > 0 で ∣ f ( x + t ) − f ( x + 0 ) ∣ ≤ C t \lvert f(x + t) - f(x + 0) \rvert \leq Ct ∣ f ( x + t ) − f ( x + 0 )∣ ≤ C t 。左側も同様で、0 < t ≤ π 0 < t \leq \pi 0 < t ≤ π では sin ( t / 2 ) ≥ t / π \sin(t/2) \geq t/\pi sin ( t /2 ) ≥ t / π なので g g g は ( 0 , π ] (0, \pi] ( 0 , π ] 上有界であり、区分的に連続なので [ 0 , π ] [0, \pi] [ 0 , π ] 上リーマン積分可能である。リーマン–ルベーグの補題(sin λ t = ( e i λ t − e − i λ t ) / ( 2 i ) \sin\lambda t = (e^{i\lambda t} - e^{-i\lambda t})/(2i) sin λ t = ( e iλ t − e − iλ t ) / ( 2 i ) に適用する)により右辺は 0 に収束する。□ \square □
収束が x x x の近くの f f f の値だけで決まることに注意する(局所化原理 )。
例 6.10 のこぎり波の級数で x = π / 2 x = \pi/2 x = π /2 とおくと、π 2 = 2 ( 1 − 1 3 + 1 5 − ⋯ ) \frac{\pi}{2} = 2\left(1 - \frac{1}{3} + \frac{1}{5} - \cdots\right) 2 π = 2 ( 1 − 3 1 + 5 1 − ⋯ ) 、すなわちライプニッツの級数 ∑ k ≥ 0 ( − 1 ) k 2 k + 1 = π 4 \sum_{k \geq 0}\frac{(-1)^k}{2k + 1} = \frac{\pi}{4} ∑ k ≥ 0 2 k + 1 ( − 1 ) k = 4 π を得る。x = π x = \pi x = π ではのこぎり波は π \pi π から − π -\pi − π に跳び、級数の値は sin n π = 0 \sin n\pi = 0 sin nπ = 0 より 0、すなわち跳びの中点である。
6.4 ギブス現象
不連続点の近くで部分和はどう振る舞うか。方形波 f = sgn x f = \operatorname{sgn}x f = sgn x の部分和 S ( x ) = 4 π ∑ k = 0 M − 1 sin ( 2 k + 1 ) x 2 k + 1 S(x) = \frac{4}{\pi}\sum_{k=0}^{M-1}\frac{\sin(2k + 1)x}{2k + 1} S ( x ) = π 4 ∑ k = 0 M − 1 2 k + 1 s i n ( 2 k + 1 ) x を考える。2 sin x ∑ k = 0 M − 1 cos ( 2 k + 1 ) x = sin 2 M x 2\sin x\sum_{k=0}^{M-1}\cos(2k + 1)x = \sin 2Mx 2 sin x ∑ k = 0 M − 1 cos ( 2 k + 1 ) x = sin 2 M x より S ′ ( x ) = 2 sin 2 M x π sin x S'(x) = \frac{2\sin 2Mx}{\pi\sin x} S ′ ( x ) = π s i n x 2 s i n 2 M x なので、S S S は x M = π 2 M x_M = \frac{\pi}{2M} x M = 2 M π で最初の極大をとる。u k = ( 2 k + 1 ) π 2 M u_k = (2k + 1)\frac{\pi}{2M} u k = ( 2 k + 1 ) 2 M π とおくと
S ( x M ) = 2 π ∑ k = 0 M − 1 sin u k u k ⋅ π M ⟶ 2 π ∫ 0 π sin u u d u ≈ 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) S ( x M ) = π 2 k = 0 ∑ M − 1 u k sin u k ⋅ M π ⟶ π 2 ∫ 0 π u sin u d u ≈ 1.1790 ( M → ∞ )
(区間 [ 0 , π ] [0, \pi] [ 0 , π ] を M M M 等分した中点でのリーマン和)。部分和は跳びの近くで、跳びの大きさ 2 のおよそ 9% だけ行き過ぎ、項を増やしてもこの行き過ぎは消えない。行き過ぎる位置 x M x_M x M が不連続点に近づくだけである。これをギブス現象 (Gibbs phenomenon) という。特に、フーリエ級数は不連続点の近くで一様収束しない。
実験フーリエ級数の部分和 項を増やすと方形波に近づく一方、不連続点の近くに残る「ギブス現象」を観察します。
この実験は JavaScript を有効にすると動きます。
6.5 フェイエールの定理
ディリクレ核は正負に大きく振動し、1 2 π ∫ ∣ D N ∣ \frac{1}{2\pi}\int\lvert D_N \rvert 2 π 1 ∫ ∣ D N ∣ は N → ∞ N \to \infty N → ∞ で発散する(問題 6.4)。これが部分和の収束を難しくしている。そこで部分和の算術平均(チェザロ平均 )
σ N f = 1 N + 1 ∑ k = 0 N S k f = f ∗ F N , F N ( t ) = 1 N + 1 ∑ k = 0 N D k ( 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) σ N f = N + 1 1 k = 0 ∑ N S k f = f ∗ F N , F N ( t ) = N + 1 1 k = 0 ∑ N D k ( t )
を考える(f ∗ F N ( x ) = 1 2 π ∫ − π π f ( x − t ) F N ( t ) d t f \ast F_N(x) = \frac{1}{2\pi}\int_{-\pi}^\pi f(x - t)F_N(t)\ dt f ∗ F N ( x ) = 2 π 1 ∫ − π π f ( x − t ) F N ( t ) d t )。F N F_N F N をフェイエール核 (Fejér kernel) という。
補題 6.11 F N ( t ) = 1 N + 1 ( sin ( N + 1 ) t 2 sin t 2 ) 2 F_N(t) = \dfrac{1}{N + 1}\left(\dfrac{\sin\frac{(N + 1)t}{2}}{\sin\frac{t}{2}}\right)^2 F N ( t ) = N + 1 1 ( sin 2 t sin 2 ( N + 1 ) t ) 2 。したがって (i) F N ≥ 0 F_N \geq 0 F N ≥ 0 、(ii) 1 2 π ∫ − π π F N = 1 \frac{1}{2\pi}\int_{-\pi}^\pi F_N = 1 2 π 1 ∫ − π π F N = 1 、(iii) 0 < δ ≤ ∣ t ∣ ≤ π 0 < \delta \leq \lvert t \rvert \leq \pi 0 < δ ≤ ∣ t ∣ ≤ π で F N ( t ) ≤ 1 ( N + 1 ) sin 2 ( δ / 2 ) → 0 F_N(t) \leq \frac{1}{(N + 1)\sin^2(\delta/2)} \to 0 F N ( t ) ≤ ( N + 1 ) s i n 2 ( δ /2 ) 1 → 0 。
証明. 2 sin t 2 sin ( k + 1 2 ) t = cos k t − cos ( k + 1 ) t 2\sin\frac{t}{2}\sin\left(k + \frac{1}{2}\right)t = \cos kt - \cos(k + 1)t 2 sin 2 t sin ( k + 2 1 ) t = cos k t − cos ( k + 1 ) t を k = 0 , … , N k = 0, \dots, N k = 0 , … , N で足すと 1 − cos ( N + 1 ) t = 2 sin 2 ( N + 1 ) t 2 1 - \cos(N + 1)t = 2\sin^2\frac{(N + 1)t}{2} 1 − cos ( N + 1 ) t = 2 sin 2 2 ( N + 1 ) t 。両辺を 2 sin 2 t 2 2\sin^2\frac{t}{2} 2 sin 2 2 t で割れば ∑ k D k ( t ) \sum_k D_k(t) ∑ k D k ( t ) の式を得る。(ii) は各 D k D_k D k の積分が 1 だから。(iii) は sin 2 ( t / 2 ) ≥ sin 2 ( δ / 2 ) \sin^2(t/2) \geq \sin^2(\delta/2) sin 2 ( t /2 ) ≥ sin 2 ( δ /2 ) から。□ \square □
(i)〜(iii) をみたす核の列を近似単位元 (approximate identity, good kernel) という。質量 1 の非負の山が原点に集中していく。
定理 6.12 (フェイエールの定理, Fejér's theorem)f f f が連続なら、σ N f \sigma_Nf σ N f は f f f に一様収束する。
証明. (ii) より σ N f ( x ) − f ( x ) = 1 2 π ∫ − π π ( f ( x − t ) − f ( x ) ) F N ( t ) d t \sigma_Nf(x) - f(x) = \frac{1}{2\pi}\int_{-\pi}^\pi\bigl(f(x - t) - f(x)\bigr)F_N(t)\ dt σ N f ( x ) − f ( x ) = 2 π 1 ∫ − π π ( f ( x − t ) − f ( x ) ) F N ( t ) d t 。ε > 0 \varepsilon > 0 ε > 0 に対し、一様連続性から ∣ t ∣ < δ \lvert t \rvert < \delta ∣ t ∣ < δ なら ∣ f ( x − t ) − f ( x ) ∣ < ε \lvert f(x - t) - f(x) \rvert < \varepsilon ∣ f ( x − t ) − f ( x )∣ < ε (x x x によらない)となる δ \delta δ をとる。∣ t ∣ < δ \lvert t \rvert < \delta ∣ t ∣ < δ の部分は (i)(ii) より ε \varepsilon ε 以下、δ ≤ ∣ t ∣ ≤ π \delta \leq \lvert t \rvert \leq \pi δ ≤ ∣ t ∣ ≤ π の部分は (iii) より 2 sup ∣ f ∣ / ( ( N + 1 ) sin 2 ( δ / 2 ) ) 2\sup\lvert f \rvert/\bigl((N + 1)\sin^2(\delta/2)\bigr) 2 sup ∣ f ∣ / ( ( N + 1 ) sin 2 ( δ /2 ) ) 以下である。よって lim sup N sup x ∣ σ N f ( x ) − f ( x ) ∣ ≤ ε \limsup_N\sup_x\lvert \sigma_Nf(x) - f(x) \rvert \leq \varepsilon lim sup N sup x ∣ σ N f ( x ) − f ( x )∣ ≤ ε 。□ \square □
系 6.13
(一意性)連続関数 f f f のフーリエ係数がすべて 0 なら f = 0 f = 0 f = 0 。
三角多項式全体は、周期 2 π 2\pi 2 π の連続関数全体の中で一様ノルムについて稠密である。
連続関数 f f f について、S N f ( x ) S_Nf(x) S N f ( x ) がある点 x x x で収束すれば、その極限は f ( x ) f(x) f ( x ) である。
証明. 1:σ N f = 0 \sigma_Nf = 0 σ N f = 0 で、f = lim σ N f = 0 f = \lim\sigma_Nf = 0 f = lim σ N f = 0 。2:σ N f \sigma_Nf σ N f は三角多項式である。3:収束する数列の算術平均は同じ極限に収束するので、lim S N f ( x ) = lim σ N f ( x ) = f ( x ) \lim S_Nf(x) = \lim\sigma_Nf(x) = f(x) lim S N f ( x ) = lim σ N f ( x ) = f ( x ) 。□ \square □
系 6.14 (ワイエルシュトラスの近似定理)[ a , b ] [a, b] [ a , b ] 上の連続関数は、多項式で一様に近似できる。
証明. 1 次変換で [ a , b ] = [ 0 , π ] [a, b] = [0, \pi] [ a , b ] = [ 0 , π ] としてよい。f f f を偶関数として [ − π , π ] [-\pi, \pi] [ − π , π ] に拡張し周期 2 π 2\pi 2 π に延長すると連続関数 g g g を得る。系 6.13 の 2 より sup ∣ g − P ∣ < ε / 2 \sup\lvert g - P \rvert < \varepsilon/2 sup ∣ g − P ∣ < ε /2 となる三角多項式 P P P がある。e i n x e^{inx} e in x のテイラー級数は [ − π , π ] [-\pi, \pi] [ − π , π ] 上で一様収束するので、P P P のテイラー展開を十分先で打ち切った多項式 Q Q Q は sup [ − π , π ] ∣ P − Q ∣ < ε / 2 \sup_{[-\pi, \pi]}\lvert P - Q \rvert < \varepsilon/2 sup [ − π , π ] ∣ P − Q ∣ < ε /2 をみたす。実部をとれば、[ 0 , π ] [0, \pi] [ 0 , π ] 上で ∣ f − Re Q ∣ < ε \lvert f - \operatorname{Re}Q \rvert < \varepsilon ∣ f − Re Q ∣ < ε 。□ \square □
6.6 L 2 L^2 L 2 理論:完全性とパーセヴァルの等式
定理 6.15 (パーセヴァルの等式, Parseval's identity)リーマン積分可能な f f f について ∥ f − S N f ∥ 2 → 0 \lVert f - S_Nf \rVert_2 \to 0 ∥ f − S N f ∥ 2 → 0 であり、
1 2 π ∫ − π π ∣ f ( x ) ∣ 2 d x = ∑ 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 2 π 1 ∫ − π π ∣ f ( x ) ∣ 2 d x = n ∈ Z ∑ ∣ f ^ ( n ) ∣ 2
が成り立つ。より一般に ⟨ f , g ⟩ = ∑ n f ^ ( n ) g ^ ( n ) ‾ \langle f, g \rangle = \sum_n\hat{f}(n)\overline{\hat{g}(n)} ⟨ f , g ⟩ = ∑ n f ^ ( n ) g ^ ( n ) 。
証明. まず f f f が連続のとき:フェイエールの定理より、ε > 0 \varepsilon > 0 ε > 0 に対し sup ∣ f − σ M f ∣ < ε \sup\lvert f - \sigma_Mf \rvert < \varepsilon sup ∣ f − σ M f ∣ < ε となる M M M がある。N ≥ M N \geq M N ≥ M なら σ M f ∈ T N \sigma_Mf \in \mathcal{T}_N σ M f ∈ T N なので、命題 6.4 より ∥ f − S N f ∥ 2 ≤ ∥ f − σ M f ∥ 2 ≤ sup ∣ f − σ M f ∣ < ε \lVert f - S_Nf \rVert_2 \leq \lVert f - \sigma_Mf \rVert_2 \leq \sup\lvert f - \sigma_Mf \rvert < \varepsilon ∥ f − S N f ∥ 2 ≤ ∥ f − σ M f ∥ 2 ≤ sup ∣ f − σ M f ∣ < ε 。
一般の f f f (∣ f ∣ ≤ B \lvert f \rvert \leq B ∣ f ∣ ≤ B )では、定理 6.6 の証明の階段関数 s s s (∣ s ∣ ≤ B \lvert s \rvert \leq B ∣ s ∣ ≤ B , ∫ ∣ f − s ∣ < ε \int\lvert f - s \rvert < \varepsilon ∫ ∣ f − s ∣ < ε )の跳びを短い線分でつなぎ、両端の値をそろえると、周期的な連続関数 h h h で ∣ h ∣ ≤ B \lvert h \rvert \leq B ∣ h ∣ ≤ B , ∫ − π π ∣ f − h ∣ < 2 ε \int_{-\pi}^\pi\lvert f - h \rvert < 2\varepsilon ∫ − π π ∣ f − h ∣ < 2 ε となるものが得られる。すると ∥ f − h ∥ 2 2 ≤ 2 B 2 π ∫ ∣ f − h ∣ < 2 B ε π \lVert f - h \rVert_2^2 \leq \frac{2B}{2\pi}\int\lvert f - h \rvert < \frac{2B\varepsilon}{\pi} ∥ f − h ∥ 2 2 ≤ 2 π 2 B ∫ ∣ f − h ∣ < π 2 B ε 。命題 6.4 より ∥ f − S N f ∥ 2 ≤ ∥ f − S N h ∥ 2 ≤ ∥ f − h ∥ 2 + ∥ h − S N h ∥ 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 ∥ f − S N f ∥ 2 ≤ ∥ f − S N h ∥ 2 ≤ ∥ f − h ∥ 2 + ∥ h − S N h ∥ 2 で、第 2 項は前半より 0 に収束する。よって ∥ f − S N f ∥ 2 → 0 \lVert f - S_Nf \rVert_2 \to 0 ∥ f − S N f ∥ 2 → 0 で、命題 6.4 の等式からパーセヴァルの等式を得る。一般形は ∣ ⋅ ∣ 2 \lvert \cdot \rvert^2 ∣ ⋅ ∣ 2 の等式を f + g f + g f + g , f + i g f + ig f + i g などに適用すれば(偏極恒等式)得られる。□ \square □
{ e n } \lbrace e_n \rbrace { e n } は「完全」である:どの関数も L 2 L^2 L 2 の意味で展開でき、係数の 2 乗和がノルムの 2 乗に等しい。これは有限次元の正規直交基底の性質そのものである。
例 6.16 (バーゼル問題)f ( x ) = x f(x) = x f ( x ) = x に適用すると、例 6.2 より ∣ f ^ ( ± n ) ∣ 2 = b n 2 / 4 = 1 / n 2 \lvert \hat{f}(\pm n) \rvert^2 = b_n^2/4 = 1/n^2 ∣ f ^ ( ± n ) ∣ 2 = b n 2 /4 = 1/ n 2 なので
1 2 π ∫ − π π x 2 d x = π 2 3 = 2 ∑ n = 1 ∞ 1 n 2 , よって ∑ n = 1 ∞ 1 n 2 = π 2 6 \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} 2 π 1 ∫ − π π x 2 d x = 3 π 2 = 2 n = 1 ∑ ∞ n 2 1 , よって n = 1 ∑ ∞ n 2 1 = 6 π 2
6.7 一様収束と項別微分
定理 6.18 f f f が連続かつ区分的に滑らかなら、∑ n ∣ f ^ ( n ) ∣ < ∞ \sum_n\lvert \hat{f}(n) \rvert < \infty ∑ n ∣ f ^ ( n )∣ < ∞ であり、S N f S_Nf S N f は f f f に一様収束する。
証明. f ′ f' f ′ (有限個の点を除いて定義される)はリーマン積分可能である。各小区間で部分積分すると、f f f の連続性と周期性から境界項は打ち消し合い、f ′ ^ ( n ) = i n f ^ ( n ) \widehat{f'}(n) = in\hat{f}(n) f ′ ( n ) = in f ^ ( n ) を得る。コーシー–シュワルツの不等式とベッセルの不等式(f ′ f' f ′ に適用)より
∑ n ≠ 0 ∣ f ^ ( n ) ∣ = ∑ n ≠ 0 ∣ f ′ ^ ( n ) ∣ ∣ n ∣ ≤ ( ∑ n ≠ 0 1 n 2 ) 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 n = 0 ∑ ∣ f ^ ( n )∣ = n = 0 ∑ ∣ n ∣ ∣ f ′ ( n )∣ ≤ n = 0 ∑ n 2 1 1/2 ( n ∑ ∣ f ′ ( n ) ∣ 2 ) 1/2 < ∞
よってワイエルシュトラスの M 判定法によりフーリエ級数は一様収束し、その和 g g g は連続である。一様収束なので項別に積分でき、g ^ ( n ) = f ^ ( n ) \hat{g}(n) = \hat{f}(n) g ^ ( n ) = f ^ ( n ) 。系 6.13 の 1 より g = f g = f g = f 。□ \square □
例 6.19 f ( x ) = x 2 f(x) = x^2 f ( x ) = x 2 (∣ x ∣ ≤ π \lvert x \rvert \leq \pi ∣ x ∣ ≤ π )の周期拡張は連続で区分的に滑らかなので、例 6.2 の級数はすべての x x x で f f f に収束する。x = π x = \pi x = π とおくと π 2 = π 2 3 + 4 ∑ 1 n 2 \pi^2 = \frac{\pi^2}{3} + 4\sum\frac{1}{n^2} π 2 = 3 π 2 + 4 ∑ n 2 1 となり、再び ∑ 1 n 2 = π 2 6 \sum\frac{1}{n^2} = \frac{\pi^2}{6} ∑ n 2 1 = 6 π 2 。x = 0 x = 0 x = 0 とおけば ∑ n ≥ 1 ( − 1 ) n + 1 n 2 = π 2 12 \sum_{n \geq 1}\frac{(-1)^{n+1}}{n^2} = \frac{\pi^2}{12} ∑ n ≥ 1 n 2 ( − 1 ) n + 1 = 12 π 2 。
滑らかさとフーリエ係数の減衰は、次のように対応する。
命題 6.20 (1) f f f が C k C^k C k 級(周期 2 π 2\pi 2 π )なら、f ( k ) ^ ( n ) = ( i n ) k f ^ ( n ) \widehat{f^{(k)}}(n) = (in)^k\hat{f}(n) f ( k ) ( n ) = ( in ) k f ^ ( n ) であり、∣ n ∣ k ∣ f ^ ( n ) ∣ → 0 \lvert n \rvert^k\lvert \hat{f}(n) \rvert \to 0 ∣ n ∣ k ∣ f ^ ( n )∣ → 0 かつ ∑ ∣ n ∣ 2 k ∣ f ^ ( n ) ∣ 2 < ∞ \sum\lvert n \rvert^{2k}\lvert \hat{f}(n) \rvert^2 < \infty ∑ ∣ n ∣ 2 k ∣ f ^ ( n ) ∣ 2 < ∞ 。(2) 連続関数 f f f について ∑ ∣ n ∣ k ∣ f ^ ( n ) ∣ < ∞ \sum\lvert n \rvert^k\lvert \hat{f}(n) \rvert < \infty ∑ ∣ n ∣ k ∣ f ^ ( n )∣ < ∞ なら、f f f は C k C^k C k 級で、フーリエ級数を k k k 回まで項別微分できる。
証明. (1) 部分積分を k k k 回行い、f ( k ) f^{(k)} f ( k ) にリーマン–ルベーグの補題とベッセルの不等式を適用する。(2) j ≤ k j \leq k j ≤ k 回項別微分した級数は M 判定法で一様収束するので、和は C k C^k C k 級で項別微分できる(微分積分学 第6章 )。和は f f f とフーリエ係数が一致する連続関数なので、系 6.13 の 1 より f f f に等しい。□ \square □
注意
項別微分はいつでも許されるわけではない。のこぎり波の級数 2 ∑ ( − 1 ) n + 1 n sin n x 2\sum\frac{(-1)^{n+1}}{n}\sin nx 2 ∑ n ( − 1 ) n + 1 sin n x を項別微分した 2 ∑ ( − 1 ) n + 1 cos n x 2\sum(-1)^{n+1}\cos nx 2 ∑ ( − 1 ) n + 1 cos n x は、一般項が 0 に収束しないのでどの点でも収束しない。周期拡張したのこぎり波は x = ± π x = \pm\pi x = ± π で跳ぶので、定理 6.18 の仮定をみたさないのである。一方、項別積分はリーマン積分可能な関数についてつねに許される。実際、コーシー–シュワルツの不等式と定理 6.15 より、∣ x ∣ ≤ π \lvert x \rvert \leq \pi ∣ x ∣ ≤ π で一様に ∣ ∫ 0 x ( f − S N f ) d t ∣ ≤ ∫ − π π ∣ f − S N f ∣ d t ≤ 2 π ∥ f − S N f ∥ 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 ∫ 0 x ( f − S N f ) d t ≤ ∫ − π π ∣ f − S N f ∣ d t ≤ 2 π ∥ f − S N f ∥ 2 → 0 である。
6.8 応用:等周不等式
長さ L L L の閉曲線で囲める面積 A A A は最大でいくらか。答えは円のときで A = L 2 / ( 4 π ) A = L^2/(4\pi) A = L 2 / ( 4 π ) である。
定理 6.21 (等周不等式, isoperimetric inequality)C 1 C^1 C 1 級の単純閉曲線(速さが 0 にならないもの)の長さを L L L 、囲む領域の面積を A A A とすると 4 π A ≤ L 2 4\pi A \leq L^2 4 π A ≤ L 2 であり、等号は円のときに限る。
証明. (フルヴィッツ)弧長に比例するパラメータで、反時計回りに z ( s ) = x ( s ) + i y ( s ) z(s) = x(s) + iy(s) z ( s ) = x ( s ) + i y ( s ) (0 ≤ s ≤ 2 π 0 \leq s \leq 2\pi 0 ≤ s ≤ 2 π )と表すと、∣ z ′ ( s ) ∣ = L / ( 2 π ) \lvert z'(s) \rvert = L/(2\pi) ∣ z ′ ( s )∣ = L / ( 2 π ) である。z z z は周期 2 π 2\pi 2 π の C 1 C^1 C 1 級関数なので、z = ∑ c n e i n s z = \sum c_ne^{ins} z = ∑ c n e in s 、z ′ ∼ ∑ i n c n e i n s z' \sim \sum inc_ne^{ins} z ′ ∼ ∑ in c n e in s 。パーセヴァルの等式より
L 2 4 π 2 = 1 2 π ∫ 0 2 π ∣ z ′ ∣ 2 d s = ∑ n n 2 ∣ c n ∣ 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 4 π 2 L 2 = 2 π 1 ∫ 0 2 π ∣ z ′ ∣ 2 d s = n ∑ n 2 ∣ c n ∣ 2
グリーンの定理より A = 1 2 ∫ 0 2 π ( x y ′ − y x ′ ) d s = 1 2 Im ∫ 0 2 π z ‾ z ′ d s A = \frac{1}{2}\int_0^{2\pi}(xy' - yx')\ ds = \frac{1}{2}\operatorname{Im}\int_0^{2\pi}\overline{z}z'\ ds A = 2 1 ∫ 0 2 π ( x y ′ − y x ′ ) d s = 2 1 Im ∫ 0 2 π z z ′ d s であり、パーセヴァルの等式の一般形より 1 2 π ∫ 0 2 π z ′ z ‾ d s = ∑ n i n c n c n ‾ \frac{1}{2\pi}\int_0^{2\pi}z'\overline{z}\ ds = \sum_n inc_n\overline{c_n} 2 π 1 ∫ 0 2 π z ′ z d s = ∑ n in c n c n なので A = π ∑ n n ∣ c n ∣ 2 A = \pi\sum_n n\lvert c_n \rvert^2 A = π ∑ n n ∣ c n ∣ 2 。よって
L 2 − 4 π A = 4 π 2 ∑ n ( n 2 − n ) ∣ c n ∣ 2 ≥ 0 L^2 - 4\pi A = 4\pi^2\sum_n(n^2 - n)\lvert c_n \rvert^2 \geq 0 L 2 − 4 π A = 4 π 2 n ∑ ( n 2 − n ) ∣ c n ∣ 2 ≥ 0
(整数 n n n について n 2 − n ≥ 0 n^2 - n \geq 0 n 2 − n ≥ 0 )。等号は n ≠ 0 , 1 n \neq 0, 1 n = 0 , 1 で c n = 0 c_n = 0 c n = 0 、すなわち z = c 0 + c 1 e i s z = c_0 + c_1e^{is} z = c 0 + c 1 e i s (円)のときに限る。□ \square □
6.9 発散する例
フェイエールの定理と定理 6.18 は、連続関数のフーリエ級数が各点で収束することまでは主張していない。実際、次が知られている(主張のみ)。
連続関数でフーリエ級数がある点で発散するものが存在する (デュ・ボア–レイモン, 1873 年)。S N f ( 0 ) = 1 2 π ∫ f D N S_Nf(0) = \frac{1}{2\pi}\int fD_N S N f ( 0 ) = 2 π 1 ∫ f D N は連続関数の空間上の有界線形汎関数で、そのノルムは 1 2 π ∫ ∣ D N ∣ \frac{1}{2\pi}\int\lvert D_N \rvert 2 π 1 ∫ ∣ D N ∣ (ルベーグ定数 )に等しく、N → ∞ N \to \infty N → ∞ で対数的に発散する(問題 6.4)。関数解析の一様有界性原理(バナッハ–シュタインハウスの定理)から、sup N ∣ S N f ( 0 ) ∣ = ∞ \sup_N\lvert S_Nf(0) \rvert = \infty sup N ∣ S N f ( 0 )∣ = ∞ となる連続関数 f f f の存在が従う(関数解析 第3章 )。
ルベーグ可積分な関数で、フーリエ級数がいたるところ発散するものが存在する(コルモゴロフ)。
一方、L 2 L^2 L 2 に属する関数のフーリエ級数はほとんどいたるところ収束する(カールソン, 1966 年)。これは非常に深い定理である。
各点収束は微妙だが、チェザロ平均による収束(フェイエール)と L 2 L^2 L 2 収束(パーセヴァル)は非常に扱いやすい。現代の解析学がフーリエ級数を L 2 L^2 L 2 の枠組みで扱う理由の一つである。
まとめ
フーリエ係数は正規直交系 e i n x e^{inx} e in x に関する内積であり、部分和 S N f S_Nf S N f は N N N 次三角多項式の中での L 2 L^2 L 2 最良近似である。ベッセルの不等式が成り立ち、係数は 0 に収束する(リーマン–ルベーグ)。
部分和はディリクレ核との畳み込みで、区分的に滑らかな関数では各点で左右の極限の平均に収束する。不連続点の近くでは約 9% の行き過ぎが残る(ギブス現象)。
チェザロ平均はフェイエール核(近似単位元)との畳み込みで、連続関数に一様収束する。ここから一意性と三角多項式の稠密性、ワイエルシュトラスの近似定理が従う。
パーセヴァルの等式 1 2 π ∫ ∣ f ∣ 2 = ∑ ∣ f ^ ( n ) ∣ 2 \frac{1}{2\pi}\int\lvert f \rvert^2 = \sum\lvert \hat{f}(n) \rvert^2 2 π 1 ∫ ∣ f ∣ 2 = ∑ ∣ f ^ ( n ) ∣ 2 が成り立ち、三角関数系は完全である。バーゼル問題 ∑ 1 / n 2 = π 2 / 6 \sum 1/n^2 = \pi^2/6 ∑ 1/ n 2 = π 2 /6 はその帰結である。
連続で区分的に滑らかな関数のフーリエ級数は絶対かつ一様に収束する。滑らかなほど係数は速く減衰する。
パーセヴァルの等式から等周不等式 4 π A ≤ L 2 4\pi A \leq L^2 4 π A ≤ L 2 が得られる。
連続関数のフーリエ級数が発散することもある。
演習問題
問題 6.1 ★ f ( x ) = ∣ sin x ∣ f(x) = \lvert \sin x \rvert f ( x ) = ∣ sin x ∣ のフーリエ級数を求め、∑ k = 1 ∞ 1 4 k 2 − 1 = 1 2 \sum_{k=1}^\infty\frac{1}{4k^2 - 1} = \frac{1}{2} ∑ k = 1 ∞ 4 k 2 − 1 1 = 2 1 を導け。
解答
偶関数なので b n = 0 b_n = 0 b n = 0 。a n = 2 π ∫ 0 π sin x cos n x d x = 1 π ∫ 0 π ( sin ( 1 + n ) x + sin ( 1 − n ) x ) d x a_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 a n = π 2 ∫ 0 π sin x cos n x d x = π 1 ∫ 0 π ( sin ( 1 + n ) x + sin ( 1 − n ) x ) d x 。n n n が偶数なら 1 ± n 1 \pm n 1 ± n は奇数で a n = 1 π ( 2 1 + n + 2 1 − n ) = 4 π ( 1 − n 2 ) a_n = \frac{1}{\pi}\left(\frac{2}{1 + n} + \frac{2}{1 - n}\right) = \frac{4}{\pi(1 - n^2)} a n = π 1 ( 1 + n 2 + 1 − n 2 ) = π ( 1 − n 2 ) 4 、n n n が奇数なら a n = 0 a_n = 0 a n = 0 (n = 1 n = 1 n = 1 も直接確かめられる)。よって ∣ sin x ∣ = 2 π − 4 π ∑ k ≥ 1 cos 2 k x 4 k 2 − 1 \lvert \sin x \rvert = \frac{2}{\pi} - \frac{4}{\pi}\sum_{k \geq 1}\frac{\cos 2kx}{4k^2 - 1} ∣ sin x ∣ = π 2 − π 4 ∑ k ≥ 1 4 k 2 − 1 c o s 2 k x (連続で区分的に滑らかなので定理 6.18 により等号)。x = 0 x = 0 x = 0 とおけば主張を得る。
問題 6.2 ★★ f ( x ) = e x f(x) = e^x f ( x ) = e x (− π < x < π -\pi < x < \pi − π < x < π )を周期 2 π 2\pi 2 π に拡張した関数のフーリエ係数を求め、ディリクレの定理を x = π x = \pi x = π で用いて ∑ n ∈ Z 1 1 + n 2 = π coth π \sum_{n \in \mathbb{Z}}\frac{1}{1 + n^2} = \pi\coth\pi ∑ n ∈ Z 1 + n 2 1 = π coth π を示せ。
解答
f ^ ( n ) = 1 2 π ∫ − π π e ( 1 − i n ) x d x = e ( 1 − i n ) π − e − ( 1 − i n ) π 2 π ( 1 − i n ) = ( − 1 ) n sinh π π ( 1 − i n ) \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)} f ^ ( n ) = 2 π 1 ∫ − π π e ( 1 − in ) x d x = 2 π ( 1 − in ) e ( 1 − in ) π − e − ( 1 − in ) π = π ( 1 − in ) ( − 1 ) n s i n h π 。x = π x = \pi x = π で左極限は e π e^\pi e π 、右極限は e − π e^{-\pi} e − π なので、定理 6.9 より S N f ( π ) → cosh π S_Nf(\pi) \to \cosh\pi S N f ( π ) → cosh π 。一方 S N f ( π ) = ∑ ∣ n ∣ ≤ N f ^ ( n ) ( − 1 ) n = sinh π π ∑ ∣ n ∣ ≤ N 1 + i n 1 + n 2 = sinh π π ∑ ∣ n ∣ ≤ N 1 1 + n 2 S_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} S N f ( π ) = ∑ ∣ n ∣ ≤ N f ^ ( n ) ( − 1 ) n = π s i n h π ∑ ∣ n ∣ ≤ N 1 + n 2 1 + in = π s i n h π ∑ ∣ n ∣ ≤ N 1 + n 2 1 (虚部は n n n と − n -n − n で打ち消し合う)。よって ∑ 1 1 + n 2 = π cosh π / sinh π \sum\frac{1}{1 + n^2} = \pi\cosh\pi/\sinh\pi ∑ 1 + n 2 1 = π cosh π / sinh π 。
問題 6.3 ★★ f ( x ) = x 2 f(x) = x^2 f ( x ) = x 2 にパーセヴァルの等式を適用して ∑ n = 1 ∞ 1 n 4 = π 4 90 \sum_{n=1}^\infty\frac{1}{n^4} = \frac{\pi^4}{90} ∑ n = 1 ∞ n 4 1 = 90 π 4 を示せ。
解答
f ^ ( 0 ) = π 2 3 \hat{f}(0) = \frac{\pi^2}{3} f ^ ( 0 ) = 3 π 2 、n ≠ 0 n \neq 0 n = 0 で f ^ ( n ) = a n 2 = 2 ( − 1 ) n n 2 \hat{f}(n) = \frac{a_n}{2} = \frac{2(-1)^n}{n^2} f ^ ( n ) = 2 a n = n 2 2 ( − 1 ) n 。1 2 π ∫ − π π x 4 d x = π 4 5 \frac{1}{2\pi}\int_{-\pi}^\pi x^4dx = \frac{\pi^4}{5} 2 π 1 ∫ − π π x 4 d x = 5 π 4 なので π 4 5 = π 4 9 + 2 ∑ n ≥ 1 4 n 4 \frac{\pi^4}{5} = \frac{\pi^4}{9} + 2\sum_{n \geq 1}\frac{4}{n^4} 5 π 4 = 9 π 4 + 2 ∑ n ≥ 1 n 4 4 、よって ∑ 1 n 4 = 1 8 ( π 4 5 − π 4 9 ) = π 4 90 \sum\frac{1}{n^4} = \frac{1}{8}\left(\frac{\pi^4}{5} - \frac{\pi^4}{9}\right) = \frac{\pi^4}{90} ∑ n 4 1 = 8 1 ( 5 π 4 − 9 π 4 ) = 90 π 4 。
問題 6.4 ★★ (ルベーグ定数)1 2 π ∫ − π π ∣ D N ( t ) ∣ d t ≥ 4 π 2 ∑ k = 1 N 1 k \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} 2 π 1 ∫ − π π ∣ D N ( t )∣ d t ≥ π 2 4 ∑ k = 1 N k 1 を示せ。特に左辺は N → ∞ N \to \infty N → ∞ で発散する。
解答
∣ sin ( t / 2 ) ∣ ≤ ∣ t ∣ / 2 \lvert \sin(t/2) \rvert \leq \lvert t \rvert/2 ∣ sin ( t /2 )∣ ≤ ∣ t ∣ /2 より ∣ D N ( t ) ∣ ≥ 2 ∣ sin ( N + 1 2 ) t ∣ ∣ t ∣ \lvert D_N(t) \rvert \geq \frac{2\lvert \sin(N + \frac{1}{2})t \rvert}{\lvert t \rvert} ∣ D N ( t )∣ ≥ ∣ t ∣ 2 ∣ s i n ( N + 2 1 ) t ∣ 。偶関数なので
1 2 π ∫ − π π ∣ D N ∣ ≥ 2 π ∫ 0 π ∣ sin ( N + 1 2 ) t ∣ t d t = 2 π ∫ 0 ( N + 1 2 ) π ∣ sin u ∣ u d u ≥ 2 π ∑ k = 1 N 1 k π ∫ ( k − 1 ) π k π ∣ sin u ∣ d u = 4 π 2 ∑ k = 1 N 1 k \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} 2 π 1 ∫ − π π ∣ D N ∣ ≥ π 2 ∫ 0 π t ∣ sin ( N + 2 1 ) t ∣ d t = π 2 ∫ 0 ( N + 2 1 ) π u ∣ sin u ∣ d u ≥ π 2 k = 1 ∑ N k π 1 ∫ ( k − 1 ) π k π ∣ sin u ∣ d u = π 2 4 k = 1 ∑ N k 1
問題 6.5 ★★ (ヴィルティンガーの不等式)f f f が周期 2 π 2\pi 2 π の C 1 C^1 C 1 級関数で ∫ 0 2 π f d x = 0 \int_0^{2\pi}f\ dx = 0 ∫ 0 2 π f d x = 0 なら ∫ 0 2 π ∣ f ∣ 2 d x ≤ ∫ 0 2 π ∣ f ′ ∣ 2 d x \int_0^{2\pi}\lvert f \rvert^2dx \leq \int_0^{2\pi}\lvert f' \rvert^2dx ∫ 0 2 π ∣ f ∣ 2 d x ≤ ∫ 0 2 π ∣ f ′ ∣ 2 d x であり、等号は f = a cos x + b sin x f = a\cos x + b\sin x f = a cos x + b sin x のときに限ることを示せ。
解答
f ^ ( 0 ) = 0 \hat{f}(0) = 0 f ^ ( 0 ) = 0 と f ′ ^ ( n ) = i n f ^ ( n ) \widehat{f'}(n) = in\hat{f}(n) f ′ ( n ) = in f ^ ( n ) およびパーセヴァルの等式より 1 2 π ∫ ∣ f ′ ∣ 2 = ∑ n ≠ 0 n 2 ∣ f ^ ( n ) ∣ 2 ≥ ∑ n ≠ 0 ∣ f ^ ( n ) ∣ 2 = 1 2 π ∫ ∣ 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 2 π 1 ∫ ∣ f ′ ∣ 2 = ∑ n = 0 n 2 ∣ f ^ ( n ) ∣ 2 ≥ ∑ n = 0 ∣ f ^ ( n ) ∣ 2 = 2 π 1 ∫ ∣ f ∣ 2 。等号は ∣ n ∣ ≥ 2 \lvert n \rvert \geq 2 ∣ n ∣ ≥ 2 で f ^ ( n ) = 0 \hat{f}(n) = 0 f ^ ( n ) = 0 、すなわち f = c 1 e i x + c − 1 e − i x f = c_1e^{ix} + c_{-1}e^{-ix} f = c 1 e i x + c − 1 e − i x のときに限る(定理 6.18 より級数は f f f に収束する)。これは第5章の最小固有値の評価の周期境界条件版である。
問題 6.6 ★★ リーマン–ルベーグの補題を用いて lim N → ∞ ∫ 0 ( N + 1 2 ) π sin u u d u = π 2 \lim_{N \to \infty}\int_0^{(N + \frac{1}{2})\pi}\frac{\sin u}{u}du = \frac{\pi}{2} lim N → ∞ ∫ 0 ( N + 2 1 ) π u s i n u d u = 2 π を示せ。広義積分 ∫ 0 ∞ sin u u d u \int_0^\infty\frac{\sin u}{u}du ∫ 0 ∞ u s i n u d u が収束することを認めて、その値を求めよ。
解答
1 2 π ∫ − π π D N = 1 \frac{1}{2\pi}\int_{-\pi}^\pi D_N = 1 2 π 1 ∫ − π π D N = 1 と偶関数性から ∫ 0 π sin ( N + 1 2 ) t sin ( t / 2 ) d t = π \int_0^\pi\frac{\sin(N + \frac{1}{2})t}{\sin(t/2)}dt = \pi ∫ 0 π s i n ( t /2 ) s i n ( N + 2 1 ) t d t = π 。h ( t ) = 1 sin ( t / 2 ) − 2 t h(t) = \frac{1}{\sin(t/2)} - \frac{2}{t} h ( t ) = s i n ( t /2 ) 1 − t 2 は ( 0 , π ] (0, \pi] ( 0 , π ] で連続で、t → 0 t \to 0 t → 0 で h ( t ) = t − 2 sin ( t / 2 ) t sin ( t / 2 ) → 0 h(t) = \frac{t - 2\sin(t/2)}{t\sin(t/2)} \to 0 h ( t ) = t s i n ( t /2 ) t − 2 s i n ( t /2 ) → 0 なので [ 0 , π ] [0, \pi] [ 0 , π ] 上の連続関数に拡張できる。リーマン–ルベーグの補題より ∫ 0 π h ( t ) sin ( N + 1 2 ) t d t → 0 \int_0^\pi h(t)\sin(N + \frac{1}{2})t\ dt \to 0 ∫ 0 π h ( t ) sin ( N + 2 1 ) t d t → 0 、よって ∫ 0 π 2 sin ( N + 1 2 ) t t d t → π \int_0^\pi\frac{2\sin(N + \frac{1}{2})t}{t}dt \to \pi ∫ 0 π t 2 s i n ( N + 2 1 ) t d t → π 。u = ( N + 1 2 ) t u = (N + \frac{1}{2})t u = ( N + 2 1 ) t と置換すれば主張を得る。広義積分が収束すれば、その値は部分列の極限 π 2 \frac{\pi}{2} 2 π に等しい。
問題 6.7 ★★★ (ワイルの一様分布定理)α \alpha α を無理数とする。周期 2 π 2\pi 2 π の任意の連続関数 f f f について lim N → ∞ 1 N ∑ k = 1 N f ( 2 π k α ) = 1 2 π ∫ 0 2 π f ( x ) d x \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 lim N → ∞ N 1 ∑ k = 1 N f ( 2 π k α ) = 2 π 1 ∫ 0 2 π f ( x ) d x を示せ。
解答
f = e m f = e_m f = e m (m ≠ 0 m \neq 0 m = 0 )なら、e 2 π i m α ≠ 1 e^{2\pi im\alpha} \neq 1 e 2 π im α = 1 (α \alpha α は無理数)なので
∣ 1 N ∑ k = 1 N e 2 π i m k α ∣ = 1 N ∣ e 2 π i m α ( e 2 π i m N α − 1 ) e 2 π i m α − 1 ∣ ≤ 2 N ∣ e 2 π i m α − 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 N 1 k = 1 ∑ N e 2 π imk α = N 1 e 2 π im α − 1 e 2 π im α ( e 2 π im N α − 1 ) ≤ N ∣ e 2 π im α − 1 ∣ 2 → 0
で、右辺の積分も 0。m = 0 m = 0 m = 0 なら両辺とも 1。線形性から三角多項式について成り立つ。一般の連続関数 f f f と ε > 0 \varepsilon > 0 ε > 0 に対し、系 6.13 の 2 より sup ∣ f − P ∣ < ε \sup\lvert f - P \rvert < \varepsilon sup ∣ f − P ∣ < ε となる三角多項式 P P P をとると、左辺の平均と右辺の積分はそれぞれ P P P の場合と ε \varepsilon ε 以内の差しかないので、lim sup \limsup lim sup で 2 ε 2\varepsilon 2 ε 以内となり主張が従う。区間の定義関数を連続関数で上下から近似すれば、k α k\alpha k α の小数部分が [ 0 , 1 ) [0, 1) [ 0 , 1 ) に一様に分布することがわかる。
問題 6.8 ★ [ 0 , π ] [0, \pi] [ 0 , π ] 上の f ( x ) = x ( π − x ) f(x) = x(\pi - x) f ( x ) = x ( π − x ) の正弦級数を求め、∑ k = 0 ∞ ( − 1 ) k ( 2 k + 1 ) 3 = π 3 32 \sum_{k=0}^\infty\frac{(-1)^k}{(2k + 1)^3} = \frac{\pi^3}{32} ∑ k = 0 ∞ ( 2 k + 1 ) 3 ( − 1 ) k = 32 π 3 を示せ。
解答
部分積分を 2 回行うと ∫ 0 π x ( π − x ) sin n x d x = 2 n 2 ∫ 0 π sin n x d x = 2 ( 1 − ( − 1 ) n ) n 3 \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} ∫ 0 π x ( π − x ) sin n x d x = n 2 2 ∫ 0 π sin n x d x = n 3 2 ( 1 − ( − 1 ) n ) なので b n = 8 π n 3 b_n = \frac{8}{\pi n^3} b n = π n 3 8 (n n n 奇数)、0(n n n 偶数)。奇関数拡張は連続で区分的に滑らかなので、定理 6.18 より x ( π − x ) = 8 π ∑ k ≥ 0 sin ( 2 k + 1 ) x ( 2 k + 1 ) 3 x(\pi - x) = \frac{8}{\pi}\sum_{k \geq 0}\frac{\sin(2k + 1)x}{(2k + 1)^3} x ( π − x ) = π 8 ∑ k ≥ 0 ( 2 k + 1 ) 3 s i n ( 2 k + 1 ) x (0 ≤ x ≤ π 0 \leq x \leq \pi 0 ≤ x ≤ π )。x = π / 2 x = \pi/2 x = π /2 とおくと π 2 4 = 8 π ∑ k ( − 1 ) k ( 2 k + 1 ) 3 \frac{\pi^2}{4} = \frac{8}{\pi}\sum_k\frac{(-1)^k}{(2k + 1)^3} 4 π 2 = π 8 ∑ k ( 2 k + 1 ) 3 ( − 1 ) k 。