Date: Wed, 28 Apr 1999 15:07:29 -0500
To: ee648on, ee648off
From: Tom Krauss
Subject: Note on Homework 7, part (d)
At 2:41 PM -0500 4/26/99, Iraj Manocheri wrote:
>Dear Professor,
>
>There is similar problem with homework 7 as was with homework 6 for Part D.
>In this case complex data is used for the impulse response of the QMF
>filters.
>In problem 6 perfect reconstruction was possible if only the imaginary part
>of the
>filter output was considered. In problem 7 with M = 4, there seems to be no
>perfect
>reconstruction for y complex, y real, or y imaginary.
>
>Please give me some explanation of the reason and under what circumstances
>is
>PR possible.
Many of us have been asking this question and with the help
of some students (I heard of a fix from Terry Charbonneau
who apparently got help from another student), I have
figured out the problem.
The issue with part (d) (in both homework 6 and this one) is
that the 2 channel QMF bank with h0 = [1 j], h1 = [1 -j],
g0 = [1 j], g1 = [-1 j] does NOT give simple perfect reconstruction,
rather it gives
y[n] = 2j * x[n-1]
We can remove this j factor by dividing it out in the g0 and
g1, so that g0 = -j*h0 and g1 = j*h1.
When we do so, the resulting nonuniform filter bank becomes perfect
reconstruction.
Regards,
Tom
PS I noticed in Zoltowski's code PRNonuniform
you need to change TWO lines to make the fix:
g0=-j*h00; g1=j*h11; <--- redefine g0,g1
...
g3=j*h11; <--- also redefine g3.