       Re: power spectrum

• To: mathgroup at smc.vnet.net
• Subject: [mg52742] Re: [mg52695] power spectrum
• From: Sseziwa Mukasa <mukasa at jeol.com>
• Date: Sat, 11 Dec 2004 05:22:05 -0500 (EST)
• References: <200412100123.UAA18915@smc.vnet.net>
• Sender: owner-wri-mathgroup at wolfram.com

```On Dec 9, 2004, at 8:23 PM, George Szpiro wrote:

> If I have a string of data - let us say a time series of 200
> observations -
> how can I find the power spectrum? I did Fourier[list] but have no
> idea what
> the output means.
>
> george at netvision.net.il)
>

If your data is real valued

ListPlot[Abs[Fourier[list]][[Range[Quotient[Length[list],2]+1]]],PlotJoi
ned->True]

will give you the power spectrum from frequencies 0 to half your sample
rate.

For complex valued observations

ListPlot[Abs[RotateRight[Fourier[list],Quotient[Length[list],2]]],PlotJo
ined->True]

will give you the power spectrum from -1/2*sample rate to 1/2*sample
rate - 1/(Length[list] * sample rate).

As for what it all means, section 13.4 of Numerical Recipes
(http://lib-www.lanl.gov/numerical/index.html) explains power spectra.

Regards,

Ssezi

```

• References:
• Prev by Date: Re: Reduce[] in version 5.1
• Next by Date: Re: Linear Programming
• Previous by thread: power spectrum
• Next by thread: Re: power spectrum