コンテンツにスキップ

フーリエ変換

出典: フリー百科事典『ウィキペディア(Wikipedia)』
上は時間領域で表現された矩形関数f(t)(左)と、周波数領域で表現されたそのフーリエ変換(ω)(右)。(ω)Sinc関数である。下は時間遅れのある矩形関数 g(t) と、そのフーリエ変換 ĝ(ω)。 時間領域における平行移動 (ディレイ)は、周波数領域では虚数部の位相シフトとして表現される。

数学においてフーリエ変換(フーリエへんかん、: Fourier transform、FT)は、変数複素または数値関数を、その周波数成分を表す別の関数ˆfに写す変換である。

工学においては、変換後の関数ˆfはもとの関数に含まれる周波数を記述していると考え、しばしばもとの関数周波数領域表現 (frequency domain representation) と呼ばれる。言い換えれば、フーリエ変換は関数を(離散的な正弦波・余弦波の和であるフーリエ級数とは対照的に)連続的な周波数成分の集まりとして分解する操作である。

フーリエ変換 (FT) は他の多くの数学的な演算と同様にフーリエ解析の主題を成す。「フーリエ変換」という言葉は関数の周波数領域表現のことを指すこともあるし、関数を周波数領域表現へ写す変換の過程・公式を言うこともある。なおこの呼称は、19世紀フランスの数学者・物理学者で次元解析の創始者とされるジョゼフ・フーリエに由来する。

概要

[編集]

フーリエ変換は、関数をその周波数成分に分解し、どの周波数がどれだけ含まれているかを示す「周波数スペクトル」に変換する数学的な操作である。この変換は、19世紀初頭にフランスの数学者・物理学者であるジョゼフ・フーリエが、熱伝導方程式の研究において、熱の伝わり方をフーリエ級数を用いて表そうとしたことに端を発する。

時間や空間に関する関数(波形など)を、その関数を構成する正弦波(または複素指数関数)の集まりとして表現し直す操作と見なすことができる。元の関数が時間領域(または空間領域)の関数 であれば、フーリエ変換後の関数 は周波数領域(または波数領域)の関数となり、元の関数 に周波数 の成分がどれだけ含まれているかを表す。

この変換は可逆であり、フーリエ逆変換と呼ばれる操作によって、周波数スペクトル から元の関数 を復元することができる。

フーリエ変換は、数学におけるフーリエ解析の中心的な道具であるだけでなく、物理学、工学(特に信号処理画像処理通信工学)、量子力学など、現代科学の非常に広範な分野で不可欠なツールとなっている。

定義

[編集]

フーリエ変換の定義には、用いる周波数の変数(周波数 [ヘルツ] か、角周波数 [ラジアン/秒] か)や、変換・逆変換の係数(正規化係数)の選び方によって、いくつかの流儀が存在する。分野によって慣習的に用いられる定義が異なるため注意が必要である。

主な定義として以下の3つが挙げられる。

よく用いられる定義のまとめ
周波数 ξ(ヘルツ) ユニタリ

角周波数 ω(ラジアン毎秒) 非ユニタリ

ユニタリ

本項では、特に断りのない限り、数学や一部の信号処理の文献で用いられることのある、周波数 (ヘルツ)を用い、変換・逆変換の係数が である定義1(ユニタリ・周波数)を採用する。この定義は、フーリエ変換と逆変換が対称形( の肩の符号が異なるだけ)である点に特徴がある。 また、この定義(および定義3)が「ユニタリ(: unitary)」と呼ばれるのは、プランシュレルの定理で示されるように、関数のノルム(信号のエネルギーに相当)が変換前後で保存される、すなわち が成り立つためであり、フーリエ変換が(空間における)ユニタリ作用素であることを意味する。定義2(非ユニタリ)では、このような単純なエネルギー保存則は成立しない(係数 が必要となる)。

絶対可積分関数に対する定義

[編集]

絶対可積分関数 f: RC のフーリエ変換は、

で定義される。ここでギリシャ文字小文字の ξ は任意の実数(周波数)である。

