EE00302C

2.
Análisis temporal de sistemas continuos y discretos 53
2. Análisis temporal de sistemas continuos y discretos
2.1 Respuesta temporal.
La respuesta temporal de un sistema lineal invariante en el tiempo puede descomponerse en dos

partes: la respuesta transitoria y la respuesta estacionaria. De este modo, si denominamos c(t) a la
expresión de la respuesta temporal:
c( t ) = ct ( t ) + css( t )
donde: ct(t) = Respuesta transitoria.
css(t) = Respuesta estacionaria.
La respuesta transitoria es originada por la propia característica dinámica del sistema y determina el
comportamiento del sistema durante la transición de algún estado inicial hasta el estado final.
La respuesta estacionaria depende fundamentalmente de la señal de excitación al sistema y, si el

sistema es estable, es la respuesta que perdura cuando el tiempo crece infinitamente.
De este modo hemos logrado determinar de un modo simple la estabilidad absoluta de un sistema; se
dice que un sistema es estable si su respuesta transistoria decae a cero cuando el tiempo tiende a
infinito.
Se define el error en estado estacionario como la diferencia entre la señal de referencia y la señal
realimentada en estado estacionario en sistemas estables. Este error coincide con el valor estacionario
de la señal originada por el detector de error.
Por otra parte, en sistemas de control, interesa minimizar la desviación de la señal de salida respecto a
la señal de entrada en estado transitorio. Por esta razón se caracteriza la respuesta transitoria respecto
a entradas típicas o estándares, conociendo que, como el sistema es lineal, la respuesta del sistema a
señales más complejas es perfectamente predecible a partir del conocimiento de la respuesta a estas
entradas de prueba más simples.
Generalmente, las entradas típicas son: función impulsional, función escalón, función en rampa y
función parabólica en el tiempo; aunque la más importante de todas ellas es, sin duda, la función
escalón.
© los autores, 2000; © Edicions UPC, 2000.

54 Teoría de control. Diseño electrónico
2.1.1 Sistemas de primer orden.

Sistema de primer orden es aquel que únicamente posee un polo en su función de transferencia.
R(s) 1 C(s)
Ts + 1
Fig. 2.1 Sistema de primer orden. T: constante de tiempo del sistema
* Respuesta al escalón:
TL 1
r ( t ) = u( t ) 
→ R (s) =
s
 −
t 
1 1 TL− 1
C(s) = ⋅ → c( t ) = 1 − e T  ⋅ u( t )
Ts + 1 s  
Obsérvese que c(t) =1 cuando t tiende a infinito si T> 0; esto implica que el polo de la función de
transferencia del sistema debe encontrarse en el semiplano izquierdo del plano transformado S. Si
T ≤ 0, el sistema no alcanza el estado estacionario, resultando, de este modo, el sistema inestable.
Gráficamente:
c(t) −t
t
c( t ) = 1 − e T T
Pendiente=1/T c( t ) = 1 − e
100 %
98.2%
86.5%
63.2%
0 1T 2T 3T 4T 5T 6T
Fig. 2.2 Respuesta al escalón
Observando la gráfica se comprueba que para t = T la señal de salida ha alcanzado el 63.2 % del valor
final, siendo esta una medida típica en la caracterización de sistemas de primer orden.

2. Análisis temporal de sistemas continuos y discretos 55
A efectos prácticos, se considera que cuando han transcurrido por lo menos cuatro constantes de
tiempo, la señal de salida ha alcanzado el valor final.
t = 4T → c( t ) = 0.982 ⇒ css ( t ) ≈ 1
De este modo se deduce que cuanto más pequeña es la constante de tiempo de un sistema de primer
orden más rápidamente alcanza el valor final.
2.1.2 Sistemas de segundo orden.

Un sistema de segundo orden se caracteriza por poseer dos polos en su función de transferencia.
Expresión normalizada:
ωn 2
G (s) = 2
s + 2ξωns + ωn 2
donde se define:
ξ= Relación de amortiguamiento.
ωn= Frecuencia natural no amortiguada.
Los polos de un sistema de segundo orden vienen determinados por la expresión:
s1,2 = −ξωn ± ωn ξ 2 − 1
* Tipos de sistemas en función del valor de ξ:
En la siguiente gráfica se visualiza la ubicación de polos de la función de transferencia en función del

valor de ξ.
jω
ξ=0
ξ<1 jωn
ξ>1 ξ>1
ξ=1 -ωn σ
ξ<1 ξ=0
Fig. 2.3 Diagrama de la ubicación de los polos
Para ξ = 0 ⇒ Sistema oscilatorio. s1,2 = ± jωn

2
Para 0 < ξ < 1 ⇒ Sistema subamortiguado . s1,2 = −ξωn ± jωn 1 − ξ = −σ ± jωd
Para ξ = 1 ⇒ Sistema con amortiguamiento crítico . s1, 2 = −ωn
2
Para ξ > 1 ⇒ Sistema sobreamortiguado . s1,2 = −ξωn ± ωn ξ − 1
* Respuesta al escalón en un sistema subamortiguado (0<ξ<1):
ωn 2 1 TL− 1  ξ 
C(s) = 2 ⋅ → c( t ) = 1 − e− ξωnt cos(ωd ⋅ t ) + ⋅ sen(ωd ⋅ t )  ; t ≥ 0 (2.1)
s + 2ξωns + ωn s
2
 1− ξ 2

1   1 − ξ2  1
c( t ) = 1 − e − ξωnt ⋅ ⋅ sen ωd ⋅ t + tg −1    = 1 − e − ξωnt ⋅ ⋅ sen[ωd ⋅ t + θ] ; t ≥ 0
1 − ξ2   ξ  1 − ξ2
  
(2.2)
 
1 − ξ2 
donde:θ = tg −1  , con lo cual θ es el ángulo que forma el polo respecto al origen; de hecho
 ξ 
 
puede expresarse: ξ = cos(θ) . En conclusión, a medida que θ aumenta la relación de amortiguamiento
disminuye.
ζ=0
ζ=0.1
ζ=0.5
ζ=1
ζ=2
Wn t
Fig. 2.4 Respuesta al escalón de un sistema de segundo orden

La respuesta del sistema para ξ = 0 es oscilatoria, como puede comprobarse en la gráfica anterior. Por
otra parte, se observa que el valor máximo de amplitud logrado en el estado transitorio decrece a
medida que aumenta ξ (0 < ξ < 1) . En concreto, para valores de ξ≥1 desaparece totalmente este
máximo, lográndose el valor máximo de la señal cuando se alcanza el estado estacionario. Esta
característica es tanto más acentuada cuanto mayor es el valor de ξ . En concreto; si ξ>>1 la respuesta
del sistema de segundo orden puede aproximarse por la respuesta de un primer orden, realizándose lo
que se conoce como aproximación de polo dominante. La respuesta, en este caso, queda caracterizada
por la raíz real de constante de tiempo mayor (la más cercana al eje imaginario en el plano S). Debe
observarse que esta aproximación no puede realizarse cuando se analiza el inicio de respuesta, donde
todavía perdura el efecto de la señal originada por la raíz no dominante.
2.2 Especificaciones de respue sta transitoria.
Generalmente, en la práctica, se especifican las características o especificaciones requeridas a un

sistema de control en términos de cantidades en el dominio del tiempo. Estas cantidades vienen
determinadas en términos de la respuesta transitoria frente a una entrada tipo escalón (u(t)). A esta
respuesta se le denomina respuesta indicial.
De este modo se caracteriza la dinámica de un sistema, aunque trabaje con otro tipo de entradas, a
través de la dinámica requerida frente a una entrada escalón. El significado de los parámetros
definidos en la caracterización de la respuesta indicial determina la forma de la respuesta transitoria de
un sistema.
En la figura siguiente puede observarse la respuesta típica de un sistema frente a una entrada escalón.
M p
T o le r a n c ia ( 5 % o 2% )
1 .0
0 .5
td
tr
tp
ts
Fig. 2.5 Respuesta de un sistema de segundo orden frente a una entrada en escalón
Definiciones:
* td: Tiempo de retardo.
Tiempo que tarda la respuesta en alcanzar el 50 % del valor final.

* tr: Tiempo de subida o crecimiento.
Sistemas subamortiguados: Tiempo que tarda la respuesta en alcanzar el 100 % del valor final.
Sistemas sobreamortiguados: Tiempo que tarda la respuesta en pasar del 10 % al 90 % del valor final.
* tp: Tiempo de pico.
Tiempo que tarda la respuesta en alcanzar el primer máximo. Solo existe en sistemas
subamortiguados.
* Mp: Máximo sobreimpulso (porcentual).
Es el valor de pico máximo de la curva de respuesta ponderado por el valor final obtenido.
c( tp) − c(∞)
Mp(%) = ⋅100 (%) (2.3)
c(∞)
El máximo sobreimpulso es un parámetro indicativo de la 'estabilidad relativa' del sistema.
* ts: Tiempo de establecimiento.
Tiempo requerido por la curva de respuesta para alcanzar y mantenerse dentro de determinado rango
(habitualmente 5 % o 2 %) del valor final. Generalmente está relacionado con las constantes de
tiempo más grandes del sistema.
Idealmente interesará lograr siempre sistemas de control que minimicen el máximo sobreimpulso y los
tiempos de respuesta transitoria, manteniendo la máxima precisión posible.
2.2.1 Particularización para sistemas de segundo orden subamortiguados.
1- Cálculo de tr:
  2
− ξωntr  ξ  1− ξ ωd
c( tr ) = 1 ⇒ 1 = 1 − e cos(ωd ⋅ tr ) + ( )
⋅ sen(ωd ⋅ tr ) ⇒ tg ωd ⋅ tr = − =− (2.4)
 2  ξ σ
 1− ξ 
1  ωd  π − θ  ωd 
tr = ⋅ arctg − = , donde θ = arctg  , con lo cual θ es el ángulo que forma el polo
ωd  σ ωd  σ
respecto al origen. En el caso de poseer raíces imaginarias puras (ξ=0), el tiempo de subida coincide
con un cuarto del periodo de señal (tr = π 2ωn ); a su vez, cuando las raíces son reales (ξ=1), el sistema
tiene una característica de amortiguamiento crítico y el tiempo de subida, según la expresión anterior,
es infinito, lo que conlleva la necesidad de variar su definición.

2- Cálculo de tp:
− ξωntp
dc( t ) ωn ⋅ e
=0⇒ ⋅ sen(ωd ⋅ tp ) = 0 ⇒ ωd ⋅ tp = k ⋅ π (2.5)
dt 2
tp 1− ξ
π Td
El primer máximo se obtiene para: tp = = , donde Td es el periodo de oscilación.
ωd 2
De este modo, a medida que aumente la parte imaginaria de los polos del sistema, el tiempo de pico
disminuirá.
3- Cálculo del máximo sobreimpulso:
 
−ξωntp  ξ − ξωnπ
Mp = c( tp ) − 1 ⇒ 1 − e cos(ωd ⋅ tp ) + ⋅ sen(ωd ⋅ tp )  = e ωd
 2 
 1− ξ 
− σπ − πξ
Mp = e ωd = e 1− ξ 2
(2.6)
Observando la expresión anterior puede trazarse una gráfica que relacione Mp en función de ξ.
%
100
−πξ
90 Mp = e 1−ξ2
80
70
60
Mp
50
40
30
20
10
0
0 0.2 0.4 0.6 0.8 1
ξ
Fig. 2.6 Curva de Mp en función de ξ
De hecho, cuanto mayor es ξ, menor es el máximo sobreimpulso obtenido en la respuesta transitoria.

Expresado de otro modo, cuanto mayor es el ángulo que forma el polo del sistema respecto al eje real
mayor es el sobreimpulso de la respuesta. En conclusión, el aumento de la parte real (en valor
absoluto) o la disminución de la parte imaginaria de las raíces conlleva la reducción del Mp. De la
figura 2.6 se puede observar que para ξ>0.6 el Mp es menor al 10%.

4- Cálculo del tiempo de establecimiento:
Para determinar las expresiones del tiempo de establecimiento en un sistema subamortiguado, se

considera que la respuesta queda comprendida entre envolventes, que permiten aproximar el tiempo
de establecimiento en función de la parte real de los polos del sistema.
Así, considerando los errores:
3 3
5 %⇒ ts ≅ = (2.7)
ξωn σ
4 4
2 %⇒ ts ≅ = (2.8)
ξωn σ
el tiempo de establecimiento es inversamente proporcional a la parte real (en valor absoluto) de las
raíces del sistema subamortiguado.
2.5
2 e − ξω t
1+
1 − ξ2
1.5
c(t)
1
0.5 e − ξωnt
1−
1− ξ2
0
-0.5
0 1 2 3 4T 4 5
Tiempo(seg)
Fig 2.7 Envolventes de respuesta al escalón unitario
En la figura 2.8 se presentan gráficas de las respuestas temporales al escalón unitario para tres
sistemas de segundo orden con diferentes posiciones de polos en el plano S. En estas figuras se
pueden verificar, mediante las expresiones 2.4, 2.5, 2.6 y 2.8, las relaciones entre las posiciones de los
polos y las características temporales de tiempo de subida, tiempo de pico, máximo sobreimpulso y
tiempo de establecimiento.

10 2
P3
8 P3 1.8
P2
6 1.6
P2 P1
4 1.4
2 1.2
Eje Imaginario
Amplitud
0 1
-2 0.8
P1*
-4 * 0.6
P2
-6 0.4
-8 * 0.2
P3
-10 0
-0.5 -0.4 -0.3 -0.2 -0.1 0 0 1 2 3 4 5 6
Eje Real Tiempo (seg)
2 2
P3 P1
1.8 1.8 P1
1.6 1.6
1.4 1.4
1.2 1.2 P2
Amplitud
1 1
Amplitud
0.8 0.8
0.6 0.6
0.4 0.4
0.2 0.2
0 0
0 1 2 3 4 5 6 0 1 2 3 4 5 6
Tiempo (seg) Tiempo (seg)
Fig 2.8 Respuestas al escalón de sistemas de segundo orden
Ejemplo 2.1
En la figura 2.9. se muestra un amplificador magnético con una baja impedancia de salida, en cascada
con un filtro paso-bajo y un preamplificador. El amplificador tiene una elevada impedancia de entrada
así como ganancia unitaria, y se usa para sumar las señales tal y como se muestra en la figura.
R=50Ω
-1 E(s) K V0(s)
V1(s)
Vin(s) 1 + τs
+1 C
Fig. 2.9 Amplificador magnético en cascada con un filtro paso-bajo y preamplificador
Calcular el valor de la capacidad C, de tal forma que la función de transferencia V0(s)/Vin(s) tenga una
relación de amortiguamiento ξ = 0.7 .

Si la constante de tiempo del amplificador magnético es igual a 1 seg. (τ = 1 seg.) y la ganancia es

