Open AccessYet another elastic plane‐wave, layer‐matrix algorithm
Author(s) -
Chapman C. H.
Publication year - 2003
Publication title -
geophysical journal international
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.302
H-Index - 168
eISSN - 1365-246X
pISSN - 0956-540X
DOI - 10.1046/j.1365-246x.2003.01958.x
Subject(s) - isotropy , diagonal , simple (philosophy) , symmetry (geometry) , plane (geometry) , matrix (chemical analysis) , algorithm , diagonal matrix , plane symmetry , stability (learning theory) , mathematical analysis , geometry , mathematics , computer science , physics , optics , materials science , composite material , philosophy , epistemology , machine learning
SUMMARY Many techniques have been developed to calculate the elastodynamic response of a stack of plane layers to a plane, spectral wave. These variations on the original propagator matrix method, commonly called the Haskell matrix method, aim to improve the efficiency and numerical stability of the algorithm. In this note, a simple variant, based on the Langer block‐diagonal decomposition and second‐order minors, is presented. The method is valid in isotropic media and anisotropic media with up–down symmetry on vertical planes of symmetry. The method is efficient and extremely simple to program—code in Matlab* is presented for the isotropic case.