Análisis de Estabilidad de Yde: Sistemas Sistemas de Control
Análisis de Estabilidad de Yde: Sistemas Sistemas de Control
sistemas y de sistemas
de control
Entrada Salida
Estabilidad: en base a la respuesta del Sistema
Intercambiador
de calor Salida
Alimentación Producto
Entrada
Vapor
Estabilidad: en base a la respuesta del Sistema
• Respuesta al impulso: será estable si la respuesta tiende a cero en un
tiempo infinito
Entrada
Salida
Estabilidad: en base a la respuesta del Sistema
Entrada
Salida
Estabilidad: en base a la respuesta del Sistema
Alimentación Producto
Entrada
Vapor
Estabilidad sobre el plano S por verificación de
los polos
Estabilidad: sobre el plano S
Sistema estable
1
𝑅 𝑠 = 𝐸(𝑠) C(𝑠) 3
𝑠 3 𝐷 𝑠 =1+ =0
𝑠 𝑠+1 𝑠+2 𝑠 𝑠+1 𝑠+2
−
𝐷(𝑠) = 𝑠 𝑠 + 1 𝑠 + 2 + 3 = 0
x 1.047j
x
-2.672 -0.164
x -1.047j
Polos de Lazo cerrado y respuesta
Sistema inestable
1
𝑅 𝑠 = 𝐸(𝑠) C(𝑠) 7
𝑠 7 𝐷 𝑠 =1+ =0
𝑠 𝑠+1 𝑠+2 𝑠 𝑠+1 𝑠+2
−
𝐷(𝑠) = 𝑠 𝑠 + 1 𝑠 + 2 + 7 = 0
1.505j x
x
-3.087 0.0434
-1.505j x
Consideraciones
𝑛𝑢𝑚(𝑠) ς𝑆𝑖=1 𝑠 + 𝛾𝑖
𝐺 𝑠 = = 2
∆(𝑠) ς𝑄 𝑠 + 𝑝𝑘 ς𝑅𝐼=1 𝑠 2 + 2𝛼𝑙 𝑠 + (𝛼𝑙 + 𝜔𝑙 2 )
𝑘=1
• Condición de estabilidad: todos los polos del sistema deben tener parte
real negativa. Es decir, todos los polos deben estar en el semiplano
izquierdo del plano s (la región estable)
𝐼𝑚
Región Región
Inestable Inestable
𝑅𝑒
Seleccionamos
una de las
entradas u(t)
• Análisis de estabilidad
Estabilidad Absoluta
Estabilidad Relativa
Comportamiento Estacionario
• Estabilidad Absoluta
Criterio de Routh
Criterio de Nyquist, etc.
• Estabilidad Relativa
Lugar de Raíces (Root Locus)
Diagramas de Bode .
Diagramas de Nyquist, etc
Condición suficiente para Inestabilidad
𝑠 + 𝑎 𝑠 − 𝑎 = 𝑠 2 − 𝑎2
𝑠 − 𝑗𝑏 𝑠 + 𝑗𝑏 = 𝑠 2 + 𝑏 2
Polinomio característico:
Ecuación característica:
∆ 𝑠 =0
Una condición necesaria, pero no suficiente
• Forma factorizada:
∆ 𝑠 = 𝑎𝑛 𝑠 − 𝑝1 𝑠 − 𝑝2 … 𝑠 − 𝑝𝑛 = 0
operando:
∆ 𝑠 = 𝑎𝑛 [𝑠 𝑛 − 𝑝1 + 𝑝2 … + 𝑝𝑛 𝑠 𝑛−1 + 𝑝1 𝑝2 + 𝑝2 𝑝3 + 𝑝1 𝑝3 … 𝑠 𝑛−2
− 𝑝1 𝑝2 𝑝3 + 𝑝1 𝑝2 𝑝4 + ⋯ 𝑠 𝑛−3 … + (−1)𝑛 𝑝1 𝑝2 …𝑝𝑛 ] = 0
∆1 𝑠 = 𝑠 3 − 4𝑠 2 + 𝑠 + 6 INESTABLE
∆2 𝑠 = 𝑠 4 + 2𝑠 3 + 3𝑠 2 + 4𝑠 + 5 ?
∆3 𝑠 = 𝑠 4 + 2𝑠 3 + 5𝑠 2 + 4𝑠 + 5 ?
• Ejemplo:
𝑅(𝑠) 𝑁 𝑠 𝐶(𝑠)
𝑎4 𝑠 4 + 𝑎3 𝑠 3 + 𝑎2 𝑠 2 + 𝑎1 𝑠 + 𝑎0
𝑠4 𝑎4 𝑎2 𝑎0
𝑠3 𝑎3 𝑎1 0 Las primeras dos filas del arreglo se
obtienen directamente de la ecuación
𝑠2 del polinomio característico.
𝑠1
𝑠0
Criterio de Routh-Hurwitz para Estabilidad
• Completando el arreglo
𝑠4 𝑎4 𝑎2 𝑎0
𝑠3 𝑎3 𝑎1 0
Cada fila subsiguiente es
𝑎4 𝑎2 𝑎4 𝑎0 𝑎 0
2 − 𝑎 𝑎 − 𝑎 − 4 completada a partir de las
𝑠 3 1
= 𝑏1 3 0
= 𝑏2
𝑎3 0
=0
𝑎3 𝑎3 𝑎3 dos filas previas. Términos
ausentes son tratados como
𝑎3 𝑎1 𝑎3 0 𝑎3 0
− 𝑏 𝑏 − − ceros.
𝑏1 0 𝑏1 0
𝑠1
1 2
= 𝑐1 =0 =0
𝑏1 𝑏1 𝑏1
𝑏1 𝑏2 𝑏1 0 𝑏1 0
− − −
𝑐1 0 𝑐1 0 𝑐1 0
𝑠0 𝑐1
= 𝑑1
𝑐1
=0
𝑐1
=0
Ejemplo 1
𝑠3 1 31 0
Cualquier fila puede ser
𝑠2 10 1 0
multiplicada por una constante
1030 103
positiva sin cambiar los valores
1 31 1 0 1 0 de las filas subsiguientes. Dos
𝑠1 −
1 103 = −72
−
1 0 =0
−
1 0 =0
1 1 1 cambios de signo en la 1era
columna indican dos polos en el
1 103 1 0 1 0
𝑠0
−
−72 0 = 103 −
−72 0 = 0
−
−72 0 = 0
semiplano derecho.
−72 −72 −72
Ejemplo 2
∆1 𝑠 = 𝑠 3 − 4𝑠 2 + 𝑠 + 6
𝑠3 1 1
𝑠2 −4 ¡ INESTABLE !
6
∆ 𝑠 = 𝑎2 𝑠 2 + 𝑎1 𝑠 + 𝑎0
∆ 𝑠 = 𝑠 6 + 4𝑠 5 +3𝑠 4 +2𝑠 3 +𝑠 2 + 4𝑠 + 4
𝑠6 1 3 1 4
¡ INESTABLE !
𝑠5 4 2 4 0
𝑠1 −76/15 0
𝑠0 4
Caso Especial 1
∆ 𝑠 = 𝑠 3 +2𝑠 2 + 𝑠 + 2
𝑠3 1 1
¡ MARGINALMENTE ESTABLE !
𝑠2 2 2
∆ 𝑠 = 𝑠 3 − 3s + 2
𝑠3 1 −3
𝑠2 0≈𝜀 2
¡ INESTABLE !
𝑠1 −3 − 2/𝜀 0
𝑠0 2
Caso Especial 1
𝐶(𝑠) 10
=
𝑅(𝑠) 𝑠 5 + 2𝑠 4 + 3𝑠 3 +6𝑠 2 + 5𝑠 + 3
𝑠5 1 3 5
𝑠4 2 6 3
𝑠3 0 𝜀 7/2 0
𝑠2 6𝜀 − 7
3 0
𝜀
𝑠1 42𝜀 − 49 − 6𝜀 2
0 0
12𝜀 − 14
𝑠0 3 0 0
Caso Especial 1
Ejemplo 2
𝐶(𝑠) 10
= 5
𝑅(𝑠) 𝑠 + 2𝑠 4 + 3𝑠 3 +6𝑠 2 + 5𝑠 + 3
Primera
columna Dos cambios de signo en la 1era
𝜀=+ 𝜀=−
columna indica dos polos en el
𝑠5 1 + + semiplano derecho.
𝑠4 2 + +
» roots([1 2 3 6 5 3])
𝑠3 0 𝜀 + − ans =
0.3429 + 1.5083i
𝑠2 6𝜀 − 7
− + 0.3429 - 1.5083i
𝜀
-1.6681
𝑠1 42𝜀 − 49 − 6𝜀 2 + + -0.5088 + 0.7020i
12𝜀 − 14
+ +
-0.5088 - 0.7020i
𝑠0 3
Caso Especial 2
• Otro caso especial ocurre cuando toda una fila es cero. Esto indica que
hay raíces imaginarias o polos reales simétricos en relación al eje
imaginario. En este caos debemos usar un polinomio auxiliar
determinado por los coeficiente de la fila anterior ala que se hace cero y
reemplazar la fila cero por los coeficientes de la derivada de este
polinomio auxiliar que también son raíces del polinomio original
𝑠5 1 11 28
𝑠4 5 23 12
𝑠3 6.4 25.6 0
𝑠2 3 12 𝐷1 𝑠 = 3𝑠 2 + 12𝑠 0
𝑠1 0 0 raíces: 𝐷1 𝑠 = 0 𝑠 = ±𝑗2
𝑑𝐷1 𝑠
𝑠1 6 0 = 6𝑠
𝑑𝑠
𝑠0 12
¡ MARGINALMENTE ESTABLE !
Criterio de Routh-Hurwitz: Ejemplo
𝑅(𝑠) 𝐾 C(𝑠)
𝑠 𝑠 + 7 𝑠 + 11
∆ 𝑠 = 𝑠 3 +18𝑠 2 + 77𝑠 + 𝐾
−
𝑠3 1 77 𝐾
𝐾 < 77 ∗ 18 𝐾 < 1386
77 − >0
18
𝑠2 18 𝐾
𝐾>0
𝑠1 77 − 𝐾/18
Rango de K para estabilidad 0 < K < 1386
𝑠0 𝐾
Probar en Simulink para valores de
K fuera y dentro del rango hallado
K 2
𝑅(𝑠) + 𝐾 C(𝑠)
𝑠2 C
s
2 s 2
s
𝑠2 R 1 K 2 1
−
s 2
s s 1
+
2 + ( K 2 s )(s 1)
s2 2
𝑠 s ( s 1) ( K 2 s )
s 2 ( K 2 s )(s 1)
3 2
1 s s 2s K
𝑠+1
Ecuación característica
𝐼𝑚
𝑠 3 +𝑠 2 + 2𝑠 + 𝐾 = 0
x
s3 1 2
s2 1 K
x 𝑅𝑒
s 1 (2 - K )
s0 K
x
K 2 s 2 K 2 s j 2
Rango de Estabilidad K 2
s3 = -1
Conclusión