FILT.spread {RSEIS}

Filter trace with a spread of filters

Description

Show a time series and a spread of user defined filters to show signal at a variety of bandwidths.

Usage

FILT.spread(x, y, dt, fl = fl, fh = fh, sfact = 1,
 WIN = NULL, PLOT = TRUE, TIT = NULL, TAPER = 0.05,
 POSTTAPER=0.05, RM=FALSE, zp=TRUE )

Arguments

x

x-axis

y

y-amplitude

dt

delta-t, sec

fl

vector of low frequency cut offs

fh

vector of high frequency cut offs

sfact

scale factor, 0,1

WIN

xlimits to constrain plotting

PLOT

logical, plotting

TIT

title

TAPER

taper data prior to filter, percent cosine, default=NULL

POSTTAPER

taper output after filter, percent cosine, default=0.05

RM

Remove mean value from trace, default=FALSE

zp

zero phase filter, default=TRUE

Details

Use the TAPER and POSTTAPER to reduce the edge effects prior to and after filtering.

Value

list:

FMAT

matrix of time series filtered

Notes

Notes for filter of each element of FMAT

Author(s)

Jonathan M. Lees<[e-mail]>

See Also

butfilt, PLOT.MATN

Examples

data(KH)
dt  <-  KH$dt[1]

y  <-   KH$JSTR[[1]]

x  <-   seq(from=0, by=dt, length=length(y))



fl <- rep(1/100, 5)
fh <- 1/c(1,2,5,10,20)

FILT.spread(x, y, dt, fl = fl, fh = fh, sfact = 1,
   WIN = NULL, PLOT = TRUE, TIT = NULL, TAPER = 0.05)