K=20. Calcular el tiempo de estabilización del sistema resultante.
4
Nota: Tiempo de estabilización ⇒ ts =
ξω n
Solución:
V0 (s)
1 
= Cs  1
V1 (s) R + 1  V0 (s) K Cs = K ⋅ 1
Cs  = ⋅
V1 (s) K  E(s) 1 + τs R + 1 1 + τs RCs + 1
= Cs
E(s) 1 + τs 

Obteniendo, entonces, el diagrama de bloques siguiente :
Vin(s) E(s) K V1(s) 1 V0(s)

+
- 1 + τs RCs + 1
V0 (s) 20 1
G (s) = = ⋅
E(s) (1 + s) (50Cs + 1)
V0 (s) G (s) 20 20
G LC (s) = = = =
Vin (s) 1 + G (s) (1 + s) ⋅ (50Cs + 1) + 20 50Cs 2 + (1 + 50C)s + 21
20
G LC (s) = 50C
2 1 + 50C 21
s + s+
50C 50C
Para calcular el tiempo de estabilización, es necesario conocer ξ y ωn, ya que ts = 4/ξωn. Si ξ debe
valer 0.7:
1 + 50C 
= 2 ξω n = 2 ⋅ 0.7 ⋅ ω n  1 + 50C 21
50C = 1.4 ⋅
21 
= ω n2  50C 50C
50C 
 0.7826 F
1 + 100C + 2500C 2 = 1.96 ⋅ 21 ⋅ 50C ; 2500C 2 − 1958C + 1 = 0 → C=
. ⋅ 10 −4 F
511

Desechamos el valor de C=0.7826 F por ser un valor físicamente muy grande, además de que provoca
un tiempo de establecimiento mayor.
Escogeremos por lo tanto: C = 5.11·10-4 F
21 21
ωn = = = 28.7 rad / s
50C . ⋅ 10 −4
50 ⋅ 511
El tiempo de estabilización quedará entonces como:
4 4
ts = = = 0.2 seg.
ξω n 0.7 ⋅ 28.7
2.2.2 Sistemas de orden superior.
Un sistema de orden superior (supuesto estable) puede caracterizarse mediante una función de
transferencia, que admite una expresión de la respuesta a una entrada escalón unitario de la forma:
m m
k ⋅ ∏ (s + z i ) k ⋅ ∏ (s + z i )
C(s) 1
= q r
i =1
→ C(s) = q r
i =1
⋅
R (s) s
∏ (s + p ) ⋅ ∏ ( s
j =1
j
k =1
2
+ 2ξ k ω k s + ω 2k ) ∏ (s + p ) ⋅ ∏ ( s
j =1
j
k =1
2
+ 2ξ k ω k s + ω 2k )
Desarrollando en fracciones parciales y aplicando la antitransformada de Laplace:
[ ]
q r
c( t ) = a ⋅ u ( t ) + ∑ a j ⋅ e + ∑ b k ⋅ e − ξ k ω k t cos(ωd k ⋅ t ) + c k ⋅ e − ξ k ω k t sen(ωd k ⋅ t )
−p j t
j =1 k =1
En la expresión aparecen términos que dependen de la parte real de los polos. Dos consideraciones
permiten aproximar la dinámica de un sistema de orden superior por la dinámica de un sistema de
primer o segundo orden.
Cancelación cero-polo: Se basa en el análisis resultante en la determinación del valor de los residuos
cuando el polo se encuentra cercano a un cero; en este caso, puede expresarse zi = pi + ε , donde
ε → 0. El cálculo del residuo resulta:

m
k⋅ ∏ (s + z )
h =1
h
ai = lim q
⋅ (s + pi)
s → − pi r
s⋅ ∏
j =1
(s + p j ) ⋅ ∏ k =1
(s 2 + 2ξ k ωk s + ω2k )
m −1
k⋅ ∏ (s + z )
h =1
h
ai = ε ⋅ lim q −1 r
s → − pi
s⋅ ∏ (s + p ) ⋅ ∏ (s
j =1
j
k =1
2
+ 2ξ k ωk s + ω2k )
en la segunda expresión no existen los términos (s+zi) en numerador ni (s+pi) en denominador. En

conclusión, en el caso en el cual el polo y el cero se encuentren muy cercanos, ε → 0 ⇒ ai → 0 ,
perdiendo importancia el término exponencial de la respuesta transitoria originado por este polo frente
al resto de términos de la respuesta total. De este modo puede simplificarse en un orden la dinámica
del sistema. Debe considerarse que se obtienen resultados análogos en el caso de poseer ceros
complejo-conjugados cercanos a polos complejo-conjugados, lo cual permite generalizar esta
aproximación a cualquier tipo de sistemas de orden superior.
Aproximación de polos dominantes: Los polos más alejados del eje imaginario poseen constantes de
tiempo menores (sean o no reales, en el caso de raíces complejo-conjugadas existe una dependencia
respecto a su parte real), de manera que puede afirmarse que las exponenciales debidas a estos polos
son importantes en el inicio de la respuesta transitoria, pero que decaen a cero mucho más
rápidamente que las exponenciales debidas a raíces con constantes de tiempo mayores. Son estas
últimas las que caracterizan plenamente la respuesta transitoria (exceptuando en el origen de la
respuesta) y permiten reducir el orden del sistema; se dice en este caso que “dominan” la respuesta del
sistema, despreciándose el efecto de las raíces con parte real mayor (en valor absoluto).
2.3 Respuesta transitoria de si stemas discretos.
2.3.1 Sistema de control discreto en lazo cerrado.

Estructura típica de un sistema de control discreto en lazo cerrado:
CONTROL
DISCRETO ZOH PLANTA
R(s) + E(s) E*(s) C(s)
Gc(z) Goh(s) G(s)
- T T
H(s)
ELEMENTO DE MEDIDA
Fig 2.10 Estructura típica de un sistema de control discreto en lazo cerrado
El control discreto es caracterizado a través de su función de transferencia en Z y es diseñado para

cumplir las especificaciones requeridas al sistema.

Se define el sistema continuo equivalente de un sistema discreto como aquel sistema que resulta de
eliminar los muestreadores y mantenedores de datos del sistema discreto.
Las respuestas temporales del sistema discreto y del sistema continuo equivalente difieren
substancialmente a medida que aumenta el periodo de muestreo, de este modo se dice que el muestreo
desvirtúa la respuesta del sistema discreto frente a la del sistema continuo equivalente. En concreto,
el muestreo tiene un efecto desestabilizador del sistema, de manera que, cuanto más desvirtuado se
halla el sistema discreto frente al sistema analógico, “peor” es su respuesta transitoria. Este efecto
conlleva una pérdida de la estabilidad relativa del sistema discreto a medida que aumenta el periodo
de muestreo.
Ejemplo 2.2
Dado los siguientes sistemas:
ZOH PLANTA PLANTA

R(s) + E*(s) C(s) R(s) + C(s)
Goh(s) G(s) G(s)
- T -
H(s) H(s)
ELEMENTO DE MEDIDA ELEMENTO DE MEDIDA
a) b)
Fig. 2.11 a) Sistema discreto
b) Sistema continuo equivalente
Con las siguientes funciones de transferencias:
1
G (s) = ; H ( s) = 1
s( s + 1)
Se puede observar el efecto desestabilizador al aumentar el periodo de muestreo en las siguientes

gráficas:

Amplitud Amplitud
1.2 1.4
1.2
1
1
0.8
0.8
0.6
0.6
0.4
0.4
0.2 0.2
0 0
0 2 4 6 8 10 0 5 10 15 20
Tiempo (seg) Nº de muestras
a) b)
Amplitud Amplitud
1.4 1.2
1.2 1
1
0.8
0.8
0.6
0.6
0.4
0.4
0.2
0.2
0 0
0 10 20 30 40 50 0 20 40 60 80 100
Nº de muestras Nº de muestras
E d)
Fig. 2.12 a) Respuesta al escalón del sistema equivalente

b) Respuesta al escalón del sistema discreto con T=1seg.
c) Respuesta al escalón del sistema discreto con T=0.25seg.
d) Respuesta al escalón del sistema discreto con T=0.1seg.
Los parámetros de medida de este proceso son, fundamentalmente, dos:
* Número de muestras por constante de tiempo:
Corresponde a la expresión: τ/T . Se considera que si se poseen más de 5 muestras por constante de
tiempo de la respuesta del sistema discreto no queda desvirtuada frente a la respuesta del sistema
continuo equivalente. Obsérvese que τ es la constante de tiempo más pequeña del sistema continuo
equivalente.
* Número de muestras por ciclo:
Corresponde a la expresión: Td/T. En este caso no existe ningún criterio al respecto, pero se considera
que son suficientes 10 muestras por ciclo para poder decir que la respuesta no queda desvirtuada.

Obsérvese que, análogamente al caso anterior, Td es el periodo de señal respuesta del sistema
continuo equivalente.
2.3.2 Correlación entre el plano S y el plano Z.
Recordando la definición de transformada Z: z = e Ts = e − σT ⋅ e ± jωT = e − σT ⋅ [cos(ωT) ± j ⋅ sin(ωT)]
Transformación de rectas de parte real constante:
jω
Im[z]
e σT 1
- 1
Re[z]
- σ1 σ2 σ
e σT
2
Fig. 2.13 Transformación de rectas de parte real constante del plano S al Z
Transformación de rectas de parte imaginaria constante:
jω
Im[z]
j ωs /2
j ω2 ω 2T
j ω1 ω 1T
1 Re[z]
σ 1
Fig. 2.14 Transformación de rectas de parte imaginaria constante del plano S al Z
Debe observarse que cuanto más pequeñas son la parte real (en valor absoluto) y la parte imaginaria
de las raíces en plano S, más cercanas se encuentran las raíces transformadas al punto z=1 en el plano
Z.
Como se ha comentado anteriormente, puede afirmarse que el sistema discreto y el sistema continuo
equivalente ofrecen dinámicas similares cuando el número de muestras por ciclo y el número de
muestras por constante de tiempo es suficientemente elevado (en este caso se dice que el sistema no
está desvirtuado). Estas indicaciones conllevan a que las raíces características en plano S sean tales
ωs T d τ 1
que mantengan elevadas las relaciones: = y o, equivalentemente, . En conclusión, un
ωd T T σT

sistema discreto no queda desvirtuado frente al sistema continuo equivalente cuando las raíces
transformadas en plano Z se encuentren en el interior del circulo de radio unidad y cercanas al punto
z=1. En estos casos, pueden asignarse al sistema discreto las mismas características dinámicas que
posee el sistema continuo equivalente. De hecho, no es posible conocer las características que presenta
la respuesta de un sistema discreto cuando las raíces en plano Z se encuentran alejadas del punto z=1,
si no es resolviendo la antitransformada Z de la señal de salida.
Mediante la observación de la expresión de la transformada Z, se puede afirmar que el aumento del

período de muestreo provocará un alejamiento de las raíces transformadas de la zona del plano Z en la
cual se garantiza una respuesta del sistema discreto no desvirtuada frente al sistema continuo
equivalente.
En la figura siguiente se muestran las respuestas obtenidas en función de la ubicación de las raíces en
plano S y en plano Z.
Eje imaginario
ωs
S Eje imaginario Z
j =20 1.5
2 8
15
1
10
0.5
5 1
2 1 65 2
7 6 4 5 3 8
0 0 7 4 3
6 5 2 1
-5 6 5 2 1
-0.5
-10
-1
-15
8
-20 -1.5
-2.5 -2 -1.5 -1 -0.5 0 0.5 1 -1.5 -1 -0.5 0 0.5 1 1.5
Eje real Eje real
1 2
500 90
80
0 70
60
-500 50
40
-1000 30
20
-1500 10
0
-2000 -10
0 5 10 15 20 25 30 0 200 400 600 800 1000
nº de muestras nº de muestras

3 4
10 7
9
6
8
7 5
6
4
5
3
4
3 2
2
1
1
0 0
0 2 4 6 8 10 0 5 10 15 20 25 30
5 6
25 12
10
20
8
15
10
4
5
2
0 0
0 5 10 15 20 25 30 0 5 10 15 20
7 8
4 2
3.5
1.5
3
1
2.5
0.5
2
0
1.5
-0.5
1
0.5 -1
0 -1.5
0 5 10 15 20 0 20 40 60 80 100 120 140 160
Fig. 2.15 Posiciones de polos correspondientes entre el plano S y el Plano Z y características de respuesta
transitoria.
2.4 Estabilidad absoluta de sis temas lineales.
2.4.1 Estabilidad en sistemas de tie mpo continuo.

Se ha comentado anteriormente que la estabilidad absoluta de un sistema lineal se logra cuando la
respuesta transitoria decae a cero al tender el tiempo a infinito. Se puede demostrar que para que un

sistema lineal sea estable es necesario que posea todos los polos de su función de transferencia en el
semiplano izquierdo del plano transformado S.
En sistemas de control, el problema fundamental es la determinación de las raíces del sistema en Lazo
Cerrado a partir del conocimiento de las raíces en Lazo Abierto. Recordando la expresión de la
función de transferencia en Lazo Cerrado.
CONTROL PLANTA
R(s) + E(s) M(s) C(s)
Gc(s) G(s)
-
B(s) H(s)
ELEMENTO DE MEDIDA
Fig. 2.16 Estructura de un sistema de control en lazo cerrado
C(s) Gc(s) ⋅ G(s)

GLC(s) = =
R(s) 1 + Gc(s) ⋅ G(s) ⋅ H(s)
La ecuación característica del sistema en Lazo Cerrado 1+Gc(s)G(s)H(s) = 0 no es fácilmente

resoluble en general; por ello aparecen métodos algorítmicos para poder determinar la estabilidad de
un sistema en lazo cerrado, el más importante de estos métodos es el criterio de estabilidad de Routh
(CER).
El criterio de estabilidad de Routh permite determinar el número de raíces de una ecuación de variable
compleja que se encuentran en el semiplano derecho, utilizándose, de este modo, para determinar si
existen polos de una función de transferencia en el semiplano derecho del plano transformado S.
Criterio de estabilidad de Routh (CER):
El criterio de estabilidad de Routh es un criterio de estabilidad absoluta. Se basa en la determinación

del número de raíces de un polinomio que se encuentran en el semiplano derecho del plano S. Para su
aplicación deben verificarse dos condiciones:
- Condición necesaria:
m m−1
C( s) b 0 ⋅ s + b1 ⋅ s +....+ bm
Dada la función de transferencia: = n n−1
, debe escribirse el denominador
R ( s) a 0 ⋅ s + a1 ⋅ s +....+ an
de la forma: a 0 ⋅ s n + a1 ⋅ s n −1 + .... + an = 0 , con an ≠ 0 (se eliminan las raíces en el eje imaginario).
Criterio: Si existe algún coeficiente negativo o cero en presencia de algún coeficiente positivo,
entonces existen una o más raíces imaginarias puras o con parte real positiva, lo cual implica que el
sistema es inestable. En otros términos, para garantizar estabilidad, a partir del primer ai ≠ 0 todos los
coeficientes deben estar presentes y ser positivos.

