Repaso Grafico LGR y Sistemas de Fase No Minima
Repaso Grafico LGR y Sistemas de Fase No Minima
Repaso Grafico LGR y Sistemas de Fase No Minima
Análisis y diseño de sistemas de control por el método del lugar de las raíces 273
EJEMPLO 6-1 Considere el sistema de la Figura 6-3. (Se supone que el valor de la ganancia K es no negativo.)
Para este sistema,
K
G(s) % , H(s) % 1
s(s ! 1)(s ! 2)
Se dibuja la gráfica del lugar de las raíces y después se determina el valor de K tal que el factor de
amortiguamiento relativo f de los polos dominantes complejos conjugados en lazo abierto sea 0.5.
Para el sistema dado, la condición de ángulo es
K
G(s) %
s(s ! 1)(s ! 2)
%. s. s!1. s!2
% u 180o(2k ! 1) (k % 0, 1, 2, ...)
La condición de magnitud es
G G
K
!G(s)! % %1
s(s ! 1)(s ! 2)
Un procedimiento común para dibujar la gráfica del lugar de las raíces es el siguiente:
1. Determinar los lugares de las raíces sobre el eje real. El primer paso al construir una
gráfica del lugar de las raíces es situar los polos en lazo abierto, s % 0, s % .1 y s % .2, en el
plano complejo. (En este sistema no hay ceros en lazo abierto.) Las localizaciones de los polos en
lazo abierto se señalan mediante cruces. (En este libro las localizaciones de los ceros en lazo abier-
to se indicarán con círculos pequeños.) Observe que los puntos iniciales de los lugares de las raíces
(los puntos que corresponden a K % 0) son los polos en lazo abierto. Los lugares de raíces indivi-
duales para este sistema son tres, que coincide con el número de polos en lazo abierto.
Para determinar los lugares de las raíces sobre el eje real, se selecciona un punto de prueba, s.
Si el punto de prueba está en el eje real positivo, entonces
s % s ! 1 % s ! 2 % 0o
Esto demuestra que no es posible satisfacer la condición de ángulo. Por tanto, no hay un lugar de
las raíces sobre el eje real positivo. A continuación, se selecciona un punto de prueba sobre el eje
real negativo entre 0 y .1. Así,
s % 180o, s ! 1 % s ! 2 % 0o
Por tanto,
. s . s ! 1 . s ! 2 % .180o
y se satisface la condición de ángulo. Así, la parte del eje real negativo entre 0 y .1 forma parte
del lugar de las raíces. Si se selecciona un punto de prueba entre .1 y .2, entonces
s % s ! 1 % 180o, s ! 2 % 0o
y
. s . s ! 1 . s ! 2 % .360o
Se observa que no se satisface la condición de ángulo. Por tanto, el eje real negativo de .1 a .2
no forma parte del lugar de las raíces. Asimismo, si se sitúa un punto de prueba sobre el eje real
negativo de .2 a .ä, se satisface la condición de ángulo. Por tanto, existen lugares de las raíces
sobre el eje real negativo entre 0 y .1 y entre .2 y .ä.
2. Determinar las asíntotas de los lugares de las raíces. Las asíntotas de los lugares de las
raíces, conforme s tiende a infinito, se determinan del modo siguiente. Si se selecciona un punto de
prueba muy lejano al origen, entonces
K K
lím G(s) % lím % lím 3
sr ä sr ä s(s ! 1)(s ! 2) sr ä s
donde A(s) y B(s) no contienen K. Observe que f (s) % 0 tiene raíces múltiples en los puntos donde
d f (s)
%0
ds
Esto se observa del modo siguiente: suponga que f (s) tiene raíces múltiples de un orden r, donde
r n 2. En este caso, f (s) se escribe como
f (s) % (s . s1)r(s . s2) ñ (s . sn)
Si se diferencia esta ecuación con respecto a s y establecemos d f (s)/ds en s % s1, se obtiene
d f (s)
ds G s%s
1
%0 (6-7)
Esto significa que múltiples raíces de f (s) satisfarán la Ecuación (6-7). A partir de la Ecuación
(6-6) se obtiene
d f (s)
% Bñ(s) ! KAñ(s) % 0 (6-8)
ds
donde
dA(s) dB(s)
Añ(s) % , Bñ(s) %
ds ds
El valor específico de K que producirá raíces múltiples de la ecuación característica se obtiene de
la Ecuación (6-8) como
Bñ(s)
K%.
Añ(s)
Si se sustituye este valor de K en la Ecuación (6-6), se obtiene
Bñ(s)
f (s) % B(s) . A(s) % 0
Añ(s)
o bien
B(s)Añ(s) . Bñ(s)A(s) % 0 (6-9)
Si se despeja la Ecuación (6-9) para s, se obtienen los puntos en los que hay raíces múltiples. Por
otra parte, a partir de la Ecuación (6-6) se obtiene
B(s)
K%.
A(s)
y
dK Bñ(s)A(s) . B(s)Añ(s)
%.
ds A2(s)
Si dK/ds se hace igual a cero, se obtiene lo mismo que en la Ecuación (6-9). Por tanto, los puntos
de ruptura se determinan sencillamente a partir de las raíces de
dK
%0
ds
Debe señalarse que no todas las soluciones de la Ecuación (6-9) o de dK/ds % 0 corresponden a
los puntos de ruptura reales. Si un punto en el cual dK/ds % 0 está sobre el lugar de las raíces, se
trata de un punto de ruptura real o un punto de ingreso. De otro modo, si en un punto en el cual
dK/ds % 0, el valor de K tiene un valor positivo real, este punto es un punto de ruptura o un punto
de ingreso real.
Capítulo 6. Análisis y diseño de sistemas de control por el método del lugar de las raíces 277
Las frecuencias en los puntos de cruce con el eje imaginario son, por tanto, u % u ∂2. El valor
de ganancia que corresponde a los puntos de cruce es K % 6.
Una aproximación alternativa es suponer que s % ju en la ecuación característica, igualar con
cero tanto la parte imaginaria como la real y después despejar u y K. Para el sistema actual, la
ecuación característica, con s % ju, es
(ju)3 ! 3(ju)2 ! 2(ju) ! K % 0
o bien
(K . 3u2) ! j(2u . u3) % 0
Si se igualan a cero tanto la parte real como la imaginaria de esta última ecuación, se obtiene
K . 3u2 % 0, 2u . u3 % 0
278 Ingeniería de control moderna
de donde
u % u ∂2, K%6 o u % 0, K%0
Por tanto, los lugares de las raíces cruzan el eje imaginario en u % u ∂2, y el valor de K en los
puntos de cruce es 6. Asimismo, una rama del lugar de las raíces sobre el eje real tocará el eje
imaginario en u % 0.
5. Seleccionar un punto de prueba en una vecindad amplia del eje ju y el origen, como se
muestra en la Figura 6-5, y aplicar la condición de ángulo. Si un punto de prueba está sobre los
lugares de las raíces, la suma de los tres ángulos, h1 ! h2 ! h3, debe ser 180o. Si el punto de prueba
no satisface la condición de ángulo, seleccione otro hasta que se cumpla tal condición. (La suma de
los ángulos en el punto de prueba indicará en qué dirección debe moverse el punto de prueba.) Conti-
núe este proceso y sitúe una cantidad suficiente de puntos que satisfagan la condición de ángulo.
6. Dibujar los lugares de las raíces, tomando como base la información obtenida en los pasos
anteriores, tal y como se muestra en la Figura 6-6.
7. Determinar un par de polos dominantes complejos conjugados en lazo cerrado tales que el
factor de amortiguamiento relativo f sea 0.5. Los polos en lazo cerrado con f % 0.5 se encuentran
sobre las líneas que pasan por el origen y forman los ángulos ucos.1 f % ucos.1 0.5 % u60o
con el eje real negativo. A partir de la Figura 6-6, tales polos en lazo cerrado con f % 0.5 se obtie-
nen del modo siguiente:
s1 % .0.3337 ! j0.5780, s2 % .0.3337 . j0.5780
El valor de K que produce tales polos se encuentra a partir de la condición de magnitud, del modo
siguiente:
K % !s(s ! 1)(s ! 2)!s%.0.3337!j0.5780
% 1.0383
Usando este valor de K, el tercer polo se encuentra en s % .2.3326.
Observe que, a partir del paso 4, se aprecia que para K%6, los polos dominantes en lazo cerra-
do se encuentran sobre el eje imaginario en s % uj ∂2. Con este valor de K, el sistema mostrará
oscilaciones sostenidas. Para K b 6, los polos dominantes en lazo cerrado se encuentran en el se-
miplano derecho del plano s, produciendo un sistema inestable.
Por último, observe que, si es necesario, se establece con facilidad la parametrización de los
lugares de las raíces en términos de K mediante la condición de magnitud. Sencillamente se elige
un punto sobre un lugar de las raíces, se miden las magnitudes de las tres cantidades complejas s,
s ! 1 y s ! 2 y se multiplican estas magnitudes; el producto es igual al valor de la ganancia K en
dicho punto, o bien
!s! · !s ! 1! · !s ! 2! % K
La parametrización del lugar de las raíces se realiza de una manera sencilla utilizando MATLAB
(véase la Sección 6-3).
EJEMPLO 6-2 En este ejemplo se dibuja la gráfica del lugar de las raíces de un sistema con polos complejos
conjugados en lazo abierto. Considere el sistema de la Figura 6-7. Para este sistema,
K(s ! 2)
G(s) % 2 , H(s) % 1
s ! 2s ! 3
donde K n 0. Se observa que G(s) tiene un par de polos complejos conjugados en
s % .1 ! j ∂2, s % .1 . j ∂2
Un procedimiento común para dibujar la gráfica del lugar de las raíces es el siguiente:
1. Determinar los lugares de las raíces sobre el eje real. Para cualquier punto de prueba s
sobre el eje real, la suma de las contribuciones angulares de los polos complejos conjugados es de
360o, como se observa en la Figura 6-8. Por tanto, el efecto neto de los polos complejos conjuga-
dos es cero sobre el eje real. La localización del lugar de las raíces sobre el eje real se determina a
partir del cero en lazo abierto sobre el eje real negativo. Una prueba sencilla revela que una sec-
ción del eje real negativo, aquella que se encuentra entre .2 y .ä, es una parte del lugar de las
raíces. Se observa que, dado que este lugar geométrico se encuentra entre dos ceros (en s % .2 y
s % .ä), es en realidad parte de dos lugares de las raíces, cada uno de los cuales empieza en uno
Figura 6-8. Determinación del lugar de las raíces sobre el eje real.
de los dos polos complejos conjugados. En otras palabras, dos lugares de las raíces ingresan en la
parte del eje real negativo entre .2 y .ä.
Como existen dos polos en lazo abierto y un cero, hay una asíntota que coincide con el eje real
negativo.
2. Determinar el ángulo de salida de los polos complejos conjugados en lazo abierto. La pre-
sencia de un par de polos complejos conjugados en lazo abierto requiere la determinación del ángulo
de salida a partir de los mismos. El conocimiento de este ángulo es importante, debido a que el lugar
de las raíces cerca de un polo complejo proporciona información con respecto a si el lugar geométri-
co que se origina en el polo complejo emigra hacia el eje real o se tiende hacia la asíntota.
En la Figura 6-9, si se elige un punto de prueba y se mueve en la vecindad misma del polo
complejo en lazo abierto en s % .p1, ocurre que la suma de las contribuciones angulares del polo en
s % p2 y el cero en s % .z1 se considera sin alteración para el punto de prueba. Si el punto de prueba
está sobre el lugar de las raíces, la suma de hñ1, .h1 y .hñ2 debe ser u180o(2k ! 1), donde k % 0, 1,
2, .... Por tanto, en este ejemplo,
hñ1 . (h1 ! hñ2) % u180o(2k ! 1)
o bien
h1 % 180o . hñ2 ! hñ1 % 180o . h2 ! h1
El valor de la ganancia K en cualquier punto sobre el lugar de las raíces se encuentra aplicando
la condición de magnitud o utilizando MATLAB (véase la Sección 6-4). Por ejemplo, el valor de K
en el cual los polos complejos conjugados en lazo cerrado tienen el factor de amortiguamiento
relativo f % 0.7 se encuentra situando las raíces, como se muestra en la Figura 6-10, y calculando
el valor de K del modo siguiente:
(s ! 1 . j ∂2)(s ! 1 ! j ∂2)
K%
G s!2 G s%.1.67!j1.70
% 1.34
o bien
tan.1
A B
p!2
. tan.1
A p!1 B . tan.1
A p!1 B % u180o(2k ! 1)
u . ∂2 u ! ∂2 u
tan.1
A p!1 B ! tan.1
A p!1 B % tan.1
A Bp!2
u 180o(2k ! 1)
Estas dos ecuaciones corresponden a los lugares de las raíces del sistema actual. Observe que la
primera ecuación, u % 0, corresponde al eje real. El eje real desde s % .2 a s % .ä correspon-
de a un lugar de las raíces para K n 0. La parte restante del eje real corresponde a un lugar de las
raíces cuando K es negativo. (En el sistema actual, K es no negativo.) La segunda ecuación para el
lugar de las raíces es una ecuación de un círculo con centro en p % .2, u % 0 y radio igual a ∂3.
Esta parte del círculo a la izquierda de los polos complejos conjugados corresponde al lugar de las
raíces para K n 0. La parte restante del círculo corresponde al lugar de las raíces cuando K es
negativo.
Es importante observar que las ecuaciones que se interpretan con facilidad para el lugar de las
raíces sólo se obtienen para sistemas sencillos. No se recomienda intentar obtener las ecuaciones
para los lugares de las raíces en sistemas complicados que tengan muchos polos y ceros. Tales
ecuaciones son muy complicadas y su configuración en el plano complejo es difícil de visualizar.
Resumen de las reglas generales para construir los lugares de las raíces. Para
un sistema complejo en lazo abierto con muchos polos y ceros, puede parecer complicado cons-
truir una gráfica del lugar de las raíces, aunque en realidad no es difícil si se aplican las reglas
para construir dicho lugar. Situando los puntos y las asíntotas específicos y calculando los ángu-
los de salida de los polos complejos y los ángulos de llegada a los ceros complejos, se puede
construir la forma general de los lugares de las raíces sin dificultad.
A continuación se resumen las reglas y el procedimiento general para construir los lugares de
las raíces del sistema de la Figura 6-11.
Primero, obtenga la ecuación característica
1 ! G(s)H(s) % 0
A continuación, vuélvase a ordenar esta ecuación para que el parámetro de interés aparezca co-
mo el factor multiplicativo, en la forma
1. Situar los polos y ceros de G(s)H(s) en el plano s. Las ramas del lugar de las raíces
empiezan en los polos en lazo abierto y terminan en los ceros (ceros finitos o ceros en infinito).
A partir de la forma factorizada de la función de transferencia en lazo abierto, sitúense los polos
y los ceros en lazo abierto en el plano s. [Obsérvese que los ceros en lazo abierto son los de
G(s)H(s), mientras que los ceros en lazo cerrado son los de G(s) y los polos de H(s).]
Obsérvese que los lugares de las raíces son simétricos con respecto al eje real del plano s,
debido a que los polos y ceros complejos sólo aparecen en pares conjugados.
Una gráfica del lugar de las raíces tendrá tantas ramas como raíces tenga la ecuación caracte-
rística. Debido a que, por lo general, el número de polos en lazo abierto es mayor que el de ceros,
el número de ramas es igual al de los polos. Si el número de polos en lazo cerrado es igual al
número de polos en lazo abierto, el número de ramas individuales del lugar de las raíces que
terminan en los ceros finitos en lazo abierto será igual al número m de ceros en lazo abierto. Las
n . m ramas restantes terminan en infinito (n . m ceros implícitos en infinito) a lo largo de las
asíntotas.
Si se incluyen los polos y los ceros en infinito, el número de polos en lazo abierto es igual al
de ceros en lazo abierto. Por tanto, siempre se puede plantear que los lugares de las raíces empie-
zan en los polos de G(s)H(s) y terminan en los ceros de G(s)H(s) conforme K aumenta de cero a
infinito, donde los polos y los ceros incluyen los finitos y los infinitos en el plano s.
2. Determinar los lugares de las raíces sobre el eje real. Los lugares de las raíces sobre el
eje real se determinan a partir de los polos y los ceros en lazo abierto que se encuentran sobre él.
Los polos y los ceros complejos conjugados de la función de transferencia en lazo abierto no
afectan a la localización de los lugares de las raíces sobre el eje real, porque la contribución del
ángulo de un par de polos o ceros complejos conjugados es 360o sobre el eje real. Cada parte del
lugar de las raíces sobre el eje real se extiende sobre un rango de un polo o cero a otro polo o
cero. Al construir los lugares sobre el eje real, selecciónese un punto en éste. Si el número total
de polos y ceros reales a la derecha de este punto de prueba es impar, este punto se encuentra en
el lugar de las raíces. Si los polos y ceros en lazo abierto son simples, el lugar de las raíces y su
forma complementaria alternan segmentos a lo largo del eje real.
3. Determinar las asíntotas de los lugares de las raíces. Si el punto de prueba s se sitúa
lejos del origen, se considera que no cambia el ángulo de cada cantidad compleja. Entonces, un
cero en lazo abierto y un polo en lazo abierto cancelan los efectos del otro. Por tanto, los lugares
de las raíces para valores de s muy grandes deben ser asintóticos para líneas rectas cuyos ángulos
(pendientes) se obtengan mediante
u180o(2k ! 1)
Ángulos de las asíntotas % (k % 0, 1, 2, ...)
n.m
Si un punto de prueba se localiza muy lejos del origen, entonces, dividiendo el denominador
entre el numerador, se puede escribir G(s)H(s) como
K
G(s)H(s) %
s n.m
! [(p1 ! p2 ! ñ ! pn) . (z1 ! z2 ! ñ ! zm)]sn.m.1 ! ñ
o bien
K
G(s)H(s) % (6-12)
(p1 ! p2 ! ñ ! pn) . (z1 ! z2 ! ñ ! zm) n.m
C s!
n.m D
La abscisa de la intersección de las asíntotas y el eje real se obtiene igualando a cero el denomi-
nador del lado derecho de la Ecuación (6-12) y despejando s, o
Figura 6-12. Construcción del lugar de las raíces. [Ángulo de salida % 180o . [h1 ! h2] ! h.]
Capítulo 6. Análisis y diseño de sistemas de control por el método del lugar de las raíces 287
8. Determinar los polos en lazo cerrado. Un punto específico de cada ramificación del lu-
gar de las raíces será un polo en lazo cerrado si el valor de K en dicho punto satisface la condi-
ción de magnitud. Por otra parte, la condición de magnitud permite determinar el valor de la
ganancia K en cualquier localización de las raíces sobre el lugar. (Si es necesario, se establece
una parametrización de los lugares de las raíces en términos de K. Los lugares de las raíces son
continuos con K.)
El valor de K que corresponde a cualquier punto s sobre el lugar de las raíces se obtiene a
partir de la condición de magnitud, o bien
producto de las longitudes entre el punto s y los polos
K%
producto de las longitudes entre el punto s y los ceros
Este valor puede calcularse de forma gráfica o analíticamente. (Se puede utilizar MATLAB
para parametrizar el lugar de las raíces con K. Véase la Sección 6-3.)
Si en este problema se da la ganancia K de la función de transferencia en lazo abierto,
entonces, aplicando la condición de magnitud se encuentra la localización correcta de los polos
en lazo cerrado para un K determinado sobre cada ramificación de los lugares de las raíces,
mediante una aproximación de prueba y error o mediante MATLAB, que se presentará en la
Sección 6-3.
Comentarios acerca de las gráficas del lugar de las raíces. Se observa que la
ecuación característica del sistema cuya función de transferencia en lazo abierto es
Cancelación de los polos G (s) con los ceros de H (s). Es importante señalar que si
el denominador de G(s) y el numerador de H(s) contienen factores típicos, los polos y ceros en
lazo abierto correspondientes se cancelarán unos con otros, reduciendo el grado de la ecuación
característica en uno o más órdenes. Por ejemplo, considérese el sistema de la Figura 6-14(a).
(Este sistema tiene una realimentación de velocidad.) Si se modifica el diagrama de bloques de la
Figura 6-14(a) para obtener el de la Figura 6-14(b), se aprecia con claridad que G(s) y H(s) tie-
nen un factor común s ! 1. La función de transferencia en lazo cerrado C(s)/R(s) es
C(s) K
%
R(s) s(s ! 1)(s ! 2) ! K(s ! 1)
La ecuación característica es
[s(s ! 2) ! K](s ! 1) % 0
Sin embargo, debido a la cancelación de los términos (s ! 1) que aparecen en G(s) y H(s), se
tiene que
K(s ! 1)
1 ! G(s)H(s) % 1 !
s(s ! 1)(s ! 2)
s(s ! 2) ! K
%
s(s ! 2)
La ecuación característica reducida es
s(s ! 2) ! K % 0
La gráfica del lugar de las raíces de G(s)H(s) no muestra todas las raíces de la ecuación caracte-
rística; sólo las raíces de la ecuación reducida.
Para obtener el conjunto completo de polos en lazo cerrado, se debe agregar el polo cancela-
do de G(s)H(s) a aquellos polos en lazo cerrado obtenidos en la gráfica del lugar de las raíces de
G(s)H(s). No debe olvidarse que el polo cancelado de G(s)H(s) es un polo en lazo cerrado del
sistema, como se observa en la Figura 6-14(c).
Obsérvese que, una vez que se ha adquirido cierta experiencia con el método, resulta fácil
evaluar los cambios en los lugares de las raíces debidos a las modificaciones en el número y
situación de los polos y ceros en lazo abierto, visualizando las gráficas de los lugares de las
raíces que se producen de las diversas configuraciones de polos y ceros.
Resumen. A partir de los análisis anteriores, es evidente que se puede dibujar un diagra-
ma razonablemente preciso del lugar de las raíces para un sistema determinado, siguiendo reglas
sencillas. (Se sugiere al lector que estudie los diversos diagramas de los lugares de las raíces que
aparecen en los problemas resueltos al final del capítulo.) En las etapas de diseño preliminares,
no se necesitan las localizaciones precisas de los polos en lazo cerrado. Con frecuencia sólo se
necesitan sus localizaciones aproximadas para hacer una estimación de la representación del sis-
tema. Por tanto, es importante que el diseñador tenga la capacidad de dibujar con rapidez los
lugares de las raíces para un sistema determinado.
Gráfica de los lugares de las raíces con MATLAB. Al dibujar los lugares de las raí-
ces con MATLAB, se utiliza la ecuación del sistema obtenida por la Ecuación (6-11), que se
escribe como
num
1!K %0
den
donde num es el polinomio del numerador y den es el polinomio del denominador. Es decir,
Observe que ambos vectores, num y den, deben escribirse en potencias decrecientes de s. Una
orden de MATLAB que se usa con frecuencia para dibujar los lugares de las raíces es
rlocus(num,den)
Con esta orden, se dibuja en la pantalla la gráfica del lugar de las raíces. El vector de ganancias
K se determina de forma automática. (El vector K contiene todos los valores de ganancias para
los cuales se van a calcular los polos en lazo cerrado.)
Para los sistemas definidos en el espacio de estados, rlocus(A,B,C,D) dibuja el lugar de
las raíces del sistema con el vector de ganancias automáticamente determinado.
Obsérvese que las órdenes
rlocus(num,den,K) y rlocus(A,B,C,D,K)
Si se quiere dibujar los lugares de las raíces con las marcas 'o' o bien 'x', es necesario usar
la orden siguiente:
r % rlocus(num,den)
plot(r,'o') o plot(r,'x')
Es instructivo dibujar los lugares de las raíces mediante las marcas 'o' o bien 'x', debido a que
cada polo en lazo cerrado calculado se muestra de forma gráfica; en alguna parte de los lugares
de las raíces estas marcas se muestran de una forma densa y en otra parte aparecen separadas.
MATLAB produce su propio conjunto de valores de ganancias que se utilizan para obtener una
gráfica del lugar de las raíces. Lo consigue mediante una rutina interna de adaptación del tamaño
de paso. Asimismo, MATLAB usa la característica automática de fijar la escala del eje de la
orden plot.
EJEMPLO 6-3 Considere el sistema de control de la Figura 6-15. Dibuje el diagrama del lugar de las raíces con
una razón de aspecto cuadrada para que una línea con una pendiente de 1 sea una línea realmente
de 45o. Para dibujar el lugar de las raíces escoja la siguiente región:
.6 m x m 6, .6 m y m 6
donde x e y son las coordenadas del eje real y del eje imaginario, respectivamente.
Con el fin de establecer la región de la gráfica en pantalla para que sea cuadrada, introduzca la
orden
v % [.6 6 .6 6]; axis (v); axis('square')
Con esta orden, una línea con una pendiente de 1 estará realmente a 45o, y no inclinada por la
forma irregular de la pantalla.
Para este problema, el denominador se obtiene como un producto de términos de primer y se-
gundo orden. Por tanto, se deben multiplicar estos términos para obtener un polinomio en s. La
multiplicación de estos términos se realiza de una manera sencilla mediante la orden de convolu-
ción, tal y como se muestra a continuación.
Defina
a % s (s ! 1): a % [1 1 0]
b % s2 ! 4s ! 16: b % [1 4 16]
r % roots(b)
r%
–2.0000 ! 3.464li
–2.0000 – 3.464li
Por tanto, el sistema tiene los siguiente ceros y polos en lazo abierto:
Ceros en lazo abierto: s % .3
Polos en lazo abierto: s % 0, s % .1, s % .2 u j3.4641
El Programa MATLAB 6-1 dibujará el lugar de las raíces para este sistema. La gráfica aparece en
la Figura 6-16.
EJEMPLO 6-4 Considere el sistema cuya función de transferencia en lazo abierto G(s)H(s) es
K
G(s)H(s) % 2
s(s ! 0.5)(s ! 0.6s ! 10)
K
%
s ! 1.1s ! 10.3s2 ! 5s
4 3
No hay ceros en lazo abierto. Los polos en lazo abierto se localizan en s % .0.3 ! j3.1480,
s % .0.3 . j3.1480, s % .0.5 y s % 0.
Si se introduce en la computadora el Programa MATLAB 6-2, se obtiene la gráfica del lugar de
las raíces de la Figura 6-17.
Observe que, en las regiones cerca de x % .0.3, y % 2.3 y x % .0.3, y % .2.3, dos lugares
tienden uno al otro. Cabe preguntarse si estas dos ramas deben tocarse o no. Para explorar esta
situación, se pueden dibujar los lugares de las raíces utilizando pequeños incrementos de K en la
región crítica.
EJEMPLO 6-5 Considere el sistema de la Figura 6-19. Las ecuaciones del sistema son
x0 % Ax ! Bu
y % Cx ! Du
u%r.y
Capítulo 6. Análisis y diseño de sistemas de control por el método del lugar de las raíces 295
En este problema de ejemplo se obtiene el lugar de las raíces del sistema definido en el espacio de
estados. Por ejemplo, se considera el caso donde las matrices A, B, C y D, son
C D C D
0 1 0 0
A% 0 0 1 , B% 1 (6-15)
.160 .56 .14 .14
C % [1 0 0], D % [0]
El lugar de las raíces para este sistema se obtiene con MATLAB mediante la siguiente orden:
rlocus(A,B,C,D)
Esta orden producirá la misma gráfica del lugar de las raíces que se obtiene mediante la orden
rlocus(num,den) , donde num y den se obtienen de
[num,den] % ss2tf(A,B,C,D)
El Programa MATLAB 6-4 generará la gráfica del lugar de las raíces de la Figura 6-20.
Figura 6-20. Gráfica del lugar de las raíces del sistema definido en el espacio de estados,
donde A, B, C y D se obtienen de la Ecuación (6-15).
Lugares de las raíces con ! constante y "n constante. Recuérdese que en el plano
complejo la razón de amortiguamiento f de un par de polos complejos conjugados se puede ex-
presar en función del ángulo h, el cual se mide desde el eje real negativo, como se muestra en la
Figura 6-21(a), con
f % cos h
Con otras palabras, las líneas con razón de amortiguamiento f constante son radiales que pasan
por el origen como muestra la Figura 6-21(b). Por ejemplo, una razón de amortiguamiento de 0.5
requiere que los polos complejos se encuentren sobre una línea que pase por el origen con ángu-
los de u60o con el eje real negativo. (Si la parte real de un par de polos complejos es positiva, lo
que significa que el sistema es inestable, el correspondiente f es negativo.) La razón de amorti-
guamiento determina la localización angular de los polos, mientras que la distancia del polo al
origen se determina mediante la frecuencia natural no amortiguada un. El lugar de las raíces para
un constante son círculos.
Para dibujar las líneas con constante f y los círculos con constante un en el lugar de las raíces
con MATLAB se utiliza la orden sgrid.
Dibujo de las rejillas polares en el lugar de las raíces. La orden
sgrid
superpone líneas de razón de amortiguamiento constante (f % 0 V 1 con incremento 0.1) y
círculos de un constante en el dibujo del lugar de las raíces. Véase el Programa MATLAB 6-5
y el gráfico resultante en la Figura 6-22.
Si sólo se desean líneas para algún f constante en particular (por ejemplo, las líneas para
f % 0.5 y f % 0.707) y círculos para alguna un constante en particular (por ejemplo, los círculos
para un % 0.5, un % 1 y un % 2), se utiliza la siguiente orden:
sgrid([0.5, 0.707], [0.5, 1, 2])
Si se desea superponer líneas de f constante y círculos de un constante como los mencionados
anteriormente en un lugar de las raíces de un sistema con
num % [0 0 0 1]
den % [1 4 5 0]
Si se desea eliminar todas las líneas de f constante o todos los círculos de un constante se
utilizan los corchetes vacíos [ ] en los argumentos de la orden sgrid. Por ejemplo, si se quiere
únicamente superponer la línea de razón de amortiguamiento f % 0.5 y ningún círculo de un
constante en el lugar de las raíces que se muestra en la Figura 6-23, se utiliza la orden
sgrid(0.5, [])
Se puede ver del diagrama del lugar de las raíces de la Figura 6-25 que este sistema es solo
estable para rangos limitados del valor de K —que es 0 a K a 12 y 73 a K a 154. El sistema se
hace inestable para 12 a K a 73 y 154 a K. (Si K toma un valor que corresponde a operación
inestable, el sistema se puede deteriorar o hacerse no lineal debido a una no linealidad de satura-
ción que pueda existir.) Tal sistema se llama condicionalmente estable.
Figura 6-26 (a) Sistema de fase no mínima; (b) gráfica del lugar de las raíces.
Capítulo 6. Análisis y diseño de sistemas de control por el método del lugar de las raíces 301
K(1 . 0.5s)
Figura 6-27. Gráfica del lugar de las raíces de G(s) % .
s(s ! 1)
Para obtener el lugar de las raíces con MATLAB, introduzca el numerador y denominador
como siempre. Por ejemplo, si T % 1 seg y Ta % 0.5 seg introduzca el siguiente numerador y
denominador en el programa:
num % [-0.5 1]
den % [1 1 0]
El programa MATLAB 6-8 proporciona el lugar de las raíces que se muestra en la Figura 6-27.
Figura 6-28 Gráficas de los lugares de las raíces de ganancia constante y de fase constante
en el plano G(s)H(s).
real negativo del plano G(s)H(s), como se aprecia en la Figura 6-28. [Obsérvese que el plano
complejo empleado aquí no es el plano s, sino el plano G(s)H(s).]
Los lugares de las raíces y los lugares de ganancia constante en el plano s son mapeos de los
lugares de G(s)H(s) % u180o(2k ! 1) y de !G(s)H(s)! % una constante, en el plano G(s)H(s).
Debido a que los lugares de fase constante y de ganancia constante en el plano G(s)H(s) son
ortogonales, los lugares de las raíces y los lugares de ganancia constante en el plano s son orto-
gonales. La Figura 6-29(a) muestra los lugares de las raíces y los lugares de ganancia constante
para el sistema siguiente:
K(s ! 2)
G(s) % 2 , H(s) % 1
s ! 2s ! 3
Figura 6-29. Gráficas de los lugares de las raíces y los lugares de ganancia constante. (a) Sistema
con G(s) % K(s ! 2)/(s2 ! 2s ! 3), H(s) % 1; (b) sistema con G(s) % K/[s(s ! 1)(s ! 2)], H(s) % 1.
Capítulo 6. Análisis y diseño de sistemas de control por el método del lugar de las raíces 303
Obsérvese que, como la configuración de polos y ceros es simétrica con respecto al eje real, los
lugares de ganancia constante también son simétricos con respecto al mismo eje.
La Figura 6-29(b) muestra los lugares de las raíces y los lugares de ganancia constante para
el sistema:
K
G(s) % , H(s) % 1
s(s ! 1)(s ! 2)
Obsérvese que, como la configuración de los polos en el plano s es simétrica con respecto al eje
real y la línea paralela hacia el eje imaginario que pasa a través del punto (p % .1, u % 0), los
lugares de ganancia constante son simétricos con respecto a la línea u % 0 (eje real) y la línea
p % .1.
De las Figuras 6-29(a) y (b), obsérvese que cada punto en el plano s tiene su correspondiente
valor de K. Si se utiliza la orden rlocfind (presentada a continuación), MATLAB dará el valor
de K de un punto determinado y los polos en lazo cerrado más próximos correspondientes a ese
valor de K.
La orden rlocfind, que debe seguir a la orden rlocus, superpone unas coordenadas x-y móviles
sobre la pantalla. Mediante el ratón, se localiza el origen de las coordenadas x-y sobre el
punto deseado del lugar de las raíces, y se pulsa el botón del ratón. A continuación MATLAB
visualiza por pantalla las coordenadas de ese punto, el valor de la ganancia en ese punto y los
polos en lazo cerrado correspondientes a ese valor de la ganancia.
Si el punto seleccionado no se encuentra sobre el lugar de las raíces, tal como el punto A en
la Figura 6-29(a), la orden rlocfind devuelve las coordenadas del punto seleccionado, el valor de
la ganancia de ese punto, tal como K % 2, y las localizaciones de los polos en lazo cerrado, tales
como los puntos B y C correspondientes a ese valor de K. [Observe que cada punto en el plano s
tiene un valor de ganancia. Véase, por ejemplo, las Figuras 6-29(a) y (b).]
C(s) G(s)
%
R(s) 1 . G(s)H(s)
La ecuación característica es
1 . G(s)H(s) % 0 (6-17)
* Referencia W-4
304 Ingeniería de control moderna
Esta ecuación se despeja de forma parecida a como se hizo el desarrollo del método del lugar de
las raíces de la Sección 6-2. Sin embargo, debe cambiarse la condición de ángulo.
La Ecuación (6-17) se escribe como
G(s)H(s) % 1
que es equivalente a las dos ecuaciones siguientes:
La suma total de todos los ángulos a partir de los polos y ceros en lazo abierto debe ser igual a
0o u k360o. Por tanto, el lugar de las raíces ocupa un lugar de 0o, en contraste con el lugar de
180o que se consideró antes. La condición de magnitud no cambia.
Para ilustrar la gráfica del lugar de las raíces para el sistema con realimentación positiva, se
utilizarán como ejemplo las siguientes funciones de transferencia G(s) y H(s).
K(s ! 2)
G(s) % , H(s) % 1
(s ! 3)(s2 ! 2s ! 2)
uk360o
Ángulos de las asíntotas % (k % 0, 1, 2, ...)
n.m
La regla 5 se modifica del modo siguiente: cuando se calcula el ángulo de salida (o el ángulo de
llegada) a partir de un polo complejo en lazo abierto (o de un cero complejo), se deben restar de
0o la suma de todos los ángulos de los vectores que parten de todos los otros polos y ceros hacia
el polo complejo (o el cero complejo) en cuestión, incluyendo los signos adecuados.
Capítulo 6. Análisis y diseño de sistemas de control por el método del lugar de las raíces 305
Las otras reglas para construir la gráfica del lugar de las raíces no cambian. Ahora se aplica-
rán las reglas modificadas para construir la gráfica del lugar de las raíces.
1. Dibuje en el plano complejo los polos (s%.1 ! j, s%.1 . j, s%.3) y cero (s%.2)
en lazo abierto. A medida que K aumenta de 0 a ä, los polos en lazo cerrado empiezan
en los polos en lazo abierto y terminan en los ceros en lazo abierto (finitos o infinitos),
igual que en el caso de los sistemas con realimentación negativa.
2. Determine los lugares de las raíces sobre el eje real. Existen lugares de las raíces sobre el
eje real entre .2 y !ä y entre .3 y .ä.
3. Determine las asíntotas de los lugares de las raíces. Para el sistema actual,
uk360o
Ángulo de la asíntota % % u180o
3.1
Esto significa simplemente que las asíntotas están sobre el eje real.
4. Determine los puntos de ruptura y de ingreso. Como la ecuación característica es
(s ! 3)(s2 ! 2s ! 2) . K(s ! 2) % 0
se obtiene
(s ! 3)(s2 ! 2s ! 2)
K%
s!2
Diferenciando K con respecto a s, se obtiene
dK 2s3 ! 11s2 ! 20s ! 10
%
ds (s ! 2)2
Obsérvese que
Figura 6-31. Gráfica del lugar de las raíces para un sistema de realimentación
negativa con G(s) % K(s ! 2)/[(s ! 3)(s2 ! 2s ! 2)], H(s) % 1.
una raíz real se introduce en el semiplano derecho del plano s. Por tanto, para valores de K mayo-
res que 3, el sistema se vuelve inestable. (Para K b 3, el sistema debe estabilizarse con un lazo
externo.)
Obsérvese que la función de transferencia en lazo cerrado para el sistema con realimentación
positiva se obtiene mediante
C(s) G(s) K(s ! 2)
% % 2
R(s) 1 . G(s)H(s) (s ! 3)(s ! 2s ! 2) . K(s ! 2)
Para comparar esta gráfica del lugar de las raíces con la del sistema con realimentación nega-
tiva correspondiente, se muestran en la Figura 6-32 los lugares de las raíces para el sistema con
realimentación negativa cuya función de transferencia en lazo cerrado es
C(s) K(s ! 2)
% 2
R(s) (s ! 3)(s ! 2s ! 2) ! K(s ! 2)
La Tabla 6-2 muestra varias gráficas del lugar de las raíces de sistemas con realimentación
negativa y positiva. Las funciones de transferencia en lazo cerrado se obtienen mediante
C G
% , para sistemas con realimentación negativa
R 1 ! GH
C G
% , para sistemas con realimentación positiva
R 1 . GH
Figura 6-32. Gráfica del lugar de las raíces para un sistema de realimentación negativa con
G(s) % K(s ! 2)/[(s ! 3)(s2 ! 2s ! 2)], H(s) % 1.
Capítulo 6. Análisis y diseño de sistemas de control por el método del lugar de las raíces 307
donde GH es la función de transferencia en lazo abierto. En la Tabla 6-2, los lugares de las raíces
para los sistemas con realimentación negativa se dibujan con líneas y curvas gruesas y los de los
sistemas con realimentación positiva se dibujan con líneas y curvas discontinuas.
Las líneas y curvas gruesas corresponden a los sistemas con realimentación negati-
va; las líneas y curvas discontinuas corresponden a los sistemas con realimentación
positiva.