La PizarraMatemática.

Matemáticas para ingeniería

Respuesta impulsional y operaciones de filtrado

Filtrar significa combinar datos con una regla: promediar para suavizar, restar para ver cambios o mantener una memoria que decae. Veremos exactamente qué muestra recibe cada peso. Después compararemos las salidas y sus frecuencias, cuidando el comienzo del registro y la memoria previa para no atribuir al filtro un borde que nosotros hemos añadido.

1. La respuesta al impulso describe el filtro LTI de estado cero

El impulso discreto δ[n] es una muestra uno en el origen y cero fuera. Para un filtro lineal e invariante en el tiempo, h[n] es la salida a ese impulso. Una entrada finita se descompone en impulsos desplazados y escalados; sumar sus respuestas produce la convolución.

La descripción es de estado cero. Si una implementación recursiva guarda memoria previa, añade una respuesta libre. Con entradas o respuestas infinitas se requieren hipótesis de convergencia; una tabla de coeficientes no elimina el dato inicial ni permite intercambiar sumas sin control.

Obtener h a partir de una regla de tres pesos

y[n]=x[n]/4+x[n−1]/2+x[n−2]/4.

  1. Impulso en cero

    En n=0 solo x[0] vale uno, así h[0]=1/4.
  2. Después

    En n=1 contribuye el peso del dato anterior: h[1]=1/2.
    En n=2 queda h[2]=1/4.
  3. Resto

    Fuera de esos índices la salida al impulso es cero.

    Así h=[1/4;1/2;1/4] desde cero.
Resultado

Cada posición de h conserva el retraso del peso correspondiente.

Lo comprobamos

La suma de pesos es uno, una propiedad que interpretaremos como conservación de constantes.

La memoria anterior añade una salida aunque no haya entrada

y[n]=(1/2)y[n−1]+(1/2)x[n], para n≥0.

  1. Reposo

    Con y[−1]=0 y entrada cero, todas las salidas son cero.

  2. Otra memoria

    Con y[−1]=4 y la misma entrada cero, y[0]=2, y[1]=1, y[2]=1/2,…
  3. Parte libre

    Esa contribución es 4(1/2)ⁿ⁺¹. No procede de convolucionar una entrada cero con h, sino del estado inicial.

Resultado

La misma recurrencia y la misma entrada no seleccionan una salida sin dato de memoria.

Lo comprobamos

Usar h de estado cero no autoriza a olvidar una y[−1] no nula.

2. Cada peso consulta un índice concreto antes de multiplicar

En y[n]=Σₖh[k]x[n−k], fijamos n y recorremos k. h[0] pesa el dato actual; h[1] el anterior; h[2] el de dos muestras antes. Primero se busca x[n−k], después se multiplica y al final se suman las aportaciones.

Índice negativo de h y peso negativo son cosas distintas. Un h[−1] no nulo consulta x[n+1], una muestra futura. Un h[1]=−1 sigue consultando el pasado y solo cambia el signo de su aportación. La causalidad se decide por posiciones, no por el signo de las amplitudes.

Usar tres datos con sus pesos correctos

h=[1/4;1/2;1/4], con x[n]=8, x[n−1]=4 y x[n−2]=0.

  1. Actual

    (1/4)·8=2.
  2. Anteriores

    (1/2)·4=2 y (1/4)·0=0.
  3. Sumar

    y[n]=2+2+0=4.

    Intercambiar ocho y cuatro cambiaría los datos asignados a pesos diferentes.

Resultado

La salida vale cuatro para esa ventana concreta.

Lo comprobamos

Una ventana constante de cuatro produciría también cuatro, pero por aportaciones distintas.

Un signo menos no implica pedir datos futuros

Comparar h₁[0]=1, h₁[1]=−1 con h₂[−1]=h₂[0]=1/2.

  1. Primera regla

    y₁[n]=x[n]−x[n−1] usa presente y pasado: es causal.
  2. Segunda regla

    y₂[n]=(x[n+1]+x[n])/2 pide el futuro: no es causal.
  3. Distinguir

    h₁ contiene un peso negativo, mientras h₂ contiene un índice negativo. Son características independientes.

Resultado

El argumento n−k permite comprobar causalidad sin adivinar.

Lo comprobamos

Un promedio puede ser no causal si la ventana se coloca hacia el futuro.

