Based on the platform of Matlab and the theory of digital signal processing, we propose a method in the cepstrum domain for dynamic load spectra identification of machinery. We demonstrate that the dynamic load spectra can be identified from the response signal of the system, based on cepstra. An ARMA model is built based on the harmonic retrieval by high-order spectra. The coefficients of a Green function are determined and the window width can be estimated. Finally the effectiveness of the method is validated by simulation results.