無限次元のピタゴラスの定理

概要

無限次元のピタゴラスの定理は、区間 $[-\pi,\pi]$ 上の連続関数 $f$ の Fourier 係数 $a_n,b_n$ について $\frac1\pi\int_{-\pi}^{\pi}f(x)^2dx=\frac{a_0^2}2+\sum_{n=1}^\infty(a_n^2+b_n^2)$ となる Parseval の等式である($f(-\pi)=f(\pi)$ は仮定しない)。関数の内積で見ると「長さの 2 乗は無限個の直交する向きの成分の 2 乗の和に等しい」ことを表し、部分和との差の 2 乗の積分が $0$ に近づくことと同値である。$f(x)=x$ に使うと $\sum_{n\ge1}\frac1{n^2}=\frac{\pi^2}6$、$f(x)=x^2$ に使うと $\sum_{n\ge1}\frac1{n^4}=\frac{\pi^4}{90}$ を得る。

$$\newcommand{C}[0]{\mathbb{C}} \newcommand{N}[0]{\mathbb{N}} \newcommand{Q}[0]{\mathbb{Q}} \newcommand{R}[0]{\mathbb{R}} \newcommand{Z}[0]{\mathbb{Z}} $$

前提知識: Fourier係数と最良近似, 三角関数の直交性(高校数学), 級数, 連続関数

高校での出発点:成分の 2 乗の和は長さの 2 乗

空間ベクトルでは、互いに直交する長さ $1$ のベクトル 3 本を座標軸にとると、成分の 2 乗の和が長さの 2 乗になる。ピタゴラスの定理を 3 次元で使ったものである。

座標軸を 3 本とも使う

$u=(1,2,3)$ とし、正規直交系 $e_1=(1,0,0)$、$e_2=\frac15(0,3,4)$、$e_3=\frac15(0,4,-3)$ をとる(どの 2 つの内積も $0$、どれも長さ $1$)。成分は
$$ \langle u,e_1\rangle=1,\qquad \langle u,e_2\rangle=\frac{6+12}5=\frac{18}5,\qquad \langle u,e_3\rangle=\frac{8-9}5=-\frac15 $$
であり、2 乗の和は $1+\frac{324}{25}+\frac1{25}=\frac{350}{25}=14$ である。一方 $\|u\|^2=1+4+9=14$ で、両者は等しい。
$e_3$ を使わずに $e_1,e_2$ だけにすると、2 乗の和は $\frac{349}{25}$ で、$14$ より $\frac1{25}$ だけ小さい。座標軸が足りないと、取り残しが出る(Fourier係数と最良近似 の Bessel の不等式)。

三角多項式では等号になる

$T(x)=1+3\cos x-2\sin2x$ の Fourier 係数は $a_0=2$(定数項 $1=\frac{a_0}2$)、$a_1=3$、$b_2=-2$ で、他はすべて $0$ である。このとき
$$ \frac{a_0^2}2+\sum_{n\ge1}(a_n^2+b_n^2)=\frac42+9+4=15 $$
であり、三角関数の直交性(高校数学) で計算した $\int_{-\pi}^{\pi}T(x)^2dx=15\pi$ を $\pi$ で割った値と等しい。三角多項式は有限個の「座標軸」だけでちょうど表せるので、取り残しがない。

関数 $f$ が三角多項式でなければ、$f$ を表すには無限個の座標軸 $1,\cos x,\sin x,\cos2x,\sin2x,\ldots$ がすべて要る。それでも「成分の 2 乗の和=長さの 2 乗」は成り立つだろうか。この記事で答える問いは次の 3 つである。

  1. 連続関数でも、Fourier 係数の 2 乗の和は長さの 2 乗に等しいか。→ thm-idp-parseval
  2. それは何の役に立つか。→ thm-idp-basel、ex-idp-fourth-powers
  3. この等式から言えないことは何か。→ ex-idp-pointwise
    高校の計算この記事の言葉大学の言葉
    $a_1^2+a_2^2+a_3^2=\vec a\cdot\vec a$Parseval の等式完全正規直交系
    座標軸が 3 本で足りる三角関数の系で足りる完全性
    長さの 2 乗の差が $0$取り残し $\int(f-S_Nf)^2dx\to0$平均二乗収束($L^2$ 収束)

