Open access peer-reviewed chapter

Analysis of Wavelet Transform Design via Filter Bank Technique

Written By

Peter Yusuf Dibal, Elizabeth Onwuka, James Agajo and Caroline Alenoghena

Submitted: 03 December 2018 Reviewed: 07 February 2019 Published: 27 September 2019

DOI: 10.5772/intechopen.85051

From the Edited Volume

Wavelet Transform and Complexity

Edited by Dumitru Baleanu

Chapter metrics overview

1,022 Chapter Downloads

View Full Metrics

Abstract

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.

Keywords

  • wavelet
  • filter bank
  • perfect reconstruction
  • orthogonality
  • paraunitary condition

1. Introduction

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.

Figure 1.

k-Channel filter bank [1, 2].

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.

Advertisement

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 y n . 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.

Figure 2.

Separation of input signals into sub-band frequencies by analysis filter bank.

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 y n can be stated as follows:

U 0 Z = Y Z H 0 Z U 1 Z = Y Z H 1 Z U 2 Z = Y Z H 2 Z U k 1 Z = Y Z H M 1 Z E1

where U i z is the z-transform of the result from the convolution operation between the z-transform of the input signal Y Z and the z-transform of the filter H i Z . The output U i z 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.

2.2 Downsampler/decimator

The downsampler shown in Figure 1 downsamples an input signal by a factor of N. This implies that it only retains all the Nth samples in a given sequence. For example, if N = 2 , then the downsampler will retain all even samples in a given sequence. Given an input signal x n , the downsampler with a factor of 2 will downsample the signal as:

x ̂ n = x 2 n , n E2

Figure 3 shows the conceptual depiction of the relationship in Eq. (2).

Figure 3.

Decimation by a factor of 2.

Mathematically, the output of the decimator in Figure 1 can be expressed as a product of the input sequence u i n and the sequence of unit impulses which are N samples apart, i.e.,

v i n = k u i n δ n kN , k E3

The relationship in Eq. (3) will only select the kNth sample of u i n , and the Fourier series expansion of the impulse series can be expressed as [5]:

k δ n kN = 1 N k = 0 N 1 e j 2 πkn / N E4

Setting W N = e j 2 π / N and n = 1 , the relationship in Eq. (4) becomes:

k δ n kN = 1 N k = 0 N 1 W N k E5

Substituting Eq. (5) into (3) yields:

v i n = 1 N k = 0 N 1 u i n W N k E6

In terms of z-transformation, the relationship in Eq. (6) can be expressed as:

V i Z = 1 N k = 0 N 1 U i Z 1 N W N k E7

Having looked at the analysis filters and downsamplers, we will now turn our attention to synthesis filter bank section of Figure 1.

Advertisement

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 w i n is as stated below [6]:

P 0 Z = W 0 Z G 0 Z P 1 Z = W 1 Z G 1 Z P 2 Z = W 2 Z G 2 Z P k 1 Z = W k 1 Z G M 1 Z E8

In Figure 2, the input to the synthesis filter bank is upsamplers or expanders. The next section gives a brief review of the upsamplers.

3.2 Upsampler/expander

The upsampler expands an input signal by a factor N. It does this by inserting zeros at every nth position in the sequence of the input signal. For example, if N = 2 , then the upsampler will insert a zero between every two adjacent samples in a given sequence as shown in Figure 4.

Figure 4.

Upsampling by a factor of 2.

Given an input signal v i n in Figure 1, an upsampler with a factor of 2 will upsample the signal using the relationship [7]:

w i n = k v i n δ n kN , k E9

Similar to the expression in Eq. (3), the z-transform of the expression in Eq. (9) which is an upsampler is stated as follows [8]:

W i Z = 1 N k = 0 N 1 V i Z N W N k E10

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.

Advertisement

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.

Figure 5.

Two-channel FIR filter bank.

The output y ̂ n is derived using Eqs. (6) and (10) as follows in terms of the signal component and aliasing component as:

Y ̂ z = signal _ component + aliasing _ component E11

where the signal_component and aliasing_component are defined as:

signal _ component = 1 2 F 0 z H 0 z + F 1 z H 1 z X z aliasing _ component = 1 2 F 0 z H 0 z + F 1 z H 1 z X z E12

To achieve perfect reconstruction, the following condition must be satisfied [1]:

F 0 z H 0 z + F 1 z H 1 z = 2 z 1 F 0 z H 0 z + F 1 z H 1 z = 0 E13

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]:

