Английская Википедия:Hartley transform
Шаблон:Short description In mathematics, the Hartley transform (HT) is an integral transform closely related to the Fourier transform (FT), but which transforms real-valued functions to real-valued functions. It was proposed as an alternative to the Fourier transform by Ralph V. L. Hartley in 1942,[1] and is one of many known Fourier-related transforms. Compared to the Fourier transform, the Hartley transform has the advantages of transforming real functions to real functions (as opposed to requiring complex numbers) and of being its own inverse.
The discrete version of the transform, the discrete Hartley transform (DHT), was introduced by Ronald N. Bracewell in 1983.[2]
The two-dimensional Hartley transform can be computed by an analog optical process similar to an optical Fourier transform (OFT), with the proposed advantage that only its amplitude and sign need to be determined rather than its complex phase.[3] However, optical Hartley transforms do not seem to have seen widespread use.
Definition
The Hartley transform of a function <math>f(t)</math> is defined by:
<math display=block> H(\omega) = \left\{\mathcal{H}f\right\}(\omega) = \frac{1}{\sqrt{2\pi}}\int_{-\infty}^\infty f(t) \operatorname{cas}(\omega t) \, \mathrm{d}t\,, </math>
where <math>\omega</math> can in applications be an angular frequency and
<math display=block> \operatorname{cas}(t) = \cos(t) + \sin(t) = \sqrt{2} \sin (t+\pi /4) = \sqrt{2} \cos (t-\pi /4)\,, </math>
is the cosine-and-sine (cas) or Hartley kernel. In engineering terms, this transform takes a signal (function) from the time-domain to the Hartley spectral domain (frequency domain).
Inverse transform
The Hartley transform has the convenient property of being its own inverse (an involution):
<math display=block>f = \{\mathcal{H} \{\mathcal{H}f \}\}\,.</math>
Conventions
The above is in accord with Hartley's original definition, but (as with the Fourier transform) various minor details are matters of convention and can be changed without altering the essential properties:
- Instead of using the same transform for forward and inverse, one can remove the <math>{1}/{\sqrt{2\pi}}</math> from the forward transform and use <math>{1}/{2\pi}</math> for the inverse—or, indeed, any pair of normalizations whose product is Шаблон:Nowrap (Such asymmetrical normalizations are sometimes found in both purely mathematical and engineering contexts.)
- One can also use <math>2\pi\nu t</math> instead of <math>\omega t</math> (i.e., frequency instead of angular frequency), in which case the <math>{1}/{\sqrt{2\pi}}</math> coefficient is omitted entirely.
- One can use <math>\cos-\sin</math> instead of <math>\cos+\sin</math> as the kernel.
Relation to Fourier transform
This transform differs from the classic Fourier transform <math>F(\omega) = \mathcal{F} \{ f(t) \}(\omega)</math> in the choice of the kernel. In the Fourier transform, we have the exponential kernel, Шаблон:Nowrap where <math>\mathrm{i}</math> is the imaginary unit.
The two transforms are closely related, however, and the Fourier transform (assuming it uses the same <math>1/\sqrt{2\pi}</math> normalization convention) can be computed from the Hartley transform via:
<math display=block>F(\omega) = \frac{H(\omega) + H(-\omega)}{2} - \mathrm{i} \frac{H(\omega) - H(-\omega)}{2}\,.</math>
That is, the real and imaginary parts of the Fourier transform are simply given by the even and odd parts of the Hartley transform, respectively.
Conversely, for real-valued functions Шаблон:Nowrap the Hartley transform is given from the Fourier transform's real and imaginary parts:
<math display=block>\{ \mathcal{H} f \} = \Re \{ \mathcal{F}f \} - \Im \{ \mathcal{F}f \} = \Re \{ \mathcal{F}f \cdot (1+\mathrm{i}) \}\,,</math>
where <math>\Re</math> and <math>\Im</math> denote the real and imaginary parts.
Properties
The Hartley transform is a real linear operator, and is symmetric (and Hermitian). From the symmetric and self-inverse properties, it follows that the transform is a unitary operator (indeed, orthogonal).
Convolution using Hartley transforms is[4] <math display=block> f(x) * g(x) = \frac{F(\omega) G(\omega) + F(-\omega) G(\omega) + F(\omega) G(-\omega) - F(-\omega) G(-\omega)}{2} </math> where <math>F(\omega) = \{\mathcal{H}f\}(\omega)</math> and <math>G(\omega) = \{\mathcal{H} g\}(\omega)</math>
Similar to the Fourier transform, the Hartley transform of an even/odd function is even/odd, respectively.
cas
The properties of the Hartley kernel, for which Hartley introduced the name cas for the function (from cosine and sine) in 1942,[1][5] follow directly from trigonometry, and its definition as a phase-shifted trigonometric function Шаблон:Nowrap For example, it has an angle-addition identity of:
<math display=block> 2 \operatorname{cas} (a+b) = \operatorname{cas}(a) \operatorname{cas}(b) + \operatorname{cas}(-a) \operatorname{cas}(b) + \operatorname{cas}(a) \operatorname{cas}(-b) - \operatorname{cas}(-a) \operatorname{cas}(-b)\,. </math>
Additionally:
<math display=block> \operatorname{cas} (a+b) = {\cos (a) \operatorname{cas} (b)} + {\sin (a) \operatorname{cas} (-b)} = \cos (b) \operatorname{cas} (a) + \sin (b) \operatorname{cas}(-a)\,, </math>
and its derivative is given by:
<math display=block> \operatorname{cas}'(a) = \frac{d}{da} \operatorname{cas} (a) = \cos (a) - \sin (a) = \operatorname{cas}(-a)\,. </math>
See also
References
- Шаблон:Cite book (NB. Also translated into German and Russian.)
- Шаблон:Cite journal
- Шаблон:Cite journal
Further reading
- Шаблон:Cite book (NB. Contains extensive bibliography.)
- ↑ 1,0 1,1 Ошибка цитирования Неверный тег
<ref>
; для сносокHartley_1942
не указан текст - ↑ Ошибка цитирования Неверный тег
<ref>
; для сносокBracewell_1983
не указан текст - ↑ Ошибка цитирования Неверный тег
<ref>
; для сносокVillasenor_1994
не указан текст - ↑ Шаблон:Cite book Equation (4.54)
- ↑ Ошибка цитирования Неверный тег
<ref>
; для сносокBracewell_1999
не указан текст