- Condición suficiente:
n n −1
Debe aplicarse el algoritmo de formación siguiente: a 0 ⋅ s + a1 ⋅ s +....+ an = 0 .
!
s n a0 a2 a4 a6
"
s n −1 a1 a 3 a 5 a 7
"
s n −2 b1 b2 b3 b4
s c c c c "
n−3
1 2 3 4
# # # # # "
s 2 e1 e2
s1 f 1
s 0 g1
Las primeras filas se obtienen directamente del polinomio característico, el resto de coeficientes se
obtienen según las expresiones:
a1 ⋅ a 2 − a 0 ⋅ a 3 a1 ⋅ a 4 − a 0 ⋅ a 5 a1 ⋅ a 6 − a 0 ⋅ a 7
b1 = ; b2 = ; b3 =
a1 a1 a1
b1 ⋅ a 3 − a1 ⋅ b2 b1 ⋅ a5 − a1 ⋅ b3 b1 ⋅ a 7 − a1 ⋅ b4
c1 = ; c2 = ; c3 =
b1 b1 b1
Como observación debe indicarse que puede multiplicarse o dividirse toda una fila por una constante
positiva.
El criterio de estabilidad de Routh determina que el número de raíces con parte real positiva del
polinomio estudiado es igual al número de cambios de signo de la primera columna del algoritmo de
formación.
De este modo la condición necesaria y suficiente para que un sistema sea estable es:
- Todos los coeficientes del polinomio característico deben existir y ser positivos.
- Todos los coeficientes de la primera columna del algoritmo de formación deben de ser positivos.
Ejemplo 2.3
a 0 ⋅ s 3 + a1 ⋅ s 2 + a 2 ⋅ s + a 3 = 0
Condición necesaria ⇒ a 0, a1, a 2, a 3 > 0

s3 a0 a2
Algoritmo: s2 a1 a3 , condición suficiente

a1 ⋅ a 2 − a 0 ⋅ a 3
s1
a1
s0 a3
a1 ⋅ a 2 − a 0 ⋅ a 3
⇒ > 0 ⇒ a1 ⋅ a 2 > a 0 ⋅ a 3
a1
Ejemplo 2.4
s 4 + 2 ⋅ s3 + 3 ⋅ s 2 + 4 ⋅ s + 5 = 0
Se verifica la condición necesaria
s4 1 3 5
s3 2 4
Algoritmo: , condición suficiente ⇒ existen dos cambios de signo ⇒
s2 1 5
s1 −6
s0 5
existen dos raíces con parte real positiva ⇒ el sistema es inestable.
* Casos especiales en el algoritmo de Routh.
1) Un término de la primera columna es cero en presencia de otros diferentes de cero.
Cuando esta situación ocurre, puede afirmarse que el sistema es inestable. Para determinar la
ubicación de las raíces que proporcionan la inestabilidad debe realizarse el siguiente procedimiento:
- Sustituir el cero por ε>0 con ε<<1.
- Aplicar el procedimiento habitual.
- Aplicar el criterio: Si los coeficientes superior e inferior a ε son de igual signo⇒existen raíces sobre
el eje imaginario; si los coeficientes superior e inferior tienen signo diferente⇒existen raíces en
semiplano derecho.

Ejemplo 2.5
s3 + 2 ⋅ s 2 + s + 2 = 0
Se verifica la condición necesaria.
s3 1 1
2
Algoritmo: s 2 2 , no hay cambio de signo ⇒ existen raíces en el eje imaginario.
s1 0≈ε
s0 2
* Ejemplo: s3 − 3 ⋅ s + 2 = 0
No se verifica la condición necesaria ⇒ el sistema es inestable.
s3 1 −3
Algoritmo: s2 0≈ε 2 , condición suficiente ⇒ existen dos cambios de signo ⇒

2
s1 −3−
ε
s0 2
existen dos raíces con parte real positiva. De hecho, s3 − 3 ⋅ s 2 + 2 = (s − 1) 2 ⋅ (s + 2) = 0
2) Todos los coeficientes de una fila son cero.
Cuando esta situación ocurre, puede afirmarse que el sistema es inestable. Existen raíces de igual
valor simétricas respecto a los ejes. Para determinar la ubicación de las raíces que proporcionan la
inestabilidad debe realizarse el siguiente procedimiento:
- Sustituir la fila de ceros por la derivada del “polinomio auxiliar” Pa(s).
Pa(s): polinomio formado por los coeficientes de la fila anterior a la de ceros. Debe indicarse que las
raíces de Pa(s) son raíces del polinomio característico P(s). P(s)=Pa(s)⋅Q(s)
- Aplicar el procedimiento habitual.
Ejemplo 2.6
s5 + 2 ⋅ s 4 + 24 ⋅ s 3 + 48 ⋅ s 2 − 25 ⋅ s − 50 = 0
No se verifica la condición necesaria.

s5 1 24 − 25
s5 1 24 − 25 s4
2 48 − 50
Algoritmo: s4
2 48 − 50 ⇒ s3 8 96 ,
s3 0 0 s2 24 − 50
s1 112.7
s0 − 50
dPa (s)
dado que Pa (s) = 2 ⋅ s 4 + 48 ⋅ s 2 − 50 ⇒ = 8 ⋅ s3 + 96 ⋅ s
ds
Existe un cambio de signo ⇒ existe una raíz con parte real positiva. Además, pueden obtenerse las
raíces simétricas calculando Pa (s) = 2 ⋅ s 4 + 48 ⋅ s 2 − 50 = 0 ⇒ s = ±1; s = ± j5 . La restante raíz se
determina dividiendo el polinomio inicial entre el polinomio auxiliar (obteniéndose s=-2).
Ejemplo 2.7
El criterio de estabilidad de Routh puede utilizarse para determinar el valor máximo de una ganancia
para el cual el sistema alcanza la estabilidad crítica o límite. Supóngase el sistema en lazo cerrado
siguiente:
R(s) + C(s)
G(s)
-
H(s)
Fig. 2.17 Diagrama de bloques de sistema en lazo cerrado
k ⋅ (s + 40) 1
G ( s) = ; H ( s) =
s ⋅ ( s + 10) s + 20
La función de transferencia en lazo cerrado es:

C( s)
=
(2
k ⋅ s + 60 ⋅ s + 800 )
R ( s) s3 + 30 ⋅ s2 + ( 200 + k ) ⋅ s + 40 ⋅ k
El rango de valores de ganancia k para el cual el sistema es estable quedará establecido al aplicar el
criterio de Estabilidad de Routh sobre el polinomio característico:
s3 + 30 ⋅ s 2 + (200 + k ) ⋅ s + 40 ⋅ k ⇒ Se verifica la condición necesaria para valores de k positiva.

s3 1 200 + k
Algoritmo: s2 30 40 ⋅ k , para que el sistema sea estable

30 ⋅ 200 − 10 ⋅ k
s1
30
s0 40 ⋅ k
⇒ k<600. En conclusión, el rango de valores de k para los que permanece el sistema estable es:
0<k<600.
Para k=600 el sistema se encuentra sobre la estabilidad límite siendo oscilatorio; puede determinarse
la frecuencia de oscilación calculando las raíces imaginarias que la caracterizan. Considerando el
polinomio auxiliar:
800
Pa (s) = 30 ⋅ s 2 + 40 ⋅ k = 0, cuando k = 600 ⇒ s = ± j 800 ⇒ fosc = (Hz)
2π
2.4.2 Estabilidad en sistemas de tie mpo discreto.

Se puede demostrar que un sistema discreto es estable cuando posee todos los polos de su función de
transferencia en el interior del circulo de radio unidad en el plano transformado Z.
La función de transferencia del sistema de control discreto de la figura se puede expresar:
CONTROL
DISCRETO ZOH PLANTA
R(s) + E(s) E*(s) C(s)
Gc(z) Goh(s) G(s)
- T T
H(s)
ELEMENTO DE MEDIDA
Fig. 2.18 Diagrama de bloques de sistema de control discreto en lazo cerrado
C(z) Gc(z) ⋅ GohG(z)

GLC(z) = =
R(z) 1 + Gc(z) ⋅ GohGH(z)
donde la ecuación característica del sistema en lazo cerrado 1+Gc(z)GohGH(z)=0 no es fácilmente

resoluble y deben buscarse métodos transformados para poder determinar la posición de sus raíces.
En este caso, la aplicación directa del criterio de estabilidad de Routh no es útil, porque determina el
número de raíces de la ecuación característica que se encuentran en semiplano derecho y no en el

exterior del circulo de radio unidad. Sin embargo, sí es posible aplicar el CER tras una transformación
que convierta el interior y el exterior del circulo de radio unidad en un semiplano izquierdo y un
semiplano derecho respectivamente; a esta transformación se le denomina transformación bilineal.
* Transformación bilineal.
La definición de la transformación bilineal no es única; sin embargo, la más conocida se realiza

T
1+ w
mediante el cambio de variable: z = 2 ; w = 2 ⋅ z −1 . (2.9)
T
1− w T z +1
2
La transformación del circulo de radio unidad ofrece el resultado:
2 z − 1 2 e jωT − 1 2 sen
ωT
2
( )
w= ⋅ = ⋅ jωT
T z +1 T e +1
=j ⋅
T cos ω T
2
( )
2  ωT 
w=j ⋅ tg   , que, efectivamente, corresponde con el eje imaginario del plano transformado W.
T  2 
jω S Z W
3
jω w
2
ω
j s
2
2
R→∞
ωs 1 2
j 3 2
4 j
1 3 1 T
σ 1 σw
Fig. 2.19 Relación entre los planos S, Z y W
ω w 100
90
80
70
60
50
40
30
20
10
0
0 0.5 1 1.5 ωs 2 2.5 3 ωs 3.5 ω
4 2
Fig. 2.20 Relación entre los ejes imaginarios del plano S y del plano W

En conclusión, el circulo de radio unidad del plano Z se transforma, mediante la transformada bilineal,
en el eje imaginario del plano W, de manera que puede aplicarse en dicho plano el criterio de
estabilidad de Routh para determinar la estabilidad de un sistema discreto.
Ejemplo 2.8
Puede analizarse la estabilidad del sistema discreto de la siguiente figura
CONTROL
DISCRETO ZOH PLANTA
R(s) + E(s) E*(s) C(s)
Gc(z) Goh(s) G(s)
- T T
H(s)
ELEMENTO DE MEDIDA
Fig. 2.21 Diagrama de bloques de sistema de control discreto en lazo cerrado
donde:
k
Gc( z ) = 1; H ( s) = 1; G (s) = ; T = 0.1 seg.
s ⋅ ( s + 1)
para ello, debe encontrarse la función de transferencia en lazo abierto como:
1 − e −Ts k  −1  k  0.0048z + 0.00468

GLA (z) = Z  ⋅  = (1 − z ) ⋅ Z 2  = k⋅
 s s ⋅ (s + 1)   s ⋅ (s + 1)  (z − 1) ⋅ (z − 0.905)
aplicando la transformada bilineal:
T
1+
w
2 = 1 + 0.05w ⇒ GLA ( w ) = k ⋅ − 0.00016w − 0.1872w + 3.81
2
z=
1 − w 1 − 0.05w 3.81w 2 + 3.8w
T
2
obteniendo la ecuación característica: 1 + GLA ( w ) = 0
(3.81 − 0.00016 ⋅ k ) ⋅ w 2 + (3.8 − 0.1872 ⋅ k ) ⋅ w + 3.81 ⋅ k = 0

El resultado, aplicando el criterio de estabilidad de Routh, indica que el sistema es estable para los
valores de k: 0<k<20.3. En concreto, para k=20.3 se pueden obtener los puntos de corte con el eje
imaginario en plano W para, posteriormente, determinar la frecuencia de oscilación del sistema
discreto. Así, considerando del polinomio auxiliar:
2
( 3.81 − 0.00016 ⋅ 20.3) ⋅ w + 3.81⋅ 20.3 = 0 ⇒ w = ± j4.5075
2 − 1 ωwT 
2 − 1 4.5075 ⋅ 0.1 
La frecuencia de oscilación es: ω=  ⋅ tg = ⋅ tg   = 4.433 rad / seg ;
T  2  0.1  2 
resultado muy parecido a la frecuencia que se obtenía en plano W debido a que la frecuencia de
muestreo es suficientemente elevada (ωs = 2 π T = 62.83 rad / seg.).
Debe observarse que, así como un sistema continuo de segundo orden con realimentación unitaria
siempre es estable (partiendo de un sistema en lazo abierto estable), un sistema discreto no tiene por
qué serlo debido al efecto desestabilizador del periodo de muestreo. En concreto, realizando el mismo
ejemplo con un periodo de muestreo mayor (por ejemplo T=1 seg), el margen de valores de ganancia
para garantizar estabilidad disminuye (0<k<2.39).
2.5 Análisis en régimen estacio nario.
Un sistema lineal estable alcanza el régimen o estado estacionario cuando, al ser excitado por una
señal de entrada, la respuesta transitoria decae a cero.
En sistemas de control, la precisión o exactitud del sistema se convierte en una de las especificaciones
más importantes que verificar; el sistema de control debe 'seguir' la señal de referencia en estado
estacionario del modo más preciso posible. Por esta razón, en sistemas de control en lazo cerrado se
obtienen las expresiones de los errores estacionarios del sistema en función del tipo de señal de
referencia introducida y de las funciones de transferencia que contiene.
Analizando el diagrama de bloques de un sistema en lazo cerrado se obtiene la señal de error como:
R ( s)
E ( s) =
1 + GLA ( s)
s ⋅ R (s )
Se define el error en estado estacionario: ess ≡ lim e( t ) = lim s ⋅ E (s) = lim
t →∞ s →0 s →0 1 + GLA(s)
Particularizando para los diversos tipos de entradas:
* Error estacionario de posición:
1
r ( t ) = u( t ) → essp = (2.10)
1 + kp

kp = lim GLA (s) → coeficiente de posición.

s →0
* Error estacionario de velocidad:
1
r ( t ) = t → essv = (2.11)
kv
kv = lim s ⋅ GLA(s) → coeficiente de velocidad.

s →0
* Error estacionario de aceleración:
t2 1
r(t) = → essa = (2.12)
2 ka
ka = lim s 2 ⋅ GLA (s) → coeficiente de aceleración.

