Estabilidad de Jury, Liapunov
Estabilidad de Jury, Liapunov
Estabilidad de Jury, Liapunov
Junio 2019.
In this document, it present the different methods that are used to evaluate the stability of
the systems in the discrete time, there is a brief difference between the systems in the continuous
time and the discrete time and the techniques necessary to apply in the discrete time, the methods
used in continuous time (Routh-Hurwitz and ). This technique consists of the bilinear
transformation where the circumference of the unit radius becomes the imaginary axis in the new
plane. It also presents the methods that are used directly on a discrete time system as the case of
the Jury´s stability criteria and the Lyapunov´s stability criteria and finally explains the basic
operation of some types of sensors that can be used in a system of control and the effect they
have on the stability of a system.
Tabla de Contenidos
Capítulo 1 Introducción...................................................................................................................1
Capítulo 2 Contenido......................................................................................................................2
Estabilidad de sistemas en lazo cerrado y lazo abierto en tiempo discreto.................................2
Criterio de Routh-Hurwitz.......................................................................................................4
Método del lugar de las raíces.................................................................................................6
Criterio de estabilidad de Jury.....................................................................................................8
Criterio de Estabilidad de Liapunov..........................................................................................11
Análisis de estabilidad de Liapunov de sistemas en tiempo discreto....................................12
Análisis de estabilidad de Liapunov de los sistemas lineales en tiempo discreto e invariantes
en el tiempo............................................................................................................................12
Sensores.....................................................................................................................................15
Capítulo 3 Conclusiones...............................................................................................................19
Referencias....................................................................................................................................20
Lista de figuras
Introducción
es la capacidad que tiene un sistema para amortiguar con el tiempo y anular las
ubicación de las raíces de la ecuación característica del sistema, es decir, para que un
sistema sea estable, los polos de la función de transferencia deben estar en el semiplano
izquierdo del plano s, o si es un sistema discreto dentro del círculo del plano z. en general
un sistema es estable si excitado con una entrada acotada, la salida permanece acotada. A
Contenido
Tanto en lazo abierto como en lazo cerrado, para que un sistema de tiempo
continuo sea estable se debe cumplir que los polos del sistema deben tener su parte real
característica:
Y (s )
G ( s )=
U (s )
Y ( s) G ( s)
=
R ( s ) 1+G ( s ) H ( s )
transferencia es:
G( z)
T ( z )=
1+GH ( z )
Para analizar la estabilidad de un sistema discreto se tiene que saber dónde deben
estar as raíces de la ecuación característica (1+GH(z)) para que el sistema tenga una
respuesta transitoria que se extinga en el tiempo. Para ello se calcula la respuesta del
ki z k z
Y ( z )= +…+ n +Y R ( z)
z −P i z−Pn
Donde Y R ( z) contiene los términos de Y(z) originados en los polos de la excitación R(z).
transitoria de Y(z). Para que el sistema sea estable, la anti transformada z de estos
ki z
de esta ecuación es Z
−1
{ }
z−Pi
=k i ( Pi )
kT
Para que este término se anule al incrementarse k, el módulo de Pi debe ser menor que la
unidad. Teniendo en cuenta que los valores Picorresponden a las raíces de la ecuación
discreto sea estable, las raíces de su ecuación característica deben estar ubicadas en el
interior del círculo de radio unitario. Las técnicas de análisis de la estabilidad de sistemas
eje imaginario jw. Estas técnicas no pueden ser aplicadas a sistemas de tiempo discreto,
2 z−1
w=
T z +1
plano w.
como Routh-Hurwitz y el método del lugar de las raíces.[ CITATION cat15 \l 6154 ]
Criterio de Routh-Hurwitz
El criterio de estabilidad de Routh dice si existen o no raíces inestables en una
ecuación polinomial, si tener que obtenerlas en realidad. Este criterio de estabilidad sólo
se aplica a los polinomios con una cantidad finita de términos. Cuando se aplica el
1. Se escribe el polinomio en w
menos un coeficiente positivo, hay una raíz o raíces imaginarias o que tienen
3. Todos los coeficientes deben ser positivos, esta condición es necesaria pero no
continuación:
w4 a0 a2 a4
w3 a1 a3
w2 b1 b2 a1 a 2−a0 a3
b 1=
w1 c1 a1
w0 d1
a1 a4−a0 a5
b 2= a 5=0 y no se calcula
a1
b1 a3−a1 b2
c 1=
b1
c 1 b2−b 1 c 2
d 1= c 2=0 y no se calcula
c1
El criterio de Routh establece que la cantidad de raíces positivas o con parte real
Puede verse que GH(z) es de la misma forma que GH(s) para sistemas continuos.
Por lo tanto, las reglas de construcción del diagrama del lugar de las raíces, desarrolladas
Para que un punto del plano z pertenezca al lugar de las raíces, GH(z) debe
Sea:
1
G ( s )=
s (s+1)
aT ( z−z 1 )
G ( z ) =K
( z−1)(z−P 1)
obtiene la ecuación de una circunferencia con centro en (z1, 0) y cuyo radio es la media
Además, fácilmente, puede verse que los puntos sobre el eje real entre el polo p1 y 1
pertenecen al lugar de raíces, lo mismo que los puntos sobre el eje real a la izquierda del
cero z1. En la figura 4 se puede observar el lugar de las raíces para un periodo de
muestreo de T=1seg
Figura 4. Lugar geométrico de las raíces con T= 1 seg
Para K = 0 las raíces coinciden con las raíces del sistema a lazo abierto. Cuando K
circunferencia para luego volver a ser reales. Por último, para K→∞ una raíz coincide
con el cero z1 y la otra está en z→∞. El sistema es estable hasta un valor de K para el
cual se produce la intersección del lugar de raíces con la circunferencia de radio unitario.
Se trata ahora de obtener una expresión del valor máximo de K para el cual el sistema
x=[ ( 1−P1 ) + z 1 ( 1+ P1 ) ] /2 z 1
Dónde
a 0 a n−k
b k=
|a n ak |
b 0 bn−1−k
ck=
|bn−1 bk |
c0 c n−2−k
dk =
|c n−2 ck |
La condición necesaria y suficiente para que el polinomio F(z) no tenga raíces
F ( 1 ) >0
(−1)n . F ( 1 ) >0
|a 0|<an
|b 0|>|b n−1|
|c 0|>|c n−2|
|d 0|>|d n−3|
Se puede notar que un sistema de segundo orden contiene una sola fila en el
arreglo. Por cada orden adicional, dos filas adiciones son agregadas. Además, para un
1
G ( s )=
s (s+1)
0.00484 z +0.00468
G ( z ) =K
( z−1)(z−0.905)
Dado que se trata de una ecuación de segundo orden, el arreglo de Jury se reduce a:
0 <K< 2.39.
z 2−0.49 z+ 1=0
Z=0.244 ∓ j 0.97
Este resultado coincide con el resultado obtenido del método del lugar geométrico
segundo método. Ambas de aplican para determinar la estabilidad de los sistemas de los
primer método está formado por procedimientos en las ecuaciones diferenciales o de las
ecuaciones en diferencias para el análisis. Por otra parte, el segundo método no requiere
de las soluciones de las ecuaciones diferenciales o en diferencias, por lo que resulta más
útil en la práctica. La ventaja de usar este método es que en diferencia los métodos
estabilidad de sistemas no lineales en los que las soluciones exactas no son posibles.
sistemas puramente matemáticos, no existe una forma simple de definir una función de
energía. Para vencer esta dificultada, Liapunov introdujo la función de Liapunov, una
función ficticia de energía. Esta idea es más general y utilizada que la de la energía. De
hecho, cualquier función escalar que satisfaga las hipótesis de los teoremas de estabilidad
∆ V ( x ( kT ) )=V ¿
x ( ( k +1 ) T ) =f ( x ( kT ))
muestreo.
x (k +1)=Gx (k )
∆ V ( x ( k ) )=V ¿
Dado que V ( x ( k ) ) se seleccionó para ser definida positiva, se requiere, para la estabilidad
∆ V ( x ( k ) )=−x¿ ( k ) Qx (k )
Dónde
De ese modo, para la estabilidad asintótica del sistema en tiempo discreto de la ecuación
especificar primero una matriz Q hermítica definida positiva (o simétrica real definida
positiva o no. Observe que una P definida positiva es una condición necesaria y
x ( k +1 )=Gx ( k )
Donde x es el vector estado (vector-n) y G es una matriz no singular constante de n*n.
Una condición necesaria y suficiente para que el estado de equilibrio x=0 sea
asintóticamente estable global es que, dada cualquier matriz Q hermítica definida positiva
(o simétrica real definida positiva) existe una matriz P hermítica definida positiva (o
x1 ( k +1) 1 x 1 (k )
= 0
[ ][
x2 ( k +1) −0.5 −1 ][ ]
x 2 (k )
ecuaciones:
p11 −2 p12=−1
De las cuales se obtiene
11 8 24
p11 = , p 12= , p 22=
5 5 5
En consecuencia
11 8
P= 5
| |
8
5
5
24
5
Observe que en vez de escoger Q como I se podría haber escogido Q como una
Q= 0 0
[ ]
0 1
¿
Siempre y cuando ∆ V ( x ( k ) )=−x ( k ) Qx (k )
∆ V ( x )=−x 22 ( k )
Para el sistema actual, x 2 ( k )idénticamente cero implica que x 1 ( k )es idénticamente cero.
origen. Por tanto, se puede escoger esta matriz Q semidefinida positiva para determinar la
3 4
P= 5
4
5
| | 5
12
5
Sensores
En un sistema de control está formado por varios elementos que realizan una
función individual, que en conjunto producen un fin para mantener estable un sistema o
de control recibe información del entorno por medio de sensores, esa información aporta
datos para que el controlador decida, si hay que realizar una acción, esta se lleva a cabo
Los sensores son dispositivos que recogen información del mundo real y la
procesar y tomar decisiones. Los sensores pueden ser analógicos, binarios y digitales
Los sensores analógicos: entregan una salida de nivel variable en función del
Los sensores binarios: entregan un todo o nada, es decir un cero o uno (1/0).
Los sensores digitales: dan información relativa a la medida con un protocolo de
pH, etc. Una magnitud eléctrica puede ser una resistencia eléctrica (un RTD), una
capacidad eléctrica (como un sensor de humedad), una tensión eléctrica (termopar), una
el sensor.
entrada es nula.
la magnitud de entrada.
ejemplos:
Posición lineal y angular: son dispositivos electrónicos utilizados para medir
posición de referencia, las tecnologías que más se utilizan son efecto Hall,
magnético por medio de corriente que circula por una bobina o un conductor,
Fuerza y par: son dispositivos que permiten obtener una señal eléctrica
por el flujo circulante. El pistón puede regresar a su posición inicial por gravedad
Cuando el pistón se mueve el imán se acerca y activa un reed switch, que cierra o
pared del sensor determina su sensibilidad, y por ende a qué caudal se activará el
sensor.
Presión: Es un dispositivo capaz de medir la presión de gases o líquido. La
tecnología como los sensores de presión capacitivos utilizan una membrana para
Wik19 \l 6154 ]
Capítulo 3
Conclusiones
sistemas de lazo abierto y lazo cerrado. Se comprobó que los métodos utilizados en los
sistemas tiempo continuo son aplicables a los sistemas de tiempo discreto haciendo
utilizar los métodos de Routh- Hurwitz y el método del lugar geométrico de las raíces.
discretos como es el caso de Jury y Liapunov, estos métodos son igual se efectivos para
obtener la estabilidad de los sistemas invariantes en el tiempo, ya que los otros métodos
comprendió estos son importantes ya que comparan el error entre la entrada y la salida y
mandan una señal para que este sea corregido, haciendo que el sistema sea realimentado