Un campo vectorial asigna un vector a cada punto del espacio. En el contexto de ecuaciones diferenciales, este vector indica la dirección y velocidad del cambio. Las partículas que "flotan" en el campo nos muestran cómo evoluciona cualquier condición inicial, revelando atractores, ciclos límite y puntos de equilibrio.
Lo que observarás en esta simulación
- Flechas del campo: vectores que indican la dirección local del flujo
- Partículas: "trazadores" que siguen las trayectorias del sistema
- Estelas: historial de posiciones que revela el retrato de fases
- Puntos fijos: donde las partículas convergen, divergen o giran
- Ciclos límite: órbitas cerradas estables en sistemas no lineales
Conceptos Fundamentales
Sistema Autónomo
Un sistema donde dx/dt y dy/dt dependen solo de (x,y), no del tiempo. Su comportamiento se representa completamente en el plano de fases.
Punto de Equilibrio
Un punto donde dx/dt = 0 y dy/dt = 0. Las partículas permanecen inmóviles allí. Su estabilidad determina el comportamiento cercano.
Retrato de Fases
El conjunto de todas las trayectorias posibles del sistema. Las estelas de las partículas revelan este retrato gradualmente.
Nulclina
Curva donde dx/dt = 0 (nulclina de x) o dy/dt = 0 (nulclina de y). Su intersección marca puntos de equilibrio.
Sistema de Ecuaciones
El campo vectorial es F(x,y) = (f(x,y), g(x,y)). En cada punto del plano, el vector F indica hacia dónde se mueve una partícula ubicada allí. La integración numérica (método RK4 en esta simulación) calcula las trayectorias.
Clasificación de Puntos Fijos
Para un punto fijo en el origen de un sistema lineal, la matriz jacobiana determina el comportamiento:
Los valores propios λ₁, λ₂ de J clasifican el punto:
(λ imaginarios puros)
(λ₁, λ₂ < 0)
(λ₁, λ₂ > 0)
(λ₁ > 0 > λ₂)
Experimento 1: Campos Lineales Clásicos
Objetivo: Identificar tipos de puntos de equilibrio
- Selecciona el campo Rotación (vórtice). Las ecuaciones son dx/dt = -y, dy/dt = x. Observa cómo las partículas orbitan el origen sin acercarse ni alejarse.
- Cambia a Sumidero (dx/dt = -x, dy/dt = -y). Todas las partículas convergen hacia el origen exponencialmente. Este es un nodo estable.
- Prueba Fuente. Es el comportamiento opuesto: todas las trayectorias escapan del origen. Nodo inestable.
- Selecciona Punto silla (dx/dt = x, dy/dt = -y). Observa cómo las partículas se acercan por un eje y se alejan por otro. Este es un punto de equilibrio inestable.
- Activa "Nulclinas" para visualizar dónde dx/dt = 0 (rosa). Nota que las nulclinas se cruzan en el origen.
Conclusión: El signo de los valores propios determina si el flujo atrae o repele. Los valores imaginarios producen rotación.
Experimento 2: Espirales y Estabilidad
Objetivo: Explorar puntos focales
- Selecciona Espiral entrante. Las ecuaciones dx/dt = -x - y, dy/dt = x - y combinan rotación con contracción.
- Aumenta la longitud de "Estela" a 150 para ver el patrón espiral claramente. Las partículas giran mientras convergen.
- Cambia a Espiral saliente. Ahora las partículas giran mientras se alejan del origen.
- Compara con el vórtice puro: la espiral tiene una componente radial además de la rotacional.
- Haz clic cerca del origen para añadir partículas. Observa cómo su distancia al centro cambia con el tiempo.
Conclusión: Los valores propios complejos con parte real negativa producen espirales estables; con parte real positiva, inestables.
Experimento 3: Van der Pol y Ciclos Límite
Objetivo: Descubrir atractores periódicos
- Selecciona Van der Pol. Este es un oscilador no lineal: dx/dt = y, dy/dt = μ(1-x²)y - x.
- Observa cómo las partículas cerca del origen son repelidas pero las lejanas son atraídas hacia una órbita cerrada.
- Añade partículas tanto cerca del origen como lejos. Todas convergen al mismo ciclo límite.
- Reduce la velocidad a 0.5x para apreciar mejor la dinámica. El ciclo no es un círculo perfecto debido a la no linealidad.
- Este ciclo límite es un atractor: trayectorias cercanas convergen hacia él desde ambos lados.
Conclusión: Los sistemas no lineales pueden tener comportamientos que no existen en sistemas lineales, como los ciclos límite estables.
Experimento 4: Lotka-Volterra (Depredador-Presa)
Objetivo: Modelar poblaciones en competencia
- Selecciona Lotka-Volterra. Las variables x e y representan poblaciones de presas y depredadores.
- Las ecuaciones son dx/dt = x(1 - 0.1y), dy/dt = y(-0.5 + 0.05x). La presa crece pero es limitada por depredadores; el depredador muere sin presas.
- Observa las órbitas cerradas. A diferencia de Van der Pol, aquí hay infinitas órbitas según la condición inicial.
- Añade partículas en diferentes posiciones. Cada una sigue su propia órbita periódica.
- Activa "Cuadrícula" para ver las escalas. Las órbitas representan ciclos de abundancia y escasez en la naturaleza.
Conclusión: Este sistema tiene un centro no lineal: infinitas órbitas cerradas que no son ni atractoras ni repulsoras. Es estructuralmente inestable.
Contexto Histórico
Henri Poincaré revolucionó el estudio de ecuaciones diferenciales al mostrar que el análisis cualitativo de sus soluciones era tan importante como las fórmulas explícitas. Su trabajo sobre el problema de tres cuerpos reveló que muchos sistemas no tienen soluciones analíticas, pero su comportamiento global puede entenderse geométricamente.
Poincaré: "Mémoire sur les courbes définies par une équation différentielle" introduce el análisis cualitativo y el retrato de fases.
Lotka: Propone ecuaciones para reacciones químicas autocatalíticas que luego Volterra aplica a poblaciones biológicas.
Van der Pol: Estudia oscilaciones en circuitos de válvulas de vacío, descubriendo el ciclo límite que lleva su nombre.
Lorenz: Descubre el atractor caótico, un nuevo tipo de comportamiento en sistemas de tres o más dimensiones.
Conexiones Interdisciplinarias
Para Explorar Más
- Strogatz, S. (2015) - "Nonlinear Dynamics and Chaos" - El texto de referencia sobre análisis cualitativo
- Hirsch, M., Smale, S. (1974) - "Differential Equations, Dynamical Systems, and Linear Algebra"
- Simulación: Atractor de Lorenz - Extensión a tres dimensiones
- Simulación: Péndulo Doble - Caos en sistemas mecánicos