s→0
Obsérvese que, cuanto mayores son los coeficientes de error, menor es el error estacionario
correspondiente cometido. El error en estado estacionario en lazo cerrado puede obtenerse a partir de
la función de transferencia en lazo abierto, de modo que, cuanto mayor es su valor en continua (para
frecuencia cero) menor es el error cometido. Para poder eliminar el error de posición es necesario que
la función de transferencia en lazo abierto contenga un polo en origen; en este caso se dice que existe
un elemento integrador en dicha función.
Por último, señalar que este error estacionario no es más que el valor que adquiere la señal de salida
del detector de error en estado estacionario. En el caso en el cual la realimentación fuese unitaria,
coincidiría con la diferencia que se tendría en estado estacionario entre la señal de referencia y la señal
de salida del sistema. Si la realimentación no es unitaria, existirá un factor de proporcionalidad entre
ambas señales que vendrá determinado por el elemento de medida.
Las expresiones de los errores estacionarios de un sistema de control discreto en lazo cerrado son:
R ( z)
Análogamente al caso anterior: Señal de error: E ( z) =
1 + GLA ( z )
(1 − z −1 ) ⋅ R (z)
Error en estado estacionario: ess ≡ lim e(kT) = lim(1 − z −1 ) ⋅ E(z) = lim
k →∞ z →1 z →1 1 + GLA ( z )
* Error estacionario de posición:
1
r ( t ) = u ( t ) → essp = (2.13)
1 + kp

kp = lim GLA (z) → coeficiente de posición.

z →1
* Error estacionario de velocidad:
1
r ( t ) = t → essv = (2.14)
kv
z −1
kv = lim ⋅ GLA (z) → coeficiente de velocidad.
z →1 T
* Error estacionario de aceleración:
t2 1
r(t) = → essa = (2.15)
2 ka
2
(z − 1)
ka = lim ⋅ GLA (z) → coeficiente de aceleración.
z →1 T2
Ejemplo 2.9
La figura 2.22 muestra los elementos de un control automático de ganancia (CAG) que se utiliza para
mantener constante el nivel de salida de un amplificador.
Ei Amplificador
Eo
de Ganancia
Variable (K)
Egc R1
ka
R
τ=RC R2
C
Filtro Pasa-Bajos - Vcc
Fig. 2.22 Diagrama de un control automático de ganancia (CAG)
donde:
Ei ≡ Amplitud de la señal de entrada (Volts).

Eo ≡ Amplitud de la señal de salida (Volts).
Egc ≡ Tensión de control de la ganancia variable (Volts).
La amplitud de la señal de entrada Ei se multiplica por la ganancia variable K (que es función de la

tensión de control de ganancia Egc) para obtener la amplitud de la señal de salida Eo (Eo=Ei⋅K).

La figura 2.23 muestra el diagrama de bloques resultante del CAG:
Ei Eo
K[Egc]
Egc
+ Er
1  R2  R1 Vcc
  ⋅ ka
1+ sτ  R1 + R 2  - R2
Fig. 2.23 Diagrama de bloques de un CAG.
donde: K[Egc] = d ⋅ e − b⋅Egc
Este sistema es altamente no lineal (debido a la ganancia variable), pero puede linealizarse utilizando
un modelo incremental para obtener el diagrama de bloques lineal de la figura 2.24.
∆(20 log Ei) (dB) + ∆Eo

0.1151 Eo (volts/dB)
+
∆(20 log Egc) (dB)
−
(
d 20 log K[Egc] ) (dB / volts)
dEgc
1  R2  - ∆Er (volts)
  ⋅ ka
1+ sτ  R1 + R 2  +
Fig. 2.24 Modelo incremental de un CAG.
Donde, en este caso, Eo es el valor deseado de la amplitud de la señal de salida ( Eo = Eo ± ∆Eo) .
Obtener:
1.- La función de transferencia en lazo abierto.

2.- La frecuencia de transición del sistema. (Corte por 0 dB del módulo de la función en lazo abierto).
3.- Las funciones de transferencia:
∆Eo
• cuando ∆Er=0
∆ (20 log Ei)

∆Eo
• cuando ∆(20logEi)=0
∆Er
Determinar el valor de ∆Eo en estado estacionario, en ambos casos, cuando las entradas son tipo
escalón.
4.- Considerando los valores de los parámetros siguientes: b=1, R2=R1, C=100nF, y que la fuente de
alimentación está perfectamente regulada, determinar los valores de la ganancia del amplificador
ka y la resistencia del filtro pasa-bajos R si se desea una amplitud de salida Eo=5 V ± 5 mV
cuando la entrada Ei varía de forma tipo escalón entre 1 mV y 1 Volt y la constante de tiempo del
sistema en lazo cerrado debe ser de 40 ns.
5.- Describir brevemente el funcionamiento de este sistema. ¿En qué aplicaciones puede utilizarse?.
Razonar la respuesta.
Solución:
1.- Recordando que no es posible obtener una caracterización dinámica de un sistema no lineal a
través de una función de transferencia, se ha obtenido una descripción lineal del sistema mediante el
modelo incremental del mismo. La función de transferencia en lazo abierto es única y es el resultado
del producto de todas las funciones de transferencia del lazo del sistema.
0.1151 ⋅ ka  d (20 log K ) 

Definiendo: ko = ⋅−  ⋅ Eo
R1  dEgc 
1+
R2
dado que:
d ( 20 log(d ⋅ e − b⋅ Egc ))
−
d ( 20 log K )
=− =−
d
[20 log d − b ⋅ Egc ⋅ 20 log e] = b ⋅ 20 log e
dEgc dEgc dEgc
0.1151 ⋅ ka
En conclusión: ko = ⋅ Eo ⋅ b ⋅ 20 log e
R1
1+
R2
ko
Función de transferencia en lazo abierto: GLA (s) =
1+ τ ⋅s
ko
2.- Dado que el sistema es de primer orden, la frecuencia de transición es: ωt =
τ
3.- El modelo incremental es lineal, como hemos mencionado anteriormente, y por ello permite
aplicar superposición para obtener la expresión de la señal de salida en función de las señales de
entrada.

Funciones de transferencia:
∆Eo
• cuando ∆Er=0
∆ (20 log Ei)
∆Eo 0.1151 ⋅ Eo 0.1151 ⋅ Eo ⋅ (1 + τ ⋅ s )

= =
∆(20 log Ei) 1 + ko ko + 1 + τ ⋅ s
1+ τ ⋅s
Valor estacionario ante entrada tipo escalón:
∆Eo 0.1151 ⋅ Eo 0.1151 ⋅ Eo

lim = ⇒ ∆Eoss = ⋅ ∆ (20 log Ei)
s →0 ∆(20 log Ei) ko + 1 ko + 1
Debe observarse que: ko ↑⇒ ∆Eoss → 0 , es decir, que para una ganancia elevada, si la fuente de
referencia no fluctúa (∆Er=0), no existirán variaciones de amplitud en la señal de salida.
∆Eo
• cuando ∆(20logEi)=0
∆Er
ko
∆Eo
= 1+ τ ⋅s =
ko
∆Er 1 + ko ko + 1 + τ ⋅ s
1+ τ⋅s
Valor estacionario ante entrada tipo escalón:
∆Eo ko ko
lim = ⇒ ∆Eoss = ⋅ ∆Erss
s → 0 ∆Er ko + 1 ko + 1
Debe observarse que: ko↑ ⇒ ∆Eoss → ∆Erss , es decir, que para una ganancia elevada, si no existen
variaciones en la amplitud de la señal de entrada (∆(20logEi)=0), podemos reducir las variaciones de
amplitud de salida únicamente a la variación debida a la fuente de referencia (∆Erss).
4.- Diseño:
- Valores de los parámetros: b=1, R2=R1 y C=100nF.

- La fuente de alimentación está perfectamente regulada, lo cual implica que ∆Er=0
- Amplitud de salida Eo=5 V ± 5 mV, lo cual implica Eo = 5V, ∆Eo = 5mV .
- Ei varía de forma tipo escalón entre 1 mV y 1 Volt, lo cual implica ∆Ei=1000.
- La constante de tiempo del sistema en Lazo Cerrado debe ser de 40 ns.

0.1151 ⋅ Eo
∆Ei = 1000 ⇒ 20 log ∆Ei = 60 dB ∆Eoss = ⋅ ∆ (20 log Ei)
ko + 1
0.1151 ⋅ 5
5 ⋅ 10−3 = ⋅ 60 ⇒ ko = 6.905 ⋅ 103
ko + 1
0.1151 ⋅ ka
ko = ⋅ Eo ⋅ b ⋅ 20 log e ⇒ ka = 2762.7
R1
1+
R2
Constante de tiempo del sistema en Lazo Cerrado:
τ
τ LC = = 40nseg = 4 ⋅ 10 −8 ⇒ τ = 2.7624 ⋅ 10 − 4
ko + 1
τ = R ⋅C ⇒ R = 2762.4Ω
5.- En conclusión, se puede deducir que este sistema es válido para obtener señales de salida dentro de
un determinado margen de valores dependientes de la buena regulación de la fuente de alimentación;
evitando el riesgo de variaciones provenientes de la señal de entrada. La ganancia del sistema se
autoregula para lograr mantener el valor final de salida aproximadamente constante, de este modo este
sistema es útil para aquellas aplicaciones en las cuales la señal de entrada sea muy variable o de
amplitud poco predecible.
2.6 El lugar geométrico de las raíces (L.G.R.).
Debido a la necesidad de conocer los polos de un sistema en lazo cerrado, ya que determinarán las
características básicas de la respuesta transitoria, se desarrolla el método del lugar geométrico de las
Raíces (también denominado Lugar de Evans). Este método permite ubicar en un gráfico los polos de
un sistema en lazo cerrado a partir del conocimiento de los polos en lazo abierto, en función de un
parámetro variable. Para ello considérese la ecuación característica de un sistema en lazo cerrado:
1+ Gc(s)G(s)H(s) = 0 ⇒ 1+ GLA(s) = 0 ⇒ GLA(s) = -1 .
La resolución de esta ecuación implica la verificación de dos condiciones:
- Condición de ángulo: fase{GLA ( s)} = ± 180 ( 2λ + 1); λ ∈ N

o
(2.16)
- Condición de módulo: GLA (s) = 1 (2.17)
Debe observarse que, de este modo, se pasa del estudio del sistema en lazo cerrado al estudio de
características del sistema en lazo abierto, lo cual debe permitir mayor facilidad en el cálculo.

Se define el lugar geométrico de las raíces como el conjunto de puntos del plano S en los que se
verifica la condición de ángulo.
En conclusión, un punto que pertenece al lugar geométrico de las raíces es un posible polo del sistema
en lazo cerrado; para ello únicamente es necesario validar la condición de módulo, y ésta se cumplirá
para un valor determinado de la ganancia del sistema en lazo abierto. Sin embargo, un punto del plano
S que no pertenezca al LGR no puede ser polo en lazo cerrado porque no verifica la condición de
ángulo, aunque varíe la ganancia del sistema en lazo abierto. De este modo, variando el parámetro k,
(s + z1).....(s + zm )
donde: GLA (s) = k ⋅ ; 0 ≤ k < ∞ , se logra trazar el lugar geométrico de las raíces que
(s + p1).....(s + pn )
proporciona los valores de los polos en lazo cerrado en función de k.
2.6.1 Reglas de construcción del LG R.

Dado que los polos y ceros complejos de la función de transferencia en lazo abierto tienen asociados
sus complejo-conjugados, el LGR será simétrico respecto al eje real.
1) Trazar el diagrama polos-ceros en lazo abierto.
Ecuación característica:
b 0 ⋅ s m + b1 ⋅ s m −1 + .... + bm (s + z1).....(s + zm )
1 + GLA (s) = 0 ⇒ 1 + n −1
= 0 ⇒ 1+ k ⋅ =0
n
a 0 ⋅ s + a1 ⋅ s + .... + an (s + p1).....(s + pn )
donde: s=-zi son ceros y s=-pi son polos en lazo abierto
2) Puntos de inicio y final del LGR.
El trazado del lugar geométrico de las raíces se inicia en los polos en lazo abierto y finaliza en los
ceros en lazo abierto (en este caso, deben considerarse los ceros infinitos). Puede demostrarse esta
sentencia resolviendo:
- Inicio en polos en lazo abierto:
(s + z1).....(s + zm ) 1  (s + z1).....(s + zm )   1
= − ⇒ lim   = klim − =∞
(s + p1).....(s + pn ) k k → 0  (s + p1).....(s + pn )  →0 
 k 
para que la expresión anterior sea cierta es necesario que s → −pi

En conclusión, los polos en lazo cerrado coinciden con los polos en lazo abierto cuando k=0. Debe
indicarse que, lógicamente, este efecto es relevante únicamente a nivel analítico, dado que no es
posible tener k=0 a nivel real.
- Final en ceros en lazo abierto:
(s + z1).....(s + zm ) 1  (s + z1).....(s + zm )   1
= − ⇒ lim   = klim − =0
(s + p1).....(s + pn ) k k → ∞  (s + p1).....(s + pn )  →∞ 
 k 
Para que la expresión anterior sea cierta es necesario que s → −zi ó s → ∞ (en el caso para el cual el
grado del denominador sea mayor que el grado del numerador).
En conclusión, los polos en lazo cerrado coinciden con los ceros en lazo abierto cuando k→∞. Así el
lugar geométrico puede tener ramas que finalicen en infinito; ahora bien, dado que el sistema es causal
nunca puede iniciarse el LGR en infinito.
El LGR se origina en los polos en lazo abierto y finaliza en los ceros en lazo abierto (finitos e
infinitos). El número de ramas del lugar geométrico indica el número de polos en lazo cerrado y
coincide con el número de polos en lazo abierto y el número de ceros en lazo abierto (finitos e
infinitos).
3) Lugar geométrico de las raíces sobre el eje real.
Los polos y ceros complejo-conjugados no afectan en la evaluación del LGR sobre el eje real, dado
que en su contribución suman múltiplos de 360o. Observando, únicamente, los polos y ceros en lazo
abierto sobre el eje real, puede aplicarse la siguiente consideración: un punto del eje real pertenece al
LGR cuando el número de total de polos y ceros a su derecha es impar (la suma angular total será un
múltiplo de 180o).
jω
j ω1
- σ1 σ
Fig. 2.25 Contribución angular de polos complejos-conjugados sobre un punto del eje real

4) Asíntotas del LGR.
El estudio asintótico se realiza para s → ∞ . En ese caso, las contribuciones angulares por parte de
todas las raíces son prácticamente iguales, y existe un efecto de cancelación de contribución angular
entre polos y ceros. De este modo, la expresión de los ángulos de las asíntotas vendrá dada por:
±180 ⋅ ( 2 λ + 1)
o
∠Asíntota = (2.18)
n−m
donde: n y m son los grados de denominador y numerador de la función de transferencia en lazo

abierto, respectivamente, y λ es un número natural.
jω
Fig. 2.26 Contribución angular sobre una asíntota
Para demostrar esta expresión debe considerarse la ecuación característica:
Gc(s)G(s)H(s) = -1 ⇒ GLA(s) = -1
(s + z1).....( s + zm)  (s + z1).....(s + zm )  k

k⋅ = −1 ; k ⋅ lim   → n −m
( s + p1).....( s + pn ) s → ∞  (s + p1).....(s + pn )  s
 k 
Recordando la condición de ángulo: fase n − m  = ±180o (2λ + 1); λ ∈ N
s 
±180 ( 2 λ + 1)
o
fase{s} = ; λ∈ N
n−m
 ±60
