Unidad 1 · Filtrado digital
El apunte de la cátedra lo dice así: truncar el sinc ideal con una ventana rectangular produce rizado de Gibbs, y con trece o con veinticinco coeficientes el rizado no baja — sólo se concentra. Medido acá, con las mismas funciones que dibujan el widget: con taps el sobrepico es 9,67%; con — ocho veces más taps — es 9,02%. Bajó 0,65 puntos, no 87,5%. Gibbs no se cura con más cómputo: converge a un número y se queda ahí.
La intuición que cualquiera trae de otros problemas numéricos es «más resolución, menos error», y por eso el apunte se toma el trabajo de marcarla como falsa. Acá no hay más resolución que agregar: el problema no es que falten muestras, es la forma de cortar.
El pasabajos ideal de corte tiene una respuesta al impulso que es un sinc, infinita hacia los dos lados y no causal: , que en vale por L'Hôpital, sin ningún salto. Para usarlo hay que cortarlo a muestras y correrlo para que empiece en . La manera más obvia de cortar es multiplicar por una ventana que vale en las muestras que se quedan y en el resto: la ventana rectangular, que no es «no usar ventana» — es la ventana más brusca de todas.
Multiplicar dos señales en el tiempo es convolucionar sus transformadas en frecuencia. La respuesta que sale, , no es la plantilla ideal: es la plantilla ideal emborronada con la transformada de la ventana. Esa transformada tiene un lóbulo principal —que ensancha la transición— y una fila de lóbulos laterales que no se apagan: son ellos los que dibujan el rizado a los dos lados de cada corte, en la banda de paso y en la de rechazo.
Antes de mover nada: si subís los taps ocho veces —de 21 a 161— ¿el sobrepico baja también a un octavo?
Subir los taps casi no baja la altura: 9,67% a 9,02% entre los dos extremos del deslizador, medido sobre los mismos filtros que dibuja el widget. Lo que sí se mueve, y mucho, es el lugar del pico: se acerca al corte cada vez que agregás taps, y la distancia se achica a la mitad cada vez que se duplica — de 0,10 a 0,05 a 0,025 a 0,0125, en unidades de , entre y . El rizado no se cura: se comprime contra el borde. Es exactamente lo que la filmina dibuja con 13 y con 25 coeficientes, sólo que acá con números.
Si más taps no bajan el rizado, ¿qué lo baja? Cambiar la ventana. La rectangular corta de golpe: salta de a en el borde, y esa discontinuidad es justo lo que le pega a su transformada los lóbulos altos. Hamming y Blackman no cortan de golpe: bajan gradualmente hasta casi cero en los dos extremos, así que su propia transformada tiene lóbulos mucho más chicos — a costa de un lóbulo principal más ancho, o sea una transición más lenta para el mismo . No hay almuerzo gratis, y el apunte lo dice tal cual: más atenuación de rizado implica transición más lenta.
| Ventana | Peor lóbulo lateral | «Primer pico» (la trampa) |
|---|---|---|
| Rectangular | −13,26 dB | −13,26 dB — coinciden |
| Hamming | −42,53 dB | −46,21 dB |
| Blackman | −58,11 dB | −96,21 dB |
En la rectangular el primer lóbulo ya es el peor: sus lóbulos decrecen todo el camino, así que medir «el primer pico» o medir «el máximo sobre todo el rango fuera de banda» dan el mismo número. Con Hamming y con Blackman no: hay un mínimo local justo después del lóbulo principal, y después los lóbulos se estabilizan en un piso más alto que ese primer valle. Medir mal —tomar el primer pico como si fuera el peor lóbulo— parece razonable y da un número que miente.
En Hamming la mentira es chica: −46,21 dB contra los −42,53 dB reales, 3,68 dB de margen que no existe. En Blackman es grande: −96,21 dB contra −58,11 dB reales, 38 dB de error — el filtro parece ochenta veces mejor de lo que es. (Ochenta y no mil: un lóbulo lateral es una razón de amplitud, así que 38 dB son , no .) Es la misma cuenta, con la misma calculadora; lo único que cambia es dónde se para a mirar. Por eso, para medir el lóbulo lateral, hay que barrer todo el rango fuera del lóbulo principal y quedarse con el máximo — nunca con el primero que aparece, que casi nunca es el más alto.
El sobrepico de Gibbs no es un problema de resolución, y por eso ningún lo arregla: es la forma de la ventana. La rectangular corta de golpe y paga −13 dB de por vida, sea o . Cambiar a Hamming o a Blackman no es agregarle más cómputo al mismo filtro: es diseñar un filtro distinto, con otro canje entre rizado y ancho de transición. Esa es la pregunta que hay que hacerse antes de mover el deslizador de taps, no después.