Archivo de la etiqueta: Sistemas gradiente

Ecuaciones Diferenciales I – Videos: Sistemas gradiente

Por Eduardo Vera Rosales

Introducción

En la entrada anterior consideramos nuevamente el sistema de ecuaciones que modelaba el movimiento de un péndulo, en este caso con fricción, de la forma $$\begin{array}{rcl} \dot{x} & = & y \\ \dot{y} & = & -\sin{x}-y \end{array}$$ y analizamos las diferencias que existen entre este sistema y el que modela al péndulo simple (sin fricción), entre ellas que el nuevo sistema ya no es hamiltoniano. Sin embargo, la función hamiltoniana $H(x,y)=\frac{1}{2}y^{2}-\cos{x}+1$ que define al sistema del péndulo simple fue útil para esbozar el plano fase del sistema con fricción.

Lo anterior nos motivó a estudiar a un conjunto de funciones que compartieran las mismas propiedades de $H$, en función de los sistemas de ecuaciones y sus puntos de equilibrio. A dichas funciones $L$ las llamamos funciones de Lyapunov. Estudiamos sus propiedades y establecimos el teorema de estabilidad de Lyapunov que nos indica la estabilidad de puntos de equilibrio cuando el sistema admite una función de Lyapunov.

En esta entrada estudiaremos un tipo particulas de sistemas para los cuales es relativamente sencillo encontrar una función de Lyapunov, bajo ciertas hipótesis que deben cumplir los puntos de equilibrio del sistema. Tales sistemas son llamados sistemas gradiente y son de la forma $$\dot{\textbf{X}}=-\nabla{G}(\textbf{X})$$ para alguna función $G: \mathbb{R}^{n} \rightarrow \mathbb{R}$, que por comodidad, la tomaremos de clase $C^{\infty}$. En términos de sistemas de dos ecuaciones, que son las hemos estado estudiando, tenemos que un sistema es gradiente si $$ \begin{array}{rcl} \dot{x} & = & -\frac{\partial{G}}{\partial{x}} \\ \dot{y} & = & – \frac{\partial{G}}{\partial{y}} \end{array}$$ para alguna función $G: \mathbb{R}^{2} \rightarrow \mathbb{R}$. Estudiaremos sus principales propiedades y como ya mencionamos veremos cuándo este tipo de sistemas admite una función de Lyapunov.

Sistemas gradiente

En el primer video definimos a los sistemas gradientes y estudiamos sus principales propiedades. Además, encontramos una función de Lyapunov para puntos de equilibrio que son mínimos locales estrictos para la función $G$ que define al sistema gradiente.

En el segundo video estudiamos un par de sistemas gradientes, así como su plano fase.

Tarea moral

Los siguientes ejercicios no forman parte de la evaluación del curso, pero te servirán para entender mucho mejor los conceptos vistos en esta entrada, así como temas posteriores.

  • Sea $\dot{\textbf{X}}=-\nabla{G(\textbf{X})}$ un sistema gradiente. Prueba que $\dot{G}(\textbf{X})=0 \iff \textbf{X}$ es un punto de equilibrio del sistema.
  • Considera el sistema gradiente $\dot{\textbf{X}}=-\nabla{G(\textbf{X})}$. Demuestra que si $(x_{0},y_{0})$ es un mínimo local estricto para $G$, entonces es asintóticamente estable.
  • Demuestra que el sistema de dos ecuaciones, obtenido por la linealización de un sistema gradiente en un punto de equilibrio, tiene únicamente valores propios reales.
  • Considera el sistema de ecuaciones $$\begin{array}{rcl} \dot{x} & = & y^{2}+2xy \\ \dot{y} & = & x^{2}+2xy \end{array}$$ que estudiamos en el primer ejemplo del segundo video. Linealiza el sistema en el origen y encuentra su estabilidad.
  • Supongamos que un sistema hamiltoniano es gradiente. Prueba que entonces la función hamiltoniana $H$ es armónica, es decir, se satisface $$\frac{\partial^{2}{H}}{\partial{x}^{2}}+\frac{\partial^{2}{H}}{\partial{y}^{2}}=0.$$ ¿El recíproco es cierto?
  • Da un ejemplo de un sistema que sea hamiltoniano y gradiente a la vez.
  • Verifica que el sistema $$\begin{array}{rcl} \dot{x} & = & x^{2}-2xy \\ \dot{y} & = & y^{2}-x^{2} \end{array}$$ es gradiente. Esboza su plano fase.
  • Considera la función $G(x,y)=x^{4}+2y^{4}$. Define un sistema gradiente con dicha función, encuentra sus puntos de equilibrio, encuentra una función de Lyapunov para el sistema (si existe) y determina la estabilidad de los puntos de equilibrio.

Más adelante

En la siguiente entrada estudiaremos un poco de bifurcaciones en sistemas de dos ecuaciones de primer orden. Es decir, estudiaremos familias de sistemas de primer orden que varían respecto a un parámetro $\lambda \in \mathbb{R}$. Estudiaremos primero familias de sistemas lineales, y posteriormente analizaremos un par de sistemas no lineales, cuyas bifurcaciones son comunes en dicha área de las ecuaciones diferenciales.

¡Hasta la próxima!

Entradas relacionadas

Agradecimientos

Trabajo realizado con el apoyo del Programa UNAM-DGAPA-PAPIME PE104522 «Hacia una modalidad a distancia de la Licenciatura en Matemáticas de la FC-UNAM – Etapa 2»