ABSTRACT This study presents a hybrid numerical technique for solving multi‐dimensional time‐fractional wave equations, combining Laplace transforms and a Chebyshev‐node‐based Lagrange pseudo‐spectral method (LPPSM). The method is designed to overcome the computational burdens of time‐stepping schemes for fractional derivatives. By transforming the temporal derivatives into the frequency domain, it enhances numerical stability and avoids time‐stepping constraints. The spatial discretization via LPPSM ensures spectral convergence. The time‐domain solution is recovered using an optimized Talbot contour inversion. The efficacy of the proposed method is demonstrated on benchmark problems, including the fractional Klein‐Gordon equation, achieving exceptional accuracy (errors on the order of ). Crucially, to bridge the gap with practical applications, the method is applied to a time‐fractional telegraph equation modeling signal propagation in a transmission line with physical parameters . The results confirm the method's high computational efficiency, superior accuracy, and potential for real‐time analysis of fractional‐order dynamics in high‐frequency electronic systems.
Kamran et al. (Sat,) studied this question.