Analytic signal

From Free net encyclopedia

In signal processing, the analytic signal, or analytic representation, of a signal <math>s(t)\,</math> is defined by:

<math>s_\mathrm{a}(t) = s(t) + j\cdot \hat s(t)\,</math>

where <math>\hat s(t)\,</math> is the Hilbert transform of <math>s(t)\,</math> and <math>j\,</math> (aka <math>i\,</math>) is the imaginary unit. The analytic representation facilitates mathematical manipulations, such as the derivation of single sideband modulation and demodulation. And as shown below, it makes certain attributes of <math>s(t)\,</math> more accessible.

The basic idea of the analytic representation is that the negative frequency components of the Fourier transform of a real-valued function are superfluous, due to spectral symmetry. They can be discarded with no loss of information, if we are willing to deal with a complex-valued function instead. And that turns out to have the advantages mentioned above. As long as the manipulated function has no negative frequency components (i.e., it is still analytic), the conversion from complex back to real is just a matter of discarding the imaginary part. And of course that also restores the suppressed negative frequency components.

  • Analytic signals are often shifted in frequency (down-converted) toward 0 Hz, which creates [non-symmetrical] negative frequency components. One motive is to minimize the highest frequency, which minimizes the necessary sample-rate for discrete-time signals. A frequency shift does not undermine the mathematical tractability of the complex signal representation. So in that sense, the down-converted signal is still "analytic". However, restoring the real-valued representation is no longer a simple matter of just extracting the real component. Up-conversion is obviously required, and if the signal has been sampled (discrete-time), interpolation (upsampling) might also be necessary to avoid aliasing.
  • The complex conjugate of an analytic signal contains only negative frequency components. In that case also, there is no loss of information or reversibility by discarding the imaginary component. Obviously the real component of the complex conjugate is the same as the real component of the analytic signal. But in this case, its extraction restores the suppressed positive frequency components.
  • Another way to achieve a spectrum of negative frequencies is to frequency-shift the analytic signal sufficiently far in the negative direction. Extracting the real component again restores the positive frequencies. But in this case their order is reversed... the low-frequency component is now the high one. This can be used to demodulate a type of single sideband signal called lower sideband or inverted sideband.


