Propiedad de un sistema dinámico donde las soluciones cercanas a un punto de equilibrio permanecen así
Se pueden discutir varios tipos de estabilidad para las soluciones de ecuaciones diferenciales o ecuaciones en diferencias que describen sistemas dinámicos . El tipo más importante es el que se refiere a la estabilidad de soluciones cercanas a un punto de equilibrio. Esto se puede discutir mediante la teoría de Aleksandr Lyapunov . En términos simples, si las soluciones que comienzan cerca de un punto de equilibrio permanecen cerca para siempre, entonces es estable de Lyapunov . Más fuertemente, si es estable de Lyapunov y todas las soluciones que comienzan cerca convergen a , entonces se dice que es asintóticamente estable (ver análisis asintótico ). La noción de estabilidad exponencial garantiza una tasa mínima de decaimiento, es decir, una estimación de qué tan rápido convergen las soluciones. La idea de estabilidad de Lyapunov se puede extender a variedades de dimensión infinita, donde se conoce como estabilidad estructural , que se refiere al comportamiento de soluciones diferentes pero "cercanas" a ecuaciones diferenciales. La estabilidad de entrada a estado (ISS) aplica las nociones de Lyapunov a sistemas con entradas.
Historia
La estabilidad de Lyapunov debe su nombre a Aleksandr Mikhailovich Lyapunov , un matemático ruso que defendió la tesis El problema general de la estabilidad del movimiento en la Universidad de Járkov en 1892. [1] AM Lyapunov fue un pionero en los esfuerzos exitosos para desarrollar un enfoque global para el análisis de la estabilidad de los sistemas dinámicos no lineales en comparación con el método local ampliamente extendido de linealizarlos sobre puntos de equilibrio. Su trabajo, publicado inicialmente en ruso y luego traducido al francés, recibió poca atención durante muchos años. La teoría matemática de la estabilidad del movimiento, fundada por AM Lyapunov, anticipó considerablemente el tiempo para su implementación en la ciencia y la tecnología. Además, Lyapunov no hizo él mismo una aplicación en este campo, su propio interés estaba en la estabilidad de las masas de fluidos en rotación con aplicación astronómica. No tuvo estudiantes de doctorado que siguieran la investigación en el campo de la estabilidad y su propio destino fue terriblemente trágico debido a su suicidio en 1918 [ cita requerida ] . Durante varias décadas, la teoría de la estabilidad se hundió en el olvido total. El matemático y mecánico ruso-soviético Nikolai Guryevich Chetaev, que trabajaba en el Instituto de Aviación de Kazán en la década de 1930, fue el primero en darse cuenta de la increíble magnitud del descubrimiento realizado por A. M. Lyapunov. La contribución de N. G. Chetaev a la teoría fue tan significativa que muchos matemáticos, físicos e ingenieros lo consideran el sucesor directo de Lyapunov y el siguiente descendiente científico en la creación y desarrollo de la teoría matemática de la estabilidad.
El interés en este campo se disparó repentinamente durante el período de la Guerra Fría , cuando se descubrió que el llamado "Segundo Método de Lyapunov" (ver más abajo) era aplicable a la estabilidad de los sistemas de guía aeroespacial que suelen contener fuertes no linealidades que no se pueden tratar con otros métodos. Apareció entonces y desde entonces una gran cantidad de publicaciones en la literatura de control y sistemas. [3] [4] [5] [6] [7]
Más recientemente, el concepto del exponente de Lyapunov (relacionado con el Primer Método de Lyapunov para analizar la estabilidad) ha recibido un amplio interés en relación con la teoría del caos . Los métodos de estabilidad de Lyapunov también se han aplicado para encontrar soluciones de equilibrio en problemas de asignación de tráfico. [8]
Definición para sistemas de tiempo continuo
Consideremos un sistema dinámico no lineal autónomo
- ,
donde denota el vector de estado del sistema , un conjunto abierto que contiene el origen, y es un campo vectorial continuo en . Supongamos que tiene un equilibrio en de modo que entonces
- Se dice que este equilibrio es estable de Lyapunov si para cada existe un tal que si entonces para cada tenemos .
- Se dice que el equilibrio del sistema anterior es asintóticamente estable si es estable de Lyapunov y existe tal que si entonces .
- Se dice que el equilibrio del sistema anterior es exponencialmente estable si es asintóticamente estable y existen tales que si entonces para todo .
Conceptualmente, los significados de los términos anteriores son los siguientes:
- La estabilidad de Lyapunov de un equilibrio significa que las soluciones que comienzan "suficientemente cerca" del equilibrio (dentro de una distancia de éste) permanecen "suficientemente cerca" para siempre (dentro de una distancia de éste). Nótese que esto debe ser cierto para cualquiera que uno desee elegir.
- La estabilidad asintótica significa que las soluciones que comienzan lo suficientemente cerca no sólo permanecen lo suficientemente cerca sino que eventualmente convergen al equilibrio.
- La estabilidad exponencial significa que las soluciones no sólo convergen, sino que de hecho convergen más rápido o al menos tan rápido como una tasa particular conocida .
La trayectoria es (localmente) atractiva si
- como
para todas las trayectorias que comienzan lo suficientemente cerca de , y globalmente atractivas si esta propiedad se cumple para todas las trayectorias.
Es decir, si x pertenece al interior de su variedad estable , es asintóticamente estable si es a la vez atractiva y estable. (Hay ejemplos que muestran que la atractividad no implica estabilidad asintótica. [9] [10] [11] Tales ejemplos son fáciles de crear usando conexiones homoclínicas ).
Si el jacobiano del sistema dinámico en equilibrio resulta ser una matriz de estabilidad (es decir, si la parte real de cada valor propio es estrictamente negativa), entonces el equilibrio es asintóticamente estable.
Sistema de desviaciones
En lugar de considerar la estabilidad sólo cerca de un punto de equilibrio (una solución constante ), se pueden formular definiciones similares de estabilidad cerca de una solución arbitraria . Sin embargo, se puede reducir el caso más general al de un equilibrio mediante un cambio de variables llamado "sistema de desviaciones". Defina , obedeciendo a la ecuación diferencial:
- .
Éste ya no es un sistema autónomo, sino que tiene un punto de equilibrio garantizado cuya estabilidad es equivalente a la estabilidad de la solución original .
Segundo método de Lyapunov para la estabilidad
Lyapunov, en su trabajo original de 1892, propuso dos métodos para demostrar la estabilidad . [1] El primer método desarrolló la solución en una serie que luego se demostró convergente dentro de límites. El segundo método, que ahora se conoce como el criterio de estabilidad de Lyapunov o el método directo, hace uso de una función de Lyapunov V(x) que tiene una analogía con la función potencial de la dinámica clásica. Se presenta de la siguiente manera para un sistema que tiene un punto de equilibrio en . Considere una función tal que
- Si y sólo si
- Si y sólo si
- para todos los valores de . Nota: para estabilidad asintótica, se requiere para .
Entonces V(x) se denomina función de Lyapunov y el sistema es estable en el sentido de Lyapunov. (Obsérvese que se requiere; de lo contrario, por ejemplo, "probaría" que es localmente estable). Se requiere una condición adicional llamada "propiedad" o "ilimitación radial" para concluir la estabilidad global. La estabilidad asintótica global (GAS) se deduce de manera similar.
Es más fácil visualizar este método de análisis pensando en un sistema físico (por ejemplo, un resorte vibrante y una masa) y considerando la energía de dicho sistema. Si el sistema pierde energía con el tiempo y esta nunca se recupera, eventualmente el sistema debe detenerse y alcanzar un estado de reposo final. Este estado final se llama atractor . Sin embargo, encontrar una función que proporcione la energía precisa de un sistema físico puede ser difícil y, para sistemas matemáticos abstractos, sistemas económicos o sistemas biológicos, el concepto de energía puede no ser aplicable.
La conclusión de Lyapunov fue que se puede demostrar la estabilidad sin necesidad de conocer la verdadera energía física, siempre que se pueda encontrar una función de Lyapunov que satisfaga las restricciones anteriores.
Definición para sistemas de tiempo discreto
La definición de sistemas de tiempo discreto es casi idéntica a la de sistemas de tiempo continuo. La definición que aparece a continuación la proporciona utilizando un lenguaje alternativo que se utiliza comúnmente en textos más matemáticos.
Sea ( X , d ) un espacio métrico y f : X → X una función continua . Se dice que un punto x en X es estable de Lyapunov si:
Decimos que x es asintóticamente estable si pertenece al interior de su conjunto estable , es decir si,
Estabilidad para modelos de espacio de estados lineales
Un modelo de espacio de estados lineal
- ,
donde es una matriz finita, es asintóticamente estable (de hecho, exponencialmente estable ) si todas las partes reales de los valores propios de son negativos. Esta condición es equivalente a la siguiente: [12]
es definida negativa para alguna matriz definida positiva . (La función de Lyapunov relevante es ).
En consecuencia, un modelo de espacio de estados lineal discreto en el tiempo
es asintóticamente estable (de hecho, exponencialmente estable) si todos los valores propios de tienen un módulo menor que uno.
Esta última condición se ha generalizado a los sistemas conmutados: un sistema de tiempo discreto conmutado lineal (regido por un conjunto de matrices ).
es asintóticamente estable (de hecho, exponencialmente estable) si el radio espectral conjunto del conjunto es menor que uno.
Estabilidad para sistemas con entradas
Un sistema con entradas (o controles) tiene la forma
donde la entrada (generalmente dependiente del tiempo) u(t) puede verse como un control , una entrada externa , un estímulo , una perturbación o una función de fuerza . Se ha demostrado [13] que cerca de un punto de equilibrio que es estable según el método de Lyapunov, el sistema permanece estable bajo pequeñas perturbaciones. Para perturbaciones de entrada mayores, el estudio de tales sistemas es el tema de la teoría de control y se aplica en la ingeniería de control . Para sistemas con entradas, se debe cuantificar el efecto de las entradas en la estabilidad del sistema. Los dos enfoques principales para este análisis son la estabilidad BIBO (para sistemas lineales ) y la estabilidad de entrada a estado (ISS) (para sistemas no lineales ).
Ejemplo
Este ejemplo muestra un sistema en el que se puede utilizar una función de Lyapunov para demostrar la estabilidad de Lyapunov, pero no se puede demostrar la estabilidad asintótica. Considere la siguiente ecuación, basada en la ecuación del oscilador de Van der Pol con el término de fricción modificado:
Dejar
de manera que el sistema correspondiente sea
El origen es el único punto de equilibrio. Elijamos como función de Lyapunov
que es claramente definida positiva . Su derivada es
Parece que si el parámetro es positivo, la estabilidad es asintótica para Pero esto es incorrecto, ya que no depende de , y será 0 en todas partes del eje. El equilibrio es estable según el método de Lyapunov pero no asintóticamente estable.
El lema de Barbalat y la estabilidad de los sistemas variables en el tiempo
Puede resultar difícil encontrar una función de Lyapunov con una derivada definida negativa, como lo exige el criterio de estabilidad de Lyapunov; sin embargo, puede estar disponible una función que sea solo semidefinida negativa. En sistemas autónomos, el teorema de conjuntos invariantes se puede aplicar para demostrar la estabilidad asintótica, pero este teorema no es aplicable cuando la dinámica es una función del tiempo. [14]
En cambio, el lema de Barbalat permite un análisis similar al de Lyapunov de estos sistemas no autónomos. El lema está motivado por las siguientes observaciones. Suponiendo que f es una función del tiempo únicamente:
- Tener no implica que tenga un límite en . Por ejemplo, .
- El hecho de haberse acercado a un límite como no implica que . Por ejemplo, .
- Al tener un límite inferior y ser decreciente ( ) implica que converge a un límite. Pero no dice si converge o no a medida que .
El lema de Barbalat dice:
- Si tiene un límite finito como y si es uniformemente continua (una condición suficiente para la continuidad uniforme es que esté acotada), entonces como . [15]
Una versión alternativa es la siguiente:
- Sea y . Si y , entonces como [16]
En la siguiente forma el Lema es verdadero también en el caso de valor vectorial:
- Sea una función uniformemente continua con valores en un espacio de Banach y supongamos que tiene un límite finito cuando . Entonces , cuando . [17]
El siguiente ejemplo está tomado de la página 125 del libro de Slotine y Li Applied Nonlinear Control . [14]
Consideremos un sistema no autónomo
Esto no es autónomo porque la entrada es una función del tiempo. Supongamos que la entrada está acotada.
Tomando da
Esto dice que, por las dos primeras condiciones, y por lo tanto y están acotados. Pero no dice nada sobre la convergencia de a cero, ya que es solo semidefinida negativa (nota: puede ser distinta de cero cuando = 0) y la dinámica no es autónoma.
Utilizando el lema de Barbalat:
- .
Esto está acotado porque , y están acotados. Esto implica que como y por lo tanto . Esto demuestra que el error converge.
Véase también
Referencias
- ^ ab Lyapunov, AM El problema general de la estabilidad del movimiento (en ruso), Tesis doctoral, Univ. Kharkov 1892 Traducciones al inglés: (1) Stability of Motion , Academic Press, Nueva York y Londres, 1966 (2) El problema general de la estabilidad del movimiento , (AT Fuller trans.) Taylor & Francis, Londres 1992. Se incluye una biografía de Smirnov y una extensa bibliografía del trabajo de Lyapunov.
- ^ Chetaev, NG Sobre trayectorias estables de la dinámica, Kazan Univ Sci Notes, vol.4 no.1 1936; La estabilidad del movimiento, publicado originalmente en ruso en 1946 por ОГИЗ. Гос. изд-во технико-теорет. лит., Москва-Ленинград. Traducido por Morton Nadler, Oxford, 1961, 200 páginas.
- ^ Letov, AM (1955). Устойчивость нелинейных регулируемых систем [ Estabilidad de sistemas de control no lineales ] (en ruso). Moscú: Gostekhizdat.Traducción al inglés de Princeton, 1961
- ^ Kalman, RE ; Bertram, J. F (1960). "Análisis y diseño de sistemas de control mediante el "segundo método" de Lyapunov: I—Sistemas de tiempo continuo". Journal of Basic Engineering . 82 (2): 371–393. doi :10.1115/1.3662604.
- ^ LaSalle, JP ; Lefschetz, S. (1961). Estabilidad por el segundo método de Lyapunov con aplicaciones . Nueva York: Academic Press.
- ^ Parks, PC (1962). "El método de Liapunov en la teoría del control automático". Control . I Nov 1962 II Dic 1962.
- ^ Kalman, RE (1963). "Funciones de Lyapunov para el problema de Lur'e en control automático". Proc Natl Acad Sci USA . 49 (2): 201–205. Bibcode :1963PNAS...49..201K. doi : 10.1073/pnas.49.2.201 . PMC 299777 . PMID 16591048.
- ^ Smith, MJ; Wisten, MB (1995). "Un modelo continuo de asignación de tráfico diario y la existencia de un equilibrio dinámico continuo de usuarios". Anales de investigación de operaciones . 60 (1): 59–79. doi :10.1007/BF02031940. S2CID 14034490.
- ^ Hahn, Wolfgang (1967). Estabilidad del movimiento. Springer. pp. 191–194, Sección 40. doi :10.1007/978-3-642-50085-5. ISBN 978-3-642-50087-9.
- ^ Braun, Philipp; Grune, Lars; Kellett, Christopher M. (2021). (In)estabilidad de inclusiones diferenciales: nociones, equivalencias y caracterizaciones similares a las de Lyapunov. Springer. págs. 19-20, Ejemplo 2.18. doi :10.1007/978-3-030-76317-6. ISBN 978-3-030-76316-9.S2CID237964551 .
- ^ Vinograd, RE (1957). "La insuficiencia del método de exponentes característicos para el estudio de ecuaciones diferenciales no lineales". Doklady Akademii Nauk (en ruso). 114 (2): 239–240.
- ^ Goh, BS (1977). "Estabilidad global en sistemas de muchas especies". The American Naturalist . 111 (977): 135–143. doi :10.1086/283144. S2CID 84826590.
- ^ Malkin IG Theory of Stability of Motion, Moscú 1952 (Gostekhizdat) Cap. II, párrafo 4 (ruso) Trad. inglesa, Language Service Bureau, Washington AEC -tr-3352; originalmente On stability under constant acting disturbies Prikl Mat 1944, vol. 8, núm. 3 241-245 (ruso); Trad. americana de matemáticas y ciencias sociales, núm. 8
- ^ ab Slotine, Jean-Jacques E.; Weiping Li (1991). Control no lineal aplicado . Nueva Jersey: Prentice Hall.
- ^ I. Barbălat, Systèmes d'équations différentielles d'oscillations non Linéaires, Rev. Math. Pures Appl. 4 (1959) 267–270, pág. 269.
- ^ B. Farkas et al., Variaciones sobre el lema de Barbălat, Amer. Math. Monthly (2016) 128, n.º 8, 825-830, DOI: 10.4169/amer.math.monthly.123.8.825, pág. 827.
- ^ B. Farkas et al., Variaciones sobre el lema de Barbălat, Amer. Math. Monthly (2016) 128, n.º 8, 825-830, DOI: 10.4169/amer.math.monthly.123.8.825, pág. 826.
Lectura adicional
- Bhatia, Nam Parshad; Szegő, Giorgio P. (2002). Teoría de la estabilidad de sistemas dinámicos . Springer. ISBN 978-3-540-42748-3.
- Chervin, Robert (1971). Estabilidad de Lyapunov y control por retroalimentación de sistemas de plasma de dos corrientes (PhD). Universidad de Columbia.
- Gandolfo, Giancarlo (1996). Dinámica económica (tercera edición). Berlín: Springer. pp. 407–428. ISBN 978-3-540-60988-9.
- Parks, PC (1992). "La teoría de estabilidad de A. M. Lyapunov: 100 años después". IMA Journal of Mathematical Control & Information . 9 (4): 275–303. doi :10.1093/imamci/9.4.275.
- Slotine, Jean-Jacques E.; Weiping Li (1991). Control no lineal aplicado . Nueva Jersey: Prentice Hall.
- Teschl, G. (2012). Ecuaciones diferenciales ordinarias y sistemas dinámicos. Providence : American Mathematical Society . ISBN 978-0-8218-8328-0.
- Wiggins, S. (2003). Introducción a los sistemas dinámicos no lineales aplicados y al caos (2.ª ed.). Nueva York: Springer Verlag . ISBN 978-0-387-00177-7.
Este artículo incorpora material de asintóticamente estable en PlanetMath , que se encuentra bajo la licencia Creative Commons Attribution/Share-Alike License .