¿Qué es una explicación intuitiva de los filtros FIR?

Los filtros generalmente se entienden formalmente en el dominio de la frecuencia. Sin embargo, dado que ha pedido una explicación intuitiva de lo mismo, intentaré hacerlo en el dominio del tiempo.

Cualquier señal, digamos audio, son datos digitales y está compuesta por un montón de muestras digitales. Cada “muestra” es en realidad solo un número (suponga un sonido de tono simple). Así que piense en un archivo de audio como un grupo de números almacenados en una matriz (aunque en realidad no lo es).

Primero, pensemos en la estructura de un filtro FIR.
Piense en el primer coeficiente como simplemente un factor de amplificación. Piense en el segundo como factor de amplificación para la muestra anterior, y el tercero como la amplificación de la segunda muestra anterior, y así sucesivamente. Entonces, todos esos coeficientes son como, tomar datos en serie, multiplicar un montón de ellos por algunos valores. Se llama filtro de línea de retardo girado.

Ahora, ¿cómo tiene sentido eso?
Imagine que desea un filtro de paso bajo. Básicamente, desea eliminar todos los componentes de alta frecuencia o los picos. Intuitivamente, eso es un poco como un promedio, ¿no? Promediar es algo así como sumar todos los componentes y dividirlos por el número de valores, ¿verdad? En los filtros FIR, usted hace algo como el promedio ponderado de todas las muestras de entrada. Entonces, eso significaría que todos los h (n) son números positivos, ya que los pesos no pueden ser negativos. Cuanto más similares son los h (n), mejor y mejor es el “promedio”, y viceversa. El promedio ponderado se reduciría a un promedio simple cuando todos los coeficientes sean iguales. Por “mejor promedio”, me refiero a una frecuencia de corte más baja del filtro de paso bajo.

Si desea hacer un filtro de paso alto o un filtro de paso de banda, notará que algunos coeficientes comienzan a ser negativos. Al hacer que algunos de ellos sean negativos, esencialmente ha comenzado a tomar diferencias. Si la diferencia es alta, el valor es alto. Si son más similares, la diferencia es cero. Así que esencialmente está reteniendo los valores “diferentes”, o rechazando las frecuencias bajas y amplificando las frecuencias altas. En el caso extremo, cada coeficiente alternativo resulta negativo, lo que se convierte en un filtro de paso alto.

¿Por qué varía el número de coeficientes?
Considere el ejemplo promedio. Cuantos más valores use para tomar el promedio, mejor será la estimación, ¿verdad? Si toma solo 2 valores anteriores y observa el promedio, y observa 20 valores anteriores y observa el promedio, el que tiene 20 obviamente sería una mejor estimación del promedio “verdadero” de toda la señal. Entonces, como regla general, cuanto mayor sea el número de coeficientes, mejor será su filtro. Sin embargo, también tomaría una mayor cantidad de tiempo calcular los “promedios”, por lo que existe una compensación.

PD: Este fue un intento de dar una explicación simple e intuitiva de los filtros FIR en el dominio del tiempo. Toma cada detalle con una pizca de sal. He simplificado demasiado muchas cosas para hacerlo más intuitivo.

Un filtro FIR no tiene equivalente en el dominio analógico. Los filtros analógicos tienen características similares a los filtros IIR.

Un filtro FIR M-tap comprende los siguientes elementos:

  • una línea de retardo intermitente muestreada a la frecuencia de muestreo [math] f_s [/ math]
  • Multiplicadores de coeficiente M
  • sumadores

La línea de retardo girado retrasa la entrada x (n) exactamente un período de muestreo [matemática] T = \ dfrac {1} {f_s} [/ matemática] por cada pulsación.

La señal x (n) y sus versiones retrasadas [matemática] x (n-1) [/ matemática], [matemática] x (n-2) [/ matemática]… [matemática] x (nM) [/ matemática] son multiplicado por los coeficientes [matemática] b_0 [/ matemática], [matemática] b_1 [/ matemática]… [matemática] b_M [/ matemática] y combinados junto con los sumadores. La combinación puede ser constructiva o destructiva.

