Guía: Ecuaciones Diferenciales

Introducción

Una ecuación diferencial es una ecuación que relaciona una función con sus derivadas. Son el lenguaje matemático del cambio: desde el crecimiento de poblaciones hasta el movimiento de planetas, las leyes de la física se expresan naturalmente como EDOs.

Esta simulación visualiza campos de pendientes (slope fields): en cada punto del plano, una pequeña línea indica la pendiente que tendría cualquier solución que pase por ahí. Haz clic para trazar curvas solución que fluyen siguiendo estas pendientes.

Notación

dy/dx = f(x, y) - EDO de primer orden
d²y/dx² = g(x, y, y') - EDO de segundo orden
y(x₀) = y₀ - Condición inicial
Isoclina - Curva donde dy/dx = constante
Nulclina - Curva donde dy/dx = 0

Ecuaciones Disponibles

Crecimiento Exponencial

dy/dx = k·y

La tasa de cambio es proporcional al valor actual. Modela poblaciones ilimitadas, decaimiento radiactivo (k < 0), interés compuesto continuo.

Solución analítica: y = Ce^(kx)

Crecimiento Logístico

dy/dx = r·y·(1 - y/K)

Crecimiento con capacidad de carga K. La población crece rápido cuando es pequeña, se frena al acercarse a K. Más realista que el exponencial.

Solución analítica: y = K / (1 + Ce^(-rx)). Punto de equilibrio en y = K.

Oscilador Armónico

d²y/dx² = -ω²·y

Sistema masa-resorte sin fricción. Produce oscilaciones sinusoidales perpetuas. El retrato de fase muestra órbitas cerradas (elipses).

Solución: y = A·cos(ωx) + B·sin(ωx). Energía se conserva.

Oscilador Van der Pol

d²y/dx² = μ(1-y²)·dy/dx - y

Oscilador no lineal con ciclo límite. Todas las soluciones (excepto el origen) convergen a una órbita periódica estable. Modela latidos cardíacos, circuitos electrónicos.

No tiene solución analítica cerrada. μ controla la no linealidad.

Lotka-Volterra

dx/dt = αx - βxy, dy/dt = δxy - γy

Modelo depredador-presa. x = presas, y = depredadores. Produce oscilaciones acopladas: cuando hay muchas presas, los depredadores crecen, reducen las presas, y luego decaen.

Las órbitas son cerradas pero no elipses. Punto de equilibrio en (γ/δ, α/β).

Péndulo No Lineal

d²θ/dt² = -(g/L)·sin(θ)

Péndulo real (no la aproximación sin(θ) ≈ θ). El retrato de fase muestra separatrices: fronteras entre oscilación y rotación completa.

Solución involucra integrales elípticas. Periodo depende de la amplitud.

Métodos Numéricos

Método de Euler

y_{n+1} = y_n + h·f(x_n, y_n)

El método más simple. Usa la pendiente actual para dar un "paso" de tamaño h. Error local O(h²), error global O(h). Puede ser inestable para h grande.

Runge-Kutta 4 (RK4)

k₁ = f(x, y)
k₂ = f(x+h/2, y+hk₁/2)
k₃ = f(x+h/2, y+hk₂/2)
k₄ = f(x+h, y+hk₃)
y_{n+1} = y_n + (h/6)(k₁+2k₂+2k₃+k₄)

Promedia 4 estimaciones de pendiente. Error local O(h⁵). Mucho más preciso que Euler para el mismo h. Método estándar en ciencia.

Experimentos Sugeridos

1. Comparar Euler vs RK4

Con crecimiento exponencial (k = 0.5), traza una solución con cada método y el mismo punto inicial.

Euler diverge ligeramente de la solución exacta
RK4 es indistinguible de la analítica
La diferencia crece con el tiempo

2. Estabilidad Logística

En la ecuación logística, traza soluciones desde diferentes condiciones iniciales.

Todas convergen a y = K (capacidad de carga)
y = 0 es equilibrio inestable
y = K es equilibrio estable (atractor)

3. Ciclo Límite de Van der Pol

Traza soluciones desde dentro y fuera del ciclo límite.

Soluciones interiores espiralan hacia afuera
Soluciones exteriores espiralan hacia adentro
Todas convergen al mismo ciclo cerrado
Aumenta μ para ver oscilaciones de "relajación"

4. Separatrices del Péndulo

En el péndulo, busca las curvas que separan regiones de comportamiento diferente.

Órbitas cerradas pequeñas: oscilación
Órbitas que dan vueltas completas: rotación
La separatriz pasa por (±π, 0): puntos de silla

5. Desfase Depredador-Presa

En Lotka-Volterra, observa cómo oscilan presas y depredadores.

Los picos de depredadores siguen a los de presas
El desfase es típicamente de ~90° (cuarto de periodo)
Las órbitas son cerradas: sistema conservativo

Conceptos Clave

Puntos de Equilibrio

Puntos donde dy/dx = 0 (y dv/dx = 0 en 2D). Las soluciones constantes. Pueden ser estables (atraen soluciones cercanas) o inestables (las repelen).

Retrato de Fase

Para sistemas de 2º orden, graficamos (y, dy/dx) en lugar de (x, y). Muestra la dinámica cualitativa: equilibrios, ciclos, separatrices, atractores.

Existencia y Unicidad

Si f(x, y) es continua y tiene derivada parcial ∂f/∂y acotada, entonces existe una única solución por cada condición inicial. Las soluciones no se cruzan (excepto en puntos de equilibrio).

Sistemas Autónomos

Cuando la ecuación no depende explícitamente de x: dy/dx = f(y). El campo de pendientes es el mismo a lo largo de líneas verticales. Las soluciones son "traslaciones" horizontales entre sí.

Conexiones

M Lorenz - Sistema caótico 3D F Péndulo Doble - Caos determinista B Ecosistemas - Depredador-Presa M Funciones - Derivadas visualizadas

Limitaciones del Modelo

⚠ Los métodos numéricos acumulan error con el tiempo. Para simulaciones largas, la solución puede desviarse significativamente.
⚠ Cerca de singularidades (donde f(x,y) explota), los métodos pueden fallar o dar resultados incorrectos.
⚠ Las isoclinas y nulclinas se calculan por muestreo discreto, pueden aparecer discontinuas o con "ruido".
⚠ Para Lotka-Volterra, los parámetros γ y δ están fijos internamente (γ = δ = 0.5) para simplificar la interfaz.

Contexto Histórico

Las ecuaciones diferenciales nacieron con el cálculo de Newton y Leibniz (c. 1680). Newton las usó para describir el movimiento planetario y la gravedad.

Euler (1768) desarrolló el primer método numérico sistemático. Runge y Kutta (1900-1901) mejoraron dramáticamente la precisión con sus métodos de orden superior.

Poincaré (c. 1890) revolucionó el campo al estudiar cualitativamente las soluciones sin resolverlas explícitamente, sentando las bases de la teoría de sistemas dinámicos y el caos.

El oscilador de Van der Pol (1920) fue uno de los primeros ejemplos estudiados de ciclo límite, surgido del análisis de circuitos de radio.