RSEIS (version 1.0-9)

mtapspec: 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}

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

Details

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

References

Lees, J. M. and Park, J., 1995: Multiple-taper spectral analysis: A stand-alone C-subroutine, {Computers & Geology, 21(2), 199-236. } [object Object]

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

misc