FACULTAD DE INGENIERÍA DE
PRODUCCIÓN Y SERVICIOS
ESCUELA PROFESIONAL DE INGENIERÍA
ELECTRÓNICA
Curso:
SISTEMAS DE CONTROL DIGITAL
Tema:
MUESTREO CON MATLAB Y SIMULINK
Elaborado por:
Roque Quispe Alejandra
Roque Ramos Ronald
Arequipa - Perú
2019
1. Muestreo con MatLab y Simulink.
Obtener la señal muestreada por MatLab y Simulink del siguiente esquema:
y(t) y(kT)
y ( t ) =sen ( wt ) f m=1 khz frecuencia de muestreo
y ( t ) =Asen ( ωt +α ) : entrada
y ( kT )= Asen ( ωkT + α ) y ( kT )= y (t=kT ): salida, k es No. entero= 0,1, 2, 3……….
a . ω=2 πf , siendo f =200 , 400 y 500 Hz , α =0 , A=1
f =200 Hz , α =0 , A=1
ω=2 πf , siendo f =500 Hz , α =0 , A=1
b. f =200 hz , α=45 ° , A=1
2. El objetivo del presente item es comparar las respuestas de dos
plantas idénticas, con la única diferencia en que la señal de entrada
(una señal sinusoidal de frecuencia ω1) en una entra directamente y
en la segunda opción la señal es previamente muestreada y luego
reconstruida a través de un bloqueador de orden cero como se muestra
en la Figura 1. Utilice un tiempo de muestreo T = 1 seg.
Analizar que ocurre al variar la frecuencia de la señal sinusoidal de entrada con:
ω1 = ωs/100, ωs/20, ωs/5, ωs/4.5, ωs/4.0, ωs/3.5, ωs/3.0, ωs/2.5, ωs/2.0 (probar para esta
frecuencia además que ocurre al variar el desfasaje), ωs/1.5 y ωs.
Frecuencia de muestreo ωs. = 2.π/T.
Explicar. ¿El tiempo de muestreo elegido es el adecuado?
ω1 = ωs/100, Frecuencia de muestreo ωs. = 2.π/. Variar el desfasaje, ωs/1.5 y
ωs
ω1 = ωs/20, Frecuencia de muestreo ωs. = 2.π/. Variar el desfasaje, ωs/1.5 y ωs
ω1 = ωs/5, Frecuencia de muestreo ωs. = 2.π/. Variar el desfasaje, ωs/1.5 y ωs
ω1 = ωs/4.5, Frecuencia de muestreo ωs. = 2.π/. Variar el desfasaje, ωs/1.5 y ωs
ω1 = ωs/4.0, Frecuencia de muestreo ωs. = 2.π/. Variar el desfasaje, ωs/1.5 y ωs
ω1 = ωs/3.5, Frecuencia de muestreo ωs. = 2.π/. Variar el desfasaje, ωs/1.5 y ωs
ω1 = ωs/3.0, Frecuencia de muestreo ωs. = 2.π/. Variar el desfasaje, ωs/1.5 y ωs
ω1 = ωs/2.5, Frecuencia de muestreo ωs. = 2.π/. Variar el desfasaje, ωs/1.5 y ωs
ω1 = ωs/2.0, Frecuencia de muestreo ωs. = 2.π/. Variar el desfasaje, ωs/1.5 y ωs
3. Utilizando el mismo esquema anterior, cambiar la fuente de señal por un
escalón unitario. ¿Se observa alguna diferencia entre ambos? ¿Por qué?
Observando la gráfica se observa que para una entrada escalón la señal reconstruída es
idéntica a la señal original. La diferencia de una entrada sinusoidal a una escalón es que
una es periódica y la otra no; por lo cual, el escalón no presenta muestras.
4. Para la planta continúa utilizada, obtener la función de transferencia discreta que representa
la misma (utilizando un bloqueador de orden cero).
Para esto MATLAB dispone de una serie de funciones para realizar las conversiones
temporales del dominio discreto al continuo y viceversa. Las funciones más importantes se
muestran en la Figura 2. Aquéllas cuyo nombre termina en m, poseen un método de
conversión a elegir por el usuario. El resto, por defecto, utilizan el sistema ZOH (bloqueador
de orden cero). En todos ellos un parámetro importante es el periodo de muestreo. Para
calcular el equivalente discreto de una planta continua con bloqueador de orden cero y
salida muestreada se emplea el comando c2d.
Figura 2. Conversión de modelos en MATLAB
Observe la relación entre los polos de la planta continua original y los polos de la función de
transferencia discreta obtenida (recuerde para eso la relación entre transformada S y
transformada Z). Compruebe que no sucede lo mismo con el cero.
Discretización de sistemas continuos con Matlab.
El paquete Matlab posee funciones que permiten discretizar un sistema continuo mediante
distintos métodos de discretización. El comando que realiza dicha tarea es c2dm. Dicho
comando tiene la siguiente sintaxis:
[numd,dend] = c2dm(num,den,Ts,'method')
De esta forma se convierte la función de transferencia continua G(s)=num(s)/den(s) en la
función de transferencia discreta G(z)=numd(z)/dend(z) a un periodo de muestreo Ts
utilizando el método ‘method’, donde num, den, numd y dend son los vectores cuyos
elementos son los coeficientes en orden decreciente de los polinomios del numerador y
denominador de las funciones de transferencia continua y discreta respectivamente.
Los distintos métodos de discretización que podemos utilizar son los siguientes:
'zoh' : discretiza asumiendo que existe un bloqueador (retenedor) de orden cero a la
entrada.
'foh' : discretiza asumiendo que existe un bloqueador (retenedor) de orden uno a la
entrada.
'tustin' : discretiza utilizando la transformación bilineal.
'prewarp' : discretiza utilizando la transformación bilineal con compensación de
frecuencias (prewarping). En este caso es necesario especificar la frecuencia de
corte y el comando quedaría:
[numd,dend]=c2dm(num,den,Ts,'prewarp',Wc), donde Wc es la frecuencia
de corte en rad./seg.
'matched' : discretiza utilizando el método de asignación de polos y ceros.
5. Realice un esquema de simulación en Simulink similar al del punto 2 pero agregando la
función transferencia discreta obtenida en el punto 4. Para esto utilice el conjunto de
bloques de la librería “Discrete” de Simulink.
Compare las respuestas de los tres esquemas de simulación.
Solucion:
CUESTIONARIO FINAL6
Realice todo el procedimiento anteriormente descrito y presente en un informe el desarrollo de
todos los pasos indicados.
OBSERVACIONES Y CONCLUSIONES
Haga sus observaciones y emita al menos cinco conclusiones en torno al trabajo realizado
Para el muestreo y reconstrucción de una señal se debe de cumplir la tasa de Nyquist
(frecuencia de muestreo >= 2 ω).
El muestreo de una señal es importante para la reconstrucción de ondas y también para
conocer los errores que puede tener la señal muestreada.
Se pudo observar como la señal reconstruida no coincide completamente con la original.
Entre menor sea el periodo de muestreo la señal reconstruida se parece más a la señal
original o de entrada.2
Pudimos observar a simple vista que el bloqueador reconstruye mejor la señal de partida o
entrada.
