Pisarenko harmonic decomposition

Pisarenko harmonic decomposition, also referred to as Pisarenko's method, is a method of frequency estimation.^[1] This method assumes that a signal, ${\displaystyle x(n)}$, consists of ${\displaystyle p}$ complex exponentials in the presence of white noise. Because the number of complex exponentials must be known a priori, it is somewhat limited in its usefulness.

Pisarenko's method also assumes that ${\displaystyle p+1}$ values of the ${\displaystyle M\times M}$ autocorrelation matrix are either known or estimated. Hence, given the ${\displaystyle (p+1)\times (p+1)}$ autocorrelation matrix, the dimension of the noise subspace is equal to one and is spanned by the eigenvector corresponding to the minimum eigenvalue. This eigenvector is orthogonal to each of the signal vectors.^[2]

The frequency estimates may be determined by setting the frequencies equal to the angles of the roots of the polynomial

${\displaystyle V_{\rm {min}}(z)=\sum _{k=0}^{p}v_{\rm {min}}(k)z^{-k}}$

or the location of the peaks in the frequency estimation function (or the pseudo-spectrum)

${\displaystyle {\hat {P}}_{\rm {PHD}}(e^{j\omega })={\frac {1}{|\mathbf {e} ^{H}\mathbf {v} _{\rm {min}}|^{2}}}}$,

where ${\displaystyle \mathbf {v} _{\rm {min}}}$ is the noise eigenvector and

${\displaystyle e={\begin{bmatrix}1&e^{j\omega }&e^{j2\omega }&\cdots &e^{j(M-1)\omega }\end{bmatrix}}^{T}}$.