![Page 1: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/1.jpg)
1Copyright © 2001, S. K. Mitra
FIR filters
ELG6163Miodrag Bolic
![Page 2: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/2.jpg)
2Copyright © 2001, S. K. Mitra
Outline
• FIR filters– Structures– Polyphase FIR filters– Parallel polyphase FIR– Decimated FIR
• Implementations of FIR filters
![Page 3: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/3.jpg)
3Copyright © 2001, S. K. Mitra
Canonic and Canonic and NoncanonicNoncanonicStructuresStructures
• A digital filter structure is said to be canonic if the number of delays in the block diagram representation is equal to the order of the transfer function
• Otherwise, it is a noncanonic structure
![Page 4: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/4.jpg)
4Copyright © 2001, S. K. Mitra
Canonic and Canonic and NoncanonicNoncanonicStructuresStructures
• The structure shown below is noncanonic as it employs two delays to realize a first-order difference equation
]1[][]1[][ 101 −++−−= nxpnxpnydny
![Page 5: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/5.jpg)
5Copyright © 2001, S. K. Mitra
Basic FIR Digital Filter Basic FIR Digital Filter StructuresStructures
• A causal FIR filter of order N is characterized by a transfer function H(z) given by
which is a polynomial in• In the time-domain the input-output relation
of the above FIR filter is given by
∑ =−= N
nnznhzH 0 ][)(
1−z
∑ = −= Nk knxkhny 0 ][][][
![Page 6: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/6.jpg)
6Copyright © 2001, S. K. Mitra
Direct Form FIR Digital Filter Direct Form FIR Digital Filter StructuresStructures
• An FIR filter of order N is characterized by N+1 coefficients and, in general, require N+1 multipliers and N two-input adders
• Structures in which the multiplier coefficients are precisely the coefficients of the transfer function are called direct formstructures
![Page 7: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/7.jpg)
7Copyright © 2001, S. K. Mitra
Direct Form FIR Digital Filter Direct Form FIR Digital Filter StructuresStructures
• A direct form realization of an FIR filter can be readily developed from the convolution sum description as indicated below for N = 4
![Page 8: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/8.jpg)
8Copyright © 2001, S. K. Mitra
Direct Form FIR Digital Filter Direct Form FIR Digital Filter StructuresStructures
• The transpose of the direct form structure shown earlier is indicated below
• Both direct form structures are canonic with respect to delays
![Page 9: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/9.jpg)
9Copyright © 2001, S. K. Mitra
Cascade Form FIR Digital Cascade Form FIR Digital Filter StructuresFilter Structures
• A higher-order FIR transfer function can also be realized as a cascade of second-order FIR sections and possibly a first-order section
• To this end we express H(z) as
where if N is even, and if N is odd, with
∏ =−− ++= K
k kk zzhzH 12
21
110 )(][)( ββ
2NK =
21+= NK
02 =Kβ
![Page 10: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/10.jpg)
10Copyright © 2001, S. K. Mitra
Cascade Form FIR Digital Cascade Form FIR Digital Filter StructuresFilter Structures
• A cascade realization for N = 6 is shown below
• Each second-order section in the above structure can also be realized in the transposed direct form
![Page 11: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/11.jpg)
11Copyright © 2001, S. K. Mitra
LinearLinear--Phase FIR StructuresPhase FIR Structures• The symmetry (or antisymmetry) property of a
linear-phase FIR filter can be exploited to reduce the number of multipliers into almost half of that in the direct form implementations
• Consider a length-7 Type 1 FIR transfer function with a symmetric impulse response:
321 3210 −−− +++= zhzhzhhzH ][][][][)(654 012 −−− +++ zhzhzh ][][][
![Page 12: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/12.jpg)
12Copyright © 2001, S. K. Mitra
LinearLinear--Phase FIR StructuresPhase FIR Structures• Rewriting H(z) in the form
we obtain the realization shown below
)]([)]([)( 516 110 −−− +++= zzhzhzH342 32 −−− +++ zhzzh ][)]([
![Page 13: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/13.jpg)
13Copyright © 2001, S. K. Mitra
LinearLinear--Phase FIR StructuresPhase FIR Structures
• A similar decomposition can be applied to a Type 2 FIR transfer function
• For example, a length-8 Type 2 FIR transfer function can be expressed as
• The corresponding realization is shown on the next slide
)]([)]([)( 617 110 −−− +++= zzhzhzH
)]([)]([ 4352 32 −−−− ++++ zzhzzh
![Page 14: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/14.jpg)
14Copyright © 2001, S. K. Mitra
LinearLinear--Phase FIR StructuresPhase FIR Structures
• Note: The Type 1 linear-phase structure for a length-7 FIR filter requires 4 multipliers, whereas a direct form realization requires 7 multipliers
![Page 15: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/15.jpg)
15Copyright © 2001, S. K. Mitra
LinearLinear--Phase FIR StructuresPhase FIR Structures
• Note: The Type 2 linear-phase structure for a length-8 FIR filter requires 4 multipliers, whereas a direct form realization requires 8 multipliers
• Similar savings occurs in the realization of Type 3 and Type 4 linear-phase FIR filters with antisymmetric impulse responses
![Page 16: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/16.jpg)
16Copyright © 2001, S. K. Mitra
PolyphasePolyphase FIR StructuresFIR Structures
• The polyphase decomposition of H(z) leads to a parallel form structure
• To illustrate this approach, consider a causal FIR transfer function H(z) with N = 8:
4321 43210 −−−− ++++= zhzhzhzhhzH ][][][][][)(8765 8765 −−−− ++++ zhzhzhzh ][][][][
![Page 17: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/17.jpg)
17Copyright © 2001, S. K. Mitra
PolyphasePolyphase FIR StructuresFIR Structures• H(z) can be expressed as a sum of two
terms, with one term containing the even-indexed coefficients and the other containing the odd-indexed coefficients:
)][][][][][()( 8642 86420 −−−− ++++= zhzhzhzhhzH)][][][][( 7531 7531 −−−− ++++ zhzhzhzh
)][][][][][( 8642 86420 −−−− ++++= zhzhzhzhh)][][][][( 6421 7531 −−−− ++++ zhzhzhhz
![Page 18: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/18.jpg)
18Copyright © 2001, S. K. Mitra
PolyphasePolyphase FIR StructuresFIR Structures
• By using the notation
we can express H(z) as
3211 7531 −−− +++= zhzhzhhzE ][][][][)(
3210 86420 4−−−− ++++= zhzhzhzhhzE ][][][][][)(
)()()( 21
120 zEzzEzH −+=
![Page 19: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/19.jpg)
19Copyright © 2001, S. K. Mitra
PolyphasePolyphase FIR StructuresFIR Structures
• In a similar manner, by grouping the terms in the original expression for H(z), we can re-express it in the form
where now)()()()( 3
223
113
0 zEzzEzzEzH −− ++=
210 630 −− ++= zhzhhzE ][][][)(
211 741 −− ++= zhzhhzE ][][][)(
212 852 −− ++= zhzhhzE ][][][)(
![Page 20: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/20.jpg)
20Copyright © 2001, S. K. Mitra
PolyphasePolyphase FIR StructuresFIR Structures
• The decomposition of H(z) in the form
or
is more commonly known as the polyphasedecomposition
)()()()( 32
231
130 zEzzEzzEzH −− ++=
)()()( 21
120 zEzzEzH −+=
![Page 21: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/21.jpg)
21Copyright © 2001, S. K. Mitra
PolyphasePolyphase FIR StructuresFIR Structures
• In the general case, an L-branch polyphasedecomposition of an FIR transfer function of order N is of the form
where
with h[n]=0 for n > N
)()( Lm
Lm
m zEzzH ∑ −=
−= 10
∑+
=
−+=LN
n
mm zmLnhzE
/)(][)(
1
0
![Page 22: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/22.jpg)
22Copyright © 2001, S. K. Mitra
PolyphasePolyphase FIR StructuresFIR Structures• Figures below show the 4-branch, 3-
branch, and 2-branch polyphaserealization of a transfer function H(z)
• Note: The expression for the polyphasecomponents are different in each case)(zEm
![Page 23: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/23.jpg)
23Copyright © 2001, S. K. Mitra
PolyphasePolyphase FIR StructuresFIR Structures
• The subfilters in the polyphaserealization of an FIR transfer function are also FIR filters and can be realized using any methods described so far
• However, to obtain a canonic realization of the overall structure, the delays in all subfilters must be shared
)( Lm zE
![Page 24: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/24.jpg)
24Copyright © 2001, S. K. Mitra
PolyphasePolyphase FIR StructuresFIR Structures• Figure below shows a canonic realization of
a length-9 FIR transfer function obtained using delay sharing
![Page 25: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/25.jpg)
25Copyright © 2001, S. K. Mitra
FIR Filter Structures Based on FIR Filter Structures Based on PolyphasePolyphase DecompositionDecomposition
• We shall demonstrate later that a parallel realization of an FIR transfer function H(z) based on the polyphase decomposition can often result in computationally efficient multirate structures
• Consider the M-branch Type I polyphasedecomposition of H(z):
)()( 10
Mk
Mk
k zEzzH ∑ −=
−=
![Page 26: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/26.jpg)
26Copyright © 2001, S. K. Mitra
FIR Filter Structures Based on FIR Filter Structures Based on PolyphasePolyphase DecompositionDecomposition
• A direct realization of H(z) based on the Type I polyphase decomposition is shown below
![Page 27: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/27.jpg)
27Copyright © 2001, S. K. Mitra
FIR Filter Structures Based on FIR Filter Structures Based on PolyphasePolyphase DecompositionDecomposition
• The transpose of the Type I polyphase FIR filter structure is indicated below
![Page 28: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/28.jpg)
28Copyright © 2001, S. K. Mitra
Computationally Efficient Computationally Efficient DecimatorsDecimators
• Consider first the single-stage factor-of-Mdecimator structure shown below
• We realize the lowpass filter H(z) using the Type I polyphase structure as shown on the next slide
M][nx )(zH ][nyv[n]
![Page 29: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/29.jpg)
29Copyright © 2001, S. K. Mitra
Computationally Efficient Computationally Efficient DecimatorsDecimators
• Using the cascade equivalence #1 we arrive at the computationally efficient decimator structure shown below on the right
Decimator structure based on Type I polyphase decomposition
y[n] y[n]x[n]x[n]
][nv
![Page 30: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/30.jpg)
30Copyright © 2001, S. K. Mitra
Computationally Efficient Computationally Efficient DecimatorsDecimators
• To illustrate the computational efficiency of the modified decimator structure, assumeH(z) to be a length-N structure and the input sampling period to be T = 1
• Now the decimator output y[n] in the original structure is obtained by down-sampling the filter output v[n] by a factor ofM
![Page 31: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/31.jpg)
31Copyright © 2001, S. K. Mitra
Computationally Efficient Computationally Efficient DecimatorsDecimators
• It is thus necessary to compute v[n] at
• Computational requirements are therefore Nmultiplications and additions per output sample being computed
• However, as n increases, stored signals in the delay registers change
...,2,,0,,2,... MMMMn −−=
)1( −N
![Page 32: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/32.jpg)
32Copyright © 2001, S. K. Mitra
Computationally Efficient Computationally Efficient DecimatorsDecimators
• Hence, all computations need to be completed in one sampling period, and for the following sampling periods the arithmetic units remain idle
• The modified decimator structure also requires N multiplications and additions per output sample being computed
)1( −N
)1( −M
![Page 33: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/33.jpg)
33Copyright © 2001, S. K. Mitra
Computationally Efficient Computationally Efficient Decimators and InterpolatorsDecimators and Interpolators
• However, here the arithmetic units are operative at all instants of the output sampling period which is 1/M times that of the input sampling period
• Similar savings are also obtained in the case of the interpolator structure developed using the polyphase decomposition
![Page 34: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/34.jpg)
34Copyright © 2001, S. K. Mitra
Computationally Efficient Computationally Efficient InterpolatorsInterpolators
• Figures below show the computationally efficient interpolator structures
Interpolator based on Type I polyphase decomposition
![Page 35: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/35.jpg)
35Copyright © 2001, S. K. Mitra
Computationally Efficient Computationally Efficient Decimators and InterpolatorsDecimators and Interpolators
• More efficient interpolator and decimator structures can be realized by exploiting the symmetry of filter coefficients in the case of linear-phase filters H(z)
• Consider for example the realization of afactor-of-3 (M = 3) decimator using alength-12 Type 1 linear-phase FIR lowpassfilter
![Page 36: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/36.jpg)
36Copyright © 2001, S. K. Mitra
Computationally Efficient Computationally Efficient Decimators and InterpolatorsDecimators and Interpolators
• The corresponding transfer function is
• A conventional polyphase decomposition ofH(z) yields the following subfilters:
54321 ]5[]4[]3[]2[]1[]0[)( −−−−− +++++= zhzhzhzhzhhzH109876 ]0[]1[]2[]3[]4[]5[ 11−−−−−− ++++++ zhzhzhzhzhzh
3210 ]2[]5[]3[]0[)( −−− +++= zhzhzhhzE
3211 ]1[]4[]4[]1[)( −−− +++= zhzhzhhzE
3212 ]0[]3[]5[]2[)( −−− +++= zhzhzhhzE
![Page 37: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/37.jpg)
37Copyright © 2001, S. K. Mitra
Computationally Efficient Computationally Efficient Decimators and InterpolatorsDecimators and Interpolators
• Note that still has a symmetric impulse response, whereas is the mirror image of
• These relations can be made use of in developing a computationally efficient realization using only 6 multipliers and 11 two-input adders as shown on the next slide
)(1 zE)(0 zE
)(2 zE
![Page 38: Digital Filter Structures - Computer Action Teamweb.cecs.pdx.edu/~mperkows/CAPSTONES/DSP1/ELG6163_FIR.pdfTitle Digital Filter Structures Author S. K. Mitra Created Date 1/18/2007 2:13:50](https://reader036.vdocument.in/reader036/viewer/2022062600/5b03ba8f7f8b9a6c0b8cc48f/html5/thumbnails/38.jpg)
38Copyright © 2001, S. K. Mitra
Computationally Efficient Computationally Efficient Decimators and InterpolatorsDecimators and Interpolators
• Factor-of-3 decimator with a linear-phase decimation filter