MathLabs

解析学

調和解析学

関数や信号を基本波に分解する分野で、フーリエ級数をより広い舞台に一般化したもの。

直観波から周波数へ

私たちが耳にする音——バイオリンの音、声、交通騒音——は実際には時間とともに変化する一つの圧力信号にすぎない。しかし耳(やグラフィックイコライザー)は、それが異なる音高を持つ多くの純音からできていることを教えてくれる。調和解析とは、その分解を扱う数学であり、関数 ff を、それぞれ固有の周波数 ξ\xi を持つ純粋な振動 e2πixξe^{2\pi i x \xi} の和や積分として書き直す。フーリエ級数は周期信号に対して離散的な周波数の列でこれを行うが、調和解析はこの発想を非周期信号、高次元、さらには他の群(円周、有限群、リー群)へと一般化する。

可変数の調和成分から構成された波形を示す、フーリエ級数のインタラクティブなウィジェット。
最初の24個のフーリエ調和成分から再構成された矩形波に近い波形。n を動かすと、ジャンプ付近のギブス現象の『リンギング』が幅は狭くなるが高さは減らないことが見える。

大学円周上のフーリエ級数:収束の問題

周期関数 ff のフーリエ級数の第 NN 部分和を SNf(x)=∑∣n∣≤Nf^(n) einxS_N f(x) = \sum_{|n|\le N} \hat f(n)\, e^{inx} と書く。ここで f^(n)=12π∫−ππf(x) e−inx dx\hat f(n) = \frac{1}{2\pi}\int_{-\pi}^{\pi} f(x)\, e^{-inx}\,dx である。L2L^2 では万事がすっきりしている:指数関数 einxe^{inx} は正規直交基底をなし、パーセバルの等式 12π∫∣f∣2=∑∣f^(n)∣2\frac{1}{2\pi}\int |f|^2 = \sum |\hat f(n)|^2 が成り立ち、SNf→fS_N f \to f が L2L^2 ノルムで成り立つ。難しいのは各点収束である:ある特定の点 xx で SNf(x)→f(x)S_N f(x) \to f(x) となるだろうか。上で示した矩形波はすでにその難しさを示唆している——ジャンプの近くでは、NN がどれほど大きくても SNfS_N f は常にジャンプの高さの約9%だけオーバーシュートする(ギブス現象であり、オーバーシュート量は半ジャンプの 2πSi(π)≈1.179\frac{2}{\pi}\mathrm{Si}(\pi) \approx 1.179 倍に近づく)。

定義: ディリクレ核

部分和はそれ自体が畳み込みである: SNf(x)=(DN∗f)(x)=12π∫−ππDN(x−y) f(y) dyS_N f(x) = (D_N * f)(x) = \frac{1}{2\pi}\int_{-\pi}^{\pi} D_N(x-y)\, f(y)\, dy。ここでディリクレ核は DN(x)=∑∣n∣≤Neinx=sin⁡((N+12)x)sin⁡(x/2)D_N(x) = \sum_{|n|\le N} e^{inx} = \dfrac{\sin\left((N+\tfrac12)x\right)}{\sin(x/2)} である。したがって収束に関するあらゆる問いは、実質的に DND_N の形状についての問いになる。

DN(x)=∑∣n∣≤Neinx=sin⁡((N+12)x)sin⁡(x/2)D_N(x) = \sum_{|n|\le N} e^{inx} = \dfrac{\sin\left((N+\tfrac12)x\right)}{\sin(x/2)}
1周期にわたるディリクレ核・フェイエール核・ポアソン核のインタラクティブなグラフで、L1ノルムのリアルタイム推定値付き。
N=9N=9 のディリクレ核 DND_N:フェイエール核やポアソン核と異なり、ゼロより下に落ち込む部分がある(サンゴ色の部分)。その L1L^1 ノルム——ルベーグ定数——は有界のままではなく (4/π2)log⁡N(4/\pi^2)\log N のように増大する。核を切り替えたり NN を大きくしたりして比べてみよう。

フェイエール核 FNF_N(D0,…,DND_0,\ldots,D_N のチェザロ平均)とポアソン核 PrP_r はどちらも非負なので、正真正銘の近似単位元である:質量 11、L1L^1 ノルム 11、そして N→∞N\to\infty(あるいは r→1r\to1)のとき質量が 00 に集中する——これだけで、任意の連続関数 ff に対して FN∗f→fF_N * f\to f と Pr∗f→fP_r * f\to f が一様に成り立つことが保証される。ディリクレ核はそうではない:その負の部分のせいで 12π∫∣DN∣\frac{1}{2\pi}\int|D_N|(ルベーグ定数)は 4π2log⁡N\frac{4}{\pi^2}\log N のように増大し、N→∞N\to\infty で有界でなくなる。バナッハ–シュタインハウスの一様有界性原理により、ルベーグ定数が有界でないことは、ある点でフーリエ級数が発散する連続関数の存在を強制する——これはまさに du Bois-Reymond が1873年に手作業で構成した現象であり、ここでは労せずして得られる。コルモゴロフは後に、単に可積分なだけの関数についてはさらに悪いことが起こりうることを示した:1923年には概発散する L1L^1 関数を構成し、1926年には至る所発散するものを構成した。

