Automática Tema 4. Análisis de la Respuesta Temporal de Sistemas

Departamento de Ingeniería
de Sistemas y Automática
Tema 4.
Análisis de la Respuesta Temporal de
Sistemas LTI
Automática
2º Curso del Grado en
Ingeniería en Tecnología Industrial
Contenido
Tema 4.- Análisis de la respuesta temporal de
sistemas LTI
4.1. Introducción.
4.2. Análisis de la respuesta transitoria de un sistema:
Departamento de Ingeniería
de Sistemas y Automática
4.2.1. Sistemas de primer orden.
4.2.2. Sistemas de segundo orden.
4.2.3. Sistemas de orden superior.
4.3. Introducción a la identificación de sistemas.
4.4. Análisis de la estabilidad de un sistema.
4.5. Análisis de la respuesta estacionaria de un sistema.
Introducción
Señales de prueba:
Departamento de Ingeniería
de Sistemas y Automática
– En el análisis y diseño de sistemas de control es necesario
tener una base para comparar los sistemas de control. Esto
se hace especificando las señales de entrada de prueba y
comparando las respuestas de varios sistemas a estas
señales de entrada.
– Las señales de prueba que se usan regularmente son
funciones escalón, rampa, parábola, impulso, senoidal, etc.
Introducción
Respuesta transitoria y respuesta estacionaria
La respuesta en el tiempo de un sistema de control consta de
dos partes: la respuesta transitoria y la respuesta estacionaria.
Departamento de Ingeniería
de Sistemas y Automática
La respuesta transitoria se
refiere a la que va del
estado inicial al estado final.
La respuesta estacionaria se
refiere a la manera en la
cual se comporta la salida
del sistema conforme el
tiempo tiende a infinito.
Introducción
• Estabilidad absoluta y error estacionario:
Departamento de Ingeniería
de Sistemas y Automática
– Un sistema de control está en equilibrio si, en ausencia de
cualquier perturbación o entrada, la salida permanece en el
mismo estado.
– Un sistema de control lineal e invariante con el tiempo es
estable si la salida termina por regresar a su estado de
equilibrio cuando el sistema está sujeto a una condición
inicial.
– Si la salida de un sistema en estado estable no coincide
exactamente con la referencia, se dice que el sistema tiene
un error estacionario (o en régimen permanente –erp-). El erp
indica la precisión del sistema.
– Al analizar un sistema de control se debe examinar el
comportamiento de la respuesta transitoria y el
comportamiento del error estacionario.
Contenido
Tema 4.- Análisis de la respuesta temporal de
sistemas LTI
4.1. Introducción.
4.2. Análisis de la respuesta transitoria de un sistema:
Departamento de Ingeniería
de Sistemas y Automática
4.2.1. Sistemas de primer orden.
4.2.2. Sistemas de segundo orden.
4.2.3. Sistemas de orden superior.
4.3. Introducción a la identificación de sistemas.
4.4. Análisis de la estabilidad de un sistema.
4.5. Análisis de la respuesta estacionaria de un sistema.
Análisis de la respuesta transitoria
Respuesta Forzada y Natural:
– Sistema continuo representado por la ecuación diferencial
con salida y(t) y entrada u(t)
Departamento de Ingeniería
de Sistemas y Automática
a n y n )  a n 1 y n 1)    a1 y  a 0 y b m u m )   b1 u b 0 u
con un conjunto de condiciones iniciales y(0), y'(0),, y n1) (0)
No nulas !!!
siendo n el orden del sistema.
– La obtención de la respuesta del sistema y(t) ante entrada u(t)
se realiza por aplicación de la transformada de Laplace.
a n ( s n Y ( s)  s n 1 y (0)  s n  2 y ' (0)  )  a n 1 ( s n 1Y ( s)  s n  2 y (0)  )  a 0Y ( s) 
b m ( s mU ( s)  s m1u(0)  ) b m1 ( s m1U ( s)  s m 2 u(0)  ) b0U ( s)
Análisis de la respuesta transitoria
Respuesta Forzada y Natural (cont.):
– Reagrupando términos:
(an sn  an1sn1 a1s  a0 )Y ( s)  (bm sm b1s  b0 )U ( s)  P( s)
Departamento de Ingeniería
de Sistemas y Automática
con P(s) polinomio que depende de las condiciones iniciales.
– La transformada de la respuesta Y(s) de un sistema continuo
se puede expresar:
Y (s) 
bm s m    b1s  b0
a n s n  a n 1s n 1    a1s  a0
Respuesta forzada
U (s) 
P( s)
a n s n  a n 1s n 1    a1s  a0
Respuesta natural
Análisis de la respuesta transitoria
Respuesta Forzada y Natural (cont.):
– Se particularizará el calculo de la respuesta transitoria para
sistemas de orden 1º, 2º y superior.
Departamento de Ingeniería
de Sistemas y Automática
– Sistema caracterizado por la respuesta forzada, asumiendo
respuesta natural nula.
U(s)
G (s)
Y ( s )  G ( s )U ( s )
Y(s)
Análisis de la respuesta transitoria
Sistemas de primer orden:
– Un sistema de primer orden (SPO) queda descrito por una
ecuación diferencial de la forma
y (t )  a0 y (t )  b0u (t )
Departamento de Ingeniería
de Sistemas y Automática
con función de transferencia
b0
G( s) 
s  a0
– La respuesta escalón de amplitud A será:
b0 A
Y ( s) 
s( s  a0 )
Análisis de la respuesta transitoria
Sistemas de primer orden (cont.):
– Descomponiendo en fracciones simples:
Departamento de Ingeniería
de Sistemas y Automática
K1
K2
b0 A 1 b0 A 1
Y ( s) 



s s  a0
a0 s a0 s  a0
– Aplicando la transformada inversa:
b0 A
y (t ) 
(1  e a0t )ue (t )
a0
respuesta de tipo exponencial.
Análisis de la respuesta transitoria
Sistemas de primer orden (cont.):
b0
1
, la constante de tiempo T 
a0
a0
como parámetros específicos de un SPO.
Departamento de Ingeniería
de Sistemas y Automática
– Se definen la ganancia K 
– Forma estándar del SPO:
b0 1
K
G (s) 

s
a0
 1 Ts  1
a0
Análisis de la respuesta transitoria
Sistemas de primer orden (cont.):
Departamento de Ingeniería
de Sistemas y Automática
– Respuesta ante una entrada escalón unitario:
Análisis de la respuesta transitoria
Sistemas de primer orden (cont.):
Departamento de Ingeniería
de Sistemas y Automática
– Respuesta ante una entrada impulso unitario:
Análisis de la respuesta transitoria
Sistemas de primer orden (cont.):
– Respuesta ante una entrada rampa pendiente unitaria:
Departamento de Ingeniería
de Sistemas y Automática
Error en régimen
permanente (erp)
Análisis de la respuesta transitoria
Sistemas de segundo orden:
– Sistema de segundo orden (SSO) queda descrito por una
ecuación diferencial:
y  a1 y  a0 y  b0 u
Departamento de Ingeniería
de Sistemas y Automática
con función de transferencia:
G (s) 
b0
s 2  a1s  a0
– La respuesta escalón de amplitud A será:
Y (s) 
b0 A
s ( s 2  a1s  a0 )
Análisis de la respuesta transitoria
Sistemas de segundo orden (cont.):
– La respuesta depende de las raíces del denominador
( s 2  a1s  a0 )  ( s  s1 )( s  s 2 )
Departamento de Ingeniería
de Sistemas y Automática
– Casos particulares:
1. Raíces reales distintas.
2. Raíces reales repetidas.
3. Raíces complejas conjugadas.
Análisis de la respuesta transitoria
Sistemas de segundo orden (cont.):
Caso 1: Raíces reales distintas.
K1
K2
K3
Y ( s) 


s s  s1 s  s2
Departamento de Ingeniería
de Sistemas y Automática
Aplicando la transformada inversa de Laplace:
y (t )  ( K1  K 2e  s1t  K 3e  s2t )ue (t )
Se denominan sistemas sobreamortiguados.
Análisis de la respuesta transitoria
Sistemas de segundo orden (cont.):
Departamento de Ingeniería
de Sistemas y Automática
Caso 1: Raíces reales distintas.
La rapidez de respuesta depende de la colocación de los polos.
Análisis de la respuesta transitoria
Sistemas de segundo orden (cont.):
Caso 2: Raíces reales repetidas.
Y (s) 
K1 K 21
K 22


s s  s1 ( s  s1 ) 2
Departamento de Ingeniería
de Sistemas y Automática
Aplicando la transformada inversa:
y (t )  ( K1  K 21e  s1t  K 22te  s1t )ue (t )
Se denominan sistemas crítico-amortiguados.
Análisis de la respuesta transitoria
Sistemas de segundo orden (cont.):
Departamento de Ingeniería
de Sistemas y Automática
Caso 2: Raíces reales repetidas.
La rapidez de respuesta depende de la colocación del polo doble.
Análisis de la respuesta transitoria
Sistemas de segundo orden (cont.):
Caso 3: Raíces complejas conjugadas.
Y (s) 
K3
K1
K2


s s    j d s    j d
Departamento de Ingeniería
de Sistemas y Automática
Reagrupando las dos fracciones complejas:
Y (s) 
K '3 d
K1
K '2 ( s   )


s ( s   ) 2  d 2 ( s   ) 2  d2
con K '2  2 Re( K 2 ) y K '3  2 Im(K 3 ).
Aplicando la transformada inversa de Laplace:
y (t )  ( K1  K '2 e t cos d t  K '3 e t sin d t )ue (t )
Se denominan sistemas subamortiguados.
Análisis de la respuesta transitoria
Sistemas de segundo orden (cont.):
Departamento de Ingeniería
de Sistemas y Automática
Caso 3: Raíces complejas conjugadas.
La forma de la respuesta depende de la colocación de los polos
( , d ).
Análisis de la respuesta transitoria
Sistemas de segundo orden (cont.):
– La respuesta de SSO admite otra representación alternativa
en función de los parámetros:
Departamento de Ingeniería
de Sistemas y Automática
• Ganancia, K
• Relación de amortiguamiento, x
• Frecuencia natural NO amortiguada, n
G (s) 
K n2
Expresión paramétrica
de los SSO
s 2  2x n s   n2
– Las raíces de la ecuación característica son:
s1 , s 2  x n  j n 1  x 2    j d
Análisis de la respuesta transitoria
Sistemas de segundo orden (cont.):
Departamento de Ingeniería
de Sistemas y Automática
– Con parámetros   x n constante de tiempo inversa y
d  n 1 - x 2 frecuencia natural amortiguada.

– Los SSO se pueden clasificar atendiendo al valor de la
constante de amortiguamiento y la ubicación de sus polos.
Análisis de la respuesta transitoria
Sistemas de segundo orden (cont.):
– x >1 (s1 y s2 reales distintos, parte real -)
SOBREAMORTIGUADO
– x =1 (s1 y s2 reales iguales, parte real -)
 CRÍTICAMENTE AMORTIGUADO
Departamento de Ingeniería
de Sistemas y Automática
– 0<x<1 (s1 y s2 conj. complejos, parte real -)
 SUBAMORTIGUADO
– x=0 (s1 y s2 sobre el eje imaginario)
 LÍMITE DE ESTABILIDAD
– -1<x<0 (s1 y s2 conj. complejos, parte real +)
 INESTABLE OSCILANTE
– x<-1 (s1 y s2 reales distintos, parte real +)
 INESTABLE NO OSCILANTE
Análisis de la respuesta transitoria
Sistemas de segundo orden (cont.):
Departamento de Ingeniería
de Sistemas y Automática
En base a la relación de amortiguamiento x constante y la
frecuencia natural no amortiguada, n , constante se establecen
los correspondientes lugares geométricos en el plano s.
Análisis de la respuesta transitoria
Sistemas de segundo orden (cont.):
Departamento de Ingeniería
de Sistemas y Automática
– Existe un conjunto de curvas normalizadas de respuesta
escalón de SSO para valores de (x , n ) .
– Para x  0 , el sistema responde con una oscilación mantenida.
Análisis de la respuesta transitoria
Sistemas de segundo orden (cont.):
Departamento de Ingeniería
de Sistemas y Automática
– Existe un conjunto de curvas normalizadas de respuesta
impulso de SSO para valores de (x , n ) .
Kn2
Y ( s)  2
s  2xn s  n2
Análisis de la respuesta transitoria
Sistemas de orden superior:
– Sistema de orden superior (SOS) queda descrito por la
función de transferencia
G( s) 
K ( s  z1 )(s  z2 ) ( s  zm )
( s  p1 )(s  p2 ) ( s  pn )
Departamento de Ingeniería
de Sistemas y Automática
con zi y pj ceros y polos en general complejos.
– La respuesta escalón de amplitud A será:
A K0 n Ki
Y (s)  G (s) 

s
s
i 1 s  pi
Análisis de la respuesta transitoria
Sistemas de orden superior (cont.):
Caso 1: Polos en general distintos.
Aplicando la transformada inversa
n
y (t )  ( K 0   K i e  pit )ue (t )
Departamento de Ingeniería
de Sistemas y Automática
i 1
• La contribución de K0 es relativa al régimen estacionario.
• La contribución de cada polo pi en la respuesta transitoria
depende la magnitud del residuo Ki y de su colocación relativa:
1. si Ki es bajo  su contribución es despreciable, y
2. si Re(pi)<0 con |Re(pi)| alto  su contribución es despreciable.
Análisis de la respuesta transitoria
Sistemas de orden superior (cont.):
Caso 1: Polos en general distintos.
Departamento de Ingeniería
de Sistemas y Automática
Forma de la respuesta no estandarizada.
Análisis de la respuesta transitoria
Sistemas de orden superior (cont.):
Caso 2: Polos en general múltiples.
La respuesta escalón de amplitud A será:
Departamento de Ingeniería
de Sistemas y Automática
Kj
A K0
K1
K2
Y (s)  G (s) 




rj
s
s s  p1 s  p2
(s  p j )
y aplicando transformada inversa:
y(t )  ( K 0  K1e
 p1t
 K 2e
 p2 t

Kj
rj  1!
t
r j 1  p j t
e
 )ue (t )
Para determinar la contribución de cada polo se sigue el mismo
razonamiento que el caso anterior.
Análisis de la respuesta transitoria
Sistemas de orden superior (cont.):
– Concepto de dominancia:
Departamento de Ingeniería
de Sistemas y Automática
• Los polos más cercanos al eje imaginario jω prevalecen, y se
denominan polos dominantes.
• Transformamos un SOS en un SPO (un único polo dominante) o
en un SSO (un par de polos dominantes).
– Criterio de dominancia:
Relación Re(pi) / Re(pd) > 5, suponiendo que no hay ceros en
cercanía de pd (efecto cancelación).
Análisis de la respuesta transitoria
Sistemas de orden superior (cont.):
Ejemplo:
Construimos la función de transferencia:
>> G = zpk([-4, -6+3.8j, -6-3.8j],[-1, -2, -6.2-4j, -6.2+4j, -10],[1])
(s+4) (s^2 + 12s + 50.44)
---------------------------------------(s+1) (s+2) (s+10) (s^2 + 12.4s + 54.44)
Pole-Zero Map
5
4
3
2
Imaginary Axis
Departamento de Ingeniería
de Sistemas y Automática
y la representamos los polos y ceros en el plano complejo:
1
0
-1
-2
-3
-4
-5
-10
-8
-6
-4
Real Axis
-2
0
Análisis de la respuesta transitoria
Sistemas de orden superior (cont.):
Ejemplo (cont.):
Aplicamos un escalón de amplitud unitaria a ambos sistemas y el
SOS y al SSO*:
Step Response
2
1.8
Original
Aproximado
1.4
Al eliminar los polos
no dominantes se
modifica la ganancia
estática del sistema.
1.2
Amplitude
Departamento de Ingeniería
de Sistemas y Automática
1.6
1
0.8
?
0.6
0.4
0.2
0
0
1
2
3
Time (sec)
4
5
6
Análisis de la respuesta transitoria
Sistemas de orden superior (cont.):
Ejemplo (cont.):
Aplicamos un escalón de amplitud unitaria a ambos sistemas y el
SOS y al SSO*:
Step Response
0.2
0.18
0.1 × G(s)*
Original
Aproximado
0.14
0.12
Amplitude
Departamento de Ingeniería
de Sistemas y Automática
0.16
0.1
0.08
0.06
0.04
0.02
0
0
1
2
3
Time (sec)
4
5
6
Análisis de la respuesta transitoria
Sistemas de orden superior (cont.):
Departamento de Ingeniería
de Sistemas y Automática
Efecto de añadir un polo/cero al sistema:
•
Los polos de G(s) afectan a los exponentes en los términos
pt
exponenciales e i (pi complejo en general) de la respuesta
transitoria.
•
Los ceros de G(s) no afectan a los exponentes en los términos
exponenciales, pero afectan las magnitudes y los signos de los
residuos.
•
Por ejemplo:
–
Si añadimos un polo real negativo  Influye con una nueva
exponencial  hace el sistema más lento y más estable
(relativamente).
–
Si añadimos un cero  Influye con un residuo  hace el sistema
más rápido y más inestable (relativamente).
Contenido
Tema 4.- Análisis de la respuesta temporal de
sistemas LTI
4.1. Introducción.
4.2. Análisis de la respuesta transitoria de un sistema:
Departamento de Ingeniería
de Sistemas y Automática
4.2.1. Sistemas de primer orden.
4.2.2. Sistemas de segundo orden.
4.2.3. Sistemas de orden superior.
4.3. Introducción a la identificación de sistemas.
4.4. Análisis de la estabilidad de un sistema.
4.5. Análisis de la respuesta estacionaria de un sistema.
Identificación de sistemas
Identificación:
– Proceso de determinación de un modelo a partir del
conocimiento previo sobre el sistema y experiencias
prácticas realizadas sobre él.
Departamento de Ingeniería
de Sistemas y Automática
– Durante el proceso de identificación el sistema es
considerado como “caja negra”, realizándose experimentos
que proporcionan pares E/S.
u1
um
SISTEMA
DINAMICO
y1
yn
– Gran variedad de métodos de identificación, particularmente
los métodos de identificación paramétricos, a través de la
obtención de coeficientes de G(s).
Identificación de sistemas
Método de análisis del transitorio:
– Consiste en la aplicación de entradas tipo (escalón, impulso,
senoidal,…) analizando la forma de la respuesta transitoria
para determinar los parámetros del modelo del sistema.
Departamento de Ingeniería
de Sistemas y Automática
u( t )
SISTEMA
y(t )
– Asume comportamiento lineal o linealizado en torno a un
punto.
– Se considerará la identificación de sistemas de primer orden,
segundo orden y orden superior ante señal de entrada
escalón.
Identificación de sistemas
Sistemas de primer orden:
Departamento de Ingeniería
de Sistemas y Automática
K
G( s) 
Ts  1
1. Determinación de K por relación de amplitud salida-entrada
en régimen estacionario.
y ()  lim y (t )
y ( )
K
t 
u (sobre
)
2. Determinación de T por inspección
0.63%
de sYy(t)
y ()  lim
(s)
s 0
estacionaria.
3. Efecto posible de retardo e
 Tr s
.
Identificación de sistemas
Circuito RC:
Departamento de Ingeniería
de Sistemas y Automática
¿¿ vi(∞) y vo(∞) ??
5.8 × 0.2V=
1.15V
K
G( s) 
Ts  1
63% de 1.15 =
0.72V
K
vo 1,15

1
vi 1,15
T  100E - 6
¿¿ T ??
0.4 × 250E-6 s.
= 100E-6 s.
Identificación de sistemas
Circuito RC (cont.):
Step Response
1.2
1
Amplitude
Departamento de Ingeniería
de Sistemas y Automática
1.4
0.8
0.6
0.4
0.2
0
0
1
2
3
Time (sec)
4
5
6
-4
x 10
Identificación de sistemas
Sistemas de segundo orden:
Departamento de Ingeniería
de Sistemas y Automática
K n2
G( s)  2
s  2x n s   n2
1. Determinación de K por relación de amplitud entrada-salida
en régimen estacionario.
2. Determinación de x y n a través de los parámetros
característicos de la respuesta transitoria de segundo orden.
3. Efecto posible de retardo e
 Tr s
.
Identificación de sistemas
Sistemas de segundo orden (cont.):
1. tiempo de subida, de 0% al 100% del valor final y(∞)
 
tr 
d
 d 

  
  tan 1 
Departamento de Ingeniería
de Sistemas y Automática
2. tiempo de pico, en el valor máximo de y(t)

tp 
d
3. sobreoscilación, definida por
SO 
y (t p )  y ()
y ( )
SO  e
100%



d
e
 t p
e
x  n t p
 e  tan
 100%
Identificación de sistemas
Sistemas de segundo orden (cont.):
4. tiempo de establecimiento, para alcanzar el régimen
permanente (y(t) incluida en banda sobre y(∞)).
Departamento de Ingeniería
de Sistemas y Automática
Dos criterios: t s 
4

(2 %) y t s 
3

(5 %)
Identificación de sistemas
Sistemas de segundo orden (cont.):
A partir de un par de valores de parámetros característicos se
determinan  y d y de ahí x y n según:
Departamento de Ingeniería
de Sistemas y Automática
  x n
x  cos   0.12
d  n 1 - x 2
– Para sistemas sobreamortiguados solo será posible computar la
función de transferencia en caso de que haya dominancia de primer
orden, ya que las expresiones anteriores son validas para el caso
subamortiguado.
– Para sistemas de orden superior no hay un método de validez
general, si bien hay métodos para el caso de sistemas con polos en
situación específica (dominancia de primero o segundo orden, …).
Identificación de sistemas
Circuito RLC:
¿¿ vi(∞) y vo(∞) ??
K n2
G( s)  2
s  2x n s   n2
1.9 × 0.5V=
0.95V
Departamento de Ingeniería
de Sistemas y Automática
tp 


 d   0.21 106
d
tp
SO 
v(t p )  v()
v ( )
 t p
¿¿ vo (tp) ??
SO  e
0.15×100E-6 s.=
15E-6 s.
3.2 × 0.5V=
1.6V
  tan 1 
v
0,95
K o 
1
vi 0,95
x  cos   0.12
¿¿ t p ??
 
 0.68
 ln( SO)
 2.57  10 4
tp
 d 
  1.45

 
  x n   n 

 2.11 10 5
x
Identificación de sistemas
Circuito RLC:
Step Response
1.6
1.2
1
Amplitude
Departamento de Ingeniería
de Sistemas y Automática
1.4
0.8
0.6
0.4
0.2
0
0
0.5
1
1.5
Time (sec)
2
2.5
-4
x 10
Contenido
Tema 4.- Análisis de la respuesta temporal de
sistemas LTI
4.1. Introducción.
4.2. Análisis de la respuesta transitoria de un sistema:
Departamento de Ingeniería
de Sistemas y Automática
4.2.1. Sistemas de primer orden.
4.2.2. Sistemas de segundo orden.
4.2.3. Sistemas de orden superior.
4.3. Introducción a la identificación de sistemas.
4.4. Análisis de la estabilidad de un sistema.
4.5. Análisis de la respuesta estacionaria de un sistema.
Estabilidad de sistemas continuos
Análisis de estabilidad:
Departamento de Ingeniería
de Sistemas y Automática
– La estabilidad es una característica del sistema que asegura
que ante cualquier entrada acotada el sistema responde con
unas salida acotada.
– La estabilidad de un sistema LTI (lineal e invariante en el
tiempo) queda asegurada si todas las raíces del polinomio
característico se encuentran en el semiplano izquierdo del
plano complejo s.
– En efecto, para un sistema definido por
P( s )
G( s) 
an s n  an 1s n 1    a1s  a0
las raíces (polos) son la solución de
Q( s)  an s n  an 1s n 1    a1s  a0  0  s1 , s2 ,, sn
Estabilidad de sistemas continuos
Análisis de estabilidad:
– Si se encuentra alguna raíz con Re(si )   i  0 el término
K i e it crecería con el
correspondiente de la respuesta
tiempo, resultando por tanto el sistema inestable.
j
s1
x
Departamento de Ingeniería
de Sistemas y Automática
s2
x
0

x
s1*
– La estabilidad no depende de la función de entrada, es una
característica del sistema. Las raíces de la entrada
contribuyen solamente en los términos de respuesta
estacionaria en la solución.
Estabilidad de sistemas continuos
Análisis de estabilidad:
Departamento de Ingeniería
de Sistemas y Automática
– El problema más importante de los sistemas de control lineal
tiene que ver con la estabilidad, sobre todo con la estabilidad
del sistema en bucle cerrado.
Estabilidad de sistemas continuos
Análisis de estabilidad:
– El problema más importante de los sistemas de control lineal
tiene que ver con la estabilidad, sobre todo con la estabilidad
del sistema en bucle cerrado.
G (s )
+
G (s )
K
Departamento de Ingeniería
de Sistemas y Automática
-
H (s )
Gbc ( s ) 
KG ( s )
1  KG ( s ) H ( s )
– De hecho, sistemas que son inestables en bucle abierto,
pueden ser estabilizados al cerrar el bucle de control (K
variable).
OJO: no confundir con la
ganancia estática del sistema
Estabilidad de sistemas continuos
Criterio de Routh-Hurwitz:
Departamento de Ingeniería
de Sistemas y Automática
– Es un método para determinar la estabilidad de un sistema
sin tener que factorizar el polinomio característico del
sistema.
– Este procedimiento no especifica la posición concreta de las
raíces, sino el número de raíces existentes en el semiplano
derecho (inestabilidad) y en el eje imaginario (estabilidad
critica).
Límite de la estabilidad
Estabilidad de sistemas continuos
Pasos del criterio de Routh-Hurwitz:
1. Escribir el polinomio característico (suponiendo a0  0 ).
Departamento de Ingeniería
de Sistemas y Automática
an s n  an 1s n 1    a1s  a0  0
2. Si cualquier ai ≤ 0 en presencia de, al menos, un aj > 0,
entonces hay una raíz o raíces que son imaginarias o con
una parte real positiva, siendo el sistema crítico-estable o
inestable.
Límite de la estabilidad
Estabilidad de sistemas continuos
Pasos del criterio de Routh-Hurwitz (cont.):
3. Si todos los coeficientes ai  0 , agrupar los coeficientes en
el siguiente arreglo:
n
n 1
a
s

a
s
   a1s  a0  0
n
n

1
n
an  2
...
a0
s n 1 an 1
s n  2 bn  2
an 3
...
a1

b0
s n  3 cn  3
 

s2
e2
e0
s1
f1
s0
g0
s
Departamento de Ingeniería
de Sistemas y Automática
bi
an

bn 2 
an 1an 2  an an 3
an 1
Estabilidad de sistemas continuos
Pasos del criterio de Routh-Hurwitz (cont.):
3. Si todos los coeficientes ai  0 , agrupar los coeficientes en
el siguiente arreglo:
n
n 1
a
s

a
s
   a1s  a0  0
n
n

1
n
Departamento de Ingeniería
de Sistemas y Automática
bi
an
an  2
...
a0
s n 1 an 1
s n  2 bn  2
an 3
...
a1

b0
s n  3 cn  3
 

s2
e2
e0
s1
f1
s0
g0
s

b0 
a n 1a0  a n a1
a n 1
De la misma forma se evalúan ci , di , ei , continuando el
proceso hasta completar la última fila.
Estabilidad de sistemas continuos
Pasos del criterio de Routh-Hurwitz (cont.):
Departamento de Ingeniería
de Sistemas y Automática
4. El criterio de estabilidad de Routh-Hurwitz establece que el
número de raíces del polinomio característico con parte real
positiva es igual al número de cambios de signo de los
coeficientes de la primera columna del arreglo.
Hay casos especiales de ceros en la primera columna o ceros en
una fila que producen sistemas inestables o crítico-estables.
Límite de la estabilidad
Estabilidad de sistemas continuos
Casos especiales del criterio de Routh-Hurwitz:
Departamento de Ingeniería
de Sistemas y Automática
a) Si el término de la primera columna de cualquier fila es cero,
pero los términos restantes no lo son, o no hay términos
restantes.
Sustituimos el cero con
un número positivo muy
pequeño e y se evalúa
normalmente.
0
Si el signo del coeficiente que está encima del cero (e) es
igual al signo que está por debajo de él, entonces hay un par
de raíces imaginarias.
Estabilidad de sistemas continuos
Casos especiales del criterio de Routh-Hurwitz:
Departamento de Ingeniería
de Sistemas y Automática
a) Si el término de la primera columna de cualquier fila es cero,
pero los términos restantes no lo son, o no hay términos
restantes.
dos raíces en el
semiplano derecho
Primer cambio
Segundo cambio
Si el signo del coeficiente que está encima del cero (e) es
opuesto al del que está abajo entonces hay un cambio de
signo.
Estabilidad de sistemas continuos
Casos especiales del criterio de Routh-Hurwitz (cont.):
Departamento de Ingeniería
de Sistemas y Automática
b) Si todos los coeficientes de cualquier fila son cero.
La evaluación del resto de la tabla
continúa mediante la formación
de un polinomio auxiliar, P(s), con
los coeficientes de la fila s4 y
sustituyendo la fila s3 con los
coeficientes de la derivada de P(s)
con respecto a s.
raíces de igual magnitud que se
encuentran radialmente opuestas en el
plano s, es decir, dos raíces con
magnitudes iguales y signos opuestos
y/o dos raíces imaginarias conjugadas.
hay dos pares de raíces de igual
magnitud y signo opuesto. Estos pares
se obtienen resolviendo la ecuación
del polinomio auxiliar P(s)=0.
Estabilidad de sistemas continuos
Casos especiales del criterio de Routh-Hurwitz (cont.):
Departamento de Ingeniería
de Sistemas y Automática
b) Si todos los coeficientes de cualquier fila son cero.
Cambio de signo
Despejando las raíces de la ecuación del polinomio auxiliar
se obtiene:
Estabilidad de sistemas continuos
Ejemplo 1:
Aplicar el criterio de estabilidad de Routh-Hurwitz al polinomio de
tercer orden:
Todos los coeficientes > 0
Departamento de Ingeniería
de Sistemas y Automática
Se construye el arreglo:
Entonces, la condición de estabilidad viene dada por:
Estabilidad de sistemas continuos
Ejemplo 2:
Aplicar el criterio de estabilidad de Routh-Hurwitz al polinomio de
tercer orden:
Todos los coeficientes > 0
Departamento de Ingeniería
de Sistemas y Automática
Se construye el arreglo:
Hay dos cambios de signo en los coeficientes de la primera
columna. Esto significa que existen dos raíces con partes reales
positivas.
Estabilidad de sistemas continuos
Criterio de Routh-Hurwitz:
Departamento de Ingeniería
de Sistemas y Automática
– El criterio de estabilidad de Routh-Hurwitz tiene una utilidad
limitada en el análisis de un sistema de control lineal, sobre
todo porque no sugiere cómo mejorar la estabilidad relativa ni
cómo estabilizar un sistema inestable.
– Sin embargo, es posible determinar los efectos de cambiar
uno o dos parámetros de un sistema si se examinan los
valores que producen inestabilidad.
Estabilidad de sistemas continuos
Ejemplo:
Departamento de Ingeniería
de Sistemas y Automática
Hallar K para que el sistema en bucle cerrado sea estable.
1. La función de transferencia en bucle cerrado es:
2. Siendo la ecuación característica la siguiente:
Estabilidad de sistemas continuos
Ejemplo (cont.):
Departamento de Ingeniería
de Sistemas y Automática
Hallar K para que el sistema en bucle cerrado sea estable.
3. Construimos el arreglo de Routh-Hurwitz:
4. Para garantizar la estabilidad del sistema todos los
elementos de la primera columna deben ser positivos.
Contenido
Tema 4.- Análisis de la respuesta temporal de
sistemas LTI
4.1. Introducción.
4.2. Análisis de la respuesta transitoria de un sistema:
Departamento de Ingeniería
de Sistemas y Automática
4.2.1. Sistemas de primer orden.
4.2.2. Sistemas de segundo orden.
4.2.3. Sistemas de orden superior.
4.3. Introducción a la identificación de sistemas.
4.4. Análisis de la estabilidad de un sistema.
4.5. Análisis de la respuesta estacionaria de un sistema.
Análisis de la respuesta estacionaria
Error en régimen permanente:
Error estacionario
Se considerarán sistemas en bucle cerrado donde la salida del
sistema C(s) tiene que seguir una consigna o referencia R(s).
R(s )
E (s)
+
G (s )
Señal referencia
-
C (s)
B(s)
Departamento de Ingeniería
de Sistemas y Automática
Señal medida
H (s )
– Se define el error en régimen permanente como:
e  lim e(t )  lim r (t )  b(t ) 
t 
t 
OJO: El erp es igual a la diferencia entre R(s)
y C(s) si y sólo si la ganancia H(s) es 1.
Análisis de la respuesta estacionaria
Error en régimen permanente (cont.):
– Aplicando la transformada de Laplace:
E(s)  R(s)  B(s)  R(s)  H (s)C(s)
Departamento de Ingeniería
de Sistemas y Automática