Parseval の等式

以下、$C([-\pi,\pi])$ は区間 $[-\pi,\pi]$ 上の実数値連続関数の全体で、$\|f\|^2=\int_{-\pi}^{\pi}f(x)^2dx$ とする。$a_n,b_n$ は $f$ の Fourier 係数、$S_Nf$ は Fourier 級数の第 $N$ 部分和である(Fourier係数と最良近似)。

Parseval の等式

$f\in C([-\pi,\pi])$ ならば
$$ \frac1\pi\int_{-\pi}^{\pi}f(x)^2\,dx=\frac{a_0^2}2+\sum_{n=1}^\infty\bigl(a_n^2+b_n^2\bigr) $$
が成り立つ。これは $\displaystyle\lim_{N\to\infty}\int_{-\pi}^{\pi}\bigl(f(x)-S_Nf(x)\bigr)^2dx=0$ と同値である。$f(-\pi)=f(\pi)$ は仮定しない。

有限個の数の列についての Parseval の等式は 離散Fourier変換と反転公式 で扱う。
左辺は長さの 2 乗の $\frac1\pi$ 倍、右辺は成分の 2 乗の和の $\frac1\pi$ 倍である(Fourier係数と最良近似 の命題「部分和は正射影である」)。Bessel の不等式の「$\le$」が、無限個の成分をすべて集めると「$=$」になる。これが「無限次元のピタゴラスの定理」である。

証明の準備

証明では、次の近似定理だけを証明せずに使う。

三角多項式による一様近似(引用)

$g\in C([-\pi,\pi])$ が $g(-\pi)=g(\pi)$ を満たすならば、任意の $\varepsilon>0$ に対して、すべての $x\in[-\pi,\pi]$ で $|g(x)-T(x)|<\varepsilon$ となる三角多項式 $T$ が存在する。

引用元

thm-idp-weierstrass は本記事では証明しない。Leb26 の 11.7 節(Stone–Weierstrass の定理)と演習 11.7.8 にあり、そこでは複素指数関数 $e^{inx}$ の 1 次結合による近似として述べられているので、その実部をとればよい。Fejér の定理からも従う。

三角多項式 $T$ はどれも $T(-\pi)=T(\pi)$ を満たす($\cos n(-\pi)=\cos n\pi$、$\sin n(\pm\pi)=0$)。したがって、$f(x)=x$ のように $f(-\pi)\ne f(\pi)$ となる関数は、三角多項式で各点で一様に近似することはできない。実際、すべての $x$ で $|f(x)-T(x)|<\varepsilon$ となる三角多項式 $T$ があったとすると、$T(\pi)=T(-\pi)$ から $|f(\pi)-f(-\pi)|\le|f(\pi)-T(\pi)|+|T(-\pi)-f(-\pi)|<2\varepsilon$ となり、$\varepsilon\le\frac12|f(\pi)-f(-\pi)|$ のとき矛盾する。そこで証明では、$f$ を端の近くだけ取り替えた関数 $g$ を使う(図 1)。
図1:端点の修正。g は右端の幅 δ の所だけ 1 次関数に取り替えたもので、両端の値が等しい 図1:端点の修正。g は右端の幅 δ の所だけ 1 次関数に取り替えたもので、両端の値が等しい

端点の修正を f(x)=x で行う

$f(x)=x$、$\delta=\frac12$ とする。$g$ を、$[-\pi,\pi-\frac12]$ では $g(x)=x$、$[\pi-\frac12,\pi]$ では点 $\bigl(\pi-\frac12,\pi-\frac12\bigr)$ と点 $(\pi,-\pi)$ を結ぶ 1 次関数とする。

  1. $g$ は $x=\pi-\frac12$ でつながっているので連続であり、$g(\pi)=-\pi=g(-\pi)$ である。
  2. $g$ の値は $-\pi$ と $\pi-\frac12$ の間にあるので、$|g|\le\pi$ である。
  3. $f$ と $g$ が違うのは $[\pi-\frac12,\pi]$ だけで、そこでは $|f(x)-g(x)|\le|f(x)|+|g(x)|\le2\pi$ なので、$\|f-g\|^2\le(2\pi)^2\cdot\frac12=2\pi^2=19.7\ldots$ である。直接積分すると $\|f-g\|^2=\frac{4\pi^2}3\cdot\frac12=6.57\ldots$ であり、確かにこの上限以下である。
    $\delta$ を小さくすれば $\|f-g\|^2$ はいくらでも小さくなる。

