doGABOR.AR: Gabor Transform with AR spectrum method In RSEIS: Seismic Time Series Analysis Tools

Description

Gabor Transform with AR spectrum method

Usage

 ```1 2``` ```doGABOR.AR(Xamp, DT = 0.008, multi = 1, scale.def = 0, TWIN = 2, TSKIP = 0.2, PCTTAP = 0.05, pord=100, PLOT=TRUE) ```

Arguments

 `Xamp` signal `DT` sample rate interval (s) `multi` Multiples of time window estimate `scale.def` scaling flag for plotting (0=raw, 1=log, 2=sqrt) `TWIN` time for window `TSKIP` time for skip `PCTTAP` percent of taper to apply to individual windows `pord` order for the AR process (default=100) `PLOT` logical, TRUE=plot to device

Details

This is a spectrogram function similar to the Gabor Transform but uses the AR method for spectrum estimation.

Value

list

 `sig` input signal `dt` deltat `numfreqs` Number of frequencies output `wpars` input parameters list(Nfft=numfreqs, Ns=Ns, Nov=Nov, fl=fl, fh=fh) `DSPEC` spectrum image `HIMAT` matrix with high values of F-test at 90 percent confidence `freqs` output frequencies (y axis) `tims` output times (x-axis)

Note

The main difference between this and other similar calls is the way the windows are determined.

Author(s)

Jonathan M. Lees<jonathan.lees@unc.edu>

References

Lees, J. M. and Park, J., 1995: Multiple-taper spectral analysis: A stand-alone C-subroutine, Computers \& Geology, 21(2), 199-236.

Percival, Donald B.,Walden, Andrew T. (1993):Spectral Analysis for Physical Applications,Cambridge University Press, Cambridge, 583p.

 ```1 2 3 4 5 6 7 8 9``` ```data(KH) ### swig(KH) Xamp <- KH\$JSTR[[1]] Xamp <- Xamp[57914:72989] EV <- doGABOR.AR(Xamp, DT = KH\$dt[1] , multi = 1, scale.def = 0, TWIN = 2, TSKIP = 0.2, PCTTAP = 0.05) ```