This invention provides a spectral analysis method for undersampled harmonic signals, comprising the following steps: acquiring multi-channel leaf tip timing signals; constructing multiple snapshot matrices based on the multi-channel leaf tip timing signals; calculating a cross-correlation matrix; calculating the following power expression of the cross-correlation matrix; calculating the steering vector; and obtaining the function beamforming output. The technical solution of this invention can be widely applied in the field of harmonicsignalprocessing technology to solve the problem of processing undersampled harmonic signals.