対象の関数における独立変数が物理量の場合、フーリエ変換は独立変数の次元をもとの逆数に移す。例えば、変換前の関数における独立変数 x時間の次元をもつとき、変換後の独立変数 ξ周波数の次元を持つ。あるいは、変換前の独立変数 x長さの次元をもつとき、変換後の独立変数 ξ波数の次元を持つ。この性質は定義より x ξ無次元量であることから従う。

適当な条件のもと、f はその変換 ˆf からフーリエ逆変換 (inverse transform)

によって復元することができる(x は任意の実数)。

緩増加超関数への拡張

[編集]

絶対可積分でない関数、例えば (定数関数)や (周期関数)、あるいはディラックのデルタ関数 のような「関数」ではない対象についてもフーリエ変換を定義したいという強い要請が物理学や工学には存在する。

このため、フーリエ変換の定義を緩増加超関数(またはシュワルツ超関数)の空間 に拡張することが行われる。

緩増加超関数 は、シュワルツ空間 (急減少関数、すなわち無限回微分可能で、自身およびその全ての導関数が で多項式の逆数よりも速く 0 に近づく滑らかな関数 の集まり)上の連続線型汎関数として定義される。つまり、 は各 に対して複素数 を対応させる写像である。

(絶対可積分関数 も、 により、 という超関数とみなすことができる。)

シュワルツ空間 はフーリエ変換 の下で閉じている(急減少関数のフーリエ変換は再び急減少関数となる)という重要な性質を持つ。

この性質を利用し、緩増加超関数 のフーリエ変換 は、任意の急減少関数 に対して、以下の関係(プランシュレルの定理の一般化)を満たす超関数 として定義される。

ここで は急減少関数 の(絶対可積分関数としての)フーリエ変換である。この定義により、超関数のフーリエ変換は再び緩増加超関数となる。

例:定数関数 を超関数 )とみなす。そのフーリエ変換 は、定義より となる。 ここで、 のフーリエ逆変換の定義 において とおくと、 という関係が得られる(この関係は、急減少関数 に対して成立する)。 したがって、 である。 一方、 という超関数 は、定義により である。 したがって、、すなわち定数関数 のフーリエ変換は となる。

例:デルタ関数 )のフーリエ変換 は、 これは、 という定数関数が定める超関数 の作用 に等しい。 したがって、デルタ関数 のフーリエ変換は定数関数 となる。

導入

[編集]

この節では、フーリエ変換の「動機」の一つとして、フーリエ級数からの類推的な導入を解説する。この方法は、フーリエ変換がどのようにして関数を周波数成分の「連続的な集まり」として表現し直すかという直観的なイメージを掴む上で、歴史的にも教育的にも重要である。

フーリエ変換を考える動機はフーリエ級数の研究に始まる。フーリエ級数の研究において、複雑な周期関数は単純な波動の数学的な表現である正弦関数余弦関数の和として表される。正弦や余弦の性質のおかげで、この和に現れる各波の量、フーリエ係数を積分によって計算することができる。

多くの場合に、オイラーの公式)を用いて、正弦関数および余弦関数の代りに基本波動 を用いた方が便利である。この場合には多くの公式が簡単化され、本項で後述するフーリエ変換のほかの類似の定式化をあたえるという点に優位性がある。この正弦・余弦から複素指数関数への移行にはフーリエ係数が複素数値であることを要する。この複素数は、関数に含まれる波動の振幅(あるいは大きさ)と、位相(あるいは初期角)の両方を与えているものと通常は解釈される。また、この移行に際して「負の周波数」も導入される。例えば、波動 および はともに周期1を持つが、複素フーリエ級数においては別々の成分として取り扱われる。したがって、周波数を単純に周期の逆数と考えることはできなくなる。

フーリエ級数を以下のようにしてフーリエ変換の動機付けに用いることができる。関数 ƒ をある区間 [−L/2, L/2] の外側で 0 となるようなものとすると、任意の TL に対して ƒ を区間 [−T /2, T /2] 上のフーリエ級数に拡張できる。ここで f のフーリエ級数に現れる波動 の係数となる で表される「量」は

