This study investigates the propagation of shear horizontal transverse waves in a functionally graded piezoelectric half-space (FGPHS), where the material properties vary linearly and quadratically. The analysis focuses on deriving and understanding the dispersion characteristics of such waves in in-homogeneous media. The WKB approximation method is employed to obtain the dispersion relation analytically, considering the smooth variation of material properties. To validate and study the wave behavior numerically, two advanced techniques were utilized: the Semi-Analytical Finite Element with Perfectly Matched Layer (SAFE-PML) and the Semi-Analytical Infinite Element (SAIFE) method incorporating a (1/r) decay model to simulate infinite media. The numerical implementation uses the Rayleigh–Ritz method to discretize the wave equation, and Gauss 3-point quadrature is applied for efficient numerical integration. The dispersion curves are plotted to illustrate the wave behavior in the graded piezoelectric medium. The results from SAFE-PML and SAIFE are in excellent agreement, indicating that these techniques effectively model the shear horizontal transverse wave propagation in such structures. This study also demonstrates that combining finite and infinite element approaches provides accurate and reliable simulation of wave phenomena in functionally graded piezoelectric materials, which has applications in sensors, actuators, and non-destructive testing.