1 / 40

Digital Signal Processing II Chapter 8: Modulated Filter Banks

Digital Signal Processing II Chapter 8: Modulated Filter Banks. Marc Moonen Dept. E.E./ESAT, K.U.Leuven marc.moonen@esat.kuleuven.be www.esat.kuleuven.be/scd/. Part-II : Filter Banks. : Preliminaries Filter bank set-up and applications

faith
Download Presentation

Digital Signal Processing II Chapter 8: Modulated Filter Banks

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Digital Signal Processing IIChapter 8: Modulated Filter Banks Marc Moonen Dept. E.E./ESAT, K.U.Leuven marc.moonen@esat.kuleuven.be www.esat.kuleuven.be/scd/

  2. Part-II : Filter Banks : Preliminaries • Filter bank set-up and applications • `Perfect reconstruction’ problem + 1st example (DFT/IDFT) • Multi-rate systems review (10 slides) : Maximally decimated FBs • Perfect reconstruction filter banks (PR FBs) • Paraunitary PR FBs : Modulated FBs • Maximally decimated DFT-modulated FBs • Oversampled DFT-modulated FBs : Special Topics • Cosine-modulated FBs • Non-uniform FBs & Wavelets • Frequency domain filtering Chapter-6 Chapter-7 Chapter-8 Chapter-9 Version 2010-2011Chapter-8: Modulated Filter Banks

  3. 3 3 F0(z) subband processing H0(z) OUT IN 3 3 F1(z) subband processing H1(z) + 3 3 F2(z) subband processing H2(z) 3 3 F3(z) subband processing H3(z) Refresh (1) General `subband processing’ set-up (Chapter-6) : PS: subband processing ignored in filter bank design synthesis bank analysis bank downsampling/decimation upsampling/expansion Version 2010-2011Chapter-8: Modulated Filter Banks

  4. u[k-3] u[k] 4 4 4 4 + 4 4 4 4 Refresh (2) Two design issues : - filter specifications, e.g. stopband attenuation, passband ripple, transition band, etc. (for each (analysis) filter!) - perfect reconstruction property (Chapter-6). PS: we are now still considering maximally decimated FB’s, i.e. Version 2010-2011Chapter-8: Modulated Filter Banks

  5. Introduction -All design procedures so far involve monitoring of characteristics (passband ripple, stopband suppression,…) of all (analysis) filters, which may be tedious. -Design complexity may be reduced through usage of `uniform’ and `modulated’ filter banks. • DFT-modulated FBs (this Chapter) • Cosine-modulated FBs (next Chapter) Version 2010-2011Chapter-8: Modulated Filter Banks

  6. H0 H1 H2 H3 uniform H0 H3 non-uniform H1 H2 H0(z) IN H1(z) H2(z) H3(z) Introduction Uniform versus non-uniform (analysis) filter bank: • N-channel uniform FB: i.e. frequency responses are uniformly shifted over the unit circle Ho(z)= `prototype’ filter (=only filter that has to be designed) Time domain equivalent is: • non-uniform = everything that is not uniform e.g. for speech & audio applications (cfr. human hearing) example: wavelet filter banks (next Chapter) Version 2010-2011Chapter-8: Modulated Filter Banks

  7. H0(z) u[k] H1(z) H2(z) H3(z) i.e. Maximally Decimated DFT-Modulated FBs Uniform filter banks can be realized cheaply based on polyphase decompositions + DFT(FFT) (hence name `DFT-modulated FB) 1. Analysis FB If (polyphase decomposition) then Version 2010-2011Chapter-8: Modulated Filter Banks

  8. i.e. Maximally Decimated DFT-Modulated FBs where F is NxN DFT-matrix(and `*’ is complex conjugate) This means that filtering with the Hi’s can be implemented by first filtering with polyphase components and then DFT Version 2010-2011Chapter-8: Modulated Filter Banks

  9. u[k] i.e. Maximally Decimated DFT-Modulated FBs conclusion: economy in… • implementation complexity (for FIR filters): N filters for the price of 1, plus DFT (=FFT) ! • design complexity: Design `prototype’ Ho(z), then other Hi(z)’s are automatically `co-designed’ (same passband ripple, etc…) ! Version 2010-2011Chapter-8: Modulated Filter Banks

  10. u[k] Ho(z) H1(z) Maximally Decimated DFT-Modulated FBs • Special case: DFT-filter bank, if all Ei(z)=1 Version 2010-2011Chapter-8: Modulated Filter Banks

  11. u[k] Ho(z) H1(z) Maximally Decimated DFT-Modulated FBs • PS: with F instead of F* (as in Chapter-6), only filter ordering is changed Version 2010-2011Chapter-8: Modulated Filter Banks

  12. u[k] u[k] 4 4 4 4 4 = 4 4 4 Maximally Decimated DFT-Modulated FBs • DFT-modulated analysis FB + maximal decimation (M=N) = efficient realization ! Version 2010-2011Chapter-8: Modulated Filter Banks

  13. y[k] + + + Maximally Decimated DFT-Modulated FBs 2. Synthesis FB phase shift added for convenience Version 2010-2011Chapter-8: Modulated Filter Banks

  14. i.e. Maximally Decimated DFT-Modulated FBs where F is NxN DFT-matrix Version 2010-2011Chapter-8: Modulated Filter Banks

  15. + + + i.e. Maximally Decimated DFT-Modulated FBs y[k] Version 2010-2011Chapter-8: Modulated Filter Banks

  16. 4 4 4 4 4 4 4 + + + + + + 4 Maximally Decimated DFT-Modulated FBs • Expansion (M=N) + DFT-modulated synthesis FB : y[k] = = efficient realization ! y[k] Version 2010-2011Chapter-8: Modulated Filter Banks

  17. 4 u[k] 4 4 4 y[k] + + + 4 4 4 4 Maximally Decimated DFT-Modulated FBs How to achieve Perfect Reconstruction (PR) with maximally decimated DFT-modulated FBs? i.e. synthesis bank polyphase components are obtained by inverting analysis bank polyphase components Version 2010-2011Chapter-8: Modulated Filter Banks

  18. 4 u[k] 4 4 4 y[k] + + + 4 4 4 4 Maximally Decimated DFT-Modulated FBs Design Procedure : 1. Design prototype analysis filter Ho(z) (see Chapter-3). 2. This determines Ei(z) (=polyphase components). 3. Assuming all Ei(z) can be inverted (?), choose synthesis filters Version 2010-2011Chapter-8: Modulated Filter Banks

  19. Maximally Decimated DFT-Modulated FBs • Will consider only FIR prototype analysis filters, leading to simple polyphase decompositions (see Chapter-6). • However, FIR Ei(z)’s generally still lead to IIR Ri(z)’s, where stability is a concern… Ri(z) ’s are stable only if Ei(z)’s have stable zeros (i.e. are `minimum-phase filters’). Example: LPC lattice filters with all |ki|<1 (see Chapter-4). The design of such minimum phase FIR filters is (significantly) more difficult.. • FIR Ri(z)’s (=guaranteed stability) are only obtained with trivial choices for the Ei(z)’s …namely Ei(z)’s with only 1 non-zero impulse response parameter. E(z) is then unimodular (see Chapter-7). Examples: see next slide. Version 2010-2011Chapter-8: Modulated Filter Banks

  20. Maximally Decimated DFT-Modulated FBs • Simple example (1) is , which leads to IDFT/DFT bank (Chapter-6) i.e. Fl(z) has coefficients of Hl(z), but complex conjugated and in reverse order (hence same magnitude response) (remember this?!) • Simple example (2) is , where wi’s are constants, which leads to `windowed’ IDFT/DFT bank, a.k.a. `short-time Fourier transform’ (see Chapter-9) Version 2010-2011Chapter-8: Modulated Filter Banks

  21. Maximally Decimated DFT-Modulated FBs Question (try to answer): Can we have paraunitary FBs here (=desirable property) ? When is maximally decimated DFT-modulated FB at the same time - PR - FIR (both analysis & synthesis) - Paraunitary ? Hint: E(z) is paraunitary only if the Ei(z)’s are all-pass filters. An FIR all-pass filter takes a trivial form, e.g. Ei(z)=1 or Ei(z)=z^{-d} Version 2010-2011Chapter-8: Modulated Filter Banks

  22. Maximally Decimated DFT-Modulated FBs • Bad news: It is seen that the maximally decimated IDFT/DFT filter bank (or trivial modifications thereof) is the only possible maximally decimated DFT- modulated FB that is at the same time... - PR - FIR (all analysis+synthesis filters) - Paraunitary • Good news: • Cosine-modulatedPR FIR FB’s (Chapter-9) • OversampledPR FIR DFT-modulated FB’s (read on) Version 2010-2011Chapter-8: Modulated Filter Banks

  23. Oversampled PR Filter Banks • So far have considered maximal decimation (M=N), where aliasing makes PR design non-trivial. • With downsampling factor (N) smaller than the number of channels (M), aliasing is expected to become a smaller problem, possibly negligible if N<<M. • Still, PR theory (with perfect alias cancellation) is not necessarily simpler ! • Will not consider PR theory as such here, only give some examples of oversampled DFT-modulated FBs that are PR/FIR/paraunitary (!) Version 2010-2011Chapter-8: Modulated Filter Banks

  24. u[k-3] u[k] 4 4 4 4 + 4 4 4 4 Oversampled PR Filter Banks • Starting point is(see Chapter-7): (delta=0 for conciseness here) where E(z) and R(z) are NxN matrices (cfr maximal decimation) • What if we try other dimensions for E(z) and R(z)…?? Version 2010-2011Chapter-8: Modulated Filter Banks

  25. u[k] 4 4 4 4 + 4 4 4 4 Oversampled PR Filter Banks ! • A more general case is : where E(z) is now MxN (`tall-thin’) and R(z) is NxM (`short-fat’) while still guarantees PR ! u[k-3] N=4 decimation M=6 channels Version 2010-2011Chapter-8: Modulated Filter Banks

  26. Oversampled PR Filter Banks • The PR condition appears to be a `milder’ requirement if M>N for instance for M=2N, we have (where Ei and Ri are NxN matrices) which does not necessarily imply that meaning that inverses may be avoided, creating possibilities for (great) DFT-modulated FBs, which can (see below) be PR/FIR/paraunitary • In the sequel, will give 2 examples of oversampled DFT-modulated FBs Version 2010-2011Chapter-8: Modulated Filter Banks

  27. Should not try to understand this… Oversampled DFT-Modulated FBs Example-1 : # channels M = 8 Ho(z),H1(z),…,H7(z) decimation N = 4 prototype analysis filter Ho(z) will consider N’-fold polyphase expansion, with Version 2010-2011Chapter-8: Modulated Filter Banks

  28. u[k] Oversampled DFT-Modulated FBs In general, it is proved that the M-channel DFT-modulated (analysis) filter bank can be realized based on an M-point DFT cascaded with an MxN `polyphase matrix’ B, which contains the (N’-fold) polyphase components of the prototype Ho(z) ps: note that if M=N, then N’=N, and then B is a diagonal matrix (cfr. supra) Example-1 (continued): N=4 decimation M=8 channels Convince yourself that this is indeed correct.. (or see next slide) Version 2010-2011Chapter-8: Modulated Filter Banks

  29. u[k] Oversampled DFT-Modulated FBs Proof is simple: Version 2010-2011Chapter-8: Modulated Filter Banks

  30. 4 u[k] 4 4 4 Oversampled DFT-Modulated FBs -With (N=) 4-fold decimation, this is… Version 2010-2011Chapter-8: Modulated Filter Banks

  31. u[k] 4 4 4 4 + 4 4 4 4 u[k-3] Oversampled DFT-Modulated FBs • Perfect Reconstruction (PR) can now be obtained based on an E(z) that is FIR and paraunitary : • If E(z )=F*.B(z) is chosen to be paraunitary, then PR is obtained with R(z)=B~(z).F (=NxM) (=DFT-modulated synthesis bank). • E(z) is paraunitary only if B(z) is paraunitary. So how can we make B(z) paraunitary ? Version 2010-2011Chapter-8: Modulated Filter Banks

  32. Oversampled DFT-Modulated FBs Example 1 (continued) : • From the structure of B(z) It follows that B(z) is paraunitary if and only if (for k=0,1,2,3) are power complementary i.e. form a lossless 1-input/2-output system (explain!) • For 1-input/2-output power complementary FIR systems, see Chapter-6 on FIR lossless lattices realizations (!)… Version 2010-2011Chapter-8: Modulated Filter Banks

  33. u[k] 4 : : 4 : Oversampled DFT-Modulated FBs Lossless 1-in/2-out • Design Procedure: Optimize parameters (=angles) of (4) FIR lossless lattices (defining polyphase components of Ho(z) ) such that Ho(z) satisfies specifications. p.30 = Version 2010-2011Chapter-8: Modulated Filter Banks

  34. Oversampled DFT-Modulated FBs • Result = oversampled DFT-modulated FB (M=8, N=4), that is PR/FIR/paraunitary !! All great properties combined in one design !! • PS: With 2-fold oversampling (M/N=2 in example-1), paraunitary design is based on 1-input/2-output lossless systems (see page 32-33). In general, with D-fold oversampling (for D=integer), paraunitary design will be based on 1-input/D-output lossless systems (see also Chapter-3 on multi-channel FIR lossless lattices). With maximal decimation (D=1), paraunitary design will then be based on 1-input/1-output lossless systems, i.e. all-pass (polyphase) filters, which in the FIR case can only take trivial forms (=page 21-22) ! Version 2010-2011Chapter-8: Modulated Filter Banks

  35. Should not try to understand this… Oversampled DFT-Modulated FBs Example-2 (non-integer oversampling) : # channels M = 6 Ho(z),H1(z),…,H5(z) decimation N = 4 prototype analysis filter Ho(z) will consider N’-fold polyphase expansion, with Version 2010-2011Chapter-8: Modulated Filter Banks

  36. u[k] Oversampled DFT-Modulated FBs DFT modulated (analysis) filter bank can be realized based on an M-point IDFT cascaded with an MxN polyphase matrix B, which contains the (N’-fold) polyphase components of the prototype Ho(z) Version 2010-2011Chapter-8: Modulated Filter Banks Convince yourself that this is indeed correct.. (or see next slide)

  37. u[k] Oversampled DFT-Modulated FBs Proof is simple: Version 2010-2011Chapter-8: Modulated Filter Banks

  38. u[k] 4 4 4 4 Oversampled DFT-Modulated FBs -With (N=) 4-fold decimation, this is… Version 2010-2011Chapter-8: Modulated Filter Banks

  39. Oversampled DFT-Modulated FBs Perfect Reconstruction by paraunitariness? - E(z) paraunitary iff B(z) paraunitary - B(z) is paraunitary if and only if submatrices are paraunitary (explain!) Hence paraunitary design based on (two) 2-input/3-output lossless systems. Such systems can again be FIR, then parameterized and optimized. Details skipped, but doable! Version 2010-2011Chapter-8: Modulated Filter Banks

  40. Conclusions • Uniform DFT-modulated filter banks are great: Economy in design- and implementation complexity • Maximally decimated DFT-modulated FBs: Sounds great, but no PR/FIR design flexibility  - Oversampled DFT-modulated FBs: Oversampling provides additional design flexibility, not available in maximally decimated case. Hence can have it all at once : PR/FIR/paraunitary!  Version 2010-2011Chapter-8: Modulated Filter Banks

More Related