で与えられ、ƒ は公式

で与えられなければならない。ξn = n/T とおき、Δξ = (n + 1)/Tn/T = 1/T とおくと、最後の和をリーマン和

として考えることができる。T → ∞ とすることにより、このリーマン和は定義節で与えられるフーリエ逆変換(積分)に収束すると期待される。このように、フーリエ変換は、関数に含まれる個々の周波数がどの程度あるかを(スペクトル密度 として)測るものと考えることができ、それらの波動を積分(「連続和」)によって再結合して元の関数を復元する操作がフーリエ逆変換に対応する。

以下の画像はフーリエ変換が特定の関数に含まれる周波数を測る方法を視覚的に現したものである。関数として、(t が秒で測られる場合には)3 ヘルツで振動し、急速に 0 になる

を描く。この関数は特に描画しやすい実フーリエ変換をもつものとして選ばれたものであり、最初の画像はそのグラフである。ˆf(3) を計算するために、e−2πi(3t)ƒ(t) を積分する。二枚目の画像はこの被積分関数の実部および虚部である。被積分関数の実部は殆ど常に正となる。これは ƒ(t) が負であるときには e−2πi(3t) の実部が同様に負となることによる。それらは同じ比率で振動するから、ƒ(t) が正であるときも同様に e−2πi(3t) の実部も正になる。

この結果、被積分関数の実部のを積分すれば、比較的大きな数値(ここでの場合 0.5)を得ることになる。

一方、(ˆf(5) を見る場合のように)含まれない周波数を測れば、被積分関数は十分に振動し、それゆえにその積分はとても小さい値となる。一般の設定ではこれよりは少し複雑になるが、それでもフーリエ変換は関数 ƒ(t) に含まれる個々の周波数がどれくらいあるかを測るものという考え方に変わりはない。

この例では、ではなく、と、不自然ながかかっているが、は、(−∞, ∞)で絶対可積分関数としてのフーリエ変換(積分)が定義できないためである。(前述の通り、超関数としては定義可能である。)一方、は、のため(−∞, ∞)で絶対可積分であり、 が成り立つ。よって、 は全体として絶対可積分関数となり、フーリエ変換(積分)が可能となっているものである。

フーリエ変換の性質

[編集]

実数直線上で定義される関数 f絶対可積分であるとは、

を満たすルベーグ可測関数であることをいう。

基本性質

[編集]

絶対可積分関数 f(x), g(x), h(x) が与えられたとき、これらのフーリエ変換をそれぞれ ˆf(ξ), ˆg(ξ), ˆh(ξ)で表す。フーリエ変換は以下の基本性質を満たす[1]

線型性

任意の複素数 a, b について h(x) = (x) + bg(x) であるならば

が成り立つ。

平行移動

任意の実数 x0 に対して h(x) = ƒ(xx0) であるならば

が成り立つ。

変調

任意の実数 ξ0 に対して h(x) = eixξ0ƒ(x) ならば

が成り立つ。

定数倍

非零実数 a に対し、h(x) = ƒ(ax) ならば

が成り立つ。a = −1 つまり h(x) = ƒ(−x) の場合には、時間反転性 (time-reversal property)

が導かれる。

複素共役

f(x)複素共役 f(x) について

が成り立つ。

微分

関数 が( で 0 に収束するなどの条件を満たし) 回微分可能で、その導関数も絶対可積分であるならば、 が成り立つ。逆に、 も成り立つ。この性質は、微分方程式を代数方程式に変換する上で重要である。

畳み込み

h(x) = (fg)(x) ならば

が成り立つ。

一様連続性とリーマン・ルベーグの補題

[編集]

絶対可積分関数のフーリエ変換は、常に成り立つというわけではない性質も持っている。絶対可積分関数 ƒ のフーリエ変換は一様連続で

を満たす[2]。絶対可積分関数のフーリエ変換は

