ハートレー変換

数学の分野におけるハートレー変換(ハートレーへんかん、英語: Hartley transform)は、フーリエ変換と非常に関係の深い、実数値関数を実数値関数へと写す積分変換である。1942年、ラルフ・ハートレーによりフーリエ変換の代替的なものとして提唱され、多くの知られているフーリエ関連変換(英語版)の内の一つとなった。フーリエ変換と比較して、ハートレー変換には実関数を実関数へと変換し、逆変換がそれ自身となるという長所がある。

1983年、ロナルド・ブレイスウェル(英語版)によりこの変換の離散版である離散ハートレー変換(英語版)が考案された。

二次元のハートレー変換は、光学フーリエ変換(英語版)と同様なあるアナログ光学処理によって計算される。その利点として、複素フェーズよりも振幅と符号のみが必要とされる、ということが提唱されている[1]。しかし、光学ハートレー変換は未だ広く利用されてはいないようである。

定義

実数値関数 f(t) のハートレー変換は

H ( ω ) = { H f } ( ω ) = 1 2 π f ( t ) cas ( ω t ) d t {\displaystyle H(\omega )=\left\{{\mathcal {H}}f\right\}(\omega )={\frac {1}{\sqrt {2\pi }}}\int _{-\infty }^{\infty }f(t)\,{\mbox{cas}}(\omega t)\mathrm {d} t}

により定義される。ここで、応用の場面での ω {\displaystyle \omega } の意味は角周波数であり、

cas ( t ) = cos ( t ) + sin ( t ) = 2 sin ( t + π / 4 ) = 2 cos ( t π / 4 ) {\displaystyle {\mbox{cas}}(t)=\cos(t)+\sin(t)={\sqrt {2}}\sin(t+\pi /4)={\sqrt {2}}\cos(t-\pi /4)\,}

は「余弦正弦(cosine-and-sine)」あるいは「ハートレー核」と呼ばれるものである。工学において、この変換は信号(関数)を時間領域からハートレースペクトル領域(周波数領域)へと写す。

逆変換

ハートレー変換は、それ自身が逆変換(対合)であるという便利な性質を持つ。すなわち

f = { H { H f } } {\displaystyle f=\{{\mathcal {H}}\{{\mathcal {H}}f\}\}}

が成立する。

慣例

上述の定義はハートレーによる元々の定義と一致するものである。しかし、様々な些細な点が慣例上の問題となっており、それらは本質的な性質を脅かすことなく変更することも出来る:

  • 元の変換と逆変換を等しいものとする代わりに、元の変換から係数 1 / 2 π {\displaystyle {1}/{\sqrt {2\pi }}} を除き、その代わりに逆変換に係数 1 / 2 π {\displaystyle {1}/{2\pi }} を付けることも出来る。あるいは、実際、積が 1 / 2 π {\displaystyle {1}/{2\pi }} となるような組み合わせならどのような係数を用いて正規化してもよい(上の例のような非対称的な正規化は、純粋数学および工学のいずれの分野においても用いられる)。
  • ω t {\displaystyle \omega t} の代わりに 2 π ν t {\displaystyle 2\pi \nu t} を用いることも出来る(すなわち、角周波数ではなく周波数を用いることも出来る)。そのような場合、係数 1 / 2 π {\displaystyle {1}/{\sqrt {2\pi }}} は完全に取り除かれる。
  • 積分核として、cos+sin の代わりに cos−sin を用いることも出来る。

フーリエ変換との関係

この変換は、積分核の選び方において従来のフーリエ変換 F ( ω ) = F { f ( t ) } ( ω ) {\displaystyle F(\omega )={\mathcal {F}}\{f(t)\}(\omega )} とは異なるものである。フーリエ変換においては、指数関数核: exp ( i ω t ) = cos ( ω t ) i sin ( ω t ) {\displaystyle \exp \left({-i\omega t}\right)=\cos(\omega t)-i\sin(\omega t)} が用いられていた。ここで i虚数単位である。

それら二つの変換は密接に関連している。特に、フーリエ変換は(正規化の系数として 1 / 2 π {\displaystyle 1/{\sqrt {2\pi }}} を用いる時)、

F ( ω ) = H ( ω ) + H ( ω ) 2 i H ( ω ) H ( ω ) 2 {\displaystyle F(\omega )={\frac {H(\omega )+H(-\omega )}{2}}-i{\frac {H(\omega )-H(-\omega )}{2}}}

によって、ハートレー変換から得ることが出来る。つまり、フーリエ変換の実部と虚部は、それぞれハートレー変換の偶部分(even part)と奇部分(odd part)によって与えられることになる。

