| Line 5: | Line 5: | ||
[[Category:lecture notes]] | [[Category:lecture notes]] | ||
| + | == [[ECE_438_Spring_2009_mboutin_Course_Notes|Course Notes]], February 23, 2009 == | ||
| + | [[ECE438_%28BoutinSpring2009%29|ECE438, Spring 2009]] | ||
| + | ---- | ||
== 1.6.3 Spectral analysis via DFT == | == 1.6.3 Spectral analysis via DFT == | ||
Latest revision as of 05:44, 16 September 2013
Contents
Course Notes, February 23, 2009
1.6.3 Spectral analysis via DFT
use DFT to approximate $ X(a) $ for a DT signal x(n)
Relationships of signals
<p>infinate duration signal: $ x(n) \Rightarrow DTFT \Rightarrow X(w) $
also infinate druation signal: $ x(n) \Rightarrow *truncation\Rightarrow \bar{x}(n) = x[n] $ finite duration
The periodic sample of $ x_p(n) = \sum_{l=-\infty}^{\infty} \bar{x}(n+lN) $
The DTFT of periodic sample is $ X(k) $
The DTFT of $ \bar{x}(n) \Rightarrow \bar{X}(w) \Rightarrow sampling\Rightarrow \bar{X}(k\frac{2\pi}{N}) = X(k) $
Example
$ x(n) = cos(w_0n) $
$ X(w) = rep_{2\pi}(\pi \delta(w-w_0) + \pi \delta(w+w_0)) $
consider $ \bar{x}(n) = x(n)w(n) $
- $ w(n) = 1 | 0 \le n \le N, 0 | else $
$ \bar{X}(w) = F(\bar{x}(n)) $
$ = \frac{1}{2\pi}X(w)*W(w) $
$ = \frac{1}{2}(W(w+w_0) + W(w-w_0)) $
Two sources of inaccuracies
- signal truncation $ \Rightarrow $ "leakage"
- frequency sampling $ \Rightarrow $ "picket fence effect"
1.6.4 FFT "Fast Fourier transform"
An algorithm (family of algo) to compute the DFT fast
Recall DFT of x(n) periodic w\ period N
$ X(k) = \sum_{n=0}^{N-1} x(n)e^{-j\frac{2\pi}{N}kn} $
$ e^{-j\frac{2\pi}{N}^r} $ is an Nth root of unity for every r=0,1,2,...,N-1
N=2
two square roots of unity
$ e^{-j\frac{2\pi}{N}*0} = 1 $
$ e^{-j\frac{2\pi}{N}*1r} = e^{-j\pi} = -1 $
Let $ W_N^r = e^{-j\frac{2\pi}{N}r} $for every r, $ W_N^r $ is an Nth root of unity
--Drestes 15:05, 23 February 2009 (UTC)