であることを述べたリーマン・ルベーグの補題をも満足する[3]。絶対可積分函数 f のフーリエ変換 ˆf は有界連続だが絶対可積分であるとは限らず、その逆変換をルベーグ積分として書くことは一般にはできない。しかしながら、ƒ および ˆf がともに絶対可積分ならば、反転公式

が殆ど全ての x において成り立つ。つまり、ƒ は右辺で定義される連続関数と殆ど至る所等しい。特に ƒ が実数直線上の連続関数として与えられたならば全ての x において等式が成り立つ。

前述の結果としてわかることは、フーリエ変換が L1(R) 上単射であることである。

プランシュレルの定理とパーセバルの定理

[編集]

f(x) および g(x) は絶対可積分であるとし、そのフーリエ変換をそれぞれ ˆf(ξ) および ˆg(ξ) と表す。f(x) および g(x) がともに自乗絶対可積分であるならばパーセバルの定理

が成立する[4]。ここで上付きバーは複素共役を表す。

パーセバルの定理と同値なプランシュレルの定理によれば

が成立する[5]。プランシュレルの定理により、L2(R) に属する関数の後述する意味でのフーリエ変換を定義することが可能になる。プランシュレルの定理は、フーリエ変換はもとの量のエネルギーを保存するという自然科学における解釈を持つ。著者によってはこれらの定理のどちらともをプランシュレルの定理あるいはパーセバルの定理と呼んでいる場合があるので注意を要する。

局所コンパクトアーベル群に関する文脈におけるフーリエ変換の概念の一般の定式化についてはポントリャーギン双対の項を参照されたい。

不確定性関係

[編集]

一般的に言って、f(x) が凝縮されればされるほどそのフーリエ変換 ˆf はより拡散される。特に、フーリエ変換のスケール性からわかることとして、関数を x において「圧搾」するならば、そのフーリエ変換は ξ において「伸展」される。したがって、関数とそのフーリエ変換の両方ともを勝手に凝縮させることはできない。

関数とそのフーリエ変換のコンパクト化のあいだの得失評価は不確定性関係の形で定式化することができる。ƒ(x) は絶対可積分かつ自乗絶対可積分であると仮定する。一般性を失うことなく関数 ƒ(x) は

に正規化されているものと仮定してよい。このとき、プランシュレルの定理により ˆf(ξ) も同様に正規化される。

x = 0 の周りでの拡散を

で定義される「0 の周りでの分散」(dispersion about zero) によって測ることにする[1]。確率の言葉で言えば、これは |f(x)|2 の 0 の周りでの二次のモーメントである。

このとき不確定性原理は、関数 ƒ(x) が絶対連続で、関数 x·ƒ(x) および ƒ′(x) が自乗絶対可積分であるならば

が成り立つことを述べる[1]。等式が成立するのは

したがって、

である場合に限る。ただし、定数 σ > 0 は任意であり、係数 C1ƒL2-正規化する定数である[1]。言い換えれば、 ƒ は 0 を中心に持つ(正規化)ガウス関数のとき等号が成り立つ。

事実として、この不等式は任意の x0, ξ0R について

が成立することをも含む[6]

量子力学において、運動量位置波動関数は(プランク定数を因子に持つ)フーリエ変換対である。プランク定数でスケールしなおせば、上述の不等式はロバートソンの不確定性関係を記述する。これは、ハイゼンベルグが構想した不確定性原理そのものではないが、深い関係がある。

ポアソン和公式

[編集]

ポアソン和公式はフーリエ変換とフーリエ級数の間の関連性を提供する。絶対可積分関数 ƒL1(Rn) が与えられたとき、ƒ周期化

によって与えられる。このとき、ポアソン和公式は f のフーリエ級数を ƒ のフーリエ変換に結びつけるもので、特に f のフーリエ級数は

で与えられることを述べるものである。ポアソン和公式を用いて、大きな次元のユークリッド球面における格子点の数に対するランダウの漸近公式を導出することができる。また、絶対可積分函数 fˆf がともにコンパクト台を持つならば ƒ = 0 を示すこともできる[1]

