A linear programming approach to the estimation of the power spectra of harmonic processes