mtapspec: MTM spectrum

In RSEIS: Seismic Time Series Analysis Tools

MTM spectrum

Description

Multi-tape Method Spectrum

Usage

mtapspec(a, dt, klen = length(a), MTP = NULL)

Arguments

a	vector time series
dt	sample rate
klen	length of fft
MTP	MTM parameters, list: kind kind of taper average nwin number of windows npi number of Pi-prolate functions inorm normalization flag

Details

MTP represent parameters that control the multi-tape pi-prolate functions used by mtapspec. See reference for details.

Value

LIST

dat	input data
dt	sample rate
spec	Estimated power spectrum
dof	degrees of freedom for each frequency
Fv	F-values for each frequency
Rspec	real part of complex spectrum
Ispec	imaginary part of complex spectrum
freq	frequencies
df	delta frequency
numfreqs	number of frequencies
klen	length used in fft
mtm	input MTM parameters, see above

Author(s)

Jonathan M. Lees<jonathan.lees.edu>

References

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

See Also

fft

Examples

data(CE1)
a <- list(y=CE1$y[CE1$x>5.443754 & CE1$x<5.615951], dt=CE1$dt)

    Mspec <-  mtapspec(a$y,a$dt, klen=4096,
             MTP=list(kind=2,nwin=5, npi=3,inorm=0)  )   

RSEIS documentation built on Sept. 13, 2024, 1:09 a.m.