Paper Menu >>
Journal Menu >>
Journal of Signal and Information Processing, 20 10 , 1, 35 -43 doi:10.4236/jsip.2010.11004 Published Online November 2010 (http://www.SciRP.org/journal/jsip) Copyright © 2010 SciRes. JSIP 35 Design of M-Band NPR Cosine-Modulated Filterbank Using IFIR Technique Ram Kumar So ni1, Alok Jain2, Rajiv Saxena3 1Department of Electronics & Communication Engineering, Samrat Ashok Technological Institute (Polytechnic), Vidisha, India, 2Department of Electronics & Instrumentation Engineering, Samrat Ashok Technological Institute, Vidisha, India, 3Department of Electronics & Communication Engineering, Jaypee University of Engineering & Technology, Raghogarh, Guna, India. Email: soniram04@rediffmail.com, alokjain6@rediffmail.com, rajiv.saxena@jiet.ac.in Received September 28th, 2010, revised November 1st, 2010; accepted November 3rd, 2010. ABSTRACT This paper presents a design of M-band near perfect reconstructed (NPR) cosine-modulated filterbank (CMFB). The prototype filter is formulated as an interpolated finite impulse response (IFIR) filter. Suitable stretched factor is used to provide considerable reduction in computation cost as well as minimum value of error parameters. Further minimiza- tion in errors has been achieved by applying a linear gradient optimization technique. The IFIR approach provides reduction in stopband energy as side-lobe-fall-off-rate (SLFOR) in magnitude response of the prototype filter is im- proved. Design examples have been included to illustrate the effectiveness of the proposed technique over the existing work. Keywords: CMFB, IFIR, NPR, Window 1. Introduction Cosine modulated filterbank find wide applications in many areas of signal processing such as data compres- sion of speech, audio and video signals, denoising, fea- ture detection and extraction and adaptive signal proc- essing [1-5]. In this approach all the analysis filters k H z and synthesis filters k F z are designed by the cosine modulation of prototype lowpass filter as shown in Figure 1. As a result, the design of filterbank becomes simple, more realizable and the computation cost reduces to the cost of one prototype filter plus modulation over- head. This design technique has been extensively studied by many prominent authors [3-7]. The use of finite im- pulse response (FIR) filter is oftenly preferred over the infinite impulse response (IIR) filter in the design of prototype lowpass filter. Such type of filter has some desirable features like linear-phase, inherent-stability, negligible-quantization noise, etc. [4]. However, the computational requirements of a FIR filter are usually greater than an IIR filter. This is essentially true in CMFB design for high stopband attenuation and narrow transition band width. The high computational require- ment of FIR filter can be reduce by exploiting the redun- dancy in filter coefficients. Neuvo et al. [8] have devel- oped a technique for FIR filter design in which the re- quired specifications are obtained by the cascade of FIR sections. The overall structure is known as interpolated finite impulse response (IFIR) filter. This approach re- quires less multipliers than the conventional single stage FIR implementation. In the proposed work, an NPR type cosine modulated filterbank is designed. The required prototype filter is formulated by IFIR technique. The comparison has been made with earlier reported work of Furado et al. [9] and Kha et al. [10]. In first one the frequency response masking (FRM) approach has been used for the design of prototype filter of a filterbank. Similarly, in second case prototype filter is designed using standard FIR approach and a convex optimization technique is used to minimize the associated errors of filterbank. The contribution of this paper is summarized below: 1) The prototype filter is designed using IFIR ap- proach to provide reduction in computation complexity at high stopband attenuation and narrow transition band width. 2) Most popular variable Kaiser window function is used to design the prototype filter. Design of M-Band NPR Cosine-Modulated Filterbank Using IFIR Technique Copyright © 2010 SciRes. JSIP 36 Figure 1. M-band cosine modulated filterbank. 3) NPR type cosine modulation approach is used for the design of filterbank. It avoids the computation of large matrix sets of perfect reconstruction (PR) condition and thereby further reduces the computational burden during the implementation. 4) The reconstructed output of NPR system is not ex- act replica of input due to reconstruction and aliasing error. A linear gradient optimization is applied in which the cutoff frequency is changed in order to approximate the power complementary property. 5) Examples are included to evaluate the performance of the proposed technique in terms of error parameters, computational complexity, group delay and stopband energy. This paper is organized as follow: The IFIR filter with their performance parameters are described in Section 2. This section is also included the design procedure of prototype lowpass filter using proposed IFIR. Section 3 covered the description of cosine modulation approach for filterbank design and parameters related with compu- tational complexity. The algorithm of linear gradient optimization technique with required objective function is provided in Section 4. Section 5 included design ex- amples of eight, sixteen and thirty two band NPR CMFB. The performance of proposed work with earlier reported publications are discussed in Section 6. Finally conclud- ing remarks are made in Section 7. 2. Design of Prototype Filter The impulse response coefficients of a causal Nth-order linear phase FIR filter p[n] using window technique is given by [4]: pnw nhn (1) where, h[n] is the impulse response of the ideal lowpass filter and is expressed as: c sin 0.5 0.5 nN hn nN (2) where, c is the cutoff frequency of the ideal low- pass filter and wn is the Kaiser window function. For the given value of passband p and stopband s frequencies the filter order is calculated by [4]: 7.95 14.36 As N (3) where, 2 sp is the normalized transi- tion bandwidth, As is the stopband attenuation and 2 csp . In the design of narrowband FIR lowpass filter the fil- ter order is very high consequently the required number of adders and multipliers are also high. Therefore, in pro- posed work an IFIR technique is used for the design of prototype lowpass filter. This technique is computation- ally more efficient than standard FIR filter design tech- nique. In IFIR technique the final prototype lowpass fil- ter can be obtained by Eq. (4) [11] L PzGz Iz (4) where, L Gz is the up-sampled version of model filter Gzwith the stretching factor of L and I(z) is the inter- polator filter. In the proposed work, both filters are de- signed by using Kaiser window function. Input parame- ters like stopband attenuation, passband ripple are taken same as desire in final prototype IFIR filter. Only the passband and stopband frequencies are different. For the model filter the values are relaxed by factor of L, i.e. , Design of M-Band NPR Cosine-Modulated Filterbank Using IFIR Technique Copyright © 2010 SciRes. JSIP 37 p L and s L , respectively. Similarly, for the in- terpolator filter the value of passband and stopband fre- quencies are p and 2 s L , respectively. The reduction in computation complexity provided by IFIR filter is very much depends upon the stretch factor. There is an optimum value of stretch factor L given in (5), which provides upto 90% reduction in computation at the cost of other performance parameters [11]. Therefore, the selection of L is depends upon the applications, for which the IFIR filter is designed. The optimum value of L is given as [11]. 2 2 opt ps sp L (5) The computational complexity is stated in terms of number of multipliers needed for implementation. The number of multiplier in case of linear phase FIR structure, polyphase structure and IFIR structure can be calculated as below [4,12] Multiplier in linear phase FIR structure = 2N (6) Multiplier in polyphase structure =2 N (7) Multiplier in IFIR structure = 22 N Ni m (8) Similarly, Adders in linear phase FIR structure = N (9) Adders in polyphase structure = 2N (10) Adders in IFIR structure = i N m N (11) where, N, Nm and Ni are the filter order for standard FIR filter, model filter and interpolator, respectively. Since, the filter order of model filter and interpolator are less because first one is stretched version of desired prototype filter and second one has wide band and large transition width. Therefore, the multiplier (multi) needed in IFIR filter is much less than standard FIR filter and achieved computational reduction is given as: () %Computational reduction F IR IFIR FIR multi multi multi (12) A FIR filter exhibits property of linear phase. There- fore, its phase response is constant and derivatives of phase with respect to frequency called the group delay ( g ) is given by [12]: 2 g dN d (13) Hence, in linear phase FIR filter the groups delay de- pends upon filter order. The proposed technique provides improvement in SLFOR, consequently reduction in stopband energy (Es). For a lowpass prototype filter stopband energy is given as [13] 2 () s s EHd (14) 3. Cosine Modulation Cosine modulation is one of the efficient technique which provides minimum computational cost in filter- bank design. In this technique all the filters of synthesis and analysis sections are obtained by cosine modulation of a linear phase lowpass prototype filter [4,12]. 2cos 211 22 4 2)cos 211 224 k k k k N hn pnkn M N fn pnkn M for 0 1, and 0kM nN (15) Here, M is the number of bands in filterbank and k hnand k f n are the impulse responses of the fil- ters of analysis and synthesis sections, respectively. The relation between the reconstructed output Y(z) and input X(z) is expressible in the z-domain as 12 0 1 M- -j πlM l l Yz TzXzTzXze (16) Here, T0(z) is the overall distortion transfer function and determines the distortion caused by the overall sys- tem for the unaliased component X(z) of the input signal and 12 0 M- jπ lM k lk k Tz FzHze (17) for l = 1,2,,M-1 is called the aliased transfer func- tion and determine how well the aliased components 2 jlM Xze of the input signal are attenuated. To can- cel aliasing and achieve PR, it is required that [4] 0,0, d Tz cznc (18) where, ndis positive integer and 0, l Tz for l = 1,2,3,, M -1 (19) For PR the criteria for the prototype filter are very strict. Fortunately, in many practical applications the PR property can be slightly relaxed, resulting in NPR [4]. 4. Optimization Technique In NPR system, the PR conditions are certainly relaxed Design of M-Band NPR Cosine-Modulated Filterbank Using IFIR Technique Copyright © 2010 SciRes. JSIP 38 by allowing small amount of error. There are three types of errors occurs at the reconstructed output, viz., ampli- tude, phase and aliasing [4,5]. The stopband attenuation plays an important role in analysis of these parameters. High stopband attenuation results smaller aliasing error but larger reconstruction error. The peak to peak recon- struction and the peak aliasing errors are define in (19-20) [10] 00 max min jj pp EMTeMTe (19) max j a EEe (20) where, 12 12 1 M jj l l EeT e In NPR system, the measures of these error parameters are approximately and are given by the following equa- tions: 0for / j Pe M (21) 0 2 / 0 1, 21 where 0 j jkM j Te M Te Pe k (22) The accuracy of the first approximation gives aliasing and accuracy of the second gives the reconstruction error. The phase error can be eliminating completely by using linear phase prototype filter. Other two distortion pa- rameters can be minimized by applying suitable optimi- zation technique. Much of work have been done in this field [5,14-16]. Initially, Johnston [16] developed a nonlinear optimization technique, later on many promi- nent authors such as Creusere et al. [14] and Lin et al. [15] and others [5] have simplified it and developed sin- gle variable linear optimization techniques with different objective functions. In the proposed work, a linear gra- dient optimization technique is used with objective func- tion given in (23). The cutoff frequency of the model filter is selected as variable parameter to optimize the objective function. 2 2 max( 1 jM j Pe Pe (23) for 0≤ ω < π / M Initially, input parameters, i.e., sampling rate, number of bands, stretch factor, passband and stopband frequen- cies, passband ripple and stopband attenuation of proto- type filter are specified. Based on these input parameters the passband and stopband frequencies of model filter and interpolator filter are calculated. Determine the cut- off frequency, transition band and filter length of model filter as well as interpolator filter. Initialize, different optimization pointers like step size, search direction, flag and initial as well as expected minimum possible values of objective function. Inside the optimization loop, de- sign the model filter, interpolated filter, up-sampled model filter and finally, obtained the prototype lowpass filter using Equation (4). Determine the other bandpass filters of analysis and synthesis sections using cosine modulation. In optimization routine cutoff frequency of model filter is gradually change as per the search direc- tion and calculate the corresponding value of objective function. Algorithm comes out from the loop as it at- tained the minimum value of objective function and ob- tained the optimized value of reconstruction error. The flowchart of the optimization algorithm is given in appendix and implemented on MATLAB 7.0 on Pentium IV processor. 5. Design Examples In this section, the performance evaluation of the pro- posed technique has been done with the help of design examples. The performance of the filterbank is examined in terms of stopband attenuation, reconstruction error, aliasing error, computational complexity, group delay and stopband energy in used prototype filter. The effec- tiveness of proposed work is compared with the earlier reported work for same input parameters. Example 1: An eight-band cosine modulated filter- bank has been designed for same specifications as given by Kha et al. [10] with stopband attenuation (As) = 35.8dB, stopband frequency s = 0.12 . Figure 2(a), shows the frequency response of the prototype filter designed by standard FIR technique. In this proposed work, the prototype filter is designed using IFIR tech- nique. For the given design parameters of prototype filter, the required model, up-sampled model and interpolator filters are obtained for stretched factor 2L. The stop- band attenuation and passband ripple for these filters are same as desire in prototype filter, however the stopband frequency is different. For model and interpolator filters these values are s L = 0.24 and 2/ s L = 0.88 , respectively. The obtained magnitude responses of these filters are shown in Figure 2(b). The magnitude response of the resultant IFIR filter is shown in Figure 2(c). Filter coefficients are optimized by the proposed gradient optimization algorithm (given in Appendix) to reduce the reconstruction error. The optimized eight- band CMFB is shown in Figure 2(d). The obtained val- ues of reconstruction error, aliasing error, group delay, stopband energy and the required number of multipliers and adders in case of Kha et al. [10] and in proposed IFIR technique are given in Table-1. The magnitude re- sponses of reconstruction and aliasing errors are shown Design of M-Band NPR Cosine-Modulated Filterbank Using IFIR Technique Copyright © 2010 SciRes. JSIP 39 in Figures 3(a) and 3(b), respectively. Example 2: In this example, sixteen-band cosine modulated filterbank is designed using same specifica- tions as given in [10]. The stopband attenuation and stopband frequency are taken as same, i.e., (As) = 45dB, s = 0.0590 . The stopband attenuation for model and interpolator are same as in prototype filter, i.e., 45dB. The stopband frequency for these filters are s L = 0.118 and 2 s L = 0.941 , respectively. The magnitude response of prototype IFIR filter is shown in Figure 4(a) and the magnitude responses of sixteen band CMFB is shown in Figure 4(b). The optimized value of reconstruction error, aliasing error and reduction in computation cost, stopband energy and group delay are given in Table 1. The reconstruction error and aliasing error are shown in Figure 4(c) and 4(d), respectively. Example 3: A thirty two-band cosine modulated fil- terbank has been designed for same specifications as given by Kha et al.[10] with stopband attenuation (As) = 100 dB, stopband frequency s =0.03125 . The same stopband attenuation, i.e., 100dB is taken for model filter as well as for interpolator. The stopband frequency for these filters are s L = 0.0625 and 2 s L = 0.968 , respectively. The magnitude response of prototype IFIR filter is shown in Figure 5(a) and the magnitude responses of first and last four bands of thirty two-band CMFB is shown in Figure 5(b). The reconstruction error and aliasing error are shown in Fig- ure 5(c) and 5(d), respectively. The optimized value of reconstruction error, aliasing error and reduction in computation cost and group delay are given in Table 1. Example 4: In this example, the performance of pro- posed technique is compared with the reported work of Furado et al. [9]. Same input parameters as reported are taken for comparison purposes. The obtained error pa- rameters from the proposed technique are given in Table 2. The plot of reconstruction error and aliasing error are given in Figure 6. Figure 2. (a) FIR lowpass filter, (b) (–––) Up-sampled model,(– – – –)model filter,(--------) Interpolator filter, (c) IFIR proto- type lowpass filter, (d) Eight-band CMF bank. Figure 3. Magnitude response of (a) Reconstruction error, (b) Aliasing error for eight-band filterbank. Design of M-Band NPR Cosine-Modulated Filterbank Using IFIR Technique Copyright © 2010 SciRes. JSIP 40 Figure 4. Magnitude response of (a) IFIR prototype filter, (b) Sixteen-bands CMF bank(c) Reconstruction error, (d) Aliasing error. Table 1. Performance comparison. As Filter order Arithmetic elements %Computa- tional saving Work Approach dB M ωs N NmNiEpp Ea E s Ad- ders Mult- liers g FIR 35.8 8 0.12π 40 - -5.50e-32.47e-31.18e-240 20 20 - FIR 45.0 160.059π 102 - -5.95e-3 3.89e-4 9.00e-4102 51 51 - FIR 100 320.031π 466 - -9.12e-4 2.38e-7 9.9e-10466 233 233 35.8 8 40 - - 20 20 45.0 16 102 - - 51 51 Kha et al.[10] Polyphase 100 32 466 - - 233 233 35.8 8 0.12π - 20065.46e-31.41e-31.0e-226 13 13 35 45.8 16 0.059π - 46062.1e-32.62e-41.0e-352 26 26 49 Pro- posed IFIR 100 320.031π - 267153.3e-31.80e-78.8e-10282 141 141 39 Design of M-Band NPR Cosine-Modulated Filterbank Using IFIR Technique Copyright © 2010 SciRes. JSIP 41 Figure 5. Magnitude response of (a) IFIR prototype filter, (b) Thirty two-band CMF bank (c) Reconstruction error, (d) aliasing error. Table 2. Performance comparison with Furado et al. [9]. work Approach As M ωs N Nm N i E pp E a Furado et al. [9] FRM 60 8 0.035π639- - - 2.1e-4 Proposed IFIR 60 8 0.035π- 12208 9.1e-3 0.92e-5 00.1 0.2 0.3 0.4 0.50.5 -5 0 5x 10-3 Magnitude Normal ized Frequ ency 00.1 0.2 0.3 0.4 0.50.5 0 0.2 0.4 0.6 0.8 11x 10 -5 Magnitude Normalized Frequency (a) (b) Figure 6. Plot of magnitude response of (a) reconstruction error (b) aliasing error. 6. Discussion In last two decades, lot of cost effective techniques have been developed for the design of the filterbank. IFIR approach is one of the techniques, which provides sig- nificant reduction in computational cost. In this proposed work, the reduction in computation is achieved by using stretch factor of two and better results are obtained. Design of M-Band NPR Cosine-Modulated Filterbank Using IFIR Technique Copyright © 2010 SciRes. JSIP 42 In the first example, 8-bands CMFB is designed. The values of resulting reconstruction and aliasing error in Kha et al. [10] and proposed work are 5.508 × 10-3; 2.477 × 10-3 and 5.46 × 10-3; 1.4 × 10-3 respectively. Hence, the proposed method offers much better per- formance in terms of reconstruction and aliasing errors. Apart from lower values of error parameters the pro- posed method provides reduction in stopband energy, less group delay and 35% reduction in computational cost in contrast to Kha et al. Similarly, in second example, 16-band CMFB is de- signed. It is quite clear from Table 1 that the value of error parameters in the proposed work is much smaller than the Kha et al. [10] with 49% reduction in computa- tional cost. In third example, 32-band CMFB is designed with the same input parameters as of Kha et al. [10]. As it is clear from the results given in Table 1 that the pro- posed technique provides smaller value of aliasing error then reported work with reduction in computational cost in terms of multipliers. The fourth example is quoted for the comparison with the work of Furado et al. [9]. In this approach the proto- type filter is designed using FRM technique and the re- quired number of filter order is very large, therefore computational complexity is also high. However, the proposed IFIR technique provides smaller value of alias- ing error using much lower value of filter order. The results in Table 1 indicate that the proposed method provides better performance in terms of compu- tational complexity then the equivalent polyphase im- plementation. It is also clear from the magnitude re- sponse of prototype filter that there is improvement in SLFOR, which provides reduction in stopband energy of prototype filter. The proposed method is not only suitable for lower range of stopband attenuation, but also appli- cable at higher value of stopband attenuations. 7. Conclusion Efficient design for M-band NPR CMFBs has been pre- sented. IFIR technique is used to design the prototype filter, which provides significant reduction in computa- tion cost as well as improvement in SLFOR of the proto- type filter. The obtained values of reconstruction and aliasing errors are much smaller than the error values in previous publications for wide range of stopband at- tenuation. Therefore, the proposed method is useful in variety of applications such as audio and image compres- sion. The improvement in SLFOR can be utilized for better suppression of cross-talk and echo cancellation. The improvement in proposed approach is possible by using multistage IFIR prototype filter. REFERENCES [1] W. P. Zhu, M. O. Ahamed and M. N. S. Swamy, “An Efficient Approach for the Design of Nearly Perfect Re- construction QMF Banks,” IEEE Transactions on Cir- cuits and System II: Analog and Digital Signal Process- ing, Vol. 45, No. 8, August 1998, pp. 1161–1165. [2] Y. J. Chen and K. S. Amaratunga, “M-Channel Lifting Factorization of Perfect Reconstruction Filter Banks and Reversible M-Band Wavelet Transforms,” IEEE Trans- actions on Circuits and System II: Analog and Digital Signal Processing, Vol. 50, No. 12, December 2003, pp. 963-976. [3] C. K. Goh and Y. C. Lim, “Novel Approach for the De- sign of Two Channel Perfect Reconstruction Linear Phase Fir Filter Banks,” IEEE Transactions on Circuits and System II: Analog and Digital Signal Processing, Vol. 45, No. 8, August 1998, pp. 1141-1146. [4] P. P. Vaidyanathan, “Multirate Systems and Filter Banks,” Prentice-Hall, Englewood Cliffs, New Jersey, 1993. [5] A. Jain, R.Saxena and S. C. Saxena, “An Improved and Simplified Design of Cosine Modulated Pseudo-QMF Filterbanks,” Digital Signal Processing, Vol. 16, No. 3, May 2006, pp. 225-232. [6] T. Q. Nguyen, “Near-Perfect-Reconstruction Pseudo-QMF Banks,” IEEE Transactions on Signal Processing, Vol. 42, No. 1, January 1994, pp. 65-76. [7] T. Q. Nguyen, “Digital Filterbank Design Quadratic-Con- strained Formulation,” IEEE Transactions on Signal Processing, Vol. 43, No. 9, September 1995, pp. 2103- 2108. [8] Y. Neuvo, C. Y. Dong and S. K. Mitra, “Interpolated Finite Impulse Response Filters,” IEEE Transactions on Acoust, Speech, Signal Processing, Vol. ASSP-32, June 1984, pp. 563-570. [9] M. B. Furado, P. S. R. Diniz, S. L. Netto and T. Saramaki, “On the Design of High-Complexity Cosine-Modulated Transmultiplexers Based on the Frequency-Response Masking Approach,” IEEE Transactions on Circuits and System, Vol. 52, No. 11, November 2005, pp. 2413-2426. [10] H. H. Kha, H. D. Tuan and T. Q. Nguyen, “Efficient De- sign of Cosine-Modulated Filterbanks via Convex Opti- mization,” IEEE Transactions on Signal Processing, Vol. 57, No. 3, March 2009, pp. 966-976. [11] A. Mehrnia and J. A. Willson, “On Optimal IFIR Filter Design,” IEEE Proceedings of the 2004 International Symposium on Circuits and Systems, Vol. 3, May 2004, pp. 133-136. [12] S. K. Mitra, “Digital Signal Processing-A computer based approach,” Tata McGraw-Hill Publishing Co. Ltd., New Delhi, 1998. [13] O. P. Sahu, M. K. Soni and I. M. Talwar, “Marquardt Optimization Method to Design Two-Channel Quadrature Mirror Filterbanks,” Digital Signal Processing, Vol. 16, No. 6, November 2006, pp. 870-879. [14] C. D. Creusere and S. K. Mitra, “A Simple Method for Designing High-Quality Prototype Filters for M-Band Design of M-Band NPR Cosine-Modulated Filterbank Using IFIR Technique Copyright © 2010 SciRes. JSIP 43 Pseudo QMF Banks,” IEEE Transactions on Signal Processing, Vol. 43, No. 4, April 1995, pp. 1005-1007. [15] Y. P. Lin and P. P. Vaidyanathan, “A Kaiser Window Approach for the Design of Prototype Filters of Co- sine-Modulated Filterbanks,” IEEE Signal Processing Letters, Vol. 5, No. 6, June 1998, pp. 132-134. [16] J. D. Johnston, “A Filter Family Designed for Use in Quadrature Mirror Filterbanks,” IEEE International Con- ference on Acoustics, Speech and Signal Processing, Denver, 1980, pp. 291-294. Appendix |