o
λ=0

±180 ( 2 λ + 1) 
o
Así, si por ejemplo: n-m=3; fase{s} =
o
= ±180 λ=1
3  o o
±300 = ±60 λ = 2

El punto de intersección de las asíntotas con el eje real, necesario para poder realizar el trazado de las
asíntotas, viene dado por la expresión :
n m
∑ p i − ∑ zj
i =1 j =1
σa = (2.19)
n−m
jω
σa σ
Fig. 2.27 Ubicación de las asíntotas
Debe observarse que en el caso de n-m=1, esto es, poseer únicamente una asíntota, no debe calcularse
el punto de intersección, dado que todo el eje real constituye la propia asíntota.
5) Puntos de ruptura.
Por definición, un punto de ruptura en el LGR corresponde a una raíz múltiple de la ecuación
característica, esto es, un punto de ruptura implica un polo en lazo cerrado múltiple. Debe resaltarse
que los puntos de ruptura pueden ser reales o complejo conjugados. Los puntos de ruptura pueden
dividirse en puntos de ruptura de dispersión (en los cuales el valor de k alcanza un máximo relativo) y
puntos de ruptura de confluencia (para los cuales k alcanza un mínimo relativo).
Para determinar el procedimiento de cálculo de puntos de ruptura, se realiza la evaluación de k cuando

aparece un punto de ruptura sobre el eje real, diferenciándose los casos:
1. LGR sobre eje real entre dos polos. En este caso, el punto de ruptura aparece cuando k alcanza un
máximo relativo, determinándose según la expresión:
(s + z1).....(s + zm ) A (s) B(s)

1 + GLA (s) = 0 ⇒ 1 + k ⋅ = 0 ⇒ 1+ k ⋅ =0⇒k =−
(s + p1).....(s + pn ) B(s) A (s )
dk
= 0 ⇒ B' (s) ⋅ A(s) − A' (s) ⋅ B(s) = 0
ds

Las soluciones de la ecuación anterior son puntos de ruptura si pertenecen al LGR y la k asociada es
real y positiva. Obviamente, las soluciones de la ecuación anterior pueden proporcionar puntos de
ruptura complejo conjugados.
2. LGR sobre eje real entre dos ceros. En este caso, el punto de ruptura aparece cuando k alcanza un
mínimo relativo, determinándose análogamente al caso anterior.
3. LGR sobre eje real entre cero y polo. En este caso, existe la posibilidad de que no aparezcan puntos
de ruptura, o bien, que existan en pares de dispersión y confluencia.
k evaluada sobre eje real k evaluada sobre eje real

k→∞ k→∞
kmáx
kmin
σ σ
k=0 k=0
k evaluada sobre eje real

kmáx
k →∞
kmin
k=0 σ
Fig. 2.28 Análisis de puntos de ruptura sobre el eje real
6) Puntos de cruce del LGR con el eje imaginario.
Métodos:
1- Sustituir s=jω en la ecuación característica, igualando parte real e imaginaria a cero.
Re(ω) = 0
1 + G (s)H (s) = 0 ⇒ 1 + G ( jω)H ( jω) = 0 ⇒ Re(ω) + j Im(ω) = 0 ⇒ 
Im(ω) = 0

2- Aplicando el algoritmo de Routh, anulando una fila de coeficientes.
Por ejemplo: s3 + b ⋅ s 2 + c ⋅ s + k ⋅ d = 0
s3 1 c
b ⋅ c − k ⋅ d b. c
Algoritmo: s2 b k ⋅d , anular filas ⇒ 
 b
= 0⇒ k =
d
b⋅c − k ⋅d  k=0
s1
b
s0 k ⋅d
Polinomio auxiliar: Pa (s) = b ⋅ s 2 + k ⋅ d = 0 ⇒ s = ± j c
7) Angulos de arranque y llegada.
Los ángulos de arranque del LGR de los polos en lazo abierto y los ángulos de llegada del LGR a los
ceros en lazo abierto se determinan a partir de la distribución del diagrama polos-ceros en lazo abierto.
Para ello, se presupone un punto perteneciente al LGR suficientemente cercano a la singularidad sobre
la que se quiere determinar el ángulo de partida o llegada como para poder considerarlo en la misma
posición que la propia singularidad; de este modo, al aplicar la condición de ángulo todas las
contribuciones angulares serán conocidas exceptuando el ángulo de arranque o llegada incógnita.
θ2? jω
P2
φ1 θ1 θ0
Z1 P1 P0 σ
θ3
P3
Fig. 2.29 Análisis de los ángulos de arranque

Ejemplo 2.10
k ⋅ (s − 2 )
Calcular el LGR de un sistema con GLA (s) = 3
s 2 ⋅ (s 2 + 10 ⋅ s + 7)
3
1) Polos y ceros en lazo abierto.
 z = 2  p1 = 0 (doble)
Ceros:  1 3 Polos:  5
z 2 = ∞ (triple) p 2,3 = − 3 ± j2.06
2) LGR en el eje real
jω
Fig. 2.30 LGR sobre el eje real
3) Asíntotas
n= 4 (nº de polos finitos) y m=1(nº de ceros finitos)
 ±60
o
λ=0

±180 ( 2 λ + 1) 
o
fase{s} =
o
= ±180 λ=1
3  o o
±300 = ±60 λ = 2
4) Intersección de las asíntotas con el eje real.
(0 + 0 − 5 3 + j2.06 − 5 3 − j2.06) − 2 3
n m
∑ pi − ∑ zj
i =1 j =1
σ = = = −1.33
a n−m 4−1
Al no existir ninguna rama entre polo y cero, no se cumple la condición necesaria para que existan
puntos de ruptura. Por tanto sólo se calcularán al final si las condiciones geométricas lo requieren.
6) Angulos de arranque y llegada

A continuación se calcula el ángulo de arranque de las ramas de los polos complejo conjugados. Para
ello, en el polo − 5 3 + j2.06 se toma un punto P de su entorno tal que P∈LGR.
[ ]
fase{P} = θ 1 − 2 ⋅ φ 1 + φ 2 + φ 3 = ±180(2 ⋅ λ + 1)
Calculo de φ1
2.06
180º − φ1 = arctg ⇒ φ1 = 129º
167
.
Calculo de θ1
2.06
180º −θ1 = arctg ⇒ θ1 = 136.8º
2.2
φ3=90º
Calculo de φ2
[ ]
± 180º ( 2 ⋅ λ + 1) = 136.8º − 2 ⋅ 129º + φ 2 + 90º ⇒ φ 2 = −312
. º
7) Puntos de cruce con el eje jω.
Para encontrar el punto de cruce del LGR con el eje jω se utiliza el criterio de estabilidad de Routh.
Se busca la ecuación característica:
10 3 2
s4 + s + 7s2 + ks − k = 0
3 3
Se aplica el algoritmo de Routh.
2
s4 1 7 − k
3
10
s3 k
3
3k 2
s2 7− − k
10 3
3k 2 83
− + k
10 9
s1
3k
7−
10
2
s0 − k
3

A continuación se busca una k que anule alguna fila.
• k=0 anula la fila s0. Pero esto indica que hay un polo en el origen que ya era conocido.
• Se busca el valor de k que anula la fila s1
3k 2 83
− + k = 0 ⇒ k = 30.7
10 9
A continuación se busca las raíces del polinomio auxiliar a esta fila.
 3k  2
 7 −  s2 − k =0
 10  3 k = 30.7
− 2.21s 2 − 20.46 = 0 ⇒ s1, 2 = ± j3.04
El LGR cortará al eje jω por ±j3.04.
• Para la fila s2 no hay un único valor de k que anule a la fila, por tanto no habrá otros posibles
puntos de corte.
Eje imaginario
4
3
φ2?
P
2 -31.2
1
θ1
60º φ1
0
-60º
-1
φ3
31.2
-2
-3
-4
-4 -3 -2 -1 0 1 2 3 4
Eje real
Fig. 2.31 Lugar geométrico de las raíces.
Ejemplo 2.11
k ⋅ (s + 2)
Calcular el LGR de un sistema con GLA (s) =
(s + 1) 2
1) Polos y ceros en lazo abierto.

 z = −2
Ceros:  1 polos: {− 1 doble
z 2 = ∞
2) LGR en el eje real
jω
Fig. 2.32 LGR sobre el eje real
3) Asíntotas
n= 2 (nº de polos finitos) y m=1(nº de ceros finitos)
± 180 o (2λ + 1)  ± 180 o λ = 0

fase{s} = =
± 540 = ±180 λ =1
o
1
La única asíntota es el eje real negativo.
4) Intersección de las asíntotas con el eje real.
Serán todo los puntos del semieje real negativo.
Existe una rama entre dos ceros, por tanto existirá como mínimo un punto de ruptura.
De la ecuación característica:
− ( s + 1)
2
k=
( s + 2)
2( s + 1)( s + 2) − ( s + 1)
2
dk
=− =0
ds ( s + 2) 2
 s r 1 = −1
s2 + 4s + 3 = 0 ⇒ 
 s r 2 = −3
( ) ( )
Los dos puntos pertenecen al LGR y poseen k s r1 = 0 y k sr 2 = 4 mayores o iguales a cero y por
tanto son puntos de ruptura.

6) Ángulos de arranque y llegada.
Ángulo de llegada a los ceros: ±180º
Ángulo de arranque de los polos:

θ 1 − 2φ 1 = ±180º ⇒ φ 1 = ±90º
θ1=0º
7) Puntos de cruce con el eje jω.
Para encontrar el punto de cruce del LGR con el eje jω se utiliza el criterio de estabilidad de Routh.
Se busca la ecuación característica:

s2 + ( k + 2)s + 2 k + 1 = 0
Se aplica el algoritmo de Routh:
s2 1 2k + 1
s1
k+2
s0 2k + 1
A continuación se busca una k que anule alguna fila.
1
• k=− anula la fila s0. Pero esto no es posible.
2
• k = −2 anula la fila s1. Pero esto no es posible.
Se comprueba que no corta con el eje jω.
Eje imaginario
3
θ1 φ1?
0
-1
-2
-3
-4 -3 -2 -1 0 1
Eje real
Fig. 2.33 Lugar geométrico de las raíces

Ejemplo 2.12
Considérese el siguiente sistema correspondiente a un intercambiador de calor descrito por la función

de transferencia G1(s)
+
e-Ts K·G1(s)
-
1
G1 ( s) =
s
Fig. 2.34 Diagrama de bloques de un intercambiador de calor.
− Ts
El término e es el retardo provocado por el sensor de temperatura que está físicamente un poco
más abajo del intercambiador, de modo que su lectura se retarda T seg.
Se desea construir el lugar geométrico de las raíces para un valor fijo de T. Para ello:
1.- Escribir la ecuación característica del sistema.

2.- Encontrar las condiciones de ángulo y módulo.
3.- La ecuación característica del sistema tiene un número infinito de raíces por ser una ecuación
transcendental. Determinar el número de ramas del lugar geométrico de raíces.
4.- Demostrar que la ecuación que determina los puntos donde el LGR cruza al eje imaginario en
frecuencias positivas es:
1 π 
ω=  + λ ⋅ 2 π  , y calcular dichos puntos para λ=0,1,2.
T 2 
Recordar que el LGR es simétrico respecto al eje real del plano s.

5.- Demostrar que los puntos del LGR para K=0 están en los polos de G1(s) y/o σ = - ∞. (s=σ+jω).
Dibujar el LGR sobre el eje real, indicando las direcciones de las ramas y la posición de los polos
en lazo cerrado en K=0.
6.- Demostrar que los puntos del LGR para K=∞ están en los ceros de G1(s) y/o σ=∞ (s = σ+jω).
Las asíntotas del LGR son infinitas en número y todas ellas paralelas al eje real en el plano S.
7.- Demostrar que los puntos de cruce de las asíntotas con el eje imaginario positivo vienen dados por
las expresiones:
2 πλ
• Asíntotas para σ=-∞: ω=
T

1
• Asíntotas para σ=∞: ω=
T
[π + 2πλ]
Calcular las asíntotas para λ=0,1,2.
Considerar: Fase[G1(s)]=-180º|s=-∞+jω
Fase[G1(s)]=0º|s=∞+jω
8.- Los puntos de ruptura del LGR deben satisfacer

[
d G1(s)K e -Ts ] = 0 . Calcular los puntos de
ds
ruptura.
9.- Dibujar el Lugar Geométrico de Raíces completo.
Solución:
K −Ts
1.- 1 + KG1(s)e −Ts = 0 ; 1+ e =0
s
2.- s=σ +jω
1 + Ke− Tσe − jωT G1(s) = 0

Ke − Tσe − jωT G1(s) = −1
Condición de módulo: Ke− Tσ G1(s) = 1
Condición de ángulo: −ωT + fase{G1(s)} = ±180( 2λ + 1)
fase{G1(s)} = ±180( 2λ + 1) + ωT
3.- Si la ecuación característica tiene infinitas raíces, habrá infinitos polos en lazo cerrado, y para cada
polo en lazo cerrado habrá una rama. Por tanto habrá infinitas ramas.
4.- Los puntos donde el LGR cruza al eje imaginario tiene que cumplir la condición de ángulo.
Por tanto aplicando la condición de ángulo a los puntos s=jω podremos encontrar ω.
fase{G1(s)} s= jω = −90 −90 = ±180( 2λ + 1) + ωT ω=

1
T
[−90 B 180( 2 λ + 1) ]
1 π 
ω=
1
T
[ T 2
]
−90 + 180( 2 λ + 1) =  + λ 2 π 


1  3π 
Descartamos la solución ω =
T
1
[ T 2
]
−90 − 180( 2λ + 1) = −  + λ 2 π  , porque hemos aplicado la

condición de ángulo a frecuencias positivas.
1 π 
Por tanto con ω = + λ2π
T  2 
π 5π 9π
λ=0 ⇒ ω = λ=1 ⇒ ω = λ=2 ⇒ ω =
2T 2T 2T
5.- Aplicamos la condición de módulo.
1
e − Tσ G1(s) =
K
Si K→ 0 ⇒ e − Tσ G1(s) = ∞
Por tanto s tiende a los polos de G1(s) ó σ tiende a -∞.
Para el eje real el LGR es:

jω
K=0
K=0
6.- Aplicamos condición de módulo:
1
e − Tσ G1(s) =
K
Si K→ ∞ ⇒ e − Tσ G1(s) = 0
Por tanto, s tiende a los ceros de G1 (s) ó σ tiende a ∞.
7.- Los puntos de s=-∞+jω deben pertenecer al LGR y por tanto cumplir la condición de ángulo.
• Puntos σ→ -∞
fase{G1(s)} = ±180( 2λ + 1) + ωT ; −180 = ±180( 2λ + 1) + ωT
Para: −180 = −180( 2λ + 1) + ωT

1 2 πλ
ω=
T
[360λ ] =
T

