doGABOR.AR: Gabor Transform with AR spectrum method

Description Usage Arguments Details Value Note Author(s) References See Also Examples

View source: R/doGABOR.AR.R

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.

See Also

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

Examples

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)

RSEIS documentation built on Aug. 2, 2021, 5:07 p.m.