next up previous print clean
Next: Prediction-error filters on synthetic Up: Time-series analysis Previous: Random numbers into the

THE ERROR FILTER FAMILY

A simple regression for a prediction filter (f1,f2) is

 
 \begin{displaymath}
\left[ 
\begin{array}
{c}
 x_2 \  
 x_3 \  
 x_4 \  
 x_5...
 ...\; \left[ 
\begin{array}
{c}
 f_1 \  
 f_2 \end{array} \right]\end{displaymath} (9)
Notice that each row in this equation says that xt fits a linear combination of x at earlier times; hence the description of f as a ``prediction" filter. The error in the prediction is simply the left side minus the right side. Rearranging the terms, we get  
 \begin{displaymath}
\left[ 
\begin{array}
{c}
 0 \  
 0 \  
 0 \  
 0 \end{ar...
 ...
\begin{array}
{c}
 1 \  
 -f_1 \  
 -f_2 \end{array} \right]\end{displaymath} (10)
We have already written programs for regressions like (9). Regressions like (10), however, often arise directly in practice. They are easier to solve directly than by transforming them to resemble (9).

Multiple reflections are predictable. It is the unpredictable part of a signal, the prediction residual, that contains the primary information. The output of the filter (1,-f1, -f2) is the unpredictable part of the input. This filter is a simple example of a ``prediction-error" (PE) filter. It is one member of a family of filters called ``error filters."

The error-filter family are filters with one coefficient constrained to be unity and various other coefficients constrained to be zero. Otherwise, the filter coefficients are chosen to have minimum power output. Names for various error filters follow:


$(1, a_1,a_2,a_3, \cdots ,a_n)$ prediction-error (PE) filter  
$(1, 0, 0, a_3,a_4,\cdots ,a_n) $		gapped PE filter with a gap of 2 
$(a_{-m}, \cdots, a_{-2}, a_{-1}, 1, a_1, a_2, a_3, \cdots ,a_n)$ interpolation-error (IE) filter 
$(a_{-m}, \cdots, a_{-4}, a_{-3}, 0,0,1,0,0, a_3, a_4, \cdots ,a_n)$ a gapped IE filter

A program for computing all the error filters will be presented after we examine a collection of examples. 1


 
next up previous print clean
Next: Prediction-error filters on synthetic Up: Time-series analysis Previous: Random numbers into the
Stanford Exploration Project
10/21/1998