<< Chapter < Page Chapter >> Page >
H 0 ( z ) = 1 + z - 1 + ... + z - ( M - 1 ) M K Q ( z ) ,

for maximal possible K . By default, every unitary scaling filter h 0 is one-regular (because n h 0 ( n ) = M - see [link] , Theorem  [link] for equivalent characterizations of K -regularity). Each of the K -identical factors in Eqn.  [link] adds an extra linear constraint on h 0 (actually, it is one linear constraint on each of the M polyphase subsequences of h 0 - see [link] ).

There is no simple relationship between the smoothness of the scaling function and K -regularity. However, the smoothness of the maximally regular scaling filter, h 0 , with fixed filter length N , tends to be an increasing function of N . Perhaps one can argue that K -regularity is an important concept independent of the smoothness of the associated wavelet system. K -regularity implies that the moments of the wavelets vanish up to order K - 1 , and therefore, functions can be better approximated by usingjust the scaling function and its translates at a given scale. Formulae exist for M -band maximally regular K -regular scaling filters (i.e., only the sequence h 0 ) [link] . Using the Householder parameterization, one can then design the remaining filters in the filter bank.

The linear constraints on h 0 that constitute K -regularity become nonexplicit nonlinear constraints on theHouseholder parameterization of the associated filter bank. However, one-regularity can be explicitly incorporated and thisgives a parameterization of all M -band compactly supportedwavelet tight frames. To see, this consider the following two factorizations of H p ( z ) of a unitary filter bank.

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

and

H p T ( z ) = k = K - 1 1 I - Q k + z - 1 Q k W 0 .

Since H p ( 1 ) = V 0 and H p T ( 1 ) = W 0 , V 0 = W 0 T . The first column of W 0 is the unit vector [ H 0 , 0 ( 1 ) , H 0 , 1 ( 1 ) , ... , H 0 , M - 1 ( 1 ) ] T . Therefore,

k = 0 M - 1 H 0 , k ( 1 ) 2 = 1 .

But since n h 0 ( n ) = H 0 ( 1 ) = M ,

k = 0 M - 1 H 0 , k ( 1 ) = M .

Therefore, for all k , H 0 , k ( 1 ) = 1 M . Hence, the first row of V 0 is [ 1 / M , 1 / M , ... , 1 / M ] . In other words, a unitary filter bank gives rise to a WTF iff the first rowof V 0 in the Householder parameterization is the vector with allentries 1 / M .

Alternatively, consider the Given's factorization of H p ( z ) for a two-channel unitary filter bank.

H p ( z ) = i = K - 1 1 cos θ i z - 1 sin θ i - sin θ i z - 1 cos θ i cos θ 0 sin θ 0 - sin θ 0 cos θ 0 .

Since for a WTF we require

1 2 1 2 - 1 2 1 2 = H 0 , 0 ( 1 ) H 1 , 0 ( 1 ) H 0 , 1 ( 1 ) H 1 , 1 ( 1 ) = cos ( Θ ) sin ( Θ ) - sin ( Θ ) cos ( Θ ) ,

we have Θ = k = 0 K - 1 θ k = π 4 . This is the condition for the lattice parameterization tobe associated with wavelets.

Modulated filter banks

Filter bank design typically entails optimization of the filter coefficients to maximize some goodness measure subject to the perfect reconstruction constraint.Being a constrained (or unconstrained for parameterized unitary filter bank design) nonlinear programming problem, numerical optimization leads to localminima, with the problem exacerbated when there are a large number of filter coefficients. To alleviate this problem one can try to impose structuralconstraints on the filters. For example, if [link] is the desired ideal response, one can impose the constraint that all analysis (synthesis)filters are obtained by modulation of a single “prototype” analysis (synthesis) filter. This is the basic idea behind modulated filter banks [link] , [link] , [link] , [link] , [link] , [link] , [link] , [link] . In what follows, we only consider the case where the number of filtersis equal to the downsampling factor; i.e., L = M .

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