OCTAVE YULEWALKER

This program computes the ALL POLE spectrum for a signal. In addition to the octave plots, a GNUPLOT (plotspec.gp) file is output along with a data file yw.dat which is readable by the plotspec.gp file. The code is interactive and one uses the mouse to pick the order of the process on the autocorrelation.

 The yulewalker.m program is run in octave.
You can run it on BSEGY, data or
you can run it on the autocorrelation of
BSEGY data.

In either case, the spectrum will be an
all pole spectrum.  You will be prompted
to use your mouse to select the maximum
lag in the autocorrelation.

If you run BXCR followed by BSTK, you can
compute an average autocorrelation for the
entire shot gather, and run yulewalker.m on
that.

Other *.m files required (bsegin.m and segyinfo.m)
must be in the same directory before yulewalker.m
is run in an octave session.

Figure 7: OCTAVE YULE WALKER: Octave program which computes the ALL POLE spectrum. Input can be either a seismic trace data or an autocorrelation of trace data (either must be in BSEGY format, *.seg file). Compare to Figure 6 FFT plot. See BXCR 12.0.16 for how to create an autocorrelation as input.