UNIVERSIDAD NACIONAL DE INGENIERA

FACULTAD DE INGENIERA MECNICA

AREA DE MECATRONICA

CURSO MT221 INGENIERIA DE CONTROL

ANLISIS DEL LUGAR GEOMTRICO DE LAS RAICES

2011
CRITERIO DE ROUTH

Es un mtodo que no nos da mucha informacin pero nos apoya en determinar el rango de
variacin de K para el cual el sistema de lazo cerrado es estable.
Sea p(s) = ansn + an-1sn-1 + an-2sn-2 +...........+ a1s + a0 = 0
Para aplicar este criterio se forma la tabla:

sn
sn-1
sn-2
sn-3
sn-4
.
.
.
s1
s0

an
an-1
b1
c1
d1

an-2
an-3
b2
c2
d2

an-4.......................ak
an-5.......................al
............... bi
0
............... cj
0
............... dk
0

0
0
0
0
0

e1
f1

e2
0

............... 0
..................0

0
0

0
0

Hallando bi, ci y di....

bi

b2

an
an1

Det

a n2
an3

( an an3 an1an2 )
an1

an 4
an5

( an an5 an1an4 )
an1

an1
an
an1

Det

an1

.
.
.
.
c1

d1

an1
b1

Det

b1

an3

b2
(an1b2 an3b1 )

b1

b1 b2

c1 c2 (b1c2 c1b2 )
c1
c1

Det

La informacin que vamos a utilizar se encuentra en la primera columna de la tabla, cada cambio
de signo en la primera columna indica una raz del polinomio en la mitad derecha del plano. Se
pueden aplicar dos reglas mas a p(s) antes de formar la tabla:
1.- Si cualquier coeficiente es negativo, entonces existen races con partes reales positivas.
2.- Si cualquier coeficiente, excepto a0, falta en el polinomio, entonces existen races con partes
reales positivas o bien races sobre el eje
imaginario.

Ejemplo 3.3.1
P(s) = s3 + 11s2 + 10s + K
La tabla de Routh es

s3

10

11

s1

b1

s0

c1

Calculando b1= -(k-110)/11

.
b 2= 0
.
b3 = K
Entonces la tabla ser

s3
1
10 0
2
s
11
K 0
s1
-(k-110)/11
0 0
0
s
K
0 0
Como las dos primeras entradas son 1 y 11 para que todas las races del polinomio
queden en la mitad izquierda del plano, las dos ultimas entradas deben ser tambin
positivas.. Por lo tanto:
-(k-110)/11>0 y K>0
Finalmente 0<K<110.
Ejemplo 3.3.2
P(s) = s3 + s2 - 4s 4 = (s+1)(s+2)(s-2)
s3
s2
s1
s0

La tabla de Routh es

1 -4
1 -4
b1 b2
c1 0

0
0
0
0

1 4

4 (4)
1 4
b1

0
1
1
1 0
Det

0 0)
1 0
b2

0
1
1
Det

Par continuar se debe reemplazar la fila de ceros por una de no ceros, se toma la fila
inmediata superior de la fila de ceros y se forma la ecuacin auxiliar s 2-4 al derivarla
obtenemos

d
( s 2 4) 2 s esto genera una nueva fila con un
ds

La nueva tabla es

s3
s2
s1

1 -4
1 -4
2 0

0
0
0

seguido de ceros

s0
1
2

Det

4
0

c1

0 ( 8)
4
2
2
La tabla es
s3
1 -4 0
s2
1 -4 0
s1
2 0
0
s0
-4 0 0
Hay un cambio de signo y por lo tanto una raz en la mitad derecha del plano.
c1

Ejemplo 3.3.2
P(s) = s4 + 3s3 + s2 +3s +2 = (s+1)(s+2)(s-2)
La tabla de Routh inicial
.

s4
s3
s2
s1
. s0

1 1
3
3
b1 b 2
c1 0
d1 0

2
0
0
0
0

0
0
0
0
0

1 1

3 3)
3 3
b1

0
3
3
1 2
Det

( 0 6)
3 0
b2

2
3
3
Como b2 0 y b1=0
Det

Reemplazamos b1 por un numero muy pequeo , donde es un numero positivo

arbitrariamente pequeo.
La tabla de Routh inicial
.

s4

