Filtro QMF

Un filtro QMF (Quadrature Mirror Filter in inglese) consente di suddividere un segnale in due segnali sottocampionati che possono essere poi ricostruiti senza aliasing, anche qualora vengano utilizzati filtri non ideali[1]. Tale filtro fu introdotto da Croisier, Esteban e Galand nel 1976[2] .

Definizione

Consideriamo il seguente schema a blocchi:


dove i due filtri H 0 {\displaystyle H_{0}} e H 1 {\displaystyle H_{1}} (rispettivamente un filtro passa basso ed un filtro passa alto) suddividono il segnale, mentre G 0 {\displaystyle G_{0}} e G 1 {\displaystyle G_{1}} operano la ricostruzione[3].
La condizione necessaria e sufficiente per la cancellazione dell'aliasing è:

G 0 ( z ) H 0 ( − z ) + G 1 ( z ) H 1 ( − z ) = 0. {\displaystyle G_{0}(z)H_{0}(-z)+G_{1}(z)H_{1}(-z)=0.}

Una soluzione a tale problema è data dalla configurazione di filtri in quadratura a specchio (QMF), che si ottiene imponendo i seguenti vincoli sui filtri d'analisi[1]:

H 1 ( z ) = H 0 ( − z ) . {\displaystyle H_{1}(z)=H_{0}(-z).}
G 0 ( z ) = H 1 ( − z ) = H 0 ( z ) . {\displaystyle G_{0}(z)=H_{1}(-z)=H_{0}(z).}
G 1 ( z ) = − H 1 ( z ) = − H 0 ( − z ) . {\displaystyle G_{1}(z)=-H_{1}(z)=-H_{0}(-z).}

Il primo vincolo stabilisce una simmetria tra H 0 {\displaystyle H_{0}} e H 1 {\displaystyle H_{1}} , determinando la denominazione di filtri a specchio. Sostituendo i vincoli soprastanti nella precedente condizione si ottiene:

H 0 ( z ) H 0 ( − z ) − H 0 ( − z ) H 0 ( z ) = 0. {\displaystyle H_{0}(z)H_{0}(-z)-H_{0}(-z)H_{0}(z)=0.}

Ciò conferma che l'aliasing è annullato per filtri QMF.
Per poter ottenere una ricostruzione perfetta, la seguente condizione dev'essere soddisfatta:

G 0 ( z ) H 0 ( z ) + G 1 ( z ) H 1 ( z ) = 2 z − l . {\displaystyle G_{0}(z)H_{0}(z)+G_{1}(z)H_{1}(z)=2z^{-l}.}

Visti i vincoli imposti, per i filtri QMF la condizione è:

H 0 2 ( z ) − H 0 2 ( − z ) = 2 z − l . {\displaystyle {H_{0}}^{2}(z)-{H_{0}}^{2}(-z)=2z^{-l}.}

A questo punto siamo in grado di ricavare tutti e quattro i filtri tramite H 0 ( z ) . {\displaystyle H_{0}(z).}
Se ci limitiamo a considerare filtri FIR (Finite Impulse Response) la precedente equazione può esser soddisfatta in modo esatto esclusivamente da filtri della forma:

H 0 ( z ) = c 0 z − 2 n 0 + c 1 z − ( 2 n 1 + 1 ) {\displaystyle H_{0}(z)=c_{0}z^{-2n_{0}}+c_{1}z^{-(2n_{1}+1)}}
H 1 ( z ) = c 0 z − 2 n 0 − c 1 z − ( 2 n 1 + 1 ) {\displaystyle H_{1}(z)=c_{0}z^{-2n_{0}}-c_{1}z^{-(2n_{1}+1)}}

dove c 0 {\displaystyle c_{0}} e c 1 {\displaystyle c_{1}} son costanti e n 0 {\displaystyle n_{0}} e n 1 {\displaystyle n_{1}} sono interi.
Un esempio di filtri QMF sono i filtri Haar della forma:

H 0 ( z ) = 1 + z − 1 . {\displaystyle H_{0}(z)=1+z^{-1}.}
H 1 ( z ) = 1 − z − 1 . {\displaystyle H_{1}(z)=1-z^{-1}.}

In tal caso c 0 = c 1 = 1 {\displaystyle c_{0}=c_{1}=1} e n 0 = n 1 = 0. {\displaystyle n_{0}=n_{1}=0.}
Maggiori informazioni sulla progettazione di filtri QMF sono reperibili in Simoncelli et al.[4].

Note

  1. ^ a b Martin Vetterli, Jelena Kovačević, Wavelets and Subband Coding. Prentice Hall PTR, Englewood Cliffs, New Jersey, 1995. Re-issued by the authors in 2007.
  2. ^ A. Croisier, D. Esteban e C. Galand, Perfect channel splitting by use of interpolation/decimation/tree decomposition techniques. International Conference on Information Sciences and Systems, pagg. 443–446, Patrasso, Grecia, Agosto 1976.
  3. ^ F. Rocca, Elaborazione numerica dei segnali Archiviato il 21 ottobre 2012 in Internet Archive.. Politecnico di Milano, 2010.
  4. ^ Simoncelli et al., Subband Transform - 4.4: Quadrature mirror filter Archiviato il 13 dicembre 2005 in Internet Archive.. Editor: John Woods. Kluwer Academic Press, 1990.

Voci correlate

  • Wavelet
  • ATRAC
  • Trasformata di Fourier

Collegamenti esterni

  • Julius O. Smith III - Spectral Audio Signal Processing - QMF. Data ultimo accesso: 28/04/2012