ff が連続で周期 2π2\pi をもつならば、チェザロ平均 σNf=1N+1∑k=0NSkf=FN∗f\sigma_N f = \frac{1}{N+1}\sum_{k=0}^N S_k f = F_N * f は N→∞N\to\infty のとき ff に一様収束する。

なぜ正しいのか?

部分和を平均することで、まさに SNfS_N f を暴れさせている振動が打ち消される:DND_N には激しい負の部分があるのに対し、その移動平均 FNF_N は滑らかで非負の山であり、常に ff の各部分を正の重みで足し合わせるだけである。そのため SNfS_N f 自体は収束しないことがあっても、その移動平均は必ず落ち着く——これは微積分において振動する部分和の列をチェザロ平均によって手なずけるのと同じ技である。

証明

ステップ1(フェイエール核は正の近似単位元である)。D0,…,DND_0,\ldots,D_N を平均し、望遠鏡型の三角関数恒等式を用いると、閉じた形 FN(x)=1N+1(sin⁡(N+12x)sin⁡(x/2))2≥0F_N(x) = \frac{1}{N+1}\left(\dfrac{\sin\left(\tfrac{N+1}{2}x\right)}{\sin(x/2)}\right)^2 \ge 0 が得られる。FN=1N+1∑k=0NDkF_N = \frac{1}{N+1}\sum_{k=0}^N D_k かつ任意の kk について 12π∫Dk=1\frac{1}{2\pi}\int D_k = 1 なので、12π∫FN=1\frac{1}{2\pi}\int F_N = 1 でもある。非負であることから 12π∫∣FN∣=1\frac{1}{2\pi}\int |F_N| = 1 でもあり、したがって L1L^1 ノルムは決して爆発しない。

ステップ2(質量が 00 に集中する)。任意に固定した δ>0\delta>0 に対し、δ≤∣x∣≤π\delta \le |x| \le \pi では分母 sin⁡(x/2)2\sin(x/2)^2 が正の定数で下から評価できるので、その範囲で一様に FN(x)=O(1/N)F_N(x) = O(1/N) となる。[−δ,δ][-\delta,\delta] の外側の質量は N→∞N\to\infty で 00 に近づく。

ステップ3(近似単位元の評価)。12π∫FN=1\frac{1}{2\pi}\int F_N = 1 を用いて σNf(x)−f(x)=12π∫−ππFN(y)(f(x−y)−f(x)) dy\sigma_N f(x) - f(x) = \frac{1}{2\pi}\int_{-\pi}^{\pi} F_N(y)\big(f(x-y)-f(x)\big)\,dy と書く。積分を ∣y∣<δ|y|<\delta と δ≤∣y∣≤π\delta\le|y|\le\pi に分ける。前半では、ff の一様連続性により δ\delta を小さくとれば ∣f(x−y)−f(x)∣<ε|f(x-y)-f(x)|<\varepsilon となり、FN≥0F_N\ge0 かつ全質量 11 よりその部分は ε\varepsilon で抑えられる。後半では FN=O(1/N)F_N=O(1/N) かつ ff が有界なので、その部分は N→∞N\to\infty のとき →0\to 0 となる。両方の評価は xx について一様なので、σNf→f\sigma_N f\to f が一様に成り立つ。

f∈L2(T)f\in L^2(\mathbb{T}) ならば、ほとんどすべての xx で SNf(x)→f(x)S_N f(x)\to f(x) となる。(Hunt は1968年にこれを p>1p>1 を満たすすべての f∈Lp(T)f\in L^p(\mathbb{T}) に拡張した。)

なぜ正しいのか?

コルモゴロフによる、至る所発散する1926年の L1L^1 関数の例を踏まえると、L2L^2——わずかに強い条件にすぎない——も同じように、あるいは少なくともほとんど至る所で失敗するだろうと考えるのはごく自然に思えた。カールソンの定理はまさに驚きである:二乗可積分というわずかな追加の仮定が、測度ゼロの集合を除いて発散を完全に排除するのである。これは1913年にルジンが予想して以来未解決だった問題に決着をつけ、20世紀解析学のもっとも深い定理の一つとみなされている。

証明

完全な証明は20世紀解析学における最も難しい論証の一つである。フェファーマンは1973年に有名な簡略化を与えており、その戦略は三つの手順で概説できる。

ステップ1(極大作用素による制御)。カールソン極大作用素 S∗f(x)=sup⁡N∣SNf(x)∣S^*f(x) = \sup_N |S_N f(x)| を L2L^2 上で評価すれば十分である。なぜなら S∗S^* に対する弱形の評価と、(収束が容易な)良い関数の稠密性を組み合わせれば、完全な概収束の主張が従うからである。

ステップ2(時間–周波数平面のタイルへの分解)。時間–周波数平面上の二進タイル——面積 ∼1\sim 1 の長方形であり、以下の『時間と周波数』の節で論じる原子そのもの——に適応した波束を用いて ff を分解する。各タイルは、時間区間と周波数帯域の両方に局在した ff の一部分を担う。

