Scribd
Upload
85 views

The Excel FFT Function v1.1

The document describes using Excel's FFT function to analyze time series data from a decaying torsion oscillator. It provides 2048 time series data points that are input to the FFT function. This returns 2048 complex Fourier coefficients corresponding to frequencies from 0-25Hz in increments of 0.25Hz. The power at each frequency is calculated from the coefficients and shows an exponential decay characteristic, with the maximum power at 0.4883Hz. The normalization used ensures the total Fourier power equals the mean squared amplitude of the original time series data.

Uploaded by

Mailson1
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
85 views

The Excel FFT Function v1.1

The document describes using Excel's FFT function to analyze time series data from a decaying torsion oscillator. It provides 2048 time series data points that are input to the FFT function. This returns 2048 complex Fourier coefficients corresponding to frequencies from 0-25Hz in increments of 0.25Hz. The power at each frequency is calculated from the coefficients and shows an exponential decay characteristic, with the maximum power at 0.4883Hz. The normalization used ensures the total Fourier power equals the mean squared amplitude of the original time series data.

Uploaded by

Mailson1
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 5

The Excel FFT Function v1.

1
P. T. Debevec
February 12, 2006
The discrete Fourier transform may be used to identify periodic structures in time
series data. Suppose that a physical process is represented by the function of time, h ( t ) .
The function is sampled at N times, tk = k t where k = 0,1, 2,..., N 1 . From these N
measurements, h k , N complex amplitudes, H n , are determined which satisfy the N
equations
N 1

Hn = h k e

ik

2 n
N

k =0

The sampled function then has the discrete Fourier expansion

hk=

N 1

1
N

Hn e

i n

2 k
N

n =0

This equation can be cast in familiar form with 2 k N = ( 2 To ) k (To N ) = o k t = o tk


hk=

N 1

1
N

H
n=0

e i notk .

The right-hand side is the discrete analogue to the complex form of the Fourier expansion

h(t ) =

ce
n =

in o t

where the complex coefficients, cn , are given by

1
cn =
To

To

0 h(t ) e

in o t

dt .

The Excel data analysis package has a Fourier analysis routine which calculates
the complex coefficients, H n , from the time series data, h k . The routine requires that
the number of samples in the time series data be a power of 2, i.e. N = 2m . The example
in this note uses N = 2048 . The Excel function is not well documented, but it is
straightforward to use. This note describes the Excel worksheet, Fourier_example.xls,
which is in the Physics 401 web site under Tutorials and Lectures, Experiment 10.

The time series data in this example are obtained from sampling a function
describing the free decay of a torsion oscillator for time t > to ,

( t ) = Ae a(t t ) sin ( 2 f ( t to ) ) .
o

The function is calculated in time steps of 0.020 s, which corresponds to sampling rate of
50.0 Hz. The time series data are shown in the Fig. 1 below.
fre e de cay dampe d oscillator
6000

raw counts

4000
2000
0
-2000
-4000
-6000
0

10

15

20

25

30

35

40

time (s)

Fig. 1 Plot of time series data decay of torsion oscillator


The data occupy cells B3 to B2050 in the data worksheet of the workbook. Click on
Tools in the Excel menu bar, and select Data Analysis. In Data Analysis select Fourier
Analysis, and a simple dialog box appears. The dialog box is shown below in Fig. 2.

Fig. 2 Fourier Analysis dialog box

For Input Range enter $B$3:$B$2050, the location of the time series data, and for output
range enter a convenient place on the worksheet, for example, $J$3:$J$2050. After
selecting the OK button, Excel returns in column J the complex coefficients, 2048 of
them. A portion of column J is shown in Fig. 3 below.

Fig. 3 Portion of Excel worksheet showing FFT output


Excel insists in displaying 15 digits for both the real and imaginary part of the coefficient,
and thus the output occupies consider space. Examination of the column of complex
numbers shows that the first entry (in cell J3) and the 1025th entry (in cell J1027) are real.
These cells represent the coefficients for the zero frequency and the folding frequency,
f c . Since the data were sampled at 50.0 Hz, the folding frequency is 50 Hz = 25 Hz.
Cells J3 through J1027 represent the complex coefficients for N 2 + 1 frequencies from 0
Hz to 25 Hz. The frequency increment is then f = f c N . It is convenient to make a
column of frequencies next to the column of coefficients. This column includes cells I3
to I1027. The Fourier Analysis routine operates only on the time series data. The
correspondence between frequency and complex coefficient must be calculated
independently. The 1023 entries from cell J1028 to cell J2050 are identical to the 1023
entries from cell J1026 to J4. These entries contain no additional information, and they
could have been omitted from the display. The Fourier Analysis of N time series data
yields the complex coefficients for only N 2 + 1 frequencies.
The power in each frequency bin is proportional to the square of the modulus of
the complex coefficient. The constant of proportionality adopted in this note has the
property that the mean squared amplitude of the time series data is equal to the total
power in the frequency domain as shown below. Excel has a number of functions for
complex arguments. Multiplication of two complex numbers is accomplished with the
function IMPRODUCT(z1, z2). Multiplication of the coefficients in column J by 1 N
produces the normalized coefficients in column N. The absolute value of a complex
number is accomplished with the function IMABS(z). Applying the IMABS function to
the normalized coefficients in column N and multiplying by a factor of 2 produces the
magnitude of the Fourier coefficients in column P. The magnitude of the Fourier
coefficients as a function of frequency is shown in Fig. 4 below. The maximum value
occurs in the frequency bin at 0.4883 Hz.

magnitude of Fourier amplitude


versus frequency
magnitude of Fourier amplitude

500
400
300
200
100
0
0

10

15

20

25

frequency (Hz)

Fig. 4 Magnitude of normalized Fourier coefficient versus frequency


Finally, the power at each frequency (except zero frequency) is the square of the
magnitude. The power is calculated in column Q and shown in Fig. 5 below. The power
at zero frequency is the square of the magnitude divided by 2. This factor of 2 also
appears in non-discrete Fourier analysis. This power distribution is characteristic of
exponential decay.
Fourier power
versus frequency
1.0E+05

Fourier power

1.0E+03

1.0E+01

1.0E-01

1.0E-03
0

10

15

frequency (Hz)

Fig. 5 Fourier power versus frequency

20

25

The normalization adopted in this note has a convenient property. Define the
mean squared amplitude of the time series data as
mean squred amplitude =

1
N

N 1

( tk ) .
k =0

The square of the amplitude is calculated in column C. The mean squared amplitude is
calculated in cell F19. For these data it has the value of 1105488.754. A pure sine wave
with an amplitude of 2 , for which an integral number of periods occur in the total
sampled time, has a mean squared amplitude of unity. Thus the mean squared amplitude
has the same role as the rms value of a periodic function. The total Fourier power, the
sum of the entries in column Q, is calculated in cell G19. For these data it also has the
value of 1105488.754, because of the normalization. If only the relative power in
different frequency bins is needed, then the normalization is not required. The relative
power is given by the square of the modulus of the Excel output directly.

You might also like