World Library  
Flag as Inappropriate
Email this Article

Analytic signal

Article Id: WHEBN0000950777
Reproduction Date:

Title: Analytic signal  
Author: World Heritage Encyclopedia
Language: English
Subject: Hilbert transform, Quadrature filter, Instantaneous phase, Phasor, Fourier analysis
Collection: Fourier Analysis, Signal Processing, Time–frequency Analysis
Publisher: World Heritage Encyclopedia
Publication
Date:
 

Analytic signal

In mathematics and signal processing, an analytic signal is a complex-valued function that has no negative frequency components.[1]  The real and imaginary parts of an analytic signal are real-valued functions related to each other by the Hilbert transform.

The analytic representation of a real-valued function is an analytic signal, comprising the original function and its Hilbert transform. This representation facilitates many mathematical manipulations. The basic idea is that the negative frequency components of the Fourier transform (or spectrum) of a real-valued function are superfluous, due to the Hermitian symmetry of such a spectrum. These negative frequency components can be discarded with no loss of information, provided one is willing to deal with a complex-valued function instead. That makes certain attributes of the function more accessible and facilitates the derivation of modulation and demodulation techniques, such as single-sideband. As long as the manipulated function has no negative frequency components (that is, it is still analytic), the conversion from complex back to real is just a matter of discarding the imaginary part. The analytic representation is a generalization of the phasor concept:[2] while the phasor is restricted to time-invariant amplitude, phase, and frequency, the analytic signal allows for time-variable parameters.

Contents

  • Definition 1
  • Examples 2
    • Example 1 2.1
    • Example 2 2.2
    • Example 3 2.3
  • Negative frequency components 3
  • Applications 4
    • Envelope and instantaneous phase 4.1
    • Complex envelope/baseband 4.2
  • Extensions of the analytic signal to signals of multiple variables 5
    • Multi-dimensional analytic signal based on an ad hoc direction 5.1
    • The monogenic signal 5.2
  • See also 6
    • Applications 6.1
  • Notes 7
  • References 8
  • Further reading 9
  • External links 10

Definition

Transfer function to create an analytic signal

If s(t) is a real-valued function with Fourier transform S(f), then the transform has Hermitian symmetry about the f = 0 axis:

S(-f) = S(f)^*,   which is the complex conjugate of S(f).

The function:

\begin{align} S_\mathrm{a}(f) &\stackrel{\mathrm{def}} \begin{cases} 2S(f), &\text{for}\ f > 0,\\ S(f), &\text{for}\ f = 0,\\ 0, &\text{for}\ f < 0 \end{cases}\\ &= \underbrace{2 \operatorname{u}(f)}_{1 + \sgn(f)}S(f) = S(f) + \sgn(f)S(f), \end{align}
where:

contains only the non-negative frequency components of S(f).  And the operation is reversible, due to the Hermitian symmetry of S(f):

\begin{align} S(f) &= \begin{cases} \frac{1}{2}S_\mathrm{a}(f), &\text{for}\ f > 0,\\ S_\mathrm{a}(f), &\text{for}\ f = 0,\\ \frac{1}{2}S_\mathrm{a}(-f)^*, &\text{for}\ f < 0\ \text{(Hermitian symmetry)} \end{cases}\\ &= \frac{1}{2}[S_\mathrm{a}(f) + S_\mathrm{a}(-f)^*]. \end{align}


The analytic signal of s(t) is the inverse Fourier transform of S_\mathrm{a}(f):