ステップ3(タイルを木にまとめて足し上げる)。カットオフ周波数 NN は xx によって変わりうるため、S∗f(x)S^*f(x) に関係するタイルは時間–周波数の包含関係で順序づけられた『木』を成す。フェファーマンの鍵となる組合せ論的評価は、すべての木が運ぶ総エネルギーを抑え、L2L^2 における S∗fS^*f を制御して証明を完成させる。

カールソンの定理をついに打ち破った手法——時間–周波数タイルを木にまとめる手法——は、現代の調和解析の他の場面でも使われる汎用的な道具であることが判明している(例えば、同じ流儀の議論が1997年にLaceyとThieleによって、双線形ヒルベルト変換に関するカルデロンの未解決問題を解決した)。

大学実数直線上のフーリエ変換

定義: フーリエ変換

実数直線 R\mathbb{R} 上の可積分関数 ff に対して、そのフーリエ変換 f^\hat f は ff の中に各純粋周波数 ξ\xi がどれだけ含まれているかを記録する。定義は f^(ξ)=∫−∞∞f(x) e−2πixξ dx\hat f(\xi) = \int_{-\infty}^{\infty} f(x)\, e^{-2\pi i x \xi}\,dx であり、指数因子 e−2πixξe^{-2\pi i x \xi} は単位長さで回転するベクトルである。 f(x)f(x) にそれを掛けて積分することで、 ff とその回転速度との相関を測る。

f^(ξ)=∫−∞∞f(x) e−2πixξ dx\hat f(\xi) = \int_{-\infty}^{\infty} f(x)\, e^{-2\pi i x \xi}\,dx

f^\hat f 自身も可積分であるとき、元の信号はそれらの純粋周波数をすべて足し戻すことで正確に復元できる: f(x)=∫−∞∞f^(ξ) e2πixξ dξf(x) = \int_{-\infty}^{\infty} \hat f(\xi)\, e^{2\pi i x \xi}\,d\xi。この反転公式こそが『関数はその周波数成分の和に等しい』という言明の厳密な意味である。

f(x)=∫−∞∞f^(ξ) e2πixξ dξf(x) = \int_{-\infty}^{\infty} \hat f(\xi)\, e^{2\pi i x \xi}\,d\xi
同じ発想の三つの側面
設定信号の定義域周波数側主要な等式
フーリエ級数円周・周期整数 Z\mathbb{Z}調和成分の離散和
フーリエ変換実数直線 R\mathbb{R}実数直線 R\mathbb{R}f(x)=∫−∞∞f^(ξ) e2πixξ dξf(x) = \int_{-\infty}^{\infty} \hat f(\xi)\, e^{2\pi i x \xi}\,d\xi
ポアソン和公式実数直線、 Z\mathbb{Z} で標本化整数 Z\mathbb{Z}∑n=−∞∞f(n)=∑k=−∞∞f^(k)\sum_{n=-\infty}^{\infty} f(n) = \sum_{k=-\infty}^{\infty} \hat f(k)

発展主要な定理

もし f∈L1(R)∩L2(R)f \in L^1(\mathbb{R}) \cap L^2(\mathbb{R}) ならば、 ∫−∞∞∣f(x)∣2 dx=∫−∞∞∣f^(ξ)∣2 dξ\int_{-\infty}^{\infty} |f(x)|^2\,dx = \int_{-\infty}^{\infty} |\hat f(\xi)|^2\,d\xi が成り立つ。言い換えれば、この正規化のもとではフーリエ変換は定数倍すら必要とせずに全エネルギーを保存する。

なぜ正しいのか?

物理的には、 ∫∣f∣2\int |f|^2 は信号のエネルギー(電圧波形が消費する電力だと考えればよい)である。プランシュレルの定理は、そのエネルギーを時間軸で信号を走査して計算しても、周波数軸でそのスペクトルを走査して計算しても同じであることを述べる——スペクトラムアナライザとオシロスコープは全電力について一致しなければならない。

証明

ステップ1(特別な場合)。まずガウス関数 f(x)=e−πx2f(x) = e^{-\pi x^2} についてこの等式を確認する。指数部の平方完成による直接計算から f^(ξ)=e−πξ2\hat f(\xi) = e^{-\pi \xi^2} が得られ、 ∫−∞∞∣f(x)∣2 dx=∫−∞∞∣f^(ξ)∣2 dξ\int_{-\infty}^{\infty} |f(x)|^2\,dx = \int_{-\infty}^{\infty} |\hat f(\xi)|^2\,d\xi の両辺は同じガウス積分 ∫e−2πx2 dx\int e^{-2\pi x^2}\,dx に等しくなる。