F 0 z = H 0 z H 1 z = H 0 z F 1 z = H 1 z E14

In Eq. (14), the synthesis filter F 0 z has the same coefficients as the analysis filter H 0 z ; the analysis filter H 1 z has the same coefficients as the analysis filter H 0 z , but every other value is negated; the synthesis filter F 1 z is a negative copy of the analysis filter H 1 z . For example, if the analysis filter H 0 z has coefficients p , q , r , s , then the filter bank in Figure 5 will assume the structure shown in Figure 6.

Figure 6.

QMF two-channel FIR filter bank.

For the CQF bank, the coefficients of the analysis filter H 1 z are a reversed version of the analysis filter H 0 z with every other value negated. The synthesis filters F 0 z and F 1 z are a reversed versions of the analysis filters H 0 z and H 1 z , respectively. These relationships can be stated mathematically as follows [10]:

H 1 z = z 1 H 0 z 1 F 0 z = H 1 z F 1 z = H 0 z E15

Based on the relationship in Eq. (15), the filter bank shown in Figure 6 for CQF will assume the structure shown in Figure 7.

Figure 7.

CQF two-channel FIR filter bank.

Based on the structure of Figures 6 or 7, the output signal y ̂ n is related to the input signal y n by the expression:

y ̂ n = pp + qq + rr + ss y n 3 E16

If we impose the condition that pp + qq + rr + ss = 1 , then Eq. (16) becomes:

y ̂ n = y n 3 E17

The relationship in Eq. (17) states that the output signal y ̂ n is delayed version of the input signal y n 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.

4.2 Orthogonality

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 [11]:

f ˜ 0 n 2 k h 1 n 2 l = 0 f ˜ 1 n 2 k h 0 n 2 l = 0 E18

where

g ˜ i n = g i n

and

g ˜ 0 n 2 k h 0 n = δ k g ˜ 1 n 2 k h 1 n = δ k E19

In Eq. (18), the inner product of the coefficients of the synthesis filter F 0 z and the analysis filter H 1 z must be zero and the inner product of the coefficients of the synthesis filter F 1 z and the analysis filter H 0 z must also be zero for the orthogonality condition to hold.

Also, the low-pass analysis filter H 0 z is related to the other three filters through the following expressions [12]:

H 1 z = cz L 1 H ˜ 0 z F 0 z = z L 1 H ˜ 0 z F 1 z = z L 1 H ˜ 1 z E20

where L denotes the length of the filter which must be even, and c is a constant with c = 1 ; H ˜ 0 z is the flipped and conjugated version of H 0 z , H ˜ 0 z is the conjugated version of H 0 z , and H ˜ 1 z is the conjugated version of H 1 z .

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 H 0 z satisfy the following power symmetry of halfband condition [8, 9]:

P z + P z = 2 E21

where P z = H 0 z H ̂ 0 z . If the low-pass filter H 0 z satisfies the required symmetry condition:

H 0 z = z L 1 H 0 z 1 E22

then P z is said to be a real filter. The implication of the constraint in Eq. (21) is that H 1 z and F 1 z be antisymmetric filters, and F 0 z 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 [13].

Figure 8.

Polyphase implementation of filter bank.

If E z in Figure 8 is type-I analysis polyphase matrix, and R z is type-II synthesis polyphase matrix, then [13]:

H 0 z H 1 z T = E z 2 1 z 1 T F 0 z F 1 z = z 1 1 R z 2 E23

The conditions in Eqs. (20)(22) hold true iff E z and R z satisfy the following:

E ̂ z E z = I R z = z k 1 E z E z = z k 1 diag 1 , 1 E z 1 J E24

where k = L / 2 , 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 [14].

Given a polyphase transfer function matrix E z , the paraunitary condition is established by the matrix iff [15]:

E H z 1 E z = I E25

where the H superscript denotes the conjugated transpose, and I denotes the identity matrix. Paraunitary filter banks also have an attractive property of losslessness, which implies that for every frequency, the total signal power is conserved [16]. From this property [17], any M × M real-coefficient lossless matrix with N 1 degree can be realized using the structure shown in Figure 9 [18].

Figure 9.

Cascade implementation of E z as FIR lossless unitary matrices separated by delays.

If the real-coefficient lossless matrix is denoted by E z ; then the matrix is said to have a special case of lossless degree of one iff it can be characterized by the relationship [18]:

E z = I vv + + z 1 vv + R E26

