02.10.2019 Views

UploadFile_6417

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

154 Chapter 5 THE DISCRETE FOURIER TRANSFORM<br />

5.3 THE DISCRETE FOURIER TRANSFORM<br />

The discrete Fourier series provides a mechanism for numerically computing<br />

the discrete-time Fourier transform. It also alerted us to a potential<br />

problem of aliasing in the time domain. Mathematics dictates that the<br />

sampling of the discrete-time Fourier transform result in a periodic sequence<br />

˜x(n). But most of the signals in practice are not periodic. They<br />

are likely to be of finite duration. How can we develop a numerically computable<br />

Fourier representation for such signals? Theoretically, we can take<br />

care of this problem by defining a periodic signal whose primary shape is<br />

that of the finite-duration signal and then using the DFS on this periodic<br />

signal. Practically, we define a new transform called the discrete Fourier<br />

transform (DFT), which is the primary period of the DFS. This DFT<br />

is the ultimate numerically computable Fourier transform for arbitrary<br />

finite-duration sequences.<br />

First we define a finite-duration sequence x(n) that has N samples<br />

over 0 ≤ n ≤ N − 1asanN-point sequence. Let ˜x(n) beaperiodic signal<br />

of period N, created using the N-point sequence x(n); that is, from (5.19)<br />

˜x(n) =<br />

∞∑<br />

r=−∞<br />

x(n − rN)<br />

This is a somewhat cumbersome representation. Using the modulo-N<br />

operation on the argument we can simplify it to<br />

˜x(n) =x(n mod N) (5.21)<br />

A simple way to interpret this operation is the following: if the argument<br />

n is between 0 and N − 1, then leave it as it is; otherwise add or subtract<br />

multiples of N from n until the result is between 0 and N − 1. Note<br />

carefully that (5.21) is valid only if the length of x(n)isN or less. Furthermore,<br />

we use the following convenient notation to denote the modulo-N<br />

operation.<br />

x((n)) N<br />

△<br />

= x(n mod N) (5.22)<br />

Then the compact relationships between x(n) and ˜x(n) are<br />

˜x(n) =x((n)) N<br />

x(n) =˜x(n)R N (n)<br />

(Periodic extension)<br />

(Window operation)<br />

(5.23)<br />

The rem(n,N) function in MATLAB determines the remainder after dividing<br />

n by N. This function can be used to implement our modulo-N<br />

Copyright 2010 Cengage Learning. All Rights Reserved. May not be copied, scanned, or duplicated, in whole or in part. Due to electronic rights, some third party content may be suppressed from the eBook and/or eChapter(s).<br />

Editorial review has deemed that any suppressed content does not materially affect the overall learning experience. Cengage Learning reserves the right to remove additional content at any time if subsequent rights restrictions require it.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!