Bitácora · DCSS · Unidad 2

Las imágenes no son un eco débil: son la copia exacta

Unidad 2 · Conversores

La nota anterior mostró que dos tonos que aliasan producen las mismas muestras: no hay eco, hay una sola lista de números. Acá la misma idea aparece en otro punto del sistema. Diezmar sin filtrar mete el alias con una amplitud idéntica a la de la banda base, no atenuada; insertar ceros para interpolar crea réplicas de la misma magnitud que la señal original. En los dos casos, lo que vuelve no es una sombra débil del original: es la copia exacta.

Un tono en ω0=0,8π diezmado por M=3 aparece en 0,4000π con una amplitud de 1,000000 — medido con un largo de señal múltiplo exacto del período del tono plegado, así que no hay fuga espectral que ensucie la lectura. No es una coincidencia de la medición: cos(0,8π·M·n)=cos(2,4π·n)=cos(0,4π·n) es el mismo coseno, escrito distinto. Lo que sigue es por qué eso es inevitable, y qué lo arregla.

01 · Diezmar estira el espectro y lo repite

Diezmar por M es quedarse con una de cada M muestras: xd[n]=x[nM]. El nuevo período es M veces más largo, y como ω=Ω·T (la nota anterior lo usa igual), un período más largo estira el eje ω por M: lo que ocupaba hasta ωN ahora ocupa hasta M·ωN. Y aparece algo más: el espectro diezmado no es una sola copia estirada, son M copias corridas cada 2π/M y promediadas:

Xd(ejω)=(1/M)·∑iX(ej(ω−2πi)/M)

— con i corriendo de 0 a M−1. Si esas M copias no se tocan, la del medio sale intacta. Si se tocan —hay contenido por encima de π/M— sus colas se suman, exactamente como en la nota anterior, y una suma no se deshace mirándola. El apunte de la cátedra lo resume en dos casos: una señal sobremuestreada al doble, diezmada por 2, cae justo en el borde crítico sin perder nada; la misma señal ya en el borde crítico, diezmada por 2, se pliega y pierde la mitad de arriba. La diferencia entre los dos casos es si π/M alcanza para lo que hay adentro, y `diezmar` por sí sola no lo sabe ni lo pregunta: sólo tira muestras.

Simulación 1 Diezmar por M = 3, con y sin filtro anti-alias — arriba, el tono antes de diezmar (y la curva del filtro, si lo activás); abajo, dónde cae después de diezmar, con la altura de referencia de la banda base.
ω0 0,1333π cae en, tras diezmar 0,4000π amplitud medida 1,000000

Con el filtro apagado, el número de la derecha no se mueve de 1,000000 en ningún lado del deslizador: `diezmar` no sabe qué hay arriba de π/M, así que no lo toca, y lo que había arriba baja intacto a ocupar el lugar de lo que había abajo. El «decimador completo» que pide el apunte —pasabajos con corte π/M y ganancia 1, antes de tirar muestras— es lo que activa la casilla: mata lo que está por encima del corte con un filtro real, no con una pared ideal, y por eso la atenuación crece de forma gradual cruzando π/M en vez de cortarse de un salto. Preferís perder esa banda, filtrada limpia, a que se pliegue encima de la que te importa.

02 · Insertar ceros comprime el espectro y lo repite

Insertar L−1 ceros entre cada muestra es la operación espejo: xe[n]=x[n/L] en los múltiplos de L, cero en todo lo demás. La DTFT de eso se escribe en tres líneas —es una sustitución directa, ningún truco de por medio— y da una identidad exacta:

Xe(ejω)=X(ejωL)

El eje se comprime por L: lo que ocupaba [0,π] ahora ocupa [0,π/L]. Pero la DTFT es periódica en 2π, así que adentro de [0,2π) entran L copias enteras del espectro original: las imágenes. No son aliasing —nada se mezcló, no se perdió información, la identidad de arriba es exacta— son copias espurias que hay que filtrar. Y son copias de verdad: la magnitud de una imagen en ωL es la misma magnitud que la banda base en ω, no una fracción de ella. Verificada contra la respuesta en frecuencia de referencia sobre una grilla de frecuencias, la identidad vale al piso de la precisión de la máquina: la diferencia máxima medida es 1,65×10−14. El «interpolador completo» —insertar ceros, después un pasabajos con corte π/L y ganancia L, no 1— necesita esa ganancia justamente porque insertar ceros diluye la amplitud media por L: la energía útil quedó en 1 de cada L muestras, y el filtro se la devuelve al mismo tiempo que borra las imágenes.

Para tener a mano

  1. Diezmar por M estira el eje ω por M y suma M copias corridas cada 2π/M. Si se pasan de π, se pliegan y sus colas se suman: es el mismo «irreversible» de la nota anterior, en otro lugar del sistema.
  2. Diezmar sin filtro no atenúa nada: es sólo un submuestreo de la secuencia, y el alias entra con la misma amplitud que la banda base. Medido: 1,000000 para un tono en 0,8π diezmado por 3.
  3. El decimador completo mata primero, con un pasabajos de corte π/M, lo que no va a sobrevivir al plegado. Lo que se pierde se pierde filtrado, no plegado.
  4. Insertar ceros comprime el eje ω por L: Xe(ejω)=X(ejωL) es una identidad exacta, no una aproximación, y las L imágenes que aparecen tienen la misma magnitud que la banda base.
  5. El interpolador completo necesita, además del pasabajos con corte π/L, una ganancia L: sin ella la energía repartida entre los ceros insertados queda diluida.
← Anterior Siguiente → Inicio