逆に、実数値関数 f(t) に対して、ハートレー変換はフーリエ変換の実部と虚部を用いることで

{ H f } = { F f } { F f } = { F f ( 1 + i ) } {\displaystyle \{{\mathcal {H}}f\}=\Re \{{\mathcal {F}}f\}-\Im \{{\mathcal {F}}f\}=\Re \{{\mathcal {F}}f\cdot (1+i)\}}

として得られる。ここで {\displaystyle \Re } および {\displaystyle \Im } は複素フーリエ変換の実部と虚部をそれぞれ表す。

性質

ハートレー変換は実線形作用素であり、対称である(そしてエルミートである)。その対称性および自己反転性により、ハートレー変換はユニタリ作用素であることが分かる(実際、直交である)。

ハートレー変換について、重畳積分定理と類似な次のような議論がある。もし二つの関数 x ( t ) {\displaystyle x(t)} および y ( t ) {\displaystyle y(t)} にそれぞれハートレー変換 X ( ω ) {\displaystyle X(\omega )} および Y ( ω ) {\displaystyle Y(\omega )} が存在するなら、それらの関数の畳み込み z ( t ) = x y {\displaystyle z(t)=x*y} にも同様にハートレー変換が存在し、次のように与えられる:

Z ( ω ) = { H ( x y ) } = 2 π ( X ( ω ) [ Y ( ω ) + Y ( ω ) ] + X ( ω ) [ Y ( ω ) Y ( ω ) ] ) / 2 {\displaystyle Z(\omega )=\{{\mathcal {H}}(x*y)\}={\sqrt {2\pi }}\left(X(\omega )\left[Y(\omega )+Y(-\omega )\right]+X(-\omega )\left[Y(\omega )-Y(-\omega )\right]\right)/2}

フーリエ変換と同様に、偶関数あるいは奇関数のハートレー変換はそれぞれ偶関数あるいは奇関数となる。

余弦正弦(cas)

余弦正弦(cas)関数の性質は、三角法およびその位相変換三角関数 cas ( t ) = 2 sin ( t + π / 4 ) {\displaystyle {\mbox{cas}}(t)={\sqrt {2}}\sin(t+\pi /4)} としての定義により従う。例えば、

2 cas ( a + b ) = cas ( a ) cas ( b ) + cas ( a ) cas ( b ) + cas ( a ) cas ( b ) cas ( a ) cas ( b ) {\displaystyle 2{\mbox{cas}}(a+b)={\mbox{cas}}(a){\mbox{cas}}(b)+{\mbox{cas}}(-a){\mbox{cas}}(b)+{\mbox{cas}}(a){\mbox{cas}}(-b)-{\mbox{cas}}(-a){\mbox{cas}}(-b)\,}

および

cas ( a + b ) = cos ( a ) cas ( b ) + sin ( a ) cas ( b ) = cos ( b ) cas ( a ) + sin ( b ) cas ( a ) {\displaystyle {\mbox{cas}}(a+b)=\cos(a){\mbox{cas}}(b)+\sin(a){\mbox{cas}}(-b)=\cos(b){\mbox{cas}}(a)+\sin(b){\mbox{cas}}(-a)\,}

などが得られ、微分は

cas ( a ) = d d a cas ( a ) = cos ( a ) sin ( a ) = cas ( a ) {\displaystyle {\mbox{cas}}'(a)={\frac {\mbox{d}}{{\mbox{d}}a}}{\mbox{cas}}(a)=\cos(a)-\sin(a)={\mbox{cas}}(-a)}

となる。

参考文献

  • Hartley, R. V. L., A more symmetrical Fourier analysis applied to transmission problems, Proc. IRE 30, 144–150 (1942).
  • Bracewell, R. N., The Fourier Transform and Its Applications (McGraw-Hill, 1965, 2nd ed. 1978, revised 1986) (also translated into Japanese and Polish)
  • Bracewell, R. N., The Hartley Transform (Oxford University Press, 1986) (also translated into German and Russian)
  • Bracewell, R. N., Aspects of the Hartley transform, Proc. IEEE 82 (3), 381-387 (1994).
  • Millane, R. P., Analytic properties of the Hartley transform, Proc. IEEE 82 (3), 413-428 (1994).
  • Villasenor, John D. (1994). “Optical Hartley transforms”. Proc. IEEE 82 (3): 391-399. doi:10.1109/5.272144. 
  • 吉村孝志、竹中博士、二宮市三:「高速ハートレ変換は高速フーリエ変換を超えるか?」、日本応用数理学会論文誌、Vol.3, No.3(1993), pp.245-255.
  1. ^ Villasenor 1994.