Power spectrum calculation and fit
- To: mathgroup at smc.vnet.net
- Subject: [mg12704] Power spectrum calculation and fit
- From: Randy Cutter <cutter at vims.edu>
- Date: Wed, 3 Jun 1998 02:20:56 -0400
- Sender: owner-wri-mathgroup at wolfram.com
If z is a data series, a 1024 record list, will the following code
correctly produce a power spectrum and the least squares fit slope to
the log transformed power spectrum versus the log of the sequence?
Thanks for any suggestions.
ftofz1 = Abs[Fourier[z]];
ftofz1 = Take[%, 512];
ftzsq = %^2;
lnftzsq = Log[%];
Fit[lnftzsq, {1, x}, Log[x]]
---------------------
Randy Cutter
Virginia Institute of Marine Science Route 1208
Gloucester Point, VA 23062 USA
phone 804 684-7365
fax 804 684-7399 or 7045
email: cutter at vims.edu
---------------------