450 likes | 697 Views
Digital Signal Processing Lecture-4 Fall 2009 Polyphase Decomposition and Mutirate Signal Processing Chap.-4. Hassan Bhatti. Motivation. Up- and down sampling combined with filtering are the usual operations in multirate systems. Polypahse approach will yield simple implementations.
E N D
Digital Signal Processing Lecture-4 Fall 2009Polyphase Decomposition and Mutirate Signal Processing Chap.-4 Hassan Bhatti
Motivation • Up- and down sampling combined with filtering are the usual operations in multirate systems. • Polypahse approach will yield simple implementations
Outlines • Two basic multirate operations • Polyphase interpolator – upsampling followed by a filter • Polyphase decimatator – a filter followed by a decimator • Two-channel filter banks • Perfect reconstruction condition • Quadrature mirror filter (QMF) filter banks • Design of two-channel filter banks with PR • Multiple-channel filter banks • Tree- structured filter banks • Octave-band filter banks
Basic Multirate Operations • Decimation and interpolation • Z-domain and Frequency domain analysis of up-and downsampled version of a signal • Polyphase decomposition • Noble Identities
Decimation and Interpolation • Decimation---down-sampling N x(n)
Decimation and Interpolation • Decimation---down-sampling N
Decimation and Interpolation • Decimation---down-sampling N y(m)
Decimation and Interpolation • Interpolation --- up-sampling N
Decimation and Interpolation • Interpolation --- up-sampling N
Decimation and Interpolation • Interpolation --- up-sampling N
H(z) N N G(z) Decimation and Interpolation A typical building block of multirate filter bank We want to know the relationships between the above signals
Decimation and Interpolation Up-sampling N
Decimation and Interpolation Upsampling when N=2
Decimation and Interpolation Downsampling followed by upsampling N N as hence
Decimation and Interpolation Downsampling followed by upsampling N N
Decimation and Interpolation Downsampling followed by upsampling N N image spectra original spectrum
Decimation and Interpolation Downsampling followed by upsampling image spectra original spectrum
Decimation and Interpolation Downsampling followed by upsampling when N=2 Image spectra
Decimation and Interpolation Downsampling N We know upsampling We know downsampling +upsampling We can get downsampling
Decimation and Interpolation Downsampling N Example, N=2
Decimation and Interpolation Downsampling when N=2 Image spectra
Decimation and Interpolation Downsampling when N=2 Image spectra
Decimation and Interpolation Downsampling when N=2 Image spectra
H(z) N N G(z) Decimation and Interpolation A typical building block of multirate filter bank We want to know the relationships between the above signals
Polyphase Decomposition • Polyphase decomposition is the decomposition of a sequence x(n) into sub-sequences x(mN+i) • There are four types of polyphase decomposition. Type-1 [ 0,1,2,3,4,5,6,7,8,9,10,11] M=3 [ 0,3,6,9] [ 1,4,7,10] [ 2,5,8,11]
Polyphase Decomposition Type-1 M M + T 2T M
Polyphase Decomposition Type-2 [ 0,1,2,3,4,5,6,7,8,9,10,11] M=3 [ 2,5,8,11] [ 1,4,7,10] [ 0,3,6,9]
Polyphase Decomposition Type-3: we want to have hence Type-3 is not very straightforward as there is a casualty problem.
Nobel Identities Identity I N G(z) G(zN) N
Nobel Identities Identity II N G(zN) G(z) N
H(z) N N G(z) Decimation and Interpolation—polyphase implementation A typical building block of multirate filter bank We want to know if there is an efficient way to implement the above system
Polyphase Interpolator N G(z) x(n) y(m) v(n) Multiplications with zeros are involved
Polyphase Intepolator • We decompose the filter into polyphase components (Type-1): • In z-domain:
Polyphase Decomposition (from the notes last week) Type-1 M M + T 2T M
N G(z) x(n) y(m) v(n) Polyphase Intepolator N x(n) y(m) v(n)
N G(z) x(n) y(m) v(n) N N N Polyphase Intepolator Using the second Noble identity: x(n) y(m)
N G(z) x(n) y(m) v(n) N x(n) N y(m) N Polyphase Intepolator For input signal of length M, and G(z) of length L, convolution of v(n) (of length NM) and g(n) (of length L) requires NML mutiplications
N x(n) N y(m) N Polyphase Intepolator • Each branch has a convolution of y(m) (of length M) with Gl(z) (of length L/N). ML/N multiplications are required; • Hence for N branches, ML multiplications are required in total; Computation is greatly reduced
Polyphase Intepolator • example: M=1024, N=2,L=64, MNL=128k multiplications are required for convolution • when using polyphase approach, only ML=64k multiplications are required.
Polyphase Decimator – polyphase implementation H(z) N y(m) x(n) v(n) The number of multiplications is: ML
Polyphase Decimator Let i=jN+k let
Polyphase Decimator N y(m) x(n) N N The number of multiplications is: N (M/N)(L/N) =ML/N. Also reduced a lot.