Lugar Geométrico de las Raíces en MATLAB
Lugar Geométrico de las Raíces en MATLAB
GL-SLS5201-L08M
LUGAR GEOMÉTRICO DE LAS RAÍCES
PROFESOR:
1. Introducción.
En ingeniería de control, el lugar de raíces o lugar de las raíces (del inglés, root locus) es el lugar geométrico de
los polos y ceros de una función de transferencia a medida que se varía la ganancia del sistema K en un
determinado intervalo. Esto es, el lugar geométrico de las raíces es la gráfica de la ubicación de las raíces (en el
plano s) del polinomio del denominador de la función de transferencia de lazo cerrado a medida se varía el
parámetro K (ganancia del sistema).
El método del lugar de raíces permite determinar la posición de los polos de la función de transferencia a lazo
cerrado para un determinado valor de ganancia K a partir de la función de transferencia a lazo abierto.
El lugar de raíces es una herramienta útil para analizar sistemas dinámicos lineales tipo SISO (single input
single output) y su estabilidad (BIBO stability). (Recuérdese que un sistema es estable si todos sus polos se
encuentran en el semiplano izquierdo del plano s (en el caso de sistemas continuos) o dentro del círculo unitario
del plano z (para sistemas discretos).)
En esta guía de laboratorio vamos a utilizar MATLAB y el Control System Toolbox para obtener la traza del
Lugar de las raíces y vamos a analizar los datos obtenidos a partir de esta gráfica.
2. Objetivos.
3. Duración.
4. Prerrequisitos
Ninguno
5. Bibliografía.
Escuela de Ingeniería
Página 1 de 34
GUÍA DE LABORATORIO
6. Marco teórico.
Sabemos que los polos de la función de transferencia se pueden desplegar gráficamente en el plano s. Vamos
a presentar un método analítico para desplegar la ubicación de los polos de la función de transferencia de lazo
cerrado (las raíces de la ecuación característica).
G
1 GH
Como una función del factor de ganancia K de la función de transferencia de lazo abierto. Este método llamado
Lugar de las raíces, solo requiere conocer la ubicación de los polos y ceros de la función de transferencia de
lazo abierto y no requiere la factorización del polinomio de la ecuación característica.
Escuela de Ingeniería
Página 2 de 34
GUÍA DE LABORATORIO
De esta forma las raíces de la ecuación característica 1 KGs H s 0 , dependen del parámetro factor de
ganancia K.
F s
s 1
ss 2
En s = -3+j4 , se obtiene;
F 3 j 4
3 j 4 1
3 j 4 3 j 4 2
20116,6
F 3 j 4
5126,9 17104,0
F 3 j 4
20
116,6 126,9 104,0
5 17
F 3 j 4 0,217 114,3
De esta forma una función de transferencia es una función compleja y se puede evaluar en un punto del plano
s, entregando un valor complejo que se puede escribir en forma polar, es decir, tiene un módulo y un ángulo.
1 KGs H s 0
O bien;
KGs H s 1
Escuela de Ingeniería
Página 3 de 34
GUÍA DE LABORATORIO
KG p1 H p1 1
O de otra forma;
N p1
K 1
D p1
N p1
K 1
D p1
A esta condición se le llama condición de magnitud. Por otro lado se debe cumplir que el ángulo de la función
compleja evaluada en el punto p1 como;
N p1
K 1 j 0
D p1
N p1 1 0
K tan 180 k 360
D p1 1
Donde k es un entero positivo. De otra forma el ángulo de fase debe ser múltiplo impar de 180°, esto es;
N p1
K 2k 1 180
D p1
Esta condición va a ser útil cuando se apliquen las reglas para la construcción del lugar de las raíces.
Específicamente la regla 10; Ángulo de despegue y ángulo de arribo.
Escuela de Ingeniería
Página 4 de 34
GUÍA DE LABORATORIO
Observe que las raíces de la ecuación característica dependen del factor de ganancia K.
s 2 10s K 0
Como la ecuación es de segundo grado tenemos dos soluciones (dos raíces) para cada valor de ganancia. Si
evaluamos para distintos valores de ganancia;
Observe que cuando K = 0 las raíces de la ecuación característica (polos de la función de transferencia de lazo
cerrado) son los polos de la función de transferencia de lazo abierto. Observe que desde K = 0 hasta K = 25, las
raíces son reales (se encuentran en el eje real). Recordamos que si las raíces de un sistema de segundo orden
son reales y diferentes la respuesta al escalón es exponencial sobreamortiguada. Si las raíces son reales e
idénticas la respuesta es críticamente amortiguada.
Escuela de Ingeniería
Página 5 de 34
GUÍA DE LABORATORIO
Observe que para todo valor de ganancia superior 25 las raíces son complejas conjugadas ( se encuentran en
el plano complejo) y por lo tanto la respuesta es subamortiguada.
Los polos de la función de transferencia de lazo abierto se muestran con una cruz negra. Y las raíces obtenidas
en la tabla anterior para cada valor de ganancia se muestran con una cruz azul.
Observe que cuando K = 0 (se inicia el lugar de las raíces) se tiene los polos de la función de transferencia de
lazo abierto. En otras palabras el lugar de las raíces se inicia en los polos de la función de transferencia de lazo
abierto. Observe además que para los polos complejos conjugados (en este caso) la parte real se mantiene
constante y aumenta la parte imaginaria.
Podemos extrapolar de inferir que a medida más aumenta la ganancia más aumentará la parte imaginaria de
los polos complejos conjugados. Finalmente unimos todos los polos para formar el lugar geométrico de las
raíces del sistema;
Escuela de Ingeniería
Página 6 de 34
GUÍA DE LABORATORIO
El punto donde el lugar geométrico de las raíces se desprende del eje real, es decir donde las raíces dejan de
ser reales y comienzan a ser complejas, se denomina punto de desprendimiento.
Existe un conjunto de reglas para la construcción de la traza del lugar de las raíces sin resolver el polinomio de
la ecuación característica para cada valor de ganancia y estas son las reglas que ahora vamos enunciar.
s 3 2s 2 3s K
Tiene tres ramas.
2. Los lugares geométricos de las raíces de un sistema con una ecuación característica real son
simétricos son simétricos respecto del eje real. Esto se debe a que las raíces complejas se presentan
en pares conjugados de la forma; s j
3. Los lugares geométricos de las raíces comienzan en los n polos donde K=0.
Escuela de Ingeniería
Página 7 de 34
GUÍA DE LABORATORIO
4. Los lugares geométricos de las raíces finalizan en los m ceros del sistema donde K = infinito. Si hay
más polos que ceros, que es el caso más común, entonces m lugares geométricos terminan en los m
ceros finitos y los (m-n) lugares geométricos restantes terminan en infinito.
5. Las porciones del eje real son secciones de los lugares geométricos de las raíces si el número de polos
y ceros a la derecha de dicha porción es impar.
6. Aquellos lugares geométricos que terminan en infinito tienden hacia las asíntotas que forman ángulos
respecto al eje real positivo de
,
3
,
5
, .....
2n m 1
nm nm nm nm
7. Las asíntotas intersectan sobre el eje real en un punto llamado centro de gravedad o centroide de las
asíntotas dado por,
Escuela de Ingeniería
Página 8 de 34
GUÍA DE LABORATORIO
Escuela de Ingeniería
Página 9 de 34
GUÍA DE LABORATORIO
La secuencia útil que se puede seguir, con la ayuda de las reglas anteriores, para construir el lugar geométrico
de las raíces se puede resumir como;
Escuela de Ingeniería
Página 10 de 34
GUÍA DE LABORATORIO
Ejemplo 1:
Dada el sistema de la figura, trazar el lugar de las raíces:
Solución:
Gs
K
s 1s 2s 3 K
Y la ecuación característica será;
s 1s 2s 3 K 0
2.- Cuando K=0 la ecuación característica se convierte en:
s 1s 2s 3 0
y los polos en lazo abierto están en -1, -2, y -3. No hay ceros
j plano s
-4 -3 -2 -1 1 2 3 4
3.- La ecuación es de tercer grado y, por lo tanto, se tendrán tres lugares geométricos de las raíces.
4.- Las porciones del lugar geométrico de las raíces sobre el eje real estarán entre -1 y -2 y desde -3 hasta
menos infinito (regla N°5)
j
plano s
-4 -3 -2 -1 1 2 3 4
Escuela de Ingeniería
Página 11 de 34
GUÍA DE LABORATORIO
5.- Tenemos tres ramas que parten en los tres polos de lazo abierto. No hay ceros en la función de
transferencia de lazo abierto, por lo tanto, el lugar de las raíces se inicia en estos tres polos y terminarán en tres
ceros que se encuentran en el infinito. Así, tenemos tres asíntotas; Los ángulos de las asíntotas, se determinan
mediante la ecuación dada en la regla 6;
,
3
,
5
, .....
2n m 1
nm nm nm nm
Donde n es el número de polos y m es el número de ceros. En este caso n = 3 y m = 0, es decir, n – m = 3. Y
por lo tanto, los ángulos de las tres asíntotas serán π/3 o 60°, 3 π/3 o 180° y 5 π/3 o 300° (o -60°).
6.- El punto de intersección de las asíntotas con el eje real, se determina usando la ecuación dada en la regla
número 7;
p
1 2 3 0 6 2
30 3
Podemos ahora graficar las asíntotas, ya que sabemos donde intersectan el eje real y cuales son sus ángulos;
(las asíntotas se muestran con líneas rojas segmentadas)
j
plano s
60.0°
180° -4 -3 -2 -1 1 2 3 4
p 2 -60.0°
j 1 j 2 j 3 K 0
j 3 2 6 j11 6 K 0
K 6 6 j
2 3
11 0 j 0
Escuela de Ingeniería
Página 12 de 34
GUÍA DE LABORATORIO
j
11 3.32
plano s
60.0°
180° -4 -3 -2 -1 1 2 3 4
p 2 -60.0°
s 1s 2s 3 K 0
K s 1s 2s 3
dK
3s 2 12s 11 0
ds
12 144 132
s 2 0.58
6
Solo es factible la solución s 2 0.58 ya que la otra queda fuera del lugar de las raíces sobre el eje real, es
decir, no queda dentro de una porción del lugar de las raíces sobre el eje real.
Escuela de Ingeniería
Página 13 de 34
GUÍA DE LABORATORIO
j
11 3.32
1.42
plano s
60.0°
180° -4 -3 -2 -1 1 2 3 4
p 2 -60.0°
Escuela de Ingeniería
Página 14 de 34
GUÍA DE LABORATORIO
j
11 3.32
1.42
plano s
60.0°
180° -4 -3 -2 -1 1 2 3 4
p 2 -60.0°
La gráfica del lugar de las raíces de este sistema usando MATLAB se obtiene mediante:
Zero/pole/gain:
1
-----------------
(s+1) (s+2) (s+3)
>> rlocus(G)
Escuela de Ingeniería
Página 15 de 34
GUÍA DE LABORATORIO
>> sisotool(G)
Escuela de Ingeniería
Página 16 de 34
GUÍA DE LABORATORIO
Seleccionando View puede desactivar y activar opciones de gráficos. Desactivamos Bode para apreciar
solamente el lugar de la raíces.
Escuela de Ingeniería
Página 17 de 34
GUÍA DE LABORATORIO
Y obtenemos:
Escuela de Ingeniería
Página 18 de 34
GUÍA DE LABORATORIO
Los cuadrados rosados sobre la gráfica indican la ubicación de los polos de lazo cerrado para la ganancia
actual del compensador (C(s) =1). Moviendo estos puntos se puede apreciar los valores de la ganancia de la
función de transferencia de lazo abierto.
Escuela de Ingeniería
Página 19 de 34
GUÍA DE LABORATORIO
Ejemplo 2:
Dada el sistema de la figura, trazar el lugar de las raíces:
K s 1
s 2 j3s 2 j3
K s 1
Gs
s 2 j3s 2 j3 K s 1
La ecuación característica será;
s 2 j3s 2 j3 K s 1 0
2.- Cuando K=0 la ecuación característica se convierte en:
s 2 j3s 2 j3 0
Y los polos en lazo abierto están en s 2 j3 y un cero en s 1.
j
3
plano s
2
1
-4 -3 -2 -1 1 2 3 4
-1
-2
-3
3.- La ecuación es de segundo grado y, por lo tanto, se tendrán dos lugares geométricos de las raíces.
Escuela de Ingeniería
Página 20 de 34
GUÍA DE LABORATORIO
4.- Las porciones del lugar geométrico de las raíces sobre el eje real estarán entre -1 y -2 y desde -3 hasta
menos infinito (regla N°5)
j
3
plano s
2
1
-4 -3 -2 -1 1 2 3 4
-1
-2
-3
5.- Puesto que n – m = 1, los ángulos de las asíntotas serán π/1 o 180°. La única asíntota se muestra con línea
roja segmentada.
j
3
plano s
2
1
180° -4 -3 -2 -1 1 2 3 4
-1
-2
-3
6.- Puesto que las asíntotas están sobre el eje real, el punto de intersección de las asíntotas con el eje real no
tiene interpretación.
7.- Las intersecciones con el eje imaginario se pueden determinar haciendo s j en la ecuación
característica, con lo que se obtiene.
j 2 j3 j 2 j3 K j 1 0
2 j 4 K 13 K 0
Escuela de Ingeniería
Página 21 de 34
GUÍA DE LABORATORIO
Al igualar las partes reales da por resultado 13 K y así; 4.1 . Debido a que la ganancia no puede
2
s 2 j3s 2 j3 K s 1 0
s 2 4 s 13 K s 1 0
s 2 4 s 13
K
s 1
dK
2 s 4 s 1 s 2 4 s 13
ds s 12
dK
2 s 2 6 s 4 s 2 4 s 13
ds s 12
dK
s 2 2s 9
ds s 12
dK s 2 2 s 9
2
ds s 2s 1
Al igualar esta ecuación a cero,
s 2 2s 9 0
Las raíces de esta ecuación son;
2 4 36
s 1 3.16
2
Y el solo el valor -4.16 es factible, el valor 2.16 no cae dentro de una porción del lugar de las raíces sobre el eje
real.
Escuela de Ingeniería
Página 22 de 34
GUÍA DE LABORATORIO
j
4.16 3
plano s
2
1
180° -4 -3 -2 -1 1 2 3 4
-1
-2
-3
9.- Se tiene un par de polos complejos y no hay ceros complejos. Los ángulos de partida del lugar geométrico
en K=0 de un polo complejo se determinan al usar la condición de ángulo y haciendo que s sea un punto sobre
el lugar geométrico de las raíces muy cercano al polo -2 + j3.
s j
θ2
3
plano s
2
θ1 1
-4 -3 -2 -1 1 2 3 4
-1
-2
θ3
-3
Puesto que hay un ángulo θ1, al cero que es 90° más tan 1 1 3 o 108,4°, un ángulo θ3, al otro polo de 90° y
un ángulo al polo de -2 + j3 de θ2, entonces
Escuela de Ingeniería
Página 23 de 34
GUÍA DE LABORATORIO
198.4°
j
4.16 3
plano s
2
1
180° -4 -3 -2 -1 1 2 3 4
-1
-2
161.6°
-3
10.- La figura muestra la gráfica completa de los lugares geométricos de las raíces
198.4°
j
4.16 3
1
180° -4 -3 -2 -1 1 2 3 4
-1
-2
161.6°
-3
Escuela de Ingeniería
Página 24 de 34
GUÍA DE LABORATORIO
1
Imag Axis
-1
-2
-3
-6 -5 -4 -3 -2 -1 0 1 2
Real Axis
Utilizando la herramienta SISO
>> sisotool(G)
Escuela de Ingeniería
Página 25 de 34
GUÍA DE LABORATORIO
Nuevamente podemos mover las raíces de la función de transferencia de lazo cerrado sobre la gráfica del lugar
geométrico de las raíces para obtener, por ejemplo, el punto de desprendimiento y la ganancia crítica de
estabilidad.
Escuela de Ingeniería
Página 26 de 34
GUÍA DE LABORATORIO
7. Actividades a realizar
Escuela de Ingeniería
Página 27 de 34
GUÍA DE LABORATORIO
a. Equipos requeridos
Un computador.
c. Instrumentos requeridos.
d. Herramientas requeridas.
Escuela de Ingeniería
Página 28 de 34
GUÍA DE LABORATORIO
e. Descripción y procedimiento.
Para completar esta actividad UD. deberá realizar con éxito los siguientes procedimientos:
Las siguientes instrucciones le solicitarán utilizar MATLAB y las funciones del control system toolboox para
obtener el lugar geométrico de las raíces para tres sistemas diferentes.
Paso Acción
1 SISTEMA 1
Para el sistema mostrado en la figura se pide obtener la traza del lugar geométrico de las raíces
utilizando el sisotool de MATLAB
U s 2 Y s
K
3s 1s 1
+
-
Controlador Planta
0.5
Sensor
Respuesta:
1. Obtener la ganancia para la cual el sistema tiene una respuesta (al escalón unitario)
críticamente amortiguada.
2. Obtener la ganancia para la cual el sistema es marginalmente estable.
Escuela de Ingeniería
Página 29 de 34
GUÍA DE LABORATORIO
2 SISTEMA 2
Para el sistema mostrado en la figura se pide obtener la traza del lugar geométrico de las raíces
utilizando el sisotool de MATLAB
U s 2 Y s
K
3s 1s 1
+
-
Controlador Planta
0.5
0.5s 1
Sensor
Respuesta:
1. Obtener la ganancia para la cual el sistema tiene una respuesta (al escalón unitario)
críticamente amortiguada.
2. Obtener la ganancia para la cual el sistema es marginalmente estable.
Escuela de Ingeniería
Página 30 de 34
GUÍA DE LABORATORIO
3 SISTEMA 3
Para el sistema mostrado en la figura se pide obtener la traza del lugar geométrico de las raíces
utilizando el sisotool de MATLAB
U s K 1 0.2s 2 Y s
+
-
3s 1s 1
Controlador
Planta
0.5
Sensor
Respuesta:
1. Obtener la ganancia para la cual el sistema tiene una respuesta (al escalón unitario)
críticamente amortiguada.
2. Obtener la ganancia para la cual el sistema es marginalmente estable.
Escuela de Ingeniería
Página 31 de 34
GUÍA DE LABORATORIO
ACTIVIDAD SOLICITADA
Paso Acción
1 Una planta tiene una función de transferencia de lazo abierto descrita por:
K s 2
G0 s
s 4s 13
2
a)
b)
c)
d)
e)
Escuela de Ingeniería
Página 32 de 34
GUÍA DE LABORATORIO
1. Cual es el comando del Control System ToolBox para desplegar el lugar de las raíces?
_______________________________________________________________________________
2. Cual es el comando del Control System Toolbox para llamara a la herramienta SISO?
_______________________________________________________________________________
_______________________________________________________________________________
_______________________________________________________________________________
Escuela de Ingeniería
Página 33 de 34
GUÍA DE LABORATORIO
Rut Nota
Alumno
Descripción Traza del lugar de las raíces para el análisis de sistemas lineales usando MATLAB
60% Habilidades
Nota % Descripción
Obtiene la gráfica del lugar de las raíces y los parámetros solicitados
Sistema 1 20% para el sistema 1.
Obtiene la gráfica del lugar de las raíces y los parámetros solicitados
Sistema 2 20% para el sistema 2.
Obtiene la gráfica del lugar de las raíces y los parámetros solicitados
Sistema 3 20% para el sistema 3.
Nota: N1 – Descuento.
Escuela de Ingeniería
Página 34 de 34