Sign in

username:

password:



Not a member?

Search Online Books



Search tips

Free Online Books

Ads

Chapters

Chapter Contents:

Search Spectral Audio Signal Processing

  

Book Index | Global Index


Would you like to be notified by email when Julius Orion Smith III publishes a new entry into his blog?

  

Quadrature Mirror Filterbanks (QMF)

The well studied subject of Quadrature Mirror Filters (QMF) is entered by imposing the following symmetry constraint on the analysis filters:

$\displaystyle H_1(z) = H_0(-z)\quad \hbox{(QMF Symmetry Constraint)} \protect$ (12.6)

That is, the filter for channel 1 is constrained to be a $ \pi$-rotation of filter 0 along the unit circle. In the time domain, $ h_1(n) = (-1)^n h_0(n)$, i.e., all odd-index coefficients are negated. If $ H_0$ is a lowpass filter cutting off near $ \omega=\pi/2$ (as is typical), then $ H_1$ is a complementary highpass filter. The exact cut-off frequency can be adjusted along with the roll-off rate to provide a maximally constant frequency-response sum.

Two-channel QMFs have been around since at least 1976 [47], and appear to be the first critically sampled perfect reconstruction filter banks. Historically, the term QMF applied only to two-channel filter banks having the QMF symmetry constraint (11.6). Today, the term ``QMF filter bank'' may refer to more general PR filter banks with any number of channels and not obeying (11.6) [264].

Combining the QMF symmetry constraint with the aliasing-cancellation constraints, given by

\begin{eqnarray*}
F_0(z) &=& \quad\! H_1(-z) = \quad\! H_0(z)\\ [0.1in]
F_1(z) &=& -H_0(-z) = -H_1(z),
\end{eqnarray*}

the perfect reconstruction requirement reduces to

$\displaystyle \hbox{constant}$ $\displaystyle =$ $\displaystyle H_0(z)F_0(z) + H_1(z)F_1(z) = H_0^2(z) - H_0^2(-z)$  
    $\displaystyle \hbox{(QMF Perfect Reconstruction Constraint)}
\protect$ (12.7)

Now, all four filters are determined by $ H_0(z)$.

It is easy to show using the polyphase representation of $ H_0(z)$ (see [264]) that the only causal FIR QMF analysis filters yielding exact perfect reconstruction are two-tap FIR filters of the form

\begin{eqnarray*}
H_0(z) &=& c_0 z^{-2n_0} + c_1 z^{-(2n_1+1)}\\
H_1(z) &=& c_0 z^{-2n_0} - c_1 z^{-(2n_1+1)}
\end{eqnarray*}

where $ c_0$ and $ c_1$ are constants, and $ n_0$ and $ n_1$ are integers. Therefore, only weak channel filters are available in the QMF case ( $ H_1(z)=H_0(-z)$), as we saw in the amplitude-complementary case. On the other hand, very high quality IIR solutions are possible. See [264, pp. 201-204] for details. In practice, approximate ``pseudo QMF'' filters are more practical, which only give approximate perfect reconstruction. We'll return to this topic in §11.7.1.

The Haar filters, which we saw gave perfect reconstruction in the amplitude-complementary case, are also examples of a QMF filter bank:

\begin{eqnarray*}
H_0(z) &=& 1 + z^{-1}\\ [0.1in]
H_1(z) &=& 1 - z^{-1}
\end{eqnarray*}

In this example, $ c_0=c_1=1$, and $ n_0=n_1=0$.


Order a Hardcopy of Spectral Audio Signal Processing

Previous: Polyphase Decomposition of Haar Example
Next: Linear Phase Quadrature Mirror Filter Banks

written by Julius Orion Smith III
Julius Smith's background is in electrical engineering (BS Rice 1975, PhD Stanford 1983). He is presently Professor of Music and Associate Professor (by courtesy) of Electrical Engineering at Stanford's Center for Computer Research in Music and Acoustics (CCRMA), teaching courses and pursuing research related to signal processing applied to music and audio systems. See http://ccrma.stanford.edu/~jos/ for details.


Comments


No comments yet for this page


Add a Comment
You need to login before you can post a comment (best way to prevent spam). ( Not a member? )