HILBERT TRANSFORM

In mathematics and in signal processing, the 'Hilbert transform' of a real-valued function, $s(t),,$ is another real-valued function in the same domain.  Together they form the very useful analytic representation of $s(t).,$  The Hilbert transform is named after the renowned mathematician David Hilbert.

Contents

Introduction

Inverse Hilbert transform

Hilbert transform examples

Narrowband model

Analytic representation

Practical considerations

Discrete Hilbert transform

Notes

See also

References

External links

Introduction

The Hilbert transform occurs in practice as the imaginary part of an analytic signal that has been created or manipulated. But it is also computed when necessary.
In the time domain, it is the convolution of  $s(t),,$  with the function  $h(t)\; =\; 1/(pi\; t),,$  which is the impulse response of a ''Hilbert filter''':'
:$$
egin{align}
widehat s(t) stackrel{mathrm{def}}{=} h(t)
★ s(t) & = int_{-infty}^{infty}h( au)cdot s(t- au), d au \
& = int_{-infty}^{infty}s( au)cdot h(t- au), d au.,
end{align}
    ^[1]
Another notation for this convolution is $mathcal\{H\}\{s(t)\}.,$
In the frequency domain, the Fourier transform of $widehat\; s(t),$ is the product of the transforms of $h(t),$ and $s(t),$':'
:$$
egin{align}
mathcal{F}{widehat s(t)} & = mathcal{F}{h(t)}cdot mathcal{F}{s(t)} \
& = H(omega) cdot S(omega),
end{align}

where':'

★
e^{+ipi/2} = 0+i, & mbox{for } omega < 0,\
0, & mbox{for } omega = 0,\
e^{-ipi/2} = 0-i, & mbox{for } omega > 0,
end{cases}

★ $mathcal\{F\}$ denotes the Fourier transform,

★ ''i'' (sometimes denoted by ''j'' ) is the imaginary unit, and

★ $omega\; ,$ is the angular frequency.
In other words, the Hilbert transform has the effect of shifting the negative frequency components of $s(t),$ by +90° ($pi/2$ radians) and the positive frequency components by −90°.

Inverse Hilbert transform

We also note that $H^2(omega\; )\; =\; e^\{pm\; i\; pi\}\; =\; -1.,$  So multiplying the above equation by $-H(omega\; ),$ gives':'
:$S(omega)\; =\; -H(omega\; )cdot\; mathcal\{F\}\{widehat\; s(t)\}$
which is a convolution in the time domain':'
:$s(t)\; =\; -\; (h(t)$
★ widehat s(t)) = -mathcal{H}{widehat s(t)} = -mathcal{H}{mathcal{H}{s(t)}}.,

Hilbert transform examples

{| class="wikitable"
|-
! Signal
$s(t),$ !! Hilbert transform
$mathcal\{H\}\{s(t)\}$
|-
| align="center"| $sin(t),$ || align="center"| $-cos(t),$
|-
| align="center"| $cos(t),$ || align="center"| $sin(t),$
|-
| align="center"| $1\; over\; t^2\; +\; 1$ || align="center"| $t\; over\; t^2\; +\; 1$
|-
| align="center"| 'Sinc function'
$sin(t)\; over\; t$|| align="center"| $1-\; cos(t)over\; t$
|-
| align="center"| 'Rectangular function'
$sqcap(t)$|| align="center"| $\{1\; over\; pi\}\; ln\; left\; |\; \{t+\{1\; over\; 2\}\; over\; t-\{1\; over\; 2\}\}$
ight |
|-
| align="center"| 'Dirac delta function'
$delta(t)\; ,$ || align="center"| $\{1\; over\; pi\; t\}$
|}
'Notice:' Some authors, e.g., Bracewell, use our $-mathcal\{H\}$ as their definition of the forward transform. A consequence is that the right column of this table would be negated.

Narrowband model

Many signals can be accurately modeled as the product of a bandlimited "message" waveform, $s\_m(t),$, and a sinusoidal "carrier"':'
:$$
s(t) = s_m(t) cdot cos(omega t + arphi),

When $s\_m(t),$ has no frequency content above the carrier frequency, Hz, then':'
:$$
widehat{s}(t) = s_m(t) cdot sin(omega t + arphi)

So, the Hilbert transform may be as simple as a circuit that produces a 90° phase shift at the carrier frequency. Furthermore':'
:
ight),
from which one can reconstruct the carrier waveform. Then the message can be extracted from  $s(t),$ by coherent demodulation.

Analytic representation

The analytic representation of a signal is defined in terms of the Hilbert transform':'
:$s\_a(t)\; =\; s(t)\; +\; icdot\; widehat\; s(t),$
E.g., for the narrowband model [above], the analytic representation is':'
:{|
|$s\_a(t),$
|
|-
|
|
ight],
|-
|
|   (by Euler's formula)
|}
This complex heterodyne operation shifts all the frequency components of $s\_m(t),$ above 0 Hz. In that case, the imaginary part of the result is a Hilbert transform of the real part. This is an indirect way to produce Hilbert transforms.
While the analytic representation of a signal is not necessarily an analytic function, there is a connection to analytic functions here, which is in fact how the Hilbert transform arose historically. The idea is as follows. Starting with a function
:$f:mathbb\{R\}\; omathbb\{R\}$
one may extend this to a harmonic function on $mathbb\{R\}^2\_+$, the upper half plane, by convolving with the Poisson kernel. Every harmonic function is the real part of some analytic function. We then consider the imaginary part of this analytic function, specifically its values along the boundary. It turns out that the boundary values are $mathcal\{H\}(f)$. It then follows that the analytic function may be described as the Poisson integral of $f+imathcal\{H\}(f)$.

Practical considerations

The function ''h'' with ''h''(''t'') = 1/(π ''t'') is a non-causal filter and therefore cannot be implemented as is, if $s$ is a time-dependent signal. If $s$ is a function of a non-temporal variable, e.g., spatial, the non-causality might not be a problem. The filter is also of infinite support which may be a problem in certain applications. Another issue relates to what happens with the zero frequency (DC), which can be avoided by assuring that $s$ does not contain a DC-component.
A practical implementation in many cases implies that a finite support filter, which in addition is made causal by means of a suitable delay, is used to approximate the computation. The approximation may also imply that only a specific frequency range is subject to the characteristic phase shift related to the Hilbert transform. See also quadrature filter.

Discrete Hilbert transform

If the signal $s(t),$ is bandlimited, then $widehat\; s(t)$ is bandlimited in the same way. Consequently, both these signals can be sampled according to the sampling theorem, resulting in the discrete signals $s[n],$ and $widehat\{s\}[n]$.  The relation between the two discrete signals is then given by the convolution':'
:$widehat\{s\}[n]\; =\; h[n]$
★ s[n],
where
:$$
h[n]=
egin{cases}
0, & mbox{for }nmbox{ even},\
rac2{pi n} & mbox{for }nmbox{ odd}
end{cases}

which is non-causal and has infinite duration. In practice, a shortened and time-shifted approximation is used. The usual filter design tradeoffs apply (e.g. filter-order and latency vs. frequency-response). Also notice, that $h[n],$ is not just a sampled version of the ''Hilbert filter'' $h(t),$, defined above. Rather it is a sequence with this discrete-time Fourier transform':'
:$$
H(omega) =
egin{cases}
e^{+ipi/2}, & -pi leq omega < 0 \
e^{-ipi/2}, & 0 leq omega < pi
end{cases}


We note that a sequence similar to $h[n],$ can be generated by sampling $H(omega),$ and computing the inverse discrete Fourier transform. The larger the transform (i.e. more samples per $2\; pi$ radians), the better the agreement (for a given value of the abscissa, n). The figure shows the comparison for a 512-point transform. (Due to odd-symmetry, only half the sequence is actually plotted.)

But that is not the actual point, because it is easier and more accurate to generate $h[n],$ directly from the formula. The point is that many applications choose to avoid the convolution by doing the equivalent frequency-domain operation':'  simple multiplication of the signal transform with $H(omega),$, made even easier by the fact that the real and imaginary components are 0 and ±1 respectively. The attractiveness of that approach is only apparent when the actual Fourier transforms are replaced by samples of the same, i.e., the DFT, which is an approximation and introduces some distortion. Thus, after transforming back to the time-domain, those applications have indirectly generated (and convolved with) not $h[n],$, but the DFT approximation to it, which is shown in the figure.
Notes on ''fast convolution''':'

