Filtro biquad digital

Filtro lineal digital recursivo de segundo orden

En el procesamiento de señales , un filtro biquad digital es un filtro lineal recursivo de segundo orden , que contiene dos polos y dos ceros . "Biquad" es una abreviatura de " biquadratic ", que se refiere al hecho de que en el dominio Z , su función de transferencia es la relación de dos funciones cuadráticas :

${\displaystyle \ H(z)={\frac {b_{0}+b_{1}z^{-1}+b_{2}z^{-2}}{a_{0}+a_{1}z^{-1}+a_{2}z^{-2}}}}$

Los coeficientes suelen normalizarse de forma que 0 ₌ 1:

${\displaystyle \ H(z)={\frac {b_{0}+b_{1}z^{-1}+b_{2}z^{-2}}{1+a_{1}z^{-1}+a_{2}z^{-2}}}}$

Los filtros de respuesta al impulso infinito de orden superior pueden ser muy sensibles a la cuantificación de sus coeficientes y pueden volverse inestables fácilmente . Esto es un problema mucho menor con los filtros de primer y segundo orden; por lo tanto, los filtros de orden superior se implementan típicamente como secciones biquad en cascada en serie (y un filtro de primer orden si es necesario). Los dos polos del filtro biquad deben estar dentro del círculo unitario para que sea estable. En general, esto es cierto para todos los filtros discretos, es decir, todos los polos deben estar dentro del círculo unitario en el dominio Z para que el filtro sea estable.

Implementación

Forma directa 1

La implementación más sencilla es la forma directa 1, que tiene la siguiente ecuación diferencial :

${\displaystyle \ y[n]={\frac {1}{a_{0}}}\left(b_{0}x[n]+b_{1}x[n-1]+b_{2}x[n-2]-a_{1}y[n-1]-a_{2}y[n-2]\right)}$

o, si se normaliza:

${\displaystyle \ y[n]=b_{0}x[n]+b_{1}x[n-1]+b_{2}x[n-2]-a_{1}y[n-1]-a_{2}y[n-2]}$

Aquí los coeficientes , y determinan ceros, y , determinan la posición de los polos. ${\displaystyle b_{0}}$ ${\displaystyle b_{1}}$ ${\displaystyle b_{2}}$ ${\displaystyle a_{1}}$ ${\displaystyle a_{2}}$

Diagrama de flujo del filtro biquad en forma directa 1:

Cuando estas secciones se conectan en cascada para filtros de orden mayor que 2, se puede mejorar la eficiencia de la implementación al notar que el retraso de la salida de una sección se clona en la entrada de la siguiente sección. Se pueden eliminar dos componentes de retraso de almacenamiento entre secciones. ${\displaystyle z^{-1}}$

Forma directa 2

La forma directa 2 implementa la misma función de transferencia normalizada que la forma directa 1, pero en dos partes:

${\displaystyle \ y[n]=b_{0}w[n]+b_{1}w[n-1]+b_{2}w[n-2],}$

y usando la ecuación diferencial :

${\displaystyle \ w[n]=x[n]-a_{1}w[n-1]-a_{2}w[n-2].}$

Diagrama de flujo del filtro biquad en forma directa 2:

La implementación de la forma directa 2 solo necesita N unidades de retardo, donde N es el orden del filtro, potencialmente la mitad que la forma directa 1. La derivación de la forma directa 1 normalizada es la siguiente:

${\displaystyle \ y[n]=b_{0}x[n]+b_{1}x[n-1]+b_{2}x[n-2]-a_{1}y[n-1]-a_{2}y[n-2]}$

Supongamos la sustitución:

${\displaystyle \ y[n]=b_{0}w[n]+b_{1}w[n-1]+b_{2}w[n-2]}$

Lo que resulta en:

${\displaystyle \ y[n]=b_{0}x[n]+b_{1}x[n-1]+b_{2}x[n-2]-a_{1}(b_{0}w[n-1]+b_{1}w[n-2]+b_{2}w[n-3])-a_{2}(b_{0}w[n-2]+b_{1}w[n-3]+b_{2}w[n-4])}$

Aislando los coeficientes , y : ${\displaystyle b_{0}}$ ${\displaystyle b_{1}}$ ${\displaystyle b_{2}}$

