15.01.2013 Views

U. Glaeser

U. Glaeser

U. Glaeser

SHOW MORE
SHOW LESS

Create successful ePaper yourself

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

frequency response involves some aliasing:<br />

where H a(s) is the Laplace transform system function of the analog filter. The aliasing effect usually causes<br />

only a slight perturbation of the digital filter with respect to the analog filter. The system function of the<br />

analog filter can be expressed in partial fraction form<br />

After sampling the digital filter has a frequency response that is also a rational form:<br />

where b[n] and a[n]are the coefficients of the designed filter. Impulse invariance is equivalent to a linear<br />

mapping of the analog frequency range [−π /T d, π /T d] into the digital frequency range [−π,π].<br />

Bilinear Transform<br />

On the other hand, the bilinear transformation performs a nonlinear mapping of the whole analog<br />

frequency range [−∞, ∞]into the finite digital frequency range [−π, π]. The mapping of the s-plane to<br />

the z-plane is done by the bilinear transform:<br />

The resulting correspondence between the analog and digital frequency domains is a tangent function:<br />

Despite the nonlinear nature of the mapping, it is relatively easy to turn the digital design specification<br />

into an analog design specification. The resulting filter is IIR and the filter coefficients can be computed<br />

with an algebraic form. The bilinear transform method is usually applied to four classical analog filter<br />

frequency selective filters: Butterworth, Chebyshev-I, Chebyshev-II, and elliptic filters. All these are well<br />

known for their frequency-selective behavior as lowpass, bandpass, or highpass filters. When using the<br />

bilinear mapping, elliptic IIR filters turn out to have the best magnitude response for given filter order,<br />

but elliptic filters have severe phase distortion, which can be a significant problem in advanced DSP<br />

applications such as telecommunications.<br />

Windowing<br />

IIR filter designs have poor phase response, so interest in FIR filters has always been strong. If the<br />

coefficients of an FIR filter are real and symmetric b[k] = b∗ [M − k] then the filter will have perfectly<br />

linear phase. The first attempt to design FIR filters in the 1960s was to truncate the inverse DTFT of the<br />

ideal frequency response (which is the impulse response h[n]of the ideal filter), so that the filter is<br />

1<br />

symmetric and linear-phase. This requires the ideal filter to have linear-phase with slope −-- M, where M<br />

© 2001 by CRC Press LLC<br />

H(ω) Ha j ω<br />

∞<br />

= ∑ Td k=−∞<br />

Ha() s<br />

=<br />

---- j 2π ⎛ + -----k ⎞<br />

⎝ ⎠<br />

N<br />

∑<br />

k=1<br />

A k<br />

--------s–<br />

sk T<br />

DTFT: H(ω)<br />

dAk 1 e Td s N<br />

b[ n]e<br />

∑-------------------------<br />

k – jω<br />

k=1 – e jωn –<br />

k=0<br />

1 a[ n]e<br />

jωn –<br />

= = --------------------------------------<br />

N<br />

+<br />

s<br />

ω a<br />

1<br />

2 1 z<br />

----<br />

Td –<br />

–<br />

1 z 1 –<br />

= ⎛--------------- ⎞<br />

⎝ + ⎠<br />

=<br />

2<br />

---- tan<br />

Td ω ⎛--- ⎞<br />

⎝2⎠ T d<br />

∑<br />

N−1<br />

∑<br />

k=1<br />

2

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

Saved successfully!

Ooh no, something went wrong!