もう 1 つ、長さについての簡単な不等式を用意する。

和の長さの 2 乗の上限

内積空間の元 $u,v$ について $\|u+v\|^2\le2\|u\|^2+2\|v\|^2$ が成り立つ。

2 つの展開を足す

内積の双線形性と対称性から
$$ \|u+v\|^2=\|u\|^2+2\langle u,v\rangle+\|v\|^2,\qquad \|u-v\|^2=\|u\|^2-2\langle u,v\rangle+\|v\|^2 $$
である(積分で定める内積 の展開公式)。2 式を足すと $\|u+v\|^2+\|u-v\|^2=2\|u\|^2+2\|v\|^2$ となる。$\|u-v\|^2\ge0$ なので、$\|u+v\|^2\le2\|u\|^2+2\|v\|^2$ である。$\square$

和の長さの 2 乗の上限を数で確かめる
  1. $\mathbb{R}^2$ で $u=(1,0)$、$v=(1,1)$ とする。$u+v=(2,1)$ なので $\|u+v\|^2=4+1=5$、$2\|u\|^2+2\|v\|^2=2\cdot1+2\cdot2=6$ で、確かに $5\le6$ である。差 $6-5=1$ は、証明で捨てた $\|u-v\|^2=\|(0,-1)\|^2=1$ に等しい。
  2. 等号になる例:$u=v=(1,1)$ とすると $\|u+v\|^2=\|(2,2)\|^2=8$、$2\|u\|^2+2\|v\|^2=2\cdot2+2\cdot2=8$ で、等号が成り立つ。$u-v=0$ なので、捨てた $\|u-v\|^2$ が $0$ だからである。

証明

端点を直して一様近似し、最良近似と比べる

方針:$\|f-S_Nf\|\to0$ を示す。まず $f$ に近い三角多項式 $T$ を 1 つ見つけ(段 3〜5)、$S_Nf$ は $T$ 以上に $f$ に近いこと(段 2)を使う。
段 1(同値性):Fourier係数と最良近似 の系「Fourier 係数の Bessel の不等式と最良近似」の 1 から
$$ \int_{-\pi}^{\pi}(f-S_Nf)^2dx=\int_{-\pi}^{\pi}f^2dx-\pi\Bigl(\frac{a_0^2}2+\sum_{n=1}^N(a_n^2+b_n^2)\Bigr) $$
である。$N\to\infty$ で左辺が $0$ に近づくことと、右辺の括弧が $\frac1\pi\int f^2dx$ に近づくこと(等式が成り立つこと)は同じである。以下、左辺が $0$ に近づくことを示す。
段 2($N$ を大きくしても悪くならない):$T$ が次数 $N_0$ 以下の三角多項式なら、$M\ge N_0$ のとき $T$ は次数 $M$ 以下でもある。同じ系の 3($S_Mf$ が次数 $M$ 以下でいちばん近い)から、$\|f-S_Mf\|\le\|f-T\|$ である。
段 3(端点の修正):$f$ は閉区間で連続なので、$|f|$ は最大値をもつ。それを $K$ とする。$0<\delta\le\pi$ とし、$g$ を、$[-\pi,\pi-\delta]$ では $f$ に等しく、$[\pi-\delta,\pi]$ では点 $(\pi-\delta,f(\pi-\delta))$ と点 $(\pi,f(-\pi))$ を結ぶ 1 次関数とする。

  • $g$ は $x=\pi-\delta$ でつながるので連続で、$g(\pi)=f(-\pi)=g(-\pi)$ である。
  • 1 次関数の部分の値は、両端の値 $f(\pi-\delta)$ と $f(-\pi)$ の間にあり、両端の絶対値は $K$ 以下なので $|g|\le K$ である。
  • $f$ と $g$ が違うのは $[\pi-\delta,\pi]$ だけで、そこでは $|f-g|\le|f|+|g|\le2K$ である。よって
    $$ \|f-g\|^2=\int_{\pi-\delta}^{\pi}\bigl(f(x)-g(x)\bigr)^2dx\le(2K)^2\delta=4K^2\delta $$
    である。
    段 4(一様近似):$\varepsilon>0$ を与える。$4K^2\delta<\varepsilon$ となるように $\delta$ を選ぶ($K=0$ ならどの $\delta$ でもよい)。$g$ に thm-idp-weierstrass を使い、すべての $x$ で $|g(x)-T(x)|<\sqrt{\frac\varepsilon{2\pi}}$ となる三角多項式 $T$ をとる。すると $(g-T)^2<\frac\varepsilon{2\pi}$ なので、長さ $2\pi$ の区間で積分して $\|g-T\|^2<\frac\varepsilon{2\pi}\cdot2\pi=\varepsilon$ である。
    段 5(組み合わせる):$f-T=(f-g)+(g-T)$ に lem-idp-sum-bound を使うと
    $$ \|f-T\|^2\le2\|f-g\|^2+2\|g-T\|^2<2\varepsilon+2\varepsilon=4\varepsilon $$
    である。
    段 6(結論):$T$ の次数を $N_0$ とすると、段 2 から $M\ge N_0$ のすべての $M$ で $\|f-S_Mf\|^2\le\|f-T\|^2<4\varepsilon$ である。$\varepsilon>0$ は任意なので、$\|f-S_Nf\|^2\to0$($N\to\infty$)である。段 1 から Parseval の等式が成り立つ。$\square$