3. FIR e IIR clasifican la duración de h, no solo la forma de una implementación

FIR significa respuesta impulsional con soporte finito; IIR, una cola que no termina exactamente. Un FIR causal puede usar una ventana finita de entradas anteriores. Un IIR puede describirse con pocos coeficientes recursivos y ser estable si su cola es absolutamente sumable.

El criterio BIBO es Σ|h[k]|<∞ para el filtro de convolución discreto habitual. La recursividad de un programa no decide por sí sola FIR o IIR: una recurrencia puede representar un FIR bajo sus condiciones iniciales. La respuesta h efectiva y los estados internos se revisan por separado.

Una ventana finita deja de responder al impulso

h=[1/4;1/2;1/4] desde cero.

  1. Duración

    Solo hay tres posiciones no nulas: cero, uno y dos.

  2. Estabilidad

    La suma absoluta es uno, así |y[n]|≤M para entradas |x|≤M.

  3. Bloque finito

    Si la entrada ocupa N posiciones desde cero, la convolución completa tiene N+2 posiciones posibles. Algunas pueden ser cero por los datos.

Resultado

Es un FIR estable con memoria de dos muestras anteriores.

Lo comprobamos

No se recorta la cola de la convolución solo porque el registro original termine antes.

Una memoria que decae puede continuar para siempre

y[n]=(1/2)y[n−1]+(1/2)x[n], con reposo antes del origen.

  1. Impulso

    h[0]=1/2, h[1]=1/4, h[2]=1/8 y en general h[k]=(1/2)ᵏ⁺¹ para k≥0.
  2. Cola y suma

    No termina en un índice finito, pero Σ|h|=(1/2)/(1−1/2)=1.
  3. Escalón

    La salida desde reposo es 1−(1/2)ⁿ⁺¹, acercándose a uno sin alcanzarlo exactamente en un n finito del modelo ideal.

Resultado

Es un IIR causal y BIBO estable.

Lo comprobamos

La duración infinita no es una prueba de inestabilidad.

Una recurrencia puede calcular un FIR

y[n]=y[n−1]+x[n]−x[n−2], con y[−1]=0 y x[n]=0 para n<0.

  1. Sumar incrementos

    y[n]=Σ desde k=0 hasta n de [x[k]−x[k−2]].
  2. Cancelar

    Todos los términos interiores se cancelan y queda y[n]=x[n]+x[n−1].

  3. Respuesta

    h=[1;1], finita, aunque la implementación use y[n−1].

    Una memoria inicial distinta añadiría una parte libre que no se debe ignorar.

Resultado

La clasificación FIR se aplica a la respuesta de estado cero efectiva.

Lo comprobamos

Una cancelación en una realización recursiva merece revisar errores de memoria, además de simplificar la relación externa.

4. Conservar constantes exige suma uno; promediar además exige pesos no negativos

Para x[n]=c en todos los índices consultados, y[n]=cΣh[k], cuando la suma existe. La suma de pesos es la ganancia de frecuencia cero. Si debe conservar cualquier constante, esa suma tiene que ser uno. Esto no promete conservar variaciones rápidas.

Para un promedio convexo de los datos, además los pesos son no negativos. Entonces la salida queda entre el menor y mayor dato de la ventana. Con suma uno y pesos de signos distintos se preservan constantes, pero pueden aparecer valores fuera de ese intervalo.

Dos ventanas distintas conservan el mismo nivel constante

Entrada seis en todos los enteros; comparar pesos [1/3;1/3;1/3] y [1/4;1/2;1/4].

  1. Primera ventana

    6/3+6/3+6/3=6.
  2. Segunda ventana

    6/4+6/2+6/4=6.
  3. Razón

    Ambas suman uno y son no negativas. No por ello tendrán la misma respuesta ante un cambio o una oscilación.

Resultado

La conservación de continua no determina todo el filtro.

Lo comprobamos

La entrada es constante también en la prehistoria; no se ha extendido por cero al inicio de un registro.

Suma uno no basta para ser un promedio convexo

h=[2;−1].

  1. Constantes

    Para x=c, y=2c−c=c.
  2. Una ventana distinta

    Con dato actual uno y anterior cero, y=2·1−0=2.

  3. Intervalo

    Dos queda fuera de [0,1], aunque los pesos sumen uno. El peso negativo permite esa salida.

