Recover the spectrum of covariance matrix: a non-asymptotic iterative method