★ Implied in the technique described above is the concept of dividing a long signal into segments of arbitrary size. The signal is filtered piecewise, and the outputs are subsequently pieced back together.

★ The segment size is an important factor in controlling the amount of distortion. As the size increases, the DFT becomes more dense and is a better approximation to the underlying Fourier transform. In the time-domain, the same distortion is manifested as "aliasing", which results in a type of convolution called circular. It is as if the same segment is repeated periodically and filtered, resulting in distortion that is worst at either or both edges of the original segment. Increasing the segment size reduces the number of edges in the pieced-together result and therefore reduces overall distortion.

★ Another mitigation strategy is to simply discard the most badly distorted output samples, because data loss can be avoided by overlapping the input segments. When the filter's impulse response is less than the segment length, this can produce a distortion-free (non-circular) convolution. That of course requires an FIR filter, which the Hilbert transform is not. So yet another technique is to design an FIR approximation to a Hilbert transform filter. That moves the source of distortion from the convolution to the filter, where it can be readily characterized in terms of imperfections in the frequency response.

★ Failure to appreciate or correctly apply these concepts is probably one of the most common mistakes made by non-experts in the digital signal processing field.

Notes

1. To avoid the singularities in $h(t),,$ the integral is considered as a Cauchy principal value. It can be shown that if $sin\; L^p(mathbb\{R\})$, then $widehat\; s(t)$ is defined and in $L^p(mathbb\{R\})$ for .

See also

★ Analytic signal

★ Single-sideband signal

★ Harmonic conjugate

★ Kramers-Kronig relations

References

★ The Fourier Transform and Its Applications, Bracewell, R., , , McGraw-Hill, 1986,

★ Communication Systems, Carlson, Crilly, and Rutledge, , , , 2002,

External links

★ Hilbert transform

★ another exposition — Hilbert transform properties

★ Mathworld Hilbert transform — Contains a table of transforms

★ Analytic Signals and Hilbert Transform Filters

★ Easy Fourier Analysis hints to compute Hilbert transform in Time domain.