R: MTM spectrum

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` | 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 \& 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 3.7-4 Index]