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.

214 Chapter 6 IMPLEMENTATION OF DISCRETE-TIME FILTERS<br />

6.2 IIR FILTER STRUCTURES<br />

The system function of an IIR filter is given by<br />

H(z) = B(z) ∑ M<br />

A(z) = n=0 b nz −n<br />

∑ N<br />

n=0 a nz = b 0 + b 1 z −1 + ···+ b M z −M<br />

−n 1+a 1 z −1 + ···+ a N z −N ; a 0 =1<br />

(6.1)<br />

where b n and a n are the coefficients of the filter. We have assumed without<br />

loss of generality that a 0 =1. The order of such an IIR filter is called N if<br />

a N ̸=0.The difference equation representation of an IIR filter is expressed<br />

as<br />

M∑<br />

N∑<br />

y(n) = b m x(n − m) − a m y(n − m) (6.2)<br />

m=0<br />

m=1<br />

Three different structures can be used to implement an IIR filter:<br />

1. Direct form: In this form the difference equation (6.2) is implemented<br />

directly as given. There are two parts to this filter, namely the moving<br />

average part and the recursive part (or equivalently, the numerator<br />

and denominator parts). Therefore this implementation leads to two<br />

versions: direct form I and direct form II structures.<br />

2. Cascade form: In this form the system function H(z) inequation<br />

(6.1) is factored into smaller 2nd-order sections, called biquads. The<br />

system function is then represented as a product of these biquads. Each<br />

biquad is implemented in a direct form, and the entire system function<br />

is implemented as a cascade of biquad sections.<br />

3. Parallel form: This is similar to the cascade form, but after factorization,<br />

a partial fraction expansion is used to represent H(z) asasum<br />

of smaller 2nd-order sections. Each section is again implemented in a<br />

direct form, and the entire system function is implemented as a parallel<br />

network of sections.<br />

We will briefly discuss these forms in this section. IIR filters are generally<br />

described using the rational form version (or the direct form structure)<br />

of the system function. Hence we will provide MATLAB functions for<br />

converting direct form structures to cascade and parallel form structures.<br />

6.2.1 DIRECT FORM<br />

As the name suggests, the difference equation (6.2) is implemented as<br />

given using delays, multipliers, and adders. For the purpose of illustration,<br />

let M = N =4.Then the difference equation is<br />

y(n) =b 0 x(n)+b 1 x(n − 1) + b 2 x(n − 2) + b 3 x(n − 3) + b 4 x(n − 4)<br />

− a 1 y(n − 1) − a 2 y(n − 2) − a 3 y(n − 3) − a 4 y(n − 4)<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!