Метод Пикара¶

Метод Пикара — итерационный метод приближённого решения обыкновенных дифференциальных уравнений, основанный на последовательных приближениях и интегральном представлении решения. Метод назван в честь французского математика Шарля Эмиля Пикара, который в конце XIX века использовал его для доказательства существования и единственности решения задачи Коши. Помимо теоретического значения, метод применяется для численного построения решения в виде ряда или последовательности функций.

¶Суть метода

Рассматривается задача Коши для обыкновенного дифференциального уравнения первого порядка:

\[ y' = f(x, y), \quad y(x_0) = y_0. \]

Интегрируя обе части уравнения от \(x_0\) до \(x\), получают эквивалентное интегральное уравнение:

\[ y(x) = y_0 + \int_{x_0}^{x} f(t, y(t)) \, dt. \]

Метод Пикара состоит в построении последовательности функций \(\{y_n(x)\}\), где нулевое приближение выбирается как начальное условие:

\[ y_0(x) = y_0, \]

а каждое следующее приближение вычисляется по рекуррентной формуле:

\[ y_{n+1}(x) = y_0 + \int_{x_0}^{x} f(t, y_n(t)) \, dt. \]

Если правая часть \(f(x, y)\) непрерывна и удовлетворяет условию Липшица по переменной \(y\) в некоторой окрестности точки \((x_0, y_0)\), то последовательность \(\{y_n(x)\}\) равномерно сходится к точному решению на некотором интервале \(|x - x_0| \le h\). Этот факт составляет основу классической теоремы Пикара — Линделёфа о существовании и единственности решения.

¶Сходимость и оценка погрешности

Скорость сходимости метода характеризуется оценкой:

\[ |y_n(x) - y(x)| \le \frac{M L^n h^{n+1}}{(n+1)!}, \]

где \(M\) — верхняя граница модуля функции \(f\), \(L\) — константа Липшица, \(h\) — радиус интервала. Благодаря факториалу в знаменателе погрешность быстро убывает с ростом номера итерации, однако на практике вычисление интегралов быстро усложняется, и метод применяют лишь для первых нескольких приближений.

¶Практическое применение

Метод Пикара используется в следующих случаях:

• для получения приближённого аналитического выражения решения в виде многочлена или ряда (например, при решении уравнения \(y' = x^2 + y^2\));
• для доказательства теоремы существования и единственности в курсах дифференциальных уравнений;
• в учебных целях для демонстрации принципа сжимающих отображений.

Для численного интегрирования дифференциальных уравнений на компьютере метод Пикара применяется редко из-за трудоёмкости вычисления интегралов и необходимости символьных преобразований. В вычислительной практике предпочтение отдаётся методам Рунге — Кутты, Адамса и др.

¶Пример

Рассмотрим задачу Коши:

\[ y' = y, \quad y(0) = 1. \]

Точное решение — \(y = e^x\). Последовательные приближения по методу Пикара:

\[ y_0 = 1, \]

\[ y_1 = 1 + \int_0^x 1 \, dt = 1 + x, \]

\[ y_2 = 1 + \int_0^x (1 + t) \, dt = 1 + x + \frac{x^2}{2}, \]

\[ y_3 = 1 + x + \frac{x^2}{2} + \frac{x^3}{6}. \]

Таким образом, \(n\)-е приближение совпадает с частичной суммой ряда Тейлора для \(e^x\), и в пределе получается точное решение.

¶Ограничения

Метод Пикара требует, чтобы функция \(f(x, y)\) была достаточно гладкой. При наличии разрывов или неединственности решения метод может не сходиться или сходиться к неверному результату. Кроме того, сходимость гарантирована лишь в локальной окрестности начальной точки; за пределами интервала сходимости последовательность приближений может расходиться.

¶Источники

• Эльсгольц Л. Э. «Дифференциальные уравнения и вариационное исчисление».
• Тихонов А. Н., Васильева А. Б., Свешников А. Г. «Дифференциальные уравнения».
• Петровский И. Г. «Лекции по теории обыкновенных дифференциальных уравнений».