Жёсткие системы

Симуляция реакции горения: тысячи химических компонентов с временными масштабами от наносекунд до секунд. Явный RK4 требует 10⁹ шагов. BDF/Radau решают ту же задачу за сотни шагов - разница в несколько порядков времени вычисления.

• **Химическая кинетика:** реакции горения, синтез полимеров - жёсткие системы с временными масштабами от 10⁻¹² до 10³ секунд; LSODA/BDF стандарт
• **Энергосистемы:** симуляция электросети - жёсткая система с быстрыми переходными процессами и медленной динамикой; неявные методы обязательны
• **Жёсткие Neural ODE:** некоторые задачи управления и физики порождают жёсткие НейронОДУ; adjoint+BDF в torchdiffeq

Предварительные знания

• Numerical Methods for ODEs

Что такое жёсткость

Жёсткая система ОДУ - система, в которой присутствуют компоненты, изменяющиеся с существенно разными скоростями. Формально: отношение наибольшего и наименьшего собственных значений матрицы Якоби большое. Для явных методов это требует очень малого шага ради устойчивости, даже если точность этого не требует.

**Число жёсткости:** stiffness ratio = max|Re(λᵢ)| / min|Re(λᵢ)| где λᵢ - собственные значения Якоби f(t, y) по y. **Устойчивость явного Эйлера:** метод устойчив при h·|λ| ≤ 2 Для жёсткой системы с λ_max = 10⁶: - Явный Эйлер требует h < 2·10⁻⁶ (миллионы шагов) - Неявный Эйлер: любой h устойчив (A-устойчивость) **Пример:** dy/dt = −1000y + cos(t)

Типичные признаки жёсткости: solve_ivp с RK45 делает тысячи маленьких шагов, скорость интегрирования неожиданно низкая. Решение: использовать method='Radau' или method='BDF' в solve_ivp.

Неявный Эйлер и A-устойчивость

Неявный Эйлер вычисляет yₙ₊₁ через **уравнение на новом шаге**, а не через известный xₙ. Это требует решения нелинейной системы (методом Ньютона) на каждом шаге, но даёт A-устойчивость: метод устойчив при любом h > 0.

**Неявный Эйлер:** yₙ₊₁ = yₙ + h·f(tₙ₊₁, yₙ₊₁) На каждом шаге решаем: G(yₙ₊₁) = yₙ₊₁ − yₙ − h·f(tₙ₊₁, yₙ₊₁) = 0 → Метод Ньютона: O(n²) или O(n³) на шаг **A-устойчивость:** метод устойчив для всех λ с Re(λ) < 0 при любом h Область устойчивости явного Эйлера: |1 + hλ| ≤ 1 Область неявного: |1/(1 − hλ)| ≤ 1 → всё левое полуплоскость!

BDF методы и LSODA

BDF (Backward Differentiation Formula) - семейство многошаговых неявных методов для жёстких ОДУ. Они используют несколько предыдущих значений yₙ, yₙ₋₁, ... для более высокой точности. LSODA (VODE) - адаптивный солвер, автоматически переключающийся между жёсткой и нежёсткой стратегиями.

**BDF порядка k:** Σⱼ₌₀ᵏ aⱼ yₙ₊₁₋ⱼ = h·β₀·f(tₙ₊₁, yₙ₊₁) BDF1 = Неявный Эйлер: O(h) BDF2: O(h²), BDF4: O(h⁴), BDF6: O(h⁶) **Ограничение:** BDF порядка > 6 нестабильны (теорема Дальквиста) **LSODA:** автоматически выбирает Adams (нежёсткий) или BDF (жёсткий) метод, определяя жёсткость на лету. `solve_ivp(f, ..., method='LSODA')` или `method='BDF'`

Ключевые идеи

• **Жёсткость:** большое соотношение max|λ|/min|λ| → явные методы требуют h → 0 для устойчивости
• **Неявный Эйлер:** A-устойчив (любой h); дорогой шаг (Ньютон для нелинейных); O(h) точность
• **BDF методы:** неявные многошаговые; порядки 1-6; теорема Дальквиста ограничивает > 6
• **Практика:** method='Radau' для жёстких задач; method='LSODA' когда жёсткость неизвестна

Связанные темы

Жёсткость связывает ОДУ с линейной алгеброй и приложениями:

• Численные методы для ОДУ — Явные методы (RK45) - стартовая точка; жёсткость объясняет, почему они неэффективны в некоторых задачах
• Метод конечных разностей для PDE — Дискретизация PDE порождает жёсткие ОДУ-системы; CFL-условие - условие устойчивости явных методов

Вопросы для размышления

• Как определить, является ли система ОДУ жёсткой, без явного вычисления собственных значений Якобиана?
• Почему для Radau IIA (неявный Рунге-Кутта) нет ограничения Дальквиста, несмотря на высокий порядок?
• Что происходит с Neural ODE при жёстких динамиках: как это влияет на обучение?

Связанные уроки

• la-01-vectors-intro