This paper presents the recursive asymptotic hybrid matrix method for acoustic waves in multilayered piezoelectric media. The hybrid matrix method preserves the numerical stability and accuracy across large and small thicknesses. For discussion and comparison, the scattering matrix method is also presented in physics-based form and coherent form. The latter form resembles closely that of hybrid matrix method and helps to highlight their relationship and distinction. For both scattering and hybrid matrix methods, their formulations in terms of eigenwaves solution are provided concisely. Making use of the hybrid matrix, the recursive asymptotic method without eigenwaves solution is described and discussed. The method bypasses the intricacies of eigenvalue-eigenvector approach and requires only elementary matrix operations along with thin- layer asymptotic approximation. It can be used to determine Green’s function matrix readily and facilitates the trade-off between computation efficiency and accuracy.