<p>The problem you are trying to solve is one generally known as spectral estimation. In addition to the practical references provided by jstarek, I suggest you also read some fundamentals. </p>

<p>The book by S.M.Kay <a href="https://rads.stackoverflow.com/amzn/click/0130151599" rel="nofollow noreferrer">Modern Spectral Estimation</a>, although almost 20 years old, is totally up to date and is a very well explained (easy to read) introduction. </p>

<p>You will see that the folding algorithm is just a flavour of the Welch periodogram, and that your cross-correlation approach fits as an MLE estimation scheme (nowadays also known as signal subspace techniques).</p>

CC BY-SA 3.0
2018-01-17T15:56:54.683Z
2018-01-17T15:56:54.683Z
Ela
<p>To understand the basics of algorithms like that, I find it's often best to look at what amateur practicioners are writing (as opposed to diving directly into scientific papers). The description of the <a href="http://neutronstar.joataman.net/technical/epoch_folding.html" rel="nofollow noreferrer">folding algorithm</a> at joataman.net is thorough and pragmatic.</p>

<p>A great in-depth book on the topic is Lorimer/Kramer: <a href="http://www.jb.man.ac.uk/pulsarhandbook/" rel="nofollow noreferrer">Handbook of Pulsar Astronomy</a>, Oxford University Press.</p>

<p>For the actual work, I would suggest the <a href="http://sigproc.sourceforge.net/" rel="nofollow noreferrer">sigproc</a> package by Duncan Lorimer. Its home is at Sourceforge, but there are forks all over the place because the original software is a bit quirky when run on modern Linuxes.</p>

CC BY-SA 3.0
2017-12-18T13:27:30.600Z
2017-12-18T13:27:30.600Z
jstarek
<p>I've a data file with 2 columns, each containing time sampled voltage output (a sequence of real numbers) from a radio telescope, so there is a pair of telescope which are measuring the signal here. The source of the data is a pulsar and I've to find its period and dispersion measure.</p>

<p>I'm an amateur trying to understand what is really happening in the method. The data is time sampled at the Nyquist sampling frequency (2B) and the signal is measured at central frequency f and broadband B.</p>

<p>I studied the cross correlation between the two voltages V1 and V2 from essential radio astronomy, but didn't quite understand the method. Do we have to evaluate the cross correlation function of each reading and then plot it with time to get the pulse profile? Or do we have to take fourier transform of each pair of voltage(how?) and then plot it to get the pulse profile? We find period of a pulsar from its average pulse profile, right?</p>

<p>It is also written that to find pulsar timing, we've to fold data from many pulses modulo instantaneous pulse period (how can we do that if we don't know the pulse period in the first place?). This sentence is everywhere and I'm not able to understand it each time.</p>

<p>I'm confused about the appropriate series of steps to be done here to get the pulsar period and dispersion measure. Could you give me some insight?</p>

<p>Thanks.</p>

CC BY-SA 4.0
2020-01-03T20:47:36.670Z
if you solve your problem can you share detail with me i am also facing same problem
user28958
<radio-astronomy><data-analysis><observational-astronomy><pulsar><radio-telescope>
How to find period of a pulsar from a given time sampled voltage data?
apooravc
# How to find period of a pulsar from a given time sampled voltage data?

(Asked by: username_0 on 2017-05-26)

I've a data file with 2 columns, each containing time sampled voltage output (a sequence of real numbers) from a radio telescope, so there is a pair of telescope which are measuring the signal here. The source of the data is a pulsar and I've to find its period and dispersion measure.

I'm an amateur trying to understand what is really happening in the method. The data is time sampled at the Nyquist sampling frequency (2B) and the signal is measured at central frequency f and broadband B.

I studied the cross correlation between the two voltages V1 and V2 from essential radio astronomy, but didn't quite understand the method. Do we have to evaluate the cross correlation function of each reading and then plot it with time to get the pulse profile? Or do we have to take fourier transform of each pair of voltage(how?) and then plot it to get the pulse profile? We find period of a pulsar from its average pulse profile, right?

It is also written that to find pulsar timing, we've to fold data from many pulses modulo instantaneous pulse period (how can we do that if we don't know the pulse period in the first place?). This sentence is everywhere and I'm not able to understand it each time.

I'm confused about the appropriate series of steps to be done here to get the pulsar period and dispersion measure. Could you give me some insight?

Thanks.

**Tags:** radio-astronomy, data-analysis, observational-astronomy, pulsar, radio-telescope

> **Comments:**
> **username_1**: if you solve your problem can you share detail with me i am also facing same problem  

---

## 2 answers

### Answer by username_2 on 2018-01-17. Score: 2

The problem you are trying to solve is one generally known as spectral estimation. In addition to the practical references provided by username_3, I suggest you also read some fundamentals.

The book by S.M.Kay [Modern Spectral Estimation](https://rads.stackoverflow.com/amzn/click/0130151599), although almost 20 years old, is totally up to date and is a very well explained (easy to read) introduction.

You will see that the folding algorithm is just a flavour of the Welch periodogram, and that your cross-correlation approach fits as an MLE estimation scheme (nowadays also known as signal subspace techniques).


---

### Answer by username_3 on 2017-12-18. Score: 2

To understand the basics of algorithms like that, I find it's often best to look at what amateur practicioners are writing (as opposed to diving directly into scientific papers). The description of the [folding algorithm](http://neutronstar.joataman.net/technical/epoch_folding.html) at joataman.net is thorough and pragmatic.

A great in-depth book on the topic is Lorimer/Kramer: [Handbook of Pulsar Astronomy](http://www.jb.man.ac.uk/pulsarhandbook/), Oxford University Press.

For the actual work, I would suggest the [sigproc](http://sigproc.sourceforge.net/) package by Duncan Lorimer. Its home is at Sourceforge, but there are forks all over the place because the original software is a bit quirky when run on modern Linuxes.


---