G (s)

E ( s )  R ( s )1  H ( s )
1  G (s) H (s) 

E (s) 
1
R( s)
1  G (s) H (s)
– El erp depende de la señal de referencia, r(t), y del “tipo” N de
la planta-sensor G(s)H(s), definido por:
G (s) H (s) 
K ( 1c s  1)...( mc s  1)
s N ( 1 p s  1)...( np s  1)
– Aplicando teorema del valor final:
s
R( s)
s 0 1  G ( s ) H ( s )
e  lim e(t )  lim sE ( s )  lim
t 
s 0
OJO: Sólo si e(t) está
acotado sistema estable.
Análisis de la respuesta estacionaria
Error en régimen permanente (cont.):
Error de posición
1. erp ante señal de entrada escalón:
s
1
s 0 1  G ( s ) H ( s ) s
e  lim
e 
1
1  lim G ( s ) H ( s )
s0
Departamento de Ingeniería
de Sistemas y Automática
Definiendo la constante de error escalón Kp
K p  lim G ( s ) H ( s )  G (0) H (0)
s 0
e 
1
1 K p
Cte. estática de error
de posición
El erp ante una entrada escalón será nulo cuando Kp∞ lo cual
se produce cuando el sistema es de tipo N ≥ 1.
Número de
integradores ≥1
Análisis de la respuesta estacionaria
Error en régimen permanente (cont.):
Error de velocidad
1. erp ante señal de entrada rampa:
s
1
s 0 1  G ( s ) H ( s ) s 2
1
1

