<< Chapter < Page Chapter >> Page >

and H is unitary iff H ˘ is unitary. [link] induces a factorization of H ˘ (and hence H ). If V 0 = d i a g ( V 0 ) and

V i = P i I - P i 0 0 P i I - P i , for i 1 , ... , K - 1 ,
H ˘ = i = K - 1 0 V i .

The factors V i , with appropriate modifications, will be used as fundamental building blocks for filter banks forfinite-length signals.

Now consider a finite input signal x = [ x ( 0 ) , x ( 1 ) , , x ( L - 1 ) ] , where L is a multiple of M and let x ˘ = [ x ( M - 1 ) , , x ( 0 ) , x ( M ) , , x ( L - 1 ) , , x ( L - M ) ] . Then, the finite vector d (the output signal) is given by

d = H ˘ x ˘ = def h p ( K - 1 ) ... h p ( 0 ) 0 ... ... ... 0 0 h p ( K - 1 ) ... h p ( 0 ) ... ... 0 h p ( K - 1 ) ... h p ( 0 ) 0 0 ... ... ... 0 h p ( K - 1 ) ... h p ( 0 ) x ˘ .

H ˘ is an ( L - N + M ) × L matrix, where N = M K is the length of the filters. Now since the rows of H ˘ are mutually orthonormal (i.e., has rank L ), one has to append N - M = M ( K - 1 ) rows from the orthogonal complement of H ˘ to make the map from x to an augmented d unitary. To get a complete description of these rows, we turn to thefactorization of H p ( z ) . Define the L × L matrix V 0 = d i a g ( V 0 ) and for i 1 , ... , K - 1 the ( L - M i ) × ( L - M i + M ) matrices

V i = P i I - P i 0 ... ... 0 0 P i I - P i ... ... 0 0 P i I - P i 0 0 ... ... 0 P i I - P i .

then H ˘ is readily verified by induction to be i = K - 1 0 V i . Since each of the factors (except V 0 ) has M more columns than rows, they can be made unitary by appending appropriate rows. Indeed, B i V i C i is unitary where, B j = Υ j T ( I - P j ) 0 ... 0 , and C j = 0 ... 0 Ξ j T P j . . Here Ξ is the δ i × M left unitary matrix that spans the range of the P i ; i.e., P i = Ξ Ξ T , and Υ is the ( M - δ i ) × M left unitary matrix that spans the range of the I - P i ; i.e., I - P i = Υ Υ T . Clearly Υ i Ξ i is unitary. Moreover, if we define T 0 = V 0 and for i 1 , ... , K - 1 ,

T i = I ( M - δ i ) ( i - 1 ) 0 0 0 B i 0 0 V i 0 0 C i 0 0 0 I δ i ( i - 1 )

then each of the factors T i is a square unitary matrix of size L - N + M and

i = K - 1 0 T i = B 1 V 0 H ˘ C 1 V 0 ,

is the unitary matrix that acts on the data. The corresponding unitary matrix that acts on x (rather than x ˘ ) is of the form U H W , where U has M K - M - Δ rows of entry filters in ( K - 1 ) sets given by [link] , while W has Δ rows of exit filters in ( K - 1 ) given by [link] :

Υ j ( I - P j ) h p j ( j - 1 ) J h p j ( j - 2 ) J ... h p j ( 0 ) J ,
Ξ j P j h p j ( j - 1 ) J h p j ( j - 2 ) J ... h p j ( 0 ) J ,

where J is the exchange matrix (i.e., permutation matrix of ones along the anti-diagonal) and

H p j ( z ) = i = j - 1 1 I - P i + z - 1 P i V 0 = def i = 0 j - 1 h p j ( i ) z - i .

The rows of U and W form the entry and exit filters respectively. Clearly they are nonunique.The input/output behavior is captured in

u d w = U H W x .

For example, in the four-coefficient Daubechies' filters in [link] case, there is one entry filter and exit filter.

0 . 8660 0 . 5000 0 0 - 0 . 1294 0 . 2241 0 . 8365 0 . 4830 - 0 . 4830 0 . 8365 - 0 . 2241 - 0 . 1294 0 0 - 0 . 5000 0 . 8660 .

If the input signal is right-sided (i.e., supported in 0 , 1 , ... ), then the corresponding filter bank would only have entry filters. Ifthe filter bank is for left-sided signals one would only have exit filters. Based on the above, we can consider switching betweenfilter banks (that operate on infinite extent input signals). Consider switching from a one-channel to an M channel filter bank. Until instant n = - 1 , the input is the same as the output. At n = 0 , one switches into an M -channel filter bank as quickly as possible. The transition is accomplished by the entry filters(hence the name entry) of the M -channel filter bank. The input/output of this time-varying filter bank is

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