L. Braile, Earth and Atmospheric Sciences

Purdue University         braile@purdue.edu

 

 

 

 

Return to L. Braile’s Home Page

 

EAS 553 Home Page:

 

http://web.ics.purdue.edu/~braile/eas553/eas553home.htm

 

Prof. Lawrence W. Braile

Department of Earth and Atmospheric Sciences

Purdue University

West Lafayette, IN 47907-1397

braile@purdue.edu

 

Last updated: March 3, 2003

 

Matlab Codes (click on any of the following files to download; save in your Matlab work folder; see EAS 553 notes for additional information, program listings and examples):

 

Refraction:

 

refrin.m  Refraction (head wave) inversion for plane, homogeneous dipping layers using slope-intercept method.

 

refrdat.m  Sample data set for refrin.m.

 

txcurv.m  Calculate refraction and reflection times for n plane, homogeneous, dipping layers.

 

rfltim.m  Reflection times – needed by txcurv.

 

rfrtim.m  Refraction times – needed by txcurv.

 

angle.m  Snell’s law angle calculations – needed by txcurv.

 

Reflection:

 

txrefl.m  Reflection travel times for a layered model – exact times and one-layer approximations using Vave and Vrms.

 

rfltim.m  Reflection travel times (exact) using Snell’s law ray trace through layered model.

 

t2x2.m  T-squared, X-squared program.

 

txdat.m  Sample data for t2x2 program.

 

linfit.m  Least squares straight line fit program for t2x2.m.

 

txline.m  Line drawing function for t2x2.m.

 

seis1d.m  1-D synthetic seismograms.

 

ricker.m  Ricker wavelet calculation; needed for seis1d.m.

 

signal.m  Signal wavelet calculation; needed for seis1d.m.

 

seis1ddataStandardModel.m  Sample data for seis1d.m; rename as seis1ddata.m to use.

 

seis1ddataGradientModel.m  Sample data for seis1d.m; rename as seis1ddata.m to use.

 

EAS 553 Data Source Links:

 

 

 

Return to Braile’s Home Page