Abstract: Based on the characteristics of the complex permittivity of three kinds of general dispersive media models,i.e. Debye model, the Lorentz model and the Drude model, may be described by rational polynomial function of j? ., theshift operator finite difference time domain method is developed to deal with the hybrid dispersive media model. Thenumerical results are in good agreement with the analysis results and show the feasibility of the presented method.