равнения существует, когда оно единственно, решается так называемыми теоремами существования и единственности. Эти теоремы очень важны, как для самой теории, так и для практики. Они гарантируют законность применения качественных методов теории дифференциальных уравнений для решения задач естествознания и техники. Численному интегрированию дифференциального уравнения обязательно должно предшествовать обращение к теоремам существования и единственности. И это необходимо делать для того, чтобы избежать недоразумений или вообще неправильных выводов.
Теорема существования. Если в уравнении
![](http://www.rae.ru/use/i/2012/4/131.jpg)
, y(x0) = y0, (2)
определенное на некотором интервале, содержащем точку x0.
Теорема существования и единственности. Если в уравнении функция f определена и непрерывна в некоторой ограниченной области D плоскости (x,y), причем она удовлетворяет в области D условию Липшица по переменной y, т.е.
Теорема о продолжении. При выполнении условий теоремы существования или теоремы существования и единственности всякое решение уравнения с начальными данными (x0, y0) Î D может быть продолжено до точки, сколь угодно близкой к границе области D. При этом в первом случае продолжение, вообще говоря, будет не обязательно единственным, во втором же случае оно единственно.
Для иллюстрации «недоразумений» возникающих при использовании численных решений дифференциальных уравнений без учета теорем существования рассмотрим несколько примеров.
Пример 1. Требуется, используя численный метод интегрирования Эйлера с итерационной схемой
c шагом h =0,1, решить начальную задачу
y(-1) = 0,21 (3)
на отрезке [-1, 3].
Решение. (c помощью пакета Mathcad)
Обратимся теперь к теореме существования. Для исследуемой начальной задачи (3) функция f, определяемая равенством, определена и непрерывна во всей плоскости (x, y) за исключением точек оси абсцисс.
Таким образом, в соответствии с теоремой существования существует решение y(x) начальной задачи (3), определенное на некотором интервале, содержащем точку x0 = -1, и это решение по теореме о продолжении может быть продолжено до значения, близкого к значению y(x) = 0. В результате численного интегрирования получаем решение начальной задачи (3) на некотором интервале (a, b), где a < -1; 1,3 < b < 1,32. Однако, учитывая, что это уравнение с разделяющимися переменными, можно аналитически найти частное решение, удовлетворяющее начальной задаче (3)
Интегрируя, получаем, что
Отсюда следует, что решение начальной задачи (3) существует только для
Оказывается, обращение к теореме существования (и к теореме о продолжении) позволило отсечь отрезок (приблизительно [1,315; 5]), на котором решение исходной начальной задачи (3) заведомо не существует. Одно же только численное интегрирование приводит к ошибочному результату. Дело здесь в том, что при приближении решения y = y(x) к оси Ox угол наклона кривой приближается к 90°. Поэтому пока аргумент x изменяется на величину 0,1 значение y успевает перескочить ось Ox, и мы попадаем на интегральную кривую, отличную от исходной. А это происходит потому, что метод Эйлера учитывает угол наклона только в текущей точке.
Пример 2. Используя метод Эйлера, а затем метод Эйлера-Коши, с шагом h = 0,1 и итерационной схемой
,
где
,
решить начальную задачу
y(-1) = -1, (4)
на промежутке [-1, 1].
Решение. На базе Mathcad методом Эйлера, а затем методом Эйлера-Коши будем иметь:
Рис. 2
Получили чертеж (рис. 3) отличный от чертежа, изображенного на (рис. 2). Чтобы лучше разобраться в причине расхождения в результатах, проинтегрируем исходную начальную задачу. Разделяя переменные, имеем
или, окончательно, .
Становится понятно, что решение по методу Эйлера приближает функцию y1(x) = x3, а по улучшенному методу Эйлера - функцию
При этом как, y1 так и y2 являются решениями начальной задачи (4), а значит, для рассматриваемой на промежутке [-1; 1] начальной задачи имеет место неединственность.
Обращаясь теперь к теореме существования и единственности, отметим, что, так как функция f, заданная равенством, непрерывна во всей плоскости (x, y), то из теоремы существования следует, что существует решение начальной задачи (4), определенное на некотором промежутке, содержащем точку x0 = -1, и это решение по теореме о продолжении может быть продолжено на любой промежуток. Далее, поскольку
,
то функция
удовлетворяет условию Липшица по переменной y в любой области, не содержащей точки оси Ox. Если же область содержит точку оси Ox, то нетрудно показать, что в ней указанная выше функция условию Липшица не удовлетворяет. Поэтому из теоремы существования и единственности (и теоремы о продолжении) следует, что в данном случае решение начальной задачи может быть продолжено единственным образом, по крайней мере, до оси Ox. Но поскольку прямая y = 0 является особой интегральной прямой для дифференциального уравнения
,
то, как только y станет равным нулю решение начальной задачи (4) не может быть единственным образом продолжено за точку O(0, 0).
Рис. 3
Итак, обращение в данном случае к теореме существования и единственности (и теореме о продолжении) позволило разобраться в результатах численного интегрирования. Если речь идет о единственном на промежутке [-1; 1] решении начальной задачи (4), то оно существует и определено лишь на отрезке [-1; 0]. В общем же случае таких решений несколько.
Список литературы
1. Roberts C.E. Jr. Why teach existence and uniqueness theorems in the first course in ordinary differential equations? // Int. J. Math. Educ. Sci. Technol. - 1976. - Vol. 7, № 1. - P. 41-44.