132 IEEE SIGNAL PROCESSING LETTERS, VOL. 5, NO. 6, JUNE 1998
A Kaiser Window Approach for the Design of
Prototype Filters of Cosine Modulated Filterbanks
Yuan-Pei Lin,
Member, IEEE,and P. P. Vaidyanathan,
Fellow, IEEEAbstract— The traditional designs for the prototype filters of cosine modulated filterbanks usually involve nonlinear optimiza-tions. We propose limiting the search of the prototype filters to the class of filters obtained using Kaiser windows. The design process is reduced to the optimization of a single parameter. An example will be given to show that very good designs can be obtained in spite of the limit of search.
Index Terms—Cosine-modulated filterbanks, filterbank design, Kaiser window.
I. INTRODUCTION
T
HE THEORY and design of -channel cosine modulated filterbanks have been studied extensively in the past (see, e.g., [1] for a list of references). The analysis and synthesis filters are cosine modulated versions of a lowpass prototype filter (as shown in Fig. 1). The design of the whole filterbank thus reduces to that of the prototype filter. The filterbank has perfect reconstruction if the polyphase components of the pro-totype satisfy some pairwise power complementary conditions [1]. The optimization of the prototype coefficients in this case is highly nonlinear. In an approximate reconstruction cosine modulated filterbank, aliasing is canceled approximately and the distortion is a delay only approximately. Suppose the prototype filter has linear phase. The conditions for approximate reconstruction can be stated in terms of as follows [1], [2]:for (1)
where
(2) The accuracy of the first approximation will in turn give a measure of the aliasing error, while the accuracy of the second approximation gives us a measure of the distortion error.
Traditional design approaches of the prototype usually in-volve nonlinear optimizations [3]. In [4], Creusere and Mitra proposed a very efficient prototype design method without using nonlinear optimizations. Instead of a full search, the search in [4] is limited to the class of filters obtained using the
Manuscript received December 24, 1997. This work was supported in part by the National Science Foundation under Grant MIP 0703755, and by Tektronix, Inc. The associate editor coordinating the review of this manuscript and approving it for publication was Prof. A. Tewfik.
Y.-P. Lin is with the Department of Electrical and Control Engineering, National Chiao Tung University, Hsinchu, Taiwan, R.O.C.
P. P. Vaidyanathan is with the Department of Electrical Engineering, California Institute of Technology, Pasadena, CA 91125 USA.
Publisher Item Identifier S 1070-9908(98)04338-7.
Fig. 1. Typical magnitude response of the prototype for anM-channel cosine modulated filterbank.
Parks–McClellan algorithm. As a result, the optimization can be reduced to that of a single parameter. The design examples provided therein have shown that very good prototypes that give close approximations of (1) and (2) can be obtained despite of the limit of search. Inspired by this, we show how the design method can be improved further by using the Kaiser window method.
In this letter, we consider limiting the search to the class of filters obtained using Kaiser window designs. Using this method, the optimization is also reduced to that of a single parameter. But the optimization procedure is even simpler than in [4] due to the fact that the filters designed by using a Kaiser window can be expressed in closed form. We will also propose the use of a new objective function that significantly reduces the amount of computation for the objective function. Experiments show that the objective function is a convex function of the parameter being optimized and that very good designs can be obtained. Furthermore, the resulting filterbank can be described by only four parameters regardless of the number of channels in the filterbank and the order of the prototype. This technique has the additional advantage that we have direct control over the stopband attenuation of the prototype filter and, hence, the stopband attenuation of the analysis and synthesis filters.
II. REVIEW OF PROTOTYPE DESIGNS USING THE PARKS–MCCLELLEN ALGORITHM [4]
In the design of a lowpass filter using the Parks–McClellan algorithm, there are four parameters to be determined, the passband edge the stopband edge passband-stopband error ratio and filter length The first condition for approximate reconstruction (1) can be achieved by letting the stopband edge Now consider the second condition. Observing the expression of in (2) we see that is periodic with period and we only need to consider
in the interval With we can
LIN AND VAIDYANATHAN: COSINE MODULATED FILTERBANKS 133
approximate by
for Let
(3) For a given filter length and an appropriately chosen passband-stopband error ratio the passband edge is opti-mized to minimize Such a search reduces the design to only the optimization of the passband edge Experimentally, the objective function is shown to be a convex function of
in the interval [4].
III. KAISER WINDOW APPROACH
A filter of length designed through Kaiser
windows is of the form [5]
where
is the impulse response of the ideal filter with cutoff frequency and is a Kaiser window. The window is defined by
otherwise where is the zeroth-order modified Bessel function, which can be computed using
The parameter can be determined by the desired stopband attenuation [5]. So a filter designed using Kaiser windows can be specified by the three parameters , and
Let us consider designing the prototype filter using Kaiser windows. Given a desired stopband attenuation and an appropriately chosen transition bandwidth the order of the window can be estimated by the following formula developed by Kaiser:
(4) In this case, the window is completely determined. The cutoff frequency is the only parameter left in the Kaiser window design. Now we can optimize such that the objective function in (3) is minimized as in [4]. In what follows, we propose the use of a different objective function that is much easier to compute than that in (3).
Define a filter given by
Observe that the condition in (2) means that is a
(a)
(b)
(c)
(d)
Fig. 2. (a) Objective function new versus the cutoff frequency !c: (b) Magnitude ofP0(ej!) (dB). (c) Zoom plot of (b) for the interval [0 2=M]:
(d) Magnitude ofT (ej!): As T (ej!) is periodic =M; only the part in the interval[0 =M] is shown.
approximately a Nyquist filter. That is to say,
is approximately where is given by
, otherwise.
Using this, we can choose a simpler objective function (5) We can adjust the parameter to find the best that yields
the smallest Experiments show that is a convex
function of Very good design can be obtained, as will be demonstrated in the next section.
IV. DESIGN EXAMPLE
We use the proposed Kaiser window approach to design a cosine modulated filterbank with specifications similar to that of a Motion Pictures Expert Group (MPEG) audio coder. In
134 IEEE SIGNAL PROCESSING LETTERS, VOL. 5, NO. 6, JUNE 1998
this case, the filterbank has 32 channels and the prototype filter has stopband attenuation around 100 dB (pp. 397, [6]). So in the design we choose to be 100 dB. The transition bandwidth is chosen to be slightly smaller than The value of the transition bandwidth used in this example is
The order of the prototype estimated
by (4), is 466. Fig. 2(a) gives the plot of the objective function versus the cutoff frequency From the plot
we see that is a convex function of when
the objective function has a minimum of 5.63 10 . Fig. 2(b) shows the magnitude response (in dB) of the optimized prototype filter. A zoom plot of Fig. 2(b) for the interval is given in Fig. 2(c). The magnitude response of the distortion function in this case is given in Fig. 2(d).
The design takes only a few seconds using a PC with a Pentium processor.1 A similar example using nonlinear
optimizations will take approximately 5 min on a Sparc 10 workstation [3].
V. CONCLUSIONS
In this letter, we propose designing the prototype filters of cosine modulated filterbanks using Kaiser windows. This technique is simple conceptually and computationally. The design of the prototype filter is formulated as a problem of 1A Matlab program code of prototype designs using Kaiser windows is
available at http://multirate.cn.nctu.edu.tw.
optimizing the cutoff frequency in Kaiser window designs. We further reduce the design complexity by introducing a new objective function Experiments show that is a con-vex function of We have demonstrated very good designs can be obtained, even when the desired stop attenuation is large and the filter order is high. Design examples have also shown that the design complexity of the Kaiser window approach is much less than that in nonlinear optimizations. On the other hand, because in Kaiser designs the stopband attenuation is determined by the window, we can have direct control over the stopband attenuation. Furthermore the coefficients of the prototype can be expressed in closed form in terms of only three parameters, filter order and cutoff frequency So the whole cosine modulated filterbank can be specified by these three parameters and the number of channels
REFERENCES
[1] P. P. Vaidyanathan, Multirate Systems and Filter Banks. Englewood Cliffs, NJ: Prentice-Hall, 1993.
[2] J. H. Rothweiler, “Polyphase Quadrature Filters—A new subband cod-ing technique,” in Proc. IEEE Int. Conf. Acoustics, Speech, and Signal
Processing, Boston, MA, Apr. 1983, pp. 1280–1283.
[3] G. Strang and T. Nguyen, Wavelets and Filter Banks. Wellesley, MA: Wellesley-Cambridge, 1996.
[4] C. D. Creusere and S. K. Mitra, “A simple method for designing high-quality prototype filters for M-band pseudo QMF banks,” IEEE Trans.
Signal Processing, vol. 43, pp. 1005–1007, Apr. 1994.
[5] A. V. Oppenheim and R. W. Schafer, Discrete-Time Signal Processing. Englewood Cliffs, NJ: Prentice-Hall, 1989.
[6] M. A. Vetterli and J. Kovacevic, Wavelets and Subband Coding. En-glewood Cliffs, NJ: Prentice-Hall, 1995.