ステップ2(線形性と畳み込み定理による拡張)。良い性質を持つ(シュワルツ)関数 ff に対して、 f~(x)=f(−x)‾\tilde f(x) = \overline{f(-x)} として g=f∗f~g = f * \tilde f と置く。すると g(0)=∫∣f(x)∣2 dxg(0) = \int |f(x)|^2\,dx であり、畳み込み定理 f∗g^=f^⋅g^\widehat{f * g} = \hat f \cdot \hat g とフーリエ反転公式により g(0)=∫g^(ξ) dξ=∫∣f^(ξ)∣2 dξg(0) = \int \hat g(\xi)\,d\xi = \int |\hat f(\xi)|^2\,d\xi となる。これによりガウス関数を含むすべての滑らかで急減少なシュワルツ関数について等式が証明される。

ステップ3(稠密性の議論)。シュワルツ関数は L2(R)L^2(\mathbb{R}) の中で稠密である: 任意の f∈L2f \in L^2 は L2L^2 ノルムに関してシュワルツ関数の列 fn→ff_n \to f の極限として得られる。ステップ2によりフーリエ変換はこの稠密な部分空間上で等長写像であるから、それは ∫−∞∞∣f(x)∣2 dx=∫−∞∞∣f^(ξ)∣2 dξ\int_{-\infty}^{\infty} |f(x)|^2\,dx = \int_{-\infty}^{\infty} |\hat f(\xi)|^2\,d\xi を満たしたまま L2(R)L^2(\mathbb{R}) 全体への有界作用素として一意に拡張される——この拡張こそが、定義積分 f^(ξ)=∫−∞∞f(x) e−2πixξ dx\hat f(\xi) = \int_{-\infty}^{\infty} f(x)\, e^{-2\pi i x \xi}\,dx が絶対収束しないことがある L2L^2 関数の『フーリエ変換』の意味である。

シュワルツ関数 ff に対して、 ∑n=−∞∞f(n)=∑k=−∞∞f^(k)\sum_{n=-\infty}^{\infty} f(n) = \sum_{k=-\infty}^{\infty} \hat f(k) が成り立つ:すべての整数にわたる関数の総和は、そのフーリエ変換のすべての整数にわたる総和に等しい。

なぜ正しいのか?

これは信号を整数点で標本化することと、そのスペクトルを周期化することとの間の橋渡しである:左辺は標本 f(n)f(n) を足し合わせた結果であり、右辺はスペクトルから得られる結果である。この一つの等式が標本化定理、結晶学における格子和、そしてテータ関数やリーマンゼータ関数の関数等式の基礎となっている。

証明

ステップ1(周期化)。シュワルツ関数 ff に対して F(x)=∑n=−∞∞f(x+n)F(x) = \sum_{n=-\infty}^{\infty} f(x+n) と定める。 ff の急減少性によりこの和は絶対かつ一様に収束し、 FF は周期1の滑らかな関数となるので、円周上で独自のフーリエ級数を持つ。

ステップ2(周期化されたもののフーリエ係数を計算する)。 FF の第 kk フーリエ係数は ck=∫01F(x)e−2πikx dx=∑n∫01f(x+n)e−2πikx dxc_k = \int_0^1 F(x) e^{-2\pi i k x}\,dx = \sum_{n} \int_0^1 f(x+n) e^{-2\pi i k x}\,dx である。各項で y=x+ny = x+n と置換し、 e−2πikxe^{-2\pi i k x} の nn に関する周期性を使うと、断片は R\mathbb{R} 全体にわたる一つの積分にまとまる: ck=∫−∞∞f(y)e−2πiky dy=f^(k)c_k = \int_{-\infty}^{\infty} f(y) e^{-2\pi i k y}\,dy = \hat f(k)。

ステップ3( x=0 で評価する)。 FF は滑らかなので、そのフーリエ級数は各点でそれ自身に収束し、特に x=0x=0 では F(0)=∑kcke0F(0) = \sum_k c_k e^{0}、すなわち ∑n=−∞∞f(n)=∑k=−∞∞f^(k)\sum_{n=-\infty}^{\infty} f(n) = \sum_{k=-\infty}^{\infty} \hat f(k) が成り立つ。左辺は定義よりまさに F(0)=∑nf(n)F(0) = \sum_n f(n) であり、これで証明が完了する。

発展ヒルベルト変換:最初の特異積分

ポアソン核やフェイエール核のような、単に有界で可積分な核に対する畳み込みはおとなしい。解析学における最も重要な作用素のいくつかは、それとは違い、原点に本物の特異性を持つ核——y=xy=x の近くで可積分ではなく、その正の部分と負の部分の間の繊細な打ち消し合いによってのみ救われる核——に対して畳み込みを行う。その原型が実数直線上のヒルベルト変換 Hf(x)=p.v. 1π∫−∞∞f(y)x−y dyHf(x) = \text{p.v.}\,\frac{1}{\pi}\int_{-\infty}^{\infty}\frac{f(y)}{x-y}\,dy であり、主値として理解される: y=xy=x における特異点は対称に取り除かれてから極限が取られる。周波数側では、この特異な畳み込みは驚くほど単純な乗算になる: (Hf^)(ξ)=−i sgn(ξ) f^(ξ)(\widehat{Hf})(\xi) = -i\,\mathrm{sgn}(\xi)\,\hat f(\xi)——ヒルベルト変換はあらゆる周波数に対する純粋な 90∘90^\circ の位相回転であり、大きさを変えずに ii の符号を反転させる。

