The pseudo-excitation method combined with the integral transform method (PEM-ITM) is presented to investigate the ground vibration of a coupled track-soil system induced by moving random loads. Commonly in the track model, the rail, sleepers, rail pads, and ballast are modelled as an infinite Euler beam, discretely distributed masses, discretely distributed vertical springs, and a viscoelastic layer, respectively. The soil is regarded as a homogenous isotropic half-space coupled with the track using the boundary condition at the surface of the ground. By introducing a pseudo-excitation, the random vibration analysis of the coupled system is converted into a harmonic analysis. The analytical form of evolutionary power spectral density responses of the simplified coupled track-soil system under a random moving load is derived in the frequency/wavenumber domain by PEM-ITM. In the numerical examples, the effects of different parameters, such as the moving speed, the soil properties, and the coherence of moving loads, on the ground response are investigated.