Fourier 係数による一意性

$f\in C([-\pi,\pi])$ の Fourier 係数がすべて $0$ ならば、$f$ は恒等的に $0$ である。したがって、Fourier 係数がすべて等しい 2 つの連続関数は等しい。

長さが 0 になる

thm-idp-parseval の右辺が $0$ なので $\int_{-\pi}^{\pi}f^2dx=0$ である。関数の内積の正定値性(積分で定める内積 の定理「関数の内積は内積である」。ここで $f$ の連続性を使う)から $f=0$ である。後半は、2 つの関数の差に前半を使えばよい(差の Fourier 係数は係数の差である)。$\square$

一意性を使って等式を確かめる

$f(x)=\cos^2x$ の Fourier 係数を計算すると、$a_0=\frac1\pi\int_{-\pi}^{\pi}\cos^2x\,dx=1$、$a_2=\frac1\pi\int_{-\pi}^{\pi}\cos^2x\cos2x\,dx=\frac12$ で、他はすべて $0$ である(たとえば $a_1=0$、$a_4=0$、$b_n=0$)。一方 $g(x)=\frac12+\frac12\cos2x$ の Fourier 係数も $a_0=1$、$a_2=\frac12$、他は $0$ である。係数がすべて一致するので、cor-idp-uniqueness により 2 つの連続関数は等しく、$\cos^2x=\frac12+\frac12\cos2x$ と結論できる。ただし、$f$ の係数 $\int_{-\pi}^{\pi}\cos^2x\cos nx\,dx$ などは、ふつう $\cos^2x$ を半角の公式や積和の公式で和に直して計算する。したがってこれは半角の公式の新しい証明ではなく、答えが分かっている式で、一意性の使い方(係数を比べれば関数が決まる)を確かめ直したものである。

応用:平方数の逆数の和

平方数の逆数の和

$$ \sum_{n=1}^\infty\frac1{n^2}=1+\frac14+\frac19+\cdots=\frac{\pi^2}6 . $$

f(x)=x に Parseval の等式を使う

$f(x)=x$ は $[-\pi,\pi]$ で連続である($f(-\pi)\ne f(\pi)$ だが、thm-idp-parseval はそれを仮定しない)。
段 1(左辺):$\frac1\pi\int_{-\pi}^{\pi}x^2\,dx=\frac1\pi\cdot\frac{2\pi^3}3=\frac{2\pi^2}3$ である。
段 2(右辺):Fourier係数と最良近似 で計算したとおり、$a_n=0$、$b_n=\frac{2(-1)^{n+1}}n$ なので $b_n^2=\frac4{n^2}$ である。右辺は $\sum_{n=1}^\infty\frac4{n^2}$ である。
段 3(比べる):Parseval の等式は $\frac{2\pi^2}3=4\sum_{n=1}^\infty\frac1{n^2}$ となり、両辺を $4$ で割って $\sum_{n=1}^\infty\frac1{n^2}=\frac{\pi^2}6$ を得る。$\square$

