In this paper, we propose a new random vibration analysis model for the vertical vehicle-track coupling system, which combines the symplectic method and the spectral element method (SEM) to enhance calculation efficiency and accuracy. The model assumed that the vehicle is composed of a multi-rigid body with 10 DOFs, and its motion equation is derived using the traditional frequency analysis method. Tracks are considered periodic structures, and the part between sleeper support points is considered a substructure. The motion equation of the periodic track-substructure is established by spectral element method (SEM), and the propagation of vibration wave between the substructures is calculated by the symplectic method. The frequency responses of the vehicle-track coupling system are calculated by using the German low-interference irregularity spectrum and short-wave irregularity spectrum as excitation. Then, the rationality of this method is verified by comparing it with the results from the traditional symplectic method and simulations using Universal Mechanism software. The comparison shows that this method surpasses the traditional symplectic method in high-frequency accuracy while maintaining high computational efficiency. Although the finite element method also has the capability for high-frequency analysis, this method can calculate the high-frequency vibration of an infinitely long track using just one element, making it more suitable for the efficient analysis of vehicle-track coupling systems. Finally, the model is applied to analyze the distribution characteristics of the vibration response of the coupled system in the frequency domain.