MECATÓNICA 2020

Teorema del Muestreo

PRACTICA 3
González García Itzel
Igg1707@[Link]
Instituto Tecnológico Superior de Huauchinango

Resumen: En esta práctica se realizara la

demostración del teorema del muestreo al
aplicarlo en un sistema digital esto por medio del
software de apoyo MATLAB y Simulink, a esto se
agregara un bloque muestreador y un retenedor Los coeficientes vienen dados por las
de orden cero expresiones:

Palabras clave: muestreo, MATLAB, sistema,

digita, Simulink, retenedor, muestreador

Abstract: In this practice, the sampling theorem Definición 2: Por Snf(x) denotaremos la suma
will be demonstrated by applying it in a digital parcial de orden n de la serie de Fourier de f
system, this using the MATLAB and Simulink
support software, this is agreed in a sampling
block and a zero order retainer

Keywords: Sampling, MATLAB, System, Digita, Donde ak y bk vienen dados por la formula
Simulink, Retainer, Sampler anterior. Ver Fig. 1 Ejemplo de una serie de
Fourier.
I. INTRODUCCIÓN

Fourier comienza a partir de principios físicos

deduciendo la ecuación que gobierna la difusión del
calor. Después resuelve el problema de la
distribución de temperatura en un tiempo dado a
partir de la distribución en el instante inicial, en
varios casos. Para ello, utiliza el método de
separación de variables, que introdujo d’Alembert
en su estudio de la ecuación de ondas. Para escribir Fig. 1 Serie de Fourier de x2
la solución, necesita escribir la función que da el
dato inicial como suma de una serie trigonométrica. II. DESARROLLO
Esta es precisamente la principal herramienta que a) Realizamos el diagrama en Simulink.
usaremos en nuestro trabajo y su teoría a día de hoy Ver Fig. 2
se conoce como análisis de Fourier, que en su
momento generó el análisis moderno.

Definición: Dada una función f de cuadrado

integrable en [−π, π] y 2π-periódica en [−π, π]
definimos la serie trigonométrica de Fourier como:

Asociación Mexicana de Mecatrónica A.C.

 MECATÓNICA 2020

Fig. 2 Diagrama de Simulink Fig. 3 Configuración de parámetros para el bloque “Sine Wave”

Para este diagrama ocupamos las herramientas:

muestreador y un retenedor de orden cero.
El bloque de retención de orden cero muestra y
retiene su entrada durante el período de muestra
especificado. El bloque acepta una entrada y
genera una salida, las cuales pueden ser escalares
o vectoriales. Si la entrada es un vector, todos los
elementos del vector se mantienen durante el
mismo período de muestra.
El bloque de onda sinusoidal proporciona una
sinusoide. El bloque puede funcionar en modo
basado en tiempo o en modo muestra.

b) Configuramos los tiempos de Fig.4 Configuración de parámetros del bloque “Pulse

muestreo de la simulación a “ode45 Generator”
(Dormand-Prince)”. También
configuramos los parámetros de los Por ultimo para los bloques de ZOH, que están
bloques “Sine Wave” y “Pulse definidos como sintonizadores colocamos el
Generator”. Ver las Fig. 3 Y 4. parámetro “rTs” y el bloque que está definido como
ZOH cambiamos el parámetro a “Ts”

c) Después en un Scrip de Matlad

escribimos el siguiente código
w1=1;
Tend = 100;
Tmax = 20;
Ts=1;
rTs=0.01;
ws=(2*pi)/Ts;
dw1=2*w1;
fN =pi/Ts;
sim('Practica1CD',Tend)
figure()

Asociación Mexicana de Mecatrónica A.C.

 MECATÓNICA 2020

subplot 311
stem(t,p)
xlim([0 Tmax])
title('Muestreador')
subplot 312
plot(t,xt)
hold on
stem (t,xs)
hold off
xlim([0 Tmax])
title('Señales originales y muestreadas')
subplot 313
tz = 0:Ts:Tend;
xz2 = interp1(tz,xz,t,'previous');
plot(t,xz2)
xlim([0 Tmax])
title('Señal retenida') Fig. 6 Grafica para un Ts=0.5 (alto)

Al correr el programa las gráficas en cuatro

escenarios, muestreo suficiente, muestreo
alto, muestreo insuficiente y frecuencia
múltiplo; son las siguientes. Ver Fig. 5, 6,7
y 8.

Fig. 7 Grafica para un Ts=4 (Insuficiente)

Fig. 5 Graficas para un Ts =1 (muestreo suficiente)

Fig. 8 Grafica para un Ts=3 (frec. Múltiplo)

d) Al agregar esta parte al código anterior,

nos genera el intervalo real. Ver Fig. 9

N = length(xt)
k=1:N;

Asociación Mexicana de Mecatrónica A.C.

 MECATÓNICA 2020

w=2*pi*(k-1)/(rTs*N); f) Al final agregamos este código para graficar

w=w'; en el dominio de la frecuencia para los
cuatro escenarios antes mencionados. Ver
Fig. 10, 11, 12 y 13.

