Теорема Пеано и нелипшицев случай

Существование без единственности

Теорема Пикара требует условия Липшица. Что если это условие нарушено? Теорема Пеано отвечает на вопрос о существовании в более общем случае.

Теорема Пеано: Если f(x, y) непрерывна в прямоугольнике R = {|x − x₀| ≤ a, |y − y₀| ≤ b}, то задача Коши y' = f(x, y), y(x₀) = y₀ имеет хотя бы одно решение на |x − x₀| ≤ h = min(a, b/M).

Ключевое отличие от теоремы Пикара: утверждается только существование, но не единственность. Без условия Липшица через одну точку может проходить несколько решений.

Доказательство использует теорему Арцела–Асколи: из последовательности приближений выбирается равностепенно непрерывная подпоследовательность, которая сходится (компактность в пространстве непрерывных функций).

Примеры неединственности

Пример 1: y' = 2√|y|, y(0) = 0.

Правая часть f = 2√|y| непрерывна, но ∂f/∂y = 1/√|y| → ∞ при y → 0 — условие Липшица нарушено.

Решения задачи Коши с y(0) = 0:

• y ≡ 0 (нулевое решение),
• y = (x − a)² при x > a и y = 0 при x ≤ a для любого a ≥ 0 (бесконечно много решений!).

Таким образом, через начало координат проходит целое однопараметрическое семейство различных решений. Начальное условие не определяет будущую траекторию.

Пример 2 (взрыв за конечное время): y' = y², y(0) = 1.

Решение: y(t) = 1/(1 − t). При t → 1⁻ решение стремится к +∞. Это «взрыв за конечное время» — типичное поведение квадратично нелинейных уравнений.

Физический смысл: Уравнение вида ẏ = y² моделирует самоусиливающийся процесс. Например, если скорость распространения пожара пропорциональна квадрату площади пожара, он «взорвётся» за конечное время. Это важно при моделировании цепных реакций, лавинных процессов, плазмы в термоядерных реакторах.

Теорема сравнения Чаплыгина

Мощный инструмент для получения оценок решений без явной формулы.

Формулировка: Пусть y' ≤ f(x, y) и z' = f(x, z), y(x₀) ≤ z(x₀). Тогда y(x) ≤ z(x) для всех x ≥ x₀ (при условии, что f монотонно возрастает по y).

Смысл: «Если начать ниже, то и останешься ниже». Если мы умеем решать «оценочное» уравнение z' = f(x, z), то можем ограничить сверху решение неравенства y' ≤ f(x, y).

Применение: Для квадратичного уравнения y' = y² + 1 оценим решение сверху. Поскольку y² + 1 ≤ y² + y² + 1 = 2y² + 1 при y ≥ 1, из уравнения z' = z² (начиная с тех же начальных данных z₀ = y₀ > 1) получаем, что y ≤ z = 1/(C − x). Это даёт оценку момента взрыва сверху.

Лемма Гронуолла–Беллмана

Одна из самых используемых лемм в теории ДУ, позволяющая получать экспоненциальные оценки.

Формулировка: Если u(x) ≥ 0 и u(x) ≤ α + β ∫_{x₀}^x u(t) dt для всех x ≥ x₀, то u(x) ≤ α · e^(β(x − x₀)).

Доказательство: Обозначим V(x) = α + β ∫_{x₀}^x u(t) dt. Тогда V' = βu ≤ βV и V(x₀) = α. Из линейного сравнительного уравнения V' = βV получаем V ≤ α e^(β(x−x₀)), откуда u ≤ V ≤ α e^(β(x−x₀)).

Стандартное применение: Оценка зависимости от начальных данных. Если δ = |y(x) − z(x)| (разность двух решений) удовлетворяет δ' ≤ L·δ (из условия Липшица), то δ(x) ≤ δ(x₀) · e^(L|x−x₀|). Вся теорема Пикара о непрерывной зависимости от начальных данных — следствие леммы Гронуолла.

Численные методы и ошибки аппроксимации

Теорема Пеано и лемма Гронуолла — теоретическая основа анализа численных методов. Метод Эйлера: yₙ₊₁ = yₙ + h·f(xₙ, yₙ). Погрешность на каждом шаге O(h²), а за N = T/h шагов накапливается как O(h) — глобальная ошибка первого порядка. Лемма Гронуолла позволяет оценить рост накопленной погрешности через те же экспоненциальные множители.

Жёсткие системы ОДУ и неявные методы

Жёсткое уравнение (stiff ODE) — это система, в которой одновременно присутствуют быстро и медленно изменяющиеся компоненты решения: λ₁ ≪ λ₂ (оба отрицательные). Метод Эйлера и явные методы Рунге–Кутты требуют шага h < 2/|λ₂|, что катастрофически мало. Неявный метод Эйлера: yₙ₊₁ = yₙ + h·f(xₙ₊₁, yₙ₊₁) — условно безусловно устойчив при любом h. Цена — решение нелинейного уравнения на каждом шаге (метод Ньютона). Метод трапеций Кранка–Николсон — золотой стандарт для жёстких задач. В химической кинетике жёсткость возникает всегда: скорости быстрых реакций на 6–8 порядков выше медленных. Без специальных методов (LSODE, VODE) численное интегрирование таких систем невозможно.

Вопрос для размышления: Лемма Гронуолла показывает, что погрешность растёт не быстрее экспоненты. Означает ли это, что численное интегрирование всегда «работает» на достаточно малых шагах? Что происходит при численном решении уравнений с λ > 0 (хаотических систем)?