Time–frequency representation

Last updated

A time–frequency representation (TFR) is a view of a signal (taken to be a function of time) represented over both time and frequency. [1] Time–frequency analysis means analysis into the time–frequency domain provided by a TFR. This is achieved by using a formulation often called "Time–Frequency Distribution", abbreviated as TFD.


TFRs are often complex-valued fields over time and frequency, where the modulus of the field represents either amplitude or "energy density" (the concentration of the root mean square over time and frequency), and the argument of the field represents phase.

Background and motivation

A signal, as a function of time, may be considered as a representation with perfect time resolution. In contrast, the magnitude of the Fourier transform (FT) of the signal may be considered as a representation with perfect spectral resolution but with no time information because the magnitude of the FT conveys frequency content but it fails to convey when, in time, different events occur in the signal.

TFRs provide a bridge between these two representations in that they provide some temporal information andsome spectral information simultaneously. Thus, TFRs are useful for the representation and analysis of signals containing multiple time-varying frequencies.

Formulation of TFRs and TFDs

One form of TFR (or TFD) can be formulated by the multiplicative comparison of a signal with itself, expanded in different directions about each point in time. Such representations and formulations are known as quadratic or "bilinear" TFRs or TFDs (QTFRs or QTFDs) because the representation is quadratic in the signal (see Bilinear time–frequency distribution). This formulation was first described by Eugene Wigner in 1932 in the context of quantum mechanics and, later, reformulated as a general TFR by Ville in 1948 to form what is now known as the Wigner–Ville distribution, as it was shown in [2] that Wigner's formula needed to use the analytic signal defined in Ville's paper to be useful as a representation and for a practical analysis. Today, QTFRs include the spectrogram (squared magnitude of short-time Fourier transform), the scaleogram (squared magnitude of Wavelet transform) and the smoothed pseudo-Wigner distribution.

Although quadratic TFRs offer perfect temporal and spectral resolutions simultaneously, the quadratic nature of the transforms creates cross-terms, also called "interferences". The cross-terms caused by the bilinear structure of TFDs and TFRs may be useful in some applications such as classification as the cross-terms provide extra detail for the recognition algorithm. However, in some other applications, these cross-terms may plague certain quadratic TFRs and they would need to be reduced. One way to do this is obtained by comparing the signal with a different function. Such resulting representations are known as linear TFRs because the representation is linear in the signal. An example of such a representation is the windowed Fourier transform (also known as the short-time Fourier transform) which localises the signal by modulating it with a window function, before performing the Fourier transform to obtain the frequency content of the signal in the region of the window.

Wavelet transforms

Wavelet transforms, in particular the continuous wavelet transform, expand the signal in terms of wavelet functions which are localised in both time and frequency. Thus the wavelet transform of a signal may be represented in terms of both time and frequency.

The notions of time, frequency, and amplitude used to generate a TFR from a wavelet transform were originally developed intuitively. In 1992, a quantitative derivation of these relationships was published, based upon a stationary phase approximation. [3]

Linear canonical transformation

Linear canonical transformations are the linear transforms of the time–frequency representation that preserve the symplectic form. These include and generalize the Fourier transform, fractional Fourier transform, and others, thus providing a unified view of these transforms in terms of their action on the time–frequency domain.

See also

Related Research Articles

Digital signal processing (DSP) is the use of digital processing, such as by computers or more specialized digital signal processors, to perform a wide variety of signal processing operations. The digital signals processed in this manner are a sequence of numbers that represent samples of a continuous variable in a domain such as time, space, or frequency. In digital electronics, a digital signal is represented as a pulse train, which is typically generated by the switching of a transistor.

<span class="mw-page-title-main">Fourier analysis</span> Branch of mathematics

In mathematics, Fourier analysis is the study of the way general functions may be represented or approximated by sums of simpler trigonometric functions. Fourier analysis grew from the study of Fourier series, and is named after Joseph Fourier, who showed that representing a function as a sum of trigonometric functions greatly simplifies the study of heat transfer.