wmax=5*pi
figure()
subplot 311
plot(w,abs(XT))
xlim([0 wmax])
title('Señal oriinal en el dominio de la
frecuencia')
subplot 312
plot(w,abs(XS))
xlim([0 wmax])
Fig. 9 Intervalo Real generado title('Señal muestreada en el dominio de la
frecuencia')
subplot 313
e) Después con este código obtuvimos la plot(w,abs(XZ))
transformada de Fourier. Ver tabla 1 xlim([0 wmax])
XT=fft(xt) title('Señal retenida en el dominio de la
XS=fft(xs) frecuencia')
XZ=fft(xz2)
Tabla 1. Datos obtenidos de la Trasformada
-0.0000 - 0.0001i
-0.0002 - 0.0006i -0.0000 - 0.0001i
-0.0002 - 0.0005i
0.0000 - 0.0001i
-0.0002 - 0.0005i
0.0000 - 0.0001i
-0.0002 - 0.0005i
0.0000 - 0.0001i
-0.0001 - 0.0004i
0.0000 - 0.0001i
-0.0001 - 0.0004i 0.0000 - 0.0001i
-0.0001 - 0.0004i
0.0000 - 0.0001i
-0.0001 - 0.0004i
0.0000 - 0.0001i
-0.0001 - 0.0004i
0.0000 - 0.0001i
-0.0001 - 0.0003i
0.0000 - 0.0001i
-0.0001 - 0.0003i
0.0000 - 0.0001i
-0.0001 - 0.0003i
0.0000 - 0.0001i
-0.0001 - 0.0003i
0.0000 - 0.0001i
-0.0001 - 0.0003i
0.0000 - 0.0001i
-0.0001 - 0.0003i
0.0000 + 0.0001i
-0.0001 - 0.0003i
0.0000 + 0.0001i
-0.0001 - 0.0003i
0.0257 - 0.0028i
-0.0001 - 0.0003i
-0.0000 - 0.0002i
-0.0000 - 0.0002i Fig. 10 Grafica en dominio de frecuencia para Ts=1
-0.0000 - 0.0002i
-0.0000 - 0.0002i
-0.0000 - 0.0002i
-0.0000 - 0.0002i
-0.0000 - 0.0002i
-0.0000 - 0.0002i
-0.0000 - 0.0002i
-0.0000 - 0.0002i
-0.0000 - 0.0002i
-0.0000 - 0.0002i
-0.0000 - 0.0002i
-0.0000 - 0.0002i
-0.0000 - 0.0001i
-0.0000 - 0.0001i
-0.0000 - 0.0001i
-0.0000 - 0.0001i
-0.0000 - 0.0001i
-0.0000 - 0.0001i
-0.0000 - 0.0001i
-0.0000 - 0.0001i

Fig. 11 Grafica en dominio de frecuencia para Ts= 0.5

Asociación Mexicana de Mecatrónica A.C.

 MECATÓNICA 2020

matemático subyacente con anterioridad, pero

en magnífica soledad.

Referencias

[1] Renato Alvarez Nodarse. Introducción al Máxima

CAS con algunas aplicaciones. Universidad de
Sevilla, 2018.
[Link]
lcurso/ [Link].
[2] Renato Álvarez Nodarse. Polinomios
hipergeométricos y q-polinomios. Monografías
del Seminario Maten ático García de Galeano,
Fig. 12 Grafica en dominio de frecuencia para Ts=4 número 26. Prensas Universitarias de Zaragoza,
2003.
[3] J. Duoandikoetxea. Lecciones sobre las series y
transformadas de Fourier. UNAM-Mangua,
2003.
[Link]
cas/
analisisdefourier/[Link].
[4] José Antonio Facenda Aguirre y Francisco José
Freniche Ibáñez. Integración de funciones de
varias variables. PIRAMIDE, 2002. ´
[5] Deguang Han, Keri Kornelson, David Larson, Eric
Weber. Frames for Undergraduates. American
Mathematical Society, 2007.
[6] Jesús Hernández Alonso. 200 años de
convergencia de las series de Fourier. La Gaceta
Fig.13 Grafica en dominio de frecuencia pata Ts=3
de la RSME, Vol. 10.3, 2007.
[Link]
C3%B1os_ convergencia_series_Fourier.pdf.

III. CONCLUSIONES
El teorema del muestreo para funciones tipo
pasa-bajos es vital en la ingeniería
de comunicaciones, vinculando las señales
continuas con las discretas.

Los distintos nombres que el teorema ha

recibido Shannon, Nyquist, Kotelnikov,
Whittaker, Someya, etc. indica
la independencia de los diferentes enfoques.
Esta historia también pone de manifiesto
un proceso muchas veces aparente en
los problemas teóricos en física o tecnología:
en primer lugar los experimentadores
establecen un conjunto de reglas empíricas,
luego los teóricos desarrollan la solución
general, y finalmente alguien descubre que los
matemáticos habían resuelto el problema

Asociación Mexicana de Mecatrónica A.C.