Review DTFT DTFT Properties Examples Summary Lecture 9: Discrete-Time Fourier Transform Mark Hasegawa-Johnson ECE 401: Signal and Image Analysis, Fall 2020 Review DTFT DTFT Properties Examples Summary 1 Review: Frequency Response 2 Discrete Time Fourier Transform 3 Properties of the DTFT 4 Examples 5 Summary Review DTFT DTFT Properties Examples Summary Outline 1 Review: Frequency Response 2 Discrete Time Fourier Transform 3 Properties of the DTFT 4 Examples 5 Summary Review DTFT DTFT Properties Examples Summary What is Signal Processing, Really? When we process a signal, usually, we're trying to enhance the meaningful part, and reduce the noise. Spectrum helps us to understand which part is meaningful, and which part is noise. Convolution (a.k.a. filtering) is the tool we use to perform the enhancement. Frequency Response of a filter tells us exactly which frequencies it will enhance, and which it will reduce. Review DTFT DTFT Properties Examples Summary Review: Convolution A convolution is exactly the same thing as a weighted local average. We give it a special name, because we will use it very often. It's defined as: X X y[n] = g[m]f [n − m] = g[n − m]f [m] m m We use the symbol ∗ to mean \convolution:" X X y[n] = g[n] ∗ f [n] = g[m]f [n − m] = g[n − m]f [m] m m Review DTFT DTFT Properties Examples Summary Review: DFT & Fourier Series Any periodic signal with a period of N samples, x[n + N] = x[n], can be written as a weighted sum of pure tones, N−1 1 X x[n] = X [k]ej2πkn=N ; N k=0 which is a special case of the spectrum for periodic signals: 2π radians 1 cycles N seconds samples !0 = ; F0 = ; T0 = ; N = ; N sample T0 second Fs cycle cycle and N−1 X X [k] = x[n]e−j2πkn=N : n=0 Review DTFT DTFT Properties Examples Summary Tones in ! Tones out Suppose I have a periodic input signal, N−1 1 X x[n] = X [k]ej2πkn=N ; N k=0 and I filter it, y[n] = h[n] ∗ x[n]; Then the output is a sum of pure tones, at the same frequencies as the input, but with different magnitudes and phases: N−1 1 X y[n] = Y [k]ej2πkn=N : N k=0 Review DTFT DTFT Properties Examples Summary Frequency Response Suppose we compute y[n] = x[n] ∗ h[n], where N−1 1 X x[n] = X [k]ej2πkn=N ; and N k=0 N−1 1 X y[n] = Y [k]ej2πkn=N : N k=0 The relationship between Y [k] and X [k] is given by the frequency response: Y [k] = H(k!0)X [k] where 1 X H(!) = h[n]e−j!n n=−∞ Review DTFT DTFT Properties Examples Summary Outline 1 Review: Frequency Response 2 Discrete Time Fourier Transform 3 Properties of the DTFT 4 Examples 5 Summary Review DTFT DTFT Properties Examples Summary Aperiodic An \aperiodic signal" is a signal that is not periodic. Periodic acoustic signals usually have a perceptible pitch frequency; aperiodic signals sound like wind noise, or clicks. Music: strings, woodwinds, and brass are periodic, drums and rain sticks are aperiodic. Speech: vowels and nasals are periodic, plosives and fricatives are aperiodic. Images: stripes are periodic, clouds are aperiodic. Bioelectricity: heartbeat is periodic, muscle contractions are aperiodic. Review DTFT DTFT Properties Examples Summary Periodic The spectrum of a periodic signal is given by its Fourier series, or equivalently in discrete time, by its discrete Fourier transform: N−1 1 X j 2πkn x[n] = X [k]e N N k=0 N−1 X −j 2πkn X [k] = x[n]e N n=0 Review DTFT DTFT Properties Examples Summary Aperiodic The spectrum of an aperiodic signal we will now define to be exactly the same as that of a periodic signal except that, since it never repeats itself, its period has to be N = 1: N−1 1 X j 2πkn x[n] ≈ lim X [k]e N N!1 N k=0 N−1 X −j 2πkn X [k] ≈ lim x[n]e N N!1 n=0 Review DTFT DTFT Properties Examples Summary An Aperiodic Signal is like a Periodic Signal with Period= 1 Review DTFT DTFT Properties Examples Summary Aperiodic The spectrum of an aperiodic signal we will now define to be exactly the same as that of a periodic signal except that, since it never repeats itself, its period has to be N = 1: N−1 1 X j 2πkn x[n] ≈ lim X [k]e N N!1 N k=0 N−1 X −j 2πkn X [k] ≈ lim x[n]e N N!1 n=0 2πk But what does that mean? For example, what is N ? Let's try this definition: allow k ! 1, and force ! to remain constant, where 2πk ! = N Review DTFT DTFT Properties Examples Summary Aperiodic Let's start with this one: N−1 1 X j 2πkn x[n] ≈ lim X [k]e N N!1 N k=0 Imagine this as adding up a bunch of tall, thin rectangles, each 2π with a height of X [k], and a width of d! = N . In the limit, as N ! 1, that becomes an integral: N−1 1 X 2π j 2πkn x[n] ≈ lim X [k]e N N!1 2π N k=0 1 Z 2π = X (!)ej!nd!; 2π !=0 where we've used X (!) = X [k] just because, as k ! 1, it makes more sense to talk about X (!). Review DTFT DTFT Properties Examples Summary Approximating the Integral as a Sum Review DTFT DTFT Properties Examples Summary Periodic Now, let's go back to periodic signals. Notice that ej2π = 1, and j 2πk(n+N) j 2πk(n−N) j 2πkn for that reason, e N = e N = e N . So in the DFT, we get exactly the same result by summing over any complete period of the signal: N−1 X −j 2πkn X [k] = x[n]e N n=0 N X −j 2πkn = x[n]e N n=1 N−4 X −j 2πkn = x[n]e N n=−3 N−1 2 X −j 2πkn = x[n]e N (N−1) n=− 2 X −j 2πkn = x[n]e N n=any complete period Review DTFT DTFT Properties Examples Summary Aperiodic Let's use this version, because it has a well-defined limit as N ! 1: N−1 2 X −j 2πkn X [k] = x[n]e N (N−1) n=− 2 The limit is: N−1 X2 X (!) = lim x[n]e−j!n N!1 (N−1) n=− 2 1 X = x[n]e−j!n n=−∞ Review DTFT DTFT Properties Examples Summary Discrete Time Fourier Transform (DTFT) So in the limit as N ! 1, 1 Z π x[n] = X (!)ej!nd! 2π −π 1 X X (!) = x[n]e−j!n n=−∞ X (!) is called the discrete time Fourier transform (DTFT) of the aperiodic signal x[n]. Review DTFT DTFT Properties Examples Summary Outline 1 Review: Frequency Response 2 Discrete Time Fourier Transform 3 Properties of the DTFT 4 Examples 5 Summary Review DTFT DTFT Properties Examples Summary Properties of the DTFT In order to better understand the DTFT, let's discuss these properties: 0 Periodicity 1 Linearity 2 Time Shift 3 Frequency Shift 4 Filtering is Convolution Property #4 is actually the reason why we invented the DTFT in the first place. Before we discuss it, though, let's talk about the others. Review DTFT DTFT Properties Examples Summary 0. Periodicity The DTFT is periodic with a period of 2π. That's just because ej2π = 1: X X (!) = x[n]e−j!n n X X X (! + 2π) = x[n]e−j(!+2π)n = x[n]e−j!n = X (!) n n X X X (! − 2π) = x[n]e−j(!−2π)n = x[n]e−j!n = X (!) n n In fact, we've already used this fact. I defined the inverse DTFT in two different ways: 1 Z π 1 Z 2π x[n] = X (!)ej!nd! = X (!)ej!nd! 2π −π 2π 0 Those two integrals are equal because X (! + 2π) = X (!). Review DTFT DTFT Properties Examples Summary 1. Linearity The DTFT is linear: z[n] = ax[n] + by[n] $ Z(!) = aX (!) + bY (!) Proof: X Z(!) = z[n]e−j!n n X X = a x[n]e−j!n + b y[n]e−j!n n n = aX (!) + bY (!) Review DTFT DTFT Properties Examples Summary 2. Time Shift Property Shifting in time is the same as multiplying by a complex exponential in frequency: −j!n0 z[n] = x[n − n0] $ Z(!) = e X (!) Proof: 1 X −j!n Z(!) = x[n − n0]e n=−∞ 1 X −j!(m+n0) = x[m]e (where m = n − n0) m=−∞ = e−j!n0 X (!) Review DTFT DTFT Properties Examples Summary 3. Frequency Shift Property Shifting in frequency is the same as multiplying by a complex exponential in time: j!0n z[n] = x[n]e $ Z(!) = X (! − !0) Proof: 1 X Z(!) = x[n]ej!0ne−j!n n=−∞ 1 X = x[n]e−j(!−!0)n n=−∞ = X (! − !0) Review DTFT DTFT Properties Examples Summary 4. Convolution Property Convolving in time is the same as multiplying in frequency: y[n] = h[n] ∗ x[n] $ Y (!) = H(!)X (!) Proof: Remember that y[n] = h[n] ∗ x[n] means that P1 y[n] = m=−∞ h[m]x[n − m]. Therefore, 1 1 ! X X Y (!) = h[m]x[n − m] e−j!n n=−∞ m=−∞ 1 1 X X = (h[m]x[n − m]) e−j!me−j!(n−m) m=−∞ n=−∞ 1 ! 0 1 1 X −j!m X −j!(n−m) = h[m]e @ x[n − m]e A m=−∞ (n−m)=−∞ = H(!)X (!) Review DTFT DTFT Properties Examples Summary Outline 1 Review: Frequency Response 2 Discrete Time Fourier Transform 3 Properties of the DTFT 4 Examples 5 Summary Review DTFT DTFT Properties Examples Summary Impulse and Delayed Impulse For our examples today, let's consider different combinations of these three signals: f [n] = δ[n] g[n] = δ[n − 3] h[n] = δ[n − 6] Remember from last time what these mean: ( 1 n = 0 f [n] = 0 otherwise ( 1 n = 3 g[n] = 0 otherwise ( 1 n = 6 h[n] = 0 otherwise Review DTFT DTFT Properties Examples Summary DTFT of an Impulse First, let's find the DTFT of an impulse: ( 1 n = 0 f [n] = 0 otherwise 1 X F (!) = f [n]e−j!n n=−∞ = 1 × e−j!0 = 1 So we get that f [n] = δ[n] $ F (!) = 1.