Bessel の不等式だけからは $\sum\frac1{n^2}\le\frac{\pi^2}6$ しか出ない。等号を保証するのが Parseval の等式である。同じ値の、三角関数の不等式による初等的な別証明は Basel問題 にある。

取り残しの数値

段 1 の式を $f(x)=x$ に使うと、取り残しは
$$ \int_{-\pi}^{\pi}(x-S_Nf)^2dx=\frac{2\pi^3}3-\pi\sum_{n=1}^N\frac4{n^2}=4\pi\Bigl(\frac{\pi^2}6-\sum_{n=1}^N\frac1{n^2}\Bigr) $$
である。数値で確かめる($\frac{\pi^2}6=1.6449340\ldots$)。

$N$$\sum_{n=1}^N\frac1{n^2}$$\frac{\pi^2}6-\sum_{n=1}^N\frac1{n^2}$$\int_{-\pi}^{\pi}(x-S_Nf)^2dx$
$1$$1$$0.6449\ldots$$8.104\ldots$
$10$$1.5497\ldots$$0.09516\ldots$$1.195\ldots$
$100$$1.63498\ldots$$0.009950\ldots$$0.1250\ldots$
$1000$$1.643934\ldots$$0.0009995\ldots$$0.01256\ldots$

$N=1$ の値 $8.104\ldots$ は、Fourier係数と最良近似 で求めた「$c\sin x$ の形でいちばん近いもの」までの距離の 2 乗と同じである。$N$ を $10$ 倍すると取り残しはほぼ $\frac1{10}$ になる(図 2)。

図2:取り残し ∫(x − S_N f)² dx と N の関係(両対数目盛)。N を 10 倍すると約 10 分の 1 になる 図2:取り残し ∫(x − S_N f)² dx と N の関係(両対数目盛)。N を 10 倍すると約 10 分の 1 になる

4 乗の逆数の和
  1. $f(x)=x^2$ に使う。$f$ は偶関数なので $b_n=0$ であり、$a_0=\frac1\pi\int_{-\pi}^{\pi}x^2dx=\frac{2\pi^2}3$ である。$n\ge1$ では部分積分を 2 回行って
    $$ \int_{-\pi}^{\pi}x^2\cos nx\,dx=\Bigl[\frac{x^2\sin nx}n\Bigr]_{-\pi}^{\pi}-\frac2n\int_{-\pi}^{\pi}x\sin nx\,dx=0-\frac2n\cdot\frac{2(-1)^{n+1}\pi}n=\frac{4(-1)^n\pi}{n^2} $$
    なので $a_n=\frac{4(-1)^n}{n^2}$ である($\int_{-\pi}^{\pi}x\sin nx\,dx=\pi b_n$ は $f(x)=x$ の計算から)。左辺は $\frac1\pi\int_{-\pi}^{\pi}x^4dx=\frac{2\pi^4}5$ なので、Parseval の等式は
    $$ \frac{2\pi^4}5=\frac12\Bigl(\frac{2\pi^2}3\Bigr)^2+\sum_{n=1}^\infty\frac{16}{n^4}=\frac{2\pi^4}9+16\sum_{n=1}^\infty\frac1{n^4} $$
    となる。$\frac{2\pi^4}5-\frac{2\pi^4}9=\frac{8\pi^4}{45}$ を $16$ で割って $\sum_{n=1}^\infty\frac1{n^4}=\frac{\pi^4}{90}$ を得る。
  2. $f(x)=|x|$ に使う。Fourier係数と最良近似 で計算したとおり $a_0=\pi$、奇数の $n$ で $a_n=-\frac4{\pi n^2}$、それ以外の係数は $0$ である。左辺は $\frac1\pi\int_{-\pi}^{\pi}x^2dx=\frac{2\pi^2}3$ なので
    $$ \frac{2\pi^2}3=\frac{\pi^2}2+\frac{16}{\pi^2}\sum_{n\text{ は奇数}}\frac1{n^4},\qquad \sum_{n\text{ は奇数}}\frac1{n^4}=1+\frac1{3^4}+\frac1{5^4}+\cdots=\frac{\pi^4}{96} $$
    である。これは 1 と矛盾しない:偶数の項の和は $\sum\frac1{(2m)^4}=\frac1{16}\cdot\frac{\pi^4}{90}$ なので、奇数の項の和は $\frac{15}{16}\cdot\frac{\pi^4}{90}=\frac{\pi^4}{96}$ となり、一致する。

