Fórmula de interpolación de Whittaker-Shannon

Algoritmo de (re)construcción de señales

La fórmula de interpolación de Whittaker-Shannon o interpolación sinc es un método para construir una función limitada en banda de tiempo continuo a partir de una secuencia de números reales. La fórmula se remonta a los trabajos de E. Borel en 1898 y ET Whittaker en 1915, y fue citada en los trabajos de JM Whittaker en 1935 y en la formulación del teorema de muestreo de Nyquist-Shannon por Claude Shannon en 1949. También se llama comúnmente fórmula de interpolación de Shannon y fórmula de interpolación de Whittaker . ET Whittaker, que la publicó en 1915, la llamó serie Cardinal .

Definición

En la figura de la izquierda, la curva gris muestra una función f(t) en el dominio del tiempo que se muestrea (los puntos negros) a velocidades de muestreo en constante aumento y se reconstruye para producir la curva dorada. En la figura de la derecha, la curva roja muestra el espectro de frecuencia de la función original f(t), que no cambia. La frecuencia más alta del espectro es la mitad del ancho de todo el espectro. El sombreado rosa que aumenta constantemente representa el espectro de frecuencia de la función reconstruida, que gradualmente llena más espectro de frecuencia de la función original a medida que aumenta la frecuencia de muestreo. Cuando el espectro de frecuencia de la función reconstruida abarca todo el espectro de frecuencia de la función original, es dos veces más ancho que la frecuencia más alta, y es entonces cuando la forma de onda reconstruida coincide con la muestreada.

Dada una secuencia de números reales, x [ n ], la función continua

${\displaystyle x(t)=\sum _{n=-\infty }^{\infty }x[n]\,{\rm {sinc}}\left({\frac {t-nT}{T}}\right)\,}$

(donde "sinc" denota la función sinc normalizada ) tiene una transformada de Fourier , X ( f ), cuyos valores distintos de cero se limitan a la región | f | ≤ 1/(2T ) . Cuando el parámetro T tiene unidades de segundos, el límite de banda , 1/(2 T ), tiene unidades de ciclos/seg ( hercios ). Cuando la secuencia x [ n ] representa muestras de tiempo, en el intervalo T , de una _función continua, la cantidad fs = 1/ T se conoce como frecuencia de muestreo , y fs /2 es la frecuencia de Nyquist _{correspondiente} . Cuando la función muestreada tiene un límite de banda, B , menor que la frecuencia de Nyquist, x ( t ) es una reconstrucción perfecta de la función original. (Consulte el teorema de muestreo ). De lo contrario, los componentes de frecuencia por encima de la frecuencia de Nyquist se "pliegan" en la región sub-Nyquist de X ( f ), lo que produce distorsión. (Ver Alias ).

Formulación equivalente: filtro de convolución/paso bajo

La fórmula de interpolación se deriva del artículo sobre el teorema de muestreo de Nyquist-Shannon , que señala que también se puede expresar como la convolución de un tren de impulsos infinito con una función sinc :

${\displaystyle x(t)=\left(\sum _{n=-\infty }^{\infty }T\cdot \underbrace {x(nT)} _{x[n]}\cdot \delta \left(t-nT\right)\right)\circledast \left({\frac {1}{T}}{\rm {sinc}}\left({\frac {t}{T}}\right)\right).}$

Esto equivale a filtrar el tren de impulsos con un filtro de paso bajo ideal ( pared de ladrillo ) con una ganancia de 1 (o 0 dB) en la banda de paso. Si la frecuencia de muestreo es suficientemente alta, esto significa que la imagen de banda base (la señal original antes del muestreo) pasa sin cambios y el filtro de pared de ladrillo elimina las otras imágenes.

Convergencia

La fórmula de interpolación siempre converge absoluta y localmente de manera uniforme siempre que

${\displaystyle \sum _{n\in \mathbb {Z} ,\,n\neq 0}\left|{\frac {x[n]}{n}}\right|<\infty .}$

Por la desigualdad de Hölder esto se cumple si la secuencia pertenece a cualquiera de los espacios con 1 ≤  p  < ∞, es decir ${\displaystyle (x[n])_{n\in \mathbb {Z} }}$ ${\displaystyle \ell ^{p}(\mathbb {Z} ,\mathbb {C} )}$

${\displaystyle \sum _{n\in \mathbb {Z} }\left|x[n]\right|^{p}<\infty .}$

Esta condición es suficiente, pero no necesaria. Por ejemplo, la suma generalmente convergerá si la secuencia de muestras proviene del muestreo de casi cualquier proceso estacionario , en cuyo caso la secuencia de muestras no es sumable al cuadrado y no está en ningún espacio. ${\displaystyle \ell ^{p}(\mathbb {Z} ,\mathbb {C} )}$

Procesos aleatorios estacionarios

Si x [ n ] es una secuencia infinita de muestras de una función muestral de un proceso estacionario de sentido amplio , entonces no es miembro de ningún espacio L p , con probabilidad 1; es decir, la suma infinita de muestras elevada a una potencia p no tiene un valor esperado finito. Sin embargo, la fórmula de interpolación converge con probabilidad 1. La convergencia puede demostrarse fácilmente calculando las varianzas de los términos truncados de la sumatoria y demostrando que la varianza puede hacerse arbitrariamente pequeña eligiendo un número suficiente de términos. Si la media del proceso es distinta de cero, entonces es necesario considerar pares de términos para mostrar también que el valor esperado de los términos truncados converge a cero. ${\displaystyle \ell ^{p}}$

Dado que un proceso aleatorio no tiene transformada de Fourier, la condición bajo la cual la suma converge a la función original también debe ser diferente. Un proceso aleatorio estacionario tiene una función de autocorrelación y, por tanto, una densidad espectral según el teorema de Wiener-Khinchin . Una condición adecuada para la convergencia a una función de muestra del proceso es que la densidad espectral del proceso sea cero en todas las frecuencias iguales y superiores a la mitad de la frecuencia de muestreo.

Ver también

• Aliasado , Filtro antialiasing , Antialiasing espacial
• función rectangular
• Muestreo (procesamiento de señales)
• Señal (electrónica)
• Función Sinc , filtro Sinc
• Remuestreo de Lanczos