# Autocorrelation - book

Series | Geophysical References Series |
---|---|

Title | Digital Imaging and Deconvolution: The ABCs of Seismic Exploration and Processing |

Author | Enders A. Robinson and Sven Treitel |

Chapter | 7 |

DOI | http://dx.doi.org/10.1190/1.9781560801610 |

ISBN | 9781560801481 |

Store | SEG Online Store |

The *autocorrelation* of a signal is defined to be

**(**)

This formula gives the autocorrelation coefficient for each integer *k*. For a finite-length causal wavelet , the expression becomes

**(**)

The superscript asterisk indicates the complex conjugate. The superscript asterisk can be omitted except in cases in which the wavelet coefficients are complex.

Let us examine wavelets, all of which have the same autocorrelation. Figures 14 through 17 depict representative members of a suite of wavelets, all with the same autocorrelation.

Energy buildup | ||||
---|---|---|---|---|

Wavelet type | ||||

Minimum-delay wavelet | 4 | 5 | 5.25 | 5.3125 |

Mixed-delay four-length wavelet | 1 | 5 | 5.0625 | 5.3125 |

Maximum-delay four-length wavelet | 0.0625 | 0.3125 | 1.3125 | 5.3125 |

Next we shall show that the autocorrelation is the convolution of a signal with its zero-point reverse. To autocorrelate a wavelet, we merely convolve the wavelet with its reverse. For example, let us find the autocorrelation of the wavelet . Its reverse is the anticausal wavelet . Therefore, the folding table for such an operation is

**(**)

We see that the autocorrelation for time index 0 is equal to the sum on the main southwest-northeast diagonal:

**(**)

For each summand, the first index minus the second index is 0. The value gives the Energy of the wavelet. Working up from the diagonal, we get the for the negative *k*s; that is,

**(**)

For each summand, the first index minus the second index is –1. In addition, we see that

**(**)

The first index minus the second index is –2. Working down from the diagonal, we get the for positive *k*; that is,

**Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \begin{align} &r_{{ l}} =b_{2}b^*_{1}+b_{1}b^{*}_0. \end{align}}****(**)

For each summand, the first index minus the second is 1. In addition, we see that

**(**)

The first index minus the second index is 2. The wavelet **Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle b=\left\{b_0{,\ }b_{1}{\ ,\ }b_{2}\right\} }**
is called the *component signal* of the autocorrelation. In other words, an autocorrelation is given by the convolution of a component signal and its zero-point reverse. A component signal for an autocorrelation is not unique.

Now let us show that a wavelet and its reverse have the same autocorrelation. Convolution is described as a folding operation between two signals. The autocorrelation of a signal is defined as the convolution of the signal with its zero-point reverse. Consider the wavelet (2, 1) in which the value 2 occurs at time 0. Its zero-point reverse is (1, 2) in which the value 2 occurs at time 0, so its autocorrelation table is

**Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \begin{align} &\ \ \ \underline{1}\,\,\,\underline{2}\\ &2|2\,\,\,\,4.\\ &1|1\,\,\,\,2 \end{align}}****(**)

Its autocorrelation is

**(**)

On the other hand, consider the wavelet (1, 2), in which the 1 occurs at time 0. Its zero-point reverse is (2, 1), in which the 1 occurs at time 0. The autocorrelation table is

**(**)

and its autocorrelation is the same. Thus, the wavelet (2, 1) has the same autocorrelation as does its reverse (1, 2). This is always so. That is, the autocorrelation of any wavelet is the same as the autocorrelation of the reverse wavelet.

Now we shall illustrate an all-pass filter. The foregoing example illustrates something further. Let be a two-length wavelet. Then the only other two-length wavelet with the same autocorrelation as that of **Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \left\{b_0{,\ }b_{1}\right\} }**
is its (midpoint) reverse . (Of course, when we speak about the wavelet and its reverse **Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \left(b^{*}_{{ 1}}{,\ }b^{*}_0\right) }**
, we mean the equivalence class of **Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle \left\{cb_O{,\ }cb_{1}\right\} }**
and for any constant *c* where . Here, we take *c* = 1.) Thus, the only other two-length wavelet with the same autocorrelation as that of (2, 1) is (1, 2). It turns out that there are wavelets of greater length with the same autocorrelation, but there is no other wavelet of length 2 with the same autocorrelation. Figure 18 shows the wavelet (2, 1) and a wavelet of greater length, both of which have the same autocorrelation. The wavelet of greater length is obtained by convolving the wave-let (2, 1) with the all-pass filter shown in the figure.

A causal all-pass filter is defined as a causal filter with an amplitude spectrum equal to unity for all frequencies. The transfer function of a causal all-pass filter is

**(**)

where *A*(*Z*) is the *Z*-transform of a finite-length minimum-delay wavelet *a*, and is the *Z*-transform of the reverse . The amplitude spectrum of this all-pass filter is unity because both *a* and its reverse have the same amplitude spectrum. This reverse is, of course, maximum delay. Because a causal all-pass wavelet has a unit amplitude spectrum, it follows that its inverse is equal to its reverse; that is, , or

**(**)

Thus the inverse is noncausal and one-sided to the past.

Generally, a causal all-pass wavelet has infinite length. However, in two special cases, a causal all-pass wavelet reduces to a finite-length wavelet. In the first special case, *A*(*Z*) is a constant. As a result, the transfer function of the causal all-pass wavelet is equal to one, so the causal all-pass wavelet is the unit spike {1, 0, 0, 0, …}, where the 1 occurs at time index 0. Such a causal all-pass wavelet is called trivial because the unit spike convolved with any signal does not change the signal. In the second special case, where *n* is a positive integer. In this case, , so the causal all-pass wavelet produces a pure delay of *n* time units.

## Continue reading

Previous section | Next section |
---|---|

Energy | Canonical representation |

Previous chapter | Next chapter |

Frequency | Synthetics |

## Also in this chapter

- Wavelets
- Fourier transform
- Z-transform
- Delay: Minimum, mixed, and maximum
- Two-length wavelets
- Illustrations of spectra
- Delay in general
- Energy
- Canonical representation
- Zero-phase wavelets
- Symmetric wavelets
- Ricker wavelet
- Appendix G: Exercises