意味 u のヒルベルト変換は、u ( t ) と関数h ( t ) = 1 / π t (コーシー核 として知られる) の畳み込み と考えることができます。1 / tは t = 0 で積分 できないため、畳み込みを定義する積分は必ずしも収束しません。代わりに、ヒルベルト変換はコーシー主値(ここでは pv と表記)を使用して定義されます。具体的には、関数(または信号) u ( t ) のヒルベルト変換は次のように与えられます。
H ( u ) ( t ) = 1 π p 。 v 。 ∫ − ∞ + ∞ u ( τ ) t − τ d τ 、 {\displaystyle \operatorname {H} (u)(t)={\frac {1}{\pi }}\,\operatorname {pv} \int _{-\infty }^{+\infty }{\frac {u(\tau )}{t-\tau }}\,\mathrm {d} \tau ,}
この積分が主値として存在する場合、これはまさにuと 緩増加分布 p.v. 1 / π t の畳み込みである。[ 1 ] あるいは、変数変換によって、主値積分は明示的に次のように書くことができる。[ 2 ]
H ( u ) ( t ) = 2 π リム ε → 0 ∫ ε ∞ u ( t − τ ) − u ( t + τ ) 2 τ d τ 。 {\displaystyle \operatorname {H} (u)(t)={\frac {2}{\pi }}\,\lim _{\varepsilon \to 0}\int _{\varepsilon }^{\infty }{\frac {u(t-\tau )-u(t+\tau )}{2\tau }}\,\mathrm {d} \tau .}
関数u にヒルベルト変換を 2 回連続して適用すると、結果は次のようになります。
H ( H ( u ) ) ( t ) = − u ( t ) 、 {\displaystyle \operatorname {H} {\bigl (}\operatorname {H} (u){\bigr )}(t)=-u(t),}
ただし、両方の反復を定義する積分が適切な意味で収束するものとします。特に、逆変換は − H {\displaystyle -\operatorname {H} } この事実は、 u ( t ) のフーリエ変換 に対するヒルベルト変換の影響を考察することで最も容易に理解できます(下記の§ フーリエ変換との関係を 参照)。
上半平面 における解析関数 に対して、ヒルベルト変換は境界値の実部と虚部の関係を表します。つまり、f ( z ) が複素平面{ z : Im{ z } > 0} で解析関数であり、u ( t ) = Re{ f ( t + 0· i )} である場合、このヒルベルト変換が存在する限り、 Im{ f ( t + 0· i )} = H( u )( t ) は加法定数を除いて成り立ちます。
表記法 信号処理 において、u ( t ) のヒルベルト変換は一般的に次のように表される。u ^ ( t ) \displaystyle {\hat {u}}(t)} [ 3 ] しかし、数学では、この表記法はu ( t ) のフーリエ変換を表すために既に広く用いられている。[ 4 ] ヒルベルト変換は、場合によっては次のように表されることもある。u ~ ( t ) \displaystyle {\tilde {u}}(t)} さらに、多くの文献では、ヒルベルト変換をここで定義したものの負として定義している。[ 5 ]
ヒルベルト変換は乗算演算子で ある。[ 14 ] H の乗算演算子はσ H ( ω ) = − i sgn( ω ) であり、sgnは 符号関数 である。したがって、次のようになる。
F ( H ( u ) ) ( ω ) = − 私 サイン ( ω ) ⋅ F ( u ) ( ω ) 、 {\displaystyle {\mathcal {F}}{\bigl (}\operatorname {H} (u){\bigr )}(\omega )=-i\operatorname {sgn}(\omega )\cdot {\mathcal {F}}(u)(\omega ),}
どこF \displaystyle {\mathcal {F}}} はフーリエ変換 を表す。sgn ( x ) = sgn(2 π x ) であることから、この結果は の 3 つの一般的な定義に適用されることがわかる。F \displaystyle {\mathcal {F}}} 。
オイラーの公式 によれば、 σ H ( ω ) = { 私 = e + 私 π / 2 もし ω < 0 0 もし ω = 0 − 私 = e − 私 π / 2 もし ω > 0 {\displaystyle \sigma _{\operatorname {H} }(\omega )={\begin{cases}~~i=e^{+i\pi /2}&{\text{if }}\omega <0\\~~0&{\text{if }}\omega =0\\-i=e^{-i\pi /2}&{\text{if }}\omega >0\end{cases}}}
したがって、H( u )( t ) は u ( t ) の負の周波数 成分の位相を+90° ( π ⁄ 2 ラジアン) だけシフトし、正の周波数成分の位相を −90° だけシフトする効果があり、i ·H( u )( t ) は正の周波数成分を復元しつつ、負の周波数成分をさらに +90° だけシフトして、それらの負の値 (つまり −1 倍) をもたらす効果があります。
ヒルベルト変換を2回適用すると、 u ( t ) の負の周波数成分と正の周波数成分の位相はそれぞれ+180°と-180°だけシフトされ、これは等量です。信号は反転します。つまり、H(H( u )) = -u となります。
( σ H ( ω ) ) 2 = e ± 私 π = − 1 のために ω ≠ 0. {\displaystyle \left(\sigma _{\operatorname {H} }(\omega )\right)^{2}=e^{\pm i\pi }=-1\quad {\text{for }}\omega \neq 0.}
以下の表では、周波数 パラメータω {\displaystyle \omega } 本当です。
注記
↑ 一部の著者(例えば、Bracewell)は、順変換の定義として我々の−H を使用しています。その結果、この表の右列が反転されます。 1 2 sin 関数と cos 関数のヒルベルト変換は、積分の主値を無限大で取ることによって定義できます。この定義は、ヒルベルト変換を分布的に定義した結果と一致します。 ヒルベルト変換の詳細な表が利用可能です。 定数のヒルベルト変換はゼロであることに注意してください。
定義領域 ヒルベルト変換がそもそも適切に 定義されていることは決して自明ではない。なぜなら、それを定義する広義積分は適切な意味で収束しなければならないからである。しかし、ヒルベルト変換は、以下の関数の広いクラスに対して適切に定義されている。L p ( R ) {\displaystyle L^{p}(\mathbb {R} )} 1 < p < ∞ の場合。
より正確には、u がL p ( R ) {\displaystyle L^{p}(\mathbb {R} )} 1 < p < ∞ の場合、広義積分を定義する極限は次のようになる。
H ( u ) ( t ) = 2 π リム ε → 0 ∫ ε ∞ u ( t − τ ) − u ( t + τ ) 2 τ d τ {\displaystyle \operatorname {H} (u)(t)={\frac {2}{\pi }}\lim _{\varepsilon \to 0}\int _{\varepsilon }^{\infty }{\frac {u(t-\tau )-u(t+\tau )}{2\tau }}\,d\tau }
ほぼすべての t に対して存在します。極限関数もL p ( R ) {\displaystyle L^{p}(\mathbb {R} )} そして、それは実際には広義積分の平均値の極限でもある。つまり、
2 π ∫ ε ∞ u ( t − τ ) − u ( t + τ ) 2 τ d τ → H ( u ) ( t ) {\displaystyle {\frac {2}{\pi }}\int _{\varepsilon }^{\infty }{\frac {u(t-\tau )-u(t+\tau )}{2\tau }}\,\mathrm {d} \tau \to \operatorname {H} (u)(t)}
ε → 0 のとき、 L p ノルムにおいて、またティッチマーシュの定理 によりほとんど至るところで点ごとに、 となる。
p = 1 の場合、ヒルベルト変換はほとんど至るところで点ごとに収束しますが、局所的にも積分可能にならない場合があります。特に、この場合、平均収束は一般には起こりません。しかし、L 1関数のヒルベルト変換は L 1 弱収束し、ヒルベルト変換はL 1から L 1,w への有界作用素です。 (特に、ヒルベルト変換はL 2 上の乗算作用素でもあるため、マルチンキェヴィッチ補間と双対性議論により、 Hが L p 上で有界であることの別の証明が得られます。)
物件
境界性 1 < p < ∞ の場合、ヒルベルト変換はL p ( R ) {\displaystyle L^{p}(\mathbb {R} )} は有界線形作用素 であり、すなわち、ある定数C pが存在して、
‖ H u ‖ p ≤ C p ‖ u ‖ p {\displaystyle \left\|\operatorname {H} u\right\|_{p}\leq C_{p}\left\|u\right\|_{p}}
すべての人々のためにu ∈ L p ( R ) {\displaystyle u\in L^{p}(\mathbb {R} )} [ 19 ]
最良の定数C p {\displaystyle C_{p}} [ 20 ] で与えられるC p = { タン π 2 p もし 1 < p ≤ 2 ベビーベッド π 2 p もし 2 < p < ∞ {\displaystyle C_{p}={\begin{cases}\tan {\frac {\pi }{2p}}&{\text{if}}~1<p\leq 2\\[4pt]\cot {\frac {\pi }{2p}}&{\text{if}}~2<p<\infty \end{cases}}}
最高のものを見つける簡単な方法C p {\displaystyle C_{p}} のためにp {\displaystyle p} 2のべき乗であることは、いわゆるコトラーの恒等式によって証明される。( H f ) 2 = f 2 + 2 H ( f H f ) {\displaystyle (\operatorname {H} f)^{2}=f^{2}+2\operatorname {H} (f\operatorname {H} f)} すべての実数値関数f に対して、同じ最適定数が周期ヒルベルト変換にも適用されます。
ヒルベルト変換の有界性は、L p ( R ) {\displaystyle L^{p}(\mathbb {R} )} 対称部分和演算子の収束 S R f = ∫ − R R f ^ ( ξ ) e 2 π 私 x ξ d ξ {\displaystyle S_{R}f=\int _{-R}^{R}{\hat {f}}(\xi )e^{2\pi ix\xi }\,\mathrm {d} \xi }
に L p ( R ) {\displaystyle L^{p}(\mathbb {R} )} [ 21 ]
反自己随伴性 ヒルベルト変換は、双対ペアリングに関して反自己共役演算子である。 L p ( R ) {\displaystyle L^{p}(\mathbb {R} )} そしてデュアルスペースL q ( R ) {\displaystyle L^{q}(\mathbb {R} )} ここで、 p とqは ヘルダー共役で あり、1 < p , q < ∞で ある。記号的には、
⟨ H u 、 v ⟩ = ⟨ u 、 − H v ⟩ {\displaystyle \langle \operatorname {H} u,v\rangle =\langle u,-\operatorname {H} v\rangle }
のためにu ∈ L p ( R ) {\displaystyle u\in L^{p}(\mathbb {R} )} そしてv ∈ L q ( R ) {\displaystyle v\in L^{q}(\mathbb {R} )} [
ヒルベルト変換は反対合である 、つまり
H ( H ( u ) ) = − u {\displaystyle \operatorname {H} {\bigl (}\operatorname {H} \left(u\right){\bigr )}=-u}
各変換が明確に定義されている限り、H は 空間を保存する。L p ( R ) {\displaystyle L^{p}(\mathbb {R} )} これは特に、ヒルベルト変換が可逆であることを意味する。 L p ( R ) {\displaystyle L^{p}(\mathbb {R} )} そして、
H − 1 = − H {\displaystyle \operatorname {H} ^{-1}=-\operatorname {H} }
分化 厳密に言えば、ヒルベルト変換の導関数は導関数のヒルベルト変換であり、つまりこれら2つの線形演算子は可換である。
H ( d u d t ) = d d t H ( u ) {\displaystyle \operatorname {H} \left({\frac {\mathrm {d} u}{\mathrm {d} t}}\right)={\frac {\mathrm {d} }{\mathrm {d} t}}\operatorname {H} (u)}
この同一性を反復すると、
H ( d k u d t k ) = d k d t k H ( u ) {\displaystyle \operatorname {H} \left({\frac {\mathrm {d} ^{k}u}{\mathrm {d} t^{k}}}\right)={\frac {\mathrm {d} ^{k}}{\mathrm {d} t^{k}}}\operatorname {H} (u)}
これは、u とその最初のk 階微分が以下に属する場合に厳密に真である。L p ( R ) {\displaystyle L^{p}(\mathbb {R} )} [ 24 ] 周波数領域では微分がωによる乗算になるので、これは簡単に確認できます 。
畳み込み ヒルベルト変換は、形式的には緩増加分布 との畳み込み として実現できる
h ( t ) = p 。 v 。 1 π t {\displaystyle h(t)=\operatorname {p.v.} {\frac {1}{\pi \,t}}}
したがって正式には、
H ( u ) = h * u {\displaystyle \operatorname {H} (u)=h*u}
しかし、これは先験的に コンパクトサポート の分布u に対してのみ定義できる。コンパクトサポート関数 (当然ながら 分布である) はL p で稠密であるため、これに対してある程度厳密に作業することが可能である。あるいは、 h ( t ) が関数log| t |/ π の分布微分 であるという事実を利用することもできる。すなわち、
H ( u ) ( t ) = d d t ( 1 π ( u * ログ | ⋅ | ) ( t ) ) {\displaystyle \operatorname {H} (u)(t)={\frac {\mathrm {d} }{\mathrm {d} t}}\left({\frac {1}{\pi }}\left(u*\log {\bigl |}\cdot {\bigr |}\right)(t)\right)}
ほとんどの運用上の目的においては、ヒルベルト変換は畳み込みとして扱うことができます。例えば、形式的には、畳み込みのヒルベルト変換は、ヒルベルト変換をいずれかの因子にのみ適用した畳み込みです。
H ( u * v ) = H ( u ) * v = u * H ( v ) {\displaystyle \operatorname {H} (u*v)=\operatorname {H} (u)*v=u*\operatorname {H} (v)}
これは、 u とvが コンパクトサポート分布である場合に厳密に真である。なぜなら、その場合、
h * ( u * v ) = ( h * u ) * v = u * ( h * v ) {\displaystyle h*(u*v)=(h*u)*v=u*(h*v)}
適切な極限を取ることにより、u ∈ L p かつv ∈ L q の場合にも、以下のことが成り立つ。
1 < 1 p + 1 q {\displaystyle 1<{\frac {1}{p}}+{\frac {1}{q}}}
ティッチマーシュによる定理から。
定義の範囲を拡張する
さらに、ヒルベルト変換を特定の分布 空間に拡張することも可能です( Pandey 1996 、第3章) 。ヒルベルト変換は微分と可換であり、 L p 上の有界作用素であるため、Hを制限すると、 ソボレフ空間 の逆極限 で連続変換が得られます。
D L p = リム ⟵ n → ∞ W n 、 p ( R ) {\displaystyle {\mathcal {D}}_{L^{p}}={\underset {n\to \infty }{\underset {\longleftarrow }{\lim }}}W^{n,p}(\mathbb {R} )}
ヒルベルト変換は、双対空間上で定義できる。D L p {\displaystyle {\mathcal {D}}_{L^{p}}} 、と表記されるD L p ′ {\displaystyle {\mathcal {D}}_{L^{p}}'} L p 分布から構成される。これは双対ペアリングによって実現される 。u ∈ D L p ′ {\displaystyle u\in {\mathcal {D}}'_{L^{p}}} 、 定義する:
H ( u ) ∈ D L p ′ = ⟨ H u 、 v ⟩ ≜ ⟨ u 、 − H v ⟩ 、 すべての人々のために v ∈ D L p 。 {\displaystyle \operatorname {H} (u)\in {\mathcal {D}}'_{L^{p}}=\langle \operatorname {H} u,v\rangle \ \triangleq \ \langle u,-\operatorname {H} v\rangle ,\ {\text{for all}}\ v\in {\mathcal {D}}_{L^{p}}.}
ゲルファントとシロフによるアプローチ[ 29 ] によって、緩増加分布 の空間上でヒルベルト変換を定義することも可能だが、積分に特異点があるため、より注意が必要である。
ヒルベルト変換は、以下の関数に対して定義できる。L ∞ ( R ) {\displaystyle L^{\infty }(\mathbb {R} )} 同様だが、いくつかの修正と注意点が必要となる。正しく理解すれば、ヒルベルト変換はL ∞ ( R ) {\displaystyle L^{\infty }(\mathbb {R} )} 有界平均振動 (BMO)クラスのバナッハ空間 へ。
素朴に解釈すると、有界関数のヒルベルト変換は明らかに定義不良である。例えば、u = sgn( x )の場合、 H( u ) を定義する積分はほとんど至るところで±∞に発散する。このような困難を軽減するために、 L ∞ 関数のヒルベルト変換は、積分の次の正則化された形式で定義される。
H ( u ) ( t ) = p 。 v 。 ∫ − ∞ ∞ u ( τ ) { h ( t − τ ) − h 0 ( − τ ) } d τ {\displaystyle \operatorname {H} (u)(t)=\operatorname {p.v.} \int _{-\infty }^{\infty }u(\tau )\left\{h(t-\tau )-h_{0}(-\tau )\right\}\,\mathrm {d} \tau }
ここで、上記のようにh ( x ) = 1 / πx で あり、
h 0 ( x ) = { 0 もし | x | < 1 1 π x もし | x | ≥ 1 {\displaystyle h_{0}(x)={\begin{cases}0&{\text{if}}~|x|<1\\{\frac {1}{\pi \,x}}&{\text{if}}~|x|\geq 1\end{cases}}}
修正された変換H は、Calderón と Zygmund による一般的な結果から、コンパクト サポートの関数については、加法定数を除いて元の変換と一致します。[ 30 ] さらに、結果として得られる積分は、ほとんど至るところで点ごとに収束し、BMO ノルムに関して、有界平均振動の関数になります。
フェファーマンの研究[ 31 ] の重要な成果は 、関数が有界平均振動を持つのは、ある関数がf + H( g )の形である場合に限るということである。 f 、 g ∈ L ∞ ( R ) {\displaystyle f,g\in L^{\infty }(\mathbb {R} )} 。
共役関数 ヒルベルト変換は、関数f ( x ) とg ( x ) のペアで理解することができ、関数は F ( x ) = f ( x ) + 私 g ( x ) {\displaystyle F(x)=f(x)+i\,g(x)} は、上半平面における正則関数 F ( z ) の境界値です。 このような状況下では、f とg が 十分に積分可能であれば、一方が他方のヒルベルト変換となります。
仮にf ∈ L p ( R ) 。 {\displaystyle f\in L^{p}(\mathbb {R} ).} すると、ポアソン積分 の理論により、f は 上半平面への一意の調和拡張を許容し、この拡張は次のように与えられる。
u ( x + 私 y ) = u ( x 、 y ) = 1 π ∫ − ∞ ∞ f ( s ) y ( x − s ) 2 + y 2 d s {\displaystyle u(x+iy)=u(x,y)={\frac {1}{\pi }}\int _{-\infty }^{\infty }f(s)\;{\frac {y}{(x-s)^{2}+y^{2}}}\;\mathrm {d} s}
これは、 f とポアソンカーネル の畳み込みである。
P ( x 、 y ) = y π ( x 2 + y 2 ) {\displaystyle P(x,y)={\frac {y}{\pi \,\left(x^{2}+y^{2}\right)}}}
さらに、上半平面で定義された一意の調和関数v が存在し、 F ( z ) = u ( z ) + iv ( z ) は正則であり、 リム y → ∞ v ( x + 私 y ) = 0 {\displaystyle \lim _{y\to \infty }v\,(x+i\,y)=0}
この調和関数は、 fを 共役ポアソンカーネル と畳み込むことによって得られる。
Q ( x 、 y ) = x π ( x 2 + y 2 ) 。 {\displaystyle Q(x,y)={\frac {x}{\pi \,\left(x^{2}+y^{2}\right)}}.}
したがって v ( x 、 y ) = 1 π ∫ − ∞ ∞ f ( s ) x − s ( x − s ) 2 + y 2 d s 。 {\displaystyle v(x,y)={\frac {1}{\pi }}\int _{-\infty }^{\infty }f(s)\;{\frac {x-s}{\,(x-s)^{2}+y^{2}\,}}\;\mathrm {d} s.}
実際、コーシー核の実部と虚部は 私 π z = P ( x 、 y ) + 私 Q ( x 、 y ) {\displaystyle {\frac {i}{\pi \,z}}=P(x,y)+i\,Q(x,y)}
コーシーの積分公式 により、F = u + iv は 正則である。
このようにしてu から得られる関数v は、 u の調和共役と呼ばれます。 y → 0 のときのv ( x , y ) の (非接線方向の) 境界極限は、f のヒルベルト変換です。 したがって、簡潔に言うと、 H ( f ) = リム y → 0 Q ( − 、 y ) ⋆ f {\displaystyle \operatorname {H} (f)=\lim _{y\to 0}Q(-,y)\star f}
ティッチマーシュの定理ティッチマーシュの定理(1937年の著作にこの定理を含めたECティッチマーシュにちなんで名付けられた)は、上半平面における正則関数の境界値とヒルベルト変換の関係を厳密に定義する。 この定理は、実数直線上の複素数値の二乗可積分 関数F ( x )が、 上半平面U における正則関数のハーディ空間 H2 ( U ) の関数の境界値となるための必要十分条件を与える。
この定理は、複素数値の二乗可積分関数に対して以下の条件が成り立つことを述べている。F : R → C {\displaystyle F:\mathbb {R} \to \mathbb {C} } 同等です。
F ( x ) は、上半平面における正則関数F ( z )の z → x の極限であり∫ − ∞ ∞ | F ( x + 私 y ) | 2 d x < K {\displaystyle \int _{-\infty }^{\infty }|F(x+i\,y)|^{2}\;\mathrm {d} x<K} F ( x ) の実部と虚部は互いにヒルベルト変換の関係にある。フーリエ変換 F ( F ) ( x ) {\displaystyle {\mathcal {F}}(F)(x)} x < 0 の場合、消滅する。 p > 1 の場合、クラスL p の関数に対しては、より弱い結果が成り立つ。具体的には、 F ( z ) が正則関数である場合、
∫ − ∞ ∞ | F ( x + 私 y ) | p d x < K {\displaystyle \int _{-\infty }^{\infty }|F(x+i\,y)|^{p}\;\mathrm {d} x<K}
すべてのy に対して、複素数値関数F ( x ) が存在する。L p ( R ) {\displaystyle L^{p}(\mathbb {R} )} y → 0のとき、 L p ノルムでF ( x + iy ) → F ( x ) となる(また、ほとんど至るところで 点ごとに成り立つ)。さらに、
F ( x ) = f ( x ) + 私 g ( x ) {\displaystyle F(x)=f(x)+i\,g(x)}
ここで、f は実数値関数であり、L p ( R ) {\displaystyle L^{p}(\mathbb {R} )} gは f のヒルベルト変換 (クラスL p )です。
これはp = 1 の 場合には当てはまりません。実際、L 1 関数f のヒルベルト変換は、平均的には別のL 1 関数に収束するとは限りません。しかしながら、 f のヒルベルト変換は、ほとんど至るところで有限関数g に収束し、
∫ − ∞ ∞ | g ( x ) | p 1 + x 2 d x < ∞ {\displaystyle \int _{-\infty }^{\infty }{\frac {|g(x)|^{p}}{1+x^{2}}}\;\mathrm {d} x<\infty }
この結果は、アンドレイ・コルモゴロフ による円盤内のハーディ関数に関する結果と直接類似している。 通常はティッチマーシュの定理と呼ばれているが、この結果はハーディ、ペイリー、ウィーナー(ペイリー・ウィーナーの定理を 参照)を含む他の多くの研究、およびリース、ヒル、タマルキンの研究を統合したものである。[ 37 ]
リーマン・ヒルベルト問題リーマン・ヒルベルト問題 の1つの形式は、上 半平面で正則で あり、下半平面で正則であるような関数 F + と F − のペアを特定しようと する もので 、 実 軸に沿った x に対して、 F + ( x ) − F − ( x ) = f ( x ) {\displaystyle F_{+}(x)-F_{-}(x)=f(x)}
ここで、f ( x ) は、ある与えられた実数値関数である。x ∈ R {\displaystyle x\in \mathbb {R} } この式の左辺は、F ± の 適切な半平面からの極限の差、または超関数 分布として理解できます。この形式の2つの関数は、リーマン・ヒルベルト問題の解です。
形式的には、F ± が リーマン・ヒルベルト問題を解く 場合f ( x ) = F + ( x ) − F − ( x ) {\displaystyle f(x)=F_{+}(x)-F_{-}(x)}
すると、 f ( x ) のヒルベルト変換は[ 38 ] で与えられる。H ( f ) ( x ) = − 私 ( F + ( x ) + F − ( x ) ) 。 {\displaystyle H(f)(x)=-i{\bigl (}F_{+}(x)+F_{-}(x){\bigr )}.}
周期関数 f に対して、円形ヒルベルト変換は次のように定義される。
f ~ ( x ) ≜ 1 2 π p 。 v 。 ∫ 0 2 π f ( t ) ベビーベッド ( x − t 2 ) d t {\displaystyle {\tilde {f}}(x)\triangleq {\frac {1}{2\pi }}\operatorname {p.v.} \int _{0}^{2\pi }f(t)\,\cot \left({\frac {x-t}{2}}\right)\,\mathrm {d} t}
円形ヒルベルト変換は、ハーディ空間の特徴付けやフーリエ級数 における共役関数の研究に用いられる。核は、 ベビーベッド ( x − t 2 ) {\displaystyle \cot \left({\frac {x-t}{2}}\right)} これは、ヒルベルト変換が最初に研究された形式であるため、ヒルベルト核 として知られています。
ヒルベルト核(円形ヒルベルト変換の場合)は、コーシー核1 ⁄ x を 周期的にすることで得られる。より正確には、x ≠ 0の場合、
1 2 ベビーベッド ( x 2 ) = 1 x + ∑ n = 1 ∞ ( 1 x + 2 n π + 1 x − 2 n π ) {\displaystyle {\frac {1}{\,2\,}}\cot \left({\frac {x}{2}}\right)={\frac {1}{x}}+\sum _{n=1}^{\infty }\left({\frac {1}{x+2n\pi }}+{\frac {1}{\,x-2n\pi \,}}\right)}
この対応関係から得られるヒルベルト変換に関する結果から、円形ヒルベルト変換に関する多くの結果を導き出すことができる。
もう一つのより直接的なつながりは、ケイリー変換 C ( x ) = ( x – i ) / ( x + i ) によって提供される。これは実数直線を円に、上半平面を単位円盤 に写像する。これはユニタリ写像を誘導する。
U f ( x ) = 1 ( x + 私 ) π f ( C ( x ) ) {\displaystyle U\,f(x)={\frac {1}{(x+i)\,{\sqrt {\pi }}}}\,f\left(C\left(x\right)\right)}
L 2 ( T ) を L 2 ( R ) 。 {\displaystyle L^{2}(\mathbb {R} ).} 演算子U はハーディ空間H 2 ( T ) をハーディ空間に運びますH 2 ( R ) {\displaystyle H^{2}(\mathbb {R} )} [
ベドロシアンの定理ベドロシアンの定理 によれば、スペクトルが重ならないローパス信号とハイパス信号の積のヒルベルト変換は、ローパス信号とハイパス信号のヒルベルト変換の積で与えられる。
H ( f LP ( t ) ⋅ f HP ( t ) ) = f LP ( t ) ⋅ H ( f HP ( t ) ) 、 {\displaystyle \operatorname {H} \left(f_{\text{LP}}(t)\cdot f_{\text{HP}}(t)\right)=f_{\text{LP}}(t)\cdot \operatorname {H} \left(f_{\text{HP}}(t)\right),}
ここで、f LP とf HP は それぞれ低域通過信号と高域通過信号である。 このことが適用される通信信号のカテゴリは、狭帯域信号モデルと呼ばれる。 そのカテゴリの1つは、高周波正弦波「搬送波」の振幅変調である。
u ( t ) = u m ( t ) ⋅ コス ( ω t + φ ) 、 {\displaystyle u(t)=u_{m}(t)\cdot \cos(\omega t+\varphi ),}
ここで、u m ( t ) は、音声や音楽などの狭帯域の「メッセージ」波形です。ベドロシアンの定理により、次のようになります。
H ( u ) ( t ) = { + u m ( t ) ⋅ 罪 ( ω t + φ ) もし ω > 0 − u m ( t ) ⋅ 罪 ( ω t + φ ) もし ω < 0 {\displaystyle \operatorname {H} (u)(t)={\begin{cases}+u_{m}(t)\cdot \sin(\omega t+\varphi )&{\text{if }}\omega >0\\-u_{m}(t)\cdot \sin(\omega t+\varphi )&{\text{if }}\omega <0\end{cases}}}
解析的表現 共役関数 の特定の種類は次のとおりです。
u 1 ( t ) ≜ u ( t ) + 私 ⋅ H ( u ) ( t ) 、 {\displaystyle u_{a}(t)\triangleq u(t)+i\cdot H(u)(t),}
解析的表現 として知られるu ( t ) 。 {\displaystyle u(t).} その名前は、主にオイラーの公式 による数学的な扱いやすさを反映している。狭帯域モデルにベドロシアンの定理を適用すると、解析的表現は次のようになる。 [ 42 ]
フーリエ変換の性質から、この複素ヘテロダイン演算によって、 u m ( t ) の負の周波数成分すべてを0 Hz より上にシフトできることが分かります。この場合、結果の虚部は実部のヒルベルト変換となります。これは、ヒルベルト変換を生成する間接的な方法です。
角度(位相/周波数)変調形式: [ 43 ]
u ( t ) = A ⋅ コス ( ω t + φ m ( t ) ) {\displaystyle u(t)=A\cdot \cos(\omega t+\varphi _{m}(t))}
は角度変調 と呼ばれ、位相変調 と周波数変調の 両方を含みます。瞬時周波数 は ω + φ m ′ ( t ) 。 {\displaystyle \omega +\varphi _{m}^{\prime }(t).} ω が十分に大きい場合、φ m ′ {\displaystyle \varphi _{m}^{\prime }} :
H ( u ) ( t ) ≈ A ⋅ 罪 ( ω t + φ m ( t ) ) {\displaystyle \operatorname {H} (u)(t)\approx A\cdot \sin(\omega t+\varphi _{m}(t))} そして: u 1 ( t ) ≈ A ⋅ e 私 ( ω t + φ m ( t ) ) 。 {\displaystyle u_{a}(t)\approx A\cdot e^{i(\omega t+\varphi _{m}(t))}.}
単側波帯変調(SSB)式1 のu m ( t ) が(メッセージ波形の)解析的表現でもある場合、すなわち:
u m ( t ) = m ( t ) + 私 ⋅ m ^ ( t ) {\displaystyle u_{m}(t)=m(t)+i\cdot {\widehat {m}}(t)}
その結果、単側波帯 変調となる。
u 1 ( t ) = ( m ( t ) + 私 ⋅ m ^ ( t ) ) ⋅ e 私 ( ω t + φ ) {\displaystyle u_{a}(t)=(m(t)+i\cdot {\widehat {m}}(t))\cdot e^{i(\omega t+\varphi )}}
送信成分は次のとおりです。[ 44 ] [ 45 ]
u ( t ) = 再 { u 1 ( t ) } = m ( t ) ⋅ コス ( ω t + φ ) − m ^ ( t ) ⋅ 罪 ( ω t + φ ) {\displaystyle {\begin{aligned}u(t)&=\operatorname {Re} \{u_{a}(t)\}\\&=m(t)\cdot \cos(\omega t+\varphi )-{\widehat {m}}(t)\cdot \sin(\omega t+\varphi )\end{aligned}}}
因果関係 機能h ( t ) = 1 / ( π t ) {\displaystyle h(t)=1/(\pi t)} 畳み込み演算の実践的な実装において、因果関係に基づく2つの課題が提示される(0における値が未定義であることに加えて)。
図 1 :周波数応答がナイキスト周波数の約95%に帯域制限されたフィルタ図 2 :ハイパス周波数応答を持つヒルベルト変換フィルタ図 3 。図 4。cos ( ωt ) のヒルベルト変換はsin( ωt ) です。この図は、sin(ωt) と、MATLABライブラリ関数hilbert()によって計算された2つの近似ヒルベルト変換を示しています。 図 5. 区分的畳み込みを用いたコサイン関数の離散ヒルベルト変換離散関数の場合、u [ n ] 、 {\displaystyle u[n],} 離散時間フーリエ変換 (DTFT)を用いて、U ( ω ) {\displaystyle U(\omega )} 離散ヒルベルト変換 u ^ [ n ] 、 {\displaystyle {\widehat {u}}[n],} DTFTu ^ [ n ] {\displaystyle {\widehat {u}}[n]} − π < ω < π の領域では、次のように与えられます。
DTFT ( u ^ ) = U ( ω ) ⋅ ( − 私 ⋅ サイン ( ω ) ) 。 {\displaystyle \operatorname {DTFT} ({\widehat {u}})=U(\omega )\cdot (-i\cdot \operatorname {sgn}(\omega )).} 畳み込み定理 を用いた逆DTFTは次のようになります。 [ 46 ] [ 47 ]
u ^ [ n ] = D T F T − 1 ( U ( ω ) ) * D T F T − 1 ( − 私 ⋅ サイン ( ω ) ) = u [ n ] * 1 2 π ∫ − π π ( − 私 ⋅ サイン ( ω ) ) ⋅ e 私 ω n d ω = u [ n ] * 1 2 π [ ∫ − π 0 私 ⋅ e 私 ω n d ω − ∫ 0 π 私 ⋅ e 私 ω n d ω ] ⏟ h [ n ] 、 {\displaystyle {\begin{aligned}{\widehat {u}}[n]&={\scriptstyle \mathrm {DTFT} ^{-1}}(U(\omega ))\ *\ {\scriptstyle \mathrm {DTFT} ^{-1}}(-i\cdot \operatorname {sgn}(\omega ))\\&=u[n]\ *\ {\frac {1}{2\pi }}\int _{-\pi }^{\pi }(-i\cdot \operatorname {sgn}(\omega ))\cdot e^{i\omega n}\,\mathrm {d} \omega \\&=u[n]\ *\ \underbrace {{\frac {1}{2\pi }}\left[\int _{-\pi }^{0}i\cdot e^{i\omega n}\,\mathrm {d} \omega -\int _{0}^{\pi }i\cdot e^{i\omega n}\,\mathrm {d} \omega \right]} _{h[n]},\end{aligned}}} どこ
h [ n ] ≜ { 0 、 もし n 平 2 π n もし n 奇数 {\displaystyle h[n]\ \triangleq \ {\begin{cases}0,&{\text{if }}n{\text{ even}}\\{\frac {2}{\pi n}}&{\text{if }}n{\text{ odd}}\end{cases}}} これは無限インパルス応答(IIR)です。
実務上の考慮事項 [ 48 ]
方法1: ストリーミングの直接畳み込みu [ n ] {\displaystyle u[n]} FIR近似を用いたデータh [ n ] 、 {\displaystyle h[n],} 我々はこれを次のように指定するh ~ [ n ] 。 {\displaystyle {\tilde {h}}[n].} 切り捨ての例h [ n ] {\displaystyle h[n]} 図1と図2に示されています。 図1は 反対称係数の数が奇数で、タイプ IIIと呼ばれます。[ 49 ] このタイプは、周波数0とナイキストで振幅がゼロの応答を本質的に示し、結果としてバンドパスフィルタの形状になります。[ 50 ] [ 51 ] タイプ IV設計(反対称係数の数が偶数)を図2 に示します。[ 52 ] [ 53 ] これはハイパス周波数応答を持ちます。[ 54 ] タイプIIIが一般的な選択肢です。[ 55 ] [ 56 ] 理由は次のとおりです。
典型的な(つまり、適切にフィルタリングされ、サンプリングされた)u [ n ] {\displaystyle u[n]} このシーケンスには、ナイキスト周波数において有用な成分は存在しない。 タイプ IVのインパルス応答には1 2 {\displaystyle {\tfrac {1}{2}}} サンプルシフトh [ n ] {\displaystyle h[n]} これにより、図 2 に示すように、ゼロ値の係数がゼロ以外の値になります。したがって、タイプ III設計はタイプIV設計の2倍の効率を持つ可能性があります 。 タイプ III 設計の群遅延 はサンプルの整数値であり、これによりアライメントが容易になります。u ^ [ n ] {\displaystyle {\widehat {u}}[n]} とu [ n ] {\displaystyle u[n]} 解析信号を 生成するため。タイプIVの群遅延は 、2つのサンプルの中間値です。 突然の中断h [ n ] {\displaystyle h[n]} フラットな周波数応答に波紋(ギブス効果)が生じます。これは、ウィンドウ関数を使用してテーパーをかけることで軽減できます。h ~ [ n ] {\displaystyle {\tilde {h}}[n]} ゼロに。[ 57 ]
方法 2: 区分的畳み込み。直接畳み込みは、畳み込み定理を介して高速フーリエ変換の効率性にアクセスできるオーバーラップセーブ などの方法よりも計算負荷がはるかに高いことはよく知られています。 [ 58 ] 具体的には、セグメントの離散フーリエ変換(DFT) は、 u [ n ] {\displaystyle u[n]} は、DFT を各ポイントごとに乗算します。h ~ [ n ] {\displaystyle {\tilde {h}}[n]} シーケンス。積に対して逆離散フーリエ変換(DFT)を実行し、セグメントの先頭と末尾の一時的なアーティファクトを破棄します。入力セグメントが重なり合うことで、出力ストリームにギャップが生じません。同等の時間領域記述は、長さのセグメントです。N {\displaystyle N} (任意のパラメータ)は周期関数と畳み込まれる。
h ~ N [ n ] ≜ ∑ m = − ∞ ∞ h ~ [ n − m N ] 。 {\displaystyle {\tilde {h}}_{N}[n]\ \triangleq \sum _{m=-\infty }^{\infty }{\tilde {h}}[n-mN].} 非ゼロ値の継続時間がh ~ [ n ] {\displaystyle {\tilde {h}}[n]} はM < N 、 {\displaystyle M<N,} 出力シーケンスには以下が含まれるN − M + 1 {\displaystyle N-M+1} サンプルu ^ 。 {\displaystyle {\widehat {u}}.} M − 1 {\displaystyle M-1} 各ブロックからの出力は破棄されますN 、 {\displaystyle N,} そして、入力ブロックはその分だけ重なり合うことで、隙間が生じないようにする。
方法3: 方法2と同じだが、DFTはh ~ [ n ] {\displaystyle {\tilde {h}}[n]} はサンプルに置き換えられます− 私 サイン ( ω ) {\displaystyle -i\operatorname {sgn} (\omega )} 分布(実数部と虚数部はすべて単に0 {\displaystyle 0} または ± 1. {\displaystyle \pm 1.} ) それはu [ n ] {\displaystyle u[n]} 周期的な総和 を伴う: [ A ]
h N [ n ] ≜ ∑ m = − ∞ ∞ h [ n − m N ] 、 {\displaystyle h_{N}[n]\ \triangleq \sum _{m=-\infty }^{\infty }h[n-mN],} [ B ] [ C ] 任意のパラメータに対して、N 。 {\displaystyle N.} h [ n ] {\displaystyle h[n]} これはFIRフィルターではないため、エッジ効果は変換全体に及びます。削除する対象とそれに伴うオーバーラップの量は、アプリケーションに依存する設計上の問題です。
図3は、 方法2と方法3の違いを示しています。反対称インパルス応答の半分のみが表示され、非ゼロ係数のみが表示されています。青色のグラフは方法2に対応しており、h [ n ] {\displaystyle h[n]} テーパーではなく、矩形窓関数によって切り捨てられます。これは、Matlab 関数hilb(65) によって生成されます。その過渡的な影響は正確にわかっており、容易に無視できます。関数引数によって決定される周波数応答は、アプリケーションに依存する唯一の設計上の問題です。
赤いグラフはh 512 [ n ] 、 {\displaystyle h_{512}[n],} 方法3に対応する。これは、の逆DFTである。− 私 サイン ( ω ) {\displaystyle -i\operatorname {sgn} (\omega )} 分布。具体的には、それは、あるセグメントと畳み込まれた関数です。u [ n ] {\displaystyle u[n]} MATLAB 関数hilbert(u,512) によって。[ 61 ] 出力シーケンスの実部は元の入力シーケンスなので、複素出力は解析的表現 である。u [ n ] 。 {\displaystyle u[n].}
入力が純粋なコサインのセグメントである場合、2つの異なる値に対する結果の畳み込みはN {\displaystyle N} 図4 (赤と青のプロット)に示されている。エッジ効果により、結果は純粋な正弦関数(緑のプロット)にはならない。h N [ n ] {\displaystyle h_{N}[n]} FIRシーケンスではないため、理論上の効果範囲は出力シーケンス全体です。ただし、正弦関数との差はエッジからの距離とともに減少します。パラメータN {\displaystyle N} は出力シーケンスの長さです。入力シーケンスの長さを超えると、入力はゼロ値の要素を追加することによって変更されます。ほとんどの場合、これによりエッジ歪みの大きさが軽減されます。しかし、その持続時間は、固有の立ち上がり時間と立ち下がり時間によって決まります。h [ n ] {\displaystyle h[n]} 衝動反応。
図5は 、方法2(青)と方法3(赤点)の両方を使用した区分的畳み込みの例です。正弦関数は、4つの重なり合うセグメントで処理され、再び結合されたコサイン関数の離散ヒルベルト変換を計算することによって作成されます。FIRの結果(青)が示すように、IIRの結果(赤)に現れる歪みは、h [ n ] {\displaystyle h[n]} そしてh N [ n ] {\displaystyle h_{N}[n]} (図3 の緑と赤)。h N [ n ] {\displaystyle h_{N}[n]} テーパー状(窓状 )になっていることは、この文脈では実際には役立ちます。本当の問題は、窓状になっていないことです。実際には、M = N 、 {\displaystyle M=N,} 一方、オーバーラップセーブ法ではM < N 。 {\displaystyle M<N.}
数論的ヒルベルト変換は、離散ヒルベルト変換を適切な素数を法とする整数に拡張したものです 。この点において、離散フーリエ変換の 数論的変換への一般化に従います。数論的ヒルベルト変換は、直交離散シーケンスの集合を生成するために使用できます
参考文献 Bargmann, V. (1947). 「ローレンツ群の既約ユニタリ表現」. Ann. of Math . 48 (3): 568–640 . doi : 10.2307/1969129 . JSTOR 1969129 . Bedrosian, E. (1962年12月).ヒルベルト変換の積定理(PDF) (報告書). Rand Corporation. RM-3439-PR. Bitsadze, AV (2001) [1994]、「解析関数論の境界値問題」、数学百科事典 、EMS Press ブレイスウェル、R. (2000).フーリエ変換とその応用 (第3 版). マグロウヒル. ISBN 0-07-116043-4 。ブランドウッド、デイビッド(2003)。レーダーと信号処理におけるフーリエ変換 。ボストン:アーテックハウス。ISBN 9781580531740 。 Calderón, AP ; Zygmund, A. (1952). 「ある種の特異積分の存在について」 Acta Mathematica . 88 (1): 85– 139. doi : 10.1007/BF02392130 .Carrick, Matt; Jaeger, Doug; harris, fred (2011).デジタル受信機におけるヒルベルト変換器の設計と応用 (PDF) . Chantilly, VA: Proceedings of the SDR 11 Technical Conference and Product Exposition, Wireless Innovation Forum . 2024年6月5日 取得 . Duoandikoetxea, J. (2000).フーリエ解析 . アメリカ数学会. ISBN 0-8218-2172-5 。 Duistermaat, JJ ;コルク、JAC (2010)。ディストリビューション 。ビルクホイザー。土井 : 10.1007/978-0-8176-4675-2。ISBN 978-0-8176-4672-1 。Duren, P. (1970). Theory of H^p Spaces . New York, NY: Academic Press. Fefferman, C. (1971). "有界平均振動の特性" . Bulletin of the American Mathematical Society . 77 (4): 587– 588. doi : 10.1090/S0002-9904-1971-12763-5 . MR 0280994 . Fefferman, C.; Stein, EM (1972). "H^p 多変数空間" . Acta Mathematica . 129 : 137–193 . doi : 10.1007/BF02392215 . MR 0447953 . Franks , LE (1969年9月)。Thomas Kailath (編)。信号理論 。情報理論。Englewood Cliffs, NJ: Prentice Hall。ISBN 0138100772 。Gel'fand, IM ; Shilov, GE (1968).一般化関数 . 第 2 巻. Academic Press. pp. 153–154 . ISBN 0-12-279502-4 。グラファコス、 ルーカス(2004)。古典的および現代的フーリエ解析 。ピアソン・エデュケーション。pp. 253–257。ISBN 0-13-035399-X 。ハーディ、GH ;リトルウッド、JE ;ポリア、G. (1952)。不等式 。ケンブリッジ、英国:ケンブリッジ大学出版局。ISBN 0-521-35880-9 。デイヴィッド・ヒルベルト (1953) [1912]。Grundzüge einer allgemeinen Theorie der lineen Integralgleichungen [ 線形積分方程式の一般理論のフレームワーク ] (ドイツ語)。ライプツィヒとベルリン、デラウェア州 (1912 年)。ニューヨーク州ニューヨーク州 (1953 年): BG Teubner (1912 年)。チェルシーパブ。社(1953)。ISBN 978-3-322-00681-3 . OCLC 988251080 . 2020年12月18日 取得– archive.org経由。 場所(リンク) Isukapalli, Yogananda. 「線形位相FIRフィルタの種類」(PDF) . p. 18. 2024年6月8日 取得 。 Johansson, Mathias. 「ヒルベルト変換、修士論文」(PDF) 。 2012年2月5日にオリジナル(PDF) からアーカイブされました。 ; また、http://www.fuchs-braun.com/media/d9140c7b3d5004fbffff8007ffffffff0.pdfも参照。2021年5月1日にWayback Machine にアーカイブ済み。カク、サブハッシュ (1970)。 「離散ヒルベルト変換」。手順IEEE 。58 (4): 585–586 .土井 : 10.1109/PROC.1970.7696。Kak, Subhash (2014). "数論的ヒルベルト変換". Circuits, Systems and Signal Processing . 33 (8): 2539–2548 . arXiv : 1308.1688 . doi : 10.1007 /s00034-014-9759-8 . S2CID 21226699 . Khvedelidze, BV (2001) [1994]、「ヒルベルト変換」、数学百科事典 、EMS Press King, Frederick W. (2009a).ヒルベルト変換 . 第 1巻. ケンブリッジ、英国: ケンブリッジ大学出版局. キング、フレデリック・W. (2009b).ヒルベルト変換 . 第 2巻. ケンブリッジ、英国: ケンブリッジ大学出版局. p. 453. ISBN 978-0-521-51720-1 。 クレス、ライナー (1989)。線形積分方程式 。ニューヨーク州ニューヨーク州: Springer-Verlag。 p. 91.ISBN 3-540-50616-0 。 ラング、セルジュ (1985)。SL (2、R {\displaystyle \mathbb {R} } ) 。数学の大学院テキスト。 Vol. 105. ニューヨーク州ニューヨーク州: Springer-Verlag。ISBN 0-387-96198-4 。Oppenheim, Alan V. ; Schafer, Ronald W. ; Buck, John R. (1999). "11.4.1".離散時間信号処理 (第2 版). Upper Saddle River, NJ: Prentice Hall. ISBN 0-13-754920-2 タイプIIIおよびIVのFIRシステムの位相の正確さは、ヒルベルト変換器の近似にそれらを使用する説得力のある動機となる 。オズグッド、ブラッド、「フーリエ変換とその応用 (PDF)」 、スタンフォード大学、 2021年4月30日 取得 Pandey, JN (1996).シュワルツ分布のヒルベルト変換とその応用 . Wiley-Interscience. ISBN 0-471-03373-1 。 Pichorides, S. (1972). 「リース、ジグムント、コルモゴロフの定理における定数の最適値について」 . Studia Mathematica . 44 (2): 165–179 . doi : 10.4064/sm-44-2-165-179 .ラビナー、ローレンス・R.、ゴールド、バーナード(1975)。デジタル信号処理の理論と応用 。ニュージャージー州エングルウッド・クリフス:プレンティス・ホール。ISBN 0-13-914101-4 。 マルセル・リース (1928)。 「シュール・レ・フォンション・コンジュゲ」。Mathematische Zeitschrift (フランス語)。27 (1): 218–244 。土井 : 10.1007/BF01171098。S2CID 123261514。 ローゼンブラム、マービン;ロヴニャク、ジェームズ(1997)。ハーディクラスと作用素論 。ドーバー出版。ISBN 0-486-69536-0 。 シュワルツ、ローラン (1950)。分布理論 。パリ、フランス:ヘルマン。Schreier, P.; Scharf, L. (2010).複素数値データの統計的信号処理:不適切信号と非円形信号の理論 . ケンブリッジ、英国:ケンブリッジ大学出版局。 Smith, JO (2007). "Analytic Signals and Hilbert Transform Filters, in Mathematics of the Discrete Fourier Transform (DFT) with Audio Applications" (第 2 版) . 2021-04-29 に取得 。 ; また、https://www.dsprelated.com/freebooks/mdft/Analytic_Signals_Hilbert_Transform.html も参照してください。スタイン、エリアス (1970)。特異積分と関数の微分可能性 。プリンストン大学出版局。ISBN 0-691-08079-8 。スタイン、エリアス ;ワイス、グイド(1971)。ユークリッド空間におけるフーリエ解析入門 。プリンストン大学出版局。ISBN 0-691-08078-X 。杉浦光夫(1990)『ユニタリ表現と調和解析入門』 ノースホランド数学ライブラリー第 44巻(第2版)エルゼビア社 。ISBN 0444885935 。 ティッチマーシュ、E. (1986) [1948].フーリエ積分理論入門 (第2 版). オックスフォード、英国:クラレンドン・プレス. ISBN 978-0-8284-0324-5 。トレッター、スティーブン A. (1995). RWLucky (編). DSP アルゴリズムを用いた通信システム設計 . ニューヨーク: Springer. ISBN 0306450321 。 ジグムント、アントニ (1988) [1968].三角関数シリーズ (第 2 版). ケンブリッジ、英国: ケンブリッジ大学出版局. ISBN 978-0-521-35885-9 。
さらに読む ベネデット、ジョン・J. (1996).調和解析とその応用 . フロリダ州ボカラトン: CRC Press. ISBN 0849378796 。 Carlson、Crilly、Rutledge (2002)。コミュニケーションシステム (第4 版)。McGraw- Hill。ISBN 0-07-011127-8 。 Gold, B.; Oppenheim, AV; Rader, CM (1969). 「離散ヒルベルト変換の理論と実装」. 1969年ブルックリン工科大学シンポジウム議事録 . ニューヨーク. Grafakos, Loukas (1994). 「離散ヒルベルト変換の二乗総和可能性の初等的証明」. American Mathematical Monthly . 101 (5). Mathematical Association of America: 456–458 . doi : 10.2307/2974910 . JSTOR 2974910 . ティッチマーシュ、E. (1926)。 「級数と積分を含む逆数の公式」。数学的ツァイシュリフト 。25 (1): 321–347 。土井 : 10.1007/BF01283842。S2CID 186237099。
外部リンク ヒルベルト変換の有界性の導出 Mathworldのヒルベルト変換— 変換表が含まれています ワイススタイン、エリック・W. 「ティッチマーシュの定理」 . MathWorld .「GS256 講義3:ヒルベルト変換」(PDF) 。 2012年2月27日にオリジナル(PDF) からアーカイブされました。 ヒルベルト変換の入門レベル。