problem with fft periodogram

Travis E. Oliphant oliphant.travis at ieee.org
Wed Oct 25 06:31:55 EDT 2006


aitor69gonzalez at gmail.com wrote:
> Hello,
> 
> I am ploting a fft periodogram of my data using a script (found in
> internet: http://linuxgazette.net/115/andreasen.html ) that gave me
> good results before. There should be a periodicity of >160 in the data
> as can be seen by eye. However, this script now says that there is
> periodicity of ~9. Can somebody explain me what's wrong? Thank you in
> advance.

The units are wrong on your period axis.  Right now you have them in 
units of "sample-spacing".  So, the plot is telling you that you have a 
periodicity of ~9 sample spacings.   To get it in unites of minutes you 
need to multiply period by the difference in minutes

period_in_minutes = period * (minutes[1] - minutes[0])

Then, plot period_in_minutes versus power.  I see a peak around 180 
minutes in your data.



-Travis




More information about the Python-list mailing list