¿Por qué se usa la ecuación de prueba y′ = λy?

Cerca de un punto estacionario muchos sistemas se linealizan en modos de la forma e^{λt}. El método numérico aplica a cada modo su factor R(hλ), así que el test escalar muestra rápido si un modo continuo que decae también decae en el esquema discreto.

¿Por qué RK4 no es estable para todo λ real negativo?

En los métodos explícitos de Runge-Kutta las regiones de estabilidad son acotadas. Sobre el eje real negativo RK4 llega más o menos hasta z ≈ −2,785. Si λ es muy negativo, hλ sale de la región si no se achica h. Es el tope típico del paso explícito ante un decaimiento rígido.

¿Estar en la región de estabilidad garantiza exactitud?

No. Estabilidad solo quiere decir que el modo discreto no crece de forma artificial. Con un h demasiado grande la solución puede ser estable pero inexacta; la exactitud la siguen fijando el error de aproximación y la suavidad del problema.

Regiones de estabilidad de Runge-Kutta

Estabilidad absoluta en el plano z = hλ para Euler explícito, RK2 y RK4: se mueven λ y h, se comparan |R(z)| y se ve por qué los modos rígidos piden pasos explícitos chicos.

Sobre el modelo

Las regiones de estabilidad de Runge-Kutta para la ecuación de prueba lineal y′ = λy salen del avance de un paso fijo: yₙ₊₁ = R(z) yₙ, con z = hλ. El método es absolutamente estable en z si los pasos repetidos no amplifican el modo: |R(z)| ≤ 1. En el plano complejo z se pintan las regiones de Euler explícito (R = 1 + z), RK2 (R = 1 + z + z²/2) y el RK4 clásico (R = 1 + z + z²/2 + z³/6 + z⁴/24). Los controles mueven λ y h, así que el punto z = hλ recorre el mapa; un segundo punto multiplica la parte real por el coeficiente de rigidez y muestra cómo un modo rápido que decae puede salir de la región aunque el modo lento siga seguro. Es el vínculo entre las fórmulas de Runge-Kutta, el factor de amplificación y el tope de paso en EDO rígidas y semidiscretizaciones de ecuaciones en derivadas parciales.

Para quién: Métodos numéricos: estabilidad, EDO rígidas, esquemas explícitos de Runge-Kutta.

Conceptos clave

método de Runge-Kutta
estabilidad absoluta
función de estabilidad
método de Euler
RK2
RK4
EDO rígida
ecuación de prueba

Cómo funciona

Para la ecuación de prueba y′ = λy el método de un paso multiplica la solución por R(z), con z = hλ. Se pintan las regiones |R(z)| ≤ 1 de Euler, RK2 y RK4; el punto naranja es el hλ actual, el rosa un modo más rígido. Si el punto queda fuera de la región, el paso explícito amplificará un modo que en continuo decae.

Fórmulas principales

Test equation y′ = λy, z = hλ. Stability functions: R_E = 1+z, R_RK2 = 1+z+z²/2, R_RK4 = 1+z+z²/2+z³/6+z⁴/24. Stable when |R(z)| ≤ 1.

Preguntas frecuentes

¿Por qué se usa la ecuación de prueba y′ = λy?: Cerca de un punto estacionario muchos sistemas se linealizan en modos de la forma e^{λt}. El método numérico aplica a cada modo su factor R(hλ), así que el test escalar muestra rápido si un modo continuo que decae también decae en el esquema discreto.
¿Por qué RK4 no es estable para todo λ real negativo?: En los métodos explícitos de Runge-Kutta las regiones de estabilidad son acotadas. Sobre el eje real negativo RK4 llega más o menos hasta z ≈ −2,785. Si λ es muy negativo, hλ sale de la región si no se achica h. Es el tope típico del paso explícito ante un decaimiento rígido.
¿Estar en la región de estabilidad garantiza exactitud?: No. Estabilidad solo quiere decir que el modo discreto no crece de forma artificial. Con un h demasiado grande la solución puede ser estable pero inexacta; la exactitud la siguen fijando el error de aproximación y la suavidad del problema.

Siga con temas cercanos de esta categoría — o todos los 85 de «Visualización matemática».

Toda la categoría →

NuevoUniversidad / investigación

Euler y RK4 (péndulo simple)

Abrir el simulador

Una misma ecuación no lineal del péndulo simple y un mismo paso h: el método de Euler explícito junto al de Runge–Kutta de cuarto orden.

NuevoEscuela

Lotka-Volterra

Abrir el simulador

Modelo clásico depredador–presa: N′ = αN − βNP, P′ = δNP − γP. Plano de fases, integración de Runge–Kutta, punto de equilibrio (γ/δ, α/β).

NuevoUniversidad / investigación

Ecuación del calor: diferencias finitas

Abrir el simulador

Ecuación unidimensional u_t = αu_xx: FTCS explícito y Euler implícito hacia atrás; se cambia r = αΔt/Δx², se ve el estallido del esquema explícito si r > 1/2 y se compara la difusión numérica.

NuevoUniversidad / investigación

Retrato de fase (EDO 2D)

Abrir el simulador

Campo de direcciones y trayectorias al hacer clic: armónico, amortiguamiento, silla, nodos, focos, péndulo (RK4).

NuevoEscuela

Modelo SIR de epidemia

Abrir el simulador

S + I + R = 1: βSI y γI; ℛ₀ ≈ β/γ, umbral de inmunidad colectiva 1 − 1/ℛ₀; gráfica en el tiempo con RK4.

NuevoEscuela

SEIR / SEIRS: epidemia

Abrir el simulador

S + E + I + R = 1: la latencia σE retrasa la infecciosidad, recuperación γI; opcional ωR → S; ℛ₀ = β/γ; gráfica RK4.

Regiones de estabilidad de Runge-Kutta

Sobre el modelo

Cómo funciona

Fórmulas principales

Preguntas frecuentes

Simuladores afines

Euler y RK4 (péndulo simple)

Lotka-Volterra

Ecuación del calor: diferencias finitas

Retrato de fase (EDO 2D)

Modelo SIR de epidemia

SEIR / SEIRS: epidemia

Regiones de estabilidad de Runge-Kutta

Sobre el modelo

Cómo funciona

Fórmulas principales

Preguntas frecuentes