Coherence (signal processing)

Contents

Background

The spectral coherence is a statistic that can be used to examine the relation between two signals or data sets. It is commonly used to estimate the power transfer between input and output of a linear system. If the signals are ergodic, and the system function linear, it can be used to estimate the causality between the input and output.

Definition and Formulation

The coherence (sometimes called magnitude-squared coherence) between two signals x(t) and y(t) is a real-valued function that is defined as [1][2]:

C_{xy} = \frac{|G_{xy}|^2}{G_{xx} G_{yy}}

where Gxy is the cross-spectral density between x and y, and Gxx and Gyy the autospectral density of x and y respectively. The magnitude of the spectral density is denoted as |G|. Given the restrictions noted above (ergodicity, linearity) the coherence function estimates the extent to which y(t) may be predicted from x(t) by an optimum linear least squares function.

Values of coherence will always satisfy 0\le C_{xy}\le 1. For a constant parameter linear system with a single input x(t) and single output y(t), the coherence will be equal to one. If the system's impulse response is h(t) (y(t) = h(t) * x(t), where * denotes convolution), then Y(f) = H(f)X(f), where Y(f) is the Fourier transform of y(t). Since Gyy = | H(f) | 2Gxx(f) and Gxy = H(f)Gxx(f), and since Gxx(f) is real,

C_{xy} = \frac{|H(f)G_{xx}(f)|^2}{G_{xx}G_{yy}} = \frac{|H(f)G_{xx}(f)|^2}{G_{xx}^2|H(f)|^2} = \frac{|G_{xx}|^2}{G_{xx}^2} = 1.

If x(t) and y(t) are completely unrelated the coherence will be zero. If Cxy is less than one but greater than zero it is an indication that either: noise is entering the measurements, that the assumed function relating x(t) and y(t) is not linear, or that y(t) is producing output due to input x(t) as well as other inputs.

The coherence of a linear system therefore represents the fractional part of the output signal power that is produced by the input at that frequency. We can also view the quantity 1 − Cxy as an estimate of the fractional power of the output that is not contributed by the input at a particular frequency. This leads naturally to definition of the coherent output spectrum:

Gvv = CxyGyy

Gvv provides a spectral quantification of the output power that is uncorrelated with noise or other inputs.

Example

Figure 1: Barometric pressure (black), ocean water levels (red), and groundwater well level (blue) near Lake Worth Florida during hurricane Frances.

Consider the two signals shown in the lower portion of figure 1. There appears to be a close relationship between the ocean surface water levels and the groundwater well levels. It is also clear that the barometric pressure has an effect on both the ocean water levels and groundwater levels.

Figure 2. shows the autospectral density of ocean water level over a long period of time.

Figure 2: Autospectral density of ocean water level. The tidal components are labeled.

As expected, most of the energy is centered on the well-known tidal frequencies. Likewise, the autospectral density of groundwater well levels are shown in figure 3.

Figure 3: Autospectral density of groundwater well level.

It is clear that variation of the groundwater levels have significant power at the ocean tidal frequencies. To estimate the extent at which the groundwater levels are influenced by the ocean surface levels, we compute the coherence between them. Let us assume that there is a linear relationship between the ocean surface height and the groundwater levels. We further assume that the ocean surface height controls the groundwater levels so that we take the ocean surface height as the input variable, and the groundwater well height as the output variable.

The computed coherence is shown in figure 4 (note that coherence is often denoted as γ2.)

Figure 4: Coherence between ocean water level (input) and groundwater well level (output).

This result indicates that at most of the major ocean tidal frequencies the variation of groundwater level at this particular site is over 90% due to the forcing of the ocean tides. However, one must exercise caution in attributing causality. If the relation (transfer function) between the input and output is nonlinear, then values of the coherence can be erroneous. Another common mistake is to assume a causal input/output relation between observed variables, when in fact the causative mechanism is not in the system model. For example, it is clear that the atmospheric barometric pressure induces a variation in both the ocean water levels and the groundwater levels, but the barometric pressure is not included in the system model as an input variable. We have also assumed that the ocean water levels drive or control the groundwater levels. In reality it is a combination of hydrological forcing and tidal potential that are driving both the observed input and output signals. Additionally, noise introduced in the measurement process, or by the spectral signal processing can contribute to the coherence.

Extension to non-stationary signals

If the signals are non-stationary, (and therefore not ergodic), the above formulations may not be appropriate. For such signals, the concept of coherence has been extended by using the concept of time-frequency distributions to represent the time-varying spectral variations of non-stationary signals in lieu of traditional spectra. For more details, see [3].

References

  1. ^ J. S. Bendat, A. G. Piersol, Random Data, Wiley-Interscience, 1986
  2. ^ http://www.fil.ion.ucl.ac.uk/~wpenny/course/course.html, chapter 7
  3. ^ L.B. White and B. Boashash, "Cross Spectral Analysis of Non-Stationary Processes", IEEE Transactions on Information Theory, Vol. 36, No. 4, pp. 830-835, July 1990.



Wikimedia Foundation. 2010.

Look at other dictionaries:

  • Coherence — Coherence, coherency, or coherent can refer to: Contents 1 In physics 2 In mathematics 3 In financial economics 4 In philosophy …   Wikipedia

  • Coherence (statistics) — In probability theory and statistics, coherence can have two meanings. When dealing with personal probability assessments, or supposed probabilities derived in nonstandard ways, it is a property of self consistency across a whole set of such… …   Wikipedia

  • Signal-to-noise ratio — For signal to noise ratio in statistics, see Cohen s d. Signal to noise ratio (often abbreviated SNR or S/N) is a measure used in science and engineering that compares the level of a desired signal to the level of background noise. It is defined… …   Wikipedia

  • Coherence (physics) — In physics, coherence is a property of waves that enables stationary (i.e. temporally and spatially constant) interference. More generally, coherence describes all properties of the correlation between physical quantities of a wave. When… …   Wikipedia

  • Optical coherence tomography — Intervention Optical Coherence Tomography (OCT) image of a sarcoma MeSH …   Wikipedia

  • Mutual coherence (linear algebra) — In linear algebra, the coherence[1] or mutual coherence[2] of a matrix A is defined as the maximum absolute value of the cross correlations between the columns of A. Formally, let be the columns of the matrix A, which are assumed to be normalized …   Wikipedia

  • Central processing unit — CPU redirects here. For other uses, see CPU (disambiguation). An Intel 80486DX2 CPU from above An Intel 80486DX2 from below …   Wikipedia

  • Degree of coherence — In optics, correlation functions are used to characterize the statistical and coherence properties of an electromagnetic field. The degree of coherence is the normalized correlation of electric fields. In its simplest form, termed g(1), it is… …   Wikipedia

  • Cross-correlation — In signal processing, cross correlation is a measure of similarity of two waveforms as a function of a time lag applied to one of them. This is also known as a sliding dot product or sliding inner product. It is commonly used for searching a long …   Wikipedia

  • optics — /op tiks/, n. (used with a sing. v.) the branch of physical science that deals with the properties and phenomena of both visible and invisible light and with vision. [1605 15; < ML optica < Gk optiká, n. use of neut. pl. of OPTIKÓS; see OPTIC,… …   Universalium

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”

We are using cookies for the best presentation of our site. Continuing to use this site, you agree with this.