← 戻る ライブラリ › 解析学 › 現代解析学 解析学
調和解析学 関数や信号を基本波に分解する分野で、フーリエ級数をより広い舞台に一般化したもの。
直観 波から周波数へ 私たちが耳にする音——バイオリンの音、声、交通騒音——は実際には時間とともに変化する一つの圧力信号にすぎない。しかし耳(やグラフィックイコライザー)は、それが異なる音高を持つ多くの純音からできていることを教えてくれる。調和解析とは、その分解を扱う数学であり、関数 f f f を、それぞれ固有の周波数 ξ \xi ξ を持つ純粋な振動 e 2 π i x ξ e^{2\pi i x \xi} e 2 π i x ξ の和や積分として書き直す。フーリエ級数は周期信号に対して離散的な周波数の列でこれを行うが、調和解析はこの発想を非周期信号、高次元、さらには他の群(円周、有限群、リー群)へと一般化する。
最初の24個のフーリエ調和成分から再構成された矩形波に近い波形。n を動かすと、ジャンプ付近のギブス現象の『リンギング』が幅は狭くなるが高さは減らないことが見える。 大学 円周上のフーリエ級数:収束の問題 周期関数 f f f のフーリエ級数の第 N N N 部分和を S N f ( x ) = ∑ ∣ n ∣ ≤ N f ^ ( n ) e i n x S_N f(x) = \sum_{|n|\le N} \hat f(n)\, e^{inx} S N f ( x ) = ∑ ∣ n ∣ ≤ N f ^ ( n ) e in x と書く。ここで f ^ ( n ) = 1 2 π ∫ − π π f ( x ) e − i n x d x \hat f(n) = \frac{1}{2\pi}\int_{-\pi}^{\pi} f(x)\, e^{-inx}\,dx f ^ ( n ) = 2 π 1 ∫ − π π f ( x ) e − in x d x である。L 2 L^2 L 2 では万事がすっきりしている:指数関数 e i n x e^{inx} e in x は正規直交基底をなし、パーセバルの等式 1 2 π ∫ ∣ f ∣ 2 = ∑ ∣ f ^ ( n ) ∣ 2 \frac{1}{2\pi}\int |f|^2 = \sum |\hat f(n)|^2 2 π 1 ∫ ∣ f ∣ 2 = ∑ ∣ f ^ ( n ) ∣ 2 が成り立ち、S N f → f S_N f \to f S N f → f が L 2 L^2 L 2 ノルムで成り立つ。難しいのは各点収束である:ある特定の点 x x x で S N f ( x ) → f ( x ) S_N f(x) \to f(x) S N f ( x ) → f ( x ) となるだろうか。上で示した矩形波はすでにその難しさを示唆している——ジャンプの近くでは、N N N がどれほど大きくても S N f S_N f S N f は常にジャンプの高さの約9%だけオーバーシュートする(ギブス現象であり、オーバーシュート量は半ジャンプの 2 π S i ( π ) ≈ 1.179 \frac{2}{\pi}\mathrm{Si}(\pi) \approx 1.179 π 2 Si ( π ) ≈ 1.179 倍に近づく)。
定義: ディリクレ核
部分和はそれ自体が畳み込みである: S N f ( x ) = ( D N ∗ f ) ( x ) = 1 2 π ∫ − π π D N ( x − y ) f ( y ) d y S_N f(x) = (D_N * f)(x) = \frac{1}{2\pi}\int_{-\pi}^{\pi} D_N(x-y)\, f(y)\, dy S N f ( x ) = ( D N ∗ f ) ( x ) = 2 π 1 ∫ − π π D N ( x − y ) f ( y ) d y 。ここでディリクレ核は D N ( x ) = ∑ ∣ n ∣ ≤ N e i n x = sin ( ( N + 1 2 ) x ) sin ( x / 2 ) D_N(x) = \sum_{|n|\le N} e^{inx} = \dfrac{\sin\left((N+\tfrac12)x\right)}{\sin(x/2)} D N ( x ) = ∑ ∣ n ∣ ≤ N e in x = sin ( x /2 ) sin ( ( N + 2 1 ) x ) である。したがって収束に関するあらゆる問いは、実質的に D N D_N D N の形状についての問いになる。
D N ( x ) = ∑ ∣ n ∣ ≤ N e i n x = sin ( ( N + 1 2 ) x ) sin ( x / 2 ) D_N(x) = \sum_{|n|\le N} e^{inx} = \dfrac{\sin\left((N+\tfrac12)x\right)}{\sin(x/2)} D N ( x ) = ∣ n ∣ ≤ N ∑ e in x = sin ( x /2 ) sin ( ( N + 2 1 ) x ) N = 9 N=9 N = 9 のディリクレ核 D N D_N D N :フェイエール核やポアソン核と異なり、ゼロより下に落ち込む部分がある(サンゴ色の部分)。その L 1 L^1 L 1 ノルム——ルベーグ定数——は有界のままではなく ( 4 / π 2 ) log N (4/\pi^2)\log N ( 4/ π 2 ) log N のように増大する。核を切り替えたり N N N を大きくしたりして比べてみよう。フェイエール核 F N F_N F N (D 0 , … , D N D_0,\ldots,D_N D 0 , … , D N のチェザロ平均)とポアソン核 P r P_r P r はどちらも非負なので、正真正銘の近似単位元 である:質量 1 1 1 、L 1 L^1 L 1 ノルム 1 1 1 、そして N → ∞ N\to\infty N → ∞ (あるいは r → 1 r\to1 r → 1 )のとき質量が 0 0 0 に集中する——これだけで、任意の連続関数 f f f に対して F N ∗ f → f F_N * f\to f F N ∗ f → f と P r ∗ f → f P_r * f\to f P r ∗ f → f が一様に成り立つことが保証される。ディリクレ核はそうではない:その負の部分のせいで 1 2 π ∫ ∣ D N ∣ \frac{1}{2\pi}\int|D_N| 2 π 1 ∫ ∣ D N ∣ (ルベーグ定数 )は 4 π 2 log N \frac{4}{\pi^2}\log N π 2 4 log N のように増大し、N → ∞ N\to\infty N → ∞ で有界でなくなる。バナッハ–シュタインハウスの一様有界性原理により、ルベーグ定数が有界でないことは、ある点でフーリエ級数が発散する連続関数の存在を強制する——これはまさに du Bois-Reymond が1873年に手作業で構成した現象であり、ここでは労せずして得られる。コルモゴロフは後に、単に可積分なだけの関数についてはさらに悪いことが起こりうることを示した:1923年には概発散する L 1 L^1 L 1 関数を構成し、1926年には至る所 発散するものを構成した。
f f f が連続で周期 2 π 2\pi 2 π をもつならば、チェザロ平均 σ N f = 1 N + 1 ∑ k = 0 N S k f = F N ∗ f \sigma_N f = \frac{1}{N+1}\sum_{k=0}^N S_k f = F_N * f σ N f = N + 1 1 ∑ k = 0 N S k f = F N ∗ f は N → ∞ N\to\infty N → ∞ のとき f f f に一様収束する。
なぜ正しいのか? 部分和を平均することで、まさに S N f S_N f S N f を暴れさせている振動が打ち消される:D N D_N D N には激しい負の部分があるのに対し、その移動平均 F N F_N F N は滑らかで非負の山であり、常に f f f の各部分を正の重みで足し合わせるだけである。そのため S N f S_N f S N f 自体は収束しないことがあっても、その移動平均は必ず落ち着く——これは微積分において振動する部分和の列をチェザロ平均によって手なずけるのと同じ技である。
証明 ステップ1(フェイエール核は正の近似単位元である)。D 0 , … , D N D_0,\ldots,D_N D 0 , … , D N を平均し、望遠鏡型の三角関数恒等式を用いると、閉じた形 F N ( x ) = 1 N + 1 ( sin ( N + 1 2 x ) sin ( x / 2 ) ) 2 ≥ 0 F_N(x) = \frac{1}{N+1}\left(\dfrac{\sin\left(\tfrac{N+1}{2}x\right)}{\sin(x/2)}\right)^2 \ge 0 F N ( x ) = N + 1 1 ( sin ( x /2 ) sin ( 2 N + 1 x ) ) 2 ≥ 0 が得られる。F N = 1 N + 1 ∑ k = 0 N D k F_N = \frac{1}{N+1}\sum_{k=0}^N D_k F N = N + 1 1 ∑ k = 0 N D k かつ任意の k k k について 1 2 π ∫ D k = 1 \frac{1}{2\pi}\int D_k = 1 2 π 1 ∫ D k = 1 なので、1 2 π ∫ F N = 1 \frac{1}{2\pi}\int F_N = 1 2 π 1 ∫ F N = 1 でもある。非負であることから 1 2 π ∫ ∣ F N ∣ = 1 \frac{1}{2\pi}\int |F_N| = 1 2 π 1 ∫ ∣ F N ∣ = 1 でもあり、したがって L 1 L^1 L 1 ノルムは決して爆発しない。
ステップ2(質量が 0 0 0 に集中する)。任意に固定した δ > 0 \delta>0 δ > 0 に対し、δ ≤ ∣ x ∣ ≤ π \delta \le |x| \le \pi δ ≤ ∣ x ∣ ≤ π では分母 sin ( x / 2 ) 2 \sin(x/2)^2 sin ( x /2 ) 2 が正の定数で下から評価できるので、その範囲で一様に F N ( x ) = O ( 1 / N ) F_N(x) = O(1/N) F N ( x ) = O ( 1/ N ) となる。[ − δ , δ ] [-\delta,\delta] [ − δ , δ ] の外側の質量は N → ∞ N\to\infty N → ∞ で 0 0 0 に近づく。
ステップ3(近似単位元の評価)。1 2 π ∫ F N = 1 \frac{1}{2\pi}\int F_N = 1 2 π 1 ∫ F N = 1 を用いて σ N f ( x ) − f ( x ) = 1 2 π ∫ − π π F N ( y ) ( f ( x − y ) − f ( x ) ) d y \sigma_N f(x) - f(x) = \frac{1}{2\pi}\int_{-\pi}^{\pi} F_N(y)\big(f(x-y)-f(x)\big)\,dy σ N f ( x ) − f ( x ) = 2 π 1 ∫ − π π F N ( y ) ( f ( x − y ) − f ( x ) ) d y と書く。積分を ∣ y ∣ < δ |y|<\delta ∣ y ∣ < δ と δ ≤ ∣ y ∣ ≤ π \delta\le|y|\le\pi δ ≤ ∣ y ∣ ≤ π に分ける。前半では、f f f の一様連続性により δ \delta δ を小さくとれば ∣ f ( x − y ) − f ( x ) ∣ < ε |f(x-y)-f(x)|<\varepsilon ∣ f ( x − y ) − f ( x ) ∣ < ε となり、F N ≥ 0 F_N\ge0 F N ≥ 0 かつ全質量 1 1 1 よりその部分は ε \varepsilon ε で抑えられる。後半では F N = O ( 1 / N ) F_N=O(1/N) F N = O ( 1/ N ) かつ f f f が有界なので、その部分は N → ∞ N\to\infty N → ∞ のとき → 0 \to 0 → 0 となる。両方の評価は x x x について一様なので、σ N f → f \sigma_N f\to f σ N f → f が一様に成り立つ。
f ∈ L 2 ( T ) f\in L^2(\mathbb{T}) f ∈ L 2 ( T ) ならば、ほとんどすべての x x x で S N f ( x ) → f ( x ) S_N f(x)\to f(x) S N f ( x ) → f ( x ) となる。(Hunt は1968年にこれを p > 1 p>1 p > 1 を満たすすべての f ∈ L p ( T ) f\in L^p(\mathbb{T}) f ∈ L p ( T ) に拡張した。)
なぜ正しいのか? コルモゴロフによる、至る所 発散する1926年の L 1 L^1 L 1 関数の例を踏まえると、L 2 L^2 L 2 ——わずかに強い条件にすぎない——も同じように、あるいは少なくともほとんど至る所で失敗するだろうと考えるのはごく自然に思えた。カールソンの定理はまさに驚きである:二乗可積分というわずかな追加の仮定が、測度ゼロの集合を除いて発散を完全に排除するのである。これは1913年にルジンが予想して以来未解決だった問題に決着をつけ、20世紀解析学のもっとも深い定理の一つとみなされている。
証明 完全な証明は20世紀解析学における最も難しい論証の一つである。フェファーマンは1973年に有名な簡略化を与えており、その戦略は三つの手順で概説できる。
ステップ1(極大作用素による制御)。カールソン極大作用素 S ∗ f ( x ) = sup N ∣ S N f ( x ) ∣ S^*f(x) = \sup_N |S_N f(x)| S ∗ f ( x ) = sup N ∣ S N f ( x ) ∣ を L 2 L^2 L 2 上で評価すれば十分である。なぜなら S ∗ S^* S ∗ に対する弱形の評価と、(収束が容易な)良い関数の稠密性を組み合わせれば、完全な概収束の主張が従うからである。
ステップ2(時間–周波数平面のタイルへの分解)。時間–周波数平面上の二進タイル——面積 ∼ 1 \sim 1 ∼ 1 の長方形であり、以下の『時間と周波数』の節で論じる原子そのもの——に適応した波束を用いて f f f を分解する。各タイルは、時間区間と周波数帯域の両方に局在した f f f の一部分を担う。
ステップ3(タイルを木にまとめて足し上げる)。カットオフ周波数 N N N は x x x によって変わりうるため、S ∗ f ( x ) S^*f(x) S ∗ f ( x ) に関係するタイルは時間–周波数の包含関係で順序づけられた『木』を成す。フェファーマンの鍵となる組合せ論的評価は、すべての木が運ぶ総エネルギーを抑え、L 2 L^2 L 2 における S ∗ f S^*f S ∗ f を制御して証明を完成させる。
カールソンの定理をついに打ち破った手法——時間–周波数タイルを木にまとめる手法——は、現代の調和解析の他の場面でも使われる汎用的な道具であることが判明している(例えば、同じ流儀の議論が1997年にLaceyとThieleによって、双線形ヒルベルト変換に関するカルデロンの未解決問題を解決した)。
大学 実数直線上のフーリエ変換 定義: フーリエ変換
実数直線 R \mathbb{R} R 上の可積分関数 f f f に対して、そのフーリエ変換 f ^ \hat f f ^ は f f f の中に各純粋周波数 ξ \xi ξ がどれだけ含まれているかを記録する。定義は f ^ ( ξ ) = ∫ − ∞ ∞ f ( x ) e − 2 π i x ξ d x \hat f(\xi) = \int_{-\infty}^{\infty} f(x)\, e^{-2\pi i x \xi}\,dx f ^ ( ξ ) = ∫ − ∞ ∞ f ( x ) e − 2 π i x ξ d x であり、指数因子 e − 2 π i x ξ e^{-2\pi i x \xi} e − 2 π i x ξ は単位長さで回転するベクトルである。 f ( x ) f(x) f ( x ) にそれを掛けて積分することで、 f f f とその回転速度との相関を測る。
f ^ ( ξ ) = ∫ − ∞ ∞ f ( x ) e − 2 π i x ξ d x \hat f(\xi) = \int_{-\infty}^{\infty} f(x)\, e^{-2\pi i x \xi}\,dx f ^ ( ξ ) = ∫ − ∞ ∞ f ( x ) e − 2 π i x ξ d x f ^ \hat f f ^ 自身も可積分であるとき、元の信号はそれらの純粋周波数をすべて足し戻すことで正確に復元できる: f ( x ) = ∫ − ∞ ∞ f ^ ( ξ ) e 2 π i x ξ d ξ f(x) = \int_{-\infty}^{\infty} \hat f(\xi)\, e^{2\pi i x \xi}\,d\xi f ( x ) = ∫ − ∞ ∞ f ^ ( ξ ) e 2 π i x ξ d ξ 。この反転公式こそが『関数はその周波数成分の和に等しい』という言明の厳密な意味である。
f ( x ) = ∫ − ∞ ∞ f ^ ( ξ ) e 2 π i x ξ d ξ f(x) = \int_{-\infty}^{\infty} \hat f(\xi)\, e^{2\pi i x \xi}\,d\xi f ( x ) = ∫ − ∞ ∞ f ^ ( ξ ) e 2 π i x ξ d ξ 同じ発想の三つの側面 設定 信号の定義域 周波数側 主要な等式 フーリエ級数 円周・周期 整数 Z \mathbb{Z} Z 調和成分の離散和 フーリエ変換 実数直線 R \mathbb{R} R 実数直線 R \mathbb{R} R f ( x ) = ∫ − ∞ ∞ f ^ ( ξ ) e 2 π i x ξ d ξ f(x) = \int_{-\infty}^{\infty} \hat f(\xi)\, e^{2\pi i x \xi}\,d\xi f ( x ) = ∫ − ∞ ∞ f ^ ( ξ ) e 2 π i x ξ d ξ ポアソン和公式 実数直線、 Z \mathbb{Z} Z で標本化 整数 Z \mathbb{Z} Z ∑ n = − ∞ ∞ f ( n ) = ∑ k = − ∞ ∞ f ^ ( k ) \sum_{n=-\infty}^{\infty} f(n) = \sum_{k=-\infty}^{\infty} \hat f(k) ∑ n = − ∞ ∞ f ( n ) = ∑ k = − ∞ ∞ f ^ ( k )
発展 主要な定理 もし f ∈ L 1 ( R ) ∩ L 2 ( R ) f \in L^1(\mathbb{R}) \cap L^2(\mathbb{R}) f ∈ L 1 ( R ) ∩ L 2 ( R ) ならば、 ∫ − ∞ ∞ ∣ f ( x ) ∣ 2 d x = ∫ − ∞ ∞ ∣ f ^ ( ξ ) ∣ 2 d ξ \int_{-\infty}^{\infty} |f(x)|^2\,dx = \int_{-\infty}^{\infty} |\hat f(\xi)|^2\,d\xi ∫ − ∞ ∞ ∣ f ( x ) ∣ 2 d x = ∫ − ∞ ∞ ∣ f ^ ( ξ ) ∣ 2 d ξ が成り立つ。言い換えれば、この正規化のもとではフーリエ変換は定数倍すら必要とせずに全エネルギーを保存する。
なぜ正しいのか? 物理的には、 ∫ ∣ f ∣ 2 \int |f|^2 ∫ ∣ f ∣ 2 は信号のエネルギー(電圧波形が消費する電力だと考えればよい)である。プランシュレルの定理は、そのエネルギーを時間軸で信号を走査して計算しても、周波数軸でそのスペクトルを走査して計算しても同じであることを述べる——スペクトラムアナライザとオシロスコープは全電力について一致しなければならない。
証明 ステップ1(特別な場合)。まずガウス関数 f ( x ) = e − π x 2 f(x) = e^{-\pi x^2} f ( x ) = e − π x 2 についてこの等式を確認する。指数部の平方完成による直接計算から f ^ ( ξ ) = e − π ξ 2 \hat f(\xi) = e^{-\pi \xi^2} f ^ ( ξ ) = e − π ξ 2 が得られ、 ∫ − ∞ ∞ ∣ f ( x ) ∣ 2 d x = ∫ − ∞ ∞ ∣ f ^ ( ξ ) ∣ 2 d ξ \int_{-\infty}^{\infty} |f(x)|^2\,dx = \int_{-\infty}^{\infty} |\hat f(\xi)|^2\,d\xi ∫ − ∞ ∞ ∣ f ( x ) ∣ 2 d x = ∫ − ∞ ∞ ∣ f ^ ( ξ ) ∣ 2 d ξ の両辺は同じガウス積分 ∫ e − 2 π x 2 d x \int e^{-2\pi x^2}\,dx ∫ e − 2 π x 2 d x に等しくなる。
ステップ2(線形性と畳み込み定理による拡張)。良い性質を持つ(シュワルツ)関数 f f f に対して、 f ~ ( x ) = f ( − x ) ‾ \tilde f(x) = \overline{f(-x)} f ~ ( x ) = f ( − x ) として g = f ∗ f ~ g = f * \tilde f g = f ∗ f ~ と置く。すると g ( 0 ) = ∫ ∣ f ( x ) ∣ 2 d x g(0) = \int |f(x)|^2\,dx g ( 0 ) = ∫ ∣ f ( x ) ∣ 2 d x であり、畳み込み定理 f ∗ g ^ = f ^ ⋅ g ^ \widehat{f * g} = \hat f \cdot \hat g f ∗ g = f ^ ⋅ g ^ とフーリエ反転公式により g ( 0 ) = ∫ g ^ ( ξ ) d ξ = ∫ ∣ f ^ ( ξ ) ∣ 2 d ξ g(0) = \int \hat g(\xi)\,d\xi = \int |\hat f(\xi)|^2\,d\xi g ( 0 ) = ∫ g ^ ( ξ ) d ξ = ∫ ∣ f ^ ( ξ ) ∣ 2 d ξ となる。これによりガウス関数を含むすべての滑らかで急減少なシュワルツ関数について等式が証明される。
ステップ3(稠密性の議論)。シュワルツ関数は L 2 ( R ) L^2(\mathbb{R}) L 2 ( R ) の中で稠密である: 任意の f ∈ L 2 f \in L^2 f ∈ L 2 は L 2 L^2 L 2 ノルムに関してシュワルツ関数の列 f n → f f_n \to f f n → f の極限として得られる。ステップ2によりフーリエ変換はこの稠密な部分空間上で等長写像であるから、それは ∫ − ∞ ∞ ∣ f ( x ) ∣ 2 d x = ∫ − ∞ ∞ ∣ f ^ ( ξ ) ∣ 2 d ξ \int_{-\infty}^{\infty} |f(x)|^2\,dx = \int_{-\infty}^{\infty} |\hat f(\xi)|^2\,d\xi ∫ − ∞ ∞ ∣ f ( x ) ∣ 2 d x = ∫ − ∞ ∞ ∣ f ^ ( ξ ) ∣ 2 d ξ を満たしたまま L 2 ( R ) L^2(\mathbb{R}) L 2 ( R ) 全体への有界作用素として一意に拡張される——この拡張こそが、定義積分 f ^ ( ξ ) = ∫ − ∞ ∞ f ( x ) e − 2 π i x ξ d x \hat f(\xi) = \int_{-\infty}^{\infty} f(x)\, e^{-2\pi i x \xi}\,dx f ^ ( ξ ) = ∫ − ∞ ∞ f ( x ) e − 2 π i x ξ d x が絶対収束しないことがある L 2 L^2 L 2 関数の『フーリエ変換』の意味である。
シュワルツ関数 f f f に対して、 ∑ n = − ∞ ∞ f ( n ) = ∑ k = − ∞ ∞ f ^ ( k ) \sum_{n=-\infty}^{\infty} f(n) = \sum_{k=-\infty}^{\infty} \hat f(k) ∑ n = − ∞ ∞ f ( n ) = ∑ k = − ∞ ∞ f ^ ( k ) が成り立つ:すべての整数にわたる関数の総和は、そのフーリエ変換のすべての整数にわたる総和に等しい。
なぜ正しいのか? これは信号を整数点で標本化することと、そのスペクトルを周期化することとの間の橋渡しである:左辺は標本 f ( n ) f(n) f ( n ) を足し合わせた結果であり、右辺はスペクトルから得られる結果である。この一つの等式が標本化定理、結晶学における格子和、そしてテータ関数やリーマンゼータ関数の関数等式の基礎となっている。
証明 ステップ1(周期化)。シュワルツ関数 f f f に対して F ( x ) = ∑ n = − ∞ ∞ f ( x + n ) F(x) = \sum_{n=-\infty}^{\infty} f(x+n) F ( x ) = ∑ n = − ∞ ∞ f ( x + n ) と定める。 f f f の急減少性によりこの和は絶対かつ一様に収束し、 F F F は周期1の滑らかな関数となるので、円周上で独自のフーリエ級数を持つ。
ステップ2(周期化されたもののフーリエ係数を計算する)。 F F F の第 k k k フーリエ係数は c k = ∫ 0 1 F ( x ) e − 2 π i k x d x = ∑ n ∫ 0 1 f ( x + n ) e − 2 π i k x d x c_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 c k = ∫ 0 1 F ( x ) e − 2 π ik x d x = ∑ n ∫ 0 1 f ( x + n ) e − 2 π ik x d x である。各項で y = x + n y = x+n y = x + n と置換し、 e − 2 π i k x e^{-2\pi i k x} e − 2 π ik x の n n n に関する周期性を使うと、断片は R \mathbb{R} R 全体にわたる一つの積分にまとまる: c k = ∫ − ∞ ∞ f ( y ) e − 2 π i k y d y = f ^ ( k ) c_k = \int_{-\infty}^{\infty} f(y) e^{-2\pi i k y}\,dy = \hat f(k) c k = ∫ − ∞ ∞ f ( y ) e − 2 π ik y d y = f ^ ( k ) 。
ステップ3( x=0 で評価する)。 F F F は滑らかなので、そのフーリエ級数は各点でそれ自身に収束し、特に x = 0 x=0 x = 0 では F ( 0 ) = ∑ k c k e 0 F(0) = \sum_k c_k e^{0} F ( 0 ) = ∑ k c k e 0 、すなわち ∑ n = − ∞ ∞ f ( n ) = ∑ k = − ∞ ∞ f ^ ( k ) \sum_{n=-\infty}^{\infty} f(n) = \sum_{k=-\infty}^{\infty} \hat f(k) ∑ n = − ∞ ∞ f ( n ) = ∑ k = − ∞ ∞ f ^ ( k ) が成り立つ。左辺は定義よりまさに F ( 0 ) = ∑ n f ( n ) F(0) = \sum_n f(n) F ( 0 ) = ∑ n f ( n ) であり、これで証明が完了する。
発展 ヒルベルト変換:最初の特異積分 ポアソン核やフェイエール核のような、単に有界で可積分な核に対する畳み込みはおとなしい。解析学における最も重要な作用素のいくつかは、それとは違い、原点に本物の特異性を持つ核——y = x y=x y = x の近くで可積分ではなく、その正の部分と負の部分の間の繊細な打ち消し合いによってのみ救われる核——に対して畳み込みを行う。その原型が実数直線上のヒルベルト変換 H f ( x ) = p.v. 1 π ∫ − ∞ ∞ f ( y ) x − y d y Hf(x) = \text{p.v.}\,\frac{1}{\pi}\int_{-\infty}^{\infty}\frac{f(y)}{x-y}\,dy H f ( x ) = p.v. π 1 ∫ − ∞ ∞ x − y f ( y ) d y であり、主値として理解される: y = x y=x y = x における特異点は対称に取り除かれてから極限が取られる。周波数側では、この特異な畳み込みは驚くほど単純な乗算になる: ( H f ^ ) ( ξ ) = − i s g n ( ξ ) f ^ ( ξ ) (\widehat{Hf})(\xi) = -i\,\mathrm{sgn}(\xi)\,\hat f(\xi) ( H f ) ( ξ ) = − i sgn ( ξ ) f ^ ( ξ ) ——ヒルベルト変換はあらゆる周波数に対する純粋な 90 ∘ 90^\circ 9 0 ∘ の位相回転であり、大きさを変えずに i i i の符号を反転させる。
H f ( x ) = p.v. 1 π ∫ − ∞ ∞ f ( y ) x − y d y Hf(x) = \text{p.v.}\,\frac{1}{\pi}\int_{-\infty}^{\infty}\frac{f(y)}{x-y}\,dy H f ( x ) = p.v. π 1 ∫ − ∞ ∞ x − y f ( y ) d y ( H f ^ ) ( ξ ) = − i s g n ( ξ ) f ^ ( ξ ) (\widehat{Hf})(\xi) = -i\,\mathrm{sgn}(\xi)\,\hat f(\xi) ( H f ) ( ξ ) = − i sgn ( ξ ) f ^ ( ξ ) これを見るもう一つの古典的な方法がある:実数直線上で定義された f f f に対して、 u = P y ∗ f u = P_y * f u = P y ∗ f をその上半平面への調和拡張(前節で扱ったポアソン核との畳み込み)とし、 v = Q y ∗ f v = Q_y * f v = Q y ∗ f を、共役ポアソン核 Q y Q_y Q y から作られるその調和共役として、F = u + i v F=u+iv F = u + i v が x + i y x+iy x + i y の正則関数となるようにする。 y → 0 + y\to 0^+ y → 0 + のとき、 v v v の境界値はまさに f f f のヒルベルト変換を復元する。これはこの理論全体の複素解析的な祖先である:境界上の特異積分作用素とは、実は一つ上の次元に住む正真正銘の正則関数の影にほかならない。
ヒルベルト変換 H H H は、すべての 1 < p < ∞ 1<p<\infty 1 < p < ∞ に対して L p ( R ) L^p(\mathbb{R}) L p ( R ) 上で有界である:あらゆる f f f について ∥ H f ∥ L p ≤ C p ∥ f ∥ L p \|Hf\|_{L^p} \le C_p\|f\|_{L^p} ∥ H f ∥ L p ≤ C p ∥ f ∥ L p を満たす定数 C p C_p C p が存在する。(マルセル・リース、1927年)
なぜ正しいのか? L 2 L^2 L 2 の場合は実質ただで得られる:乗数 − i s g n ( ξ ) -i\,\mathrm{sgn}(\xi) − i sgn ( ξ ) はほとんどいたるところで絶対値がちょうど 1 1 1 であるため、プランシュレルの定理により H H H は文字通り L 2 L^2 L 2 上の等長写像である——あらゆる周波数の 90 ∘ 90^\circ 9 0 ∘ 位相回転は全エネルギーを何も変えない。しかし有界性を他の p p p に拡張することは本質的に別の、はるかに難しい問題である:L 2 L^2 L 2 の外にはパーセバルの等式の類似物が存在しないため、代数的な近道は使えず、証明は H f Hf H f が大きくなる集合の大きさと幾何に関する実変数的な評価に頼らざるを得ない。
証明 ステップ1(L 2 L^2 L 2 、直接)。プランシュレルの定理と乗数の公式により、ξ ≠ 0 \xi\ne0 ξ = 0 で ∣ s g n ( ξ ) ∣ = 1 |\mathrm{sgn}(\xi)|=1 ∣ sgn ( ξ ) ∣ = 1 であることを用いて ∥ H f ∥ L 2 2 = ∫ ∣ H f ^ ( ξ ) ∣ 2 d ξ = ∫ ∣ s g n ( ξ ) ∣ 2 ∣ f ^ ( ξ ) ∣ 2 d ξ = ∥ f ∥ L 2 2 \|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 ∥ H f ∥ L 2 2 = ∫ ∣ H f ( ξ ) ∣ 2 d ξ = ∫ ∣ sgn ( ξ ) ∣ 2 ∣ f ^ ( ξ ) ∣ 2 d ξ = ∥ f ∥ L 2 2 が成り立つ(一点は積分に影響しない)。したがって H H H は L 2 ( R ) L^2(\mathbb{R}) L 2 ( R ) 上の等長写像であり、特に C 2 = 1 C_2=1 C 2 = 1 で有界である。
ステップ2(一般の 1 < p < ∞ 1<p<\infty 1 < p < ∞ 、後方参照)。指数の全範囲は、ステップ1からいかなる易しい議論によっても従わない。これは実変数コルデロン–ジグムント理論の正真正銘の定理である:ヒルベルト変換の核 1 / ( π x ) 1/(\pi x) 1/ ( π x ) は、コルデロン–ジグムント核の滑らかさと打ち消し条件をまさに満たしており、そのため次節で展開される一般的な特異積分の機構がそのまま適用できる。その機構は、ステップ1で証明した L 2 L^2 L 2 の評価から弱形 ( 1 , 1 ) (1,1) ( 1 , 1 ) の評価を生み出し、弱-( 1 , 1 ) (1,1) ( 1 , 1 ) 評価と L 2 L^2 L 2 評価の間のマルツィンキェヴィッチ補間(そして p > 2 p>2 p > 2 に対する双対性の議論)によって、すべての 1 < p < ∞ 1<p<\infty 1 < p < ∞ について証明が完成する——この議論の後半部分は、ここで繰り返す代わりに、すぐ下のコルデロン–ジグムントの節で完全に実行される。
例: 指示関数のヒルベルト変換
[ − 1 , 1 ] [-1,1] [ − 1 , 1 ] の指示関数 f = 1 [ − 1 , 1 ] f = \mathbf{1}_{[-1,1]} f = 1 [ − 1 , 1 ] を考える:有界でコンパクトな台を持つ、まったく取り立てて特徴のない関数である。H f Hf H f を求めよ。
解答 定義より H f ( x ) = 1 π p.v. ∫ − 1 1 d y x − y Hf(x) = \frac{1}{\pi}\,\text{p.v.}\int_{-1}^{1} \frac{dy}{x-y} H f ( x ) = π 1 p.v. ∫ − 1 1 x − y d y である。[ − 1 , 1 ] [-1,1] [ − 1 , 1 ] の外側では取り除くべき特異点はなく、1 / ( x − y ) 1/(x-y) 1/ ( x − y ) の y y y に関する原始関数は − log ∣ x − y ∣ -\log|x-y| − log ∣ x − y ∣ であるから、∫ − 1 1 d y x − y = [ − log ∣ x − y ∣ ] y = − 1 y = 1 = log ∣ x + 1 x − 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 x − y d y = [ − log ∣ x − y ∣ ] y = − 1 y = 1 = log x − 1 x + 1 となる(( − 1 , 1 ) (-1,1) ( − 1 , 1 ) の内部の x x x についても、y = x − y=x^- y = x − と y = x + y=x^+ y = x + における2つの発散が打ち消し合うため、同じ計算が正真正銘の主値として成り立つ)。これより H f ( x ) = 1 π log ∣ x + 1 x − 1 ∣ Hf(x) = \frac{1}{\pi}\log\left|\frac{x+1}{x-1}\right| H f ( x ) = π 1 log x − 1 x + 1 が得られる。
際立った特徴は x = ± 1 x=\pm1 x = ± 1 で起こることである:f f f 自身はそこで完全に有界である(単に 1 1 1 から 0 0 0 へ跳ぶだけ)にもかかわらず、H f Hf H f はまさにその2点で対数的に発散する。有界でコンパクトな台を持つ、教科書的に素性の良い関数が、H H H によって2つの正真正銘の特異点を持つ関数に変えられてしまう——これがあらゆる特異積分作用素に共通する局所的な特徴であり、まさにこれゆえに H H H を制御することは、素性の良い有界な核との畳み込みを制御するよりもはるかに繊細な評価を必要とするのである。
よくある誤り. H H H は L 1 ( R ) L^1(\mathbb{R}) L 1 ( R ) でも L ∞ ( R ) L^\infty(\mathbb{R}) L ∞ ( R ) でも有界ではない——M.リースの定理は両端点の手前で厳密に止まっており、これは証明上の隙間ではなく、本当の破綻である。上の例はすでにその仕組みを示している:f = 1 [ − 1 , 1 ] f=\mathbf{1}_{[-1,1]} f = 1 [ − 1 , 1 ] は L 1 ∩ L ∞ L^1\cap L^\infty L 1 ∩ L ∞ 関数として考えうる限り素性が良いにもかかわらず、H f Hf H f は x = ± 1 x=\pm1 x = ± 1 に除去不能な対数的特異点を持ち、それはまさに H f Hf H f をそれらの点の近くで L ∞ L^\infty L ∞ から締め出し、高次元における類似の現象が H f Hf H f の他の自然なノルムを無限大に押しやる原因となる、そのタイプの発散である。安易な工夫で易しい L 2 L^2 L 2 の評価を p = 1 p=1 p = 1 や p = ∞ p=\infty p = ∞ にまで拡張しようとしてもうまくいかない——代わりに本当に新しい実変数の機構(次に扱うコルデロン–ジグムント理論)が必要であり、その機構をもってしても端点 p = 1 p=1 p = 1 では、より弱い弱形 ( 1 , 1 ) (1,1) ( 1 , 1 ) の代替物しか回復できない。発展 ハーディ–リトルウッド極大関数とコルデロン–ジグムント分解 ヒルベルト変換だけでなく、一般の特異積分作用素についてM.リースの定理を完成させるには、解析学者たちはフーリエ変換について何も語らない、純粋に実変数的な道具を必要とした。ハーディとリトルウッドは1930年にそれを導入した:極大関数 M f ( x ) = sup r > 0 1 2 r ∫ ∣ y − x ∣ < r ∣ f ( y ) ∣ d y Mf(x) = \sup_{r>0} \frac{1}{2r}\int_{|y-x|<r} |f(y)|\,dy M f ( x ) = sup r > 0 2 r 1 ∫ ∣ y − x ∣ < r ∣ f ( y ) ∣ d y であり、これは x x x を中心とする任意の区間にわたる ∣ f ∣ |f| ∣ f ∣ のありうる最大の平均値である。これは、x x x の近くで f f f に対して行われうるあらゆる平均化操作を、一度にまとめて制御するための、無骨だが極めて有効な道具である。
M M M は弱形 ( 1 , 1 ) (1,1) ( 1 , 1 ) である:普遍定数 C C C (常に C = 5 C=5 C = 5 ととれる)に対して ∣ { x : M f ( x ) > λ } ∣ ≤ C λ ∥ f ∥ 1 |\{x : Mf(x) > \lambda\}| \le \frac{C}{\lambda}\|f\|_1 ∣ { x : M f ( x ) > λ } ∣ ≤ λ C ∥ f ∥ 1 が成り立つ。したがって補間により、M M M はすべての 1 < p ≤ ∞ 1<p\le\infty 1 < p ≤ ∞ に対して L p ( R ) L^p(\mathbb{R}) L p ( R ) 上でも有界である——しかし M M M 自身は L 1 ( R ) L^1(\mathbb{R}) L 1 ( R ) 上では決して有界にならない。
なぜ正しいのか? M f ( x ) Mf(x) M f ( x ) は、x x x の周りの区間上でとりうるあらゆる f f f の平均——あらゆる尺度での移動平均——を一挙に制御する。この一つの量を制御することが、まさにルベーグ微分定理(縮小していく x x x の周りの区間にわたる f f f の平均が、ほとんどすべての x x x で f ( x ) f(x) f ( x ) に収束する)を証明するために必要なことであることが判明する:M f Mf M f がほとんどいたるところ有限であるとわかれば、短くて易しい議論によって完全な微分の主張へと格上げできる。これは、易しい L 2 L^2 L 2 評価に対してプランシュレルの等式が果たした役割を、以下のコルデロン–ジグムント理論に対して果たす、実変数的な原動力である。
証明 ステップ1(良い球による被覆)。λ > 0 \lambda>0 λ > 0 を固定する。M f ( x ) > λ Mf(x)>\lambda M f ( x ) > λ を満たす各 x x x に対し、上限の定義より 1 ∣ B x ∣ ∫ B x ∣ f ∣ > λ \frac{1}{|B_x|}\int_{B_x}|f| > \lambda ∣ B x ∣ 1 ∫ B x ∣ f ∣ > λ を満たす x x x を中心とする区間 B x B_x B x が存在する。これらの球 { B x } \{B_x\} { B x } は集合 { M f > λ } \{Mf>\lambda\} { M f > λ } を被覆する。
ステップ2(ヴィタリの 5 r 5r 5 r 被覆補題)。半径が有界な球の任意の族から、可算な互いに素な部分族 { B i } \{B_i\} { B i } を常に取り出すことができ、5 5 5 倍に拡大した球 { 5 B i } \{5B_i\} { 5 B i } が元の族全体の和集合を依然として被覆するようにできる。これを { B x } \{B_x\} { B x } に適用して、{ M f > λ } ⊆ ⋃ i 5 B i \{Mf>\lambda\}\subseteq\bigcup_i 5B_i { M f > λ } ⊆ ⋃ i 5 B i を満たす互いに素な部分族 { B i } \{B_i\} { B i } を得る。
ステップ3(互いに素な部分を足し合わせる)。各 B i B_i B i は構成により λ ∣ B i ∣ < ∫ B i ∣ f ∣ \lambda|B_i| < \int_{B_i}|f| λ ∣ B i ∣ < ∫ B i ∣ f ∣ を満たす。B i B_i B i は互いに素であるから、i i i について足し合わせると λ ∑ i ∣ B i ∣ < ∑ i ∫ B i ∣ f ∣ ≤ ∫ R ∣ f ∣ = ∥ f ∥ 1 \lambda\sum_i|B_i| < \sum_i\int_{B_i}|f| \le \int_{\mathbb{R}}|f| = \|f\|_1 λ ∑ i ∣ B i ∣ < ∑ i ∫ B i ∣ f ∣ ≤ ∫ R ∣ f ∣ = ∥ f ∥ 1 となり、∑ i ∣ B i ∣ < ∥ f ∥ 1 / λ \sum_i|B_i| < \|f\|_1/\lambda ∑ i ∣ B i ∣ < ∥ f ∥ 1 / λ が得られる。最後に ∣ { M f > λ } ∣ ≤ ∑ i ∣ 5 B i ∣ = 5 ∑ i ∣ B i ∣ < 5 λ ∥ f ∥ 1 |\{Mf>\lambda\}| \le \sum_i|5B_i| = 5\sum_i|B_i| < \dfrac{5}{\lambda}\|f\|_1 ∣ { M f > λ } ∣ ≤ ∑ i ∣5 B i ∣ = 5 ∑ i ∣ B i ∣ < λ 5 ∥ f ∥ 1 となり、これはまさに ∣ { x : M f ( x ) > λ } ∣ ≤ 5 λ ∥ f ∥ 1 |\{x : Mf(x) > \lambda\}| \le \frac{5}{\lambda}\|f\|_1 ∣ { x : M f ( x ) > λ } ∣ ≤ λ 5 ∥ f ∥ 1 である。
定義: コルデロン–ジグムント分解
f ∈ L 1 ( R ) f\in L^1(\mathbb{R}) f ∈ L 1 ( R ) とある高さ α > 0 \alpha>0 α > 0 が与えられたとき、コルデロンとジグムント(1952年)は、直線を f f f と α \alpha α に適応した互いに素な『停止』二進区間に分割する方法を示した:一つの大きな二進区間から出発し、それを再帰的に二等分していき、二進区間 Q j Q_j Q j の平均 1 ∣ Q j ∣ ∫ Q j ∣ f ∣ \frac{1}{|Q_j|}\int_{Q_j}|f| ∣ Q j ∣ 1 ∫ Q j ∣ f ∣ が α \alpha α を超えたまさに最初の時点でそれを停止させ保持する(その親の平均はまだ ≤ α \le\alpha ≤ α だったので、二等分してもたかだか2倍にしかならない)。これにより、互いに素な停止区間の族 { Q j } \{Q_j\} { Q j } が生まれ、α < 1 ∣ Q j ∣ ∫ Q j ∣ f ∣ ≤ 2 α \alpha < \frac{1}{|Q_j|}\int_{Q_j}|f| \le 2\alpha α < ∣ Q j ∣ 1 ∫ Q j ∣ f ∣ ≤ 2 α を満たし、全長は ∑ ∣ Q j ∣ ≤ ∥ f ∥ 1 / α \sum|Q_j|\le \|f\|_1/\alpha ∑ ∣ Q j ∣ ≤ ∥ f ∥ 1 / α であり、⋃ j Q j \bigcup_j Q_j ⋃ j Q j の外側に残る部分ではほとんどいたるところ ∣ f ∣ ≤ α |f|\le\alpha ∣ f ∣ ≤ α となる。
[0,1] 上の固定されたでこぼこな |f|(x≈0.18 付近に細い尖った峰、x≈0.62 付近にほぼ特異的な盛り上がり、x≈0.86 付近に幅広い盛り上がりを持つ)を示し、この水準でのコルデロン–ジグムント停止区間が赤く塗られている:塗られた各区間が一つの停止立方体 Q_j であり、その上に引かれた太い赤線分がその平均値を示す——それは常に破線で示された α と 2α の基準線の間に収まる。スライダーを動かして α を変えると、停止区間が縮んだり、広がったり、分裂したり、併合したりする様子が見える。 この分解は f f f を『良い』部分と無限個の『悪い』断片に分ける、f = g + ∑ j b j f = g + \sum_j b_j f = g + ∑ j b j :良い部分 g g g は、すべての Q j Q_j Q j の外側では f f f に等しく、保持された各 Q j Q_j Q j の上では(だいたい α \alpha α 程度の大きさの)f f f の平均値に等しいので、いたるところ ∣ g ∣ ≤ 2 α |g|\le 2\alpha ∣ g ∣ ≤ 2 α であり、g g g は易しい L 2 L^2 L 2 理論で扱える。各悪い断片 b j b_j b j はそれ自身の立方体 Q j Q_j Q j 上に台を持ち、そこで平均ゼロである、∫ Q j b j d x = 0 \int_{Q_j} b_j\,dx = 0 ∫ Q j b j d x = 0 ——この打ち消し合いこそが、Q j Q_j Q j から十分離れれば滑らかな核による b j b_j b j の寄与が小さいままでいられる理由である。なぜなら核は Q j Q_j Q j の全体にわたってほぼ一定に見え、ほぼ一定の核は平均ゼロの関数を積分するとほとんど何も残らないからである。この良い–悪いの分割を、弱形 ( 1 , 1 ) (1,1) ( 1 , 1 ) のハーディ–リトルウッド極大評価と組み合わせ、マルツィンキェヴィッチ補間を通すことこそが、コルデロンとジグムントがヒルベルト変換だけでなく一般の特異積分作用素についてM.リースの定理を完成させ、L p L^p L p の有界性をすべての 1 < p < ∞ 1<p<\infty 1 < p < ∞ に拡張するために使った実変数的な機構であり、これによって上の定理へと話がつながって輪が閉じる。
大学 実世界での応用と具体例 フーリエ変換は『時間的な形』を『周波数的な内容』に変えるため、技術者や科学者が振動を分離・フィルタリング・計数する必要のあるあらゆる場面で標準的な道具となっている:オーディオイコライザー、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) f ( t ) = cos ( 2 π ⋅ 440 t ) + cos ( 2 π ⋅ 3000 t ) を記録した:440 Hz の音符(A4)と 3000 Hz の電子的なヒス雑音が混じっている。技術者は ∣ ξ ∣ > 1000 |\xi| > 1000 ∣ ξ ∣ > 1000 のときは常に f ^ ( ξ ) = 0 \hat f(\xi) = 0 f ^ ( ξ ) = 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) f ( t ) = cos ( 2 π ⋅ 440 t ) + cos ( 2 π ⋅ 3000 t ) のエネルギーは ξ = ± 440 \xi = \pm 440 ξ = ± 440 と ξ = ± 3000 \xi = \pm 3000 ξ = ± 3000 にのみ集中している( f ^ \hat f f ^ 中のスパイクとして)。
フィルタは ∣ ξ ∣ ≤ 1000 |\xi| \le 1000 ∣ ξ ∣ ≤ 1000 を満たすものをすべて保持し、それ以外を除去する。 440 ≤ 1000 440 \le 1000 440 ≤ 1000 だが 3000 > 1000 3000 > 1000 3000 > 1000 なので、 ± 3000 \pm 3000 ± 3000 のスパイクは除去され、 ± 440 \pm 440 ± 440 のスパイクはそのまま残る。
残ったものに反転公式を適用すると、生き残った周波数が正確に再構成される:出力は cos ( 2 π ⋅ 440 t ) \cos(2\pi \cdot 440\,t) cos ( 2 π ⋅ 440 t ) であり、ヒス雑音の消えたきれいな音符である。これはグラフィックイコライザーのローパスつまみが実際に行っていることそのものである。
例: 熱が広がる理由:熱方程式を解く
ある棒の初期温度分布が u ( x , 0 ) = f ( x ) u(x,0) = f(x) u ( x , 0 ) = f ( x ) であり、熱方程式 u t = u x x u_t = u_{xx} u t = u xx に従う。フーリエ変換( x x x について)を用いて、 t > 0 t>0 t > 0 における u ( x , t ) u(x,t) u ( x , t ) を求めよ。
解答 x x x についてフーリエ変換をとると、各空間微分 ∂ x \partial_x ∂ x は 2 π i ξ 2\pi i \xi 2 π i ξ による乗算になるので、 u x x u_{xx} u xx は − 4 π 2 ξ 2 u ^ -4\pi^2\xi^2 \hat u − 4 π 2 ξ 2 u ^ になる。方程式 u t = u x x u_t = u_{xx} u t = u xx は、固定した各 ξ \xi ξ に対する t t t についての常微分方程式 ∂ t u ^ ( ξ , t ) = − 4 π 2 ξ 2 u ^ ( ξ , t ) \partial_t \hat u(\xi,t) = -4\pi^2\xi^2\, \hat u(\xi,t) ∂ t u ^ ( ξ , t ) = − 4 π 2 ξ 2 u ^ ( ξ , t ) になる——難しい偏微分方程式が、互いに独立な常微分方程式の族という扱いやすいものになる大きな簡略化である。
この常微分方程式の解は u ^ ( ξ , t ) = f ^ ( ξ ) e − 4 π 2 ξ 2 t \hat u(\xi,t) = \hat f(\xi)\, e^{-4\pi^2 \xi^2 t} u ^ ( ξ , t ) = f ^ ( ξ ) e − 4 π 2 ξ 2 t であり、高周波数 ξ \xi ξ を急速に消し去る指数減衰因子である:細かい空間的な詳細は急速に滑らかになり、これは鋭い温度スパイクが最初に平らになるという日常的事実と一致する。
変換を逆にすると(畳み込み定理を逆向きに使えば、周波数領域での積は空間領域での畳み込みになる) u ( x , t ) = ( f ∗ K t ) ( x ) , K t ( x ) = 1 4 π t e − x 2 / ( 4 t ) u(x,t) = (f * K_t)(x),\ \ K_t(x) = \dfrac{1}{\sqrt{4\pi t}}\, e^{-x^2/(4t)} u ( x , t ) = ( f ∗ K t ) ( x ) , K t ( x ) = 4 π t 1 e − x 2 / ( 4 t ) が得られる:時刻 t t t における温度は、初期分布が広がるガウス型の山 K t K_t K t によってぼかされた(畳み込まれた)ものである——これはまさに熱が拡散する様子の数学的な描像である。
定義: 不確定性原理
関数は時間と周波数の両方で同時に鋭く局在することはできない。 ∣ f ∣ 2 |f|^2 ∣ f ∣ 2 と ∣ f ^ ∣ 2 |\hat f|^2 ∣ f ^ ∣ 2 の広がりをそれぞれ Δ x \Delta x Δ x と Δ ξ \Delta \xi Δ ξ で測ると、 Δ x ⋅ Δ ξ ≥ 1 4 π \Delta x \cdot \Delta \xi \ge \dfrac{1}{4\pi} Δ x ⋅ Δ ξ ≥ 4 π 1 が成り立つ。パルスを時間的に押しつぶす(小さい Δ x \Delta x Δ x )と、そのスペクトルは広がらざるを得ず(大きい Δ ξ \Delta \xi Δ ξ )、逆もまた然りである——これは量子力学におけるハイゼンベルクの不確定性原理を導くのと同じトレードオフである。
Δ x ⋅ Δ ξ ≥ 1 4 π \Delta x \cdot \Delta \xi \ge \dfrac{1}{4\pi} Δ x ⋅ Δ ξ ≥ 4 π 1 よくある誤り. 『コンパクトな台を持つ』信号(有限の時間窓でのみ非ゼロ)は、コンパクトな台を持つスペクトル(有限の帯域でのみ非ゼロ)も持てると考えたくなるかもしれない。しかし不確定性原理は、ゼロでない任意の関数についてこれを排除する: f f f と f ^ \hat f f ^ が同時にコンパクトな台を持つことは決してない。これが、(上の例のローパスフィルタのような)理想的な『壁』フィルタが数学的な理想化にすぎない理由である——本当に時間的に制限された信号は、遠く離れたところでは振幅がごく小さいとはいえ、必ずすべての周波数に漏れ出す。 歴史的ノート
1854年の三角級数に関する教授資格論文の中で、ベルンハルト・リーマンはどのような関数がそのような級数で表現できるかを厳密に研究し、(厳密な積分の概念、および可積分な f f f について ∣ ξ ∣ → ∞ |\xi| \to \infty ∣ ξ ∣ → ∞ のとき f ^ ( ξ ) → 0 \hat f(\xi) \to 0 f ^ ( ξ ) → 0 となるというリーマン–ルベーグの補題という)基礎を築いた。それは後にフーリエ解析および実数直線上の調和解析の厳密な理論に不可欠なものとなった。
ベルンハルト・リーマン ジョゼフ・フーリエ
研究の最前線 2026年時点
調和解析には依然として未解決の鋭い問題が多く残っており、進展は突然訪れることがある。制限予想(restriction conjecture)—— f ^ \hat f f ^ の情報を球面のような曲面へ制限することでどれだけ復元できるか——は次元2でのみ完全に解決されており(Fefferman、Carleson–Sjölin、Zygmund、1970年代)、デカップリング法(Bourgain–Demeter、2015年)による深い部分的結果があるにもかかわらず n ≥ 3 n \ge 3 n ≥ 3 では未解決のままである。これと密接に関係する幾何学的な問題であるカケヤ予想(あらゆる方向に単位長さの線分を含む集合は小さくなり得るか)は、何十年もの挑戦に耐えた末、2025年にホン・ワンとジョシュア・ザールによって三次元の場合にのみ解決された。 n ≥ 4 n \ge 4 n ≥ 4 の場合は両予想ともまだ未解決である。これらは孤立したパズルではない:カケヤ型集合に関する評価は制限予想やボホナー–リース評価に直接影響し、それらの評価がフーリエ級数の収束や分散型偏微分方程式の適切性を制御している。 これらの結果は、Bochner、Riesz、Steinらによって発見された一連の含意の連鎖を成している:カケヤ予想における進展はフーリエ制限予想に対する評価に直接寄与し、それがひいては一次元より高い次元でフーリエ積分を足し合わせるボホナー–リース予想を制御する。
規約 f ^ ( ξ ) = ∫ f ( x ) e − 2 π i x ξ d x \hat f(\xi) = \int f(x) e^{-2\pi i x\xi}\,dx f ^ ( ξ ) = ∫ f ( x ) e − 2 π i x ξ d x のもとで、ガウス関数 f ( x ) = e − π x 2 f(x) = e^{-\pi x^2} f ( x ) = e − π x 2 のフーリエ変換は:
f ^ ( ξ ) = e − π ξ 2 \hat f(\xi) = e^{-\pi \xi^2} f ^ ( ξ ) = e − π ξ 2 f ^ ( ξ ) = 0 \hat f(\xi) = 0 f ^ ( ξ ) = 0 f ^ ( ξ ) = e ξ \hat f(\xi) = e^{\xi} f ^ ( ξ ) = e ξ f ^ ( ξ ) = π e − π ξ 2 \hat f(\xi) = \pi e^{-\pi \xi^2} f ^ ( ξ ) = π e − π ξ 2 プランシュレルの定理 ∫ − ∞ ∞ ∣ f ( x ) ∣ 2 d x = ∫ − ∞ ∞ ∣ f ^ ( ξ ) ∣ 2 d ξ \int_{-\infty}^{\infty} |f(x)|^2\,dx = \int_{-\infty}^{\infty} |\hat f(\xi)|^2\,d\xi ∫ − ∞ ∞ ∣ f ( x ) ∣ 2 d x = ∫ − ∞ ∞ ∣ f ^ ( ξ ) ∣ 2 d ξ を最もよく表す記述は:
時間領域と周波数領域の間での全エネルギーの保存 フーリエ変換の周期性 跳躍不連続点におけるフーリエ級数の収束 原始関数の一意性
MRI装置は空間信号のフーリエ変換の標本(『kスペース』)を測定する。低周波数の標本( ∣ ξ ∣ |\xi| ∣ ξ ∣ が小さいもの)のみを収集した場合、再構成される画像は:
ぼやけて細部を欠く(輪郭は高周波数に存在する) すべての周波数を使う場合より鮮明になる 低周波数がすべての情報を含むため同一になる ノイズのみを含む
不確定性原理 Δ x ⋅ Δ ξ ≥ 1 4 π \Delta x \cdot \Delta \xi \ge \dfrac{1}{4\pi} Δ x ⋅ Δ ξ ≥ 4 π 1 が意味するのは:
ゼロでない関数は、fとその変換の両方がコンパクトな台を持つことはできない すべての関数は周期的なフーリエ変換を持つ フーリエ変換は線形ではない 畳み込みと乗算は同じ演算である
フェイエールの定理は、チェザロ平均 σ N f \sigma_N f σ N f が任意の連続関数 f f f に一様収束することを述べている。カールソンの定理はこれとは別の、はるかに深い、次のことについての主張である:
生の部分和 S N f S_N f S N f 自体がほとんどいたるところ収束するが、それは f f f が L 2 L^2 L 2 かそれ以上であると仮定した場合に限られる チェザロ平均 σ N f \sigma_N f σ N f も、任意の可積分な f f f に対して一様収束する ディリクレ核 D N D_N D N はフェイエール核と異なり非負である すべての連続関数はそのフーリエ級数が一様収束する M.リースの定理により、ヒルベルト変換 H H H が L p ( R ) L^p(\mathbb{R}) L p ( R ) 上で有界となるのはちょうど次の場合である:
例外なく、すべての 1 ≤ p ≤ ∞ 1\le p\le\infty 1 ≤ p ≤ ∞ p = 2 p=2 p = 2 の場合のみすべての 1 < p < ∞ 1<p<\infty 1 < p < ∞ だが、p = 1 p=1 p = 1 や p = ∞ p=\infty p = ∞ ではない H H H が特異積分であるため、p = 1 p=1 p = 1 の場合のみ