The technique of filter banks has been extensively applied in signal processing in the last three decades. It provides a very efficient way of signal decomposition, characterization, and analysis. It is also the main driving idea in almost all frequency division multiplexing technologies. With the advent of wavelets and subsequent realization of its wide area of application, filter banks became even more important as it has been proven to be the most efficient way a wavelet system can be implemented. In this chapter, we present an analysis of the design of a wavelet transform using the filter bank technique. The analysis covers the different sections which make up a filter bank, i.e., analysis filters and synthesis filters, and also the upsamplers and downsamplers. We also investigate the mathematical properties of wavelets, which make them particularly suitable in the design of wavelets. The chapter then focuses attention to the particular role the analysis and the synthesis filters play in the design of a wavelet transform using filter banks. The precise procedure by which the design of a wavelet using filter banks can be achieved is presented in the last section of this chapter, and it includes the mathematical techniques involved in the design of wavelets.
- filter bank
- perfect reconstruction
- paraunitary condition
Filter banks can be defined as the cascaded arrangement of filters, i.e., low-pass, high-pass, and band-pass filters connected by sampling operators in such a manner as to achieve the decomposition and recomposition of a signal from a spectrum perspective. The sampling operators could either be downsamplers or upsamplers. The downsamplers are called decimators while the upsamplers are called expanders. The technique of filter banks plays an important role in most digital systems that rely on signal processing for their operations. Using this technique, any signal feature can be reliably extracted and analyzed; hence filter banks have wide applications in digital signal processing systems. A filter bank as shown in Figure 1 [1, 2] consists of different parts, which collectively execute a desired function.
As can be seen in Figure 1, the filter bank is made of two sections: the analysis filter bank section (composed of analysis filters and downsamplers), and the synthesis filter bank section (composed of upsamplers and synthesis filters). In this chapter, we will discuss the analysis and synthesis filter bank sections, their responses to incoming signals, and how they work together in the derivation of a wavelet transform function.
2. Analysis filter bank section
The analysis filter bank section is made up of the analysis filter banks, and downsamplers or decimators which together act on an input signal to perform a desired function through decomposition of the signal. In this section, we will analyze the mathematical relationship that exists between these two components. To have a thorough understanding of this relationship, it is important to briefly discuss these components separately.
2.1 Analysis filter bank
The filters that make up the analysis filter banks could either be low-pass filters, or high-pass filters. Each of these filters, as shown in Figure 2, allows the passage of only a particular frequency component of the input signal . Thus, specific features of the input signal embedded at different frequencies can be individually extracted and investigated using the analysis filter bank [3, 4]. The k-channel filter bank in Figure 2 separates the frequencies of the input signal in the manner presented.
It can be seen from the frequency responses that the output of the filters overlap each other. This is because in practice, the filters are not ideal. However, the overlapping condition can be improved through an optimized design of the filters. Mathematically, the effect of each of the filters in the filter bank on the input signal can be stated as follows:
where is the z-transform of the result from the convolution operation between the z-transform of the input signal and the z-transform of the filter . The output in Figure 2 is fed into the corresponding downsampler of Figure 1. In the next section, we will analyze the downsampler and state the mathematical operation it performs on a given signal.
The downsampler shown in Figure 1 downsamples an input signal by a factor of
Figure 3 shows the conceptual depiction of the relationship in Eq. (2).
Mathematically, the output of the decimator in Figure 1 can be expressed as a product of the input sequence
and the sequence of unit impulses which are
The relationship in Eq. (3) will only select the
Setting and , the relationship in Eq. (4) becomes:
Substituting Eq. (5) into (3) yields:
In terms of z-transformation, the relationship in Eq. (6) can be expressed as:
Having looked at the analysis filters and downsamplers, we will now turn our attention to synthesis filter bank section of Figure 1.
3. Synthesis filter bank section
The synthesis filter bank section is made of the upsamplers and synthesis filter banks. These components work together to perform the opposite operation performed by the analysis filter bank section shown in Figure 1. In this section, we will make an analysis of the mathematical relationship that governs the operation of the synthesis filters and upsamplers.
3.1 Synthesis filter bank
Similar to the analysis filter bank, the synthesis filter bank is made of low-pass and high-pass filters. The output of these filters as shown in Figure 1, are summed to a common output. In typical filter bank applications, the frequency responses of these filters are typically matched to those of the analysis filters shown in Figure 2. The mathematical expression for the effect each of these filters has on the corresponding input signal is as stated below :
In Figure 2, the input to the synthesis filter bank is upsamplers or expanders. The next section gives a brief review of the upsamplers.
The upsampler expands an input signal by a factor
Given an input signal in Figure 1, an upsampler with a factor of 2 will upsample the signal using the relationship :
Similar to the expression in Eq. (3), the z-transform of the expression in Eq. (9) which is an upsampler is stated as follows :
To be useful in wavelet designs, filter banks must be designed to have certain characteristics which guarantee that a signal at the input of a filter bank will be received accurately at the output of the filter bank. In the next section, we will examine the properties of filter banks and how these properties influence the design of wavelets.
4. Properties of filter banks for wavelet design
In wavelet designs, filter banks are required to possess three important properties which are fundamental to the realization of a wavelet function. These properties include: perfect reconstruction, orthogonality, and paraunitary condition.
4.1 Perfect reconstruction
This property guarantees that the signal at the output of a given filter bank is a delayed version of the signal at the input of the filter bank. Perfect reconstruction is an important property of a filter bank because it cancels the effect of aliasing of the input signal at the output, caused by the downsamplers and upsamplers. To understand this point, consider a two-channel finite impulse response FIR filter bank shown in Figure 5.
The output is derived using Eqs. (6) and (10) as follows in terms of the signal component and aliasing component as:
where the signal_component and aliasing_component are defined as:
To achieve perfect reconstruction, the following condition must be satisfied :
The relationships in Eqs. (11) and (13) are possible when the filter bank is constructed as a QMF (quadrature mirror filter) filter bank or CQF (conjugate quadrature filter) filter bank. Both QMF and CQF banks provide a mechanism by which complete cancellation of the aliasing component in Eq. (11) can be accomplished. Using QMF, aliasing cancellation can be achieved by constructing the filters in Figure 5 based on the following relationships [4, 5]:
In Eq. (14), the synthesis filter has the same coefficients as the analysis filter ; the analysis filter has the same coefficients as the analysis filter , but every other value is negated; the synthesis filter is a negative copy of the analysis filter . For example, if the analysis filter has coefficients , then the filter bank in Figure 5 will assume the structure shown in Figure 6.
For the CQF bank, the coefficients of the analysis filter are a reversed version of the analysis filter with every other value negated. The synthesis filters and are a reversed versions of the analysis filters and , respectively. These relationships can be stated mathematically as follows :
Based on the relationship in Eq. (15), the filter bank shown in Figure 6 for CQF will assume the structure shown in Figure 7.
Based on the structure of Figures 6 or 7, the output signal is related to the input signal by the expression:
If we impose the condition that , then Eq. (16) becomes:
The relationship in Eq. (17) states that the output signal is delayed version of the input signal by three samples. We leave the verification of Eq. (16) as an exercise for the reader.
Having looked at perfect reconstruction as a necessary property for a filter bank in wavelet design, we now look at orthogonality as also an essential property for a filter bank in the design of wavelets.
Orthogonality in a filter bank is a situation in which the synthesis filter bank is a transpose of the analysis filter bank. This is a useful property in the sense that it allows for the energy preservation of the signal being processed. This important property is achieved through the imposition of the orthogonality condition on both the analysis and filter bank sections while at the same time preserving the perfect reconstruction condition of the filter bank. The imposition of the orthogonality condition in a filter bank (see Figure 5) occurs when the following relationships are satisfied :
In Eq. (18), the inner product of the coefficients of the synthesis filter and the analysis filter must be zero and the inner product of the coefficients of the synthesis filter and the analysis filter must also be zero for the orthogonality condition to hold.
Also, the low-pass analysis filter is related to the other three filters through the following expressions :
The condition in Eq. (20) also describe the necessary requirement for a filter bank to be paraunitary (which we shall examine in the next section), i.e., the low-pass filter satisfy the following power symmetry of halfband condition [8, 9]:
where . If the low-pass filter satisfies the required symmetry condition:
then is said to be a real filter. The implication of the constraint in Eq. (21) is that and be antisymmetric filters, and is a symmetric filter. The relationships in Eqs. (20)–(22) give the necessary and sufficient condition for the characterization of a filter bank with orthogonality and symmetry.
The orthogonality condition for a filter bank can also be examined from a polyphase perspective. Consider the polyphase representation of the filter bank in Figure 5 as illustrated in Figure 8 .
If in Figure 8 is type-I analysis polyphase matrix, and is type-II synthesis polyphase matrix, then :
The conditions in Eqs. (20)–(22) hold true iff and satisfy the following:
where , with the first and second condition in Eq. (24) relating to the filter bank orthogonality condition, and the last represents the filter bank symmetry.
We now look at the paraunitary condition of a filter bank, which is also a necessary property in filter bank implementation of wavelets.
4.3 Paraunitary condition
In the filter bank implementation of a wavelet transform, the paraunitary condition plays the critical role of guaranteeing the generation of orthonormal wavelets, and also perfect recovery of a decomposed signal. The paraunitary condition guarantees that recovered signal will suffer no phase or aliasing effect if a filter bank satisfies the paraunitary condition .
Given a polyphase transfer function matrix , the paraunitary condition is established by the matrix iff :
where the H superscript denotes the conjugated transpose, and
If the real-coefficient lossless matrix is denoted by ; then the matrix is said to have a special case of lossless degree of one iff it can be characterized by the relationship :
where is an arbitrary unitary matrix and is an column vector with unit norm. From Eq. (26), the paraunitary condition for a filter bank is obtained as follows :
Having looked at the filter bank and its three important properties for the design of a wavelet, we will in the next section examine the application of these properties in the design of a wavelet.
5. Filter bank design of a wavelet transform
The filter bank design of a wavelet transform is usually implemented from the analysis filter bank segment to the synthesis filter bank segment.
5.1 Analysis filter bank in wavelet transform design
Given that the expression for a scaling function is the series sum of the shifted versions of , then according to [15, 16], can be represented as:
where denotes the scaling coefficients. If is transformed such that , then the relationship in Eq. (28) becomes :
which translates into:
In a similar consideration to Eq. (28), the wavelet function can be represented as :
where denotes the wavelet coefficients. Also, if is transformed such that , then the relationship in Eq. (31) becomes :
which translates into:
5.2 Synthesis filter bank in wavelet transform design
In the synthesis filter bank, the reconstruction of the original coefficients of a signal can be achieved through the combination of the scaling and wavelet function coefficients at a coarse level of resolution. Given a signal at scaling space , then according to [16, 17], the reconstruction is derived as follows:
For the next scale, Eq. (34) becomes:
Substituting Eqs. (28) and (31) into Eq. (35) and after algebraic manipulations yields :
6. Wavelet transform design procedure using filter banks
In the design of a wavelet system using filter banks, it is of utmost importance that the filters which will execute the filter bank system as shown in Figure 1, possess the properties discussed in Section 4. Owing to the fact that in a filter bank, all the filters can be derived from an initial filter as described in Eq. (13), then this initial filter must be designed in such a manner that the relationships in Sections 5.1 and 5.2 are realized. To this end, the following steps as shown in the state diagram in Figure 10 are necessary.
In the first state in Figure 10, the design problem formulation which can be achieved using trigonometric polynomial, takes the following into consideration :
Compact support which guarantees that the wavelet is characterized by finite non-zero coefficient.
Paraunitary condition which guarantees the generation of orthonormal wavelets.
Flatness/k-regularity which guarantees the smoothness of the wavelet in both time and frequency domains.
The second state which involves conditioning the problem as a tractable problem involves, if necessary, transforming a non-linear formulation of the problem to a linear formulation, and then optimizing the problem using techniques like convex optimization. The generation of the filter coefficients using solvers in the third state of the machine involves techniques like spectral factorization. Through simulation in the fourth state of the chart, the generated coefficients can be verified whether or not they meet the design constraints. Using the QMF or CQF relationships in Eqs. (13) and (14), the other filters in the filter bank are generated in the fifth state of the chart.
In this chapter, we have presented an analysis of the design of wavelets using filter bank technique. The chapter looked at the two major components of a filter bank which the analysis and the synthesis components. The properties of filter banks which are desirable in the design of wavelets were also investigated, alongside the mathematical description of these properties. The chapter also gave a brief mathematical description of the role the analysis and the synthesis filter banks play in the design of wavelets. Finally, the required general procedure for the design of wavelets was presented, showing the necessary steps to take in order to achieve an effective design.
The major contribution of this chapter is the provision of a step by step analysis and procedure for the design of filter banks in a precise and concise manner.
Mertins A. Signal Analysis: Wavelets, Filter Banks, Time-Frequency Transforms and Applications. Sussex: Wiley; 1999
Mon AT, Aye SM, Tun HM, Naing ZM, Moe WK. Analysis of multichannel filter banks-based tree-structured design for communication system. Software Engineering. 2018; 6(2):37-46
Mallat S. A Wavelet Tour of Signal Processing: The Sparse Way. 3rd ed. Burlington: Academic Press; 2009
Tay DB, Ortega A, Anis A. Cascade and lifting structures in the spectral domain for bipartite graph filter banks. In: APSIPA Annual Summit and Conference. Cambridge University Press; 2018. pp. 1141-1147
Chu E. Discrete and Continuous Fourier Transforms: Analysis, Applications and Fast Algorithms. Boca Raton: CRC Press; 2008
Abbaszadeh MH, Khalaj BH, Haghbin A. Optimum low complexity filter bank for generalized orthogonal frequency division multiplexing. EURASIP Journal on Wireless Communications and Networking. 2018; 2018(12):1-14
Vantita S, Rajesh M. Rational sampling rate converter design analysis using symmetric technique. International Journal of Computer Trends and Technology. 2015; 27(2):116-120
Devis T, Manuel M. Multirate and filterbank approaches in digital hearing aid design: A review. In: IOP Materials Science and Engineering Conference. IOP Publishing; 2018. pp. 1-8
Agrawal SK, Sahu OP. Two channel quadrature mirror filter bank: An overview. Hindawi—ISRN Signal Processing. 2013; 2013:1-10
Giron-Sierra JM. Digital Signal Processing with MATLAB Examples. Singapore: Springer Singapore; 2017
Vetterli M, Herley C. Wavelets and filter banks: Theory and design. IEEE Transactions on Signal Processing. 1992; 40(9):2207-2231
Vaidyanathan PP. Multirate Systems and Filter Banks. New Jersey: Prentice-Hall; 1993
Gao XQ, Nguyen TQ, Strang G. A study of two-channel complex-valued filterbanks and wavelets with orthogonality and symmetry properties. IEEE Transactions on Signal Processing. 2002; 50(4):824-833
Nikookar H. Wavelet Radio: Adaptive and Reconfigurable Wireless Systems Based on Wavelets. Cambridge: Cambridge University Press; 2013
Gao X, Nguyen TQ, Strang G. On factorization of M-channel paraunitary filterbanks. IEEE Transactions on Signal Processing. 2001; 49(7):1433-1466
Redif S, McWhirter JG, Weiss S. Design of FIR paraunitary filter banks for subband coding using a polynomial eigenvalue decomposition. IEEE Transactions on Signal Processing. 2011; 59(11):5253-5264
Vaidyanathan PP, Dongata Z, Nguyen TQ. More results on the perfect reconstruction problem in M-band parallel QMF bank. In: IEEE International Symposium on Circuits and Systems. IEEE; 1987. pp. 847-850
Vaidyanathan PP, Nguyen TQ, Doganata Z, Saramaki T. Improved technique for design of perfect reconstruction FIR QMF banks with lossless polyphase matrices. IEEE Transactions on Acoustics, Speech, and Signal Processing. 1989; 37(7):1042-1056
Burrus CS, Gopinath RA, Guo H. Introduction to Wavelets and Wavelet Transform: A Primer. Upper Saddle River: Prentice Hall; 1998