\begin{align} s_\mathrm{a}(t) &\stackrel{\mathrm{def}} \mathcal{F}^{-1}[S_\mathrm{a}(f)]\\ &= \mathcal{F}^{-1}[S (f)+ \sgn(f) \cdot S(f)]\\ &= \underbrace{\mathcal{F}^{-1}\{S(f)\}}_{s(t)} + \overbrace{ \underbrace{\mathcal{F}^{-1}\{\sgn(f)\}}_{j\frac{1}{\pi t}} * \underbrace{\mathcal{F}^{-1}\{S(f)\}}_{s(t)} }^{convolution}\\ &= s(t) + j\underbrace{\left_{\operatorname{\mathcal{H}}[s(t)]}\\ &= s(t) + j\hat{s}(t), \end{align}

where

Examples

Example 1

s(t) = \cos(\omega t),   where  \omega > 0.

Then:

\hat{s}(t) = \cos(\omega t - \pi/2) = \sin(\omega t),
s_\mathrm{a}(t) = s(t) + j\hat{s}(t) = \cos(\omega t) + j\sin(\omega t) = e^{j\omega t}.  The third equality is Euler's formula.


A corollary of Euler's formula is  \cos(\omega t) = \tfrac{1}{2} (e^{j\omega t} + e^{j (-\omega) t}).  In general, the analytic representation of a simple sinusoid is obtained by expressing it in terms of complex-exponentials, discarding the negative frequency component, and doubling the positive frequency component. And the analytic representation of a sum of sinusoids is the sum of the analytic representations of the individual sinusoids.

Example 2

Here we use Euler's formula to identify and discard the negative frequency.

s(t) = \cos(\omega t+\theta) = \tfrac{1}{2} (e^{j (\omega t+\theta)} + e^{-j (\omega t+\theta)})

Then:

s_\mathrm{a}(t) = \begin{cases} e^{j(\omega t + \theta)} \ \ = \ e^{j |\omega| t}\cdot e^{j\theta} , & \text{if} \ \omega > 0, \\ e^{-j(\omega t + \theta)} = \ e^{j |\omega| t}\cdot e^{-j\theta} , & \text{if} \ \omega < 0. \end{cases}

Example 3

This is another example of using the Hilbert transform method to remove negative frequency components. We note that nothing prevents us from computing s_\mathrm{a}(t) for a complex-valued s(t). 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 s(t).

s(t) = e^{-j\omega t}, where \omega > 0.

Then:

\hat{s}(t) = je^{-j\omega t},
s_\mathrm{a}(t) = e^{-j\omega t} + j^2 e^{-j\omega t} = e^{-j\omega t} - e^{-j\omega t} = 0.

Negative frequency components

Since s(t) = \operatorname{Re}[s_\mathrm{a}(t)], restoring the negative frequency components is a simple matter of discarding \operatorname{Im}[s_\mathrm{a}(t)] which may seem counter-intuitive. We can also note that the complex conjugate s_\mathrm{a}^*(t) comprises only the negative frequency components. And therefore \operatorname{Re}[s_\mathrm{a}^*(t)] restores the suppressed positive frequency components.

Applications

Envelope and instantaneous phase

A function in blue and the magnitude of its analytic representation in red, showing the envelope effect.

An analytic signal can also be expressed in terms of its time-variant magnitude and phase (polar coordinates):

s_\mathrm{a}(t) = s_\mathrm{m}(t)e^{j\phi(t)},

where:

  • s_\mathrm{m}(t) \stackrel{\mathrm{def}} |s_\mathrm{a}(t)| is called the instantaneous amplitude or the envelope;
  • \phi(t) \stackrel{\mathrm{def}} \arg\!\left[s_\mathrm{a}(t)\right] is called the instantaneous phase.

In the accompanying diagram, the blue curve depicts s(t) and the red curve depicts the corresponding s_\mathrm{m}(t).

The time derivative of the unwrapped instantaneous phase has units of radians/second, and is called the instantaneous angular frequency:

\omega(t) \stackrel{\mathrm{def}} \frac{d\phi}{dt}(t).

The instantaneous frequency (in hertz) is therefore:

f(t)\stackrel{\mathrm{def}} \frac{1}{2\pi}\omega(t).  [3]

The instantaneous amplitude, 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.

Complex envelope/baseband

Analytic signals are often shifted in frequency (down-converted) toward 0 Hz, possibly creating [non-symmetrical] negative frequency components:

\underline{s_\mathrm{a}}(t) \stackrel{\mathrm{def}} s_\mathrm{a}(t)e^{-j\omega_0 t} = s_\mathrm{m}(t)e^{j(\phi(t) - \omega_0 t)},

where \omega_0 is an arbitrary angular frequency.[2]

This function goes by various names, such as complex envelope and complex baseband. The complex envelope is not unique; it is determined by the choice of \omega_0. This concept is often used when dealing with passband signals. If s(t) is a modulated signal, \omega_0 might be equated to its carrier frequency.

In other cases, \omega_0 is selected to be somewhere in the middle of the desired passband. Then a simple low-pass filter with real coefficients can excise the portion of interest. Another motive is to reduce the highest frequency, which reduces the minimum rate for alias-free sampling. 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 may be required, and if the signal has been sampled (discrete-time), interpolation (upsampling) might also be necessary to avoid aliasing.

If \omega_0 is chosen larger than the highest frequency of s_\mathrm{a}(t), then \underline{s_\mathrm{a}}(t) has no positive frequencies. In that case, extracting the real component restores them, but in reverse order; the low-frequency components are now high ones and vice versa. This can be used to demodulate a type of single sideband signal called lower sideband or inverted sideband.

Other choices of \omega_0

Sometimes \omega_0 is chosen to minimize

\int_{0}^{+\infty}(\omega - \omega_0)^2|S_\mathrm{a}(\omega)|^2\, d\omega.

Alternatively,[4] \omega_0 can be chosen to minimize the mean square error in linearly approximating the unwrapped instantaneous phase \phi(t):

\int_{-\infty}^{+\infty}[\omega(t) - \omega_0]^2 |s_\mathrm{a}(t)|^2\, dt

or another alternative (for some optimum \theta):

\int_{-\infty}^{+\infty}[\phi(t) - (\omega_0 t + \theta)]^2\, dt.

In the field of time-frequency signal processing, it was shown that the analytic signal was needed in the definition of the Wigner–Ville distribution so that the method can have the desirable properties needed for practical applications.[5]

Sometimes the complex envelope is identified as synonymous with complex amplitude;[1][2] other times it is presented as a time-dependent generalisation.[3] Their relationship is not unlike that in the real-valued case: varying envelope generalizing constant amplitude.

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 unit vector \boldsymbol \hat{u} in the Fourier domain and label any frequency vector \boldsymbol \xi as negative if \boldsymbol \xi \cdot \boldsymbol \hat{u} < 0. The analytic signal is then produced by removing all negative frequencies and multiply the result by 2, in accordance to the procedure described for the case of one-variable signals. However, there is no particular direction for \boldsymbol \hat{u} which must be chosen unless there are some additional constraints. Therefore, the choice of \boldsymbol \hat{u} 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 (n+1)-dimensional vector-valued function for the case of n-variable signals.

See also

Applications

Notes

  1. ^ "the complex envelope (or complex amplitude)"[6]
  2. ^ "the complex envelope (or complex amplitude)", p.586 [7]
  3. ^ "Complex envelope is an extended interpretation of complex amplitude as a function of time." p.85[8]

References

  1. ^ ``Mathematics of the Discrete Fourier Transform (DFT), with Audio Applications --- Second Edition, by Julius O. Smith III, W3K Publishing, 2007, ISBN 978-0-9745607-4-8. Copyright © 2014-04-21 by Julius O. Smith III Center for Computer Research in Music and Acoustics (CCRMA), Stanford University, https://ccrma.stanford.edu/~jos/r320/Analytic_Signals_Hilbert_Transform.html[7/16/2014 1:07:57 PM]
  2. ^ a b Bracewell, Ron. The Fourier Transform and Its Applications. McGraw-Hill, 1965. p269
  3. ^ B. Boashash, "Estimating and Interpreting the Instantaneous Frequency of a Signal-Part I: Fundamentals", Proceedings of the IEEE, Vol. 80, No. 4, pp. 519-538, April 1992
  4. ^ [5]
  5. ^ B. Boashash, “Notes on the use of the Wigner distribution for time frequency signal analysis”, IEEE Trans. on Acoustics, Speech, and Signal Processing , vol. 26, no. 9, 1987
  6. ^ Time-Frequency Analysis edited by Franz Hlawatsch, François Auger [6]
  7. ^ Encyclopedia of Optical Engineering, Volume 1, edited by Ronald G. Driggers [7]
  8. ^ Global Environment Remote Sensing edited by Kenʼichi Okamoto [8]

Further reading

  • Leon Cohen, Time-frequency analysis, Prentice Hall, Upper Saddle River, 1995.
  • Frederick W. King, Hilbert Transforms, vol. II, Cambridge University Press, Cambridge, 2009.
  • B. Boashash, Time-Frequency Signal Analysis and Processing: A Comprehensive Reference, Elsevier Science, Oxford, 2003.

External links

  • Analytic Signals and Hilbert Transform Filters
This article was sourced from Creative Commons Attribution-ShareAlike License; additional terms may apply. World Heritage Encyclopedia content is assembled from numerous content providers, Open Access Publishing, and in compliance with The Fair Access to Science and Technology Research Act (FASTR), Wikimedia Foundation, Inc., Public Library of Science, The Encyclopedia of Life, Open Book Publishers (OBP), PubMed, U.S. National Library of Medicine, National Center for Biotechnology Information, U.S. National Library of Medicine, National Institutes of Health (NIH), U.S. Department of Health & Human Services, and USA.gov, which sources content from all federal, state, local, tribal, and territorial government publication portals (.gov, .mil, .edu). Funding for USA.gov and content contributors is made possible from the U.S. Congress, E-Government Act of 2002.
 
Crowd sourced content that is contributed to World Heritage Encyclopedia is peer reviewed and edited by our editorial staff to ensure quality scholarly research articles.
 
By using this site, you agree to the Terms of Use and Privacy Policy. World Heritage Encyclopedia™ is a registered trademark of the World Public Library Association, a non-profit organization.
 


Copyright © World Library Foundation. All rights reserved. eBooks from Project Gutenberg are sponsored by the World Library Foundation,
a 501c(4) Member's Support Non-Profit Organization, and is NOT affiliated with any governmental agency or department.