畳み込み定理

[編集]

フーリエ変換は、関数の畳み込みと関数の(点毎の)積とを相互に変換する。ƒ(x) および g(x) が絶対可積分関数であるとし、そのフーリエ変換をそれぞれ ˆf(ξ) および ˆg(ξ) で表す。さらに ƒg との畳み込みが存在して絶対絶対可積分であるならば、この畳み込みのフーリエ変換はフーリエ変換 ˆf(ξ)ˆg(ξ) との積で与えられる(ただし、フーリエ変換の定義の仕方によっては定数因子が現れる場合もある)。

これを式で表せば、∗ を畳み込みとして

と表されるとき、

が成立することを意味する。線型時不変 (LTI) 系理論において、f(x) を単位インパルスで置き換えたものが h(x) = g(x) を与えることから、通例 g(x) は、入力 ƒ(x) と出力 h(x) に関する LTI 系のインパルス応答として解釈される。この場合、ˆg(ξ) はこの系の周波数応答を表す。

逆に、ƒ(x) がふたつの自乗絶対可積分函数 p(x) および q(x) の積に分解されるならば、 ƒ(x) のフーリエ変換は、各因子のフーリエ変換 ˆp(ξ) および ˆq(ξ) の畳み込みで与えられる。

相互相関定理

[編集]

同様の方法で、h(x) が ƒ(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]ƒ(x) = ƒ0(|x|)P(x) (P(x) ∈ ) と表される関数のフーリエ変換は

となる。ただし、

であり、J(n + 2k − 2)/2 は次数 (n + 2k − 2)/2 の第一種ベッセル関数である。k = 0 のとき、これは動径関数のフーリエ変換に対する有用な公式を与える[8]

一般化

[編集]

他の函数空間上のフーリエ変換

[編集]

フーリエ変換の定義を他の函数空間に対するものへ拡張することができる。コンパクト台を持つ滑らかな函数は絶対可積分で、その全体は L2(R) において稠密であるから、プランシュレルの定理を用いて、L2(R) の一般の函数にまで(コンパクト台をもつ滑らかな函数によって近似して)フーリエ変換の定義を拡張することができる。さらに

ユニタリ作用素である[9]。フーリエ変換の多くの性質はこの場合にもそのまま成立する。ハウスドルフ・ヤング不等式を用いて 1 ≤ p ≤ 2 に対する Lp(R) の函数を含むようにフーリエ変換の定義を拡張することができる。

だが、さらなる拡張はもっと技巧的である。2 < p < ∞ の範囲でのLp に属する函数のフーリエ変換には超函数の研究が必要である[2]。事実として、p > 2 に関する Lp に属する函数のフーリエ変換は函数としては定義できないことを示すことができる [3]

多次元版

[編集]

フーリエ変換は勝手な次元 n において考えることができる。1-次元の場合と同様にさまざまな流儀があるが、本項では絶対可積分函数 ƒ(x) に対して、

をフーリエ変換の定義とする。ここで、x および ξn-次元ベクトルであり、x · ξ はベクトルの点乗積である。点乗積はしばしば <x, ξ> とも書き表される。

プランシュレルの定理やパーセバルの定理がそうであるように、上述の基本性質は n-次元フーリエ変換においても成立する。函数が絶対可積分であるとき、フーリエ変換はやはり一様連続であり、リーマン・ルベーグの補題が成立する[3]

より高い次元ではフーリエ変換の制限問題の研究が興味深いものになる。絶対可積分函数のフーリエ変換は連続で、この函数の任意の集合への制限が定義される。しかし自乗絶対可積分函数のフーリエ変換は自乗絶対可積分函数の一般の類を成す。そのような L2(Rn)-函数のフーリエ変換の制限は測度 0 の集合上では定義することができない。1 ≤ p ≤ 2 に対する Lp における制限問題の理解はいまだ活発な研究の行われる領域である。驚くべきことに、集合 S の曲率が非零であるようないくつかの場合には、フーリエ変換の S への制限を定義することができる。SRn における単位球面であるときが特に興味深い。この場合に、トマス-ステインの制限定理によれば、フーリエ変換の Rn における単位球面への制限は 1 ≤ p ≤ (2n + 2)/(n + 3) に対する Lp 上で有界作用素である。

