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.

Decimation by a Factor D 483<br />

The output of the filter is a sequence v(n) given as<br />

v(n) △ =<br />

∞∑<br />

h(k)x(n − k) (9.21)<br />

k=0<br />

which is then downsampled by the factor D to produce y(m). Thus,<br />

y(m) =v(mD) =<br />

∞∑<br />

h(k)x(mD − k) (9.22)<br />

k=0<br />

Although the filtering operation on x(n) islinear and time invariant, the<br />

downsampling operation in combination with the filtering results also in<br />

a time-variant system.<br />

The frequency-domain characteristics of the output sequence y(m) obtained<br />

through the filtered signal v(n) can be determined by following the<br />

analysis steps given before—i.e., by relating the spectrum of y(m) tothe<br />

spectrum of the input sequence x(n). Using these steps, we can show that<br />

or that<br />

Y (z) = 1 D<br />

Y (ω y )= 1 D<br />

D−1<br />

∑<br />

k=0<br />

H ( e −j2πk/D z 1/D) X ( e −j2πk/D z 1/D) (9.23)<br />

D−1<br />

∑<br />

k=0<br />

(<br />

ωy − 2πk<br />

H<br />

D<br />

)<br />

X<br />

( )<br />

ωy − 2πk<br />

D<br />

(9.24)<br />

With a properly designed filter H D (ω), the aliasing is eliminated and,<br />

consequently, all but the first term in (9.24) vanish. Hence,<br />

Y (ω y )= 1 D H D<br />

( ωy<br />

)<br />

X<br />

D<br />

( ωy<br />

)<br />

D<br />

= 1 (<br />

D X ωy<br />

)<br />

D<br />

(9.25)<br />

for 0 ≤|ω y |≤π. The spectra for the sequences x(n), h(n), v(n), and<br />

y(m) are illustrated in Figure 9.7.<br />

MATLAB Implementation MATLAB provides the function y =<br />

decimate(x,D) that resamples the sequence in array x at 1/D times<br />

the original sampling rate. The resulting resampled array y is D times<br />

shorter—i.e., length(y) = length(x)/D. Anideal lowpass filter given<br />

in (9.20) is not possible in the MATLAB implementation; however, fairly<br />

accurate approximations are used. The default lowpass filter used in the<br />

function is an 8th-order Chebyshev type-I lowpass filter with the cutoff<br />

frequency of 0.8π/D. Using additional optional arguments, the filter order<br />

can be changed or an FIR filter of specified order and cutoff frequency<br />

can be used.<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!