Для линейного тестового уравнения y′ = λy любой одношаговый метод с фиксированным шагом даёт переход yₙ₊₁ = R(z)yₙ, где z = hλ. Метод абсолютно устойчив в этой точке, если повторные шаги не усиливают моду: |R(z)| ≤ 1. Симулятор закрашивает области устойчивости на комплексной плоскости z для явного Эйлера (R = 1 + z), RK2 (R = 1 + z + z²/2) и классического RK4 (R = 1 + z + z²/2 + z³/6 + z⁴/24). Ползунки двигают λ и h, поэтому оранжевая точка z = hλ перемещается по карте; розовая точка умножает вещественную часть на коэффициент жёсткости и показывает, как быстрая затухающая мода может выйти из области устойчивости, даже если медленная мода выглядит безопасной. Это геометрическая связка между формулами Рунге–Кутты, коэффициентами усиления и ограничениями шага в жёстких ОДУ и полу-дискретизациях PDE.
Для кого: Студенты численных методов, scientific computing и вычислительной физики, изучающие функции устойчивости, ограничения шага и трудности явных RK-методов на жёстких задачах.
Ключевые понятия
Метод Рунге–Кутты
Абсолютная устойчивость
Функция устойчивости
Метод Эйлера
RK2
RK4
Жёсткое ОДУ
Тестовое уравнение
Как это работает
Для тестового y′ = λy одношаговый метод умножает решение на R(z), где z = hλ. Закрашены области |R(z)| ≤ 1 для Euler, RK2, RK4; оранжевая точка — текущий hλ, розовая — более жёсткая мода. Если точка вне области, явный шаг будет усиливать затухающую моду.
Основные формулы
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.
Часто задаваемые вопросы
Почему используют тестовое уравнение y′ = λy?
Около стационарной точки многие системы линеаризуются на моды вида e^{λt}. Численный метод применяет к каждой такой моде свой множитель R(hλ), поэтому скалярный тест быстро показывает, будет ли затухающая непрерывная мода также затухать в дискретной схеме.
Почему RK4 не устойчив для всех отрицательных вещественных λ?
У явных RK-методов области устойчивости ограничены. На отрицательной вещественной оси RK4 доходит примерно до z ≈ −2,785. Если λ очень отрицательно, то hλ покидает область, если h не сделать малым. Это типичное ограничение явного шага для жёсткого затухания.
Попадание в область устойчивости гарантирует точность?
Нет. Устойчивость означает только, что дискретная мода не растёт искусственно. При слишком большом h решение может быть устойчивым, но неточным; точность всё равно определяется ошибкой аппроксимации и гладкостью задачи.