<< Chapter < Page Chapter >> Page >

where θ 0 = π 3 and θ 1 = - π 12 . The fact that the filter bank is associated with waveletsis precisely because θ 0 + θ 1 = π 4 . More generally, for a filter bank with filters of length 2 K to be associated with wavelets, k = 0 K - 1 θ k = π 4 . This is expected since for filters of length 2 K to be associated with wavelets we have seen (from the Householder factorization)that there are K - 1 parameters v k . Our second example belongs to a class of unitary filter banks called modulated filter banks , which is described in a following section. A Type 1 modulated filter bank with filters of length N = 2 M and associated with a wavelet orthonormal basis is defined by

h i ( n ) = 1 2 M sin π ( i + 1 ) ( n + . 5 ) M - ( 2 i + 1 ) π 4 - sin π i ( n + . 5 ) M - ( 2 i + 1 ) π 4 ,

where i 0 , ... , M - 1 and n 0 , ... , 2 M - 1 [link] , [link] . Consider a three-band example with length six filters.In this case, K = 2 , and therefore one has one projection P 1 and the matrix V 0 . The projection is one-dimensional and given by the Householder parameter

v 1 T = 1 6 1 - 2 1 and V 0 = 1 3 1 1 1 - 3 + 1 2 1 3 - 1 2 3 - 1 2 1 - 3 + 1 2 .

The third example is another Type 1 modulated filter bank with M = 4 and N = 8 . The filters are given in [link] . H p ( z ) had the following factorization

H p ( z ) = I - P 1 + z - 1 P 1 V 0 ,

where P 1 is a two-dimensional projection P 1 = v 1 v 1 T + v 2 v 2 T (notice the arbitrary choice of v 1 and v 2 ) given by

v 1 = 0 . 41636433418450 - 0 . 78450701561376 0 . 32495344564406 0 . 32495344564406 , v 2 = 0 . 00000000000000 - 0 . 14088210492943 0 . 50902478635868 - 0 . 84914427477499

and

V 0 = 1 2 1 1 1 1 - 2 0 2 0 0 2 0 - 2 1 - 1 1 - 1 .

Notice that there are infinitely many choices of v 1 and v 2 that give rise to the same projection P 1 .

M -band wavelet tight frames

In [link] , Theorem 7 , while discussing the properties of M -band wavelet systems, we saw that the lowpass filter h 0 ( h in the notation used there) must satisfy the linear constraint n h 0 ( n ) = M . Otherwise, a scaling function with nonzero integral could not exit.It turns out that this is precisely the only condition that an FIR unitary filter bank has to satisfy in order for it to generatean M -band wavelet system [link] , [link] . Indeed, if this linear constraint is not satisfied the filter bank does not generate a wavelet system. This single linear constraint (for unitary filter banks) also impliesthat h h i ( n ) = 0 for i 1 , 2 , ... , M - 1 (because of Eqn.  [link] ). We now give the precise result connecting FIR unitary filter banks and wavelet tight frames.

Theorem 46 Given an FIR unitary filter bank with n h 0 ( n ) = M , there exists an unique, compactly supported, scaling function ψ 0 ( t ) L 2 ( ) (with support in [ 0 , N - 1 M - 1 ] , assuming h 0 is supported in [ 0 , N - 1 ] ) determined by the scaling recursion :

ψ 0 ( t ) = M k h 0 ( k ) ψ 0 ( M t - k ) .

Define wavelets, ψ i ( t ) ,

ψ i ( t ) = M k h i ( k ) ψ 0 ( M t - k ) i 1 , 2 , ... , M - 1 ,

and functions, ψ i , j , k ( t ) ,

ψ i , j , k ( t ) = M j / 2 ψ i ( M j t - k ) .

Then ψ i , j , k forms a tight frame for L 2 ( ) . That is, for all f L 2 ( )

f ( t ) = i = 1 M - 1 j , k = - f , ψ i , j , k ψ i , j , k ( t ) .

Also,

f ( t ) = k f , ψ 0 , 0 , k ψ 0 , 0 , k ( t ) + i = 1 M - 1 j = 1 k = - f , ψ i , j , k ψ i , j , k ( t ) .

Remark: A similar result relates general FIR (not necessarily unitary) filter banks and M -band wavelet frames [link] , [link] , [link] .

Starting with [link] , one can calculate the scaling function using either successive approximationor interpolation on the M -adic rationals—i.e., exactly as in the two-band case in Chapter  [link] . Equation  [link] then gives the wavelets in terms of the scaling function. As in the two-band case, the functions ψ i ( t ) , so constructed, invariably turn out highly irregular and sometimes fractal. Thesolution, once again, is to require that several moments of the scaling function (or equivalently the moments of the scaling filter h 0 ) are zero. This motivates the definition of K -regular M -band scaling filters: A unitary scaling filter h 0 is said to be K regular if its Z -transform can be written in the form

Get Jobilize Job Search Mobile App in your pocket Now!

Get it on Google Play Download on the App Store Now




Source:  OpenStax, Wavelets and wavelet transforms. OpenStax CNX. Aug 06, 2015 Download for free at https://legacy.cnx.org/content/col11454/1.6
Google Play and the Google Play logo are trademarks of Google Inc.

Notification Switch

Would you like to follow the 'Wavelets and wavelet transforms' conversation and receive update notifications?

Ask