DFT, 1 Sinusoid and related functions, Range and Resolution.

Cuthbert Nyack
In this applet a section of a sinusoid from T1 to T2 is sampled and the DTFT and DFT calculated. An important property of the DFT is the relation between the spectrum and the beginning and end of the sample interval of a periodic signal. The variation of the spectra as T1, T2 are changed is illustrated by the applet. The real and imaginary parts and the magnitude of the DFT are also shown.

eg parameters:-
(3.0, 0.5, 4.21, 6.2, 20, 0, 15, 10.0, 15.0, 0.5) in this case one period is sampled and the DFT contains 1 line at the sinusoid frequency. The sampling frequency is ~ 60rad/s so there are 20 samples per period and the line in the DFT spectrum corresponds to m = 1.
(3.0, 0.5, 2.11, 6.19, 40, 0, 15, 10.0, 15.0, 0.25) in this case two periods are sampled and the DFT contains 1 line at the sinusoid frequency. The sampling frequency is ~ 60rad/s so there are 40 samples in 2 periods and the line in the DFT spectrum corresponds to m = 2.
(3.0, 0.5, 2.11, 10.27, 40, 0, 15, 10.0, 20.0, 0.25) in this case four periods are sampled and the DFT contains 1 line at the sinusoid frequency. The sampling frequency is ~ 30rad/s so there are 40 samples in 4 periods and the line in the DFT spectrum corresponds to m = 4.

(6.0, 0.5, 0.01, 8.49, 10, 0, 10, 10.0, 15.0, 1.0) in this case the Nyquist rate is violated since the sampling frequency(6.668rad/s) is less than twice the sinusoid frequency(12rad/s). The DFT spectrum between 0 and ws/2 contains 1 line at 0.668rad/s (ie ws - w). In this case the samples appear to trace out a sinusoid whose period is equal to 9 periods of the sinusoid at 6rad/s ie with a frequency of 6/9 = 0.667rad/s.

When an integral number of periods are sampled the DFT frequencies outside of the main peak falls on the zeros of the DTFT. When (n + ½) periods are sampled, the DFT frequencies fall near the maxima of the DTFT and the spectrum contains contains several lines. The effect is illustrated by the eg parameters below:-
(3.0, 0.5, 2.11, 5.12, 30, 0, 20, 10.0, 15.0, 0.5) in this case 1.5 periods are sampled and the DFT spectrum contains several lines. The sampling frequency is ~ 60.53rad/s, the DFT line separation is 2.017rad/s so the sinusoid frequency of 3rad/s corresponds to m ~ 1.5. Whenever the sinusoid frequency does not correspond to one of the DFT frequencies, then the spectrum is spread out among all the values of m. In this case the 2 largest frequency components correspond to m = 1 and m = 2.
(3.0, 0.5, 0.01, 11.32, 60, 0, 12, 10.0, 30.0, 0.2) in this case 5.5 periods are sampled and the DFT spectrum is a PSK spectrum rather than a sinusoid spectrum.

When n periods are sampled starting from the zero of the sine, the spectrum is imaginary. For (n + ½) periods, it is real and is an almost mixture for (n + ¼) and (n + ¾) periods. When Sampling starts from the peak of the sinusoid, this is equivalent to sampling a cosine and the opposite applies. eg (3.0, 0.5, 2.61, 6.65, 25, 0, 16, 10.0, 15.0, 0.4). For (n + ½) periods, it is not completely imaginary.

Range and Resolution.
When using the DFT it is necessary to know how to control the range and resolution of the spectrum. The range is determined by the sampling time and is equal to ws/2. The resolution is determined by the length of the sampling interval ~ 2p/NTs = ws/N .
eg parameters (3.0, 1.0, 0.01, 62.2, 100, 0, 100, 225.0, 65.5, 0.1)
The sampling frequency is 10.0rad/s (2p(Ns - 1)/(T2 - T1) and the range is from 0 to 5.0rad/s. The Resolution is shown on the applet as the Line Sep = 0.1rad/s(ws/Ns = 10.0/100). There is only 1 line in the spectrum of significance with magnitude ~50 which corresponds to a sinusoid of amplitude 50 * (2/100) = 1. The frequency of the spectral line is 3.0rad/s.

Reducing the sampling interval while keeping the number of samples constant increases both the range and the resolution.
eg parameters (3.0, 1.0, 0.01, 24.89, 100, 0, 100, 90.0, 65.5, 0.1)
The sampling frequency is 25.0rad/s, the range is 12.5rad/s and the resolution is 0.25rad/s.

eg parameters (3.0, 1.0, 0.01, 30.79, 50, 0, 50, 225.0, 65.5, 0.2)
show the range is 5.0rad/s and the resolution is 0.20rad/s.

Bad choice of parameters can result in a range which is less than the frequency of interest.
eg parameters (3.0, 1.0, 0.01, 63.32, 40, 0, 40, 168.0, 65.5, 0.2)
show the range is 1.935rad/s and the resolution is 0.096rad/s. The frequency shown in the range between 0 and ws/2 is at the aliased frequency 3.87 - 3.0 ~ 0.87rad/s. The correct frequency of 3rad/s occurs in the range ws/2 to ws





The applet also allows the DFT spectra of other functions to be examined by varying the Wfm parameter.

Some of the options available are:-
Wfm = 1, a half wave rectified sinusoid.
Wfm = 2, a full wave rectified sinusoid.
Wfm = 3, a saturated sinusoid, saturation level set by a.
Wfm = 4, a half wave saturated sinusoid, saturation level set by a.
Wfm = 5, a full wave rectified saturated sinusoid, saturation level set by a.
Wfm = 6, an SCR controlled waveform, set by a.
Wfm = 7, a Triac controlled waveform, set by a.
Wfm = 8, a parabolic approximation to a sinusoid.
Wfm = 9, a cubic approximation to a sinusoid.

Return to main page
Return to page index
COPYRIGHT 2007 Cuthbert Nyack.