<< Chapter < Page Chapter >> Page >

Although the DFT filterbanks are widely used, there is a problem with aliasing in the decimated channels. At first glance, onemight think that this is an insurmountable problem and must simply be accepted. Clearly, with FIR filters and maximaldecimation, aliasing will occur. However, a simple example will show that it is possible to exactly cancel out aliasing under certain conditions!!!

Consider the following trivial filterbank system, with two channels. ( [link] )

Note x ^ n x n with no error whatsoever, although clearly aliasing occurs in both channels! Note that the overall data rate isstill the Nyquist rate, so there are clearly enough degrees of freedom available to reconstruct the data, if the filterbank isdesigned carefully. However, this isn't splitting the data into separate frequency bands, so one questions whether somethingother than this trivial example could work.

Let's consider a general two-channel filterbank, and try to determine conditions under which aliasing can be cancelled, andthe signal can be reconstructed perfectly ( [link] ).

Let's derive x ^ n , using z-transforms, in terms of the components of this system. Recall ( [link] ) is equivalent to Y z H z X z Y ω H ω X ω
and note that ( [link] ) is equivalent to Y z m x m z L m x z L Y ω X L ω
and ( [link] ) is equivalent to Y z 1 M k 0 M 1 X z 1 M W M k Y ω 1 M k 0 M 1 X ω M 2 k M
Y z is derived in the downsampler as follows: Y z m x M m z m Let n M m and m n M , then Y z n x n p δ n M p z n M

Now

x n p δ n M p IDFT x ω 2 M k 0 M 1 δ ω 2 k M IDFT 2 M k 0 M 1 X ω 2 k M W M 2 M 1 M k 0 M 1 X n W M - n k
so
Y z n 1 M k 0 M 1 x n W M - n k z n M 1 M k 0 M 1 x n W M + k z 1 M n 1 M k 0 M 1 X z 1 M W M k
Armed with these results, let's determine X ^ z x ^ n . ( [link] )
Note U 1 z X z H 0 z U 2 z 1 2 k 0 1 X z 1 2 2 k 2 H 0 z 1 2 k 1 2 X z 1 2 H 0 z 1 2 1 2 X z 1 2 H 0 z 1 2 U 3 z 1 2 X z H 0 z 1 2 X z H 0 z U 4 z 1 2 F 0 z H 0 z X z 1 2 F 0 z H 0 z X z and L 4 z 1 2 F 1 z H 1 z X z 1 2 F 1 z H 1 z X z 1 2 F 1 z H 1 z X z 1 2 F 1 z H 1 z X z Finally then,
X ^ z U 4 z L 4 z 1 2 H 0 z F 0 z X z H 0 z F 0 z X z H 1 z F 1 z X z H 1 z F 1 z X z 1 2 H 0 z F 0 z H 1 z F 1 z X z 1 2 H 0 z F 0 z H 1 z F 1 z X z
Note that the X z X ω corresponds to the aliasing terms!

There are four things we would like to have:

  1. No aliasing distortion
  2. No phase distortion (overall linear phase → simple time delay)
  3. No amplitude distortion
  4. FIR filters

No aliasing distortion

By insisting that H 0 z F 0 z H 1 z F 1 z 0 , the X z component of X ^ z can be removed, and all aliasing will be eliminated! There may be many choices for H 0 , H 1 , F 0 , F 1 that eliminate aliasing, but most research has focused on the choice F 0 z H 1 z : F 1 z H 0 z We will consider only this choice in the following discussion.

Phase distortion

The transfer function of the filter bank, with aliasing cancelled, becomes T z 1 2 H 0 z F 0 z H 1 z F 1 z , which with the above choice becomes T z 1 2 H 0 z H 1 z H 1 z H 0 z . We would like T z to correspond to a linear-phase filter to eliminate phase distortion: Call P z H 0 z H 1 z Note that T z 1 2 P z P z Note that P z 1 n p n , and that if p n is a linear-phase filter, 1 n p n is also (perhaps of the opposite symmetry). Also note that the sum of two linear-phase filters of the same symmetry( i.e. , length of p n must be odd ) is also linear phase, so if p n is an odd-length linear-phase filter, there will be no phase distortion. Also note that Z p z p z p n 1 n p n 2 p n n is odd 0 n is even means p n 0 , when n is even. If we choose h 0 n and h 1 n to be linear phase, p n will also be linear phase. Thus by choosing h 0 n and h 1 n to be FIR linear phase, we eliminate phase distortion and get FIR filters as well (condition 4).

Amplitude distortion

