I have non-uniformly sampled data and I'm trying to look at its Power Spectral Density (PSD). For that, I looked at the Periodogram PSD estimate of an interpolated (with uniform sampling period) version of my data and compared it to the Lomb-Scargle PSD estimate of my original data.

enter image description here

enter image description here

I'm surprised to see that the resulting PSDs are different, especially for high frequencies, given that the interpolated data seems to follow quite well the original data. Now, I'm not sure which one should I trust!

I would appreciate it if anyone can explain why this difference and which one to use.

P.S: I found this related question For non-uniformly sampled data, how to decide whether interpolation or Lomb-Scargle periodogram is better? but no answers/comments were given.

  • $\begingroup$ How did you interpolate? Is the underlying signal bandlimited? $\endgroup$
    – Gillespie
    Jan 27 at 5:09
  • $\begingroup$ @Gillespie, interpolated using 'interp1' matlab function with the linear method. I don't have a clear idea on the spectrum of the underlying signal, that's exactly why I'm trying to look at the PSD. $\endgroup$
    – Likely
    Jan 27 at 7:04
  • $\begingroup$ @Gillespie Also, this particular data I included is supposed to contain only noise. $\endgroup$
    – Likely
    Jan 27 at 7:14
  • $\begingroup$ Given that it's not bandlimited, there's no telling how aliased the signal was in the sampling process. And interpolation of nonuniformly sampled signals is complicated even when they are bandlimited: dsp.stackexchange.com/a/84232/55647. So I would say it's hard (impossible?) to tell which PSD is more accurate. $\endgroup$
    – Gillespie
    Jan 27 at 13:13
  • $\begingroup$ @Gillespie, Does this mean that we can't say, for sure, that Lomb-Scargle PSD gives a more accurate PSD, although it's supposed to be adequate for dealing with nonuniform sampling. $\endgroup$
    – Likely
    Jan 27 at 13:24


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.