next up previous print clean
Next: Multidimensional deconvolution breakthrough Up: FILTERING ON A HELIX Previous: FILTERING ON A HELIX

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 $A(Z) = 1 + a_1Z + a_2Z^2 + \cdots$.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  
 \begin{displaymath}
y_k \quad=\quad x_k + \sum_{i=1}^{N_a} a_i x_{k-i}\end{displaymath} (1)
Convolution computes yk from xk whereas deconvolution (also called back substitution) does the reverse. Rearranging (1) we get  
 \begin{displaymath}
x_k \quad=\quad y_k - \sum_{i=1}^{N_a} a_i x_{k-i}\end{displaymath} (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 up previous print clean
Next: Multidimensional deconvolution breakthrough Up: FILTERING ON A HELIX Previous: FILTERING ON A HELIX
Stanford Exploration Project
4/27/2004