Teorema de muestreo

[gL.edu] Este artículo recoge contribuciones de Verónica Velasco López, Benson Nketia y J.M. Díaz Nafría, elaboradas en el contexto de la Clarificación conceptual en "teoría de la señal y la comunicación", bajo la supervisión de J.M. Díaz Nafría.

Definiciones

El teorema de muestreo establece que si una señal continua en el tiempo, $x_a(t)$, tiene un ancho de banda limitado, B, es posible reconstruirla de forma exacta a partir de muestras de la señal original muestreada a una frecuencia $F_m\ge 2B$. En otras palabras, el teorema de muestreo nos indica las condiciones del muestreo para evitar que se produzca el problema de solape en frecuencia (aliasing) y así no perder información durante el proceso de muestreo.

Formalmente el teorema de muestreo puede expresarse así: Sea $x(t)$ una señal de banda limitada con $X(F)=0$ para $|F|>F_{max}$ entonces se determina unívocamente mediante sus muestras $x(nT)$ donde $n=\pm 1,\pm 2,\pm 3,...,$ si la frecuencia de muestreo $F_m>2F_{max}$, donde $F_m=1/T$, siendo T el periodo de muestreo.

Muestreo de señales periódicas

Supongamos una señal analógica compuesta por un conjunto finito de sinusoides de diferentes amplitudes, frecuencias y fases:

$x_a(t)=\sum _{i=0}^N{A}_i\cos (2\pi F_{i}t+{\theta }_i)$

Donde N es el número de componentes armónicas de la señal analógica y las frecuencias se encuentran acotadas en el intervalo $[0,F_{max}]$. De acuerdo al teorema de muestreo, para evitar pérdida de información en el proceso de muestreo de la señal original, sería necesario que la frecuencia de muestreo, $F_m$, cumpla: $F_m>2F_{max}$

Seleccionando esta frecuencia de muestreo, cualquier componente de frecuencia $F_i<F_{max}$ se corresponde con una sinusoide discreta en el tiempo de frecuencia normalizada: $f_i=\frac{F_i}{F_m}\le \frac{1}{2}\mathrm{\forall }i$. De este modo, todas las componentes de frecuencia de la señal analógica están representadas en la forma muestreada sin ambigüedad y puede reconstruirse sin distorsión a partir de los valores de las muestras, usando la función de interpolación ideal:

$g(t)=\frac{sin\left(2\pi Bt\right)}{2\pi Bt}$

De modo que: $x_a(t)=\sum _{n=-\mathrm{\infty }}^{n=\mathrm{\infty }}{x}_a\left(\frac{n}{F_m}\right)g(t-\frac{n}{F_m})$

La frecuencia mínima de muestreo que nos permite recuperar sin pérdidas la señal de ancho banda limitado es $F_N=2B=2F_{max}$, denominada Frecuencia de Nyquist.

Señales periódicas en relación armónica

Si las componentes sinusoidales de $x_a(t)$ se encontraran en relación armónica con una frecuencia fundamental $F_0=1/T$ de modo que $F_i=i{F}_0=i/T$ con $i=0,1,...N$, entonces T sería el periodo de la señal analógica y bastaría tomar muestras dentro de este periodo para poder determinar unívocamente todas las componentes armónicas. El resto de muestras serían repetitivas y, por tanto, redundantes. Esto es, si aplicamos el criterio de Nyquist antes referido ($f_m\ge 2f_{max}=2N{f}_0=2N/T$) sería suficiente un número finito de muestras para preservar toda la información de la señal original: Nº muestras $\ge T{f}_N+1=T\cdot 2N{f}_0+1=T\frac{2N}{T}+1=2N+1$ (donde el valor en exceso respecto a 2N se debe a que se toma una muestra al inicio y otra al final del periodo). La necesidad de esta cantidad de datos para la reconstrucción de la señal original se puede entender teniendo en cuenta que, de esta forma, dentro del periodo del armónico de mayor frecuencia, N/T, tendríamos dos valores que nos permitirían determinar $A_N$ y ${\theta }_N$. Por otra parte, en conjunto habrá 2N+1 valores $\{A_i{\}}_{N+1}$ y $\{{\theta }_i{\}}_N$ que deben determinarse para la reconstrucción exacta de la señal original, $x_a(t)$.

Demostración del teorema de muestreo

Antes nos hemos referido a una señal periódica representada en términos de una seria de Fourier, sin embargo, el resultado es extensible a cualquier señal sea o no periódica. La transformada de Fourier sería precisamente la expresión en el dominio de la frecuencia de la señal analógica.

La operación de muestreo podemos representarla idealmente como el producto de la señal original por un tren de deltas de Dirac separadas Tm. Pues bien, basta con tener en cuenta que el la transformada de Fourier de ese tren de deltas en el dominio del tiempo es un otro tren de deltas en el dominio de la frecuencia situadas en k/Tm= k fm, para observar -aplicando el teorema de la convolución- que: la transformada de Fourier de la señal muestreada es periódico con periodo fm, repitiendo el espectro de la señal original. Es decir,

$X_m(f)=\sum _{n=-\mathrm{\infty }}^{n=\mathrm{\infty }}{X}_a(f-n{F}_m)$

Desde este punto de vista, se observa que la única manera de que no se produzca un solape espectral es haciendo que B sea inferior a $F_m/2$, así como que bastará con aplicar un filtro paso bajo ideal con frecuencia de corte $f_m/2$ cuya respuesta impulsional es precisamente $g(t)$, la función de interpolación antes referida. En efecto la aplicación de dicho filtro hace que

$X_a(f)=X_m(f)\cdot G(f)=X_m(f)\cdot \mathrm{\Pi }(\frac{f}{f_m})=\{\sum _{n=-\mathrm{\infty }}^{n=\mathrm{\infty }}{X}_a(f-n{F}_m)\}\cdot \mathrm{\Pi }(\frac{f}{f_m})=X_a(f)\leftrightarrow x_a(t)$[qed]

Referencias

Proakis, J.G. y Manolakis, D.G. (2007). Tratamiento Digital de Señales. Madrid: Pearson Education.