Next: EXAMPLES OF SIMPLE 2-D Up: Claerbout: Recursion via the Previous: CONVOLUTION ON A HELIX

# BASIC REVIEW OF 1-D RECURSIVE FILTERS

Convolution is the operation we do on polynomial coefficients when we multiply polynomials. Deconvolution is likewise for polynomial division. Often these ideas are described as polynomials in the variable Z. Take X(Z) to denote the polynomial whose coefficients are samples of input data, and let A(Z) likewise denote the filter. The convention I adopt here is that the first coefficient of the filter has the value +1, so the filter's polynomial is .To see how to convolve, we now identify the coefficient of Zk in the product Y(Z)=A(Z)X(Z). The usual case (k larger than the number Na of filter coefficients) is
 (1)
Deconvolution is a complicated process in seismology that includes estimating A(Z). More simply, let us take the filter A(Z) as known, the output Y(Z) as known, and we want to go back to find the input X(Z)=Y(Z)/A(Z). Again we simply identify the coefficient of each power of Zk in Y(Z)=A(Z)X(Z), but now we seek to recursively find xk instead of yk. Rearranging (1) we get
 (2)
where now we are finding the output xk from its past outputs xk-i and from the present input yk. We see that the deconvolution process is essentially the same as the convolution process, except that the filter coefficients are used with opposite polarity; and they are applied to the past outputs instead of the past inputs. That is why deconvolution must be done sequentially while convolution can be done in parallel.

Next: EXAMPLES OF SIMPLE 2-D Up: Claerbout: Recursion via the Previous: CONVOLUTION ON A HELIX
Stanford Exploration Project
10/14/1997