Para: −180 = +180( 2λ + 1) + ωT
1 2 π(1 + λ )
ω=−
T
[360 + 360λ ] = − T
Como calculamos sobre el eje positivo, descartamos frecuencias negativas.
2 πλ
Tendremos asíntotas para σ→ -∞ en ω =
T
2π 4π
λ=0 ⇒ ω = 0 λ=1 ⇒ ω = λ=2 ⇒ ω =
T T
• Puntos σ→ ∞
fase{G1(s)} = ±180( 2λ + 1) + ωT
0 = ±180( 2λ + 1) + ωT
Para: 0 = −180( 2λ + 1) + ωT
ω=
1
T
[ 1
]
180( 2 λ + 1) = [π + 2 πλ ]
T
Para: 0 = +180( 2λ + 1) + ωT
ω=−
1
T
[ 1
]
180( 2 λ + 1) = − [π + 2 πλ ]
T
Como calculamos sobre el eje positivo, descartamos frecuencias negativas.
1
Tendremos asíntotas para σ→ -∞ en ω = −
T
[π + 2πλ ]
π 3π 5π
λ=0 ⇒ ω = λ=1 ⇒ ω = λ=2 ⇒ ω =
T T T
d Ke − Ts 1
8.- =0 e − Ts ( Ts + 1) = 0 s=−
ds s T

9.-
jω
5π/Τ
4π/Τ
3π/Τ
2π/Τ
π/Τ
σ
1/Τ
Fig. 2.35 Diagrama del lugar geométrico de raíces.
2.6.2 Evaluación de Ceros en lazo c errado.
La determinación de los ceros en lazo cerrado debe realizarse cuidadosamente, dado que no aparecen
de un modo directo en el LGR y su influencia es importante en la respuesta del sistema. Observando
las topologías que se muestran en la figura 2.36, puede decirse que los ceros en lazo cerrado son los
ceros en lazo directo conjuntamente con los polos del elemento de medida.
R(s) + k(s+2) C(s) R(s) + k(s+2) 1 C(s)
- (s+1) - (s+1) (s+1)
1
(s+1)
Fig. 2.36 Evaluación de los ceros en lazo cerrado
Las funciones de transferencia en lazo cerrado de los sistemas descritos son:
C( s) k ⋅ ( s + 1) ⋅ ( s + 2 ) C( s) k ⋅ (s + 2)
= 2
y = 2
R ( s) ( s + 1) + k ⋅ ( s + 2 ) R ( s) ( s + 1) + k ⋅ ( s + 2 )
Los dos sistemas tienen la misma ecuación característica, lo cual implica que tienen los mismos polos
en lazo cerrado para el mismo valor de k. El trazado del lugar geométrico de las raíces es el mismo,
pero los ceros en lazo cerrado son diferentes debido al efecto de la función de transferencia del

elemento de medida. De este modo, para conocer los ceros en lazo cerrado es necesario observar la
estructura concreta del sistema.
Interesará minimizar el efecto del cero de lazo cerrado introducido por el elemento de medida; para
ello, deberá alejarse substancialmente del eje imaginario del plano S respecto a los polos en lazo
cerrado. En conclusión, los polos del elemento de medida deben tener constantes de tiempo mucho
menores que las características de la planta, esto es, el elemento de medida debe “procesar” muy
rápidamente la información respecto a la respuesta transitoria de la planta.
2.6.3 Aspectos importantes de cons trucción del LGR.
1. Efectos de adición de polos y ceros.
La adición de un polo a la función de transferencia en lazo abierto GLA(s) en el semiplano izquierdo

del plano S tiene el efecto de “empujar” el LGR hacia el semiplano derecho. De este modo,
manteniendo la ganancia del sistema, su dinámica tenderá a empeorar.
Eje imaginario Eje imaginario
2 4
1.5 3
1 2
0.5 1
0 0
-0.5 -1
-1 -2
-1.5 -3
-2 -4
-2 -1.5 -1 -0.5 0 0.5 1 1.5 2 -4 -3 -2 -1 0 1 2 3 4
Eje real Eje Real
Eje Imaginario
6
-2
-4
-6
-6 -4 -2 0 2 4 6
Eje Real
Fig. 2.37 Efecto sobre el LGR de la adición de polos
La adición de un cero a la función de transferencia en lazo abierto GLA(s) en el semiplano izquierdo

del plano S tiene el efecto de “atraer” el LGR hacia el semiplano izquierdo. En este caso, la dinámica
sufre un efecto contrario al anterior.

Eje Imaginario
4
-1
-2
-3
-4
-4 -3 -2 -1 0 1 2 3 4
Eje Real
Fig. 2.38 Efecto sobre el LGR de la adición de un cero.
2. Efectos de movimientos de polos y ceros.
A medida que un polo en lazo abierto se acerca al eje imaginario, las ramas del LGR se acercan más al
semiplano derecho. Observaríamos el efecto contrario en el caso del movimiento de un cero en lazo
abierto.
Eje Imaginario Eje Imaginario

5 5
4 4
3 3
2 2
1 1
0 0
-1 -1
-2 -2
-3 -3
-4 -4
-5 -5
-14 -12 -10 -8 -6 -4 -2 0 -14 -12 -10 -8 -6 -4 -2 0
Eje real Eje Real
Eje Imaginario Eje Imaginario
5 2
4
1.5
3
1
2
0.5
1
0 0
-1 -0.5
-2
-1
-3
-1.5
-4
-5 -2
-6 -5 -4 -3 -2 -1 0 1 -2 -1.5 -1 -0.5 0 0.5 1
Eje Real Eje Real
Fig. 2.39. Efecto sobre el LGR del movimiento de un polo hacia el semiplano derecho.

Ejemplo 2.13
Se desea dibujar el lugar geométrico de las raíces del siguiente sistema:
R(s) (s + 0.2)(s + 1) C(s)

k G (s) =
s ( s + 0.04) ( s + 40)
+ 2 2
-
Calcular:
1.- Lugar geométrico de las raíces sobre el eje real.
2.- Asíntotas, sus puntos de cruce por el eje real y sus puntos de intersección con el eje imaginario.
3.- Puntos de cruce del LGR con el eje imaginario, utilizando para ello el diagrama de Bode de la
figura 2.40, que representa la evaluación de la función de transferencia G(s) sobre el eje
imaginario.
4.- Rango de valores de k para el cual el sistema en lazo cerrado es estable.
5.- Puntos de ruptura del LGR y valores de k asociados a los mismos, utilizando para ello la gráfica de
la figura 2.41, que representa la evaluación sobre el eje real del módulo de la función de
transferencia G(s).
6.- Dibujar el lugar geométrico de las raíces resultante.
Ganancia (dB)
50
-50
-100
-150
-3 -2 -1 0 1 2
10 10 10 10 10 10
Fase (grados) Frecuencia (rad/seg)
-50
-100
-150
-180
-200
-250
-3 -2 -1 0 1 2
10 10 10 10 10 10
Frecuencia (rad/seg)
Fig. 2.40 Diagrama de Bode en módulo y fase de G(s)

10·log|G(s)|
60
40
20
-20
-40
-60
-80
-100
-120
-3 -2 -1 0 1 2
10 10 10 10 10 10
Semieje real negativo plano S
Fig. 2.41 Módulo de G(s) evaluada sobre el eje real
Solución:
1.- LGR sobre el eje real:
Polos LA: sp1 = 0 sp2,3 = -40 (doble) sp4,5 = -0.04 (doble)
Ceros LA: sc1 = -1 sc2 = -0.2
jω
σ
-40 -1 -0.2 -0.04
2.- Asíntotas:
±180 ⋅ ( 2 λ + 1) ±180 ⋅ ( 2λ + 1)  ±60º


θa = = =
n−m 5− 2  180º
σa =
∑ pi − ∑ z j =
−2 ⋅ 40 − 2 ⋅ 0.04 + 1 + 0.2
= −26.3
n−m 5− 2
Punto de intersección de las asíntotas con el eje imaginario: ωa = σa·tg 60º = 45.55

3.- Puntos de cruce con el eje imaginario:
Sabemos que el diagrama de Bode de la figura 2.40 representa el módulo y la fase de la función de
transferencia G(s) evaluada sobre el eje imaginario, por tanto conocemos el valor que toman estas dos
magnitudes para cada punto sobre dicho eje.
Sabemos también que los puntos de cruce del LGR por el eje imaginario son puntos que deben
pertenecer al LGR, y para ello deben cumplir la condición de ángulo:
G(sa) ∈ LGR ⇔ Fase[G(sa)] = ±180º
En consecuencia, para obtener los puntos del LGR que se encuentran sobre el eje imaginario (puntos
de intersección del LGR con el eje jω), únicamente deberemos considerar las frecuencias del diagrama
de Bode para las cuales tengamos ±180º de fase.
Observando el diagrama de Bode de la figura 2.40, comprobamos que las frecuencias a las cuales la
fase de G(s) es de ±180º, son:
ω = 0.055 ω = 0.35 ω = 40
que se corresponderán con los puntos de cruce del LGR con el eje imaginario.
4.- Margen de valores de k para el cual el sistema es estable:
Sabiendo que el límite de estabilidad se da cuando los polos en lazo cerrado se sitúan sobre el eje
imaginario, podemos determinar los valores de k asociados a los puntos de intersección del LGR sobre
dicho eje, obtenidos en el apartado anterior, para poder obtener los valores de k que situarán al sistema
en el límite de la estabilidad.
Para s = ±j0.055 ⇒ k = 1.9592
Para s = ±j0.35 ⇒ k = 162.73
Para s = ±j40 ⇒ k = 127958.55
Conociendo la forma del LGR, podemos decir que:
0 < k <1.9592 ⇒ Sistema estable
1.9592 ≤ k ≤ 162.73 ⇒ Sistema inestable
162.73 < k < 127958.55 ⇒ Sistema estable
k ≥ 127958.55 ⇒ Sistema inestable

5.- Puntos de ruptura del LGR:
dk d ( G (s))
Se puede demostrar que: =0 ⇒ =0
ds ds
1
Para ello, de la ecuación característica: 1 + k ⋅ G (s) = 0 ⇒ k = −
G (s)
d(G (s) )
dk
=0 ⇒
d  1  ds = 0 ⇒ d(G (s) ) = 0
− =
ds ds  G (s)  (G (s) )2 ds
En conclusión: los puntos de ruptura se corresponden con los máximos y mínimos de G(s).
De este modo, para obtener los puntos de ruptura sobre el eje real, podremos hacerlo observando los
máximos y mínimos de la gráfica del módulo de G(s) evaluado sobre el eje real, que tenemos en la
figura 2.41. Ahora bien, deberemos tener cuidado, porque los puntos que se corresponden con los
polos y ceros en lazo abierto provocan valores infinitos en la gráfica, debido a que la escala es
logarítmica. Así, consideraremos únicamente los máximos y mínimos relativos.
De la gráfica obtenemos que los puntos de ruptura son:
σ ≈ 1.5·10-2 ⇒ Punto de ruptura de dispersión.
σ ≈ 2.2 ⇒ Punto de ruptura de confluencia.
σ ≈ 12 ⇒ Punto de ruptura de dispersión.
Nota: Podemos observar, según la figura 2.41, que en la zona del semieje real comprendida
entre σ = 0.2 y σ = 1 existe otro máximo relativo, pero, a diferencia de los demás, éste
no se corresponde con un punto de ruptura, ya que los valores de σ correspondientes a
esta zona no pertenecen al LGR, como podemos comprobar en el apartado 1.
Los valores de ganancia k asociados a dichos puntos de ruptura, calculados mediante la sustitución del
valor de σ obtenido en la ecuación que nos define k, serán los siguientes:
Para s = -1.5·10-2 ⇒ k = 0.08225
Para s = -2.2 ⇒ k = 6110.8
Para s = -12 ⇒ k = 10367.7

6.- Dibujar el LGR:

40
30 k=120840
s=j40
20
k=10375
10 s=-12.44
-10
-20
-30
-40
-50 -40 -30 -20 -10 0 10
Eje Real
Fig. 2.42 Lugar geométrico de las raíces de G(s)
2 0.4
1.5 0.3
k=144.46
s=j0.33
1 0.2
k=2.01
0.5 0.1 s=j0.055
k=6110
s=-2.22
0 0
-0.5 -0.1
k=0.0825
-2
s=-1.41·10
-1 -0.2
-1.5 -0.3
-2 -0.4
-2 -1 0 -0.2 -0.1 0
Eje Real Eje Real
Fig. 2.43 Ampliación de la zona del LGR próxima al origen
2.6.4 LGR en sistemas discretos.
En un sistema discreto pueden evaluarse las raíces de la ecuación característica mediante el LGR
trazado en el plano Z. Ello es posible debido a que la ecuación característica de un sistema discreto
lineal invariante en el tiempo es una función racional de polinomios en Z. Por lo tanto, puede aplicarse

el mismo conjunto de reglas de trazado del LGR que en sistemas analógicos, con la salvedad de que,
además, deben obtenerse los puntos de cruce del LGR con el circulo de radio unidad en el plano Z.
Ejemplo 2.14
Calcular el LGR del siguiente sistema
ZOH PLANTA
R(s) + E*(s) C(s)
Goh(s) G(s)
- T
H(s)
ELEMENTO DE MEDIDA
k
siendo H ( s) = 1; G (s) = ; T = 1 seg.
s ⋅ (s + 1)
1 − e − TS k  0.368k ⋅ ( z + 0.717)
G LA ( z) = Z =
 s s ⋅ ( s + 1)  z 2 − 1368
. z + 0.368
• Polos en: zp1=1; zp2=0.368 y ceros en zc1=-0.717

• Asíntotas: ±180º
• Puntos de ruptura: z1=0.67 ; Kz1=0.196
z2=-2.11; Kz2=15.01
• Corte con el eje imaginario: z3=±j1.161
• Corte con el circulo unitario. Para calcular estos puntos se aplica el criterio de Routh sobre el
plano transformado W, ya que el circulo de radio unidad del plano Z se transforma mediante la
transformada bilineal en el eje imaginario del plano W.
G LA ( w ) = G LA ( z) z = 1+0.5w
1− 0.5w
(1 − 0.038k ) w 2 + (0.924 − 0.386k ) w + 0.924 k

0 = 1 + G LA ( w ) =
w( w + 0.924)
Si a continuación se aplica el criterio de Routh, se encuentra que el sistema será estable para
0<k<2.39. Si se sustituye el valor límite de k en la ecuación característica, se encuentra la
posición del LGR sobre el circulo unitario, z=0.245±j0.97 (modulo=1; fase=±75.8º).

