数学および信号処理におけるヒルベルト変換(ヒルベルトへんかん、: Hilbert transform)は、実変数関数 u(t) を別の実変数関数 H(u)(t) へ写すある特定の線型作用素を言う。具体的にこの作用素は 1πt との畳み込み: で与えられる。ただし、現れる広義積分コーシー主値の意味でとる。このヒルベルト変換は周波数領域において特に単純な表現 —引数となる函数の各フーリエ成分に π/2(90°) の位相ずれ (phase shift) を生じさせる— を持つ。例えば、余弦函数 cos(ωt) (ω > 0) のヒルベルト変換は cos(ωtπ/2) となる。

信号処理におけるヒルベルト変換は、それが実数値信号 u(t)解析的表現英語版を導くという点において重要である。具体的に、u のヒルベルト変換を v とすれば、vu調和共軛英語版となる。すなわち、v は実変数 t の函数であって、複素数値函数 u+ivコーシー–リーマン方程式を満足するように複素上半平面まで延長可能となる。この設定でヒルベルト変換を最初に導入したのはダフィット・ヒルベルトで、解析函数に対するリーマン–ヒルベルト問題英語版の特別の場合を解決するためであった。

導入

編集

u のヒルベルト変換は u(t)コーシー核英語版と呼ばれる函数 h(t) ≔ 1/πt との畳み込みと考えることができる。h(t)可積分でないから畳み込みの定義積分は収束しないが、代わりにコーシー主値(ここでは p.v. と書く)を用いることでヒルベルト変換は定義される。陽に書けば、函数(あるいは信号)u(t) のヒルベルト変換は   で—上記の積分が主値の意味で存在する限りにおいて—定義される。これはちょうど、u緩増加超函数 p.v. 1/πt との畳み込みになっている[注 1]。あるいはまた、変数変換を施すことにより、上記の主値積分を   と明示的に書くことができる[2]

函数 u にヒルベルト変換 H を続けて二回施すとき、定義に現れる二度の積分が適当な意味で収束する限りにおいて、得られる結果は u の符号反転:   である。特に、H の逆変換は H になる。この事実を確認するには、u のフーリエ変換のもとでのヒルベルト変換の振る舞い(後述)を見ることがもっとも簡便である。

上半平面上の解析函数に対し、ヒルベルト変換は境界値の実部と虚部との間の関係を記述する。つまり、複素函数 f(z) が平面 ℑm z > 0 上で解析的かつ u(t) ≔ ℜe f(t + 0⋅i) と書くとき、u のヒルベルト変換が存在する限りにおいて ℑm f(t + 0⋅i) = H(u)(t) が定数を加える違いを除いて成立する。

記法について

編集

信号処理において、u(t) のヒルベルト変換は一般的に ˆu(t) と書かれる[3]。しかし、数学においてこの記法は、広く一般に u(t) のフーリエ変換を表すものとして既に用いられている[4]。稀に、ヒルベルト変換が ~u(t) と書かれているかもしれない。さらにいえば、本項で定義したのとは符号が逆のものをヒルベルト変換と定義する文献も多い[5]

フーリエ変換との関係

編集

ヒルベルト変換はフーリエ乗算作用素英語版(すなわち、フーリエ変換のもと定数を乗じる操作として働く作用素)である[6]。その乗率は、符号函数 sgn を用いて σH(ω) ≔ −i sgn(ω) で与えられる。すなわち、フーリエ変換を F と書くとき   が成り立つ。フーリエ変換 F は、—適当な定数を掛けるだけの違いだが— 異なる定義がよくもちいられるものでも三種類あるが、符号函数は引数を正数倍しても変わらず sgn(x) = sgn(2πx) が成り立つから、上記の結果はどの定義のフーリエ変換でも変わらず適用できる。

オイラーの公式を適用すれば、  と書けるから、したがって H(u)(t)u(t)負の周波数英語版成分に +90° (π/2 rad) および正の周波数成分に −90° の位相ずらし (phase shft) を引き起こす。また i⋅H(u)(t) は正の周波数成分をそのままに、負の周波数成分をさらに +90°(計 +180° のずれ、つまり符号反転)を引き起こす。

ヒルベルト変換を二回反復適用するとき、u の負および正の周波数成分は、それぞれ +180° および −180° ずれて、これらのずれは一致する。よって、信号自体は符号が反転する: H(H(u)) = −u。これは   による。

ヒルベルト変換表

編集

次の表では、周波数パラメータ は実数である。

信号  ヒルベルト変換[fn 1] 
正弦[fn 2]   
余弦[fn 2]   
   
   
   (ドーソン積分英語版を参照)
sinc関数    
矩形関数    
ディラックのデルタ関数    
特性関数    
Note
  1. ^ 文献によっては(例えば Bracewell[要文献特定詳細情報])、本項で言うと H にあたるものをヒルベルト順変換の定義とするものもある。その場合には、一覧表の右列はすべて符号が逆になる。
  2. ^ a b sincos のヒルベルト変換は無限遠点において積分の主値をとることで定義することができる(シュヴァルツ超函数の意味で定義したものとも一致する)。

幅広いヒルベルト変換の一覧表が利用可能 (King 2009b). 定数のヒルベルト変換は 0 となることに注意。

注釈

編集
  1. ^ Schwartz (1950) による[1]

出典

編集
  1. ^ Pandey 1996, Chapter 3.
  2. ^ Zygmund 1968, §XVI.1.
  3. ^ 例えば Brandwood 2003, p. 87
  4. ^ 例えば Stein & Weiss 1971
  5. ^ 例えば Bracewell 2000, p. 359
  6. ^ Duoandikoetxea 2000, Chapter 3.

参考文献

編集

関連項目

編集

外部リンク

編集