R: MTMdisp

MTMdisp {RSEIS} | R Documentation

## MTMdisp

### Description

Display MTM displacement spectrum.

### Usage

    MTMdisp(a, f1 = f1, f2 = f2, len2 = 1024, PLOT = FALSE)

### Arguments

`a` | seismic velocity trace, as a ts structure (list(y=trace, dt=sample rate)
`f1` | low frequency
`f2` | high frequency
`len2` | length of fft
`PLOT` | logical, TRUE=plot

### Details

Uses Multi-taper estimate of spectrum and divides the spectrum by 1/(2*pi*f) to get integration of velocity seismogram.

### Value

Returns displacement spectrum. Graphical Side effect.

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

mtapspec

### Examples

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

    len2 <- next2(length(xvel$y))
    Spec <- MTMdisp(xvel, f1=.01, f2=25, len2=len2, PLOT=FALSE )

[Package RSEIS version 3.7-4 Index]