🔄 Retrato de Fase
Visualiza el plano de fase de sistemas dinámicos 2D. Elige entre Lotka-Volterra, el oscilador de van der Pol, el péndulo, la bifurcación de Hopf y más. Haz clic para trazar trayectorias, superponer nulclinas y localizar puntos fijos.
Acerca de esta simulación
Esta simulación renderiza el plano de fase de un sistema dinámico 2D enteramente en la GPU usando un shader GLSL de estilo Line-Integral-Convolution (LIC): el ruido deformado por el dominio se difumina a lo largo de la dirección local del campo vectorial, produciendo líneas de flujo continuas en lugar de flechas discretas. El color codifica la velocidad local |(f, g)|, y un brillo suave marca los puntos fijos donde la velocidad se anula. Hay cinco sistemas incorporados — una matriz lineal 2×2 ajustable, el oscilador de van der Pol, un péndulo amortiguado, el modelo depredador-presa de Lotka-Volterra y el oscilador de doble pozo de Duffing — cada uno revela un tipo distinto de comportamiento a largo plazo, desde espirales y sillas hasta un ciclo límite genuinamente no lineal.
🔬 Qué muestra
El campo de velocidad (dx/dt, dy/dt) = (f(x,y), g(x,y)) evaluado en cada píxel, visualizado como líneas de flujo animadas cuyo color pasa de azul→cian→amarillo→rojo con velocidad creciente. Los puntos fijos, donde el campo es cero, brillan suavemente. Para el sistema lineal, una lectura en vivo clasifica el equilibrio — silla, nodo, espiral o centro — directamente a partir de la traza y el determinante de la matriz.
🎮 Cómo usarlo
Elige un sistema de la lista y luego arrastra Parámetro a y Parámetro b para remodelar su dinámica: en el caso lineal desplazan el equilibrio entre espirales estables/inestables, nodos, una silla y un centro; para van der Pol, el parámetro a es el amortiguamiento no lineal μ; para el péndulo es el amortiguamiento y el par motor; para Lotka-Volterra son las tasas de crecimiento de presas y muerte de depredadores; para Duffing es el amortiguamiento. La Velocidad de flujo controla la rapidez con la que se animan las líneas de flujo, y Reiniciar restaura los valores por defecto del sistema elegido.
💡 ¿Sabías que…?
El oscilador de van der Pol, inventado en la década de 1920 para modelar circuitos de válvulas de vacío, se asienta en un ciclo límite — un único bucle cerrado hacia el que convergen en espiral las trayectorias tanto desde dentro como desde fuera. Ningún sistema lineal, por mucho que se ajuste su matriz, puede producir este comportamiento: los ciclos límite son un fenómeno genuinamente no lineal, y la ecuación de van der Pol se convirtió en uno de los ejemplos fundacionales en el estudio de las oscilaciones autosostenidas.
Preguntas frecuentes
¿Qué es un retrato de fase y qué representan las líneas de flujo?
Un retrato de fase es una imagen del espacio de estados de un sistema dinámico que muestra cómo evoluciona con el tiempo cada punto de partida posible. En esta simulación, cada punto del plano tiene un vector de velocidad (dx/dt, dy/dt) determinado por las ecuaciones del sistema elegido, y las líneas de flujo son una representación de tipo Line-Integral-Convolution de ese campo vectorial — ruido difuminado a lo largo de la dirección local del flujo — de modo que puedes ver el patrón general de movimiento de un vistazo en lugar de seguir una trayectoria a la vez.
¿Cómo se clasifica un punto fijo como silla, nodo, espiral o centro?
Para el sistema lineal con matriz [[a, b], [-b, a]], la traza es T = 2a y el determinante es D = a² + b². Si D < 0 el origen es una silla (inestable en una dirección, estable en otra); si D > 0 y T² < 4D es una espiral, estable cuando T < 0 e inestable cuando T > 0; si D > 0 y T² > 4D es un nodo; y si T = 0 con D > 0 es un centro, rodeado de órbitas cerradas. El panel de lectura calcula esta clasificación en vivo a medida que mueves los deslizadores.
¿Qué diferencia al oscilador de van der Pol de un oscilador amortiguado simple?
Un oscilador lineal amortiguado estándar siempre converge en espiral hacia un único punto de equilibrio y permanece allí. La ecuación de van der Pol ẋ = y, ẏ = μ(1−x²)y − x tiene un amortiguamiento que cambia de signo según la amplitud: energiza las oscilaciones pequeñas y amortigua las grandes, de modo que las trayectorias tanto desde dentro como desde fuera convergen hacia la misma órbita cerrada, llamada ciclo límite. Esta amplitud de oscilación autocorregida es la razón por la que el modelo de van der Pol se usa para describir latidos del corazón, circuitos de válvulas de vacío y otros osciladores autosostenidos.
¿Qué representan las órbitas cerradas en el sistema de Lotka-Volterra?
Las ecuaciones de Lotka-Volterra ẋ = a·x − x·y, ẏ = x·y − b·y modelan un par depredador-presa: las presas (x) crecen exponencialmente a una tasa a en ausencia de depredadores, y los depredadores (y) mueren a una tasa b sin presas que comer. Las órbitas cerradas que rodean el punto de coexistencia (b, a) significan que las poblaciones oscilan periódicamente en lugar de estabilizarse en un valor fijo — un aumento de presas alimenta un aumento de depredadores, que luego agota a las presas, lo que hace morir de hambre a los depredadores, y el ciclo se repite, sin converger nunca al equilibrio.
¿Por qué el oscilador de Duffing tiene dos espirales estables separadas?
El sistema de Duffing ẋ = y, ẏ = −b·y + x − x³ describe el movimiento en un potencial de doble pozo, físicamente como una bola que puede asentarse en cualquiera de dos hondonadas separadas por una joroba. El término cúbico −x³ crea dos equilibrios estables simétricos en x = ±1 (cada uno una espiral cuando el amortiguamiento b > 0) y una silla inestable en el origen que separa sus cuencas de atracción, de modo que en qué espiral termina una trayectoria depende sensiblemente de su posición inicial respecto a la silla.
Ajusta los parámetros de un sistema dinámico 2D y observa cómo las trayectorias del campo vectorial giran en espiral, orbitan o se alejan de los puntos fijos.
3D · Renderizador Three.js / WebGL · Objetivo de 60 FPS · funciona completamente en el navegador, sin instalación