Los coeficientes se calculan de tal manera que la combinación es constructiva en la banda de paso y destructiva en la banda de parada. De manera similar a la interferencia en las ondas, cuando la señal y las versiones retrasadas se suman en fase, la salida alcanza la amplitud máxima para esa frecuencia. Por el contrario, cuando las señales y las versiones retrasadas están desfasadas, la amplitud de salida es mínima para esa frecuencia.

Consideremos algunos filtros FIR muy simples para comprender cómo las salidas del retraso derivado pueden combinarse de manera constructiva o destructiva. Los filtros son filtros de paso bajo con M = 4, 8 y 16 coeficientes. Los coeficientes de cada filtro son los mismos.

[matemáticas] b_0 = b_1 =… = b_M = \ dfrac {1} {M} [/ matemáticas]

La gráfica representa la respuesta de amplitud de los filtros.

[matemática] {\ color {azul} {\ text {M = 4 ___}}} [/ matemática] [matemática] {\ color {verde} {\ text {M = 8 ___}}} [/ matemática] [matemática] { \ color {rojo} {\ text {M = 16 ___}}} [/ matemáticas]

El eje horizontal representa la frecuencia normalizada [matemática] \ dfrac {f} {f_s} [/ matemática], donde [matemática] f_s [/ matemática] es la frecuencia de muestreo.

Supongamos que [math] f_s = 8 [/ math] [math] kHz [/ math].

Vamos a probar los filtros con una entrada constante y una entrada de onda sinusoidal con unidad de amplitud y frecuencia [math] f_0 [/ math].

Para la entrada de onda sinusoidal, la salida de los grifos viene dada por

[matemáticas] x (n) = sen 2 \ pi \ dfrac {f_0} {f_s} n [/ matemáticas]
[matemáticas] x (n-1) = sen 2 \ pi \ dfrac {f_0} {f_s} (n-1) [/ matemáticas]

[matemáticas] x (nk) = sen 2 \ pi \ dfrac {f_0} {f_s} (nk) [/ matemáticas]

[matemáticas] x (nM) = sen 2 \ pi \ dfrac {f_0} {f_s} (nM) [/ matemáticas]

FIR con M = 4

El resultado es:

[matemáticas] y (n) = \ dfrac {1} {4} [x (n) + x (n-1) + x (n-2) + x (n-3)] [/ matemáticas]

Para una entrada constante = 1

[matemáticas] x (n) = x (n-1) = x (n-2) = x (n-3) = 1 [/ matemáticas]

[matemáticas] y (n) = 1 [/ matemáticas]

Este es un ejemplo muy simple de combinación constructiva. Como la entrada es constante, la línea de retardo girado está en fase y la salida es máxima.

Consideremos ahora una entrada de onda sinusoidal con unidad de amplitud y frecuencia [matemática] f_0 = 2 [/ matemática] [matemática] kHz, [/ matemática] correspondiente al primer cero de la función de transferencia.

[matemáticas] x (n) = sin \ dfrac {\ pi} {2} n [/ matemáticas]

[matemáticas] x (n-1) = sin \ dfrac {\ pi} {2} (n-1) [/ matemáticas]

[matemáticas] x (n-2) = sin \ dfrac {\ pi} {2} (n-2) = sin (\ dfrac {\ pi} {2} n – \ pi) = – x (n) [/ matemáticas]

[matemáticas] x (n-3) = sin \ dfrac {\ pi} {2} (n-3) = sin (\ dfrac {\ pi} {2} (n-1) – \ pi) = – x ( n-1) [/ matemáticas]

Como podemos ver [math] x (n) [/ math] y [math] x (n-2) [/ math] están 180 [math] ^ \ circ [/ math] fuera de fase. Lo mismo es cierto para [matemáticas] x (n-1) [/ matemáticas] y [matemáticas] x (n-3) [/ matemáticas]. Como consecuencia de [matemática] f_0 = 2 [/ matemática] [matemática] kHz [/ matemática]