s 0 s  sG ( s ) H ( s )
lim sG ( s ) H ( s )
e  lim
e  lim
s 0
Departamento de Ingeniería
de Sistemas y Automática
Definiendo la constante de error rampa Kv
K v  lim sG ( s ) H ( s )
s 0
e 
1
Kv
Cte. estática de
error de velocidad
El erp ante una entrada rampa será nulo cuando Kv∞ lo cual se
produce cuando el sistema es de tipo N ≥ 2.
Número de
integradores ≥2
Análisis de la respuesta estacionaria
Error en régimen permanente (cont.):
Error de aceleración
1. erp ante señal de entrada parábola:
s
1
s 0 1  G ( s ) H ( s ) s 3
e  lim
e  lim
s 0
1
s 2  sG ( s ) H ( s )

1
lim s 2 G ( s ) H ( s )
s 0
Departamento de Ingeniería
de Sistemas y Automática
Definiendo la constante de error parábola Ka
K a  lim s 2 G ( s ) H ( s )
s 0
e 
1
Ka
Cte. estática de error
de aceleración
El erp ante una entrada parábola será nulo cuando Ka∞ lo cual
se produce cuando el sistema es de tipo N ≥ 3.
Número de
integradores ≥3
Análisis de la respuesta estacionaria
Error en régimen permanente (cont.):
Departamento de Ingeniería
de Sistemas y Automática
Tabla resumen de errores en régimen permanente de sistemas
continuos: