In this paper, an efficient Transmission Line Matrix (TLM) approach based on the shift operator (SO) has been developed to model electromagnetic wave interactions with gyroelectric media. The main idea of this technique is to formulate the electric current density vector components by introducing the equivalence between time differential operator ∂ ∂t and discrete time shift operator z. A concise formulation of voltage sources modeling the frequency dispersive properties of gyroelectric media is then deduced and implemented. Numerical simulations illustrate the Faraday rotation phenomenon in time domain, and in the frequency domain, reflection and transmission coefficients of left hand circular polarization and right-hand circular polarization waves are also calculated. A comparison of SO-TLM scheme with five other approaches according to the criteria of accuracy and CPU time is presented. Numerical experiments show that SO-TLM provides the most accurate and fastest results.