<< Chapter < Page Chapter >> Page >

In this appendix we develop most of the results on scaling functions, wavelets and scaling and wavelet coefficients presented in [link] and elsewhere. For convenience, we repeat [link] , [link] , [link] , and [link] here

φ ( t ) = n h ( n ) 2 φ ( 2 t - n )
n h ( n ) = 2
φ ( t ) φ ( t - k ) d t = E δ ( k ) = E i f k = 0 0 o t h e r w i s e

If normalized

φ ( t ) φ ( t - k ) d t = δ ( k ) = 1 i f k = 0 0 o t h e r w i s e

The results in this appendix refer to equations in the text writtenin bold face fonts.

Equation [link] is the normalization of [link] and part of the orthonormal conditions required by [link] for k = 0 and E = 1 .

Equation [link] If the φ ( x - k ) are orthogonal, [link] states

φ ( x + m ) φ ( x ) d x = E δ ( m )

Summing both sides over m gives

m φ ( x + m ) φ ( x ) d x = E

which after reordering is

φ ( x ) m φ ( x + m ) d x = E .

Using [link] , [link] , and [link] gives

φ ( x ) d x A 0 = E

but φ ( x ) d x = A 0 from [link] , therefore

A 0 2 = E

If the scaling function is not normalized to unity, one can show the more general result of [link] . This is done by noting that a more general form of [link] is

m φ ( x + m ) = φ ( x ) d x

if one does not normalize A 0 = 1 in [link] through [link] .

Equation [link] follows from summing [link] over m as

m φ ( x + m ) φ ( x ) d x = φ ( x ) 2 d x

which after reordering gives

φ ( x ) m φ ( x + m ) d x = φ ( x ) 2 d x

and using [link] gives [link] .

Equation [link] is derived by applying the basic recursion equation to its own right hand side to give

φ ( t ) = n h ( n ) 2 k h ( k ) 2 φ ( 2 ( 2 t - n ) - k )

which, with a change of variables of = 2 n + k and reordering of operation, becomes

φ ( t ) = n h ( n ) h ( - 2 n ) 2 φ ( 4 t - ) .

Applying this j times gives the result in [link] . A similar result can be derived for the wavelet.

Equation [link] is derived by defining the sum

A J = φ ( 2 J )

and using the basic recursive equation [link] to give

A J = n h ( n ) 2 φ ( 2 2 J - n ) .

Interchanging the order of summation gives

A J = 2 n h ( n ) φ ( 2 J - 1 - n )

but the summation over is independent of an integer shift so that using [link] and [link] gives

A J = 2 2 n h ( n ) φ ( 2 J - 1 ) = 2 A J - 1 .

This is the linear difference equation

A J - 2 A J - 1 = 0

which has as a solution the geometric sequence

A J = A 0 2 J .

If the limit exists, equation [link] divided by 2 J is the Riemann sum whose limit is the definition of the Riemann integral of φ ( x )

lim J A J 1 2 J = φ ( x ) d x = A 0 .

It is stated in [link] and shown in [link] that if φ ( x ) is normalized, then A 0 = 1 and [link] becomes

A J = 2 J .

which gives [link] .

Equation [link] shows another remarkable property of φ ( x ) in that the bracketed term is exactly equal to the integral, independent of J . No limit need be taken!

Equation [link] is the “partitioning of unity" by φ ( x ) . It follows from [link] by setting J = 0 .

Equation [link] is generalization of [link] by noting that the sum in [link] is independent of a shift of the form

φ ( 2 J - L 2 M ) = 2 J

for any integers M J and L . In the limit as M , L 2 M can be made arbitrarily close to any x , therefore, if φ ( x ) is continuous,

φ ( 2 J - x ) = 2 J .

This gives [link] and becomes [link] for J = 0 . Equation [link] is called a “partitioning of unity" for obvious reasons.

The first four relationships for the scaling function hold in a generalized form for the more general defining equation [link] . Only [link] is different. It becomes

k φ ( k M J ) = M J

for M an integer. It may be possible to show that certain rational M are allowed.

Equations [link] , [link] , and [link] are the recursive relationship for the Fourier transform of the scaling function and areobtained by simply taking the transform [link] of both sides of [link] giving

Φ ( ω ) = n h ( n ) 2 φ ( 2 t - n ) e - j ω t d t

which after the change of variables y = 2 t - n becomes

Φ ( ω ) = 2 2 n h ( n ) φ ( y ) e - j ω ( y + n ) / 2 d y

and using [link] gives

Φ ( ω ) = 1 2 n h ( n ) e - j ω n / 2 φ ( y ) e - j ω y / 2 d y = 1 2 H ( ω / 2 ) Φ ( ω / 2 )

which is [link] and [link] . Applying this recursively gives the infinite product [link] which holds for any normalization.

Equation [link] states that the sum of the squares of samples of the Fourier transform of the scaling function is one if the samples areuniform every 2 π . An alternative derivation to that in Appendix A is shown here by taking the definition of the Fourier transform of φ ( x ) , sampling it every 2 π k points and multiplying it times its complex conjugate.

Φ ( ω + 2 π k ) Φ ( ω + 2 π k ) ¯ = φ ( x ) e - j ( ω + 2 π k ) x d x φ ( y ) e j ( ω + 2 π k ) y d y

Summing over k gives

k | Φ ( ω + 2 π k ) | 2 = k φ ( x ) φ ( y ) e - j ω ( x - y ) e - j 2 π k ( x - y ) d x d y
= φ ( x ) φ ( y ) e j ω ( y - x ) k e j 2 π k ( y - x ) d x d y
= φ ( x ) φ ( x + z ) e j ω z k e j 2 π k z d x d z

but

k e j 2 π k z = δ ( z - )

therefore

k | Φ ( ω + 2 π k ) | 2 = φ ( x ) φ ( x + ) e - j ω d x

which becomes

φ ( x ) φ ( x + ) d x e j ω

Because of the orthogonality of integer translates of φ ( x ) , this is not a function of ω but is | φ ( x ) | 2 d x which, if normalized, is unity as stated in [link] . This is the frequency domain equivalent of [link] .

Equations [link] and [link] show how the scaling function determines the equation coefficients. This is derived by multiplying bothsides of [link] by φ ( 2 x - m ) and integrating to give

φ ( x ) φ ( 2 x - m ) d x = n h ( n ) φ ( 2 x - n ) φ ( 2 x - m ) d x
= 1 2 n h ( n ) φ ( x - n ) φ ( x - m ) d x .

Using the orthogonality condition [link] gives

φ ( x ) φ ( 2 x - m ) d x = h ( m ) 1 2 | φ ( y ) | 2 d y = 1 2 h ( m )

which gives [link] . A similar argument gives [link] .

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