where R is an arbitrary M × M unitary matrix and v is an M × 1 column vector with unit norm. From Eq. (26), the paraunitary condition for a filter bank is obtained as follows [18]:

I v k v k + + v k v k + z E k z = E k 1 z E27

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.

Advertisement

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 φ n is the series sum of the shifted versions of φ 2 n , then according to [15, 16], φ n can be represented as:

φ n = k h k 2 φ 2 n k , k E28

where h k denotes the scaling coefficients. If n is transformed such that n 2 α n β , then the relationship in Eq. (28) becomes [14]:

φ 2 α n β = k h k 2 φ 2 2 α n β k E29

which translates into:

φ 2 α n β = m = 2 β + k h m 2 β 2 φ 2 α + 1 n m E30

when k = m 2 β .

In a similar consideration to Eq. (28), the wavelet function ψ n can be represented as [19]:

ψ n = k g k 2 φ 2 n k , k E31

where g k denotes the wavelet coefficients. Also, if n is transformed such that n 2 α n β , then the relationship in Eq. (31) becomes [14]:

φ 2 α n β = k g k 2 φ 2 2 α n β k E32

which translates into:

ψ 2 α n β = m = 2 β + k g m 2 β 2 φ 2 α + 1 n m E33

when k = m 2 β .

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 α + 1 scaling space f n V α + 1 , then according to [16, 17], the reconstruction is derived as follows:

f n = 1 M β = λ α + 1 , β φ α + 1 , β n = 1 M β = λ α + 1 , β 2 α + 1 ψ 2 α + 1 n β E34

For the next scale, Eq. (34) becomes:

f n = 1 M β λ α , β 2 α / 2 φ 2 α n β + β γ α , β 2 α / 2 ψ 2 α n β E35

Substituting Eqs. (28) and (31) into Eq. (35) and after algebraic manipulations yields [14]:

λ α + 1 , β = m λ α , β h β 2 m + m γ α , β g β 2 m E36
Advertisement

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 H 0 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.

Figure 10.

State chart for wavelet design procedure.

In the first state in Figure 10, the design problem formulation which can be achieved using trigonometric polynomial, takes the following into consideration [14]:

  1. Compact support which guarantees that the wavelet is characterized by finite non-zero coefficient.

  2. Paraunitary condition which guarantees the generation of orthonormal wavelets.

  3. 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.

Advertisement

7. Conclusion

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.

References

  1. 1. Mertins A. Signal Analysis: Wavelets, Filter Banks, Time-Frequency Transforms and Applications. Sussex: Wiley; 1999
  2. 2. 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
  3. 3. Mallat S. A Wavelet Tour of Signal Processing: The Sparse Way. 3rd ed. Burlington: Academic Press; 2009
  4. 4. 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
  5. 5. Chu E. Discrete and Continuous Fourier Transforms: Analysis, Applications and Fast Algorithms. Boca Raton: CRC Press; 2008
  6. 6. 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
  7. 7. 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
  8. 8. 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
  9. 9. Agrawal SK, Sahu OP. Two channel quadrature mirror filter bank: An overview. Hindawi—ISRN Signal Processing. 2013;2013:1-10
  10. 10. Giron-Sierra JM. Digital Signal Processing with MATLAB Examples. Singapore: Springer Singapore; 2017
  11. 11. Vetterli M, Herley C. Wavelets and filter banks: Theory and design. IEEE Transactions on Signal Processing. 1992;40(9):2207-2231
  12. 12. Vaidyanathan PP. Multirate Systems and Filter Banks. New Jersey: Prentice-Hall; 1993
  13. 13. 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
  14. 14. Nikookar H. Wavelet Radio: Adaptive and Reconfigurable Wireless Systems Based on Wavelets. Cambridge: Cambridge University Press; 2013
  15. 15. Gao X, Nguyen TQ, Strang G. On factorization of M-channel paraunitary filterbanks. IEEE Transactions on Signal Processing. 2001;49(7):1433-1466
  16. 16. 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
  17. 17. 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
  18. 18. 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
  19. 19. Burrus CS, Gopinath RA, Guo H. Introduction to Wavelets and Wavelet Transform: A Primer. Upper Saddle River: Prentice Hall; 1998

Written By

Peter Yusuf Dibal, Elizabeth Onwuka, James Agajo and Caroline Alenoghena

Submitted: 03 December 2018 Reviewed: 07 February 2019 Published: 27 September 2019