<< Chapter < Page Chapter >> Page >
A ( ω ) = φ ( x ) e - j ω x d x φ ( y ) e j ω y d y
= Φ ( ω ) Φ ( - ω ) = | Φ ( ω ) | 2

If we look at [link] as being the inverse Fourier transform of [link] and sample a ( t ) at t = k , we have

a ( k ) = 1 2 π - | Φ ( ω ) | 2 e j ω k d ω
= 1 2 π 0 2 π | Φ ( ω + 2 π ) | 2 e j ω k d ω = 1 2 π 0 2 π | Φ ( ω + 2 π ) | 2 e j ω k d ω

but this integral is the form of an inverse discrete-time Fourier transform (DTFT) which means

, a ( k ) e j ω k = | Φ ( ω + 2 π ) | 2 .

If the integer translates of φ ( t ) are orthogonal, a ( k ) = δ ( k ) and we have our result

| Φ ( ω + 2 π ) | 2 = 1 .

If the scaling function is not normalized

| Φ ( ω + 2 π ) | 2 = | φ ( t ) | 2 d t

which is similar to Parseval's theorem relating the energy in the frequency domain to the energy in the time domain.

Proof 6 Equation [link] states a very interesting property of the frequency response of an FIR filter with the scaling coefficients as filtercoefficients. This result can be derived in the frequency or time domain. We will show the frequency domain argument. The scaling equation [link] becomes [link] in the frequency domain. Taking the squared magnitude of both sides of a scaled version of

Φ ( ω ) = 1 2 H ( ω / 2 ) Φ ( ω / 2 )

gives

| Φ ( 2 ω ) | 2 = 1 2 | H ( ω ) | 2 | Φ ( ω ) | 2

Add k π to ω and sum over k to give for the left side of [link]

k | Φ ( 2 ω + 2 π k ) | 2 = K = 1

which is unity from [link] . Summing the right side of [link] gives

k 1 2 | H ( ω + k π ) | 2 | Φ ( ω + k π ) | 2

Break this sum into a sum of the even and odd indexed terms.

k 1 2 | H ( ω + 2 π k ) | 2 | Φ ( ω + 2 π k ) | 2 + k 1 2 | H ( ω + ( 2 k + 1 ) π ) | 2 | Φ ( ω + ( 2 k + 1 ) π ) | 2
= 1 2 | H ( ω ) | 2 k | Φ ( ω + 2 π k ) | 2 + 1 2 | H ( ω + π ) | 2 k | Φ ( ω + ( 2 k + 1 ) π ) | 2

which after using [link] gives

= 1 2 | H ( ω ) | 2 + 1 2 | H ( ω + π ) | 2 = 1

which gives [link] . This requires both the scaling and orthogonal relations but no specific normalization of φ ( t ) . If viewed as an FIR filter, h ( n ) is called a quadrature mirror filter (QMF) because of the symmetry of its frequencyresponse about π .

Proof 10 The multiresolution assumptions in [link] require the scaling function and wavelet satisfy [link] and [link]

φ ( t ) = n h ( n ) 2 φ ( 2 t - n ) , ψ ( t ) = n h 1 ( n ) 2 φ ( 2 t - n )

and orthonormality requires

φ ( t ) φ ( t - k ) d t = δ ( k )

and

ψ ( t ) φ ( t - k ) d t = 0

for all k Z . Substituting [link] into [link] gives

n h 1 ( n ) 2 φ ( 2 t - n ) h ( ) 2 φ ( 2 t - 2 k - ) d t = 0

Rearranging and making a change of variables gives

n , h 1 ( n ) h ( ) 1 2 φ ( y - n ) φ ( y - 2 k - ) d y = 0

Using [link] gives

n , h 1 ( n ) h ( ) δ ( n - 2 k - ) = 0

for all k Z . Summing over gives

n h 1 ( n ) h ( n - 2 k ) = 0

Separating [link] into even and odd indices gives

m h 1 ( 2 m ) h ( 2 m - 2 k ) + h 1 ( 2 + 1 ) h ( 2 + 1 - 2 k ) = 0

which must be true for all integer k . Defining h e ( n ) = h ( 2 n ) , h o ( n ) = h ( 2 n + 1 ) and g ˜ ( n ) = g ( - n ) for any sequence g , this becomes

h e h ˜ 1 e + h o h ˜ 1 o = 0 .

From the orthonormality of the translates of φ and ψ one can similarly obtain the following:

h e h ˜ e + h o h ˜ o = δ .
h 1 e h ˜ 1 e + h 1 o h ˜ 1 o = δ .

This can be compactly represented as

h e h o h 1 e h 1 o h ˜ e h ˜ 1 e h ˜ o h ˜ 1 o = δ 0 0 δ .

Assuming the sequences are finite length [link] can be used to show that

h e h 1 o - h o h 1 e = ± δ k ,

where δ k ( n ) = δ ( n - k ) . Indeed, taking the Z-transform of [link] we get using the notation of  Chapter: Filter Banks and Transmultiplexers H p ( z ) H p T ( z - 1 ) = I . Because, the filters are FIR H p ( z ) is a (Laurent) polynomial matrix with a polynomial matrix inverse. Therefore the determinant of H p ( z ) is of the form ± z k for some integer k . This is equivalent to [link] . Now, convolving both sides of [link] by h ˜ e we get

± h ˜ e δ k = h e h 1 o - h o h 1 e h ˜ e = h e h ˜ e h 1 o - h 1 e h ˜ e h o = h e h ˜ e h 1 o + h 1 o h ˜ o h o = h e h ˜ e + h o h ˜ o h 1 o = h 1 o .

Similarly by convolving both sides of [link] by h ˜ o we get

h ˜ o δ k = h 1 e .

Combining [link] and [link] gives the result

h 1 ( n ) = ± ( - 1 ) n h ( - n + 1 - 2 k ) .

Proof 11 We show the integral of the wavelet is zero by integrating both sides of ( [link] b) gives

ψ ( t ) d t = n h 1 ( n ) 2 φ ( 2 t - n ) d t

But the integral on the right hand side is A 0 , usually normalized to one and from [link] or [link] and [link] we know that

n h 1 ( n ) = 0

and, therefore, from [link] , the integral of the wavelet is zero.

The fact that multiplying in the time domain by ( - 1 ) n is equivalent to shifting in the frequency domain by π gives H 1 ( ω ) = H ( ω + π ) .

Get Jobilize Job Search Mobile App in your pocket Now!

Get it on Google Play Download on the App Store Now




Source:  OpenStax, Intermodular linking test collection. OpenStax CNX. Sep 09, 2015 Download for free at http://legacy.cnx.org/content/col11841/1.4
Google Play and the Google Play logo are trademarks of Google Inc.

Notification Switch

Would you like to follow the 'Intermodular linking test collection' conversation and receive update notifications?

Ask