Prev: How to Plot this Kind of Matrix
Next: Cross-correlations - how to read
From: Pat van Heerden on 2 Aug 2010 09:36
Hi

I have a data-set (variable star light-curve) that is sampled non-uniformly (avg sampling 0.548Hz) - I re-sample the data by using 'interp1' over a new time-vector specified within the same begin/end values of original (sampling 1Hz). Total time is 4.8hours.

I am looking for the rotation period of the star - which is buried in the light-curve data. What I want to do is to filter the data to remove all frequencies outside the specified region wherein the rotation period lies. (Same kind of method as 'idealfilter')

Frequencies of interest between 0.02 - 0.08 Hz. I played around with 'tstool' to do the timing analysis, but I have a bit of an issue.

I want to do datafolding or fftstacking to get rid off all the background noise - in the hope that the fundamental frequencies that I am looking for will arise, but to do that I need to work with only the relevant frequency range.

My method that I am trying to do is the following.

1.Re-sample the data to 1Hz.
2.Filter the data to relevant frequencies.
3.Break the data into 1hour segmant's
4.Get 'fft' of each segment.
5.Stack the fft's.
6.Plot the results over the frequency range.

Hope for the best.

My current problem is with section 2. (Or am I going about this in the wrong way???)

Hope someone can give me a starting point.

Thanx
Pat
 | 
Pages: 1
Prev: How to Plot this Kind of Matrix
Next: Cross-correlations - how to read