CONTROL DIGITAL
WILLIAM EDISON VANEGAS RODRIGUEZ
C.C: 5820973
Tutor: Ing. Joan Sebastián Bustos
GRUPO 203041_19
UNIVERSIDAD NACIONAL ABIERTA Y A DISTANCIA - UNAD
ESCUELA ECBTI
PROGRAMA INGENIERIA ELECTRONICA
CEAD, IBAGUE SEPTIEMBRE 2020
DESARROLLO ACTIVIDAD
Usando como base la temática de la unidad uno (1) del curso, la cual encontrará en el
entorno de aprendizaje y con el uso de Matlab y Proteus dar solución a los siguientes
puntos del proyecto:
[Link] sobre la identificación de modelos a partir de la curva de reacción de un
proceso y realizar un pequeño resumen sobre los modelos generales. (Máximo 4 páginas en
este punto).
[Link] el montaje en el software Proteus de la figura 1, analizar la curva de reacción del
Heated Oven (Lumped Model) el cual se debe configurar con los siguientes parámetros:
Temperature Ambient (°C) =27
•Thermal Resistence to Ambient (°C/W) =0.5
•Oven Time Constant (sec)= 10
•Heater Time Constant (sec)= 1
•Temperature coefficient (V/°C) =1
•Heating Power (W)= 150
A partir de la curva de reacción calcular lo siguiente:
•Establecer el modelo analítico (ecuación matemática) según la gráfica del simulador
Proteus.
•Simular el modelo analítico (ecuación matemática) en Matlab y aplicar una entrada tipo
escalón.
•Hacer el análisis de estabilidad de manera analítico y computacional y comparar resultados
El sistema resultante es de primer orden por la forma que tiene la gráfica obtenida en proteus.
t 0=2 seg
t 1=3 seg
y 1=840 ° C
u1=10 V
y 0=27 ° C
t 2=15,2 seg
τ 0=t 1−t 0 =3−2=1 seg
γ 0=t 2−t 1=15.2−3=12.2 seg
y 1− y 0 840−27 813
k 0= = = =81.3
u 1−u0 10−0 10
K 0 e−st 0 e−s 1
G ( s )= =81.3
1+ γoS 1+12.2 s
Análisis Matemático
Se tiene en cuenta que una función de transferencia de primer orden tiene la siguiente forma:
C( s) K
=
R( s) τ s +1
Donde K es la ganancia del sistema y 𝛕 la constante de tiempo del sistema
K e−Ls
En nuestro caso sería de esta forma: G ( s )=
τs+1
Donde L es el tiempo que tarda en reaccionar el sistema en el momento de aplicar el impulso,
entonces tenemos que:
e−s 1
Y ( s )=81.3
1+ 12.2 s
−t
(
y ( t ) =K ∆ u 1−e τ )
−t
y ( t ) =29.9 (1−e )τ
τ =12.2 seg
Si el tiempo es igual a 𝛕 tenemos que:
t=τ
K ∆ u=( 1−e−1 ) =0.632 ( K ∆ u )=18.96
Se tiene en cuenta la variación de temperatura la cual es ∆ T =840−27=813 °C
En la gráfica se puede apreciar que, como es de primer orden y la entrada al sistema es un paso
escalón es aproximadamente de 63.2 % entonces 8.13∗63.2=513.816 ° C
En Matlab se aprecia que en τ =12.2 s llega a 489 ° C aprox
En 2𝛕 86.5% entonces 8 . 13∗86.5=703.24 ° C
En 3𝛕 95% entonces 8.13∗95=772.35° C
En 4𝛕 98.2 % entonces 8.13∗98.2=798.366 ° C
En 5𝛕 99.3 % entonces 8.13∗99.3=807.309° C
Por definición tenemos que:
t d=1 seg tiempo retardo
t r=12.2∗2=24.4 seg tiempo subida
t p=12.2∗3=36.6 seg tiempo pico
t s=12.2∗4=48.8 seg tiempo asentamiento
M p=0
Por ser un sistema de primer orden, la sobre elongación es casi nula por lo tanto se desprecia