RSEIS (version 2.4-5)

PLOT.MATN: plot a matrix of several seismograms

Description

Matrix of several seismograms

Usage

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

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)

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

See Also

PICK.GEN

Examples

dt = 0.001

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

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

#   % sample rate = 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)