<< Chapter < Page
  Wavelets and wavelet transforms     Page 22 / 28
Chapter >> Page >

Nonsmooth local trigonometric bases

To construct local trigonometric bases we have to choose: (a) the window functions w k ( t ) ; and (b) the trigonometric functions (i.e., α , β and γ in Eq.  [link] ). If we use the rectangular window (which we know is a bad choice),then it suffices to find a trigonometric basis for the interval that the window spans. Without loss of generality, we couldconsider the unit interval ( 0 , 1 ) and hence we are interested in trigonometric bases for L 2 ( ( 0 , 1 ) ) . It is easy to see that the following four sets of functions satisfy this requirement.

  1. Φ n ( t ) = 2 cos ( π ( n + 1 2 ) t ) , n 0 , 1 , 2 , ... ;
  2. Φ n ( t ) = 2 sin ( π ( n + 1 2 ) t ) , n 0 , 1 , 2 , ... ;
  3. Φ n ( t ) = 1 , 2 cos ( π n t ) , n 1 , 2 , ... ;
  4. Φ n ( t ) = 2 sin ( π n t ) , n 0 , 1 , 2 , ... .

Indeed, these orthonormal bases are obtained from the Fourier series on ( - 2 , 2 ) (the first two) and on ( - 1 , 1 ) (the last two) by appropriately imposing symmetries and hence are readily verified to becomplete and orthonormal on ( 0 , 1 ) . If we choose a set of nonoverlapping rectangular window functions w k ( t ) such that k w k ( t ) = 1 for all t , and define χ k , n ( t ) = w k ( t ) Φ n ( t ) , then, χ k , n ( t ) is a local trigonometric basis for L 2 ( ) , for each of the four choices of p h i n ( t ) above.

Construction of smooth windows

We know how to construct orthonormal trigonometric bases for disjoint temporal bins or intervals. Now we need to construct smooth windows w k ( t ) that when applied to cosines and sines retain orthonormality. An outline of the process is as follows:A unitary operation is applied that “unfolds” the discontinuities of all the local basis functions at the boundaries ofeach temporal bin. Unfolding leads to overlapping (unfolded) basis functions. However,since unfolding is unitary, the resulting functions still form an orthonormal basis. The unfolding operator is parameterizedby a function r ( t ) that satisfies an algebraic constraint (which makes the operator unitary). The smoothness of the resulting basisfunctions depends on the smoothness of this underlying function r ( t ) .

The function r ( t ) , referred to as a rising cutoff function, satisfies the following conditions (see [link] ) :

r ( t ) 2 + r ( - t ) 2 = 1 , for all t I R ; r ( t ) = 0 , if t - 1 1 , if t 1

r ( t ) is called a rising cutoff function because it rises from 0 to 1 in the interval [ - 1 , 1 ] (note: it does not necessarily have to be monotone increasing). Multiplying a function by r ( t ) would localize it to [ - 1 , ] . Every real-valued function r ( t ) satisfying [link] is of the form r ( t ) = s i n ( θ ( t ) ) where

θ ( t ) + θ ( - t ) = π 2 for all t I R ; r ( t ) = 0 , if t - 1 . π 2 , if t 1 .

This ensures that r ( - t ) = sin ( θ ( - t ) ) = sin ( π 2 - θ ( t ) ) = cos ( θ ( t ) ) and therefore r 2 ( t ) + r 2 ( - t ) = 1 . One can easily construct arbitrarily smooth risingcutoff functions. We give one such recipe from [link] (p.105) . Start with a function

r [ 0 ] ( t ) = 0 , if t - 1 sin π 4 ( 1 + t ) , if - 1 < t < 1 1 , if t 1

It is readily verified to be a rising cutoff function. Now recursively define r [ 1 ] ( t ) , r [ 2 ] ( t ) , ... as follows:

r [ n + 1 ] ( t ) = r [ n ] ( sin ( π 2 t ) ) .

Notice that r [ n ] ( t ) is a rising cutoff function for every n . Moreover, by induction on n it is easy to show that r [ n ] ( t ) C 2 n - 1 (it suffices to show that derivatives at t = - 1 and t = 1 exist and are zero up to order 2 n - 1 ).

Folding and unfolding

Using a rising cutoff function r ( t ) one can define the folding operator, U , and its inverse, the unfolding operator U as follows:

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