${\displaystyle \ y[n]=b_{0}(x[n]-a_{1}w[n-1]-a_{2}w[n-2])+b_{1}(x[n-1]-a_{1}w[n-2]-a_{2}w[n-3])+b_{2}(x[n-2]-a_{1}w[n-3]-a_{2}w[n-4])}$

Suponiendo que se obtiene el resultado anterior: ${\displaystyle \ w[n]=x[n]-a_{1}w[n-1]-a_{2}w[n-2]}$

${\displaystyle \ y[n]=b_{0}w[n]+b_{1}w[n-1]+b_{2}w[n-2],}$

La desventaja es que la forma directa 2 aumenta la posibilidad de desbordamiento aritmético para filtros de Q alto o resonancia. ^[1] Se ha demostrado que a medida que Q aumenta, el ruido de redondeo de ambas topologías de forma directa aumenta sin límites. ^[2] Esto se debe a que, conceptualmente, la señal pasa primero a través de un filtro de todos los polos (que normalmente aumenta la ganancia en las frecuencias resonantes) antes de que el resultado de eso se sature, luego pasa a través de un filtro de todos los ceros (que a menudo atenúa mucho de lo que amplifica la mitad de todos los polos).

La implementación de la forma directa 2 se denomina forma canónica, porque utiliza la mínima cantidad de retrasos, sumadores y multiplicadores, produciendo la misma función de transferencia que la implementación de la forma directa 1.

Formas directas transpuestas

Cada una de las dos formas directas se puede transponer invirtiendo el gráfico de flujo sin alterar la función de transferencia. Los puntos de ramificación se cambian por sumadores y los sumadores se cambian por puntos de ramificación. ^[3] Estos proporcionan implementaciones modificadas que logran la misma función de transferencia que puede ser matemáticamente significativa en una implementación del mundo real donde la precisión se puede perder en el almacenamiento de estados.

Las ecuaciones diferenciales para la forma directa transpuesta 2 son:

${\displaystyle \ y[n]=b_{0}x[n]+s_{1}[n-1],}$

dónde

${\displaystyle \ s_{1}[n]=s_{2}[n-1]+b_{1}x[n]-a_{1}y[n]}$

y

${\displaystyle \ s_{2}[n]=b_{2}x[n]-a_{2}y[n].}$

Forma directa transpuesta 1

La forma directa 1se transpone a

Forma directa transpuesta 2

La forma directa 2se transpone a

Cuantificación del ruido

Cuando una muestra de n bits se multiplica por un coeficiente de m bits, el producto tiene n+m bits. Estos productos se acumulan típicamente en un registro DSP, la suma de cinco productos puede necesitar 3 bits de desbordamiento; este registro suele ser lo suficientemente grande como para contener n+m+3 bits. El z ⁻¹ se implementa almacenando un valor para un tiempo de muestra; este registro de almacenamiento suele tener n bits, el registro acumulador se redondea para que quepan n bits y esto introduce ruido de cuantificación.

En la disposición de forma directa 1, hay una única función de cuantificación/redondeo Q(z):

En la disposición de forma directa 2, también hay una función de cuantificación/redondeo para un valor intermedio. En una cascada, es posible que el valor no necesite redondeo entre etapas, pero el resultado final sí puede necesitarlo.

El DSP de punto fijo suele preferir las formas no transpuestas y tiene un acumulador con una gran cantidad de bits, y se redondea cuando se almacena en la memoria principal. El DSP de punto flotante suele preferir la forma transpuesta, cada multiplicación y potencialmente cada adición se redondean; las adiciones son un resultado de mayor precisión, cuando ambos operandos tienen una magnitud similar.

Véase también

• Filtro biquad
• Filtro digital

Referencias

1. JO Smith III, Forma directa II
2. LB Jackson, "Sobre la interacción del ruido de redondeo y el rango dinámico en filtros digitales", Bell Sys. Tech. J. , vol. 49 (febrero de 1970), reimpreso en Digital Signal Process , LR Rabiner y CM Rader, Eds. (IEEE Press, Nueva York, 1972).
3. "Formas directas transpuestas".

Enlaces externos

• Fórmulas de libro de cocina para los coeficientes del filtro biquad del ecualizador de audio
• WikiLibro sobre procesamiento de señales digitales
• Filtros digitales de segundo orden acoplados