As a non-linear time variant system, the four-quadrant line converter of an electric multiple unit (EMU) was expressed by linear time periodic functions near an operating point and modeled by a steady-state harmonic domain matrix. The components were then combined according to the circuit connection and relations of the feedback control loops to form a complete converter model. The proposed modeling method allows the study of the amplitude of harmonic impedances to explore harmonic coupling. Moreover, the proposed method helps provide a better design for the converter controllers, as well as solves the problem in coordination operation between the EMUs and the AC supply. On-site data from an actual $CRH_2$ high-speed train were used to validate the modeling principles presented in the paper.