mtapspec {RSEIS}

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)  )