EC 2314 Digital Signal Processing Filter Realization By

EC 2314 Digital Signal Processing -Filter Realization By Dr. K. Udhayakumar

Introduction n n In earlier chapters we studied theory of discrete systems in both the time and frequency domains. We will now use this theory for the processing of digital signals. To process signals, we have to design and implement systems called filters. The filter design issue is influenced by such factors as n n n The type of the filter: IIR or FIR The form of its implementation: structures Different filter structures dictate different design strategies.

Introduction n IIR filters are characterized by infinite-duration impulse response. Some of these impulse responses can be modeled by n n Rational system functions Difference equations ARMA or recursive filters We will treat FIR filter separately from IIR filters for both design and implementation purposes.

Introduction n Since our filters are LTI systems, we need the following three elements to describe digital filter structures. n n n Adder Multiplier (Gain) Delay element (shift or memory)

IIR Filter Structures n The system function of an IIR filter is given by The order of such an IIR filter is called N if a. N~=0. The difference equation representation of an IIR filter is expressed as

Three different structures can be used to implement an IIR filter: n Direct form n n n Cascade form n n n In this form, there are two parts to this filter, the moving average part and the recursive part (or the numerator and denominator parts) Two version: direct form I and direct form II The system function H(z) is factored into smaller second-order sections, called biquads. H(z) is then represented as a product of these biquads. Each biquad is implemented in a direct form, and the entire system function is implemented as a cascade of biquad sections. Parallel form n n n H(z) is represented as a sum of smaller second-order sections. Each section is again implemented in a direct form. The entire system function is implemented as a parallel network of sections.

Direct Form I Structure n n As the name suggests, the difference equation is implemented as given using delays, multipliers, and adders. For the purpose of illustration, Let M=N=2,

Direct Form II Structure The commutative law of the convolution Direct Form II structure

Matlab Implementation n In Matlab the direct from structure is described by two row vectors: b containing the {bn} coefficients and a containing the {an} coefficients. The structure is implemented by the filter function, which is discussed in Chapter 2.

Cascade Form n n In this form the system function H(z) is written as a product of second-order section with real coefficients. This is done by factoring the numerator and denominator polynomials into their respective roots and then combining either a complex conjugate root pair or any two real roots into second-order polynomials.

Cascade Form n We assume that N is an even integer. Then Where, K is equal to N/2, and Bk, 1, Bk, 2, Ak, 1, Ak, 2 are real numbers representing the coefficients of second-order section.

Biquad Section Is called the k-th biquad section. The input to the k-th biquad section is the output from the (k-1)-th section, while the output from the k-th biquad is the input to the (k+1)-th biquad. Each biquad section can be implemented in direct form II.

The entire filter is then implemented as a cascade of biquads Cascade form structure for N=4

Matlab Implementation n Given the coefficients {bn} and {an} of the direct form filter The function [b 0, B, A] = dir 2 cas(b, a) can be used to obtain the coefficients b 0, {Bk, i}, and {Ak, i}. The cascade form is implemented using a casfiltr function Function y = casfiltr(b 0, B, A, x) Function cas 2 dir converts a cascade form to a direct form. Function [b, a] = cas 2 dir(b 0, B, A) n Examples: 6. 1 n n n

Parallel Form n In this form the system function H(z) is written as a sum of second-order section using partial fraction expansion. K=N/2, and B, A are real numbers

The second-order section Is the k-th proper rational biquad section. The filter input is available to all biquad section as well as to the polynomial section if M>=N (which is an FIR part) The output from these sections is summed to form the filter output. Each biquad section can be implemented in direct form II.

Parallel form structure for N=4 (M=N=4)

Matlab Implementation n n The function dir 2 par converts the direct form coefficients {bn} and {an} into parallel form coefficients {Bk, i} and {Ak, i} Functions: n n n [C, B, A] = dir 2 par(b, a) I = cplxcomp(p 1, p 2) y = parfiltr(C, B, A, x) [b, a] = par 2 dir(C, B, A) Examples

FIR Filter Structure A finite-duration impulse response filter has a system function of the form: Hence the impulse response h(n) is And the difference equation representation is Which is a linear convolution of finite support. The order of the filter is M-1, while the length of the filter is M.

FIR Filter Structure n Direct form: n n The difference equation is implemented as a tapped delay line since there are no feedback paths. Figure 6. 10 Note that since the denominator is equal to unity, there is only one direct form structure. Matlab implementation n Function: y = filter(b, 1, x)

FIR Filter Structure n Cascade form: n n Figure 6. 11 Matlab Implementation n Function: dir 2 cas, cas 2 dir