Harmonic analysis is a branch of mathematics concerned with the representation of functions or signals as the superposition of basic waves, and the study of and generalization of the notions of Fourier series and Fourier transforms. In the past two centuries, it has become a vast subject with applications in areas as diverse as number theory, representation theory, signal processing, quantum mechanics, tidal analysis and neuroscience.

<span class="mw-page-title-main">Wavelet</span> Function for integral Fourier-like transform

A wavelet is a wave-like oscillation with an amplitude that begins at zero, increases or decreases, and then returns to zero one or more times. Wavelets are termed a "brief oscillation". A taxonomy of wavelets has been established, based on the number and direction of its pulses. Wavelets are imbued with specific properties that make them useful for signal processing.

<span class="mw-page-title-main">Fourier transform</span> Mathematical transform that expresses a function of time as a function of frequency

A Fourier transform (FT) is a mathematical transform that decomposes functions into frequency components, which are represented by the output of the transform as a function of frequency. Most commonly functions of time or space are transformed, which will output a function depending on temporal frequency or spatial frequency respectively. That process is also called analysis. An example application would be decomposing the waveform of a musical chord into terms of the intensity of its constituent pitches. The term Fourier transform refers to both the frequency domain representation and the mathematical operation that associates the frequency domain representation to a function of space or time.

<span class="mw-page-title-main">Frequency domain</span> Signal representation

In physics, electronics, control systems engineering, and statistics, the frequency domain refers to the analysis of mathematical functions or signals with respect to frequency, rather than time. Put simply, a time-domain graph shows how a signal changes over time, whereas a frequency-domain graph shows how much of the signal lies within each given frequency band over a range of frequencies. A frequency-domain representation can also include information on the phase shift that must be applied to each sinusoid in order to be able to recombine the frequency components to recover the original time signal.

In signal processing, time–frequency analysis comprises those techniques that study a signal in both the time and frequency domains simultaneously, using various time–frequency representations. Rather than viewing a 1-dimensional signal and some transform, time–frequency analysis studies a two-dimensional signal – a function whose domain is the two-dimensional real plane, obtained from the signal via a time–frequency transform.

Stransform as a time–frequency distribution was developed in 1994 for analyzing geophysics data. In this way, the S transform is a generalization of the short-time Fourier transform (STFT), extending the continuous wavelet transform and overcoming some of its disadvantages. For one, modulation sinusoids are fixed with respect to the time axis; this localizes the scalable Gaussian window dilations and translations in S transform. Moreover, the S transform doesn't have a cross-term problem and yields a better signal clarity than Gabor transform. However, the S transform has its own disadvantages: the clarity is worse than Wigner distribution function and Cohen's class distribution function.

<span class="mw-page-title-main">Chirplet transform</span>

In signal processing, the chirplet transform is an inner product of an input signal with a family of analysis primitives called chirplets.

In mathematics, in the area of harmonic analysis, the fractional Fourier transform (FRFT) is a family of linear transformations generalizing the Fourier transform. It can be thought of as the Fourier transform to the n-th power, where n need not be an integer — thus, it can transform a function to any intermediate domain between time and frequency. Its applications range from filter design and signal analysis to phase retrieval and pattern recognition.

<span class="mw-page-title-main">Wigner distribution function</span>

The Wigner distribution function (WDF) is used in signal processing as a transform in time-frequency analysis.

Geophysical survey is the systematic collection of geophysical data for spatial studies. Detection and analysis of the geophysical signals forms the core of Geophysical signal processing. The magnetic and gravitational fields emanating from the Earth's interior hold essential information concerning seismic activities and the internal structure. Hence, detection and analysis of the electric and Magnetic fields is very crucial. As the Electromagnetic and gravitational waves are multi-dimensional signals, all the 1-D transformation techniques can be extended for the analysis of these signals as well. Hence this article also discusses multi-dimensional signal processing techniques.

