A note on the equivalence of three major propagator algorithms for computational stability and efficiency
A note on the equivalence of three major propagator algorithms for computational stability and efficiency
复制标题
关于计算稳定性和效率的三种主要传播器算法的等效性的说明
DOI：
--
复制
发表时间：
2012
期刊：
影响因子：
--
通讯作者：
Hui
中科院分区：
文献类型：
--
作者：
Yanlu Ma;Rongjiang Wang;Hui
本文指出，标准正交化方法、子矩阵方法和递归反射-透射矩阵方法是密切相关的，它们同样很好地解决了原始的B-H传播矩阵方法中的数值不稳定性。本文提出了一种基于正交化和Langer块对角分解的稳定而有效的方法来计算水平分层模型对平面谱波的响应。它是一种数值稳健的用于高频率、大层厚和水平慢度的B-H矩阵方法。该技术适用于计算反射-透射系数，体波接收器功能和瑞利波色散。
It is shown in this note that the three methods, the orthonormalization method, the minor matrix method and the recursive reflection-transmission matrix method are closely related and solve the numerical instability in the original Thomson-Haskell propagator matrix method equally well. Another stable and efficient method based on the orthonormalization and the Langer block-diagonal decomposition is presented to calculate the response of a horizontal stratified model to a plane, spectral wave. It is a numerically robust Thomson-Haskell matrix method for high frequencies, large layer thicknesses and horizontal slownesses. The technique is applied to calculate reflection-transmission coefficients, body wave receiver functions and Rayleigh wave dispersion.