Теорема Пикара–Линделёфа

Фундаментальный вопрос

Прежде чем решать дифференциальное уравнение, нужно убедиться, что решение вообще существует. И если оно существует — единственно ли оно? Эти вопросы не академические: от их ответа зависит, можно ли в принципе предсказывать поведение системы.

Рассмотрим физическую аналогию. Вы описываете траекторию частицы уравнением x'' = F(x)/m. Если в данной точке пространства скорость x' задана — единственна ли дальнейшая траектория? Если да, то «детерминизм Лапласа» справедлив: зная состояние системы в момент t₀, можно однозначно восстановить её прошлое и предсказать будущее. Если нет — система принципиально непредсказуема, даже без квантовой неопределённости.

Условие Липшица

Обычной непрерывности f(x, y) оказывается недостаточно для единственности. Нужно более сильное условие.

Функция f(x, y) удовлетворяет условию Липшица по y в области D с константой L, если для всех (x, y₁), (x, y₂) ∈ D:

|f(x, y₁) − f(x, y₂)| ≤ L|y₁ − y₂|.

Смысл: насколько бы близко ни взять два значения y, разность f «управляется» линейной функцией расстояния между ними. Это запрещает функции «слишком быстро расходиться» при малых изменениях y.

Как проверить: если ∂f/∂y ограничена в D (|∂f/∂y| ≤ L), то условие Липшица выполнено — это следует из теоремы Лагранжа о среднем значении.

Почему y' = y^(1/3) нарушает условие: ∂f/∂y = (1/3) y^(-2/3) → ∞ при y → 0. Действительно, через точку (0, 0) проходит бесконечно много решений: y = 0 и y = ±(2x/3)^(3/2).

Теорема Пикара–Линделёфа

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

Ограничение h = min(a, b/M) имеет простой геометрический смысл: решение не должно «вылететь» за рамку R раньше, чем мы достигнем конца отрезка по x.

Метод последовательных приближений Пикара

Конструктивное доказательство даёт алгоритм нахождения решения. Строим последовательность функций:

y₀(x) = y₀ (константа — начальное значение),

yₙ(x) = y₀ + ∫_{x₀}^x f(t, yₙ₋₁(t)) dt.

Это итерационный процесс: подставляем текущее приближение yₙ₋₁ в правую часть уравнения, интегрируем — получаем следующее приближение yₙ. Каждая итерация «учитывает» динамику системы точнее предыдущей.

Развёрнутый пример: y' = y, y(0) = 1.

y₀(x) = 1.

y₁(x) = 1 + ∫₀ˣ y₀(t) dt = 1 + ∫₀ˣ 1 dt = 1 + x.

y₂(x) = 1 + ∫₀ˣ y₁(t) dt = 1 + ∫₀ˣ (1 + t) dt = 1 + x + x²/2.

y₃(x) = 1 + ∫₀ˣ y₂(t) dt = 1 + x + x²/2 + x³/6.

Видно, что yₙ(x) = Σₖ₌₀ⁿ xᵏ/k! — частичная сумма ряда Тейлора для eˣ. В пределе: y(x) = eˣ — точное решение. ✓

Оценка скорости сходимости: ‖yₙ − y‖ ≤ ML^n h^(n+1)/(n+1)! → 0 — сходимость аналогична геометрической прогрессии, поэтому метод быстро сходится.

Теорема о максимальном интервале существования

Решение задачи Коши существует на некотором максимальном интервале (α, β). При x → β⁻ (или x → α⁺) одно из двух: либо x достигает границы области определения f, либо |y(x)| → ∞ (взрыв за конечное время).

Пример взрыва: y' = y², y(0) = 1. Решение: y = 1/(1 − x). При x → 1⁻ решение стремится к +∞. Максимальный интервал: (−∞, 1). За пределами x = 1 решение не существует.

Физически это означает: квадратичная нелинейность может приводить к «катастрофическому» нарастанию за конечное время. Это принципиально отличается от линейных уравнений, у которых решение всегда глобальное.

Пример глобального существования: Линейное уравнение y' = p(x)y + q(x) с непрерывными p, q — решение существует на всей оси.

Применение: предсказуемость физических систем

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

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

Теорема Пикара и размер области существования

Теорема гарантирует единственное решение на отрезке |x − x₀| ≤ h, где h = min(a, b/M), a — ширина прямоугольника, b — его высота, M = max|f|. Это локальная гарантия. Глобальное существование: если |f(x, y)| ≤ A + B|y| (линейный рост по y), решение существует на всём промежутке — следствие теоремы Гронуолла. Если |f| растёт быстрее линейного, решение может уйти на бесконечность за конечное время («взрыв»). Пример: y' = y², y(0) = 1 имеет решение y = 1/(1 − x) — «взрыв» при x = 1.