Resultado

Conservar constantes y suavizar como un promedio convexo son propiedades distintas.

Lo comprobamos

La suma absoluta de esos pesos es tres, no uno: su cota BIBO también refleja esa amplificación posible.

5. Suavizar reparte cambios y puede introducir un retraso de fase

El promedio causal de dos datos es y[n]=(x[n]+x[n−1])/2.

Combina el presente con el pasado. Conserva constantes y puede cancelar alternancias, pero una transición se extiende durante varias salidas. No conserva cada detalle del registro.

El promedio conserva una zona constante y extiende las transiciones
Un cambio se reparte en el tiempoEntrada: [0;0;2;2;2;0]20246Promedio: (x[n]+x[n−1])/220246Los bordes ocupan más de una muestra.
La entrada desde cero es [0,0,2,2,2,0] y se declara cero fuera del registro. El promedio causal de dos muestras produce [0,0,1,2,2,1,0] en índices cero a seis. La subida y la bajada aparecen repartidas, y el tramo donde ambas muestras consultadas valen dos conserva ese valor. Ambas gráficas usan las mismas escalas.

Para pesos reales simétricos de longitud M, la respuesta tiene un factor exp(−iΩ(M−1)/2) por una función real. En intervalos donde H no se anula y se elige una fase coherente, el retardo de grupo es (M−1)/2 muestras. Puede ser semientero: describe la pendiente de la fase, no una consulta a x en un índice fraccionario.

Ver una subida y una bajada después de promediar

x=[0;0;2;2;2;0] desde cero y cero fuera.

  1. Antes del cambio

    y[0]=y[1]=0.
  2. Al entrar y mantener

    y[2]=(2+0)/2=1; y[3]=y[4]=2.
  3. Al salir

    y[5]=(0+2)/2=1 y y[6]=0.

    La lista completa en cero a seis es [0;0;1;2;2;1;0].

Resultado

Los bordes se reparten y la zona verdaderamente constante se conserva.

Lo comprobamos

No se omite la salida que aún usa el último dato no nulo del registro.

Medio periodo de muestra no significa un índice no entero

Entrada bilateral x[n]=cos(πn/2), promedio de dos.

  1. Primeras salidas

    El patrón de salida es 1/2,1/2,−1/2,−1/2 y se repite.

  2. Forma senoidal

    Equivale a (√2/2)cos(πn/2−π/4): atenuación √2/2 y fase −π/4.

  3. Lectura temporal

    En esta frecuencia, el desfase corresponde a media muestra. La implementación sigue consultando solo x[n] y x[n−1], ambos de índice entero.

Resultado

La descripción de fase y la regla de consulta cuentan aspectos compatibles.

Lo comprobamos

No se necesita inventar el valor general de x[n−1/2] para aplicar ese promedio.

6. La diferencia elimina constantes y destaca los puntos de cambio

La primera diferencia es d[n]=x[n]−x[n−1], con pesos [1;−1].

Suman cero, así una región constante produce salida cero. Un cambio positivo o negativo deja una respuesta con su signo. Esto también puede resaltar variaciones rápidas de ruido.

Si los datos vienen de una función suave separada por T, dividir entre T da una aproximación de derivada con unidades de amplitud por tiempo. Una diferencia finita no es el límite de la derivada para un T cualquiera. El error depende de regularidad, paso y ruido.

Un pulso deja dos cambios de signo contrario

x=[0;0;2;2;2;0] desde cero y cero fuera.

  1. Subida

    d[2]=2−0=2.
  2. Tramos constantes

    d[0]=d[1]=d[3]=d[4]=0.
  3. Bajada

    d[5]=0−2=−2 y d[6]=0.

    Se obtiene [0;0;2;0;0;−2;0] en cero a seis.

Resultado

La salida marca dónde cambia el nivel, no el nivel mantenido dos.

Lo comprobamos

La suma de los dos cambios es cero porque la entrada vuelve a su nivel inicial.

El intervalo de muestreo cambia unidades y precisión

