Standard algorithms for estimating the power spectral density of finite discrete data require interpolating missing samples, which usually produces biased estimates. An unbiased estimate can be obtained by taking the Fourier transform of the unbiased estimator of the circular autocorrelation, using only the available data. With missing samples, this estimator can produce negative power spectral densities, but converges to positive values when averaged over a sufficient number of realizations.