Example 1: <math>s(t) = \cos(\omega_0 t)\,</math>, for some parameter <math>\omega_0 > 0\,</math>
<math>\hat s(t) = \cos(\omega_0 t -\begin{matrix} \frac{\pi }{2} \end{matrix}) = \sin(\omega_0 t)\,</math>
<math>s_\mathrm{a}(t) = \cos(\omega_0 t) + j\cdot \sin(\omega_0 t) = e^{j \omega_0 t}\,</math> (The 2nd equality is Euler's formula.)
  • This is a complex-valued signal with increasing phase (positive frequency).
It also follows from Euler's formula that <math>\cos(\omega t) = \begin{matrix} \frac{1}{2} \end{matrix}(e^{j \omega t}+e^{-j \omega t})\,</math>.   So <math>s(t)\,</math> comprises both positive and negative frequency components. <math>s_\mathrm{a}(t)\,</math> is just the positive portion.


Example 2: <math>s(t) = \cos(\omega _1 t) + \cos(\omega _2 t)\,</math>
<math>\hat s(t) = \sin(|\omega _1 |t) + \sin(|\omega _2 |t)\,</math>
<math>s_\mathrm{a}(t) = e^{j |\omega _1 |t}+e^{j |\omega _2 |t}\,</math>


The definition can also be expressed in terms of the Fourier transforms of <math>s(t)\,</math> and <math>s_\mathrm{a}(t)\,</math>, respectively denoted by <math>S(\omega)\,</math> and <math>S_\mathrm{a}(\omega)\,</math>:

<math>S_\mathrm{a}(\omega)\,</math> <math> = 2 S(\omega)\,</math>,   for <math>\omega > 0\,</math>
<math> = 0\,</math>,   for <math>\omega < 0\,</math>


So it is a general result that <math>s_\mathrm{a}(t)\,</math> comprises only the positive frequency components of <math>s(t)\,</math>.

Yet the analytic transforms are reversible for real-valued <math>s(t)\,</math>. I.e., the original functions can be recovered:

<math>S(-\omega)\,</math> <math>= S^{*}(\omega) \,</math> ( Hermitian property )
<math>= \begin{matrix} \frac{1}{2} \end{matrix} S_a^{*}(\omega) \,</math>,   for <math>\omega > 0\,</math>
and
<math>s(t) = Re\{s_\mathrm{a}(t)\}\,</math>.


For completeness, we note that nothing prevents us from computing <math>s_\mathrm{a}(t)\,</math> for a complex-valued <math>s(t)\,</math>. The negative-frequency components will be eliminated, as usual (see example 3). But it might not be a reversible representation, because the original spectrum is not symmetrical in general. So except for this example, the general discussion assumes real-valued <math>s(t)\,</math>.

Example 3: <math>s(t) = e^{-j \omega_0 t}\,</math>, for some parameter <math>\omega_0 > 0\,</math>
<math>\hat s(t) = j\cdot e^{-j \omega_0 t}\,</math>
<math>s_\mathrm{a}(t) = e^{-j \omega_0 t} + j^2\cdot e^{-j \omega_0 t} = 0\,</math>


Contents

Applications

Image:Analytic.png

The analytic signal can also be expressed in terms of complex polar coordinates, <math>s_\mathrm{a}(t) = A(t)e^{j\phi(t)}\,</math>, where:

  • <math>A(t) = |s_\mathrm{a}(t)| = \sqrt {s^2(t) + \hat s^2(t)}\,</math>
  • <math>\phi(t) = \arg(s_\mathrm{a}(t)) = \arctan({\hat s(t) \over s(t)})\,</math>

These functions are respectively called the amplitude envelope and instantaneous phase of the signal <math>s(t)\,</math>. In the accompanying diagram, the blue curve depicts <math>s(t)\,</math>, and the red curve depicts the corresponding <math>A(t)\,</math>.
The time derivative of the instantaneous phase is called the instantaneous frequency:

  • <math>\omega_\mathrm{inst}(t) = \phi '(t) = {d \over dt} \phi(t)\,</math>

The amplitude function, and the instantaneous phase and frequency are in some applications used to measure and detect local features of the signal. Another application of the analytic representation of a signal relates to demodulation of modulated signals. The polar coordinates conveniently separate the effects of amplitude modulation and phase (or frequency) modulation, and effectively demodulates certain kinds of signals.

Extensions of the analytic signal to signals of multiple variables

The concept of analytic signal is well-defined for signals of a single variable which typically is time. For signals of two or more variables, an analytic signal can be defined in different ways, and two approaches are presented below

Multi-dimensional analytic signal based on an ad-hoc direction

A straightforward generalization of the analytic signal can be done for a multi-dimensional signal once it is established what is meant by negative frequencies for this case. This can be done by introducing a normalized vector <math>\hat{n}</math> in the Fourier domain and label any frequency vector <math>u</math> as negative if <math>u \cdot \hat{n} < 0</math>. The analytic signal is then produced by removing all negative frequencies and multiply the result by <math>2</math>, in accordance to the procedure described for the case of one-variable signals. However, it should be noted that there is no particular direction for <math>\hat{n}</math> which must be chosen unless there are some additional constraints. Therefore, the choice of <math>\hat{n}</math> is ad-hoc, or application specific.

The monogenic signal

The real and imaginary parts of the analytic signal correspond to the two elements of the vector-valued monogenic signal, as it is defined for one-variable signals. However, the monogenic signal can be extended to arbitrary number of variables in a straightforward manner, producing an <math>n + 1</math> dimensional vector-valued function for the case of <math>n</math> variable signals.

Not to be confused with

In mathematics an analytic function refers to a function which satisfies certain properties related to differentiability. The concept of analytic signal should not be confused with analytic functions.

See also

References

  • Bracewell, R; The Fourier Transform and Its Applications, 2nd ed, 1986, McGraw-Hill.
  • Leon Cohen, "Time-frequency analysis", Prentice-Hall (1995)

External links