Assuming aliasing cancellation and elimination of phase distortion, we might also desire noamplitude distortion ( T ω 1 ). All of these conditions require T z 1 2 H 0 z H 1 z H 1 z H 0 z c z D where c is some constant and D is a linear phase delay. c 1 for T ω 1 . It can be shown by considering that the following can be satisfied! T z P z P z 2 c z D 2 p z 2 c δ n D n is odd p n anything n is even Thus we require P z n 0 N p 2 n z 2 n z D Any factorization of a P z of this form, P z A z B z can lead to a Perfect Reconstruction filter bank of the form H 0 z A z H 1 z B z [This result is attributed to Vetterli.] A well-known specialcase (Smith and Barnwell) H 1 z z 2 D 1 H 0 z Design techniques exist for optimally choosing the coefficients of these filters, under all of these constraints.

Quadrature mirror filters

H 1 z H 0 z H 1 ω H 0 ω H 0 * ω
for real-valued filters. The frequency response is "mirrored" around ω 2 . This choice leads to T z H 0 z 2 H 0 z 2 : it can be shown that this can be a perfect reconstruction system only if H 0 z c 0 z 2 n 0 c 1 z 2 n 1 which isn't a very flexible choice of filters, and not a very good lowpass! The Smith and Barnwell approach is more commonlyused today.

Questions & Answers

A golfer on a fairway is 70 m away from the green, which sits below the level of the fairway by 20 m. If the golfer hits the ball at an angle of 40° with an initial speed of 20 m/s, how close to the green does she come?
Aislinn Reply
cm
tijani
what is titration
John Reply
what is physics
Siyaka Reply
A mouse of mass 200 g falls 100 m down a vertical mine shaft and lands at the bottom with a speed of 8.0 m/s. During its fall, how much work is done on the mouse by air resistance
Jude Reply
Can you compute that for me. Ty
Jude
what is the dimension formula of energy?
David Reply
what is viscosity?
David
what is inorganic
emma Reply
what is chemistry
Youesf Reply
what is inorganic
emma
Chemistry is a branch of science that deals with the study of matter,it composition,it structure and the changes it undergoes
Adjei
please, I'm a physics student and I need help in physics
Adjanou
chemistry could also be understood like the sexual attraction/repulsion of the male and female elements. the reaction varies depending on the energy differences of each given gender. + masculine -female.
Pedro
A ball is thrown straight up.it passes a 2.0m high window 7.50 m off the ground on it path up and takes 1.30 s to go past the window.what was the ball initial velocity
Krampah Reply
2. A sled plus passenger with total mass 50 kg is pulled 20 m across the snow (0.20) at constant velocity by a force directed 25° above the horizontal. Calculate (a) the work of the applied force, (b) the work of friction, and (c) the total work.
Sahid Reply
you have been hired as an espert witness in a court case involving an automobile accident. the accident involved car A of mass 1500kg which crashed into stationary car B of mass 1100kg. the driver of car A applied his brakes 15 m before he skidded and crashed into car B. after the collision, car A s
Samuel Reply
can someone explain to me, an ignorant high school student, why the trend of the graph doesn't follow the fact that the higher frequency a sound wave is, the more power it is, hence, making me think the phons output would follow this general trend?
Joseph Reply
Nevermind i just realied that the graph is the phons output for a person with normal hearing and not just the phons output of the sound waves power, I should read the entire thing next time
Joseph
Follow up question, does anyone know where I can find a graph that accuretly depicts the actual relative "power" output of sound over its frequency instead of just humans hearing
Joseph
"Generation of electrical energy from sound energy | IEEE Conference Publication | IEEE Xplore" ***ieeexplore.ieee.org/document/7150687?reload=true
Ryan
what's motion
Maurice Reply
what are the types of wave
Maurice
answer
Magreth
progressive wave
Magreth
hello friend how are you
Muhammad Reply
fine, how about you?
Mohammed
hi
Mujahid
A string is 3.00 m long with a mass of 5.00 g. The string is held taut with a tension of 500.00 N applied to the string. A pulse is sent down the string. How long does it take the pulse to travel the 3.00 m of the string?
yasuo Reply
Who can show me the full solution in this problem?
Reofrir Reply
Got questions? Join the online conversation and get instant answers!
Jobilize.com Reply

Get Jobilize Job Search Mobile App in your pocket Now!

Get it on Google Play Download on the App Store Now




Source:  OpenStax, Digital signal processing: a user's guide. OpenStax CNX. Aug 29, 2006 Download for free at http://cnx.org/content/col10372/1.2
Google Play and the Google Play logo are trademarks of Google Inc.

Notification Switch

Would you like to follow the 'Digital signal processing: a user's guide' conversation and receive update notifications?

Ask