plot a matrix of several seismograms

Share:

Description

Matrix of several seismograms

Usage

1
2
3
PLOT.MATN(ascd, tim=1, dt=1,  WIN=c(0,1), labs="",
notes=notes, sfact=1,ampboost=0,  shift=NULL, LOG="",
COL='red', add=1, AXES=1, units=NULL, VS=FALSE)

Arguments

ascd

N by K matrix of seismograms where

tim

time values fo x-axis

dt

sample interval, seconds

WIN

vector, time window for zoom

labs

vector of labels for each panel

notes

vector of notes for each panel

sfact

scaling factor, 1=window, 2=trace

ampboost

increase each amplitude by this multiplier

shift

vector, shift each trace by these time

LOG

log x-axis

COL

vector of colors or indexes to colors

add

numeric, to existing plot. add = 1,2,3 if add=1 plot and add traces, add =2 plot, but no traces, add = 3 no plot, but add traces. DEFAULT=1

AXES

numeric, 0,1,2,3,4; default=1

units

label for units of Y-axis

VS

var-squiggle display

Details

Plots a matrix of seismograms that each have the same starting time. For the AXES argument, 0 = no axes, AXES=1 plot scale for largest amplitude band and a multiplier for all others, AXES=2 left side, AXES=3 right side, AXES=4 alternate sides

Value

Graphical side effects and,

n

n

windiv

matrix of n rows, with columns=(window Y min, window Y max, user Y min, user Y max)

Author(s)

Jonathan M. Lees<jonathan.lees@unc.edu>

See Also

swig

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
dt  <-  0.001

 t  <-  seq(0, 6, by=0.001)

thefreqs  <-  seq(from=10, to=100, by=10)
theamps  <-   runif(length(thefreqs))

#   sample rate is 1000 Hz, 0.001 seconds 601 samples
 x  <-  NULL

for(i in 1:length(thefreqs))
{
x  <-  cbind(x, theamps[i]*sin(2*pi*thefreqs[i]*t))


}

PLOT.MATN(x,  dt = dt)

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.