Статьи Соросовского Образовательного журнала в текстовом формате
В статье дается введение в теорию обыкновенных дифференциальных уравнений первого порядка на плоскости. Описываются следующие основные понятия и методы, используемые в теории таких уравнений: поле направлений, элементарные методы построения траекторий поля направлений и интегрирования дифференциальных уравнений, теорема существования и единственности решения задачи Коши, метод Эйлера.
ПОЛЯ НАПРАВЛЕНИЙ И СООТВЕТСТВУЮЩИЕ ИМ ТРАЕКТОРИИМ. И. ВИШИК
Московский государственный университет
им. М.В. Ломоносова
1. ПОЛЕ НАПРАВЛЕНИЙ И ТРАЕКТОРИИ ЭТОГО ПОЛЯ
Во многих задачах физики и механики исследуются установившиеся движения жидкости или газа. При этом представляют интерес траектории движения частиц. Обычно они представляют собой семейство кривых. Мы будем рассматривать для простоты лишь семейства кривых на плоскости R2, заполняющие некоторую область G ? R2. Направление установившегося движения жидкости вдоль кривой (изображающей движение частицы) в данной ее точке P определяется касательной к этой кривой в точке P. Для определенности мы будем это направление изображать отрезком касательной, центр которого есть точка P. Если в каждой точке семейства кривых провести такой отрезок касательной, мы получим семейство отрезков с отмеченными точками P (их серединами). Такое семейство отрезков с отмеченными точками P называется полем направлений этого семейства кривых. Можно также определить поле направлений независимо от семейства кривых.
Определение 1. Говорят, что в плоской области G задано поле направлений, если через каждую точку P k G проходит отрезок, имеющий эту точку P своей серединой (рис. 1).
Определение 2. Траекторией поля направлений (или интегральной кривой поля направлений) называется кривая L, L ? G, которая в каждой точке P k G касается отрезка поля направлений, имеющего эту точку P своей серединой (рис. 2).
В физических, геометрических и других вопросах естествознания возникает следующая задача: пусть в области G задано поле направлений (рис. 1); найти все траектории (интегральные кривые) L этого поля направлений.
Кроме того, представляет интерес еще такая задача: в области G задано поле направлений (рис. 1) и точка P0 k G ; найти траекторию этого поля направлений, которая проходит через точку P0 .
2. НЕКОТОРЫЕ ЭЛЕМЕНТАРНЫЕ МЕТОДЫ ПОСТРОЕНИЯ ТРАЕКТОРИЙ ПОЛЯ НАПРАВЛЕНИЙ
а) Рассмотрим поле направлений на плоскости с выколотой точкой (началом координат), определяемое следующим образом: отрезок поля направлений, проходящий через точку P = (x, y) ? (0, 0), перпендикулярен отрезку, соединяющему (x, y) с началом координат (0, 0). Легко видеть, что траекториями L этого поля направлений будут любые окружности с центром в начале координат x 2 + y 2 = R 2, где R > 0 - любое фиксированное число.
б) Рассмотрим поле направлений на плоскости с выколотым началом координат, состоящее из отрезков, лежащих на прямых, проходящих через начало координат. Очевидно, траекториями такого поля направлений являются любые лучи без начала, выходящие из начала координат.
Замечание. В этих примерах поля направлений не задаются в начале координат. Отметим, что эти поля направлений нельзя доопределить в точке (0, 0) таким образом, чтобы они оставались непрерывно зависящими от (x, y) в точке (0, 0).
в) Пусть дано поле направлений, заданное следующим свойством: угловой коэффициент tg a отрезка этого поля с центром в точке (x, y) определяется равенством
(на оси ординат поле вертикально; в начале координат - опять-таки не определено).
Легко видеть, что траектории этого поля - это ветви гипербол вида
xy = C,
где C - любая константа (рис. 3).
Действительно, пусть y = y(x) - одна из гипербол (2). Тогда
xy(x) ╞ C.
Взяв производную по x от обеих частей (3), получим
Отсюда
Так как - угловой коэффициент касательной к графику (x, y(x)) гиперболы, то из (4) следует, что
Таким образом, доказано, что любая гипербола, определяемая уравнением (2), является траекторией поля направлений, заданного условием (1). Отметим, что положительные и отрицательные полуоси прямых Ox и Oy также являются траекториями поля направлений.
3. НЕКОТОРЫЕ ЭЛЕМЕНТАРНЫЕ МЕТОДЫ ИНТЕГРИРОВАНИЯ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ
Весьма часто вместо поля направлений задают соотношения вида (4), которые, как мы видим, соответствуют полям направлений. (Ниже мы это объясним подробнее.) В (4) - угловой коэффициент касательной к траектории поля направлений в точке (x, y(x)), а правая часть (4) указывает, чему равно в точке (x, y(x)). Соотношения вида (4) называются дифференциальными уравнениями. При этом вместо y(x) пишут обычно y. Таким образом, (4) записывают в виде
Задача о нахождении всех функций y = y(x), удовлетворяющих такому соотношению, называется задачей интегрирования этого дифференциального уравнения. Приведем некоторые примеры дифференциальных уравнений и методы их интегрирования.
а) Рассмотрим дифференциальное уравнение вида
где f(x) - заданная непрерывная функция. Найдем решения y = y(x) = j(x) этого уравнения, то есть такие функции y = j(x), для которых
Как известно, первообразной функции f(x) называется такая дифференцируемая функция F(x), что . Отметим, что если F(x) - первообразная функции f(x), то F(x) + C также является первообразной f(x).Таким образом, решениями уравнения (6) является семейство функций y = F(x) + C или - неопределенный интеграл функции f(x).
Отметим, что уравнению (6) соответствует поле направлений, полученное по следующему правилу: угловой коэффициент tg a отрезка поля направлений в точке P = (x, y) k G, где G - область в R2, равен f(x):
Таким образом, отрезки поля направлений, лежащие на одной вертикальной прямой, то есть имеющие в точках P = (x, y) одинаковые значения абсциссы x, параллельны (так как в них согласно (7) tg a = f(x) = const при заданном x). Интегральными кривыми этого поля направлений являются траектории, полученные любым сдвигом одной из его траекторий y = F(x) в направлении оси Oy (y = F(x) + C ).
б) Рассмотрим уравнение
где g(y) - непрерывная функция переменной y. Для нахождения решений этого уравнения удобно считать y независимой переменной, а x = y(y) - искомой функцией. При этом воспользуемся следующим известным фактом: если y = j(x) - решение уравнения (8), то обратная к y = j(x) функция x = y(y) удовлетворяет уравнению
Мы найдем решения уравнения (8) или, что эквивалентно, уравнения (9) в таких областях плоскости (x, y), в которых g(y) ? 0. Уравнение (9) имеет такой же вид, как и уравнение (6), если поменять x и y местами. Поэтому решения уравнения (9) задаются формулой
где G(y) - первообразная функции 1/ g(y). Что касается исходного уравнения (8), то его решениями являются функции y = F(x - C), где y = F(x) - функция, обратная к x = G(y).
Отметим, что поле направлений, соответствующее уравнению (8), состоит из отрезков, одинаково направленных при фиксированном значении y, так как согласно (8) .
При этом, зная одно решение x = G(y) уравнения (8) (или (9)), все остальные решения получаются с помощью параллельного сдвига этого решения в направлении оси Ox: x = G(y) + C.
4. ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ ПЕРВОГО ПОРЯДКА НА ПЛОСКОСТИ И СООТВЕТСТВУЮЩИЕ ИМ ПОЛЯ НАПРАВЛЕНИЙ
Как известно, производная функции y = j(x) имеет следующий геометрический смысл: пусть L - график функции y = j(x), то есть совокупность точек (x, j(x)), a < x < b, на плоскости (x, y). Тогда совпадает с угловым коэффициентом k = tg a касательной, проведенной к кривой L в точке (x, j(x)). При этом a, 0 # a # p, - угол между касательной к L в точке (x, j(x)) и положительной полуосью Ox. Очевидно, предполагается, что существует производная в точке x.
Ряд физических и геометрических задач приводит к нахождению всех таких функций y = y(x) = = j(x), производная которых в точке (x, y) = (x, y(x)) = (x, j(x)) равна значению заданной функции f(x, y) в точке (x, y(x)) = (x, j(x)). Иными словами, задается уравнение вида
Решением этого уравнения называется любая функция y = y(x) ╞ j(x), a < x < b, производная которой совпадает со значением правой части (10), в которой вместо y подставлена функция j(x), то есть
Например, пусть дано дифференциальное уравнение вида
Его решениями являются функции вида , где C - произвольная постоянная. Действительно, для этих функций
а правая часть уравнения (11) для них
Таким образом, для
Мы доказали, что семейство решений при любом постоянном значении C является решением уравнения (11).
Уравнению (10) можно сопоставить поле направлений следующим образом: через каждую точку P = (x, y) k G, в которой задано f(x, y), провести отрезок прямой, имеющей угловой коэффициент
k = tg a = f(x, y).
При этом точка (x, y) = P - середина такого отрезка. Тогда в области G возникает поле направлений, полученное согласно (12).
Задача о нахождении решений
y = j(x), a < x < b,
уравнения (10) и задача о нахождении траекторий L поля направлений, порожденного уравнением (10), эквивалентны. Действительно, график L любого решения y = j(x), a < x < b, уравнения (10) является траекторией поля направлений этого уравнения (10).
В самом деле, согласно (10) в любой точке (x, j(x)) k L отрезки построенного по уравнению (10) поля направлений (см. рис. 4) касаются кривой L в силу условия (12), которому удовлетворяет функция j(x):
Таким образом, график L = {(x, j(x)), a < x < b}, где j(x) - решение уравнения (10), является траекторией поля направлений, порожденного уравнением (10).
С другой стороны, каждому полю направлений (см. рис. 1) можно сопоставить дифференциальное уравнение вида (10), для которого оно является его полем направлений. Действительно, в каждой точке P области G с декартовыми координатами (x, y) (P = (x, y)) можно вычислить tg a = k - угловой коэффициент отрезка поля направлений, проходящего через эту точку. Положим в точке (x, y) функцию f(x, y) равной значению tg a отрезка поля направлений, проходящего через (x, y):
f(x, y) = tg a
(tg a = k - угловой коэффициент отрезка поля направлений в точке (x, y) (рис. 1)).Теперь заданному полю направлений (см. рис. 1) сопоставим дифференциальное уравнение
где f(x, y) определено формулой (13). Легко видеть, что траектория L поля направлений (рис. 1), являющаяся графиком некоторой дифференцируемой функции y = j(x) (L = {(x, j(x)), a < x < b}), удовлетворяет уравнению (14). Действительно,
где tg a - угловой коэффициент касательной к кривой L в точке (x, j(x)). По приведенному выше построению (см. (13)) tg a = f(x, y) = f(x, j(x)).Таким образом, функция y = j(x) удовлетворяет тождеству
и, следовательно, является решением уравнения (14). Мы доказали, что между полями направлений (рис. 1) и дифференциальными уравнениями вида (10) имеется взаимно однозначное соответствие, устанавливаемое соотношением (12).
5. НАЧАЛЬНАЯ ЗАДАЧА ДЛЯ ЗАДАННОГО ДИФФЕРЕНЦИАЛЬНОГО УРАВНЕНИЯ
Многие задачи геометрии, физики, механики и других естественных наук приводят к нахождению решения y = j(x) дифференциального уравнения (10), проходящего через точку (x0 , y0) области G, то есть такого решения y = j(x), a < x < b, a < x0 < b, для которого
y(x0) = y0 .
Эту задачу называют начальной задачей, или задачей Коши.
Предполагается, что:
(I) функция f(x, y) (то есть правая часть уравнения (10)) непрерывно зависит от x, y, (x, y) k G ;
(II) производная функции f(x, y) по второму аргументу (y) при фиксированном первом аргументе (x) непрерывно зависит от точки (x, y) k G.
Приведем основную теорему о разрешимости этой задачи.
Теорема 1. Предположим, что выполнены условия (I), (II).
1) Для всякой точки (x0 , y0) области G найдется решение y = y(x) = j(x) уравнения (10), удовлетворяющее условию (15).
2) Если два решения y = y(x) и y = j(x) уравнения (10) совпадают при x = x0 , то есть
y(x0) = j(x0),
то эти решения тождественно равны (y(x) = j(x)) для всех тех значений x, для которых они определены.
Идея доказательства утверждения 1) приводится в п. 6. Доказательство утверждения 2) можно найти в книгах [1 - 5].
Если выполнено утверждение 1), то говорят, что имеет место теорема существования решения начальной задачи (10), (15). При выполнении утверждения 2) говорят, что имеет место единственность решения этой задачи.
Как отмечалось в п. 4, каждому уравнению (10) отвечает поле направлений, состоящее из отрезков, имеющих середину в точке (x, y) и угловой коэффициент tg a (см. (14')). При этом если функция f(x, y) удовлетворяет условиям (I) и (II) настоящего пункта, то, согласно теореме 1, через каждую точку P = (x, y) k G проходит и притом единственная траектория этого поля.
6. МЕТОД ЛОМАНЫХ ЭЙЛЕРА ДЛЯ ПРИБЛИЖЕННОГО РЕШЕНИЯ НАЧАЛЬНОЙ ЗАДАЧИ (10), (15)
Нахождение явного решения y = j(x) начальной задачи (10), (15) иногда сопряжено с большими трудностями или просто невозможно. Поэтому применяют различные приближенные методы нахождения решения y = j(x) задачи (10), (15).
Опишем один из самых простых приближенных методов, а именно приближенный метод Эйлера.
Пусть дано уравнение (10) и начальная точка (x0 , y0) k G. Уравнению (10) соответствует поле направлений, отрезки которого имеют в точке (x, y) k G угловой коэффициент k, определяемый равенством (13):
k = tg a = f(x, y).
Выберем достаточно малое число h > 0. Ломаная Эйлера, отвечающая числу h, строится следующим образом (см. рис. 5): из начальной точки (x0 , y0) мы сдвигаемся вдоль отрезка поля направлений, имеющего середину в (x0 , y0), в точку (x0 + h, y1), где
y1= y0 + h tg a0 , tg a0 = f(x0 , y0).
Далее, из точки (x0 + h, y1) мы сдвигаемся вправо на шаг h вдоль отрезка поля направлений, имеющего середину в точке (x0 + h, y1). Так как в точке (x0 + h, y1) отрезок поля направлений имеет угловой коэффициент
k1 = tg a1 = f(x0 + h, y1),
где y1 задается формулой (18), то мы придем таким образом в точку (x0 + 2h, y2), где
y2 = y1 + h tg a1 = y1 + hf(x0 + h, y1),
и т.д. Мы получим ломаную y = jh(x), которая задается соотношениями
y0 + [x - x0]tg a0 , если x0 # x # x0 + h,
y1 + [x - (x0 + h)]tg a1 ,
jh(x) = если x0 + h # x # x0 + 2h,
y2 + [x - (x0 + 2h)]tg a2 ,
если x0 + 2h # x # x0 + 3h,
. . . . . . . . . . . . . . . . . . . . . . .
где
tg a0 = f(x0 , y0) ,
tg a1 = f(x0 + h, y1), y1 = y0 + hf(x0 , y0),
tg a2 = f(x0 + 2h, y2), y2 = y1 + hf(x0 + h, y1),
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
Полученная таким образом ломаная y = jh(x), x0 # x # x0 + hk называется ломаной Эйлера. Зафиксируем интервал (x0 , x0 + a) и выберем h и k такими, чтобы a = hk.
Теорема 2. Если функция f(x, y) в (10) удовлетворяет условиям (I) и (II) предыдущего пункта и число a выбрано настолько малым, что при достаточно малом h ломаные Эйлера y = jh(x) лежат в области G, то jh(x) при h 0 стремится к j(x), где y = j(x) - решение начальной задачи (10), (15).
Утверждение теоремы 2 интуитивно очевидно. Его строгое доказательство содержится в учебниках по обыкновенным дифференциальным уравнениям (см. [1 - 5]).
7. ПРИМЕР НЕЕДИНСТВЕННОСТИ РЕШЕНИЯ НАЧАЛЬНОЙ ЗАДАЧИ
Рассмотрим дифференциальное уравнение
Правая часть этого уравнения 3y2/3 ╞ f(x, y) ╞ f(y) - непрерывная функция переменных (x, y) на всей плоскости, следовательно, выполнено условие (I) пункта 5. Производная по y правой части (22)
обращается в ? при y = 0 и, следовательно, не является непрерывной функцией на оси x. Покажем, что если начальная точка (x0 , y0) лежит на оси x, то есть
y(x0) = 0 = y0 ,
то для начальной задачи (22), (23) нарушается свойство единственности ее решения.
Действительно, легко проверить, что решением задачи (22), (23) является кубическая парабола
y = (x - x0)3.
В самом деле, , а правая часть уравнения (22) равна 3[(x - x0)3]2/3. Следовательно, функция (24) удовлетворяет уравнению (22) и начальному условию (23): y(x0) = 0.
С другой стороны, прямая
y ╞ 0
также удовлетворяет уравнению (22) и начальному условию (23). Таким образом, через любую начальную точку (x0 , 0) проходят по крайней мере два решения y = (x - x0)3 и y = 0 уравнения (22). Мы видим, что вторая часть теоремы 1 (единственность) не имеет места для уравнения (22).
В том случае, когда правая часть уравнения (10) разрывна в точке (x0 , y0), может нарушаться и теорема существования, то есть первая часть теоремы 1.
ЛИТЕРАТУРА
1. Арнольд В.И. Обыкновенные дифференциальные уравнения. М., 1971.
2. Петровский И.Г. Лекции по теории обыкновенных дифференциальных уравнений. М., 1964.
3. Понтрягин Л.С. Обыкновенные дифференциальные уравнения. М., 1961.
4. Степанов В.В. Курс дифференциальных уравнений. М.; Л.: ГОНТИ, 1939.
5. Эльсгольц Л.Э. Дифференциальные уравнения. М., 1957.
* * *
Марко Иосифович Вишик, доктор физико-математических наук, профессор Московского государственного университета, главный научный сотрудник Института проблем передачи информации РАН, автор 232 научных работ и 4 монографий.