[matemáticas] y (n) = 0 [/ matemáticas]

Este es un ejemplo de combinación totalmente destructiva de las salidas de línea de retardo girado.

FIR con M = 8

El resultado es:

[matemáticas] y (n) = \ dfrac {1} {8} [x (n) + x (n-1) + x (n-2) + x (n-3) + x (n-4) + x (n-5) + x (n-6) + x (n-7)] [/ matemáticas]

Es fácil ver que para una entrada constante = 1

[matemáticas] y (n) = 1 [/ matemáticas]

Consideremos ahora una entrada de onda sinusoidal con unidad de amplitud y frecuencia [matemática] f_0 = 1 [/ matemática] [matemática] kHz, [/ matemática] correspondiente al primer cero de la función de transferencia.

[matemáticas] x (n) = sin \ dfrac {\ pi} {4} n [/ matemáticas]

[matemáticas] x (n-1) = sin \ dfrac {\ pi} {4} (n-1) [/ matemáticas]

[matemáticas] x (n-2) = sin \ dfrac {\ pi} {4} (n-2) [/ matemáticas]

[matemáticas] x (n-3) = sin \ dfrac {\ pi} {4} (n-3) [/ matemáticas]

[matemáticas] x (n-4) = sin \ dfrac {\ pi} {4} (n-4) = sin (\ dfrac {\ pi} {4} n – \ pi) = – x (n) [/ matemáticas]

[matemáticas] x (n-5) = sin \ dfrac {\ pi} {4} (n-5) = sin (\ dfrac {\ pi} {4} (n-1) – \ pi) = – x ( n-1) [/ matemáticas]

[matemáticas] x (n-6) = sin \ dfrac {\ pi} {4} (n-6) = sin (\ dfrac {\ pi} {4} (n-2) – \ pi) = – x ( n-2) [/ matemáticas]

[matemáticas] x (n-7) = sin \ dfrac {\ pi} {4} (n-7) = sin (\ dfrac {\ pi} {4} (n-3) – \ pi) = – x ( n-3) [/ matemáticas]

[matemática] x (n) [/ matemática], [matemática] x (n-1) [/ matemática], [matemática] x (n-2) [/ matemática] y [matemática] x (n-3) [ / math] están respectivamente 180 [math] ^ \ circ [/ math] fuera de fase con respecto a [math] x (n-4) [/ math], [math] x (n-5) [/ math], [matemáticas] x (n-6) [/ matemáticas] y [matemáticas] x (n-7) [/ matemáticas]. Como consecuencia de [matemáticas] f_0 = 1 kHz [/ matemáticas]

[matemáticas] y (n) = 0 [/ matemáticas]

Este es nuevamente un ejemplo de combinación totalmente destructiva de las salidas de línea de retardo intervenido.

Los ejemplos anteriores son ejemplos de combinaciones totalmente constructivas y totalmente destructivas de las salidas de línea de retardo girado.

Es intuitivo que con solo cambiar la frecuencia [matemática] f_0 [/ matemática] es posible encontrar toda una serie de combinaciones parcialmente constructivas o destructivas que producen todos los valores de la respuesta de amplitud.

Los filtros utilizados en los ejemplos son a propósito muy simples, por lo que es más fácil mostrar cómo se combinan las salidas de la línea de retardo girado. Al elegir adecuadamente los coeficientes, es posible controlar mejor la forma en que se combinan las salidas de la línea de retardo girado. Esto hace posible diseñar filtros FIR con banda de paso muy plana, bandas de transición muy estrechas y rechazos más altos en la banda de parada.

La implementación de FIR básicamente está tratando de realizar una convolución de las muestras de entrada con las muestras almacenadas de la función sinc (derivaciones de filtro).

Uno puede recordar que sinc tiene forma rectangular en el dominio de frecuencia. Además, la convolución en el dominio del tiempo puede considerarse como una multiplicación en el dominio de la frecuencia.

Entonces, mediante la operación FIR, estamos multiplicando el espectro de entrada con una ventana rectangular, realizando la operación de filtrado.