Next: INTERPOLATION ERROR
Up: CALCULATING ERROR FILTERS
Previous: CALCULATING ERROR FILTERS
Theory for
stabilizing least squares,
using equations
() and
(),
was described earlier in this book.
I installed this stabilization, along with the filter determinations
discussed in this chapter,
but as I expected,
stabilization in this highly overdetermined application
showed no advantages.
Nevertheless, it is worth seeing how stabilization is implemented,
particularly since the changes to the program
calling iner()
make for more informative plots.
The input data is modified by appending a zeropadded impulse
at the data's end.
The output will contain the filter impulse response in that region.
The spike size is chosen to be compatible with the data size,
for the convenience of the plotting programs.
The weighting function in the appended region
is scaled according to how much stabilization is desired.
Figure 12 shows the complete input and residual.
wz33
Figure 12
Data from the North Sea
(extracted from Yilmaz and Cumro dataset 33)
processed by prediction error.
Rightmost box is weighted according to the desired stabilization.
The truncation event is weighted by zero.
It also illustrates the problem that output data flows beyond the length
of the input data because of the nonzero length of the filter.
This extra output is undoubtedly affected by the truncation of
the data, and its energy should not be part of the energy minimization.
Therefore it is weighted by zero.
EXERCISES:

Given a sinusoidal function ,a threeterm recurrence relationship
predicts x_{t} from the previous two points,
namely, x_{t} = a_{1} x_{t1} + a_{2} x_{t2}.
Find a_{1} and a_{2} in terms of .HINT: See chapter .
(Notice that the coefficients depend on but not .)

Figure 9 has a separate filter for each trace.
Consider the problem of finding a single filter for all the traces.
What is the basic operator and its adjoint?
Assemble these operators using subroutine contrunc() .

Examine the filters on Figure 12.
Notice that, besides the pulse at the water depth,
another weak pulse occurs at double that depth.
Suggest a physical mechanism.
Suggest a mechanism relating to computational approximations.
Next: INTERPOLATION ERROR
Up: CALCULATING ERROR FILTERS
Previous: CALCULATING ERROR FILTERS
Stanford Exploration Project
10/21/1998