Eigenvalue decomposition of a parahermitian matrix: extraction of analytic eigenvalues

An analytic parahermitian matrix admits an eigenvalue decomposition (EVD) with analytic eigenvalues and eigenvectors except in the case of multiplexed data. In this paper, we propose an iterative algorithm for the estimation of the analytic eigenvalues. Since these are generally transcendental, we find a polynomial approximation with a defined error. Our approach operates in the discrete Fourier transform (DFT) domain and for every DFT length generates a maximally smooth association through EVDs evaluated in DFT bins; an outer loop iteratively grows the DFT order and is shown, in general, to converge to the analytic eigenvalues. In simulations, we compare our results to existing approaches.