s3

s2

s1

c1

d1

. s0

3 3

6 3 ) 3( 2)
2
c1

, entonces c1 es negativo

2
Det

(0 2c1 )
c1 0
d1

2
c1
c1
Existen 2 cambios de signo y por consiguiente dos polos en la mitad derecha del plano.
Det

PROBLEMAS PROPUESTOS
1. Grafique las races del polinomio caracterstico para las ganancias especificadas.
4

a. GH ( s ) s ( s 4) ,

K 0.1,0.51,4,20,50

K ( s 2)
, K 0.1,0.51,4,20,50
( s 4)( s 8)
K
c. GH ( s ) s ( s 2)(s 40) , K 1,5,50,100
K ( s 4)
, K 1,2,5,50,500.
d. GH ( s ) 2
s ( s 2)( s 10)

b. GH ( s )

2. Criterio de Routh. Verifique si algunas de las races queda en la mitad derecha del
plano
a. s4 + 7 s3 + 24 s2 + 58 s + 40
b. s4 + 23 s3 + 196 s2 + 624 s + 640
c. s6 + 5 s5 + 6 s4 + 10 s3 + 4 s + 1
d. s3 -2 s2 - 5 s + 6
3. Criterio de Routh :Rango de estabilidad. Determine e que rango de K el sistema de
lazo cerrado tiene polos estables, si GH es:
K ( s 1)

a. ( s 2)( s 1)
K

b. s ( s 1)( s 50)
c.

K ( s 1)
s 2 ( s 10)

6. INTRODUCCIN.- Veamos la ecuacin caracterstica

1 + G(s)H(s)

(6.1)

Esta ecuacin se analizara geomtricamente para obtener con rapidez la imagen completa de
cmo la respuesta de lazo cerrado variara con los cambio de ganancia.
Pero escribamos la ecuacin (6.1) de la siguiente manera:
GH(s) = -1 = 1-180
En coordenadas polares podemos escribir

GH(s)| GH(s) = 1 -180

De lo cual podemos escribir dos ecuaciones
|GH(s)| = 1 .....................(6.2)
GH(s) = -180............(6.3)
6.1.

REPRESENTACIN GRAFICA.-Sabemos que:

G ( s) H ( s)

g n ( s ) hn ( s ) K im1 ( s z i )

......( 6.4)
g d ( s ) hd ( s )
in1 ( s pi )

El factor K puede variar los polos de Tc(s).

Consideremos ahora un solo termino (s + ai)

Im(s)
s

+ai

s+ai

s+ai
-ai

ai

Re(s)

Representacin vectorial de s+ai

K im1 ( s zi )

in1 ( s pi )

G ( s) H ( s)
Entonces

K s zi s zi
s pi s pi

G ( s) H ( s )

K im1 s zi s zi
n
i 1 s pi s pi
K s zi
s zi s pi
s pi

Si igualamos al lado derecho que es 1 -180 tenemos las ecuaciones bases que son
las siguientes:
m

K s zi
i 1
n

s p

(6.5)

i 1

s z i

i 1

s p

180

(6.6)

i 1

Estas ecuaciones nos permiten realizar algunos anlisis.

MTODO GEOMTRICO LUGAR DE LAS RAICES.Veamos el siguiente sistema
.

.
.
.

K
( s 2)( s 4)

.
.
El ecuacin caracterstica ser 1

Tenemos:

k
0
( s 2)( s 4)

K
1 180
( s 2)( s 4)

Escribiendo en forma polar

K
1 180
( s 2 1 )( s 4 2 )

Tomamos como punto de prueba S1 (al lado derecho de s = -2)

K
K

0
( s1 2 0 )( s1 4 0 )
s1 2 s1 4

No

cumple con la condicin de ngulo pues 0 es diferente de -180.

Tomamos como punto de prueba S2 (al lado izquierdo de s = -4)
K
K
K

360
0
( s 2 2 180 )( s 2 4 180 )
s2 2 s2 4
s2 2 s 2 4

No cumple con la condicin de ngulo pues 0 es diferente de -180.

Tomamos como punto de prueba S3 ( -4< S3<-2)
K
K

180
( s3 2 180 )( s3 4 0 )
s3 2 s3 4

Por tanto

cualquier punto sobre el eje real entre los polos de GH es una solucin de la ecuacin
caracterstica.

Para el

sistema

Para
K

s3 2 s3 4
1
s3 2.05
2.05 2 2.05 4
1

0.0975 0.1

El mismo valor se obtiene para s = -3.95.

Tambin sobre la bisectriz perpendicular se cumple que 1+ 2 = 180.
La bisectriz pasa por 3.
El s = -3 es un punto de ruptura, punto en el cual el lugar geomtrico se aparta del eje real.

S4
S4+2

S4+4

S2
-5

S3

-4 S3+4 -3 S3+2

S1
-2

S1+4

-1
S1+2

S2+2

-1

S2+4
-2

GRAFICO DEL LUGAR DE LAS RAICES

En resumen, los pasos seguidos fueron:
1. Se escribe la ecuacin caracterstica de la forma GH = 1-180
2. Se grafican los polos y ceros reales de GH.
3. Se grafican los polos y ceros reales de GH con respecto a un punto s de solucin
prospectiva.
4. Si la suma de los ngulos de los ceros de GH menos la suma de los ngulos de los polos
de GH es igual a 180, entonces el punto es una solucin de la ecuacin caracterstica
con una ganancia especifica.
5. Se calcula la ganancia que coloca un polo de lazo cerrado en ese punto por medio de la
ecuacin respectiva.
6. K se calcula solamente cuando la condicin de ngulo ha sido satisfecha.
REGLAS DEL LUGAR GEOMTRICO DE LAS RAICES
Regla 1.- El lugar geomtrico de las races en el eje real se encuentra a la izquierda de una
cuenta impar de polos reales y ceros reales de GH.
Ejemplo 2.- Si GH

K ( s 1)
s ( s 4)( s 10)

Im(s)

Soluciones de GH + 1 = 0
.

X
-10

X
-4

Ejemplo 3.- Considere

GH

O
-1

Re(s)

K ( s 1) 2
s 2 ( s 10)

Soluciones de GH + 1 = 0
Im(s)

X
-10

O(2)
-1

X(2)

Regla 2.- El inicio del lugar geomtrico de las races es en los polos de GH.

Re(s)

Regla 3.- El final de las ramas en el lugar geomtrico de las races es en los ceros de GH.
n

Tenemos K

s p

cuando s p i

i 1
m

sz

cuando s z i

K 0
K

i 1

Ejemplo 4- GH

K ( s 1)( s 4)
s ( s 10)

Im(s)

Soluciones de GH + 1 = 0
X
-10

O
-4

Re(s)

-1

Regla 4.- El nmero de ramas que se extienden al infinito es igual a la diferencia entre el nmero
de polos y el nmero de ceros de GH. Este numero llamado exceso de polos sobre ceros, est
denotado por pex.
Las ramas que se extienden al infinito tienden a asntotas de lnea recta que se originan en un
punto comn. Los ngulos que las asntotas forman con el eje real se calculan con la formula.
(1 2l )(180 )
, l 0,1,... pex
pex
Regla 5.- La interseccin de las asntotas con el eje real esta dada por la formula

n
j

Re{ p i }

m
j

Re{ z i }

p ex

10

Ejemplo 1: Sea GH ( s )

K ( s 1)( s 5)
s ( s 6)( s 20)

El lugar geomtrico de las races se muestra en la figura siguiente:

Im(s)

Asntota en 180
Re(s)
-20

-6

-5

-1

La nica asntota esta en = -180. El polo de lazo cerrado que se origina a partir del polo de GH
en s = -20 viaja a lo largo de esta asntota hacia cero en s = - .
Regla 6.Los puntos en los cuales el lugar geomtrico de las races se aparta de y se acerca al eje real
pueden determinarse encontrando los puntos mximos y mnimos de la ganancia K en funcin de
s; con s restringida a valores reales.
Ejemplo 2.- Tomando el ejemplo ya analizado con
GH

K
( s 4)( s 2)

Vemos un punto de ruptura en s =-3. La figura muestra el lugar geomtrico de las raices como
una curva de K como una funcin de valores reales de s entre 2 y 4 con las ordenadas de las
dos curvas alineadas. El mximo ocurre en s = -3 con K = 1. En el punto donde K = 1, la ecuacin
caracterstica tiene una doble raz en s = -3. Esta es la ganancia mxima con la que los polos son
reales; las ganancias ms altas producen races complejas.

