PLOT.MATN: plot a matrix of several seismograms

In RSEIS: Seismic Time Series Analysis Tools

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

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<[e-mail]>

See Also

swig, matsquiggle, dowiggles, varsquiggle

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

RSEIS documentation built on Aug. 19, 2023, 5:07 p.m.