反例と注意

外した条件崩れる主張ボックス
系が全部そろっている($\cos$ も $\sin$ も)等号(取り残しが $0$)ex-idp-incomplete
「長さで測る」近さ各点での近さex-idp-pointwise
$f$ の連続性係数から $f$ が決まることex-idp-discontinuous
反例:座標軸が欠けていると等号にならない

$\sin nx$($n\ge1$)だけを使い、$\cos$ の系を使わないとする。$f(x)=\cos x$ について、$\sin$ の係数はすべて $b_n=\frac1\pi\int_{-\pi}^{\pi}\cos x\sin nx\,dx=0$ なので、その 2 乗の和は $0$ である。一方 $\frac1\pi\int_{-\pi}^{\pi}\cos^2x\,dx=1$ である。$0<1$ で、等号は成り立たない。Parseval の等式は、$1,\cos nx,\sin nx$ をすべて使ったときの主張である。ex-idp-vector で $e_3$ を落としたときと同じことが起きている。

反例:長さで近いことは各点で近いことではない

thm-idp-parseval の収束 $\int(f-S_Nf)^2dx\to0$ を平均二乗収束という。これは、各点 $x$ で $S_Nf(x)\to f(x)$ となること(各点収束)とは別の主張である。

  1. $f(x)=x$ では、$\sin n\pi=0$ からすべての $N$ で $S_Nf(\pi)=0$ であり、$f(\pi)=\pi$ に近づかない(図 3)。それでも取り残し $\int(x-S_Nf)^2dx$ は $0$ に近づく(ex-idp-table)。$-\pi< x<\pi$ の各点では $S_Nf(x)\to x$ となるが、それは平均二乗収束からは出ず、各点収束の定理(Leb26 定理 11.8.10。本記事では証明しない)による。
  2. 連続関数 $g_k(x):=\max(0,1-k|x|)$($k=1,2,3,\dots$。原点で高さ $1$、幅 $\frac2k$ の山)は、$\|g_k\|^2=2\int_0^{1/k}(1-kx)^2dx=\frac2{3k}\to0$ を満たすので、長さの意味では $0$ に近づく。しかし $g_k(0)=1$ はずっと $0$ に近づかない。

図3:x = π の近くでの部分和 S_5 f、S_20 f、S_80 f。どの N でも S_N f(π) = 0 で、f(π) = π に近づかない 図3:x = π の近くでの部分和 S_5 f、S_20 f、S_80 f。どの N でも S_N f(π) = 0 で、f(π) = π に近づかない

反例:連続でない関数では係数から関数が決まらない

区間 $[-\pi,\pi]$ で、$x=0$ のとき $1$、それ以外で $0$ をとる関数 $h$ を考える。$h$ は連続でないが、積分は定まり、$h\cos nx$ も $h\sin nx$ も 1 点を除いて $0$ なので、Fourier 係数はすべて $0$ である。しかし $h$ は $0$ でない。cor-idp-uniqueness の証明で使った正定値性($\int h^2dx=0$ なら $h=0$)が、連続でない関数では成り立たないからである。

さらに先へ

  • 2 乗が積分できる関数まで広げ(ほとんどいたるところ等しい関数を同じとみなす)、極限をとっても外に出ないようにした空間 $L^2$ は Hilbert空間 の代表例であり、三角関数系はその「完全正規直交系」になる。Parseval の等式はその言いかえである(本記事では扱わない)。
  • thm-idp-weierstrass は、Stone–Weierstrassの定理 の特別な場合である。
  • $\sum\frac1{n^{2k}}$ はどれも $\pi^{2k}$ の有理数倍になる(Basel問題、Riemannゼータ関数)。

関連項目

参考文献

Mathpediaは寄付と、参考文献の書籍リンク(Amazonアソシエイト)の紹介料で運営されています。 支援について / 寄付する