Im(s)

11

Re(s)
-4

-3

-2

Ganancia mxima es s=-3

-4
K=s2+6s+8

-2
dk/ds = 2s +6 = 0

mximo es s = -3

Ejemplo 3.- Consideremos el caso con:

GH

K ( s 1)
s ( s 4)( s 10)

En este caso, los polos de lazo cerrado comienzan en s = -10 y s = -4 y emigran uno hacia el
otro. Como en le caso del ejemplo habr una ganancia mxima en algn lugar a lo largo del eje
real entre estos dos puntos de inicio. La tabla proporciona la ganancia en la regin de s = -7,
aproximadamente en un punto intermedio entre los polos de GH en s igual a -10 y s =-4. Si GH
tuviera solo dos polos en s = -10 y s = -4, el alejamiento habra sido precisamente en s = -7 .
Sin embargo, el polo de GH en s = 0 y el cero en s = -1 hacen que el punto de alejamiento este
exactamente a la derecha de s = -7. Incluso as, el punto de alejamiento aun esta muy cerca de s
= -7 y el punto medido del intervalo constituye un excelente lugar para iniciar la bsqueda del
punto de alejamiento.
Tabla: Valores de ganancia a lo largo del eje real cerca del punto de alejamiento

12

s
K

-7.2
10.41

-7.1
10.46

-7.0
10.5

-6.9
10.51

Como pex = 2, habr dos asntotas en = 90. Las asntotas se cortan en

i = (-10-4-(-1))/2 = -6.5
Se muestra
un bosquejo de LGR.

-10

-4

-1

13

-6.8
10.50

-6.7
10.47

Ejemplo 4.Sea GH ( s)

K ( s 1)
2

s ( s 2)( s 20)

Hallar el L.G.R:

Solucin:
El primer paso es encontrar el L.G.R. sobre el eje real:

En el intervalo 2 < s < -1 ( la suma es 3, 2 polos mas 1 cero)

En el intervalo - < s < -10 ( 5, 4 polos mas 1 cero)

Pex = 4 1 = 3 que ser tambin el numero de asntotas

Las asntotas estarn en :
Considerando la formula l

(1 2l )(180 )
,
Pex

l 0,1,2..., Pex

(1 0)(180 )
60
3
(1 2)(180 )
2
180
3
(1 4)(180 )
3
300
3

14

(0 0 ( 2) ( 20)) ( 1) 21

7
4 1
3

En el origen se aparta 90.

Asntotas

-20

-7

-2

-1

Regla 7.- El lugar geomtrico de las races es simtrico con respecto al eje del plano s
Regla 8.- El ngulo de partida de una rama del lugar geomtrico de las races desde un polo de
GH, o el ngulo de arribo de una rama en un cero de GH puede determinarse satisfaciendo la
condicin del ngulo sobre un circulo de radio pequeo centrado en el polo o cero en cuestin.
Como este crculo se hace mas pequeo, todos los ngulos de los dems polos y ceros tienden a
valores exactos que son fciles de calcular. Solo el ngulo con respecto al polo o cero alrededor
del cual se construye este crculo de radio pequeo es entonces desconocido.
Ejercicio.- Considerando el sistema
GH ( s )

K ( s 2)
s ( s 1)( s 40)

15

Siguiendo los pasos anteriores delinear el LGR para que quede como el mostrado con matlab.

PROGRAMAS A USAR EN EL LABORATORIO

RL01
num=1;
den=[1 6 8];
rlocus(num,den);
axis([-10 0 -10 10]);
sgrid;
RL02
num=[1 1];
den=[1 14 40 0];
rlocus(num,den);
axis([-20 0 -100 100]);
sgrid;
RL03
num=[1 2 1];
den=[1 10 0 0];
rlocus(num,den);
axis([-15 0 -15 15]);
sgrid;
RL04
num=[1 1];
den=[1 14 40 0];
rlocus(num,den);
axis([-15 0 -15 15]);
sgrid;
for k=11:30
ec=[1 14 40+k k];
r=roots(ec)
k
end

16

RL05
num=[1 1];
den=[1 22 40 0 0];
rlocus(num,den);
axis([-16 16 -16 16]);
sgrid;

17