🌊 Диффуры🔢 Численные методы

Жёсткие и нежёсткие ОДУ: почему явный метод Эйлера ломается

Q: Как распознать жёсткую задачу, не решая её?

Посмотрите на якобиан правой части: если собственные значения лежат в левой полуплоскости, а их модули различаются в сотни раз, задача жёсткая. Косвенный тест дешевле счёта: решите задачу двумя методами с одинаковым шагом. Явный взорвался или зашёлся пилой, неявный спокоен — перед вами жёсткость, и лечится она переходом на неявные схемы.

Q: Почему неявные методы считаются дороже явных?

Явному хватает одного вычисления правой части на шаг. Неявный решает уравнение относительно нового значения: для нелинейной задачи это метод Ньютона — несколько вычислений $f$ и якобиана плюс решение системы линейных уравнений на каждой итерации. Плату отбивает шаг: его можно поднять на порядки, и суммарная стоимость счёта падает.

Q: Нужны ли жёсткие решатели всегда?

Нет. Если собственные значения якобиана одного порядка, шаг всё равно выбирается ради точности, и неявная схема лишь добавит возни с Ньютоном. Жёсткость становится проблемой на длинных интервалах, где условие $h < 2/|\lambda_{\text{макс}}|$ связывает руки именно тогда, когда решение давно вышло на плавную траекторию и быстрого процесса больше не видно.

Q: Как связаны шаг и точность при решении жёстких задач?

У неявного метода шаг ограничен только точностью, поэтому его выбирает контроль локальной ошибки: погрешность шага Эйлера порядка $h^2$. Рабочий приём — считать с шагом $h$ и $h/2$ и сравнивать решения: разность прикидывает ошибку. Если при укрупнении шага решение не портится, устойчивость вас больше не диктует — счёт вышел на точностный режим.

Жёсткие ОДУНежёсткие ОДУ

Коротко: Жёсткость — свойство задачи, а не метода; лечится переходом на неявные схемы. Явный берёт шаг из устойчивости (h < 2/|λ_макс|), неявный — из точности: A-устойчивость снимает запрет на шаг, L-устойчивость — следы быстрых мод. Быстрый диагноз: два метода, один шаг, разные ответы — задача жёсткая.

План на семестр Вердикт: с чего начать

