6. Решение задачи коши для обыкновенных дифференциальных уравнений 1-го порядка Справочная информация

Обыкновенное дифференциальное уравнение 1-го порядка

с

Рис.1.

вязывает независимую переменнуюx, искомую функцию y и её производную. Решение дифференциального уравнения заключается в отыскании функции y = y(x), обращающей это уравнение в тождество на конечном или бесконечном интервале (a, b). Различают общее и частные решения дифференциального уравнения. Общее решение уравнения имеет вид y = y(x, C), где C – произвольная постоянная интегрирования. Его графическим отображением является семейство кривых (см. рис.1), называемых интегральными. Каждая интегральная кривая является отображением частного решения, которое определяется своим значением постоянной C. Для выделения частного решения из множества кривых общего решения необходимо задать начальное условие

.

Такая постановка задачи отыскания решения дифференциальных уравнений называется задачей Коши (A.L.Cauchy, 1789–1857). Для существования единственного решения задачи Коши необходимо и достаточно существование и ограниченность правой части дифференциального уравнения f(x, y) и её частной производной f(x, y)/y в некоторой окрестности начальной точки (x₀, y₀).

Для численного решения задачи Коши существует множество методов, которые условно делятся на две группы: одношаговые и многошаговые. Все эти методы позволяют получить искомое решение дифференциального уравнения в виде таблично заданной функции, в той или иной мере согласую­щееся с истинным частным решением (см. рис.2). Эти группы методов различаются объёмом информации, которая используется для вычисления координат очередной точки табличной функции. Одношаговые методы используют значения функции и её производной только в одной предыдущей точке, в то время как многошаговые – в нескольких. К одношаговым методам решения задачи Коши относятся м

Рис.2.

етод Эйлера, модифицированный метод Эйлера, методы Рунге–Кутта и другие.

Метод Эйлера (L.Euler, 1768)

Он является старейшим методом решения задачи Коши и заключается в последовательном применении следующих формул

,

, ,

г

Рис.3.

еометрическая интерпре­тация которых приk = 0 представлена на рис.3. В точке (x₀, y₀) вычисляется значение производной dy/dx через правую часть f(x, y) исходного уравнения. Оно определяет тангенс угла наклона касательной к графику точного решения задачи Коши. Следующая точка численного решения определяется как точка на этой касательной с абсциссой x₁ = x₀ + h. В компактном виде для k = 0, 1, 2,… эти соотношения записываются следующим образом

, .

Метод Эйлера относится к методам первого порядка точности, поскольку его решение совпадает с истинным только в том случае, когда последнее является линейной функцией y = a₁+ a₂x. Его абсолютная погрешность ε_абс(x_k+1, h) на каждом шаге пропорциональна величине h². Это обусловлено тем, что в качестве направления, определяющего положение следующей точки численного решения, используется касательная в левой точке каждого отрезка [x_k, x_k+1]. На рис.3 видно, что для получения более точного численного решения недостаточно знания параметров функции в единственной левой точке отрезка [x_k, x_k+1]. Требуется собрать дополнительную информацию о её поведении на отрезке интегрирования для отыскания решения при x = x_k+1 с меньшей погрешностью. Для этого можно использовать некоторые промежуточные направления, определяемые касательными к графику неизвестного точного решения в характерных точках рассматриваемого отрезка (крайние точки, середина отрезка и т.д.).