Discrete Fourier Transform Example

or

The Transform of a Sine Wave and its Frequency Spectrum
Showing how Kaiser-Bessel Windowing Reduces Ringing
from the Jump Caused by the Periodic Continuation of the
Endpoints

Sine wave with known frequency.

Number of points in the DFT.

[Graphics:Images/index_gr_1.gif]
[Graphics:Images/index_gr_2.gif]

Sampling frequency.

[Graphics:Images/index_gr_3.gif]
[Graphics:Images/index_gr_4.gif]

Mapping from index to frequency.

[Graphics:Images/index_gr_5.gif]

Maximal frequency at index = n/2.

[Graphics:Images/index_gr_6.gif]
[Graphics:Images/index_gr_7.gif]

Frequency of sine wave in Hz.

[Graphics:Images/index_gr_8.gif]
[Graphics:Images/index_gr_9.gif]

Sine wave points from 0 to n-1.

[Graphics:Images/index_gr_10.gif]

Import values from a file and readjust the length and sampling frequency.

[Graphics:Images/index_gr_11.gif]
[Graphics:Images/index_gr_12.gif]
[Graphics:Images/index_gr_13.gif]
[Graphics:Images/index_gr_14.gif]
[Graphics:Images/index_gr_15.gif]

Plot the discretized sine wave.

[Graphics:Images/index_gr_16.gif]

[Graphics:Images/index_gr_17.gif]

[Graphics:Images/index_gr_18.gif]

Kaiser-Bessel Window.

[Graphics:Images/index_gr_19.gif]

[Graphics:Images/index_gr_20.gif]
[Graphics:Images/index_gr_21.gif]
[Graphics:Images/index_gr_22.gif]
[Graphics:Images/index_gr_23.gif]
[Graphics:Images/index_gr_24.gif]
[Graphics:Images/index_gr_25.gif]

[Graphics:Images/index_gr_26.gif]

[Graphics:Images/index_gr_27.gif]

Fourier Transform.

Compute a Fourier transform [Graphics:Images/index_gr_28.gif][Graphics:Images/index_gr_29.gif]and plot its power spectrum.

[Graphics:Images/index_gr_30.gif]
[Graphics:Images/index_gr_31.gif]
[Graphics:Images/index_gr_32.gif]

[Graphics:Images/index_gr_33.gif]

[Graphics:Images/index_gr_34.gif]
[Graphics:Images/index_gr_35.gif]
[Graphics:Images/index_gr_36.gif]
[Graphics:Images/index_gr_37.gif]
[Graphics:Images/index_gr_38.gif]
[Graphics:Images/index_gr_39.gif]
[Graphics:Images/index_gr_40.gif]

Windowed Fourier Transform.

Window the time domain first.

[Graphics:Images/index_gr_41.gif]
[Graphics:Images/index_gr_42.gif]

[Graphics:Images/index_gr_43.gif]

[Graphics:Images/index_gr_44.gif]

Now compute the power spectrum.  Note it is much cleaner.

[Graphics:Images/index_gr_45.gif]
[Graphics:Images/index_gr_46.gif]
[Graphics:Images/index_gr_47.gif]

[Graphics:Images/index_gr_48.gif]

[Graphics:Images/index_gr_49.gif]


Converted by Mathematica      January 28, 2007