Статьи Соросовского Образовательного журнала в текстовом формате
Излагается один из методов интегрирования обыкновенных дифференциальных уравнений - метод сравнения. На его основе решается задача об устойчивости решений. Изложение элементарно и доступно.
ИНТЕГРИРОВАНИЕ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ МЕТОДОМ СРАВНЕНИЯЕ. В. ВОСКРЕСЕНСКИЙ
Мордовский государственный университет им. Н.П. Огарева, Саранск
ВВЕДЕНИЕ
В статье рассказано о некоторых методах решения дифференциальных уравнений вида
где функция f определена и непрерывна в некоторой области D плоскости R 2. Пусть решение уравнения (1), удовлетворяющее начальным данным (t0 , x0) единственно, и его будем обозначать так: x(t: t0 , x0). Итальянский математик Дж. Пеано в 1890 году доказал, что при любых t0 и x0 , таких, что точка (t0 , x0) k k D, всегда существует решение x(t: t0 , x0) уравнения (1). Как найти его? Оказалось, найти решение в виде аналитических формул, как это удается сделать для уравнений с разделяющимися переменными и некоторых других, в общем случае невозможно. Ж. Лиувилль в 1841 году впервые указал уравнение, которое не интегрируется в квадратурах. Более того, большинство уравнений в этом смысле решению принципиально не поддается. Так, решение уравнения
нельзя выразить через элементарные функции и интегралы от них. Другими словами, уравнение (2) является нерешаемым в классе элементарных функций. И это несмотря на то, что на основании теоремы Пеано решения этого уравнения существуют. Вскоре обнаружилось, что этот факт не является трагическим. Оказалось, что более важное значение имеет знание свойств решений, чем его задание формулой, а свойства можно узнать без использования формулы, если даже она имеется. Поэтому в современной математике термин "интегрирование" дифференциального уравнения обращен к изучению свойств решений и в зависимости от методов изучения свойств имеются различные методы интегрирования: численные, качественные, асимптотические и т.д.
Важное значение имеет знание оценок вида
| x(t: t0 , x0) | # y(t),
где y - известная функция. Если она определена на полуоси [T, + ?), то из неравенства (3) можно получить важные асимптотические свойства решений. Основным методом получения оценок вида (3) в настоящее время является метод сравнения, суть которого заключается в следующем.
Рассматриваются два уравнения вида (1)
причем свойства решений уравнения (5) известны и функция
R (t, x) = | f1(t, x) - f2(t, x) |
в каком-нибудь смысле является малой в области D. Например:
1) R (t, x) # m, где m - достаточно малое положительное число, или
2) сходится и т.д.
Малость в каждой конкретной задаче индивидуальна: в некоторых случаях она вида 1), в других - вида 2). Могут быть и другие виды малости. Оказалось, что при подходящей малости R многие важные свойства решений уравнения (4) наследуются от свойств решений уравнения (5). В этом суть метода сравнения. Уравнение (5) называется уравнением сравнения, а уравнение (4) - исследуемым уравнением, сам процесс определения характера малости R и нахождения списка наследуемых свойств решениями называется интегрированием дифференциальных уравнений методом сравнения.
НЕРАВЕНСТВО ВАЖЕВСКОГО
Уравнения сравнения подбираются так, чтобы они не были слишком сложными. Малость R часто задается при помощи скалярной функции, которая определяется на множестве, когда переменная t меняется в пределах полуоси, а вторая переменная z $ $ 0. Пусть уравнение (5) имеет простейший вид
и
| f (t, x) | = R (t, x) # l(t, | x |)
при всех x и t > T, а l - непрерывная неотрицательная функция, такая, что произвольные начальные данные (t0 , z0), t0 $ 0, z0 $ 0 определяют единственное решение уравнения
Тогда
| f (t, x) | # l(t, | x |)
и малость R определяется свойствами функции l.
В основе метода сравнения находится простая теорема [1], доказанная Г. Важевским в 1950 году.
Теорема (Важевский). Пусть функция V такая, что:
1) она определена и непрерывна на множестве I0 = = [t0 , b], V(t0) # z0 ;
2) на множестве I0 .
Тогда если z(t : t0 , z0) определено на I0 , то V(t) # # z(t : t0 , z0) при всех t из I0 .
В дальнейшем будем считать, что D состоит из точек плоскости R 2, когда t > T, а x - произвольное число и любые начальные данные однозначно определяют решение x(t: t0 , x0). С учетом неравенства (7) для решения x(t: t0 , x0) уравнения (1) получим
Тогда на основании теоремы Важевского имеем
| x(t: t0 , x0) | # z(t: t0 , z0)
для всех | x0 | # z0 и t k I0 .
Пример 1. Рассмотрим уравнение
Тогда
и уравнение (6) здесь имеет вид
Поэтому
Заметим, что при любых z0 и любом t0 > 1 решение уравнения (9) принимает малые значения при всех t $ t0 . Это свойство решения x = 0 уравнения (9) называется устойчивостью решения x = 0.
УСТОЙЧИВОСТЬ РЕШЕНИЙ
Перейдем к точной формулировке понятия устойчивости решения. Пусть x = 0 есть решение уравнения (1).
Определение. Решение x = 0 называется устойчивым, если для любой окрестности U точки x = 0 существует такая меньшая окрестность U0 K U, что любое решение, проходящее в момент времени t0 через точку x0 k U0 , остается в U при всех t $ t0 .
Понятие устойчивости отнесено здесь лишь к решению x = 0. Однако оно легко переносится и на другие решения. Пусть x(t: t0 , x0) - исследуемое решение. Выполним для уравнения (1) замену
y = x - x(t: t0 , x0),
тогда уравнение (1) перейдет в новое уравнение
для которого y = 0 является решением. Если это решение устойчиво, то решение x(t: t0 , x0) называется устойчивым, если y = 0 неустойчиво, то x(t: t0 , x0) - неустойчивым. Поэтому впредь считаем, что замена (10) выполнена, и поэтому всегда задача об устойчивости решения x(t: t0 , x0) сводится к изучению устойчивости решения y = 0 уравнения (11), а это уравнение такого же типа, что и уравнение (1). По этой причине исследование всегда ведется лишь относительно решения y = 0 уравнения (1). Не будем отступать здесь от традиций.
Неравенство (8) и пример 1 показывают, что свойство устойчивости решения x = 0 уравнения сравнения наследуется уравнением (1) для решения x = 0. Таким образом, эта задача для уравнения (9) получила очень простое решение, тогда как оно не интегрируется в квадратурах.
ЛЯПУНОВСКИЕ ПРЕОБРАЗОВАНИЯ
Прежде чем решать задачу об устойчивости решения x = 0 уравнения (1), полезно уравнение (1) преобразовать так, чтобы полученное уравнение было простым. Если преобразование имеет вид
x = j(t, y)
и | x | # k1| y |, | y | # k2| x |, где k1 , k2 > 0, то оно переводит уравнение (1) в уравнение, для которого решение y = 0 имеет тот же характер, что и решение x = 0 уравнения (1). Другими словами, эти решения либо одновременно устойчивы, либо одновременно неустойчивы. Такие преобразования называются ляпуновскими [2].
Предположим, что f непрерывно дифференцируема по x и
| f (t, x) | # y(t)| x |,
где y - непрерывная функция, сходится.
Пусть . Тогда l(t, | x |) = y(t)| x | и на основании теоремы Важевского имеем
| x(t: t0 , x0) | # k1| x0 |, k1 = ek, t $ t0 .
Рассмотрим замену
Тогда | x | # c1 | y |, c1 = kk1 . Легко проверяется неравенство | y | # c2 | x |, c2 > 0. Следовательно, преобразование (13) ляпуновское и переводит уравнение (1) в уравнение
Уравнение (14) является простейшим, и все его решения являются устойчивыми. Заметим, что и уравнение (9) ляпуновским преобразованием переводится в уравнение (14).
ВОЗМУЩЕННЫЕ УРАВНЕНИЯ
Теперь рассмотрим уравнения более общего вида
которые получаются из уравнения (1) путем добавления к правой части непрерывной функции f0 . Такие уравнения называются возмущенными.
Предположим, что для уравнения (1) существует ляпуновское преобразование j, переводящее уравнение (1) в уравнение (14). Тогда этим же преобразованием j уравнение (15) преобразуется в уравнение
Еще предположим
а функция l обладает такими же свойствами, что и l из уравнения (6). Тогда если решение уравнения (6) z = 0 устойчиво, то и решение уравнения (15) x = 0 также устойчиво.
Пример 2. Исследовать характер решения x = 0 уравнения
где функция f0 удовлетворяет всем перечисленным выше условиям. Так как ляпуновским преобразованием уравнение (9) переводится в уравнение (14), то характеры решений x = 0 уравнения (18) и y = 0 уравнения (16) совпадают. Из неравенства (17) следует устойчивость решения x = 0 уравнения (18), если решение z = 0 уравнения (6) устойчиво.
ЗАКЛЮЧЕНИЕ
В статье сделана попытка изложить один из важнейших современных методов интегрирования дифференциальных уравнений [3, 4]. При помощи этого метода решены важнейшие прикладные задачи, относящиеся к исследованию космоса, небесной механики, экологии и многих других областей знаний. Более подробное изложение требует привлечения серьезного аппарата. Однако изначальные идеи метода требуют лишь некоторого минимального запаса знаний из математики. Его дальнейшее уменьшение не позволило бы изложить суть дела.
ЛИТЕРАТУРА
1. Руш Н., Абетс П., Лалуа М. Прямой метод Ляпунова в теории устойчивости. М.: Мир, 1980. 300 с.
2. Воскресенский Е.В. Ляпуновские группы преобразований // Изв. вузов. Математика. 1994. ╧ 7. С. 1-7.
3. Матросов В.М. Принцип сравнения с вектор-функцией Ляпунова // Дифференц. уравнения. 1968. Т. 4, ╧ 8. С. 1374-1386.
4. Воскресенский Е.В. Методы сравнения в нелинейном анализе. Изд-во Сарат. ун-та, 1990. 224 с.
* * *
Евгений Викторович Воскресенский, доктор физико-математических наук, профессор, зав. кафедрой прикладной математики Мордовского государственного университета, заслуженный деятель науки Республики Мордовия, президент Средневолжского математического общества. Область научных интересов - обыкновенные дифференциальные уравнения, математическая теория управления, численные методы. Автор более 150 работ, в том числе трех монографий.