ヒルベルト変換

提供: testwiki
ナビゲーションに移動 検索に移動

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

信号処理におけるヒルベルト変換は、それが実数値信号 テンプレート:Mathテンプレート:Ill2を導くという点において重要である。具体的に、テンプレート:Mvar のヒルベルト変換を テンプレート:Mvar とすれば、テンプレート:Mvarテンプレート:Mvarテンプレート:Ill2となる。すなわち、テンプレート:Mvar は実変数 テンプレート:Mvar の函数であって、複素数値函数 テンプレート:Mathコーシー–リーマン方程式を満足するように複素上半平面まで延長可能となる。この設定でヒルベルト変換を最初に導入したのはダフィット・ヒルベルトで、解析函数に対するテンプレート:Ill2の特別の場合を解決するためであった。

導入

テンプレート:Mvar のヒルベルト変換は テンプレート:Mathテンプレート:Ill2と呼ばれる函数 テンプレート:Math との畳み込みと考えることができる。テンプレート:Math可積分でないから畳み込みの定義積分は収束しないが、代わりにコーシー主値(ここでは テンプレート:Math と書く)を用いることでヒルベルト変換は定義される。陽に書けば、函数(あるいは信号)テンプレート:Math のヒルベルト変換は H(u)(t):=1πp.v.+u(τ)tτdτ で—上記の積分が主値の意味で存在する限りにおいて—定義される。これはちょうど、テンプレート:Mvar緩増加超函数 テンプレート:Math との畳み込みになっているテンプレート:Efn2。あるいはまた、変数変換を施すことにより、上記の主値積分を H(u)(t)=2πlimε0εu(t+τ)u(tτ)2τdτ と明示的に書くことができるテンプレート:Sfn

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

上半平面上の解析函数に対し、ヒルベルト変換は境界値の実部と虚部との間の関係を記述する。つまり、複素函数 テンプレート:Math が平面 テンプレート:Math 上で解析的かつ テンプレート:Math と書くとき、テンプレート:Mvar のヒルベルト変換が存在する限りにおいて テンプレート:Math が定数を加える違いを除いて成立する。

記法について

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

フーリエ変換との関係

ヒルベルト変換はテンプレート:Ill2(すなわち、フーリエ変換のもと定数を乗じる操作として働く作用素)であるテンプレート:Sfn。その乗率は、符号函数 テンプレート:Math を用いて テンプレート:Math で与えられる。すなわち、フーリエ変換を テンプレート:Mathcal と書くとき (H(u))(ω)=(isgn(ω))(u)(ω) が成り立つ。フーリエ変換 テンプレート:Mathcal は、—適当な定数を掛けるだけの違いだが— 異なる定義がよくもちいられるものでも三種類あるが、符号函数は引数を正数倍しても変わらず テンプレート:Math が成り立つから、上記の結果はどの定義のフーリエ変換でも変わらず適用できる。

オイラーの公式を適用すれば、σH(ω)={i=e+iπ/2,for ω<00,for ω=0i=eiπ/2,for ω>0 と書けるから、したがって テンプレート:Mathテンプレート:Mathテンプレート:Ill2成分に テンプレート:Math (テンプレート:Math rad) および正の周波数成分に テンプレート:Math の位相ずらし (phase shft) を引き起こす。また テンプレート:Math は正の周波数成分をそのままに、負の周波数成分をさらに テンプレート:Math(計 テンプレート:Math のずれ、つまり符号反転)を引き起こす。

ヒルベルト変換を二回反復適用するとき、テンプレート:Mvar の負および正の周波数成分は、それぞれ テンプレート:Math および テンプレート:Math ずれて、これらのずれは一致する。よって、信号自体は符号が反転する: テンプレート:Math。これは (σH(ω))2=e±iπ=1for ω0 による。

ヒルベルト変換表

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

信号u(t) ヒルベルト変換[fn 1]H(u)(t)
正弦[fn 2]sin(ωt) sgn(ω)sin(ωtπ2)=sgn(ω)cos(ωt)
余弦[fn 2]cos(ωt) sgn(ω)cos(ωtπ2)=sgn(ω)sin(ωt)
eiωt sgn(ω)ei(ωtπ2)=isgn(ω)eiωt
1t2+1 tt2+1
et2 2π1/2F(t)(テンプレート:Illを参照)
sinc関数 sin(t)t 1cos(t)t
矩形関数 (t) 1πln|t+12t12|
ディラックのデルタ関数 δ(t) 1πt
特性関数 χ[a,b](t) 1πln|tatb|
Note
  1. 文献によっては(例えば テンプレート:Harvnbテンプレート:Full)、本項で言うと テンプレート:Math にあたるものをヒルベルト順変換の定義とするものもある。その場合には、一覧表の右列はすべて符号が逆になる。
  2. 2.0 2.1 テンプレート:Mathテンプレート:Math のヒルベルト変換は無限遠点において積分の主値をとることで定義することができる(シュヴァルツ超函数の意味で定義したものとも一致する)。

幅広いヒルベルト変換の一覧表が利用可能 (テンプレート:Harvnb). 定数のヒルベルト変換は テンプレート:Math となることに注意。

テンプレート:脚注ヘルプ

注釈

テンプレート:Notelist2

出典

テンプレート:Reflist

参考文献

関連項目

外部リンク

テンプレート:Commonscat