mtapspec {RSEIS}R Documentation

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


[Package RSEIS version 2.1-4 Index]