Почему используют тестовое уравнение y′ = λy?

Около стационарной точки многие системы линеаризуются на моды вида e^{λt}. Численный метод применяет к каждой такой моде свой множитель R(hλ), поэтому скалярный тест быстро показывает, будет ли затухающая непрерывная мода также затухать в дискретной схеме.

Почему RK4 не устойчив для всех отрицательных вещественных λ?

У явных методов Рунге–Кутты области устойчивости ограничены. На отрицательной вещественной оси RK4 доходит примерно до z ≈ −2,785. Если λ очень отрицательно, то hλ покидает область, если h не сделать малым. Это типичное ограничение явного шага для жёсткого затухания.

Попадание в область устойчивости гарантирует точность?

Нет. Устойчивость означает только, что дискретная мода не растёт искусственно. При слишком большом h решение может быть устойчивым, но неточным; точность всё равно определяется ошибкой аппроксимации и гладкостью задачи.

Области устойчивости Рунге–Кутты

Абсолютная устойчивость в плоскости z = hλ для явного Эйлера, RK2 и RK4: двигайте λ и h, сравнивайте |R(z)| и смотрите, почему жёсткие моды требуют малых явных шагов.

О модели

Для линейного тестового уравнения y′ = λy одношаговый метод с фиксированным шагом даёт переход yₙ₊₁ = R(z) yₙ, где z = hλ. Метод абсолютно устойчив в точке z, если повторные шаги не усиливают моду: |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λ перемещается по карте; вторая точка умножает вещественную часть на коэффициент жёсткости и показывает, как быстрая затухающая мода может выйти из области, даже если медленная мода ещё безопасна. Это связь формул Рунге–Кутты, коэффициента усиления и ограничения шага в жёстких ОДУ и полудискретизациях уравнений в частных производных.

Для кого: Численные методы: устойчивость, жёсткие ОДУ, явные схемы Рунге–Кутты.

Ключевые понятия

метод Рунге–Кутты
абсолютная устойчивость
функция устойчивости
метод Эйлера
RK2
RK4
жёсткое ОДУ
тестовое уравнение

Как это работает

Для тестового y′ = λy одношаговый метод умножает решение на R(z), где z = hλ. Закрашены области |R(z)| ≤ 1 для Эйлера, 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 не устойчив для всех отрицательных вещественных λ?: У явных методов Рунге–Кутты области устойчивости ограничены. На отрицательной вещественной оси RK4 доходит примерно до z ≈ −2,785. Если λ очень отрицательно, то hλ покидает область, если h не сделать малым. Это типичное ограничение явного шага для жёсткого затухания.
Попадание в область устойчивости гарантирует точность?: Нет. Устойчивость означает только, что дискретная мода не растёт искусственно. При слишком большом h решение может быть устойчивым, но неточным; точность всё равно определяется ошибкой аппроксимации и гладкостью задачи.

Продолжите с близкими темами в этой категории — или все 85 в «Визуализация математики».

Вся категория →

НовоеУниверситет / научные

Эйлер и РК4 (математический маятник)

Запустить симулятор

Одно нелинейное уравнение математического маятника и один шаг h: явный метод Эйлера рядом с методом Рунге–Кутты четвёртого порядка.

НовоеШкольные

Лотка–Вольтерра

Запустить симулятор

Классическая модель хищник–жертва: N′ = αN − βNP, P′ = δNP − γP. Фазовая плоскость, интегрирование Рунге–Кутты, точка равновесия (γ/δ, α/β).

НовоеУниверситет / научные

Уравнение теплопроводности: конечные разности

Запустить симулятор

Одномерное u_t = αu_xx: явная FTCS и неявная схема Эйлера назад; меняйте r = αΔt/Δx², наблюдайте взрыв явной схемы при r > 1/2 и сравнивайте численную диффузию.

НовоеУниверситет / научные