1-次元の場合と多次元の場合とで、フーリエ変換の大きな違いは部分和作用素に関係する。与えられた絶対可積分函数 ƒ に対し

で定義される函数 ƒR を考える。さらに ƒLp(Rn) に属すると仮定する。n = 1 で 1 < p < ∞ とし、SR = (−R, R) と置くと、ヒルベルト変換の有界性から ƒRR を無限大に飛ばす極限で ƒLp 内で収束する。素朴に n > 1 の場合にも同様であることを期待するかもしれない。SR を一辺の長さが R の立方体とするならば、確かに部分和作用素はもとの函数に収束する。別の自然な候補としてユークリッド球体 SR = {ξ : |ξ| < R} をとると、部分和作用素が収束するためには単位球体に対するマルチプライヤーが Lp(Rn) において有界である必要がある。n ≥ 2 に対しては、単位球体に対するマルチプライヤーは p = 2 でない限り有界にはならないというよく知られたチャールズ・フェファーマンの定理がある[7]。事実として、p ≠ 2 のときには ƒRƒLp 内で収束しないだけではなく、函数 ƒLp(Rn) であっても ƒRLp の元でさえないようなものまでが存在する。

フーリエ・スティルチェス変換

[編集]

Rn 上の有限ボレル測度 μ のフーリエ変換は

によって与えられる[1]。この変換は絶対可積分函数のフーリエ変換がもつ多くの性質を引き続き満足する。大きな違いの一つに、測度に関してリーマン・ルベーグの補題が成り立たないことが挙げられる[2]dμ = ƒ(x)dx の場合には上述の定義式を f の通常のフーリエ変換の定義に簡約化することができる。

このフーリエ変換を用いて連続測度の特徴づけを与えることができる。ボホナーの定理はそのような函数を測度のフーリエ・スティルチェス変換として得られるものとして特徴付ける[2]

さらに言えば、ディラックのデルタ函数は函数ではないが有限ボレル測度であり、そのフーリエ変換は定数函数となる(特殊値は用いるフーリエ変換の形に依存する)。

局所コンパクトアーベル群

[編集]

フーリエ変換を任意の局所コンパクトアーベル群に対して一般化することができる。局所コンパクトアーベル群とは、抽象アーベル群であると同時に局所コンパクトハウスドルフ空間であって、なおかつその位相に関して群演算が連続となるものである。G が局所コンパクトアーベル群ならば、Gハール測度と呼ばれる平行移動不変な測度 μ を持つ。また、局所コンパクトアーベル群 G に対して、その位相を指標全体の成す集合 ˆG へ移行することができて、ˆG 自身も局所コンパクトアーベル群の構造を持つ。L1(G) に属する函数 f に対して、そのフーリエ変換を

によって定義することができる[2]

この一般化を概周期函数に適用した理論や、準周期函数に適用した理論が知られている。

応用

[編集]

フーリエ変換は、関数を周波数成分に分解するというその性質から、理論的な解析から実践的なデータ処理まで、科学・工学の極めて広い分野で応用されている。

微分方程式の解析学

[編集]

フーリエ変換は、微分作用素を(周波数変数に関する)単純な掛け算に変換する性質を持つ。例えば、関数 の導関数 のフーリエ変換は となる(定義による定数の違いあり)。この性質を利用することで、線型微分方程式を、周波数領域における代数方程式に変換して解くことができる。

特に、熱伝導方程式波動方程式シュレーディンガー方程式のような偏微分方程式の解法において強力な道具となる。フーリエがこの変換を導入した直接の動機も、熱伝導方程式の解析であった。

信号処理

[編集]