Hf(x)=p.v. 1π∫−∞∞f(y)x−y dyHf(x) = \text{p.v.}\,\frac{1}{\pi}\int_{-\infty}^{\infty}\frac{f(y)}{x-y}\,dy
(Hf^)(ξ)=−i sgn(ξ) f^(ξ)(\widehat{Hf})(\xi) = -i\,\mathrm{sgn}(\xi)\,\hat f(\xi)

これを見るもう一つの古典的な方法がある:実数直線上で定義された ff に対して、 u=Py∗fu = P_y * f をその上半平面への調和拡張(前節で扱ったポアソン核との畳み込み)とし、 v=Qy∗fv = Q_y * f を、共役ポアソン核 QyQ_y から作られるその調和共役として、F=u+ivF=u+iv が x+iyx+iy の正則関数となるようにする。 y→0+y\to 0^+ のとき、 vv の境界値はまさに ff のヒルベルト変換を復元する。これはこの理論全体の複素解析的な祖先である:境界上の特異積分作用素とは、実は一つ上の次元に住む正真正銘の正則関数の影にほかならない。

ヒルベルト変換 HH は、すべての 1<p<∞1<p<\infty に対して Lp(R)L^p(\mathbb{R}) 上で有界である:あらゆる ff について ∥Hf∥Lp≤Cp∥f∥Lp\|Hf\|_{L^p} \le C_p\|f\|_{L^p} を満たす定数 CpC_p が存在する。(マルセル・リース、1927年)

なぜ正しいのか?

L2L^2 の場合は実質ただで得られる:乗数 −i sgn(ξ)-i\,\mathrm{sgn}(\xi) はほとんどいたるところで絶対値がちょうど 11 であるため、プランシュレルの定理により HH は文字通り L2L^2 上の等長写像である——あらゆる周波数の 90∘90^\circ 位相回転は全エネルギーを何も変えない。しかし有界性を他の pp に拡張することは本質的に別の、はるかに難しい問題である:L2L^2 の外にはパーセバルの等式の類似物が存在しないため、代数的な近道は使えず、証明は HfHf が大きくなる集合の大きさと幾何に関する実変数的な評価に頼らざるを得ない。

証明

ステップ1(L2L^2、直接)。プランシュレルの定理と乗数の公式により、ξ≠0\xi\ne0 で ∣sgn(ξ)∣=1|\mathrm{sgn}(\xi)|=1 であることを用いて ∥Hf∥L22=∫∣Hf^(ξ)∣2 dξ=∫∣sgn(ξ)∣2∣f^(ξ)∣2 dξ=∥f∥L22\|Hf\|_{L^2}^2 = \int |\widehat{Hf}(\xi)|^2\,d\xi = \int |\mathrm{sgn}(\xi)|^2|\hat f(\xi)|^2\,d\xi = \|f\|_{L^2}^2 が成り立つ(一点は積分に影響しない)。したがって HH は L2(R)L^2(\mathbb{R}) 上の等長写像であり、特に C2=1C_2=1 で有界である。

ステップ2(一般の 1<p<∞1<p<\infty、後方参照)。指数の全範囲は、ステップ1からいかなる易しい議論によっても従わない。これは実変数コルデロン–ジグムント理論の正真正銘の定理である:ヒルベルト変換の核 1/(πx)1/(\pi x) は、コルデロン–ジグムント核の滑らかさと打ち消し条件をまさに満たしており、そのため次節で展開される一般的な特異積分の機構がそのまま適用できる。その機構は、ステップ1で証明した L2L^2 の評価から弱形 (1,1)(1,1) の評価を生み出し、弱-(1,1)(1,1) 評価と L2L^2 評価の間のマルツィンキェヴィッチ補間(そして p>2p>2 に対する双対性の議論)によって、すべての 1<p<∞1<p<\infty について証明が完成する——この議論の後半部分は、ここで繰り返す代わりに、すぐ下のコルデロン–ジグムントの節で完全に実行される。

例: 指示関数のヒルベルト変換

[−1,1][-1,1] の指示関数 f=1[−1,1]f = \mathbf{1}_{[-1,1]} を考える:有界でコンパクトな台を持つ、まったく取り立てて特徴のない関数である。HfHf を求めよ。

解答

定義より Hf(x)=1π p.v.∫−11dyx−yHf(x) = \frac{1}{\pi}\,\text{p.v.}\int_{-1}^{1} \frac{dy}{x-y} である。[−1,1][-1,1] の外側では取り除くべき特異点はなく、1/(x−y)1/(x-y) の yy に関する原始関数は −log⁡∣x−y∣-\log|x-y| であるから、∫−11dyx−y=[−log⁡∣x−y∣]y=−1y=1=log⁡∣x+1x−1∣\int_{-1}^{1} \frac{dy}{x-y} = \big[-\log|x-y|\big]_{y=-1}^{y=1} = \log\left|\frac{x+1}{x-1}\right| となる((−1,1)(-1,1) の内部の xx についても、y=x−y=x^- と y=x+y=x^+ における2つの発散が打ち消し合うため、同じ計算が正真正銘の主値として成り立つ)。これより Hf(x)=1πlog⁡∣x+1x−1∣Hf(x) = \frac{1}{\pi}\log\left|\frac{x+1}{x-1}\right| が得られる。

