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
harmonic signal processing technology to solve the problem of
processing undersampled
harmonic signals.