[time-nuts] Calculate spectral content from a series of zero crossing time stamps?

[Tijd Dingen]

Consider the following scenario. We have a signal source of about 10 kHz,
with unknown phase noise. Let's for simplicity's sake assume for now that
the phase noise is large enough that it will be detectable by the following
approach.

We measure every zero crossing with lets say 1 ns accuracy. So we have a
signal with a nominal period of 100 us, and we can measure every zero crossing
to within 1 ns. This gives you ~ 10,000 data points every second.

Now how does one efficiently calculate the spectral content based on these
10,0000 zero crossings? The end result would be the spectral density, centered
around that nominal 10 kHz frequency.