信号処理は、フーリエ変換の最も代表的な応用分野の一つである。

  • スペクトル解析: 音声信号や電波などの時間信号 をフーリエ変換することで、その信号に含まれる周波数成分 (周波数スペクトル)を知ることができる。これにより、信号の特性分析、ノイズの同定、特定の周波数成分の検出などが可能となる。
  • フィルタリング: 信号から特定の周波数成分を除去(例:ノイズ除去)したり、特定の周波数成分のみを通過させたりする操作(フィルタリング)は、周波数領域で容易に行うことができる。例えば、不要な高周波成分を除去するローパスフィルタは、周波数領域で に特定の関数(高周波域で0になる関数)を掛け合わせ、それを逆フーリエ変換することで実現される。これは、畳み込み定理に基づいている。
  • 非定常信号の解析: 実際の音声や音楽のように、時間と共に周波数成分が変化する信号(非定常信号)を解析するためには、短時間フーリエ変換 (STFT) やウェーブレット変換が用いられる。これらは、時間領域の局所的な区間(窓関数を適用した区間)に対してフーリエ変換を繰り返し適用する手法である。

画像処理

[編集]

2次元フーリエ変換は、デジタル画像処理において広く用いられる。画像は空間的な輝度(または色)の分布であり、そのフーリエ変換は空間周波数(画像の模様の細かさや向き)のスペクトルを表す。

  • 画像フィルタリング: 信号処理と同様に、周波数領域での操作により、画像のノイズ除去、エッジ強調、ぼかし処理などが行われる。
  • 画像圧縮: JPEGなどの画像圧縮規格は、画像を小さなブロックに分割し、離散コサイン変換(フーリエ変換の変種)を用いて空間周波数成分に変換し、人間の視覚が敏感でない高周波成分の情報を削減することでデータ量を圧縮している。

量子力学

[編集]

量子力学において、粒子の状態は波動関数 で記述される。位置 に関する波動関数 と、運動量 に関する波動関数 は、互いにフーリエ変換(の定数倍)の関係にある。

この関係は、粒子の位置と運動量を同時に正確に決定することはできないという不確定性原理の数学的な表現(不確定性関係)と密接に関連している。

離散フーリエ変換

[編集]

現実の信号処理や画像処理では、コンピュータを用いて数値的にフーリエ変換を計算する必要がある。この目的で用いられるのが、連続的な関数ではなく、離散的な(サンプリングされた)データ列を扱う離散フーリエ変換: Discrete Fourier Transform, DFT)である。

DFTは、有限個のデータ点 から、有限個の周波数成分 を計算する変換であり、その計算を劇的に高速化するアルゴリズムが高速フーリエ変換: Fast Fourier Transform, FFT)である。

FFTの登場により、膨大なデータに対しても実用的な時間でフーリエ変換(およびその逆変換)が計算可能となり、上記のような信号処理や画像処理の応用が爆発的に普及する基盤となった。

フーリエ変換の定義域と値域のまとめ

[編集]