Eje imaginario
2
1.5
0.5
-0.5
-1
-1.5
-2
-3 -2 -1 0 1 2
Eje real
Fig. 2.44 Lugar geométrico de raíces discreto
Por último, debe indicarse que disminuir el periodo de muestreo conlleva modificar el LGR, de
manera que aumenta el margen de valores de k para el cual el sistema es estable, y se produce un
acercamiento de las ramas del LGR a z=1, indicativo de la menor desvirtuación del sistema discreto
frente al sistema continuo equivalente.
Ejemplo 2.15
Si se repite el ejemplo anterior con T=0.1 seg se obtiene el siguiente LGR.
Eje imaginario
2
1.5
0.5
-0.5
-1
-1.5
-2
-3 -2 -1 0 1 2
Eje real
Fig. 2.45 Lugar geométrico de raíces discreto
Donde el rango de valores de k para el cual el sistema es estable es 0<k<20.3

2.7 Problemas.
Problema 1
Dado el sistema de la figura :
R(s) +
k e −sT C(s)
Goh(s)
- T s +1
1 − e − sT
Donde Goh(s) =
s
1.- Obtener la función de transferencia de lazo abierto y la función de transferencia de lazo cerrado del
sistema.
2.- Dibujar el LGR del sistema en función de k sabiendo que el periodo de muestreo es T.
3.- Calcular la relación entre el periodo de muestreo T y la ganancia k, que proporciona el límite de
estabilidad. Dibujar la curva que relaciona k y T, calculando para ello los valores en los puntos T=
0, 0.1, 0.5, 1, 2 y 10. Razonar el resultado obtenido y compararlo con el sistema continuo
equivalente.
Solución:
− sT

1.- G ( z) = TZ1 − e e − sT  −1 −1  1   1  1 1  z z
⋅  = (1 − z ) ⋅ z ⋅ TZ ; TZ  = TZ − = −
 s s + 1   s ⋅ (s + 1)   s ⋅ ( s + 1)   s s + 1 z − 1 z − e−T
La función de transferencia de lazo abierto es :
1 − e−T
FTLA ( z) = k ⋅ G ( z) = k ⋅
z( z − e − T )
Y la función de transferencia de lazo cerrado :
k (1 − e − T )
FTLC ( z) = 2
z − z ⋅ e − T + k (1 − e − T )
2.- En lazo abierto hay dos polos y teniendo en cuenta que uno de los polos tiende a z=1 cuando T
tiende a cero, y tiende a z=0 cuando T tiende a infinito, habrá por tanto tres posibles LGR
dependiendo del valor de T.

T→0 T→∞ T
0.5 e-T e-T

2
3.- La ganancia límite de estabilidad corresponde a la situación de las raíces de lazo cerrado cuando
estas cruzan el círculo de radio unidad. La forma más sencilla de calcular esta ganancia es aplicar la
condición de módulo :
1
1
z ⋅ z − e−T
k ⋅ G ( z) = 1 k=
1 − e −T
e-T e-T
2
De la gráfica se observa que z = 1; z − e−T = 1 .
1
Así que : k=
1 − e−T
La gráfica que se obtiene de esta relación es :
10
9
8
k 5
0
0 0.5 1 1.5 2 2.5 3 3.5 4
T
De la gráfica se observa que, a medida que el periodo de muestreo aumenta, el margen de estabilidad
disminuye.

Problema 2
Determinar el lugar geométrico de las raíces de un sistema cuya ganancia de lazo viene dada por:
G (s) =
(
k ⋅ s 2 + 0.1056 )
(
s ⋅ s +12
)
Especificar para ello:
* Polos y ceros en lazo abierto.
* Asíntotas y su cruce por el eje real.
* Puntos de corte con el eje jω.
* Ángulos de salida y llegada de polos y ceros.
* Puntos de ruptura. Especificar si son de confluencia o dispersión.
Solución:
1.- Polos y ceros en lazo abierto:
Polos: 0, ±j
Ceros: ±0.325j
2.- LGR en el eje real:
jω
j
0.325j
σ
-0.325j
-j

3.- Asíntotas:
±180 ⋅ ( 2 λ + 1)
θa =
n−m
n = número de polos finitos ⇒ n = 3
m = número de ceros finitos ⇒ m = 2
±180 ⋅ ( 2λ + 1)  λ = 0 θ a = 180º

θa = ⇒
1  λ = 1 θ a = 540º = 180º
La única asíntota es el eje real.
4.- Puntos de corte con el eje imaginario:
1+
(
k ⋅ s 2 + 0.1056 ) = 0; s 3 + k ⋅ s 2 + s + 0.1056 ⋅ k = 0
(
s ⋅ s +1
2
)
Algoritmo de Routh:
s3 1 1
s 2
k 0.1056 ⋅ k
1
s 0.8944
s 0.1056 ⋅ k
0
No existe ningún valor k > 0 que anule una fila. Los únicos puntos de corte con el eje imaginario se
darán para k = 0, es decir, son los polos en lazo abierto ya conocidos.
5.- Ángulos de salida y entrada de polos y ceros:

jω
φ1
j
θ1
0.325j
φ2
σ
θ2
-0.325j
-j φ3
Para conocer los ángulos de arranque y llegada de los polos y ceros deberemos aplicar la condición de
ángulo del LGR.

θ1 + θ2 - ( φ1 + φ2 + φ3 ) = ±180º
Para φ1: 90º + 90º - φ1 - 90º - 90º = ±180º ⇒ φ1 = 180º
Para φ2: -90º + 90º + 90º -φ2 - 90º = ±180º ⇒ φ2 = 180º
Para θ1: θ1 + 90º + 90º -90º -90º = ±180º ⇒ θ1 = 180º
Como que el LGR es simétrico respecto al eje real, podemos decir que:
θ2 = 180º φ3 = 180º
6.- Determinación de los puntos de ruptura:
1+
(
k ⋅ s 2 + 0.1056 ) = 0; k=
− s3 − s
(
s ⋅ s +12
) (s 2
+ 0.1056 )
 s = ±0.4871
dk 
=0 ⇒ 
ds  s = ±0.6677
Los puntos s = 0.4871 y s = 0.6677 no pertenecen al LGR, por tanto no pueden ser puntos de ruptura.
Los únicos posibles puntos de ruptura son s = -0.4871 y s = -0.6677.
Comprobaremos si en esos puntos k > 0.
k s=−0.4871 = 17577
. >0
k s=−0.6677 = 1.7569 > 0
Por tanto tendremos dos puntos de ruptura en el eje real situados en s = -0.4871 y en s = -0.6677 .
Para determinar si son de confluencia o de dispersión, lo podemos hacer de dos maneras:
d2k
* Buscar la segunda derivada y ver si los puntos son máximos o mínimos.
ds 2
* Sabiendo que las ramas salen de los polos con ángulos de 180º y que llegan a los ceros con
ángulos de 180º, y que además, estas dos ramas no pueden cruzarse, ya que provocarían
puntos de ruptura adicionales inexistentes, la única solución posible es que s = -0.6677 sea
de confluencia y s = -0.4871 sea de dispersión.
Por tanto el gráfico del LGR quedará de la forma:

1.5
0.5
Eje Imaginario
0
-0.5
-1
-1.5
-1.5 -1 -0.5 0 0.5 1 1.5
Eje Real
Eje Real
Problema 3
Determinar el lugar geométrico de las raíces de un sistema cuya ganancia de lazo viene dada por la
expresión:
k
G (s) =
(s 2
)(
+ 2s + 2 s 2 + 2s + 5 )
Especificar para ello:
* Polos y ceros en lazo abierto.
* Asíntotas y su cruce por el eje real.
* Puntos de corte con el eje jω.
* Ángulos de salida de los polos.

dk
* Puntos de ruptura sabiendo que una raíz de la ecuación = 0 es s = -1.
ds
Solución:
1.- Polos y ceros en lazo abierto:
sp1,2 = -1 ± j sp3,4 = -1 ± 2j
sc = ∞ (cuarto orden)
2.- LGR sobre el eje real:

Como en el eje real no existe ninguna singularidad, el eje real no puede pertenecer al LGR.
3.- Asíntotas:
±180 ⋅ ( 2 λ + 1) ±180 ⋅ ( 2λ + 1)
θa = = = ±45º
n−m 4−0
σa =
∑ i −∑zj
p
=
( −1 + j) + ( −1 − j) + ( −1 + 2 j) + ( −1 − 2 j)
= −1
n−m 4−0
4.- Puntos de cruce con el eje imaginario:
1 + G(s) = 0 ⇒ s4 + 4s3 + 11s2 + 14s + 10 +k = 0
Aplicando el algoritmo de Routh y anulando la fila s1 , obtenemos el valor de k que lo permite:
k = 16.25
Para este valor de k, y a partir del polinomio auxiliar, conseguimos los puntos de cruce con jω.
s1,2 = ±j1.87
5.- Ángulos de arranque:
Aplicando la condición de ángulo sobre cada uno de los polos complejos conjugados en LA, se
obtiene que los ángulos de salida de cada uno de ellos son:
θs1 = 90º θs2 = -90º θs3 = -90º θs4 = 90º
6.- Puntos de ruptura:
De la ecuación característica 1 + G(s) = 0, obtenemos

dk
ds
(
= − 4s 3 + 12s 2 + 22s + 14 = 0 .)
Sabiendo que s = -1 es una raíz de la ecuación anterior, podemos aplicar Rufini para calcular las otras
dos raíces:
dk
= ( s + 1)( s + 1 + j158
. )( s + 1 − j158
. )=0
ds
Para comprobar si las raíces de la ecuación corresponden o no a puntos de ruptura, deberemos

verificar en primer lugar si pertenecen o no al LGR.
Evidentemente s = -1 no puede ser un punto de ruptura porque es un punto que no pertenece al LGR,
ya que no existe LGR sobre el eje real, como hemos comprobado anteriormente.

En cuanto a las raíces complejas conjugadas, podemos comprobar su pertenencia al LGR mediante la
condición de ángulo.
jω
θ3
2j
-1+j1.58
θ1
j
-1 σ
θ2
-j
-1-j1.58
θ4 -2j
Respecto al punto s = -1 +j1.58: -(θ1 + θ2 + θ3 + θ4) = -(90º + 90º -90º + 90º) = -180º. Es decir, que
cumple la condición de ángulo y por tanto pertenece al LGR. Por simetría, su complejo conjugado
también pertenecerá.
Podemos comprobar que los valores de k asociados a estos puntos son positivos, en consecuencia son
efectivamente puntos de ruptura.
El lugar geométrico de las raíces quedará entonces de la forma siguiente:
1
Eje Imaginario
-1
-2
-3
-4
-4 -3 -2 -1 0 1 2
E je R e a l
Eje Real

Problema 4
Dado el sistema de la figura siguiente:
H2(s)
H1(s)
R(s)
K 1 − e − Ts 1 C(s)
+
- T T s s+1
Se pretende observar cómo depende la estabilidad del sistema de la elección del periodo de muestreo
T. Para ello, encontrar una relación entre T y el valor de K que proporciona el límite de estabilidad.
Solución:
C ( z) H1 ( z ) H 2 ( z )
H ( z) = =
R ( z ) 1 + H1 ( z ) H 2 ( z )
H1 ( z) = K
− Ts
 1 − e 1   1 
H 2 ( z) = TZ
 s
⋅ ( −1
 = 1 − z ⋅ TZ
s + 1
) −1
(
1
 = 1 − z ⋅ TZ −
 s( s + 1) 
)
1 

 s s + 1
z−1  z z  1 − e− T
H 2 ( z) = ⋅ − =
z  z − 1 z − e− T  z − e− T
1 − e− T
H ( z) =
K
z − e− T =
K 1 − e− T ( )
1 − e− T z − e− T (1 + K) + K
1+ K
z − e− T
La ecuación característica vendrá determinada por el denominador de la función de transferencia:
z − e− T (1 + K) + K = 0
1 − e− T 1
−T
=−
z− e K
Dibujaremos el LGR para observar la dependencia de la estabilidad del sistema con K.
La función de transferencia en lazo abierto viene dada por la expresión: H1 ( z) H 2 ( z) =

(
K 1 − e− T )
z − e− T

Podemos observar que sólo existe un polo en lazo abierto situado en z=e-T. Este polo será el punto
inicial del LGR. En consecuencia, el LGR tendrá la forma siguiente:
Im(z)
Kosc
1
-T
e Re(z)
Para T=1, el polo en lazo abierto está situado en z = e-1 = 0.368. Para este valor de T, el sistema deja
de ser estable cuando la rama del LGR cruza el círculo de radio unidad:
z = e −1 (1 + K) − K = −1 ⇒ K osc = 2.2
Ahora trataremos de encontrar una relación entre el periodo de muestreo T y el factor de ganancia K
que hace oscilar el sistema.
e − T (1 + K osc ) − K osc = −1
1+ e −T
K osc =
1− e −T
Este relación puede expresarse gráficamente como se puede observar en la figura:
Kosc
4
2.1
1.3
0.5 1 2 3 4 T
Con lo cual podemos comprobar que, a medida que T aumenta, el rango de estabilidad del sistema
disminuye.

Problema 5
Dado el diagrama de bloques de un sistema discreto como el de la figura:
R(s) C(s)
D(z) ZOH H(s)
+ T T
-
donde: T = 2 ⋅ ln 2 seg. D ( z) =
( )
2 3 − z −1
H ( s) =
k
3(1 − z )−1 2s + 1
C ( z)
1.- Determinar, para k=1.5, la función de transferencia .
R ( z)
2.- Dibujar el lugar geométrico de las raíces del sistema, para k desde 0 a ∞, especificando para ello:
* LGR sobre el eje real.
* Asíntotas.
* Puntos de ruptura.
* Puntos de corte con el círculo unitario.
* Puntos de corte con el eje imaginario.
3.- ¿Para qué valores de k se vuelve inestable el sistema?
Nota: La ecuación característica en el plano transformado W tiene la forma:
−0.44444 ⋅ ( w + 0.721347) ⋅ ( w − 1442695

. )⋅k
1+ =0
w ⋅ ( w + 0.4808984 )
Solución:
1.-
C ( z) D ( z) Z oh H ( z)
=
R ( z) 1 + D ( z) Z oh H ( z)

 15 . 
 2s + 1  z − 1  15  z − 1 3
(
Z oh H ( z) = 1 − z −1
) ⋅ TZ
 s
=
 z
⋅ TZ
 (
s 2
.
s + 1) 
=
z
⋅ ⋅
4 z − 1 ⋅
z
( ) ( z − 0.5)
 
2( 3z − 1) z z−1
⋅ ⋅
C ( z) 3( z − 1) 4( z − 1) ⋅ ( z − 0.5) z 3z − 1
= =
R ( z) 2( 3z − 1) z z−1 2z 2
1+ ⋅ ⋅
3( z − 1) 4( z − 1) ⋅ ( z − 0.5) z
2.-
z−1  k  z − 1  0.5  0.5 ⋅ k

TZ{Z oh H ( s)} = ⋅ TZ = ⋅ k ⋅ TZ =
z  s( 2s + 1)  z  s( s + 0.5)  z − 0.5
D ( z) ⋅ TZ{Z oh H ( s)} =
2(3z − 1) 0.5 ⋅ k
⋅ =
k ⋅ z − 13 ( )
3( z − 1) z − 0.5 ( z − 1) ⋅ ( z − 0.5)
* LGR sobre el eje real:
Im(z)
1/3 1/2 1 Re(z)
* Asíntotas:
±180 ⋅ ( 2 λ + 1) ±180 ⋅ ( 2λ + 1)  λ = 0 ⇒ ± 180º