Задача одной строкой: $y^{'} = - 1000 (y - cos t)$ , $y (0) = 1$ . Решение образцово-приличное: за первые несколько тысячных секунды оно прилипает к квазистационарной кривой $y \approx cos t + 0, 001 sin t$ и дальше просто следит за косинусом с микроскопическим запаздыванием. А теперь явный Эйлер с шагом $h = 0, 0021$ : на первой секунде лёгкая рябь, к пятой — числа порядка $1 0^{91}$ . Уравнение не меняли, метод не меняли — шаг перевалил через $0, 002$ , и вычисление сгорело.

Это и есть жёсткость: в задаче смешаны процессы сильно разных скоростей, и самый быстрый из них, давно умерев в решении, продолжает диктовать шаг явному методу. Ниже — вывод границы $h < 2/∣ λ ∣$ , определение жёсткости через собственные значения якобиана, честный эксперимент с двумя методами на двух шагах и практическая диагностика. Механика явного метода — в уроке метод Эйлера, фон по устойчивости — в уроке устойчивость решений, вычислительная кухня — в уроке численное решение ОДУ.

Критерий	Жёсткие ОДУ	Нежёсткие ОДУ
Собственные значения якобиана	$∣ λ_{макс} ∣/∣ λ_{мин} ∣ ≫ 1$ — сотни и тысячи	один порядок величин
Масштабы времени	быстрая мода умирает за $1/∣ λ_{макс} ∣$ , но навязывает шаг всему счёту	одно характерное время жизни
Явный Эйлер	шаг $h < 2/∣ λ_{макс} ∣$ — доли миллисекунды	шаг диктуется точностью
Неявный Эйлер	работает при любом $h$ (A-устойчивость)	тоже работает, но Ньютон не окупается
Симптом поломки явного метода	пила или взрыв при формально разумном $h$	просто крупная ошибка при крупном $h$
Лечение	неявные схемы, L-устойчивые методы Гира/BDF	явные методы Рунге—Кутты

Жёсткость — это разделение масштабов времени, а не свойство алгоритма

Почему явный Эйлер взрывается уже при h = 0,0021?#

Достаточно посмотреть на тестовое уравнение $y^{'} = λ (y - c)$ с $λ < 0$ . Один шаг явного Эйлера умножает отклонение от равновесия на фиксированный множитель — и ошибка обязана гаснуть, а не расти:

y_{n + 1} = y_{n} + hλ (y_{n} - c) = (1 + hλ) y_{n} - hλ c

множитель шага

1 + hλ

: при

∣1 + hλ ∣ > 1

ошибка растёт геометрически, при

∣1 + hλ ∣ < 1

— гаснет

∣ 1 + hλ ∣ \leq 1 ⟹ h \leq \frac{2}{∣ λ ∣}

условие устойчивости явного Эйлера: для

λ = - 1000

получаем

h \leq 0, 002

— шаг, продиктованный не точностью, а выживанием

Проверяем на нашем уравнении ( $λ = - 1000$ ), интегрируя до $t = 5$ : при $h = 0, 0019$ множитель $- 0, 9$ — решение ограничено; при $h = 0, 002$ множитель ровно $- 1$ — незатухающая пила; при $h = 0, 0021$ множитель $- 1, 1$ , и максимум решения за пять секунд достигает порядка $1 0^{91}$ . Двухпроцентное превышение шага — девяносто порядков на выходе. При этом неявный Эйлер на том же отрезке спокойно работает хоть с $h = 0, 1$ : его множитель $1/ (1 - hλ)$ по модулю всегда меньше единицы.

Как определить жёсткость через собственные значения якобиана?#

Для системы $y^{'} = f (t, y)$ матрица частных производных — якобиан — несёт собственные значения $λ_{1}, \dots, λ_{n}$ ; как их находят — в уроке собственные значения и собственные векторы. Задача жёсткая, если все $λ$ лежат в левой полуплоскости, а их модули разделены в разы:

s = \frac{max _{i} ∣ λ _{i} ∣}{min _{i} ∣ λ _{i} ∣} ≫ 1

показатель жёсткости: для системы с

λ_{1} = - 1

λ_{2} = - 1000

получаем

s = 1000

: быстрая мода затухает в тысячу раз быстрее медленной

Физика за формулой: каждой моде отвечает своя постоянная времени $1/∣ λ ∣$ — время затухания в $e$ раз (см. термин постоянная времени). Быстрая мода $e^{- 1000 t}$ умирает за пять тысячных секунды и дальше невооружённым глазом не видна, но явному методу она выставляет счёт за шаг на всём интервале. Решение тем временем прилипает к медленной траектории: во вступительной задаче медленный процесс — косинус, а быстрый — релаксация с $λ = - 1000$ ; частное решение вида $a cos t + b sin t$ для таких линейных уравнений строится в уроке линейные неоднородные ДУ.

Что показывают числа: как ведут себя оба метода на двух шагах?#

Эксперимент: $y^{'} = - 20 (y - cos t)$ , $y (0) = 0$ — сознательно стартуем не на квазистационарной кривой, даём начальный толчок величиной около единицы. Точное решение известно: $y (t) = \frac{400}{401} cos t + \frac{20}{401} sin t - \frac{400}{401} e^{- 20 t}$ . Гоняем явный и неявный Эйлер и меряем максимальную ошибку против него:

Метод	Шаг $h$	Множитель шага	max\|ошибка\| на $[0{,}5;\, 5]$	Что видит глаз
явный Эйлер	0,05	$∣1 - 20 h ∣ = 0$	≈ 0,001	толчок убит за один шаг, решение гладкое
явный Эйлер	0,1	$∣1 - 20 h ∣ = 1$	≈ 1,0	пила: ошибка не гаснет, решение дёргается с размахом около 2
явный Эйлер	0,11	$∣1 - 20 h ∣ = 1, 2$	≈ 3 700	ошибка умножается на 1,2 за шаг: $1, 2^{45} \approx 3, 7 \cdot 1 0^{3}$
неявный Эйлер	0,1	$1/ (1 + 20 h) = 1/3$	≈ 0,006	толчок гаснет втрое за шаг, дальше гладко
неявный Эйлер	0,05	$1/ (1 + 20 h) = 1/2$	≈ 0,002	то же, вдвое точнее

Сетка

t \in [0; 5]

; значения экспериментальные и зависят от сетки и начального условия, но картина воспроизводится всегда

Три симптома в одной таблице. Шаг $h = 0, 05$ : множитель нулевой — начальный толчок истреблён мгновенно, живёт только ошибка аппроксимации. Шаг $h = 0, 1$ : ровно граница устойчивости, $∣1 - 20 h ∣ = 1$ — толчок вспоминается бесконечно пилой, хотя за 0,5 секунды он в решении давно умер. Шаг $h = 0, 11$ : рост в 1,2 раза за шаг, и через 45 шагов ошибка — тысячи. Неявный метод при любом из этих шагов спокоен: его множитель $1/ (1 + 20 h)$ по модулю меньше единицы всегда, а при большом $h$ близок к нулю — сверхбыстрые моды он не гасит, а хоронит.

Тот же счёт для вступительной задачи с $λ = - 1000$ : явный жив при $h < 0, 002$ , у границы пилит, за границей улетает в $1 0^{91}$ ; неявный при $h = 0, 01$ ведёт решение с максимальной ошибкой порядка $5 \cdot 1 0^{- 6}$ — шаг в пять раз крупнее всего допустимого бюджета явного метода. Вот цена вопроса: не «точнее/грубее», а «возможно/невозможно».

Чем неявная схема платит за устойчивость: что такое A- и L-устойчивость?#

Неявный Эйлер подставляет правую часть в неизвестную точку:

y_{n + 1} = y_{n} + h f (t_{n + 1}, y_{n + 1})

y_{n + 1}

сидит в обеих частях уравнения — за устойчивость платят решением уравнения на каждом шаге

Для нашей задачи уравнение линейное, и шаг сводится к одной строке: $y_{n + 1} = (y_{n} + 20 h cos t_{n + 1}) / (1 + 20 h)$ . Для общей нелинейной $f$ на каждом шаге решают нелинейное уравнение методом Ньютона — итерационным методом с двумя–пятью вычислениями $f$ и якобиана, а для систем — ещё и системой линейных уравнений (техника — в уроке численные методы для СЛАУ). Шаг дорожает в разы; окупается тем, что шаг можно поднять на порядки.

Классификация устойчивости короткая. A-устойчивость: область устойчивости содержит всю левую полуплоскость — любой шаг годится, пока $Re λ < 0$ ; неявный Эйлер ею обладает. L-устойчивость — усиление: множитель метода стремится к нулю при $hλ \to - \infty$ , то есть бесконечно быстрые моды уничтожаются за один шаг:

R (z) = \frac{1}{1 - z}, R (z) z \to - \infty 0

функция устойчивости неявного Эйлера: A-устойчив и L-устойчив; у метода трапеций

R (\infty) = - 1

— A-устойчив, но пилу быстрых мод не хоронит

Разница между A и L не академична: трапеции (Кранк—Николсон) при крупном шаге на жёсткой задаче дают незатухающую рябь с множителем $- 1$ — устойчиво по букве закона и некрасиво по существу. Поэтому в жёстких решателях стоят именно L-устойчивые схемы: неявный Эйлер, методы Гира (BDF), неявные Рунге—Кутты типа Radau.

Как диагностировать жёсткость на практике?#

Шаг 1: решите задачу двумя методами с одинаковым шагом — явным и неявным. Явный взорвался или зашёлся пилой, а неявный спокоен? Перед вами жёсткая задача.
Шаг 2: выпишите якобиан и посчитайте собственные значения: отношение максимального модуля к минимальному больше сотни — готовьтесь к неявным схемам.
Шаг 3: оцените бюджет явного счёта: сколько шагов продиктует условие $h < 2/∣ λ_{макс} ∣$ на всём интервале. Миллионы однотипных шагов — время менять метод.
Шаг 4: проверьте, нужна ли быстрая мода вообще: если процесс умирает за $1/∣ λ_{макс} ∣$ и дальше не виден, L-устойчивая схема похоронит его честно.
Шаг 5: в готовых пакетах жёсткость включается выбором решателя (ode15s, Radau, BDF) — но диагноз «жёсткая/нежёсткая» всё равно остаётся за вами: пакет угадает не всегда.

Свойство	Явный Эйлер	Неявный Эйлер
Формула шага	$y_{n + 1} = y_{n} + h f (t_{n}, y_{n})$	$y_{n + 1} = y_{n} + h f (t_{n + 1}, y_{n + 1})$
Устойчивость	$h \leq 2/∣ λ ∣$ при $λ < 0$ — шаг под диктовку быстрой моды	A-устойчив: любой $h$ при $Re λ < 0$
Быстрые моды	навязывают шаг всему счёту	L-устойчивость: гасятся за один шаг
Стоимость шага	одно вычисление $f$	линейно — одна формула; в общем случае Ньютон: $f$ , якобиан, СЛАУ
Шаг ограничен	устойчивостью	только точностью
Когда применять	нежёсткие задачи, дешёвая $f$ , короткие интервалы	жёсткие задачи, длинные интервалы, кинетика, цепи

Явный и неявный Эйлер: цена шага против цены устойчивости

Проверь себя+15 XP

$y^{'} = - 1000 (y - cos t)$ . При каких шагах явный Эйлер останется устойчивым?

Проверь себя+15 XP

На задаче $y^{'} = - 20 (y - cos t)$ , $y (0) = 0$ : явный Эйлер при $h = 0, 1$ выдаёт пилу с размахом около 2, неявный при том же шаге — гладкое решение. Диагноз?

Частые вопросы

Как распознать жёсткую задачу, не решая её?

Посмотрите на якобиан правой части: если собственные значения лежат в левой полуплоскости, а их модули различаются в сотни раз, задача жёсткая. Косвенный тест дешевле счёта: решите задачу двумя методами с одинаковым шагом. Явный взорвался или зашёлся пилой, неявный спокоен — перед вами жёсткость, и лечится она переходом на неявные схемы.

Почему неявные методы считаются дороже явных?

Явному хватает одного вычисления правой части на шаг. Неявный решает уравнение относительно нового значения: для нелинейной задачи это метод Ньютона — несколько вычислений $f$ и якобиана плюс решение системы линейных уравнений на каждой итерации. Плату отбивает шаг: его можно поднять на порядки, и суммарная стоимость счёта падает.

Нужны ли жёсткие решатели всегда?

Нет. Если собственные значения якобиана одного порядка, шаг всё равно выбирается ради точности, и неявная схема лишь добавит возни с Ньютоном. Жёсткость становится проблемой на длинных интервалах, где условие $h < 2/∣ λ_{макс} ∣$ связывает руки именно тогда, когда решение давно вышло на плавную траекторию и быстрого процесса больше не видно.

Как связаны шаг и точность при решении жёстких задач?

У неявного метода шаг ограничен только точностью, поэтому его выбирает контроль локальной ошибки: погрешность шага Эйлера порядка $h^{2}$ . Рабочий приём — считать с шагом $h$ и $h /2$ и сравнивать решения: разность прикидывает ошибку. Если при укрупнении шага решение не портится, устойчивость вас больше не диктует — счёт вышел на точностный режим.

Следующий шаг: прогоните эксперимент из таблицы руками — пять строк кода на явный и неявный Эйлер, и граница устойчивости перестанет быть абстракцией. Порядок действий при решении уравнений — в разборе, как решать ДУ, задачи для тренировки — в задачах по дифференциальным уравнениям, счёт — в тренажёре по дифурам; формулы под рукой — в шпаргалке по ДУ. Ошибки аппроксимации и округления, которые добавляются поверх устойчивости, разобраны в уроке источники погрешностей; сравнение двух явных методов — в сравнении РК4 против Эйлера, а исходные определения — в терминах дифференциальное уравнение и в уроке системы ДУ.