-# Copyright (C) 2008-2010 W. Trevor King <wking@drexel.edu>
+# Copyright (C) 2010-2012 W. Trevor King <wking@tremily.us>
#
# This file is part of Hooke.
#
-# Hooke is free software: you can redistribute it and/or
-# modify it under the terms of the GNU Lesser General Public
-# License as published by the Free Software Foundation, either
-# version 3 of the License, or (at your option) any later version.
+# Hooke is free software: you can redistribute it and/or modify it under the
+# terms of the GNU Lesser General Public License as published by the Free
+# Software Foundation, either version 3 of the License, or (at your option) any
+# later version.
#
-# Hooke is distributed in the hope that it will be useful,
-# but WITHOUT ANY WARRANTY; without even the implied warranty of
-# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
-# GNU Lesser General Public License for more details.
+# Hooke is distributed in the hope that it will be useful, but WITHOUT ANY
+# WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
+# A PARTICULAR PURPOSE. See the GNU Lesser General Public License for more
+# details.
#
-# You should have received a copy of the GNU Lesser General Public
-# License along with Hooke. If not, see
-# <http://www.gnu.org/licenses/>.
+# You should have received a copy of the GNU Lesser General Public License
+# along with Hooke. If not, see <http://www.gnu.org/licenses/>.
"""Wrap :mod:`numpy.fft` to produce 1D unitary transforms and power spectra.
Notes
-----
-
If the units on your data are Volts,
and your sampling frequency is in Hz,
then `freq_axis` will be in Hz,
def unitary_avg_power_spectrum(data, freq=1.0, chunk_size=2048,
overlap=True, window=window_hann):
- """Compute the unitary avgerage power spectrum of `data`.
+ """Compute the unitary average power spectrum of `data`.
See Also
--------
Notes
-----
-
The expected return values are [#dft]_:
.. math:: X_k = \sum_{m=0}^{n-1} x_m \exp^{-2\pi imk/n}
Notes
-----
-
Which is:
.. math:: \sum_{m=0}^{n-1} |x_m|^2 dt = \sum_{k=0}^{n-1} |X_k|^2 df
Notes
-----
-
Analytic result:
.. math:: \rfft(\rect(at)) = 1/|a|\cdot\sinc(f/a)
Notes
-----
-
Analytic result:
.. math::