θa = = =
n−m 2 −1  λ = 1 ⇒ ± 540º = ±180º

* Puntos de ruptura:
 ( z − 1) ⋅ ( z − 0.5)   z2 − 3 z + 1 
k = −  = − 2 2

 z − (
1
3

 ) 
 z − 1
3


dk
=
( )( )
2 z − 3 2 ⋅ z − 13 − z 2 − 3 2 z + 12
=
(z2 − 23 z
=0
)
dz z − 13
2
( ) z − 13
2
( )
 z=0 ∈ LGR
2 
z −2 3z=0 ⇒ 
 z = 0.6666 ∈ LGR
k z= 0 = 15
. ≥ 0 ⇒ Es punto de ruptura
k z= 0.6666 = 0.16666 ≥ 0 ⇒ Es punto de ruptura
* Puntos de corte con el círculo unitario:
Deberemos aplicar el criterio de Routh modificado.
Podemos observar, según el LGR sobre el eje real, que cortará al círculo unitario en z = 1 y en z = -1,
pero debemos saber con qué valor de k cortará en estos puntos, y si existen o no otros puntos de corte.
0 = 1 + D( z) ⋅ TZ{Z oh H (s)} = 1 +
(
k ⋅ z − 13 )
T
( z − 1) ⋅ ( z − 0.5) z= 1+ 2 w
T
1− w
2
De la nota del enunciado sabemos que la transformada bilineal de la ecuación característica es:
−0.44444 ⋅ ( w + 0.721347) ⋅ ( w − 1442695

. )⋅k
1+ =0
w ⋅ ( w + 0.4808984 )
(1 − 0.4444 k ) ⋅ w 2 + (0.4808984 + 0.32059 k ) ⋅ w + 0.462521⋅ k = 0
w2 1 − 0.4444 ⋅ k 0.462521 ⋅ k
w1 0.4808984 + 0.32059 ⋅ k
w0 0.462521 ⋅ k

0.462521·k > 0 ⇒ k > 0
0.4808984 + 0.32059·k > 0 ⇒ k > -1.5
1-0.4444·k > 0 ⇒ k < 2.25
Margen de estabilidad: 0 < k < 2.25
Sustituyendo en la ecuación característica:
1+
(
k ⋅ z − 13 ) =0 ⇒
 z 1 = 0.25

z 2 + 0.75z − 0.25 = 0 ⇒ 
( z − 1) ⋅ ( z − 0.5) k =2.25  z 2 = −1
El sistema deja de ser estable cuando el LGR corta el círculo unitario. Para el caso que nos ocupa, esto
se da en z = -1 para un valor de k = 2.25.
Para el otro posible valor límite de k (k = 0):
1+
(
k ⋅ z − 13 ) =0 ⇒ 2
 z1 = 1

z − 2 z + 0.5 = 0 ⇒ 
3
( z − 1) ⋅ ( z − 0.5) k =2.25  z 2 = 0.5
Es decir, para k = 0, estaremos situados en z = 1 y z = 0.5, siendo z = 1 otro de los márgenes de

estabilidad. En consecuencia, el LGR no cortará al círculo unitario en otros puntos que no sean z = 1 y
z = -1.
* Puntos de corte con el eje imaginario:
Para ello aplicaremos el criterio de estabilidad de Routh directamente al polinomio resultante de la

ecuación característica 1+G(z) = 0, igualando una fila a 0 para hallar las raíces imaginarias.
1+
(
k ⋅ z − 13 ) =0
( z − 1) ⋅ ( z − 0.5)
( )
z 2 + k − 3 2 ⋅ z + 0.5 −
k
3
=0
k
z2 1 0.5 −
3
z1 k − 3 2
k
z 0 0.5 −
3

k
0.5 − = 0 ⇒ k = 15
.
3
k − 3 2 = 0 ⇒ k = 15.
Para k = 1.5 en la fila z1: ( )

z 2 − 0.5 − 1.5 3 = 0 ⇒ z 2 = 0 ⇒ z = 0
1
0.8
0.6
0.4
Eje Imaginario
0.2
k = 2.25
0
-0.2
-0.4
-0.6
-0.8
-1
-1.5 -1 -0.5 0 0.5 1 1.5

Eje Real
Eje Real
Problema 6
Dado el sistema de la figura 1:
R(s) k C(s)
+ 2
s + 2s + 3.25
-
1+ b ⋅ s
Figura 1
1.- Diseñar k y b para tener un error estacionario de posición del 10%, así como un coeficiente de
amortiguamiento de los polos en lazo cerrado de 0.707.
2.- Manteniendo el valor de k calculado anteriormente, dibujar el lugar geométrico de las raíces en
función de b. Para ello calcular:
* LGR sobre el eje real.
* Ángulos de las asíntotas y su punto de intersección con el eje real.
* Puntos de ruptura.
* Ángulos de arranque de los polos y llegada a los ceros.

3.- En la figura 2 se muestra una topología alternativa a la anterior, que posee los mismos polos en
lazo cerrado para los mismos valores de k y b.
R(s) k C(s)
+ 1+ b ⋅ s 2
- s + 2s + 3.25
Figura 2
Relacionar cada una de las topologías presentadas con las respuestas temporales en lazo cerrado al
escalón que se pueden ver en las gráficas siguientes. Razonar la respuesta.
1.2 1.2
1 1
0.8 0.8
0.6 0.6
0.4 0.4
0.2 0.2
0 0
0 1 2 0 1 2
Tiempo (seg) Tiempo (seg)
Gráfica 1 Gráfica 2
Solución:
1.- Valor de k y b.
La función de transferencia en lazo abierto tiene la forma:
k ⋅ (1 + b ⋅ s)
G LA (s) =
s 2 + 2s + 3.25
Se desea un error estacionario de posición del 10%, luego:
1
e ssp = = 01
. ⇒ Kp = 9
1+ K p
k ⋅ (1 + b ⋅ s) k
K p = lim G LA (s) = lim = ⇒ k = 3.25 ⋅ 9 = 29.25
s→0 s→0 s 2 + 2s + 3.25 3.25

Por otra parte, se exige un coeficiente de amortiguamiento de los polos en lazo cerrado de 0.707.
k ⋅ (1 + b ⋅ s)
La ecuación característica: 1 + G LA (s) = 0 ⇒ 1 + =0
s 2 + 2s + 3.25
s 2 + ( 2 + 29.25 ⋅ b ) ⋅ s + 32.5 = 0
ω n 2 = 32.5 ⇒ ω n = 5.7 rad / s
2 ξω n = 2 + 29.25 ⋅ b ⇒ b = 0.207
2.- Lugar geométrico de las raíces en función de b.
Arreglando la ecuación característica:
29.25 ⋅ b ⋅ s
1+ =0
2
s + 2s + 32.5
Raíces:
s1,2 = -1 ± j5.6125 ⇒ Polos complejos conjugados.
s = 0 ⇒ Cero en origen.
* LGR sobre el eje real: Todo el semieje real negativo.
* Ángulos de las asíntotas y su punto de intersección con el eje real:
Tenemos el eje real como única asíntota, en consecuencia no es necesario calcular el punto de
intersección de la misma.
* Puntos de ruptura:
s 2 + 2s + 3.25
b=−
29.25 ⋅ s
db
= 0 ⇒ s 2 − 32.5 = 0 ⇒ s1,2 = ±5.7
ds
s = 5.7 ∉ LGR ⇒ No es punto de ruptura.
s = -5.7 ∈ LGR; b s=−5.7 = 0.3214 ≥ 0 ⇒ s = -5.7 es punto de ruptura.
* Ángulos de arranque de los polos y llegada a los ceros.
- Ángulo de llegada al cero en origen: 180º
- Ángulo de arranque de los polos complejos conjugados:

5.6125
180º − arctg − 90º −θ = 180º
1
θ = −159.8974º
El lugar geométrico de las raíces quedará de la forma:

10
4
Eje Imaginario
-2
-4
-6
-8
-10
-10 -5 0 5 10
Eje Real
Eje Real
3.- Identificación de las respuestas con las topologías correspondientes.
Para ello deberán observarse las funciones de transferencia en lazo cerrado:
k 29.25
Fig. 1 ⇒ G LC (s) = 2 = 2
s + 2s + 3.25 + k ⋅ (1 + b ⋅ s) s + 8.054s + 32.5
k ⋅ (1 + b ⋅ s) 29.25 ⋅ (1 + 0.207s)
Fig. 2 ⇒ G LC (s) = 2 = 2
s + 2s + 3.25 + k ⋅ (1 + b ⋅ s) s + 8.054s + 32.5
Aplicando la expresión del máximo sobreimpulso para un sistema de segundo orden subamortiguado:
− πξ
1− ξ 2
Mp = e = 4.32% ⇒ la gráfica 1 corresponde con la topología de la figura 1.
Problema 7
Determinar los parámetros k, a y b del sistema mostrado en la figura siguiente, conociendo que,
cuando el sistema tiene una rampa de entrada, la señal de salida sigue la señal de entrada con un error
en estado estacionario finito, y que, cuando la ganancia se dobla a 2k, la señal de salida a un impulso
de entrada es una señal senoidal con un periodo de 0.314 segundos.

R(s) k C(s)
+
- ( s + 40)(s + a )(s + b)
Solución:
Error finito a la rampa ⇒ essp = 0 ⇒ kp = ∞
k k
k p = lim G (s) = lim = =∞ ⇒ aob=0
s→0 s→0 ( s + 40)( s + a )( s + b) 40 ⋅ a ⋅ b
k
1+ =0
(s + 40)(s + a )(s + b)
(s + 40)(s + a)(s + b) + k = 0
Suponiendo b = 0:
(s + 40)(s + a)s + k = 0 ⇒ s 3 + ( 40 + a )s 2 + 40as + k = 0
Aplicando el criterio de estabilidad de Routh:
s3 1 40a
s2 40 + a k
(40 + a ) ⋅ 40a − k
s1
40 + a
s0 k
Si la ganancia pasa de k a 2k, el sistema es oscilatorio ⇒ Raíces sobre eje jω ⇒ Fila del algoritmo de
Routh se hace cero:
(40+a)·40a-2k = 0
Las raíces sobre el eje jω se corresponden con las raíces del polinomio auxiliar:
Pa(s) = (40+a)s2 + 2k = 0

2k 2π 2π
s1,2 = ± j = ± jω d = ± = ±j = ± j20 ⇒ 2 k = 400 ⋅ ( 40 + a )
40 + a Td 0.314
Sustituyendo:
 10

(40+a)·40a-2k = 0 ⇒ (40+a)·40a-400·(40+a) = 0 ⇒ a = 
−40 ⇒ Sist. inestable
Así: a= 10 ⇒ k = 200·(40+a) = 200·50 = 10000
10000
G (s) =
s( s + 10)( s + 40)

EE00302C

Cargado por

Copyright:

Formatos disponibles

EE00302C

Cargado por

Información del documento

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

EE00302C

Cargado por

Copyright:

Formatos disponibles

2.

Análisis temporal de sistemas continuos y discretos 53

2. Análisis temporal de sistemas continuos y discretos

2.1 Respuesta temporal.

La respuesta temporal de un sistema lineal invariante en el tiempo puede descomponerse en dos

donde: ct(t) = Respuesta transitoria.

css(t) = Respuesta estacionaria.

La respuesta estacionaria depende fundamentalmente de la señal de excitación al sistema y, si el

© los autores, 2000; © Edicions UPC, 2000.

2.1.1 Sistemas de primer orden.

Fig. 2.1 Sistema de primer orden. T: constante de tiempo del sistema

Fig. 2.2 Respuesta al escalón

© los autores, 2000; © Edicions UPC, 2000.

2.1.2 Sistemas de segundo orden.

ωn= Frecuencia natural no amortiguada.

Los polos de un sistema de segundo orden vienen determinados por la expresión:

* Tipos de sistemas en función del valor de ξ:

En la siguiente gráfica se visualiza la ubicación de polos de la función de transferencia en función del

Fig. 2.3 Diagrama de la ubicación de los polos

Para ξ = 0 ⇒ Sistema oscilatorio. s1,2 = ± jωn

© los autores, 2000; © Edicions UPC, 2000.

Para ξ = 1 ⇒ Sistema con amortiguamiento crítico . s1, 2 = −ωn

* Respuesta al escalón en un sistema subamortiguado (0<ξ<1):

© los autores, 2000; © Edicions UPC, 2000.

2.2 Especificaciones de respue sta transitoria.

Generalmente, en la práctica, se especifican las características o especificaciones requeridas a un

* td: Tiempo de retardo.

Tiempo que tarda la respuesta en alcanzar el 50 % del valor final.

© los autores, 2000; © Edicions UPC, 2000.

* tr: Tiempo de subida o crecimiento.

* tp: Tiempo de pico.

* Mp: Máximo sobreimpulso (porcentual).

El máximo sobreimpulso es un parámetro indicativo de la 'estabilidad relativa' del sistema.

* ts: Tiempo de establecimiento.

2.2.1 Particularización para sistemas de segundo orden subamortiguados.

© los autores, 2000; © Edicions UPC, 2000.

3- Cálculo del máximo sobreimpulso:

De hecho, cuanto mayor es ξ, menor es el máximo sobreimpulso obtenido en la respuesta transitoria.

© los autores, 2000; © Edicions UPC, 2000.

4- Cálculo del tiempo de establecimiento:

Para determinar las expresiones del tiempo de establecimiento en un sistema subamortiguado, se

Así, considerando los errores:

© los autores, 2000; © Edicions UPC, 2000.

Fig 2.8 Respuestas al escalón de sistemas de segundo orden

Fig. 2.9 Amplificador magnético en cascada con un filtro paso-bajo y preamplificador

© los autores, 2000; © Edicions UPC, 2000.

Si la constante de tiempo del amplificador magnético es igual a 1 seg. (τ = 1 seg.) y la ganancia es

Obteniendo, entonces, el diagrama de bloques siguiente :

Vin(s) E(s) K V1(s) 1 V0(s)

© los autores, 2000; © Edicions UPC, 2000.

Escogeremos por lo tanto: C = 5.11·10-4 F

El tiempo de estabilización quedará entonces como:

2.2.2 Sistemas de orden superior.

Desarrollando en fracciones parciales y aplicando la antitransformada de Laplace:

© los autores, 2000; © Edicions UPC, 2000.

en la segunda expresión no existen los términos (s+zi) en numerador ni (s+pi) en denominador. En

2.3 Respuesta transitoria de si stemas discretos.

2.3.1 Sistema de control discreto en lazo cerrado.

El control discreto es caracterizado a través de su función de transferencia en Z y es diseñado para

© los autores, 2000; © Edicions UPC, 2000.