The method of reassignment is a technique for sharpening a time-frequency representation by mapping the data to time-frequency coordinates that are nearer to the true region of support of the analyzed signal. The method has been independently introduced by several parties under various names, including method of reassignment, remapping, time-frequency reassignment, and modified moving-window method. In the case of the spectrogram or the short-time Fourier transform, the method of reassignment sharpens blurry time-frequency data by relocating the data according to local estimates of instantaneous frequency and group delay. This mapping to reassigned time-frequency coordinates is very precise for signals that are separable in time and frequency with respect to the analysis window.

A Modified Wigner distribution function is a variation of the Wigner distribution function (WD) with reduced or removed cross-terms.

Bilinear time–frequency distributions, or quadratic time–frequency distributions, arise in a sub-field of signal analysis and signal processing called time–frequency signal processing, and, in the statistical analysis of time series data. Such methods are used where one needs to deal with a situation where the frequency composition of a signal may be changing over time; this sub-field used to be called time–frequency signal analysis, and is now more often called time–frequency signal processing due to the progress in using these methods to a wide range of signal-processing problems.

The cone-shape distribution function, also known as the Zhao–Atlas–Marks time-frequency distribution,, is one of the members of Cohen's class distribution function. It was first proposed by Yunxin Zhao, Les E. Atlas, and Robert J. Marks II in 1990. The distribution's name stems from the twin cone shape of the distribution's kernel function on the plane. The advantage of the cone kernel function is that it can completely remove the cross-term between two components having the same center frequency. Cross-term results from components with the same time center, however, cannot be completely removed by the cone-shaped kernel.

In the field of time–frequency analysis, several signal formulations are used to represent the signal in a joint time–frequency domain.

In reflection seismology, a seismic attribute is a quantity extracted or derived from seismic data that can be analysed in order to enhance information that might be more subtle in a traditional seismic image, leading to a better geological or geophysical interpretation of the data. Examples of seismic attributes can include measured time, amplitude, frequency and attenuation, in addition to combinations of these. Most seismic attributes are post-stack, but those that use CMP gathers, such as amplitude versus offset (AVO), must be analysed pre-stack. They can be measured along a single seismic trace or across multiple traces within a defined window.

Time–frequency analysis for music signals is one of the applications of time–frequency analysis. Musical sound can be more complicated than human vocal sound, occupying a wider band of frequency. Music signals are time-varying signals; while the classic Fourier transform is not sufficient to analyze them, time–frequency analysis is an efficient tool for such use. Time–frequency analysis is extended from the classic Fourier approach. Short-time Fourier transform (STFT), Gabor transform (GT) and Wigner distribution function (WDF) are famous time–frequency methods, useful for analyzing music signals such as notes played on a piano, a flute or a guitar.

Spectroscopic optical coherence tomography (SOCT) is an optical imaging and sensing technique, which provides localized spectroscopic information of a sample based on the principles of optical coherence tomography (OCT) and low coherence interferometry. The general principles behind SOCT arise from the large optical bandwidths involved in OCT, where information on the spectral content of backscattered light can be obtained by detection and processing of the interferometric OCT signal. SOCT signal can be used to quantify depth-resolved spectra to retrieve the concentration of tissue chromophores, characterize tissue light scattering, and/or used as a functional contrast enhancement for conventional OCT imaging.


  1. E. Sejdić, I. Djurović, J. Jiang, "Time-frequency feature representation using energy concentration: An overview of recent advances," Digital Signal Processing, vol. 19, no. 1, pp. 153-183, January 2009.
  2. B. Boashash, "Note on the use of the Wigner distribution for time frequency signal analysis", IEEE Trans. on Acoust. Speech. and Signal Processing, vol. 36, issue 9, pp 1518–1521, Sept. 1988. doi : 10.1109/29.90380
  3. Delprat, N., Escudii, B., Guillemain, P., Kronland-Martinet, R., Tchamitchian, P., and Torrksani, B. (1992). "Asymptotic wavelet and Gabor analysis: extraction of instantaneous frequencies". IEEE Transactions on Information Theory. 38 (2): 644–664. doi:10.1109/18.119728.{{cite journal}}: CS1 maint: multiple names: authors list (link)