1 / 45 Designs of Orthogonal Filter Banks and Orthogonal Cosine-Modulated Filter Banks Jie Yan Department of Electrical and Computer Engineering University of Victoria April 16, 2010
2 / 45 OUTLINE 1 INTRODUCTION 2 LS DESIGN OF ORTHOGONAL FILTER BANKS AND WAVELETS 3 MIMINAX DESIGN OF ORTHOGONAL FILTER BANKS AND WAVELETS 4 DESIGN OF ORTHOGONAL COSINE-MODULATED FILTER BANKS 5 CONCLUSIONS AND FUTURE RESEARCH
3 / 45 1. INTRODUCTION A two-channel conjugate quadrature (CQ) filter bank where H 0 (z) = N 1 n=0 h nz n H 1 (z) = z (N 1) H 0 ( z 1 ) G 0 (z) = H 1 ( z) G 1 (z) = H 0 ( z) H ( z) 0 G ( z) 0 H ( z) 1 Analysis Filter Bank G ( z) 1 Synthesis Filter Bank
4 / 45 Two-Channel Orthogonal Filter Banks Perfect reconstruction (PR) condition N 1 2m n=0 h n h n+2m = δ m for m = 0, 1,..., (N 2)/2 Vanishing moment (VM) requirement: A CQ filter has L vanishing moments if N 1 ( 1) n n l h n = 0 for l = 0, 1,..., L 1 n=0
5 / 45 Two-Channel Orthogonal Filter Banks Cont d A least squares (LS) design of CQ lowpass filter H 0 (z) having L VMs minimize subject to: π H 0 (e jω ) 2 dω ω a PR condition and VM requirement The LS problem above can be expressed as minimize subject to: h T Qh N 1 2m n=0 h n h n+2m = δ m for m = 0, 1,..., (N 2)/2 N 1 ( 1) n n l h n = 0 for l = 0, 1,..., L 1 n=0
6 / 45 Two-Channel Orthogonal Filter Banks Cont d A minimax design minimizes the maximum instantaneous power of H 0 (z) over its stopband minimize subject to: maximize H 0 (e jω ) ω a ω π PR condition and VM requirement The minimax problem can be further cast as minimize subject to: η T(ω) h η for ω Ω N 1 2m n=0 h n h n+2m = δ m for m = 0, 1,..., (N 2)/2 N 1 ( 1) n n l h n = 0 for l = 0, 1,..., L 1 n=0
7 / 45 Orthogonal Cosine-Modulated Filter Banks An orthogonal cosine-modulated (OCM) filter bank [ ( π h k (n) = 2h(n) cos k + 1 ) ( n D ) + ( 1) k π ] M 2 2 4 [ ( π f k (n) = 2h(n) cos k + 1 ) ( n D ) ( 1) k π ] M 2 2 4 for 0 k M 1 and 0 n N 1 x(n) H 0 (z) M M F 0 (z) + y(n) H 1 (z) M M F 1 (z)...... H M-1 (z) M M F M-1 (z)
8 / 45 Orthogonal Cosine-Modulated Filter Banks Cont d An M-channel OCM filter bank is uniquely characterized by its prototype filter (PF) The design of the PF of an OCM filter bank can be formulated as minimize subject to: π ω s H 0 (e jω ) 2 dω PR condition As the PF has linear phase, h is symmetrical. The design problem can be reduced to minimize e 2 (ĥ) = ĥt ˆPĥ subject to: a l,n (ĥ) = ĥt ˆQl,n ĥ c n = 0 for 0 n m 1 and 0 l M/2 1 where the design variables are reduced by half to ĥ = [h 0 h 1 h N/2 1 ] T.
9 / 45 Overview and Contribution of the Thesis Overview We have formulated three nonconvex optimization problems LS design of CQ filter banks Minimax design of CQ filter banks Design of OCM filter banks Contribution of the thesis Several improved local design methods for the three problems Several strategies proposed for potentially GLOBAL solutions of the three problems
10 / 45 Global Design Method at a Glance Multiple local solutions exist for a nonconvex problem Algorithms in finding a locally optimal solution are available Start the local design algorithm from a good initial point How do we secure such a good initial point?
11 / 45 2. LS DESIGN OF ORTHOGONAL FILTER BANKS AND WAVELETS A least squares (LS) design of a conjugate quadrature (CQ) filter of length-n with L vanishing moments (VMs) can be cast as minimize subject to: h T Qh N 1 2m n=0 h n h n+2m = δ m for m = 0, 1,..., (N 2)/2 N 1 ( 1) n n l h n = 0 for l = 0, 1,..., L 1 n=0
12 / 45 Local LS Design of CQ Filter Banks An effective direct design method is recently proposed by W.-S. Lu and T. Hinamoto Based on the direct design technique, we develop two local methods Sequential convex-programming (SCP) method Sequential quadratic-programming (SQP) method Both methods produce improved local designs than the direct method
13 / 45 Local LS Design of CQ Filter Banks Cont d Sequential Convex-Programming Method Suppose we are in the kth iteration to compute δ h so that h k+1 = h k + δ h reduces the filter s stopband energy and better satisfies the constraints, then h T k+1qh k+1 = δ T h Qδ h + 2δ T h Qh k + h T k Qh k N 1 N 1 ( 1) n n l (δ h ) n = ( 1) n n l (h k ) n n=0 n=0 N 1 2m n=0 (h k ) n (δ h ) n+2m + N 1 2m n=0 N 1 2m δ m (h k ) n (h k ) n+2m n=0 (h k ) n+2m (δ h ) n
Local LS Design of CQ Filter Banks Cont d With h bounded to be small, the kth iteration assumes the form minimize subject to: δ T h Qδ h + δ T h g k A k δ h = a k Cδ h b By using SVD to remove the equality constraint, the problem is reduced to minimize subject to: x T ˆQx + xtĝ k Ĉx ˆb We modify the problem to make it always feasible as minimize subject to: x T ˆQx + xtĝ k Fx a which is a convex QP problem. 14 / 45
15 / 45 Local LS Design of CQ Filter Banks Cont d Sequential Quadratic-Programming Method The design problem is a general nonlinear optimization problem minimize subject to: f (h) a i (h) = 0 for i = 1, 2,..., p By using the first-order necessary conditions of a local minimizer, the problem can be reduced to minimize subject to: 1 2 δt h W k δ h + δ T h g k A k δ h = a k δ h is small
16 / 45 Local LS Design of CQ Filter Banks Cont d where p W k = 2 hf (h k ) (λ k ) i 2 ha i (h k ) i=1 A k = [ h a 1 (h k ) h a 2 (h k ) h a p (h k ) ] T g k = h f (h k ) a k = [ a 1 (h k ) a 2 (h k ) a p (h k ) ] T (13a) (13b) (13c) (13d) By removing the equality constraint using the SVD or QR decomposition, the problem assumes the form of a QP problem. Once the minimizer δ h is found, the next iterate is set to h k+1 = h k + δ h, λ k+1 = (A k A T k ) 1 A k (W k δ h + g k )
17 / 45 Global LS Design of Low-Order CQ Filter Banks The LS design problem is a polynomial optimization problem (POP) Two recent breakthroughs in solving POPs Global solutions of POPs are made available by Lasserre s method Sparse SDP relaxation is proposed for global solutions of POPs of relatively larger scales MATLAB toolbox SparsePOP and GloptiPoly can be used to find global solutions of POPs, but only for POPs of limited sizes
18 / 45 Global LS Design of Low-Order CQ Filter Banks Cont d Example: Design a globally optimal LS CQ filter with N = 6, L = 2 and ω a = 0.56π MATLAB toolbox GloptiPoly and SparsePOP are utilized to produce the globally optimal solution h (6,2) LS = 0.33268098788629 0.80689591454849 0.45986215652386 0.13501431772967 0.08543638600240 0.03522516035714 However, GloptiPoly and SparsePOP fail to work as long as the filter length N is greater than or equal to 18
19 / 45 Global LS Design of High-Order CQ Filter Banks A common pattern shared among globally optimal low-order impulse responses. 0.8 N = 6, L = 2 0.7 N = 8, L = 2 N = 10, L = 2 0.6 0.5 0.4 0.3 0.2 0.1 0 0.1 0.2 0 0.2 0.4 0.6 0.8 1
20 / 45 Global LS Design of High-Order CQ Filter Banks Cont d h 6 : Globally optimal impulse response when N = 6 h zp 8 : Impulse response generated by zero-padding h 6 h 8 : Globally optimal impulse response when N = 8 0.8 h 6 (N=6, L=2) zp 0.7 h 8 h 0.6 8 (N=8, L=2) 0.5 0.4 0.3 0.2 0.1 0 0.1 0.2 0 0.2 0.4 0.6 0.8 1 Generate initial point by zero-padding!
21 / 45 Global LS Design of High-Order CQ Filter Banks Cont d Global design strategy in brief: 1 Design a globally optimal CQ filter of short length, say 4, using e.g. GloptiPoly 2 Generating an impulse response for higher order design by zero-padding 3 Apply the SCP or SQP method with the zero-padded impulse response as the initial point to obtain the optimal impulse response of higher order 4 Follow this concept in an iterative way, until desired filter length is reached
22 / 45 Global LS Design of High-Order CQ Filter Banks Cont d The designs obtained are quite likely to be globally optimal because: 1 Zero-padded initial point sufficiently close to the global minimizer. 2 The local design methods are known to converge to a nearby minimizer.
23 / 45 Design Examples Potentially globally optimal design of an LS CQ filter with N = 96, L = 3 and ω = 0.56π 0 20 40 60 80 100 120 0 0.2 0.4 0.6 0.8 1 Normalized frequency
24 / 45 Design Examples Cont d Comparisons Global design Global design based on SCP Energy in stopband 1.18103e-9 Largest eq. error 1e-14 Local design Local design based on SCP Energy in stopband 3.15564e-9 Largest eq. error 1e-14
25 / 45 Design Examples Cont d Zero-pole plots Imaginary Part 1 0.8 0.6 0.4 0.2 0 0.2 0.4 0.6 0.8 1 4 95 1 0.5 0 0.5 1 Real Part Global design Imaginary Part 1 0.5 0 3 95 0.5 1 1.5 1 0.5 0 0.5 1 1.5 Real Part Local design The globally optimal LS CQ filter possesses minimum phase
26 / 45 3. MIMINAX DESIGN OF ORTHOGONAL FILTER BANKS AND WAVELETS A minimax design of a conjugate quadrature (CQ) filter of length-n with L vanishing moments (VMs) can be cast as minimize subject to: η T(ω) h η for ω Ω N 1 2m n=0 h n h n+2m = δ m for m = 0, 1,..., (N 2)/2 N 1 ( 1) n n l h n = 0 for l = 0, 1,..., L 1 n=0
27 / 45 Local Minimax Design of CQ Filter Banks Like the LS design, an effective direct design method is recently proposed by W.-S. Lu and T. Hinamoto Based on the direct design technique, we develop an improved method named the SCP-GN method The SCP-GN method can achieve convergence at a small tolerance ε, by implementing two techniques 1 Constructing Ω by locating magnitude-response peaks 2 A Gauss-Newton method with adaptively controlled weights
Local Minimax Design of CQ Filter Banks Cont d In the kth iteration, the problem assumes the form minimize subject to: η T(ω)(h k + δ h) η for ω Ω A kδ h = a k Cδ h b By using SVD of matrix A k to remove the equality constraints, the problem can be reduced to minimize subject to: η T k(ω)x + e k(ω) η for ω Ω Ĉx ˆb As a technical remedy to make the above problem to be always feasible, we modify the problem as minimize subject to: η T k(ω)x + e k(ω) η for ω Ω Fx a which is a second-order cone programming (SOCP) problem for which efficient solvers such as SeDuMi exist. 28 / 45
Global Minimax Design of Low-Order CQ Filter Banks Example: Design a globally optimal minimax CQ filter with N = 4, L = 1 and ω a = 0.56π Since the Minimax design problem is a POP, GloptiPoly and SparsePOP can be used to produce the globally optimal solution 0.48296282173531 h (4,1) minimax = 0.83651623138234 0.22414405492402 0.12940935473280 However, GloptiPoly and SparsePOP fail to work as long as the filter length N is greater than or equal to 6 29 / 45
30 / 45 Global Minimax Design of High-Order CQ Filter Banks Method 1 Globally optimal minimax impulse responses appear to exhibit a pattern similar to that in the LS case Thus, we proposed method 1 in spirit similar to that utilized in the global LS designs by passing the zero-padded impulse response as the initial point for the SCP-GN local method in each round of iteration Method 2 We simply pass the impulse response of the globally optimal LS filter as an initial point for the SCP-GN method to design an optimal minimax filter with the same design specifications
31 / 45 Design Examples Potentially globally optimal design of a minimax CQ filter with N = 96, L = 3 and ω = 0.56π 0 20 40 60 80 100 0 0.2 0.4 0.6 0.8 1 Normalized frequency
32 / 45 Design Examples Cont d Comparisons Global design Maximum instantaneous energy in stopband Largest eq. error Local design Maximum instantaneous energy in stopband Largest eq. error Global design based on Method 1 6.75750e-9 <1e-15 Local design based on SCP-GN 1.81165e-8 2.9e-14
33 / 45 Design Examples Cont d Zero-pole plots 1 0.8 1 0.6 0.4 0.5 Imaginary Part 0.2 0 0.2 0.4 4 95 Imaginary Part 0 4 95 0.6 0.5 0.8 1 1 0.5 0 0.5 1 Real Part 1 1 0.5 0 0.5 1 1.5 Real Part Global design Local design The globally optimal minimax CQ filter possesses minimum phase
34 / 45 4. DESIGN OF ORTHOGONAL COSINE-MODULATED FILTER BANKS We have formulated the design of the prototype filter (PF) of an orthogonal cosine-modulated (OCM) filter bank as minimize e 2 (ĥ) = ĥt ˆPĥ subject to: a l,n (ĥ) = ĥt ˆQl,n ĥ c n = 0 for 0 n m 1 and 0 l M/2 1
35 / 45 Local Design of OCM Filter Banks We improved an effective direct design method proposed by W.-S. Lu, T. Saramäki and R. Bregović Gauss-Newton method with adaptively controlled weights was applied for the algorithm to converge to a highly accurate solution
Local Design of OCM Filter Banks Cont d Suppose we are in the kth iteration to compute δ so that ĥk+1 = ĥk + δ reduces the PF s stopband energy and better satisfies the PR conditions. Then, ĥ T k+1ˆpĥk+1 = δ T ˆPδ + 2δ T ˆPĥ k + ĥt k ˆPĥk And the kth iteration assumes the form a l,n (ĥk + δ) a l,n (ĥk) + g T l,n (ĥk)δ = 0 for 0 n m 1 and 0 l M/2 1 minimize subject to: δ T ˆPδ + δ T b k G k δ = a k δ is small 36 / 45
37 / 45 Local Design of OCM Filter Banks Cont d The equality constraint can be eliminated via SVD of G k = UΣV as Thus, the problem can be cast as δ = V e φ + δ s (18) minimize subject to: φ T Pk φ + φ T bk φ is small The Gauss-Newton technique with adaptively controlled weights is used as a post-processing step to achieve convergence at a small tolerance.
Global Design of Low-Order OCM Filter Banks Example: Design a globally optimal OCM filter bank with M = 2, m = 1 and ρ = 1 GloptiPoly and SparsePOP can be used to produce the globally optimal solution h (2,1) = 0.235923416966353 0.440840267366581 0.440840267366581 0.235923416966353 The software was found to work only for the following cases: a) M = 2, 1 m 5; b) M = 4, 1 m 3; c) M = 6, m = 1; d) M = 8, m = 1. 38 / 45
39 / 45 Global Design of High-Order OCM Filter Banks Two observations: 1 For a fixed M, the impulse responses with different m exhibit a similar pattern and are close to each other 2 For m = 1, the impulse responses with different M also exhibit a similar shape. 0.45 0.4 0.35 0.3 m=1,m=2 m=2,m=2 m=3,m=2 m=4,m=2 m=1,m=4 m=2,m=4 0.25 0.2 0.15 0.1 0.05 0 0.05 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Global Design of High-Order OCM Filter Banks Cont d 0.35 0.3 m=1,m=4 zp h 0 m=2,m=4 0.5 0.45 m=1,m=2 int h 0 m=1,m=4 0.25 0.4 0.2 0.35 0.15 0.3 0.1 0.25 0.05 0.2 0 0.15 0.05 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0.1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 Effect of zero-padding when M = 4 Effect of linear interpolation when m = 1 40 / 45
41 / 45 Global Design of High-Order OCM Filter Banks Cont d An improvement in initial point when m = 1, by downshifting h int 0 by a constant value d computed using the Gauss-Newton method with adaptively controlled weights. 0.5 0.45 h 0 int h 0 h 0.4 d 0.35 0.3 0.25 0.2 0.15 0.1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
42 / 45 Global Design of High-Order OCM Filter Banks Cont d An order-recursive algorithm in brief 1 Obtaining a low-order global design; 2 Using zero-padding/linear interpolation in conjuction of the G-N method with adaptively controlled weights of the impulse response to produce a desirable initial point for PF of slightly increased order, and carrying out the design by a locally optimal method; 3 Repeating step 2 until the filter order reaches the targeted value.
43 / 45 Design Examples Design of an OCM filter bank with m = 20, M = 4 and ρ = 1. Shown below are the impulse responses of the PF from global and local design, respectively. 20 0 20 40 60 80 100 120 140 160 180 200 0 0.2 0.4 0.6 0.8 1 Normalized frequency 20 0 20 40 60 80 100 120 140 160 0 0.2 0.4 0.6 0.8 1 Normalized frequency Global design Local design
44 / 45 Design Examples Cont d Performance comparison for OCM filter banks with m = 20, M = 4 and ρ = 1 Global design Local design Energy in stopband 8.226e-13 6.585e-10 Largest eq. error 1.839e-15 2.297e-10 By comparing the OCM filter banks reported in the literature, the OCM filter bank designed using our proposed algorithm offers the BEST performance, because it is a globally optimal design.
45 / 45 5. CONCLUSIONS AND FUTURE RESEARCH We have investigated three design problems, 1 LS design of orthogonal filter banks and wavelets 2 Minimax design of orthogonal filter banks and wavelet 3 Design of OCM filter banks Improved local design methods for the three problems Several strategies proposed for GLOBAL designs of the three design scenarios Future research Theoretical proof of the global optimality of our proposed method