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

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

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

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

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

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

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

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

О модели

Для линейного тестового уравнения 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 решение может быть устойчивым, но неточным; точность всё равно определяется ошибкой аппроксимации и гладкостью задачи.

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

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

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

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

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

Одно и то же нелинейное уравнение маятника и шаг h; прямое сравнение методов Эйлера и РК4.

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

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

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

N′ = αN−βNP, P′ = δNP−γP; фазовая плоскость RK4; точка равновесия.

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

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

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

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

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