<< Chapter < Page Chapter >> Page >

If we assume the length of the sequence of the signal is L and the length of the sequence of scaling filter coefficients h ( n ) is N , then the number of multiplications necessary to calculate each scaling functionand wavelet expansion coefficient at the next scale, c ( J 1 - 1 , k ) and d ( J 1 - 1 , k ) , from the samples of the signal, f ( T k ) c ( J 1 , k ) , is L N . Because of the downsampling, only half are needed to calculate thecoefficients at the next lower scale, c ( J 2 - 1 , k ) and d ( J 2 - 1 , k ) , and repeats until there is only one coefficient at a scale of j = J 0 . The total number of multiplications is, therefore,

Mult = L N + L N / 2 + L N / 4 + + N
= L N ( 1 + 1 / 2 + 1 / 4 + + 1 / L ) = 2 N L - N

which is linear in L and in N . The number of required additions is essentially the same.

If the length of the signal is very long, essentially infinity, the coarsest scale J 0 must be determined from the goals of the particular signal processing problem being addressed. For this case, the number ofmultiplications required per DWT coefficient or per input signal sample is

Mult/sample = N ( 2 - 2 - J 0 )

Because of the relationship of the scaling function filter h ( n ) and the wavelet filter h 1 ( n ) at each scale (they are quadrature mirror filters), operations can be shared between them through the use of alattice filter structure, which will almost halve the computational complexity. That is developed in Chapter: Filter Banks and Transmultiplexers and [link] .

The periodic case

In many practical applications, the signal is finite in length (finite support) and can be processed as single “block," much as the FourierSeries or discrete Fourier transform (DFT) does. If the signal to be analyzed is finite in length such that

f ( t ) = 0 t < 0 0 t > P f ( t ) 0 < t < P

we can construct a periodic signal f ˜ ( t ) by

f ˜ ( t ) = n f ( t + P n )

and then consider its wavelet expansion or DWT. This creation of a meaningful periodic function can still be done, even if f ( t ) does not have finite support, if its energy is concentrated and some overlap isallowed in [link] .

Periodic Property 1: If f ˜ ( t ) is periodic with integer period P such that f ˜ ( t ) = f ˜ ( t + P n ) , then the scaling function and wavelet expansion coefficients (DWT terms) at scale J are periodic with period 2 J P .

If f ˜ ( t ) = f ˜ ( t + P ) then d ˜ j ( k ) = d ˜ j ( k + 2 j P )

This is easily seen from

d ˜ j ( k ) = - f ˜ ( t ) ψ ( 2 j t - k ) d t = f ˜ ( t + P n ) ψ ( 2 j t - k ) d t

which, with a change of variables, becomes

= f ˜ ( x ) ψ ( 2 j ( x - P n ) - k ) d x = f ˜ ( x ) ψ ( 2 j x - ( 2 j P n + k ) ) d x = d ˜ j ( k + 2 j P n )

and the same is true for c ˜ j ( k ) .

Periodic Property 2: The scaling function and wavelet expansion coefficients (DWT terms) can be calculated from the inner product of f ˜ ( t ) with φ ( t ) and ψ ( t ) or, equivalently, from the inner product of f ( t ) with the periodized φ ˜ ( t ) and ψ ˜ ( t ) .

c ˜ j ( k ) = f ˜ ( t ) , φ ( t ) = f ( t ) , φ ˜ ( t )

and

d ˜ j ( k ) = f ˜ ( t ) , ψ ( t ) = f ( t ) , ψ ˜ ( t )

where φ ˜ ( t ) = n φ ( t + P n ) and ψ ˜ ( t ) = n ψ ( t + P n ) .

This is seen from

d ˜ j ( k ) = - f ˜ ( t ) ψ ( 2 j t - k ) d t = n 0 P f ( t ) ψ ( 2 j ( t + P n ) - k ) d t = 0 P f ( t ) n ψ ( 2 j ( t + P n ) - k ) d t
d ˜ j ( k ) = 0 P f ( t ) ψ ˜ ( 2 j t - k ) d t

where ψ ˜ ( 2 j t - k ) = n ψ ( 2 j ( t + P n ) - k ) is the periodized scaled wavelet.

Periodic Property 3: If f ˜ ( t ) is periodic with period P , then Mallat's algorithm for calculating the DWT coefficients in [link] becomes

c ˜ j ( k ) = m h ( m - 2 k ) c ˜ j + 1 ( 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