x(t)=t², T=1/10 y n=3.

  1. Datos

    t₃=0,3 y t₂=0,2; x[3]=0,09 y x[2]=0,04.
  2. Diferencia y cociente

    La diferencia es 0,05. Dividir entre T da 0,5, mientras la derivada actual 2·0,3 vale 0,6.

  3. Interpretar error

    Para esta cuadrática el cociente coincide con la derivada en el punto medio 0,25. No es una igualdad universal en el instante actual; sin dividir por T, tampoco tendría las mismas unidades.

Resultado

El filtro de diferencias y una derivada exacta se distinguen conservando paso y referencia temporal.

Lo comprobamos

Un error de amplitud en muestras consecutivas puede amplificarse al restar y dividir por un T pequeño.

7. El espectro explica qué oscilaciones cambian de tamaño y fase

Para x[n]=exp(iΩn), la convolución de estado cero produce y[n]=H(Ω)exp(iΩn), con H(Ω)=Σh[k]exp(−iΩk) si la suma converge. Para FIR es una suma finita. En Ω=0 todos los factores son uno, de modo que H(0)=Σh[k].

Cada operación tiene una respuesta distinta según la frecuencia
Suavizar y destacar cambiosMódulo de H(Ω)0120π/2πAzul: promedio de dosRojo: primera diferenciaConstantes: Ω=0Alternancia: Ω=π
Entre cero y π, el promedio de dos tiene módulo cos(Ω/2), que baja de uno a cero; la primera diferencia tiene módulo 2sin(Ω/2), que sube de cero a dos. El gráfico muestra amplitud, no fase: las fases se calculan con sus fórmulas complejas y no están definidas donde la respuesta vale cero.

El módulo escala amplitudes y el argumento cambia fase. Si H=0, no se define una fase para esa frecuencia. Las gráficas de módulo por sí solas no describen el retraso ni una respuesta completa; se interpretan junto con la fórmula compleja y la prehistoria de la entrada.

El promedio de dos mantiene continua y anula alternancia

h=[1/2;1/2].

  1. Fórmula

    H=(1+exp(−iΩ))/2=exp(−iΩ/2)cos(Ω/2).
  2. Extremos

    H(0)=1; H(π)=0.

    Una constante pasa y la alternancia (−1)ⁿ se anula.

  3. Interior

    En Ω=π/2, H=(1−i)/2, módulo √2/2 y fase −π/4, coincidiendo con el ejemplo temporal.

Resultado

Una regla de promedio tiene un efecto diferente según la frecuencia.

Lo comprobamos

No se asigna −π/2 como fase de H en Ω=π cuando el resultado es cero.

La diferencia tiene el efecto contrario sobre esos extremos

h=[1;−1].

  1. Fórmula

    H=1−exp(−iΩ).
  2. Extremos

    H(0)=0 y H(π)=2: suprime constantes y duplica la amplitud de una alternancia.

  3. Interior

    En Ω=π/2, H=1+i, con módulo √2 y fase +π/4.

    El signo imaginario no coincide con el del promedio.

Resultado

Destacar cambios incluye amplificar algunas oscilaciones rápidas.

Lo comprobamos

Que el filtro sea estable no promete reducir la amplitud de cada señal.

Una ventana de tres no anula todas las frecuencias rápidas

h=[1/3;1/3;1/3].

  1. Factorizar

    H=exp(−iΩ)[1+2cosΩ]/3.
  2. Un cero

    En Ω=2π/3, cosΩ=−1/2 y H=0.
  3. Alternancia

    En Ω=π, H=1/3.

    La alternancia no se anula: queda con un tercio de amplitud. El factor real puede cambiar de signo y añadir un salto de fase de π al atravesar un cero.

Resultado

Suavizar no equivale a borrar cualquier patrón de alta frecuencia.

Lo comprobamos

El retardo de grupo uno se interpreta en intervalos sin ceros y con una rama de fase coherente.

8. El borde del registro y la discretización son parte de la definición

Una ventana puede pedir datos anteriores o posteriores al bloque conocido. Podemos calcular solo posiciones con ventana completa o declarar extensión por cero, constante, reflexión o periodicidad. Cada convención define otras señales de borde. Una constante limitada a un registro y extendida por cero ya no es constante en todo el tiempo.

En continuo, filtrar usa una integral. Al aproximarla por una suma de rectángulos se conserva su anchura T: aparecen pesos T h(kT), junto con la regla de cuadratura y sus extremos. Usar sin más las alturas h(kT) puede perder escala y unidades. El método numérico también requiere sus hipótesis y error.

