The simultaneous estimation of a time series harmonic components and covariance structure