フーリエ変換は、定義する関数の空間(定義域)によって、変換後の関数の属する空間(値域)が異なる。主な空間と変換の性質を以下にまとめる。

  • シュワルツ空間 : 急減少関数の空間。フーリエ変換は から への全単射であり、この空間で閉じている。(詳細は超関数への拡張を参照)
  • 空間 : 絶対可積分関数の空間。フーリエ変換 から (無限遠で 0 に収束する連続関数の空間)への有界線型写像である(リーマン・ルベーグの補題)。
  • 空間 : 自乗可積分関数の空間。フーリエ変換は から へのユニタリ作用素である(プランシュレルの定理)。ただし、この空間での変換は、必ずしも積分として定義されるとは限らず、極限操作によって定義される。(詳細は一般化を参照)
  • 空間): フーリエ変換は から への有界線型写像である()。(ハウスドルフ・ヤング不等式
  • 緩増加超関数 : シュワルツ空間の双対空間。フーリエ変換は から への全単射であり、閉じている。これにより、デルタ関数や定数関数など、 に属さない対象のフーリエ変換も定義可能となる。(詳細は超関数への拡張を参照)

その他の記法

[編集]

フーリエ変換の記法として ˆf(ξ) 以外によく用いられるものに

などがある。あるいはもっと他の記号を使うことも在りうる。たとえば、(f(x) と F(ξ) のように)もとの函数を表している文字の対応する大文字を用いてそのフーリエ変換を表すことは自然科学や工学においてとくによく用いられる記法である。

複素函数 ˆf(ξ) は、極座標に関してこれを表示することにより、振幅

および位相

と呼ばれるふたつの実函数 A(ξ) および φ(ξ) を用いて

なる形に解釈することができる。

このとき逆変換は ƒ(x) の周波数成分すべての再結合として

と書くことができる。各成分は振幅A(ξ) で(x = 0 における)初期位相角が φ(ξ) であるような eixξ のかたちの複素正弦曲線である。

フーリエ変換は函数空間の間の写像として考えることもできる。この写像はここでは で表し、函数 f のフーリエ変換には が用いられる。この写像 は函数空間上の線型変換とみることができ、それによって と書く代わりに、ベクトル(ここでは函数 f)の線型変換を表す線型代数学の標準的な記法で と書くこともできる。函数にフーリエ変換を施した結果は再び函数となるから、この新たな函数の ξ における値というものには意味があり、それを あるいは などと表す。前者の場合には はまず f に施されて、その後に得られた函数の ξ における値が評価されるものと暗黙に理解されているということに注意しなければならない。なお、混乱を防ぐため(そしてラプラス変換の記法と統一して)括弧の種類を変えてのように書く場合もある。

数学や多くの応用科学において、函数 f それ自身と函数 f の変数 x における値 f(x) とを峻別しなければならないことがしばしばある。このことが意味するのは、たとえば のような記法は、形式的には fx における「値」のフーリエ変換と解釈できてしまうということである。このような不具合にもかかわらず、特定の函数あるいは特定の変数の函数を頻繁に変換しなければならないような場合には、このような記法はよく用いられる。たとえば

は矩形函数のフーリエ変換が sinc-函数であることを表すために用いられることがあり、またたとえば

はフーリエ変換のシフト性を表すのに用いられることがある。最後の例は、変換される函数 fx0 のではなく x の函数であるという前提のもとでのみ正しいということに注意を要する。

主なフーリエ変換の一覧

[編集]

以下にフーリエ変換の閉じた表示に関する表を掲げる(フーリエ変換はよく用いられる三種類を挙げてある)。函数 ƒ(x) , g(x), h(x) に対して、それらのフーリエ変換をそれぞれ ˆf, ˆg, ˆh で表す。

函数の関係式

[編集]

以下の表におけるフーリエ変換は (Erdélyi 1954) あるいは (Kammler 2000) の付録に見つけることができる。

もとの函数 ユニタリ・周波に関するフーリエ変換 ユニタリ・角周波に関するフーリエ変換 非ユニタリ・角周波に関するフーリエ変換 備考

101 線型性
102 時間領域シフト
103 周波数領域シフト
102の双対
104 |a| が大きければ f(ax) は 0 の周りに集中し は平らに広がる
105 ここで、 は、それぞれの列で考えているフーリエ変換を施した結果の、変数を x に取替えたものである。
106
107 106の双対
108 fgfg との畳み込みである。この公式は畳み込み定理と呼ばれる。
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 で表される[10]

一般の n-変数函数

[編集]
もとの函数 ユニタリ・周波に関するフーリエ変換 ユニタリ・角周波に関するフーリエ変換 非ユニタリ・角周波に関するフーリエ変換 備考

501


χ[0,1] は区間 [0, 1] の指示関数、Γ(x ) はガンマ関数、Jn /2+δn /2 + δ次の第1種ベッセル関数である。n = 2 およびδ = 0とすると402を得る[11]

関連項目

[編集]

出典

[編集]

参考文献

[編集]
  • 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日).

外部リンク

[編集]