The Discrete Fourier Transform Sampling the Fourier Transform • Consider an aperiodic sequence with a Fourier transform DT FT x[n] X e j • Assume that a sequence is obtained by sampling the DTFT ~ Xk X e j X e j2 / Nk 2 / Nk • Since the DTFT is periodic resulting sequence is also periodic • We can also write it in terms of the z-transform ~ Xk Xz z e2 / Nk X e j2 / Nk • • • DSP The sampling points are shown in figure ~ Xk could be the DFS of a sequence Write the corresponding sequence 1 N1 ~ ~ x[n] Xk e j2 / Nkn N k 0 Sampling the Fourier Transform Cont’d • The only assumption made on the sequence is that DTFT exist xme Xe j ~ Xk X e j2 / Nk jm m 1 N1 ~ ~ x[n] Xk e j2 / Nkn N k 0 • Combine equation to get 1 N 1 ~ x[n] xme j2 / Nkm e j2 / Nkn N k 0 m 1 N 1 j2 / Nk n m ~ xm e x m p n m m N k 0 m • Term in the parenthesis is • So we get 1 N1 j2 / Nk n m ~ p n m e N k 0 ~ x[n] xn DSP n m rN r r r n rN xn rN Sampling the Fourier Transform Cont’d DSP Sampling the Fourier Transform Cont’d • Samples of the DTFT of an aperiodic sequence – can be thought of as DFS coefficients – of a periodic sequence – obtained through summing periodic replicas of original sequence • If the original sequence – is of finite length – and we take sufficient number of samples of its DTFT – the original sequence can be recovered by ~ x n 0 n N 1 xn else 0 • It is not necessary to know the DTFT at all frequencies – To recover the discrete-time sequence in time domain • Discrete Fourier Transform – Representing a finite length sequence by samples of DTFT DSP The Discrete Fourier Transform • Consider a finite length sequence x[n] of length N xn 0 outside of 0 n N 1 • For given length-N sequence associate a periodic sequence ~ x n xn rN r • The DFS coefficients of the periodic sequence are samples of the DTFT of x[n] • Since x[n] is of length N there is no overlap between terms of x[n-rN] and we can write the periodic sequence as ~ x n xn mod N xnN • To maintain duality between time and frequency ~ – We choose one period of Xk as the Fourier transform of x[n] ~ ~ Xk 0 k N 1 Xk Xk mod N Xk N Xk 0 else DSP The Discrete Fourier Transform Cont’d • The DFS pair ~ Xk N 1 ~ x[n]e j2 / Nkn n0 1 N1 ~ ~ x[n] Xk e j2 / Nkn N k 0 • The equations involve only on period so we can write N1 ~ x[n]e j2 / Nkn 0 k N 1 Xk n0 0 else 1 N1 ~ Xk e j2 / Nkn 0 k N 1 x[n] N k 0 0 else • The Discrete Fourier Transform Xk N 1 x[n]e n0 j2 / N kn 1 N 1 x[n] Xk e j2 / Nkn N k 0 • The DFT pair can also be written as DFT Xk x[n] DSP Example • The DFT of a rectangular pulse • x[n] is of length 5 • We can consider x[n] of any length greater than 5 • Let’s pick N=5 • Calculate the DFS of the periodic form of x[n] ~ Xk 4 j2 k / 5 n e n0 1 e j2 k 1 e j2 k / 5 5 k 0,5,10,... else 0 DSP Example Cont’d • If we consider x[n] of length 10 • We get a different set of DFT coefficients • Still samples of the DTFT but in different places DSP Properties of DFT • Linearity x1 n DFT x2 n DFT X1 k X2 k ax1 n bx2 n DFT aX1 k bX2 k • Duality xn DFT Xk Xn DFT Nx k N • Circular Shift of a Sequence xn DFT Xk xn mN 0 n N - 1 DFT Xk e j2k / Nm DSP Example: Duality DSP Symmetry Properties DSP Circular Convolution • Circular convolution of of two finite length sequences x3 n N 1 x mx n m 1 2 N m0 x3 n N 1 x mx n m 2 m0 DSP 1 N Example • Circular convolution of two rectangular pulses L=N=6 1 0 n L 1 x1 n x2 n else 0 • DFT of each sequence X1 k X2 k N 1 e j 2 kn N n0 N k 0 0 else • Multiplication of DFTs N2 k 0 X3 k X1 k X2 k 0 else • And the inverse DFT N 0 n N 1 x3 n else 0 DSP Example • We can augment zeros to each sequence L=2N=12 • The DFT of each sequence X1 k X2 k 1e j 1e 2 Lk N j 2 k N • Multiplication of DFTs 2 Lk j N 1 e X3 k 2 k j 1e N DSP 2
© Copyright 2026 Paperzz