OCTAVE cafwd3.m - Perform Kelvin-Voigt forward computation given stiffness and damping values
cafwd3 prompts for:
Choose to
compute with observed data for comparison or not (GUI Click
Box)
IF WITHOUT DATA
User enters stiffness, C1 (m^2/s^2), and damping, C2 (m^2/s)
(GUI Entry Box)
LOOPS for additional parameter choices and computes plots
of phase velocity (m/s), decay (1/m), and Q (unitless) vs
frequency.
ELSE
Assuming bvas.his and bamp.his are available (from bvas and
bamp runs)
Prompts for 2 mouse clicks selecting a frequency range (low
to high frequency)
This permits an estimate for initial C1 and C2 from the
data. (GUI Entry Box)
LOOPS for additional parameter choices and computes plots
of phase velocity (m/s), decay (1/m), and Q (unitless) vs
frequency.
END
Basic Seismic Utilities (BSU) program run in Octave or Matlab. Within an octave session, type cafwd3 and a GUI pops up as shown above. One either computes the forward problem, displayed against results from a down-hole survey (requires two files, bvas.his and bamp.his), or just computes the forward problem with no data comparisons. This computation is in the context of dispersive body waves (transmission measurement) under the Kelvin-Voigt representation.
Requires two
files from execution of bvas and bamp programs ONLY if
"with data" selected.
bvas.his
bamp.his
Plots of phase
velocity, decay, and Q, as a function of frequency under
Kelvin-Voigt representation
of transmission of body waves.
bamp(1) bvas(1) caplot(1) caplot3.m(7) cainv3.m(7)
No known bugs.
Copyright © 2017 by Paul Michaels
This program is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2 of the License, or (at your option) any later version. This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with this program; if not, write to the Free Software Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
P. Michaels, PE. <pm@cgiss.boisestate.edu>