際立った特徴は x=±1x=\pm1 で起こることである:ff 自身はそこで完全に有界である(単に 11 から 00 へ跳ぶだけ)にもかかわらず、HfHf はまさにその2点で対数的に発散する。有界でコンパクトな台を持つ、教科書的に素性の良い関数が、HH によって2つの正真正銘の特異点を持つ関数に変えられてしまう——これがあらゆる特異積分作用素に共通する局所的な特徴であり、まさにこれゆえに HH を制御することは、素性の良い有界な核との畳み込みを制御するよりもはるかに繊細な評価を必要とするのである。

発展ハーディ–リトルウッド極大関数とコルデロン–ジグムント分解

ヒルベルト変換だけでなく、一般の特異積分作用素についてM.リースの定理を完成させるには、解析学者たちはフーリエ変換について何も語らない、純粋に実変数的な道具を必要とした。ハーディとリトルウッドは1930年にそれを導入した:極大関数 Mf(x)=sup⁡r>012r∫∣y−x∣<r∣f(y)∣ dyMf(x) = \sup_{r>0} \frac{1}{2r}\int_{|y-x|<r} |f(y)|\,dy であり、これは xx を中心とする任意の区間にわたる ∣f∣|f| のありうる最大の平均値である。これは、xx の近くで ff に対して行われうるあらゆる平均化操作を、一度にまとめて制御するための、無骨だが極めて有効な道具である。

MM は弱形 (1,1)(1,1) である:普遍定数 CC(常に C=5C=5 ととれる)に対して ∣{x:Mf(x)>λ}∣≤Cλ∥f∥1|\{x : Mf(x) > \lambda\}| \le \frac{C}{\lambda}\|f\|_1 が成り立つ。したがって補間により、MM はすべての 1<p≤∞1<p\le\infty に対して Lp(R)L^p(\mathbb{R}) 上でも有界である——しかし MM 自身は L1(R)L^1(\mathbb{R}) 上では決して有界にならない。

なぜ正しいのか?

Mf(x)Mf(x) は、xx の周りの区間上でとりうるあらゆる ff の平均——あらゆる尺度での移動平均——を一挙に制御する。この一つの量を制御することが、まさにルベーグ微分定理(縮小していく xx の周りの区間にわたる ff の平均が、ほとんどすべての xx で f(x)f(x) に収束する)を証明するために必要なことであることが判明する:MfMf がほとんどいたるところ有限であるとわかれば、短くて易しい議論によって完全な微分の主張へと格上げできる。これは、易しい L2L^2 評価に対してプランシュレルの等式が果たした役割を、以下のコルデロン–ジグムント理論に対して果たす、実変数的な原動力である。

証明

ステップ1(良い球による被覆)。λ>0\lambda>0 を固定する。Mf(x)>λMf(x)>\lambda を満たす各 xx に対し、上限の定義より 1∣Bx∣∫Bx∣f∣>λ\frac{1}{|B_x|}\int_{B_x}|f| > \lambda を満たす xx を中心とする区間 BxB_x が存在する。これらの球 {Bx}\{B_x\} は集合 {Mf>λ}\{Mf>\lambda\} を被覆する。

ステップ2(ヴィタリの 5r5r 被覆補題)。半径が有界な球の任意の族から、可算な互いに素な部分族 {Bi}\{B_i\} を常に取り出すことができ、55 倍に拡大した球 {5Bi}\{5B_i\} が元の族全体の和集合を依然として被覆するようにできる。これを {Bx}\{B_x\} に適用して、{Mf>λ}⊆⋃i5Bi\{Mf>\lambda\}\subseteq\bigcup_i 5B_i を満たす互いに素な部分族 {Bi}\{B_i\} を得る。

ステップ3(互いに素な部分を足し合わせる)。各 BiB_i は構成により λ∣Bi∣<∫Bi∣f∣\lambda|B_i| < \int_{B_i}|f| を満たす。BiB_i は互いに素であるから、ii について足し合わせると λ∑i∣Bi∣<∑i∫Bi∣f∣≤∫R∣f∣=∥f∥1\lambda\sum_i|B_i| < \sum_i\int_{B_i}|f| \le \int_{\mathbb{R}}|f| = \|f\|_1 となり、∑i∣Bi∣<∥f∥1/λ\sum_i|B_i| < \|f\|_1/\lambda が得られる。最後に ∣{Mf>λ}∣≤∑i∣5Bi∣=5∑i∣Bi∣<5λ∥f∥1|\{Mf>\lambda\}| \le \sum_i|5B_i| = 5\sum_i|B_i| < \dfrac{5}{\lambda}\|f\|_1 となり、これはまさに ∣{x:Mf(x)>λ}∣≤5λ∥f∥1|\{x : Mf(x) > \lambda\}| \le \frac{5}{\lambda}\|f\|_1 である。

