Diseño de Controladores con Matlab
Diseño de Controladores con Matlab
- Diseño de Controladores
Controlador
Acción de Planta
R(s) E(s) control U(s) Y(s)
+ G ( s)
- Gc ( s )
Y1(s) Sensor
H (s)
A los controladores que usan el esquema de la figura 11.1 también se les llama controladores
basados en la señal de error actuante, la cual es la diferencia entre la referencia deseada r (t ) y
la salida medida de la planta y1 (t ) , es decir,
e(t ) = r (t ) − y1 (t ) (11.1)
o bien,
E ( s ) = R ( s ) − Y1 ( S ) (11.2)
En la mayoría de los casos se usan sensores tales que G ( s ) = 1 , por lo tanto, es común suponer
que y1 (t ) = y (t ) .
Las acciones de clásicas de control usadas en el esquema de la figura 11.1 se pueden clasificar
como sigue:
A continuación se hace una breve descripción de cada una de estas acciones de control.
Esta acción de control solamente puede producir uno de dos valores: U OFF , o U ON , los cuales
conmuta y envía a la entrada de la planta de acuerdo con el valor de la señal de error actuante. La
versión más sencilla de esta acción de control se muestra en la parte izquierda de la figura 11.2 y
está dada por
U si e(t ) > ET
u (t ) = ON (11.3)
U OFF si e(t ) < ET
Sin embargo, esta implementación produce cambios muy bruscos que pueden dañar los
elementos del actuador, por esta razón se le agrega una brecha diferencial en la forma de un
efecto de histéresis de manera que una implementación más adecuada está dada por
U ON si e(t ) > EH
u (t ) = U OFF si e(t ) < EL (11.4)
no cambia si E < e(t ) < E
L H
u(t) u(t)
VON VON
VOFF VOFF
ET e(t) EL EH
Conmutación simple e(t)
Conmutación con Histéresis
Figura 11.2.- Efecto Entrada- Salida de la Acción de control ON-OFF
Para implementar esta acción de control en Matlab se Simulink se cuenta con el bloque 'relay'
mostrado en la figura 11.4, el cual se puede configurar dando doble click para que aparezca el
cuadro de diálogo mostrado en la misma figura.
Para que la configuración mostrada en el cuadro de diálogo de la figura 11.4 coincida con el
comportamiento descrito por la ecuación (11.4), se deberán introducir los siguientes valores:
Switch on point: EH
Switch off point: EL
Output when on: U ON
Output when off: U OFF
Objetivo de control: Se desea que al cerrar el lazo con el controlador, la velocidad del
motor se mantenga a la referencia de velocidad r (t ) = constante, aún en presencia de
perturbaciones.
Scope
Perturbacion
= -70rad/seg
15.7
0.02s+1
Ref=100r/s Relay 0 a 12v Gm(s)
h=20
Figura 11.5.- Diagrama de simulación del control ON-OFF para la velocidad del motor de cd
Condiciones de la simulación:
Referencia = 100 rad/seg
Perturbación constante = -70 rad/seg a partir del instante t = 0.1 seg
parámetros del controlador ON-OFF: EH = 10 , EL = −10 , (h = 20) , U ON = 12 volts ,
U OFF = 0 volts
Los resultados de la simulación se muestran en la figura 11.6, en la cual se puede apreciar que la
brecha diferencial (h) evita un seguimiento exacto de la referencia, pero hace que la acción de
control no conmute tan rápido Si se disminuye h se tendrá un mejor seguimiento, pero a costa de
hacer más rápida la conmutación. También se observa que a partir de la perturbación la
conmutación se hace más lenta, pero la acción de control Va (t ) le cuesta más "esfuerzo", ya que
dura más tiempo en su valor máximo de 12 volts.
100 h
ω(t) (rad/seg)
80
60
40
20
0
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2
t (seg)
Va (volts)
10
5
0
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2
t (seg)
Figura 11.6.- Acción de control Va (t ) ON-OFF y respuesta de velocidad ω (t ) del motor de CD.
Esta acción de control produce una corrección proporcional al tamaño del error, es decir, el valor
de la acción de control se obtiene multiplicando el error por una constante, es decir,
u (t ) = K p e(t ) (11.6)
donde K p se denomina ganancia proporcional.
De la ecuación anterior se observa que el error en estado estable nunca es cero, pero disminuye
para valores grandes de K p . En la figura 11.8 se muestra la respuesta del control proporcional
para la velocidad del motor de CD con diferentes valores de K P .
90 ess
80
70 Referencia rad/seg
velocidad (rad/seg) Kp=0.5
60 Velocidad (rad/seg) Kp=1
Velocidad (rad/seg) Kp=5
50
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2
0
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2
Figura 11.8.- Acción de control proporcional Va (t ) y respuesta de velocidad ω (t ) del motor de CD
En la figura 11.8 se aprecia que el error de estado estable ess (el cual se marca sólo para el caso
K p = 0.5 ) disminuye conforme K p aumenta, además también disminuye el efecto de la
perturbación, sin embargo, esta mejora se logra con el costo de que la acción de control aumenta
considerablemente, sobre todo en el transitorio (La gráfica está recortada para apreciar el
transitorio de la acción de control Va (t ) , ya que los picos de voltaje llegan hasta 500 volts).
Una alternativa para disminuir el error de estado estable sin que la acción de control tenga los
picos producidos por la acción de control Integral, la cual produce una corrección proporcional a la
integral de la señal de error actuante, es decir,
1 t
u (t ) = ∫ e(τ )dτ (11.14)
Ti 0
donde Ti es una constante denominada Tiempo Integral.
la ecuación (11.14) significa que la magnitud de la acción de control producida por este
controlador es proporcional a la magnitud del error, pero también a la duración de dicho error,
por lo tanto, conforme pasa el tiempo, si el error no desaparece esta acción de control aumenta
sin parar.
Ejemplo: Si consideramos el control integral de velocidad para el motor de CD, entonces podemos
obtener la Función de transferencia del sistema de lazo cerrado de la figura 11.9 , como sigue
1
G (s)
Y (s) Ti s
= (11.16)
R( s) 1 + 1 G ( s) H ( s)
Ti s
sustituyendo la función de transferencia del motor (11.5) y además H ( s ) = 1 , se obtiene
Y (s) 15.7
= 2
(11.17)
R( s ) 0.02Ti s + Ti s + 15.7
Usando el teorema del valor final obtenemos el valor de estado estable de la respuesta en el
tiempo
1570
y (∞) = lim sY ( s ) = = 100 (11.19)
s →0 0.02Ti (0)2 + Ti (0) + 15.7
por lo tanto, el error de estado estable ess será
ess = e(∞ ) = r (∞ ) − y (∞) (11.20)
100
50
0
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2
40
Va (t) (volts)
20
0
-20
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2
t (seg)
Figura 11.10.- Velocidad del motor de CD y acción de control Va (t ) con un controlador integral.
Esta acción de control produce una corrección proporcional a la rapidez con que crece el error
actuante. Es decir,
de(t )
u (t ) = Td (11.22)
dt
U (s)
= Td s (11.23)
E ( s)
Una manera de evitar el efecto indeseable de la acción puramente derivativa es combinarla con
una acción de control proporcional, para obtener la acción de control combinada Proporcional
Derivativa (PD), la cual producirá la acción de control como la combinación lineal de una parte
Laboratorio de Electrónica “Ing. Luis García Reyes”. UMSNH-FIE 8
Elaboró: José Juan Rincón Pasaye
Laboratorio de Control Analógico I Práctica No. 11 y 12.- Diseño de Controladores
proporcional al error actuante y otra parte proporcional a la rapidez del error actuante, es decir,
de(t )
u (t ) = K p 1 + Td (11.24)
dt
Es decir, su función de transferencia es
U (s)
= K p (1 + Td s ) (11.25)
E ( s)
Una acción de control Derivativa ideal no puede ser implementada, ni conviene hacerlo, ya que el
efecto que tendría en un sistema real que contiene señales con ruido de alta frecuencia por
pequeña que sea la amplitud de este ruido, su derivada podría ser muy grande y causar problemas
en el sistema. Por esta razón, los derivadores "prácticos" son aproximaciones a los derivadores
ideales. Por ejemplo, una aproximación muy usada es el llamado "derivador sucio" de primer
orden:
s
s≈ (11.26)
Ns + 1
donde N se elige de manera que el efecto del ruido no se amplifique mucho: Entre más pequeño
es N, más se acerca el derivador al ideal y más amplifica el ruido de alta frecuencia. Entre más
grande N más se atenúa el efecto del ruido de alta frecuencia, pero también se atenúan las señales
de rápida variación en el sistema lo cual puede disminuir las ventajas de la acción derivativa.
Ejemplo: Para poner a prueba la acción de control PD para la velocidad del motor de CD se puede
utilizar el diagrama de simulación mostrado en la figura 11.12.
s Scope
0.05
0.01s+1 Perturbacion
Td Derivador = -70rad/seg
15.7
0.5
0.02s+1
Ref=100r/s Kp Gm(s)
80
70
60
50
velocidad (rad/seg)
40
Referencia (rad/seg)
30
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2
100
Acción de control (volts)
va(t)
50
0
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2
t (seg)
Figura 11.13.- Velocidad del motor de CD y acción de control Va (t ) con controlador PD.
Se puede demostrar que para este ejemplo ess es exactamente igual para el controlador
PD que para el controlador puramente proporcional y por lo tanto no depende de Td .
Como pudo observarse en la figura 11.10. la acción de control Integral pura puede producir
oscilaciones muy pronunciadas en la respuesta del sistema, poor esta razón es recomendable
combinarla con otra acción como la proporcional. La acción de control Proporcional Integral (PI) es
una combinación lineal de una parte proporcional y una parte integral, es decir,
1 t
u (t ) = K p 1 + ∫ e(τ )dτ (11.27)
T1
0
Ejemplo: Para poner a prueba la acción de control PI para la velocidad del motor de CD se puede
utilizar el diagrama de simulación mostrado en la figura 11.13.
1 Scope
1/0.01
s
Perturbacion
1/Ti Integrator = -70rad/seg
15.7
0.5
0.02s+1
Ref=100r/s Kp Gm(s)
Figura 11.15.- Diagrama de simulación del controlador PI para la velocidad del motor de CD.
En la figura 11.15 se puede observar que para la simulación se han elegido los parámetros:
K p = 0.5 , Ti = 0.01 , con los cuales se obtiene la respuesta mostrada en la figura 11.16.
Se puede demostrar que para este ejemplo con la acción de controlador PD se logra
siempre ess = 0 siguiendo el procedimiento mostrado para el controlador proporcional
mediante el teorema del valor final.
Con un poco de experimentación en simulación es posible obtener una mejor elección de los
valores de los parámetros K p y Ti de manera que el esfuerzo de la acción de control se
disminuya.
100
80
ω (rad/seg)
60
40
velocidad (rad/seg)
20 referencia (rad/seg)
0
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2
va (t) (volts)
20
0
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2
t (seg)
Figura 11.16.- Velocidad del motor de CD y acción de control Va (t ) con controlador PI.
Es posible sacar ventaja de los tres tipos de acciones de control haciendo una combinación
adecuada de cada uno de ellos. La acción de control Proporcional-Integral-Derivativa (PID)
produce una corrección cuyo valor es la suma ponderada de las tres acciones básicas, la
proporcional, la integral y la derivativa, es decir,
1 t de(t )
u (t ) = K p e(t ) +
Ti ∫ e(τ )dτ + T
0
d
dt
(11.29)
Sintonización de PID's.
Aunque el controlador PID puede mejorar la respuesta del sistema en lazo cerrado, si no se eligen
adecuadamente los valores de sus parámetros K p , Ti y Td , puede empeorar el comportamiento
del sistema que se buscaba mejorar.
Ziegler y Nichols desarrollaron en 1942 un conjunto de reglas heurísticas para lograr una
sintonización que proporciona una buena atenuación de perturbaciones, pero permitiendo un
sobreimpulso no muy pequeño. Estas reglas fueron obtenidas buscando los parámetros que
proporcionaran un adecuado amortiguamiento de la oscilación dominante en lazo cerrado y un
valor pequeño de la integral del valor absoluto del error para una entrada escalón. Las reglas son
fundamentalmente empíricas más que teóricas y no garantizan un buen funcionamiento en todos
los casos.
0
t
L T
En este método se miden sobre la curva de reacción los parámetros R y L , donde R=K/T es el valor
de la pendiente de la recta tangente de mayor pendiente que sea posible trazar sobre la curva de
la respuesta al escalón del proceso, y L es el tiempo al cual la recta tangente cruza el eje del
tiempo, ver figura 5.7.
En el método de la ganancia límite, ver figura 11.19, se hace un experimento con la planta en lazo
cerrado con un controlador proporcional que lleva al sistema al límite de su estabilidad es decir,
que lo hace oscilar de manera sostenida, los dos parámetros que se utilizan son Ku y Tu, donde
Ku es la ganancia del controlador proporcional que logra la oscilación y Tu es el período de la
oscilación resultante.
escalón + y(t)
- Ku Planta
Tu
Figura 11.19. Método de ganancia límite
Las fórmulas de sintonización dadas por Ziegler y Nichols para un PID de la forma (11.30), se
muestran en la Tabla 11.1.
Ejercicio 1:
Si se desea controlar la velocidad del motor de CD de los ejemplos anteriores mediante un PID,
ninguno de los métodos de Ziegler y Nichols se puede utilizar para hacer la sintonización. Explica
Laboratorio de Electrónica “Ing. Luis García Reyes”. UMSNH-FIE 14
Elaboró: José Juan Rincón Pasaye
Laboratorio de Control Analógico I Práctica No. 11 y 12.- Diseño de Controladores
por qué. y utiliza las recomendaciones empíricas que se hicieron la sección anterior para obtener
una sintonización adecuada del PID para lograr los siguientes objetivos de control:
Ejercicio 2:
Perturbación
τL Motor de CD
Engranes
Vr(s) E(s) va Ke τm - Km ωm 1 θm Θ(s)
+ Controlador Ka + N
- Te s + 1 Tm s + 1 s
Actuador
Vs(s)
Ks
Sensor
Figura 11.20.- Control de posición de un servomecanismo mediante un PID.
2.1.- Obtener la respuesta al escalón unitario del sistema en lazo abierto (¿tiene la forma de "s"
del primer método de Ziegler y Nichols (figura 11.18) ?
2.2.- Probar un Controlador ON-OFF y sintonizar sus parámetros (UON, UOFF, EL, EH) por prueba y
error hasta lograr si es posible los objetivos de control.
2.3.- Probar un Controlador Proporcional y sintonizar su ganancia (Kp) por prueba y error hasta
lograr si es posible los objetivos de control.
2.4.- Probar un Controlador PD y sintonizar sus parámetros (Kp, Td) por prueba y error hasta lograr
si es posible los objetivos de control.
2.5.- Probar un Controlador PI y sintonizar sus parámetros (Kp, Ti) por prueba y error hasta lograr si
es posible los objetivos de control.
2.6.- Probar un controlador PID, Utilizar el método de Ziegler y Nichols para realizar la
sintonización inicial de los parámetros (Kp,Ti, Td). Después de haber sintonizado de esta manera los
parámetros del PID, modificarlos por prueba y error hasta lograr los objetivos de control.
2.7.- Decidir cuál de los controladores anteriores es el más adecuado para lograr los objetivos de
control.
Desarrollo de la Práctica.
1. Probar todos los ejemplos propuestos por el profesor conforme los va explicando.
2. Realizar todos los ejercicios propuestos.
Reportar:
1.- El diagrama de simulación, los parámetros obtenidos y las gráficas de respuesta de velocidad
ω (t ) (rad/seg) y de acción de control va (t ) (volts) para el Ejercicio 1.
2.- Los diagramas de simulación, los parámetros obtenidos y las gráficas de respuesta de posición
angular θ (t ) (rad) y de acción de control va (t ) (volts) para todos los incisos del ejercicio 2.