Editing Quadrature mirror filter

In [[digital signal processing]], a '''quadrature mirror filter''' is a filter whose magnitude response is the mirror image around <math>\pi/2</math> of that of another filter. Together these filters, first introduced by Croisier et al., are known as the quadrature mirror filter pair.

A filter <math>H_1(z)</math> is the quadrature mirror filter of <math>H_0(z)</math> if <math>H_1(z) = H_0(-z)</math>.

The filter responses are symmetric about <math>\Omega = \pi / 2</math>:

: <math>\big|H_1\big(e^{j\Omega}\big)\big| = \big|H_0\big(e^{j(\pi - \Omega)}\big)\big|.</math>

In audio/voice codecs, a quadrature mirror filter pair is often used to implement a [[filter bank]] that splits an input [[signal processing|signal]] into two bands. The resulting high-pass and low-pass signals are often reduced by a factor of 2, giving a critically sampled two-channel representation of the original signal. The analysis filters are often related by the following formula in addition to quadrate mirror property:
: <math>\big|H_0\big(e^{j\Omega}\big)\big|^2 + \big|H_1\big(e^{j\Omega}\big)\big|^2 = 1,</math>
where <math>\Omega</math> is the [[frequency]], and the sampling rate is normalized to <math>2\pi</math>.
This is known as power complementary property.
In other words, the power sum of the high-pass and low-pass filters is equal to 1.  

Orthogonal [[wavelet]]s – the [[Haar wavelet]]s and related [[Daubechies wavelet]]s, [[Coiflet]]s, and some developed by [[Stéphane Mallat|Mallat]], are generated by [[Wavelet#Scaling_function|scaling functions]] which, with the wavelet, satisfy a quadrature mirror filter relationship.

==Relationship to other filter banks==
The earliest wavelets were based on expanding a function in terms of rectangular steps, the Haar wavelets. This is usually a poor approximation, whereas Daubechies wavelets are among the simplest but most important families of wavelets. A linear filter that is zero for “smooth” signals, given a record of <math>N</math> points <math>x_n</math> is defined as
: <math>y_n = \sum_{i=0}^{M-1} b_i x_{n-i}.</math>

It is desirable to have it vanish for a constant, so taking the order <math>m = 4</math>, for example,
: <math>b_0 \cdot 1 + b_1 \cdot 1 + b_2 \cdot 1 + b_3 \cdot 1 = 0.</math>
And to have it vanish for a linear ramp, so that
: <math>b_0 \cdot 0 + b_1 \cdot 1 + b_2 \cdot 2 + b_3 \cdot 3 = 0.</math>

A linear filter will vanish for any <math>x = \alpha n + \beta</math>, and this is all that can be done with a fourth-order wavelet. Six terms will be needed to vanish a quadratic curve, and so on, given the other constraints to be included. Next an accompanying filter may be defined as
: <math>z_n = \sum_{i=0}^{M-1} c_i x_{n-i}.</math>

This filter responds in an exactly opposite manner, being large for smooth signals and small for non-smooth signals. A linear filter is just a convolution of the signal with the filter’s coefficients, so the series of the coefficients is the signal that the filter responds to maximally. Thus, the output of the second filter vanishes when the coefficients of the first one are input into it. The aim is to have
: <math>\sum_{i=0}^{M-1} c_i b_i = 0.</math>

Where the associated time series flips the order of the coefficients because the linear filter is a convolution, and so both have the same index in this sum. A pair of filters with this property are defined as quadrature mirror filters.<ref>
{{Citation
| last1=Gershenfeld
| first1=Neil 
| title=The Nature of Mathematical Modeling
| pages=132–135
| publisher=Cambridge University Press
| location=Cambridge, England
| isbn=0521570956
| year=1998}}.</ref>
Even if the two resulting bands have been subsampled by a factor of 2, the relationship between the filters means that approximately [[filter bank#Perfect reconstruction filter banks | perfect reconstruction]] is possible. That is, the two bands can then be upsampled, filtered again with the same filters and added together, to reproduce the original signal exactly (but with a small delay). (In practical implementations, numeric precision issues in [[floating-point arithmetic]] may affect the perfection of the reconstruction.)

==Further reading==
* A. Croisier, D. Esteban, C. Galand. ''Perfect channel splitting by use of interpolation/decimation tree decomposition techniques''. First International Conference on Sciences and Systems, Patras, August 1976, pp. 443–446.
* Johnston, J. D. [http://www.info490b.ece.mcgill.ca/Data/Exp4/Johnston.pdf ''A Filter Family Designed for use in Quadrature Mirror Filter Banks''].{{dead link|date=July 2016 |bot=InternetArchiveBot |fix-attempted=yes }}, Acoustics, Speech and Signal Processing, IEEE International Conference, 5, 291–294, April, 1980.
* [[Binomial QMF]], also known as [[Daubechies wavelet]] filters.
* [http://web.njit.edu/~akansu/sym.htm NJIT Symposia on Subbands and Wavelets 1990, 1992, 1994, 1997].
* Mohlenkamp, M. J. [http://www.ohiouniversityfaculty.com/mohlenka/20044/PASIII/waveletIPAM.pdf ''A Tutorial on Wavelets and Their Applications'']. University of Colorado, Boulder, Dept. of Applied Mathematics, 2004.
* Polikar, R. [https://web.archive.org/web/20180430094742/http://users.rowan.edu/~polikar/WAVELETS/WTpart4.html ''Multiresolution Analysis: The Discrete Wavelet Transform'']. Rowan University, NJ, Dept. of Electrical and Computer Engineering.

==References==
{{reflist}}

[[Category:Digital signal processing]]
[[Category:Filter theory]]
[[Category:Wavelets]]