定義: コルデロン–ジグムント分解

f∈L1(R)f\in L^1(\mathbb{R}) とある高さ α>0\alpha>0 が与えられたとき、コルデロンとジグムント(1952年)は、直線を ff と α\alpha に適応した互いに素な『停止』二進区間に分割する方法を示した:一つの大きな二進区間から出発し、それを再帰的に二等分していき、二進区間 QjQ_j の平均 1∣Qj∣∫Qj∣f∣\frac{1}{|Q_j|}\int_{Q_j}|f| が α\alpha を超えたまさに最初の時点でそれを停止させ保持する(その親の平均はまだ ≤α\le\alpha だったので、二等分してもたかだか2倍にしかならない)。これにより、互いに素な停止区間の族 {Qj}\{Q_j\} が生まれ、α<1∣Qj∣∫Qj∣f∣≤2α\alpha < \frac{1}{|Q_j|}\int_{Q_j}|f| \le 2\alpha を満たし、全長は ∑∣Qj∣≤∥f∥1/α\sum|Q_j|\le \|f\|_1/\alpha であり、⋃jQj\bigcup_j Q_j の外側に残る部分ではほとんどいたるところ ∣f∣≤α|f|\le\alpha となる。

[0,1] 上の固定された関数のコルデロン–ジグムント分解を示すインタラクティブなグラフで、調整可能な水準スライダー α により、停止区間、その平均値、および α と 2α における基準線が表示される。
[0,1] 上の固定されたでこぼこな |f|(x≈0.18 付近に細い尖った峰、x≈0.62 付近にほぼ特異的な盛り上がり、x≈0.86 付近に幅広い盛り上がりを持つ)を示し、この水準でのコルデロン–ジグムント停止区間が赤く塗られている:塗られた各区間が一つの停止立方体 Q_j であり、その上に引かれた太い赤線分がその平均値を示す——それは常に破線で示された α と 2α の基準線の間に収まる。スライダーを動かして α を変えると、停止区間が縮んだり、広がったり、分裂したり、併合したりする様子が見える。

この分解は ff を『良い』部分と無限個の『悪い』断片に分ける、f=g+∑jbjf = g + \sum_j b_j:良い部分 gg は、すべての QjQ_j の外側では ff に等しく、保持された各 QjQ_j の上では(だいたい α\alpha 程度の大きさの)ff の平均値に等しいので、いたるところ ∣g∣≤2α|g|\le 2\alpha であり、gg は易しい L2L^2 理論で扱える。各悪い断片 bjb_j はそれ自身の立方体 QjQ_j 上に台を持ち、そこで平均ゼロである、∫Qjbj dx=0\int_{Q_j} b_j\,dx = 0——この打ち消し合いこそが、QjQ_j から十分離れれば滑らかな核による bjb_j の寄与が小さいままでいられる理由である。なぜなら核は QjQ_j の全体にわたってほぼ一定に見え、ほぼ一定の核は平均ゼロの関数を積分するとほとんど何も残らないからである。この良い–悪いの分割を、弱形 (1,1)(1,1) のハーディ–リトルウッド極大評価と組み合わせ、マルツィンキェヴィッチ補間を通すことこそが、コルデロンとジグムントがヒルベルト変換だけでなく一般の特異積分作用素についてM.リースの定理を完成させ、LpL^p の有界性をすべての 1<p<∞1<p<\infty に拡張するために使った実変数的な機構であり、これによって上の定理へと話がつながって輪が閉じる。

大学実世界での応用と具体例

フーリエ変換は『時間的な形』を『周波数的な内容』に変えるため、技術者や科学者が振動を分離・フィルタリング・計数する必要のあるあらゆる場面で標準的な道具となっている:オーディオイコライザー、MRIや電波望遠鏡のイメージング、そして物理学の熱方程式や波動方程式などである。

例: 録音からヒス雑音を除去する

マイクが f(t)=cos⁡(2π⋅440 t)+cos⁡(2π⋅3000 t)f(t) = \cos(2\pi \cdot 440\,t) + \cos(2\pi \cdot 3000\,t) を記録した:440 Hz の音符(A4)と 3000 Hz の電子的なヒス雑音が混じっている。技術者は ∣ξ∣>1000|\xi| > 1000 のときは常に f^(ξ)=0\hat f(\xi) = 0 とする理想的なローパスフィルタを適用する。出力される信号は何か。

解答

各コサインは周波数 ±ν\pm\nu における2つの純粋な回転指数関数の和である:周波数領域では、 f(t)=cos⁡(2π⋅440 t)+cos⁡(2π⋅3000 t)f(t) = \cos(2\pi \cdot 440\,t) + \cos(2\pi \cdot 3000\,t) のエネルギーは ξ=±440\xi = \pm 440 と ξ=±3000\xi = \pm 3000 にのみ集中している( f^\hat f 中のスパイクとして)。