Una constante en el registro puede bajar al inicio por extensión cero

Registro [4;4;4] desde cero; promedio causal de tres.

  1. Cero fuera

    Las salidas completas son [4/3;8/3;4;8/3;4/3] en índices cero a cuatro.

  2. Ventana completa

    Si solo usamos ventanas enteramente dentro del registro, queda la salida de n=2, igual a cuatro.

  3. Extensión constante

    Si declaramos cuatro también antes y después, todas las salidas valen cuatro. No hay contradicción con la ganancia de continua: hemos cambiado la entrada fuera del bloque.

Resultado

El filtro no perdió una constante que se hubiera definido en todos los índices.

Lo comprobamos

El descenso inicial de la primera convención viene de los ceros añadidos.

Las anchuras de una integral no desaparecen al usar muestras

h(v)=1 en [0,1], cero fuera; entrada continua constante dos.

  1. Integral

    La salida es ∫₀¹2dv=2.
  2. Rectángulos

    Con T=1/4 y regla izquierda usamos k=0,1,2,3.
    Cada peso es T h(kT)=1/4; la suma da cuatro veces (1/4)·2=2.
  3. Factor perdido

    Si se usaran cuatro pesos uno sin T, saldría ocho. Esta regla es exacta para esta constante, no una identidad para cualquier función variable.

Resultado

La suma necesita conservar la medida de cada aportación.

Lo comprobamos

Los extremos y cantidad de paneles pertenecen a la cuadratura declarada; no se añade una quinta muestra con el mismo peso por ver el extremo uno.

Ahora lo resolvemos de principio a fin.

EJEMPLO RESUELTO 1

Promediar dos muestras

y[n]=(x[n]+x[n−1])/2. Si x[0]=2,x[1]=6 y x[−1]=0, calcula y[0],y[1].

  1. Regla y borde

    y[n]=(x[n]+x[n−1])/2, con x[−1]=0 dado explícitamente.
  2. Índice cero

    y[0]=(2+0)/2=1.
  3. Índice uno

    y[1]=(6+2)/2=4.
  4. Resultados

    Las salidas pedidas son uno y cuatro.

  5. Respuesta impulsional

    h[0]=h[1]=1/2; esos dos pesos reproducen las cuentas.

  6. Control

    Si los dos datos consultados fueran ambos cuatro, el promedio sería cuatro. El primer borde usó cero porque así estaba definido el dato anterior.

Resultado

y[0]=1, y[1]=4.

Lo comprobamos

Sustituir h[0]=h[1]=1/2 en la convolución reproduce ambas cuentas. Una constante igual a cuatro en las dos posiciones daría cuatro.

EJEMPLO RESUELTO 2

Un filtro de diferencias

y[n]=x[n]−x[n−1]. Halla h y su ganancia de continua.

  1. Regla

    y[n]=x[n]−x[n−1].
  2. Entrada impulso

    Con x=δ, en n=0 la salida vale uno.
  3. Siguiente índice

    En n=1 vale menos uno, por el dato anterior.

  4. Respuesta

    h=[1;−1] desde cero, con cero fuera.
  5. Continua

    La suma de pesos es 1−1=0; toda entrada constante bilateral se anula.

  6. Frecuencia

    H=1−exp(−iΩ); en Ω=0 también da cero, coincidiendo con la cuenta temporal.

Resultado

h=[1,−1]; ganancia de continua cero.

Lo comprobamos

Su respuesta frecuencial es 1−e^(−iΩ). En Ω=0 vale 1−1=0, la misma ganancia obtenida en el tiempo.

Los símbolos que hemos usado

x[n], y[n]
Muestras de entrada y salida en el índice entero n.
δ[n]
Impulso unitario: uno en el origen y cero en los demás índices.
h[n]
Respuesta a δ[n] del filtro de estado cero.
k
Índice auxiliar de la suma de convolución.
M
Número positivo de muestras de un promedio móvil.
H(e^(iΩ))
Ganancia compleja del filtro a frecuencia angular discreta Ω.
FIR, IIR
Respuesta impulsional de duración finita o infinita, respectivamente.

Ahora te toca a ti.

Tienes cinco ejercicios. Inténtalos; las pistas te ayudan si te atascas.

Practicar este tema →

Para ampliar