get.corner {RSEIS} | R Documentation |
Search for low frequency asymptote, corner frequency, and fall off slope of seismic spectrum.
get.corner(INfreq, INspec, dt, f1, f2, PLOT = FALSE, VERBOSE = FALSE)
INfreq |
frequency vector |
INspec |
spectrum |
dt |
deltaT |
f1 |
low frequency for modeling, Hz |
f2 |
High frequency for modeling, Hz |
PLOT |
logical, TRUE=plot |
VERBOSE |
TRUE=diagnostics |
This routine does not assume any particular mathematical model. It searches for a three parameters that describe two lines that mimic the displacement spectrum. The search is done via least squares.
Model of 3 parameters, best fit.
Jonathan M. Lees<jonathan.lees.edu>
brune.doom
data(CE1) #### set frequency range for modeling f1 = 0.01 f2 = 14.0 ## set up data and parameters amp = CE1$y len2 = 2*next2(length(amp)) a = list(y=amp, dt=CE1$dt) Spec =MTMdisp(a, f1=f1, f2=f2, len2=len2, PLOT=FALSE ) lspec = Spec$displ ### get initial estimate of parameters xc = get.corner( Spec$f , lspec, dt, f1, f2, PLOT=FALSE)