フィルタは ∣ξ∣≤1000|\xi| \le 1000 を満たすものをすべて保持し、それ以外を除去する。 440≤1000440 \le 1000 だが 3000>10003000 > 1000 なので、 ±3000\pm 3000 のスパイクは除去され、 ±440\pm 440 のスパイクはそのまま残る。

残ったものに反転公式を適用すると、生き残った周波数が正確に再構成される:出力は cos⁡(2π⋅440 t)\cos(2\pi \cdot 440\,t) であり、ヒス雑音の消えたきれいな音符である。これはグラフィックイコライザーのローパスつまみが実際に行っていることそのものである。

例: 熱が広がる理由:熱方程式を解く

ある棒の初期温度分布が u(x,0)=f(x)u(x,0) = f(x) であり、熱方程式 ut=uxxu_t = u_{xx} に従う。フーリエ変換( xx について)を用いて、 t>0t>0 における u(x,t)u(x,t) を求めよ。

解答

xx についてフーリエ変換をとると、各空間微分 ∂x\partial_x は 2πiξ2\pi i \xi による乗算になるので、 uxxu_{xx} は −4π2ξ2u^-4\pi^2\xi^2 \hat u になる。方程式 ut=uxxu_t = u_{xx} は、固定した各 ξ\xi に対する tt についての常微分方程式 ∂tu^(ξ,t)=−4π2ξ2 u^(ξ,t)\partial_t \hat u(\xi,t) = -4\pi^2\xi^2\, \hat u(\xi,t) になる——難しい偏微分方程式が、互いに独立な常微分方程式の族という扱いやすいものになる大きな簡略化である。

この常微分方程式の解は u^(ξ,t)=f^(ξ) e−4π2ξ2t\hat u(\xi,t) = \hat f(\xi)\, e^{-4\pi^2 \xi^2 t} であり、高周波数 ξ\xi を急速に消し去る指数減衰因子である:細かい空間的な詳細は急速に滑らかになり、これは鋭い温度スパイクが最初に平らになるという日常的事実と一致する。

変換を逆にすると(畳み込み定理を逆向きに使えば、周波数領域での積は空間領域での畳み込みになる) u(x,t)=(f∗Kt)(x),  Kt(x)=14πt e−x2/(4t)u(x,t) = (f * K_t)(x),\ \ K_t(x) = \dfrac{1}{\sqrt{4\pi t}}\, e^{-x^2/(4t)} が得られる:時刻 tt における温度は、初期分布が広がるガウス型の山 KtK_t によってぼかされた(畳み込まれた)ものである——これはまさに熱が拡散する様子の数学的な描像である。

定義: 不確定性原理

関数は時間と周波数の両方で同時に鋭く局在することはできない。 ∣f∣2|f|^2 と ∣f^∣2|\hat f|^2 の広がりをそれぞれ Δx\Delta x と Δξ\Delta \xi で測ると、 Δx⋅Δξ≥14π\Delta x \cdot \Delta \xi \ge \dfrac{1}{4\pi} が成り立つ。パルスを時間的に押しつぶす(小さい Δx\Delta x)と、そのスペクトルは広がらざるを得ず(大きい Δξ\Delta \xi)、逆もまた然りである——これは量子力学におけるハイゼンベルクの不確定性原理を導くのと同じトレードオフである。

Δx⋅Δξ≥14π\Delta x \cdot \Delta \xi \ge \dfrac{1}{4\pi}

規約 f^(ξ)=∫f(x)e−2πixξ dx\hat f(\xi) = \int f(x) e^{-2\pi i x\xi}\,dx のもとで、ガウス関数 f(x)=e−πx2f(x) = e^{-\pi x^2} のフーリエ変換は:

プランシュレルの定理 ∫−∞∞∣f(x)∣2 dx=∫−∞∞∣f^(ξ)∣2 dξ\int_{-\infty}^{\infty} |f(x)|^2\,dx = \int_{-\infty}^{\infty} |\hat f(\xi)|^2\,d\xi を最もよく表す記述は:

MRI装置は空間信号のフーリエ変換の標本(『kスペース』)を測定する。低周波数の標本( ∣ξ∣|\xi| が小さいもの)のみを収集した場合、再構成される画像は:

不確定性原理 Δx⋅Δξ≥14π\Delta x \cdot \Delta \xi \ge \dfrac{1}{4\pi} が意味するのは:

フェイエールの定理は、チェザロ平均 σNf\sigma_N f が任意の連続関数 ff に一様収束することを述べている。カールソンの定理はこれとは別の、はるかに深い、次のことについての主張である:

M.リースの定理により、ヒルベルト変換 HH が Lp(R)L^p(\mathbb{R}) 上で有界となるのはちょうど次の場合である:

参考文献

  1. Terence Tao (2003). Recent progress on the restriction conjecture · arXiv:math/0311181
  2. Elias M. Stein, Guido Weiss (1971). Introduction to Fourier Analysis on Euclidean Spaces
  3. Loukas Grafakos (2014). Classical Fourier Analysis · DOI:10.1007/978-1-4939-1194-3