Estabilidad de Lyapunov

De Wikipedia, la enciclopedia libre
Saltar a: navegación, búsqueda

En matemáticas, la noción de estabilidad de Lyapunov se da en el estudio de los sistemas dinámicos.

De manera esquemática, diremos que un punto de equilibrio X_o de la ecuación diferencial homogénea \dot X=f(X) es estable si todas las soluciones a la ecuación que parten en un entorno de X_o se mantienen cerca de X_o para todo tiempo posterior.

Esta definición de estabilidad lleva el nombre de Aleksandr Liapunov, quien publicó en 1892 su tesis de doctorado El problema general de la estabilidad del movimiento, donde define este concepto.

Definición[editar]

Sea f:M\to TM un campo de vectores en una variedad diferenciable M. Consideremos la ecuación diferencial

\dot{x} = f(x),

p\in M tal que f(p)=0 (es decir, un punto de equilibrio de la ecuación). Diremos que p es:

  1. estable en el sentido de Lyapunov si para todo \epsilon > 0, existe \delta > 0 tal que si x es solución a la ecuación con \|x(0)-p\| < \delta, entonces para t \geq 0 tenemos \|x(t)-p\| < \epsilon.
  2. asintóticamente estable si cumple con el punto anterior y además el \delta > 0 puede elegirse de manera que \lim_{t \rightarrow +\infty} \|x(t)-p\| = 0.

Ejemplos[editar]

(1) Sea la ecuación diferencial en \R, \ \dot x= -x. El 0 es un punto de equilibrio de la ecuación. Veamos que es asintóticamente estable.

Si x_0\in\R entonces la solución a la ecuación con condición x(0)=x_0 es x(t)=x_0e^{-t}. Es fácil ver que para todo x_0\in\R tendremos que esa solución es decreciente y tiende a 0 cuando t\to+\infty.

Por lo tanto, dado \epsilon>0, tomando \delta=\epsilon se cumple: si  |x(0)| < \delta entonces |x(t)| < \epsilon para t \geq 0 y \lim_{t \rightarrow +\infty} |x(t)| = 0.


(2) Para la ecuación  \dot x= x el 0 también es un punto de equilibrio. Veamos que no es estable.

Si x_0\in\R entonces la solución a la ecuación con condición x(0)=x_0 es x(t)=x_0e^{t}.

Tomando \epsilon=1 tenemos que ningún \delta>0 sirve para la definición de estabilidad: dado \delta>0 la solución \varphi(t)=\frac\delta 2 e^t verifica |\varphi(0)|<\delta pero existe T>0 tal que \varphi(T)>1.


(3) Sea la ecuación \dot X=AX, donde A=\left(\begin{array}{rr} 0 & 1\\ -1 & 0     \end{array}  \right). Veamos que el origen es un punto de equilibrio estable pero no asintóticamente estable.

Para ello mostremos que si \varphi=(x,y) es solución a la ecuación entonces x^2+y^2 es constante: (x^2+y^2)'= 2xx'+2yy'=2xy-2yx=0 . Por lo tanto, toda solución que parte a distancia r del origen se mantendrá a distancia r siempre. Esto implica que el origen es estable pero no asintóticamente.

El caso lineal[editar]

Para el caso de ecuaciones en \R^n del tipo \dot X=AX, donde A\in M_n(\R), se conoce una clasificación completa de los casos en que el origen es un punto de equilibrio estable o asintóticamente estable, estudiando sus valores propios.

Si A tiene todos sus valores propios con parte real negativa entonces el origen es un punto de equilibrio asintóticamente estable. Si la matriz tiene algún valor propio con parte real positiva entonces el origen no es estable.

Para el caso en que A tenga valores propios con parte real nula se sabe que el origen no es asintóticamente estable. Para ver si es estable debemos estudiar las multiplicidades geométricas de dichos valores propios. Cuando la matriz tiene valores propios con parte real menor o igual a cero tendremos que: el origen es estable si y solo para todo valor propio \lambda con parte real 0 se tiene que la multiplicidad algebraica de \lambda es igual a la geométrica.

Algunos resultados[editar]

El teorema de Hartman-Grobman[editar]

Sea f:\R^n\to\R^n una función diferenciable. El teorema de Hartman-Grobman indica que para estudiar la estabilidad de un punto de equilibrio de la ecuación \dot X=f(X) puede utilizarse su aproximación lineal en algunos casos. Más en concreto: sea p\in\R^n tal que f(p)=0 y su matriz jacobiana J=J_f(p) no tiene valores propios con parte real nula, entonces p es (asintóticamente) estable si y solo si el origen es (asintóticamente) estable para la ecuación \dot X= J X.

Funciones de Lyapunov[editar]

Sea f:\R^n\to\R^n una función de clase C^1. Consideremos la ecuación \dot X=f(X). Supongamos p\in\R^n verifica f(p)=0.

Sea U un entorno de p, V:U\subseteq\mathbb{R}^n \to \mathbb{R} derivable tal que V(p) = 0 \,, V(x) > 0 \quad \forall x \in U\setminus\{0\}. A una función así la llamaremos función de Lyapunov. Para  x solución a la ecuación diferencial, la derivada de V\circ x es \dot{V}(x) = \frac{\partial V}{\partial x} \frac{dx}{dt} = \nabla V \cdot \dot{x} = \nabla V \cdot f(x).

Existen dos resultados debidos a Lyapunov que conciernen este tipo de funciones:

  1. si \dot{V}(x) \leq 0 \quad \forall x \in U entonces p es estable;
  2. si \dot{V}(x) < 0 \quad \forall x \in U\setminus\{p\} entonces p es asintóticamente estable.

Véase también[editar]

Bibliografía[editar]

  • Sotomayor, Jorge (1979). Lições de Equações Diferenciais Ordinárias (en portugués). Río de Janeiro: IMPA (Projeto Euclides). 
  • Gil, Omar (2002). «Ecuaciones diferenciales ordinarias: teoría básica». Curso introductorio a las ecuaciones diferenciales. Montevideo: IMERL (Facultad de Ingeniería, Universidad de la República). pp. 245–272. 
  • Imaz, Carlos; Vorel, Zdenek (1968). «El problema de estabilidad». Ecuaciones diferenciales ordinarias (primera edición). México D.F.: Limusa-Wiley S.A. pp. 123–156. 
  • Lyapunov, Aleksandr (1992). The general problem of stability of motion (en inglés) (primera edición). Londres: Taylor & Francis.