FIR Filter Structure n Linear-phase form n n For frequency-selective filters (e. g. , lowpass filters) it is generally desirable to have a phase response that is a linear function of frequency. That is For a causal FIR filter with impulse over [0, M-1] interval, the linear-phase conditions Symmetric impulse response vs. antisymmetric impulse response

Linear-phase form n n n Consider the difference equation with a symmetric impulse response. Figure 6. 12: M=7(odd) and M=6(even) Matlab implementation n The linear-phase structure is essentially a direct form draw differently to save on multiplications. Hence in a Matlab implementation the linearphase structure is equivalent to the direct form.

Frequency Sampling Form n In this form we use the fact that the system function H(z) of an FIR filter can be reconstructed from it samples on the unit circle It is also interesting to note that the FIR filter described by the above equation has a recursive form similar to an IIR filter because it contains both poles and zeros.

Frequency Sampling Form n n n The system function leads to a parallel structure as shown in Figure 6. 15 for M=4. One problem with the structure in Fig. 6. 15 is that it requires a complex arithmetic implementation. Using the symmetry properties of the DFT and the (WM–k) factor.

let p 1=WM-k=exp(-j*2*pi*k/M) =cos(2*pi*k/M)+jsin(2*pi*k/M) Hk=mag. Hk*exp(j*pha. Hk) =mag. Hk*(cos(pha. Hk)+jsin(pha. Hk)) Then: n Figure 6. 16

MATLAB functions and Examples n n n Notify the structure type [C, B, A]=dir 2 fs(h) of textbook A practical problem: unstable, avoid this problem by sampling H(z) on a circle |z|=r. Ex 6. 6 find frequency sample form Ex 6. 7 In this example the frequency sample form has less computational complexity than direct form.

Lattice Filter Structure n n The lattice filter is extensively used in digital speech processing and in implementation of adaptive filter. It is a preferred form of realization over other FIR or IIR filter structures because in speech analysis and in speech synthesis the small number of coefficients allows a large number of formants to be modeled in real-time. n n All-zeros lattice is the FIR filter representation of the lattice filter. The lattice ladder is the IIR filter representation.

All-zero Lattice Filters n An FIR filter of length M (or order M-1) has a lattice structure with M-1 stages. Km: reflection coefficients

All-zero lattice filters If the FIR filter is given by the direct form And if we denote the polynomial Then the lattice filter coefficients {Km} can be obtained by the following recursive algorithm

Note that the above algorithm will fail if |Km|=1 for any m. Clearly, this condition is satisfied by the linear-phase FIR filter. Therefore, linear-phase FIR filter cannot be implemented using lattice structure.

Matlab Implementation n Functions: [K] = dir 2 latc(b) [y] = latcfilt(K, x) [b] = latc 2 dir(K) n Example 6. 8 n n n

All-pole Lattice Filter n n A lattice structure for an IIR filter is restricted to an all-pole system unction. It can be developed from an FIR lattice structure. This IIR filter of order N has a lattice structure with N stages as shown in Fig. 6. 20. Each stage of the filter has an input and output that are related by the order-recursive equations.

All-pole Lattice Filter

Matlab Implementation n n Function [K] = dir 2 latc(a) Care must be taken to ignore the K 0 coefficient in the K array. Function [a] = latc 2 dir(K) K 0=1

Lattice ladder Filters n A general IIR filter containing both poles and zeros can be realized as a lattice-type structure by using an all-pole lattice as the basic building block. Consider an IIR filter with system function n Where, without loss of generality, we assume that N>=M n

Lattice ladder Filters n A lattice type structure can be constructed by first realizing an all-poles lattice with coefficient {Km} for the denominator and then adding a ladder part by taking the output as a weighted linear combination of {gm(n)} as shown in Fig. 6. 22 for M=N.

Lattice ladder Filters The output of the lattice-ladder structure is given by {Cm} are called the ladder coefficients that determine the zeros of the system function H(z).

Matlab Implementation n Function [K, C] = dir 2 ladr(b, a) n n n To use this function, N>=M. If M>N, then the numerator AN(z) should be divided into the denominator BM(z) using the deconv function to obtain a proper rational part and a polynomial part. The proper rational part can be implemented using the lattice-ladder structure, while the polynomial part is implemented using the direct structure. Function [b, a] ladr 2 dir(K, C) Function [y] = ladrfilt(K, C, x)

Readings and exercises n Textbook : pp 182~217 Chinese ref. Book: pp. 128~137, 231~234 n Exercises : n n n 1. 6. 2, 6. 3 both except (e) 2. 6. 2 e, 6. 3 e, 6. 8, 选6. 10