フーリエ変換

数学においてフーリエ変換(フーリエへんかん、英: Fourier transform、FT)は、実変数の複素または実数値関数を、その周波数成分を表す別の関数に写す変換である。
工学においては、変換後の関数はもとの関数に含まれる周波数を記述していると考え、しばしばもとの関数の周波数領域表現 (frequency domain representation) と呼ばれる。言い換えれば、フーリエ変換は関数を(離散的な正弦波・余弦波の和であるフーリエ級数とは対照的に)連続的な周波数成分の集まりとして分解する操作である。
フーリエ変換 (FT) は他の多くの数学的な演算と同様にフーリエ解析の主題を成す。「フーリエ変換」という言葉は関数の周波数領域表現のことを指すこともあるし、関数を周波数領域表現へ写す変換の過程・公式を言うこともある。なおこの呼称は、19世紀フランスの数学者・物理学者で次元解析の創始者とされるジョゼフ・フーリエに由来する。
概要
[編集]フーリエ変換は、関数をその周波数成分に分解し、どの周波数がどれだけ含まれているかを示す「周波数スペクトル」に変換する数学的な操作である。
時間や空間に関する関数(波形など)を、その関数を構成する正弦波(または複素指数関数)の集まりとして表現し直す操作と見なすことができる。元の関数が時間領域(または空間領域)の関数 (t) であれば、フーリエ変換後の関数 (ξ) は周波数領域(または波数領域)の関数となり、元の関数 f(t) に周波数 ξ の成分がどれだけ含まれているかを表す。
この変換は可逆であり、フーリエ逆変換と呼ばれる操作によって、周波数スペクトル (ξ) から元の関数 f(t) を復元することができる。
フーリエ変換は、数学におけるフーリエ解析の中心的な道具であるだけでなく、物理学、工学(特に信号処理、画像処理、通信工学)、量子力学など、現代科学の非常に広範な分野で不可欠なツールとなっている。
なお、「時間」と「周波数」(あるいは「空間座標」と「空間周波数」)という呼称がしばしば用いられるが、これらは物理的な応用に由来する慣用的表現であり、フーリエ変換が時間(あるいは空間座標)を変数としてのみ行われるということではない。物理量を含む関数にフーリエ変換を適用する場合、積分核の指数関数に含まれる変数の積(xξ など)が無次元量となる必要があるため、変換前の変数が特定の物理的次元を持つ場合、変換後の変数の次元はその逆数となる。例えば、元の関数が質量や加速度の関数であれば、変換後の関数はそれぞれ質量の逆数や加速度の逆数の関数となる。確率論において確率密度関数から特性関数を導出する際にも、このような時間・周波数以外の変数の組み合わせが用いられる。また、偏微分方程式の解法への応用でも、「両辺をフーリエ変換する」ことで解が導き出せる場合があるが、ここでも時間(あるいは空間)を変数とするフーリエ変換に限定されるものではなく、積分可能であればどの変数に関するフーリエ変換を行っても良い。数学では、物理量に依存しない形式でフーリエ変換は記述される。
定義
[編集]フーリエ変換の定義には、用いる周波数の変数(周波数 ξ [ヘルツ] か、角周波数 ω = 2πξ [ラジアン/秒] か)や、変換・逆変換の係数(正規化係数)の選び方によって、いくつかの流儀が存在する。分野によって慣習的に用いられる定義が異なるため注意が必要である。
主な定義として以下の3つが挙げられる。
| 周波数 ξ(ヘルツ) | ユニタリ | |
|---|---|---|
| 角周波数 ω(ラジアン毎秒) | 非ユニタリ | |
| ユニタリ |
本項では、特に断りのない限り、数学や一部の信号処理の文献で用いられることのある、周波数 (ヘルツ)を用い、変換・逆変換の係数が である定義1(ユニタリ・周波数)を採用する。この定義は、フーリエ変換と逆変換が対称形( の肩の符号が異なるだけ)である点に特徴がある。 また、この定義(および定義3)が「ユニタリ(英: unitary)」と呼ばれるのは、プランシュレルの定理で示されるように、関数のノルム(信号のエネルギーに相当)が変換前後で保存される、すなわち が成り立つためであり、フーリエ変換が(空間における)ユニタリ作用素であることを意味する。定義2(非ユニタリ)では、このような単純なエネルギー保存則は成立しない(係数 が必要となる)。
絶対可積分関数に対する定義
[編集]ここでは、最も基本的な 1次元の場合について定義を述べる。
絶対可積分関数 f: R → C のフーリエ変換は、
で定義される。ここでギリシャ文字小文字の ξ は任意の実数(周波数)である。
対象の関数における独立変数が物理量の場合、フーリエ変換は独立変数の次元をもとの逆数に移す。例えば、変換前の関数における独立変数 x が時間の次元をもつとき、変換後の独立変数 ξ は周波数の次元を持つ。あるいは、変換前の独立変数 x が長さの次元をもつとき、変換後の独立変数 ξ は波数の次元を持つ。この性質は定義より x ξ が無次元量であることから従う。
適当な条件のもと、f はその変換 からフーリエ逆変換 (inverse transform)
によって復元することができる(x は任意の実数)。
緩増加超関数への拡張
[編集]絶対可積分でない関数、例えば f(x)=c (定数関数)や f(x)=cos(x)(周期関数)、あるいはディラックのデルタ関数 δ(x) のような「関数」ではない対象についてもフーリエ変換を定義したいという強い要請が物理学や工学には存在する。
このため、フーリエ変換の定義を緩増加超関数(またはシュワルツ超関数)の空間 に拡張することが行われる。
緩増加超関数 は、シュワルツ空間 (急減少関数、すなわち無限回微分可能で、自身およびその全ての導関数が で多項式の逆数よりも速く 0 に近づく滑らかな関数 の集まり)上の連続線型汎関数として定義される。つまり、 は各 に対して複素数 を対応させる写像である。
(絶対可積分関数 も、 により、 という超関数とみなすことができる。)
この拡張により、定数関数や多項式など、無限遠で減衰しない関数のフーリエ変換が可能となる。しかし、 のような指数関数的に増大する関数は、緩増加超関数の枠組みでも扱うことができない。
シュワルツ空間 はフーリエ変換 の下で閉じている(急減少関数のフーリエ変換は再び急減少関数となる)という重要な性質を持つ。
この性質を利用し、緩増加超関数 T のフーリエ変換 は、任意の急減少関数 に対して、以下の関係(プランシュレルの定理の一般化)を満たす超関数 として定義される。
ここで は急減少関数 の(絶対可積分関数としての)フーリエ変換である。この定義により、超関数のフーリエ変換は再び緩増加超関数となる。
例:定数関数 を超関数 ()とみなす。そのフーリエ変換 は、定義より となる。 ここで、 のフーリエ逆変換の定義 において とおくと、 という関係が得られる(この関係は、急減少関数 に対して成立する)。 したがって、 である。 一方、 という超関数 は、定義により である。 したがって、、すなわち定数関数 のフーリエ変換は となる。
例:デルタ関数 (すなわち )のフーリエ変換 を計算する。定義より、任意の急減少関数 に対して以下が成り立つ。 ここで、 は のフーリエ変換の定義式において としたものに等しいため、 である。この右辺は、定数関数 が定める超関数 の作用 に他ならない。したがって、超関数の意味において、デルタ関数 のフーリエ変換は定数関数 となる。
フーリエ変換の性質
[編集]実数直線上で定義される関数 f が絶対可積分であるとは、
を満たすルベーグ可測関数であることをいう。
なお、ルベーグ積分においては、可測関数が絶対可積分であることの必要十分条件は、可積分であること(積分が有限な値を持つこと)であるため、単に可積分と言っても良い。
基本性質
[編集]絶対可積分関数 f(x), g(x), h(x) が与えられたとき、これらのフーリエ変換をそれぞれ (ξ), (ξ), (ξ)で表す。フーリエ変換は以下の基本性質を満たす[1]。
- 線型性
- 任意の複素数 a, b について h(x) = aƒ(x) + bg(x) であるならば
- が成り立つ。
- 平行移動
- 任意の実数 x0 に対して h(x) = ƒ(x − x0) であるならば
- が成り立つ。
- 変調
- 任意の実数 ξ0 に対して h(x) = e2πixξ0ƒ(x) ならば
- が成り立つ。
- 定数倍
- 非零実数 a に対し、h(x) = ƒ(ax) ならば
- が成り立つ。a = −1 つまり h(x) = ƒ(−x) の場合には、時間反転性 (time-reversal property)
- が導かれる。
- 複素共役
- f(x) の複素共役 f(x) について
- が成り立つ。
- 微分
- 関数 が( で 0 に収束するなどの条件を満たし) 回微分可能で、その導関数も絶対可積分であるならば、 が成り立つ。逆に、関数に を掛けたもののフーリエ変換は周波数領域での微分となり、(または )が成り立つ。
- 畳み込み
- h(x) = (f ∗ g)(x) ならば
- が成り立つ。
一様連続性とリーマン・ルベーグの補題
[編集]絶対可積分関数のフーリエ変換は、常に成り立つというわけではない性質も持っている。絶対可積分関数 f のフーリエ変換は一様連続で
を満たす[2]。絶対可積分関数のフーリエ変換は
であることを述べたリーマン・ルベーグの補題をも満足する[3]。絶対可積分函数 f のフーリエ変換 は有界連続だが絶対可積分であるとは限らず、その逆変換をルベーグ積分として書くことは一般にはできない。しかしながら、f および がともに絶対可積分ならば、反転公式
が殆ど全ての x において成り立つ。つまり、f は右辺で定義される連続関数と殆ど至る所等しい。特に f が実数直線上の連続関数として与えられたならば全ての x において等式が成り立つ。
前述の結果としてわかることは、フーリエ変換が L1(R) 上単射であることである。
プランシュレルの定理とパーセバルの定理
[編集]f(x) および g(x) は絶対可積分であるとし、そのフーリエ変換をそれぞれ (ξ) および (ξ) と表す。f(x) および g(x) がともに自乗絶対可積分であるならばパーセバルの定理
が成立する[4]。ここで上付きバーは複素共役を表す。
パーセバルの定理と同値なプランシュレルの定理によれば
が成立する[5]。プランシュレルの定理により、L2(R) に属する関数の後述する意味でのフーリエ変換を定義することが可能になる。プランシュレルの定理は、フーリエ変換はもとの量のエネルギーを保存するという自然科学における解釈を持つ。著者によってはこれらの定理のどちらともをプランシュレルの定理あるいはパーセバルの定理と呼んでいる場合があるので注意を要する。
局所コンパクトアーベル群に関する文脈におけるフーリエ変換の概念の一般の定式化についてはポントリャーギン双対の項を参照されたい。
不確定性関係
[編集]一般的に言って、f(x) が凝縮されればされるほどそのフーリエ変換 はより拡散される。特に、フーリエ変換のスケール性からわかることとして、関数を x において「圧搾」するならば、そのフーリエ変換は ξ において「伸展」される。したがって、関数とそのフーリエ変換の両方ともを勝手に凝縮させることはできない。
関数とそのフーリエ変換のコンパクト化のあいだの得失評価は不確定性関係の形で定式化することができる。f(x) は絶対可積分かつ自乗絶対可積分であると仮定する。一般性を失うことなく関数 f(x) は
に正規化されているものと仮定してよい。このとき、プランシュレルの定理により (ξ) も同様に正規化される。
x = 0 の周りでの拡散を
で定義される「0 の周りでの分散」(dispersion about zero) によって測ることにする[1]。確率の言葉で言えば、これは f(x)|2 の 0 の周りでの二次のモーメントである。
このとき不確定性原理は、関数 f(x) が絶対連続で、関数 x·f(x) および f′(x) が自乗絶対可積分であるならば
が成り立つことを述べる[1]。等式が成立するのは
したがって、
である場合に限る。ただし、定数 σ > 0 は任意であり、係数 C1 は f を L2-正規化する定数である[1]。言い換えれば、 f は 0 を中心に持つ(正規化)ガウス関数のとき等号が成り立つ。
事実として、この不等式は任意の x0, ξ0 ∈ R について
が成立することをも含む[6]。
量子力学において、運動量と位置の波動関数は(プランク定数を因子に持つ)フーリエ変換対である。プランク定数でスケールしなおせば、上述の不等式はロバートソンの不確定性関係を記述する。これは、ハイゼンベルグが構想した不確定性原理そのものではないが、深い関係がある。
ポアソン和公式
[編集]ポアソン和公式はフーリエ変換とフーリエ級数の間の関連性を提供する。絶対可積分関数 f ∈ L1(Rn) が与えられたとき、f の周期化が
によって与えられる。このとき、ポアソン和公式は f のフーリエ級数を f のフーリエ変換に結びつけるもので、特に f のフーリエ級数は
で与えられることを述べるものである。ポアソン和公式を用いて、大きな次元のユークリッド球面における格子点の数に対するランダウの漸近公式を導出することができる。また、絶対可積分函数 f と がともにコンパクト台を持つならば f = 0 を示すこともできる[1]。
畳み込み定理
[編集]フーリエ変換は、関数の畳み込みと関数の(点毎の)積とを相互に変換する。ƒ(x) および g(x) が絶対可積分関数であるとし、そのフーリエ変換をそれぞれ (ξ) および (ξ) で表す。さらに f と g との畳み込みが存在して絶対可積分であるならば、この畳み込みのフーリエ変換はフーリエ変換 (ξ) と (ξ) との積で与えられる(ただし、フーリエ変換の定義の仕方によっては定数因子が現れる場合もある)。
これを式で表せば、∗ を畳み込みとして
と表されるとき、
が成立することを意味する。線型時不変 (LTI) 系理論において、f(x) を単位インパルスで置き換えたものが h(x) = g(x) を与えることから、通例 g(x) は、入力 ƒ(x) と出力 h(x) に関する LTI 系のインパルス応答として解釈される。この場合、(ξ) はこの系の周波数応答を表す。
逆に、f(x) がふたつの自乗絶対可積分函数 p(x) および q(x) の積に分解されるならば、f(x) のフーリエ変換は、各因子のフーリエ変換 (ξ) および (ξ) の畳み込みで与えられる。
相互相関定理
[編集]同様の方法で、h(x) が f(x) と g(x) との相互相関
であるならば h(x) のフーリエ変換が
で与えられることが示される。
固有関数
[編集]L2(R) の正規直交基底の重要な一つはエルミート函数系
で与えられる。ここで Hn(x) は「確率論者の」エルミート多項式と呼ばれる、 で定義される関数である。この規約の下、フーリエ変換は
で与えられる。言い換えれば、エルミート関数系は L2(R) 上のフーリエ変換の固有関数からなる完全正規直交系を成す[1]。しかしながら、この固有関数系の選び方は一意ではなく、フーリエ変換の相異なる固有値は {±1, ±i} の 4 つしかなく、同じ固有値に属する固有関数の任意の線型結合はふたたび固有関数になる。この結果として L2(R) を 4 つの空間 H0, H1, H2, H3 で、フーリエ変換が Hk 上で単に ik-倍として作用するものの直和に分解することができる。この方法によるフーリエ変換の定義はウィーナーによる[7]。エルミート関数を選ぶのが便利なのは、それらが周波数域と時間域の両方で指数関数的に局在することと、それゆえに時間周波数解析において用いられる非整数次フーリエ変換が得られることにある [要出典]。
球面調和関数
[編集]で次数 k の斉次調和多項式全体の成す集合を表す。集合 は体球面調和関数系として知られる。高次元において体球面調和関数系はエルミート多項式と同様の役割を演じる。具体的には、 の適当な P(x) に対し、f(x) = e−π|x|2P(x) のフーリエ変換は
で与えられる。集合 を f(|x|)P(x) (P(x) ∈ ) の形の関数から作られる線型結合全体の成す集合の L2(Rn) における閉包とする。このとき、空間 L2(Rn) は空間 の直和に分解され、フーリエ変換は各空間 をそれ自身に移す。また、各空間 へのフーリエ変換の作用を特徴付けることができる[3]。f(x) = f0(|x|)P(x) (P(x) ∈ ) と表される関数のフーリエ変換は
となる。ただし、
であり、J(n + 2k − 2)/2 は次数 (n + 2k − 2)/2 の第一種ベッセル関数である。k = 0 のとき、これは動径関数のフーリエ変換に対する有用な公式を与える[8]。
フーリエ級数からの類推による導入
[編集]この節では、フーリエ変換の「動機」の一つとして、フーリエ級数からの類推的な導入を解説する。この方法は、フーリエ変換がどのようにして関数を周波数成分の「連続的な集まり」として表現し直すかという直観的なイメージを掴む上で、歴史的にも教育的にも重要である。
フーリエ変換を考える動機はフーリエ級数の研究に始まる。フーリエ級数の研究において、複雑な周期関数は単純な波動の数学的な表現である正弦関数や余弦関数の和として表される。正弦や余弦の性質のおかげで、この和に現れる各波の量、フーリエ係数を積分によって計算することができる。
多くの場合に、(オイラーの公式)を用いて、正弦関数および余弦関数の代りに基本波動 を用いた方が便利である。この場合には多くの公式が簡単化され、本項で後述するフーリエ変換のほかの類似の定式化をあたえるという点に優位性がある。この正弦・余弦から複素指数関数への移行にはフーリエ係数が複素数値であることを要する。この複素数は、関数に含まれる波動の振幅(あるいは大きさ)と、位相(あるいは初期角)の両方を与えているものと通常は解釈される。また、この移行に際して「負の周波数」も導入される。例えば、波動 および はともに周期1を持つが、複素フーリエ級数においては別々の成分として取り扱われる。したがって、周波数を単純に周期の逆数と考えることはできなくなる。
フーリエ級数を以下のようにしてフーリエ変換の動機付けに用いることができる。関数 ƒ をある区間 [−L/2, L/2] の外側で 0 となるようなものとすると、任意の T ≥ L に対して ƒ を区間 [−T /2, T /2] 上のフーリエ級数に拡張できる。ここで f のフーリエ級数に現れる波動 の係数となる で表される「量」は
で与えられ、ƒ は公式
で与えられなければならない。ξn = n/T とおき、Δξ = (n + 1)/T − n/T = 1/T とおくと、最後の和をリーマン和
として考えることができる。T → ∞ とすることにより、このリーマン和は定義節で与えられるフーリエ逆変換(積分)に収束すると期待される。このように、フーリエ変換は、関数に含まれる個々の周波数がどの程度あるかを(スペクトル密度 として)測るものと考えることができ、それらの波動を積分(「連続和」)によって再結合して元の関数を復元する操作がフーリエ逆変換に対応する。
以下の画像はフーリエ変換が特定の関数に含まれる周波数を測る方法を視覚的に現したものである。関数として、(t が秒で測られる場合には)3 ヘルツで振動し、急速に 0 になる
を描く。この関数は特に描画しやすい実フーリエ変換をもつものとして選ばれたものであり、最初の画像はそのグラフである。(3) を計算するために、e−2πi(3t)ƒ(t) を積分する。二枚目の画像はこの被積分関数の実部および虚部である。被積分関数の実部は殆ど常に正となる。これは ƒ(t) が負であるときには e−2πi(3t) の実部が同様に負となることによる。それらは同じ比率で振動するから、ƒ(t) が正であるときも同様に e−2πi(3t) の実部も正になる。
この結果、被積分関数の実部のを積分すれば、比較的大きな数値(ここでの場合 0.5)を得ることになる。
一方、((5) を見る場合のように)含まれない周波数を測れば、被積分関数は十分に振動し、それゆえにその積分はとても小さい値となる。一般の設定ではこれよりは少し複雑になるが、それでもフーリエ変換は関数 ƒ(t) に含まれる個々の周波数がどれくらいあるかを測るものという考え方に変わりはない。
この例では、ではなく、と、不自然ながかかっているが、は、(−∞, ∞)で絶対可積分関数としてのフーリエ変換(積分)が定義できないためである。(前述の通り、超関数としては定義可能である。)一方、は、のため(−∞, ∞)で絶対可積分であり、 から が成り立つ。よって、 は全体として絶対可積分関数となり、フーリエ変換(積分)が可能となっているものである。
- 3ヘルツの振動を示すもとの関数
- 3ヘルツにおけるフーリエ変換の被積分関数の実部および虚部
- 5ヘルツにおけるフーリエ変換の被積分関数の実部および虚部
- 3ヘルツおよび5ヘルツでラベル付けされたフーリエ変換
歴史
[編集]フーリエ変換の起源は、19世紀初頭のフランスの数学者・物理学者ジョゼフ・フーリエによる熱伝導の研究に遡る。
フーリエの研究
[編集]1807年、フーリエはフランス科学アカデミーに提出した論文において、熱伝導方程式の解を三角関数の無限級数(後にフーリエ級数と呼ばれる)として表現する方法を提示した[9]。この研究は、任意の周期関数が正弦波と余弦波の和として表現できるという革新的な主張を含んでいた。
当初の「不連続な関数までもが三角級数で表現できる」というフーリエの主張は、厳密な数学的証明を欠いていたため、論争を引き起こした。フーリエの論文は直ちには出版されず、1822年になってようやく『熱の解析的理論』(Théorie analytique de la chaleur)として刊行された[10]。
理論の厳密化と発展
[編集]フーリエの研究に触発され、19世紀の数学者たちはフーリエ級数の収束条件や表現可能な関数のクラスを明らかにする研究を進めた。ペーター・グスタフ・ディリクレは1829年に、区分的に単調で不連続点が有限個しかない関数についてフーリエ級数が収束することを証明し、フーリエ級数論に厳密な基礎を与えた[11]。
19世紀後半には、ベルンハルト・リーマンがリーマン積分の理論を発展させ、より広い関数のクラスに対してフーリエ級数を定義できるようになった[12]。
非周期関数への拡張
[編集]周期関数を対象とするフーリエ級数から、非周期関数を扱うフーリエ変換(フーリエ積分)への移行は、19世紀後半から20世紀初頭にかけて進展した。周期を無限大に近づける極限操作により、離散的な周波数成分の和(級数)が連続的な周波数成分の積分(フーリエ変換)に移行するという考え方が確立された。
20世紀に入ると、アンリ・ルベーグによる測度論と積分論の革新により、フーリエ変換の理論は大きく発展した[13]。1910年、ミシェル・プランシュレルは、自乗可積分関数(関数)に対してフーリエ変換がノルムを保存するユニタリ作用素であることを示す定理(プランシュレルの定理)を証明した[14]。この定理は、フーリエ変換を空間全体に拡張する数学的基礎を提供し、後の量子力学における応用にも不可欠なものとなった。
1920年代から1930年代にかけて、ノーバート・ウィーナーは一般化調和解析の理論を構築し、フーリエ変換を非周期的な信号の解析に適用する枠組みを整備した[15]。また、ローラン・シュワルツは1940年代後半に超関数の理論を確立し、デルタ関数や定数関数のような通常の関数ではない対象にもフーリエ変換を適用できるようにした[16]。
計算技術の発展
[編集]理論の発展と並行して、フーリエ変換の数値計算法も進歩した。1965年、ジェームズ・クーリーとジョン・テューキーは高速フーリエ変換(FFT)アルゴリズムを発表した[17]。FFTは、離散フーリエ変換の計算量を劇的に削減する画期的な手法であり、点のデータに対して従来のからへと計算量を減少させた。
一般化
[編集]他の関数空間上のフーリエ変換
[編集]フーリエ変換の性質は、対象とする関数が属する関数空間によって大きく異なる。数学的・物理的に重要な主要な空間における性質を以下に詳述する。
- 空間(絶対可積分関数)
- 最も基本的な定義は、絶対値の積分が有限である関数の空間 上でなされる。このとき、フーリエ変換された関数 は有界かつ一様連続であり、無限遠で0に収束する(リーマン・ルベーグの補題)。すなわち、 である。しかし、逆変換の積分が絶対収束するとは限らず、変換後の関数が再び に戻るとは限らないという非対称性がある。
- 空間(自乗可積分関数・エネルギー有限信号)
- 物理学や工学(信号処理、量子力学)において最も重要な空間である。信号のエネルギー(または粒子の存在確率)が有限であることを意味する。 関数は必ずしも ではない(例: )ため、単純な積分定義は通用しない場合がある。しかし、 が で稠密であることを利用し、極限操作によって定義を拡張できる。
- プランシュレルの定理により、 においてフーリエ変換はノルム(エネルギー)を保存するユニタリ変換(等長写像)となる。すなわち、 が成り立つ。また、フーリエ変換は から への全単射であり、情報の損失なく周波数領域への変換と復元が可能であることが保証される。
- 空間(急減少関数・シュワルツ空間)
- 理論的な解析、特に超関数の定義の基礎として用いられる、「無限回微分可能で、自身も導関数も多項式の逆数より速く減衰する関数」の空間。フーリエ変換は をそれ自身に移す全単射(自己同型)である。 や のような定義域と値域のズレが生じないため、微分などの演算とフーリエ変換の関係を議論するのに適している。
- 空間(緩増加超関数)
- シュワルツ空間の双対空間であり、通常の関数の枠組みには収まらない対象を含む。(インパルス)、定数関数 (直流成分)、正弦波 など、工学的に必須な対象のフーリエ変換を正当化するために導入される。フーリエ変換は 上の線形同相写像となり、一般化された意味での変換・逆変換が常に可能である。
- 空間()
- ハウスドルフ・ヤング不等式を用いて、フーリエ変換は から (ここで )への有界線型写像として定義される。なお、 の場合、フーリエ変換は通常の関数の意味では定義できず、超関数として扱う必要がある[3]。
多次元版
[編集]フーリエ変換は勝手な次元 n において考えることができる。1-次元の場合と同様にさまざまな流儀があるが、本項では絶対可積分函数 ƒ(x) に対して、
をフーリエ変換の定義とする。ここで、x および ξ は n-次元ベクトルであり、x · ξ はベクトルの点乗積である。点乗積はしばしば <x, ξ> とも書き表される。
プランシュレルの定理やパーセバルの定理がそうであるように、上述の基本性質は n-次元フーリエ変換においても成立する。函数が絶対可積分であるとき、フーリエ変換はやはり一様連続であり、リーマン・ルベーグの補題が成立する[3]。
より高い次元ではフーリエ変換の制限問題の研究が興味深いものになる。絶対可積分函数のフーリエ変換は連続で、この函数の任意の集合への制限が定義される。しかし自乗絶対可積分函数のフーリエ変換は自乗絶対可積分函数の一般の類を成す。そのような L2(Rn)-函数のフーリエ変換の制限は測度 0 の集合上では定義することができない。1 ≤ p ≤ 2 に対する Lp における制限問題の理解はいまだ活発な研究の行われる領域である。驚くべきことに、集合 S の曲率が非零であるようないくつかの場合には、フーリエ変換の S への制限を定義することができる。S が Rn における単位球面であるときが特に興味深い。この場合に、トマス-ステインの制限定理によれば、フーリエ変換の Rn における単位球面への制限は 1 ≤ p ≤ (2n + 2)/(n + 3) に対する Lp 上で有界作用素である。
1-次元の場合と多次元の場合とで、フーリエ変換の大きな違いは部分和作用素に関係する。与えられた絶対可積分函数 ƒ に対し
で定義される函数 ƒR を考える。さらに ƒ が Lp(Rn) に属すると仮定する。n = 1 で 1 < p < ∞ とし、SR = (−R, R) と置くと、ヒルベルト変換の有界性から ƒR は R を無限大に飛ばす極限で ƒ に Lp 内で収束する。素朴に n > 1 の場合にも同様であることを期待するかもしれない。SR を一辺の長さが R の立方体とするならば、確かに部分和作用素はもとの函数に収束する。別の自然な候補としてユークリッド球体 SR = {ξ : |ξ| < R} をとると、部分和作用素が収束するためには単位球体に対するマルチプライヤーが Lp(Rn) において有界である必要がある。n ≥ 2 に対しては、単位球体に対するマルチプライヤーは p = 2 でない限り有界にはならないというよく知られたチャールズ・フェファーマンの定理がある[7]。事実として、p ≠ 2 のときには ƒR が ƒ に Lp 内で収束しないだけではなく、函数 ƒ ∈ Lp(Rn) であっても ƒR が Lp の元でさえないようなものまでが存在する。
フーリエ・スティルチェス変換
[編集]Rn 上の有限ボレル測度 μ のフーリエ変換は
によって与えられる[1]。この変換は絶対可積分函数のフーリエ変換がもつ多くの性質を引き続き満足する。大きな違いの一つに、測度に関してリーマン・ルベーグの補題が成り立たないことが挙げられる[2]。dμ = ƒ(x)dx の場合には上述の定義式を f の通常のフーリエ変換の定義に簡約化することができる。
このフーリエ変換を用いて連続測度の特徴づけを与えることができる。ボホナーの定理はそのような函数を測度のフーリエ・スティルチェス変換として得られるものとして特徴付ける[2]。
さらに言えば、ディラックのデルタ函数は函数ではないが有限ボレル測度であり、そのフーリエ変換は定数函数となる(特殊値は用いるフーリエ変換の形に依存する)。
局所コンパクトアーベル群
[編集]フーリエ変換を任意の局所コンパクトアーベル群に対して一般化することができる。局所コンパクトアーベル群とは、抽象アーベル群であると同時に局所コンパクトなハウスドルフ空間であって、なおかつその位相に関して群演算が連続となるものである。G が局所コンパクトアーベル群ならば、G はハール測度と呼ばれる平行移動不変な測度 μ を持つ。また、局所コンパクトアーベル群 G に対して、その位相を指標全体の成す集合 へ移行することができて、 自身も局所コンパクトアーベル群の構造を持つ。L1(G) に属する函数 f に対して、そのフーリエ変換を
によって定義することができる[2]。
応用
[編集]フーリエ変換は、関数を周波数成分に分解するというその性質から、理論的な解析から実践的なデータ処理まで、科学・工学の極めて広い分野で応用されている。
微分方程式の解析学
[編集]フーリエ変換は、微分作用素を(周波数変数に関する)単純な掛け算に変換する性質を持つ。例えば、関数 の導関数 のフーリエ変換は となる(定義による定数の違いあり)。この性質を利用することで、線型微分方程式を、周波数領域における代数方程式に変換して解くことができる。
特に、熱伝導方程式や波動方程式、シュレーディンガー方程式のような偏微分方程式の解法において強力な道具となる。フーリエがこの変換を導入した直接の動機も、熱伝導方程式の解析であった。
信号処理
[編集]信号処理は、フーリエ変換の最も代表的な応用分野の一つである。
- スペクトル解析: 音声信号や電波などの時間信号 をフーリエ変換することで、その信号に含まれる周波数成分 (周波数スペクトル)を知ることができる。これにより、信号の特性分析、ノイズの同定、特定の周波数成分の検出などが可能となる。
画像処理
[編集]2次元フーリエ変換は、デジタル画像処理において広く用いられる。画像は空間的な輝度(または色)の分布であり、そのフーリエ変換は空間周波数(画像の模様の細かさや向き)のスペクトルを表す。
- 画像フィルタリング: 信号処理と同様に、周波数領域での操作により、画像のノイズ除去、エッジ強調、ぼかし処理などが行われる。
量子力学
[編集]量子力学において、粒子の状態は波動関数 で記述される。位置 に関する波動関数 と、運動量 に関する波動関数 は、互いにフーリエ変換(の定数倍)の関係にある。
この関係は、粒子の位置と運動量を同時に正確に決定することはできないという不確定性原理の数学的な表現(不確定性関係)と密接に関連している。
離散フーリエ変換
[編集]現実の信号処理や画像処理では、コンピュータを用いて数値的にフーリエ変換を計算する必要がある。この目的で用いられるのが、連続的な関数ではなく、離散的な(サンプリングされた)データ列を扱う離散フーリエ変換(英: Discrete Fourier Transform, DFT)である。
DFTは、有限個のデータ点 から、有限個の周波数成分 を計算する変換であり、その計算を劇的に高速化するアルゴリズムが高速フーリエ変換(英: Fast Fourier Transform, FFT)である。
FFTの登場により、膨大なデータに対しても実用的な時間でフーリエ変換(およびその逆変換)が計算可能となり、上記のような信号処理や画像処理の応用が爆発的に普及する基盤となった。
その他の記法
[編集]フーリエ変換の記法として (ξ) 以外によく用いられるものに
などがある。あるいはもっと他の記号を使うことも在りうる。たとえば、(f(x) と F(ξ) のように)もとの函数を表している文字の対応する大文字を用いてそのフーリエ変換を表すことは自然科学や工学においてとくによく用いられる記法である。
複素函数 (ξ) は、極座標に関してこれを表示することにより、振幅
および位相
と呼ばれるふたつの実函数 A(ξ) および φ(ξ) を用いて
なる形に解釈することができる。
このとき逆変換は ƒ(x) の周波数成分すべての再結合として
と書くことができる。各成分は振幅が A(ξ) で(x = 0 における)初期位相角が φ(ξ) であるような e2πixξ のかたちの複素正弦曲線である。
フーリエ変換は函数空間の間の写像として考えることもできる。この写像はここでは で表し、函数 f のフーリエ変換には が用いられる。この写像 は函数空間上の線型変換とみることができ、それによって と書く代わりに、ベクトル(ここでは函数 f)の線型変換を表す線型代数学の標準的な記法で と書くこともできる。函数にフーリエ変換を施した結果は再び函数となるから、この新たな函数の ξ における値というものには意味があり、それを あるいは などと表す。前者の場合には はまず f に施されて、その後に得られた函数の ξ における値が評価されるものと暗黙に理解されているということに注意しなければならない。なお、混乱を防ぐため(そしてラプラス変換の記法と統一して)括弧の種類を変えてのように書く場合もある。
数学や多くの応用科学において、函数 f それ自身と函数 f の変数 x における値 f(x) とを峻別しなければならないことがしばしばある。このことが意味するのは、たとえば のような記法は、形式的には f の x における「値」のフーリエ変換と解釈できてしまうということである。このような不具合にもかかわらず、特定の函数あるいは特定の変数の函数を頻繁に変換しなければならないような場合には、このような記法はよく用いられる。たとえば
は矩形函数のフーリエ変換が sinc-函数であることを表すために用いられることがあり、またたとえば
はフーリエ変換のシフト性を表すのに用いられることがある。最後の例は、変換される函数 f をx0 のではなく x の函数であるという前提のもとでのみ正しいということに注意を要する。
主なフーリエ変換の一覧
[編集]以下にフーリエ変換の閉じた表示に関する表を掲げる(フーリエ変換はよく用いられる三種類を挙げてある)。函数 ƒ(x) , g(x), h(x) に対して、それらのフーリエ変換をそれぞれ , , で表す。
函数の関係式
[編集]以下の表におけるフーリエ変換は (Erdélyi 1954) あるいは (Kammler 2000) の付録に見つけることができる。
| もとの函数 | ユニタリ・周波に関するフーリエ変換 | ユニタリ・角周波に関するフーリエ変換 | 非ユニタリ・角周波に関するフーリエ変換 | 備考 | |
|---|---|---|---|---|---|
|
|
|
||||
| 101 | 線型性 | ||||
| 102 | 時間領域シフト | ||||
| 103 | 周波数領域シフト 102の双対 | ||||
| 104 | |a| が大きければ f(ax) は 0 の周りに集中し は平らに広がる | ||||
| 105 | ここで、 は、それぞれの列で考えているフーリエ変換を施した結果の、変数を x に取替えたものである。 | ||||
| 106 | |||||
| 107 | 106の双対 | ||||
| 108 | f ∗ g は f と g との畳み込みである。この公式は畳み込み定理と呼ばれる。 | ||||
| 109 | 108の双対 | ||||
| 110 | 純実偶関数 | はいずれも純実偶関数 | 正弦・余弦変換も参照 | ||
| 111 | 純実奇関数 | はいずれも純虚奇関数 | |||
自乗絶対可積分函数
[編集]以下の表におけるフーリエ変換は (Campbell & Foster 1948), (Erdélyi 1954) あるいは (Kammler 2000) の付録に見つけることができる。
| もとの函数 | ユニタリ・周波に関するフーリエ変換 | ユニタリ・角周波に関するフーリエ変換 | 非ユニタリ・角周波に関するフーリエ変換 | 備考 | |
|---|---|---|---|---|---|
|
|
|
|
|||
| 201 | 矩形波と標準化されたsinc関数でsinc関数はsinc(x) = sin(πx)/(πx)で表される | ||||
| 202 | 201の双対で矩形波は理想的なローパスフィルターである。sinc関数はそのようなフィルターの非因果波応答である。 | ||||
| 203 | tri(x)は三角形関数である。 | ||||
| 204 | 203の双対 | ||||
| 205 | u(x)はヘビサイドの単位ステップ関数であり、a>0 | ||||
| 206 | これが示すものは、ガウス関数exp(−αx2)でαを選んだ場合はユニタリフーリエ変換である。 Re(α)>0で積分可能である | ||||
| 207 | a>0である | ||||
| 208 |
|
|
|
関数Jn (x)は、n次の第1種ベッセル関数である。関数Un (x)は第2種チェビシェフ多項式である。下記315と316を参照 | |
| 209 | 双曲線正割は自分自身をフーリエ変換したものである |
超函数
[編集]以下の表におけるフーリエ変換は (Erdélyi 1954) あるいは (Kammler 2000) の付録に見つけることができる。
| もとの函数 | ユニタリ・周波に関するフーリエ変換 | ユニタリ・角周波に関するフーリエ変換 | 非ユニタリ・角周波に関するフーリエ変換 | 備考 | |
|---|---|---|---|---|---|
|
|
|
|
|||
| 301 | δ(ξ) はディラックのデルタ関数 | ||||
| 302 | 301の双対 | ||||
| 303 | 103と301より導かれる。 | ||||
| 304 | 101、303とオイラーの公式:より導かれる。 | ||||
| 305 | 101、303と より導かれる。 | ||||
| 306 | |||||
| 307 | |||||
| 308 | n は自然数、 δ(n )(ξ) はディラックのデルタ関数のn 階微分。107と301より導かれる。さらに101と組み合わせることで、任意の多項式を変換できる。 | ||||
| 309 | sgn(ξ) は符号関数。1/x は超関数ではないことに注意。シュワルツ関数に対してテストするときにコーシーの主値を使用する必要がある。この規則はヒルベルト変換を研究するとき有用である。 | ||||
| 310 | 309の一般化 | ||||
| 311 | |||||
| 312 | 309の双対。積分はコーシーの主値を考える。 | ||||
| 313 | u (x ) はヘヴィサイドの階段関数。101、301および312より導かれる。 | ||||
| 314 | この関数はくし型関数といわれる。302、102および、超関数として であることから導かれる。 | ||||
| 315 | J0 (x ) は0次の第1種ベッセル関数 | ||||
| 316 | 315の一般化。Jn (x ) はn 次の第1種ベッセル関数、Tn (x ) は第1種チェビシェフ多項式。 |
二変数函数
[編集]| もとの函数 | ユニタリ・周波に関するフーリエ変換 | ユニタリ・角周波に関するフーリエ変換 | 非ユニタリ・角周波に関するフーリエ変換 | 備考 | |
|---|---|---|---|---|---|
|
|
|
|
ξx , ξy , ωx , ωy , νx , νy は実変数。積分領域は全平面である。 | ||
| 401 | 両方のガウス関数は規格化されている必要はない。 | ||||
| 402 | 元の函数は circ(r ) = 1 (0≤r ≤1), and 0 (otherwise) で定義される。これはエアリー分布であり、1次の第1種ベッセル函数 J1 で表される[18]。 |
一般の n-変数函数
[編集]| もとの函数 | ユニタリ・周波に関するフーリエ変換 | ユニタリ・角周波に関するフーリエ変換 | 非ユニタリ・角周波に関するフーリエ変換 | 備考 | |
|---|---|---|---|---|---|
|
|
|
||||
| 501 | χ[0,1] は区間 [0, 1] の指示関数、Γ(x ) はガンマ関数、Jn /2+δ はn /2 + δ次の第1種ベッセル関数である。n = 2 およびδ = 0とすると402を得る[19]。 |
ラプラス変換との関係
[編集]フーリエ変換は、ラプラス変換と密接な関係にある。数学的な観点において、フーリエ変換はラプラス変換の特殊な場合、あるいは複素平面の虚軸上に制限したものと解釈されることが多いが、工学的な応用においては「扱う時間の範囲」と「信号の安定性」に関して明確な使い分けが存在する。
両側ラプラス変換との関係
[編集]連続時間信号 の両側ラプラス変換 は、複素変数 を用いて以下のように定義される。
すなわち、ラプラス変換とは、元の関数 に実数の指数関数 を乗じて減衰(あるいは増幅)させ、その結果をフーリエ変換したものに他ならない。特に ()の場合、これは通常のフーリエ変換と一致する。
片側ラプラス変換との関係
[編集]工学(制御工学や電気回路)で一般的に用いられる「ラプラス変換」は、積分範囲を に限定した片側ラプラス変換である。 これは、元の関数 に単位ステップ関数 を乗じて を切り捨てたものの両側ラプラス変換と等価である。
フーリエ変換の視点から見ると、片側ラプラス変換への移行は以下の二つの操作を同時に行っていることになる。
- 因果律の導入:信号に を乗じる。周波数領域では、元のスペクトルと のフーリエ変換()との畳み込みが行われ、スペクトルが複素平面上で平滑化される。
- 複素領域への拡張:周波数変数を実数 から複素数 へ拡張する。
この操作により、片側ラプラス変換は「スイッチを入れた瞬間()」からの挙動を記述するのに適しており、かつ微分方程式を代数方程式に変換する際の初期値の取り扱いが容易になるという利点を持つ。
関連項目
[編集]出典
[編集]- 1 2 3 4 5 6 7 Pinsky 2002.
- 1 2 3 4 Katznelson 1976.
- 1 2 3 4 Stein & Weiss 1971.
- ↑ Rudin 1987, p. 187.
- ↑ Rudin 1987, p. 186.
- ↑ Stein & Shakarchi 2003.
- 1 2 Duoandikoetxea 2001.
- ↑ Grafakos 2004.
- ↑ Fourier, Joseph (1807), Théorie de la propagation de la chaleur dans les solides
- ↑ Fourier, Joseph (1822), Théorie analytique de la chaleur
- ↑ Dirichlet, Peter Gustav Lejeune (1829), “Sur la convergence des séries trigonométriques qui servent à représenter une fonction arbitraire entre des limites données”, Journal für die reine und angewandte Mathematik
- ↑ Riemann, Bernhard (1867), “Über die Darstellbarkeit einer Function durch eine trigonometrische Reihe”, Abhandlungen der Königlichen Gesellschaft der Wissenschaften zu Göttingen
- ↑ Lebesgue, Henri (1904), Leçons sur l'intégration et la recherche des fonctions primitives, Paris: Gauthier-Villars
- ↑ Plancherel, Michel (1910), “Contribution à l'étude de la représentation d'une fonction arbitraire par des intégrales définies”, Rendiconti del Circolo Matematico di Palermo
- ↑ Wiener, Norbert (1933), The Fourier Integral and Certain of its Applications, Cambridge University Press
- ↑ Schwartz, Laurent (1950-1951), Théorie des distributions, Paris: Hermann
- ↑ Cooley, James W.; Tukey, John W. (1965), “An algorithm for the machine calculation of complex Fourier series”, Mathematics of Computation 19: 297-301, doi:10.1090/S0025-5718-1965-0178586-1
- ↑ Stein & Weiss 1971, Thm. IV.3.3.
- ↑ Stein & Weiss 1971, Thm. 4.13.
参考文献
[編集]- Bochner, S.; Chandrasekharan, K. (1949). Fourier Transforms. Princeton University Press
- Bracewell, R. N. (2000), The Fourier Transform and Its Applications (3rd ed.), Boston: McGraw-Hill.
- Campbell, George; Foster, Ronald (1948), Fourier Integrals for Practical Applications, New York: D. Van Nostrand Company, Inc..
- Duoandikoetxea, Javier (2001), Fourier Analysis, American Mathematical Society, ISBN 0-8218-2172-5.
- Dym, H; McKean, H (1985), Fourier Series and Integrals, Academic Press, ISBN 978-0122264511.
- Erdélyi, Arthur, ed. (1954), Tables of Integral Transforms, 1, New Your: McGraw-Hill
- Grafakos, Loukas (2004), Classical and Modern Fourier Analysis, Prentice-Hall, ISBN 0-13-035399-X.
- Hörmander, L. (1976), Linear Partial Differential Operators, Volume 1, Springer-Verlag, ISBN 978-3540006626.
- James, J.F. (2002), A Student's Guide to Fourier Transforms (2nd ed.), New York: Cambridge University Press, ISBN 0-521-00428-4.
- Kaiser, Gerald (1994), A Friendly Guide to Wavelets, Birkhäuser, ISBN 0-8176-3711-7
- Kammler, David (2000), A First Course in Fourier Analysis, Prentice Hall, ISBN 0-13-578782-3
- Katznelson, Yitzhak (1976), An introduction to Harmonic Analysis, Dover, ISBN 0-486-63331-4
- Bao Luong(2009): Fourier Analysis on Finite Abelian Groups, Birkhäuser, ISBN 978-0-8176-4916-6. ※ 有限アーベル群上のフーリエ解析
- Pinsky, Mark (2002), Introduction to Fourier Analysis and Wavelets, Brooks/Cole, ISBN 0-534-37660-6
- Polyanin, A. D.; Manzhirov, A. V. (1998), Handbook of Integral Equations, Boca Raton: CRC Press, ISBN 0-8493-2876-4.
- Dinakar Ramakrishnan,and Robert J. Valenza (1999): Fourier Analysis on Number Fields, Springer (GMT 186),ISBN 978-1-4757-3085-2. ※ 数論関連
- Rudin, Walter (1987), Real and Complex Analysis (Third ed.), Singapore: McGraw-Hill, ISBN 0-07-100276-6.
- Stein, Elias; Shakarchi, Rami 新井仁之、杉本充、高木啓行、千原浩之訳 (2003), Fourier Analysis: An introduction [フーリエ解析入門], Princeton University Press, ISBN 0-691-11384-X.
- Stein, Elias; Weiss, Guido (1971), Introduction to Fourier Analysis on Euclidean Spaces, Princeton, N.J.: Princeton University Press, ISBN 978-0-691-08078-9.
- Audrey Terras (1999): Fourier Analysis on Finite Groups and Applications, Cambridge University Press, ISBN 978-0-52145718-7. ※ 有限群上のフーリエ解析
- Wilson, R. G. (1995), Fourier Series and Optical Transform Techniques in Contemporary Optics, New York: Wiley, ISBN 0471303577.
- Yosida, K. (1968), Functional Analysis, Springer-Verlag, ISBN 3-540-58654-7.
関連図書
[編集]- 高橋洋一郎:「実関数とフーリエ解析」、岩波書店、ISBN 4-00-005457-0 (2006年7月7日).
外部リンク
[編集]- Fourier Series Applet (Tip: drag magnitude or phase dots up or down to change the wave form).
- Tables of Integral Transforms at EqWorld: The World of Mathematical Equations.
- Weisstein, Eric W. “Fourier Transform”. mathworld.wolfram.com (英語).
- Fourier Transform Module by John H. Mathews
- The DFT “à Pied”: Mastering The Fourier Transform in One Day at The DSP Dimension
- FFT in Python
- Fourier-transform (mathematics) - ブリタニカ百科事典
- 日本大百科全書(ニッポニカ)『フーリエ変換』 - コトバンク