doGABOR.AR {RSEIS}R Documentation

Gabor Transform with AR spectrum method

Description

Gabor Transform with AR spectrum method

Usage

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.

See Also

evolfft, evolMTM, MTM.drive, GETARAIC, doGABOR.AR, DOsgram, doGABOR.MTM

Examples




data(KH)
###   PICK.GEN(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)







[Package RSEIS version 2.3-3 Index]