Портал | Содержание | О нас | Авторам | Новости | Первая десятка | Дискуссионный клуб | Чат Научный форум --> Первая десятка "Русского переплета" Темы дня:

Президенту Путину о создании Института Истории Русского Народа. |Нас посетило 40 млн. человек | Чем занимались русские 4000 лет назад?

| Кому давать гранты или сколько в России молодых ученых?

Описан универсальный метод исследования качественного поведения траекторий векторного поля на плоскости в окрестности как угодно сложной особой точки.

Н. Б. МЕДВЕДЕВА

Челябинский государственный университет

ВЕКТОРНОЕ ПОЛЕ

Векторное поле на плоскости с координатами (x, y) задается парой функций u1(x, y), u2(x, y). Это означает, что в точке плоскости с координатами (x, y) приложен вектор V(x, y) = (u1(x, y), u2(x, y)).

Траекторией векторного поля V называется кривая g(t) = (x(t), y(t)), t - вещественный параметр, такая, что функции x(t) и y(t) удовлетворяют системе обыкновенных дифференциальных уравнений первого порядка:

Траектория векторного поля в каждой точке (x, y) касается вектора поля V(x, y). Мгновенная скорость движения по траектории в точке (x, y) равна длине вектора V(x, y).

О векторных полях, физическом смысле связанных с ними понятий, а также о качественных методах исследования см., например, [1].

А. Пуанкаре принадлежит тезис о том, что для качественного исследования некоторого дифференциального объекта удобно привести его к более простому виду с помощью подходящей замены переменной.

Замена переменных

x = x(u, u), y = y(u, u)

в векторном поле V (или, что то же самое, в системе (1)) производится в соответствии с правилами дифференцирования сложной функции:

Отсюда

где

есть матрица Якоби замены переменных (2).

Рассмотрим задачу о построении локального фазового портрета векторного поля: требуется исследовать качественное поведение траекторий векторного поля вблизи фиксированной точки плоскости, причем нас будут интересовать только формы траекторий, но не скорость и не направление движения по ним.

Точка плоскости (x0 , y0) называется особой точкой векторного поля V или положением равновесия, если V(x0 , y0) = (0, 0). В противном случае точка (x0 , y0) называется неособой точкой векторного поля V. В некоторой окрестности неособой точки существует замена переменных, переводящая исходное векторное поле в векторное поле V0 = (1, 0). Траектории последнего векторного поля параллельны оси абсцисс. Поэтому и траектории исходного векторного поля вблизи неособой точки ведут себя подобно параллельным прямым. Поведение же траекторий в окрестности особой точки может быть весьма сложным.

Пример: линейные системы. Рассмотрим линейную систему на плоскости с невырожденной матрицей A. Начало координат является единственной особой точкой такой системы. В силу невырожденности системы оба собственных значения l1 и l2 матрицы A отличны от нуля. Поведение траекторий линейных невырожденных систем на плоскости изучают в университетском курсе обыкновенных дифференциальных уравнений [2]. Все возможные фазовые портреты в окрестности начала координат для такой системы изображены на рис. 1, а-г.

Рассмотрим гладкое векторное поле в окрестности особой точки на плоскости. Далее будем считать, что особая точка расположена в начале координат. Пусть A - матрица линейной части векторного поля в особой точке:

Тогда

где второе слагаемое стремится к нулю быстрее, чем вектор (x, y). Это означает, что векторное поле вблизи особой точки "почти" линейно.

Особая точка векторного поля называется:

- элементарной, если по крайней мере одно из собственных значений матрицы линейной части векторного поля в этой особой точке отлично от нуля;

- невырожденной, если оба собственных значения отличны от нуля. Всякая невырожденная особая точка является элементарной. В курсе обыкновенных дифференциальных уравнений доказывается, что в случае, когда вещественные части собственных значений невырожденной особой точки отличны от нуля, поведение траекторий в ее окрестности такое же, как у линейной части векторного поля, то есть возможны случаи, изображенные на рис. 1, а-в. В случае центра по линейным членам (то есть когда l1, 2 = ? ib - чисто мнимые числа) ситуация не столь проста. При одной и той же линейной части возможны и центр и фокус. Алгоритм различения центра и фокуса в случае центра по линейным членам восходит к А. Пуанкаре и А. Ляпунову и является классическим. Мы не будем на нем останавливаться, поскольку основной целью нашего исследования являются более сложные особые точки, не поддающиеся исследованию классическими методами.

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

1) вырожденный узел (пример: );

2) вырожденное седло (пример: );

3) седло-узел (пример: ).

Поведение траекторий в первых двух случаях качественно такое же, как в случае невырожденного узла и седла. Седло-узел изображен на рис. 1, д.

Элементарные особые точки являются как бы кирпичиками, на которые может быть "рассыпана" как угодно сложная особая точка. Точный смысл этого высказывания станет ясен ниже.

РАЗДУТИЕ ОСОБЕННОСТЕЙ

Для того чтобы исследовать окрестность неэлементарной особой точки, применяют так называемый метод раздутия особенностей. Существуют различные варианты этого метода. Суть его состоит в следующем. Особая точка заменяется некоторой кривой типа окружности с углами. Окрестность особой точки превращается в полуокрестность этой кривой с помощью одной или нескольких замен переменных специального вида. Векторные поля, полученные после замен переменных и определенные в различных кусках полуокрестности "ломаной окружности", имеют только элементарные особые точки. В окрестности каждой такой особой точки можно нарисовать фазовый портрет. Затем спроектировать эти картинки в окрестность исходной особой точки и получить фазовый портрет в ней (рис. 2).

Наиболее распространенными вариантами метода раздутия особенностей являются:

1) переход к полярным координатам (Ф. Дюмортье);

2) s-процесс [2];

3) метод диаграмм Ньютона [3, 4].

Первый вариант использует трансцендентные функции. Во втором картинки после раздутия расположены на двумерных неориентируемых многообразиях в трехмерном пространстве. Третий метод хорош тем, что использует степенные замены переменных, а после раздутия получаются плоские картинки, к тому же он является более быстрым, чем первые два.

Поясним метод диаграмм Ньютона на простом примере.

ДИАГРАММА НЬЮТОНА

Рассмотрим систему уравнений, зависящую от нескольких параметров:

Умножим первое уравнение системы на y, а второе - на x:

На плоскости переменных (i, j) отметим показатели всех мономов xiy j правой части: (2, 2), (1, 3), (0, 4), (6, 0) (рис. 3, а). Полученное множество, состоящее из четырех точек, называется носителем системы (3). Каждой точке носителя системы соответствует набор коэффициентов мономов правой части, который называется векторным коэффициентом. Например, точке (2, 2) соответствует векторный коэффициент (a, b). Ломаная, которая является огибающей полученного множества, называется многоугольником Ньютона системы (3), а та ее часть, которая "смотрит" на начало координат (в нашем случае состоящая из двух ребер), называется диаграммой Ньютона системы (3) (рис. 3, а).

Диаграмма Ньютона отражает сложность особой точки. Чем "меньше" диаграмма Ньютона, тем "проще" особая точка. Например, компоненты векторного поля, диаграмма Ньютона которого состоит из одной точки и изображена на рис. 3, б, имеют вид

V1(x, y) = x(l1 + f1(x, y)), V2(x, y) = y(l2 + f2(x, y)),

где f1(0, 0) = f2(0, 0) = 0. Поэтому (l1 , l2) - набор собственных значений линейной части этого поля. Так как (l1 , l2) ? (0, 0), то особая точка элементарная.

Как же "уменьшить" диаграмму Ньютона? Попытаемся это сделать с помощью степенных замен переменных.

СТЕПЕННЫЕ ЗАМЕНЫ ПЕРЕМЕННЫХ

Замены переменных вида

замечательны тем, что при таких заменах показатели степеней мономов в правой части системы преобразуются линейным образом и, значит, диаграмма Ньютона претерпевает линейные преобразования. Поэтому мы можем подобрать линейное преобразования так, чтобы новая особая точка имела "меньшую" диаграмму Ньютона, а значит, стала бы проще исходной.

Предложение 1. Если - матрица показателей замены (4), то в результате замены (4) точка носителя (i, j ) переходит в точку носителя новой системы с координатами где CT - транспонированная матрица, а (I, J ) - целочисленный вектор. Векторные коэффициенты точек носителя претерпевают линейное преобразование с матрицей C -1.

Пример. Степенные замены переменных применяют для интегрирования уравнений специального вида. Например, в однородной системе переменные разделяются с помощью замены y = ux. Система называется квазиоднородной с весами n и m переменных x и y, если ее носитель лежит на прямой ni + mj = = d, d > 0. Если m = n = 1, то квазиоднородная система является однородной. Замена переменных

x = zn, y = zmw

превращает квазиоднородную систему в систему, носитель которой - вертикальный отрезок. Такой системе соответствует уравнение с разделяющимися переменными.

Проверим утверждение предложения 1 для случая замены y = xu, которой соответствует матрица Действительно, рассмотрим мономиальную систему и сделаем в ней замену y = xu. После замены получаем

Применяя к обеим частям преобразование C -1, получаем

где

ЗАМЕНА ПЕРЕМЕННЫХ, СООТВЕТСТВУЮЩАЯ ВЕРШИНЕ ДИАГРАММЫ НЬЮТОНА

Рассмотрим вершину диаграммы Ньютона системы (3), имеющую координаты (2, 2). С помощью линейного преобразования и сдвига поместим ее в точку (1, 1) так, чтобы примыкающие к ней два ребра превратились в вертикальный и горизонтальный отрезки (см. рис. 3, б ). Например, линейное преобразование с матрицей подходит для этой цели. Соответствующая замена переменных x = = zw, y = zw 2, имеющая матрицу показателей C = C T, превращает систему (3) в систему, имеющую в начале координат (z, w) = (0, 0) элементарную особую точку. Соответствующие собственные значения можно сосчитать пользуясь предложением 1:

Предположим, что l1l2 = (2a - b)(b - a) ? 0, то есть особая точка (z, w) = (0, 0) не вырождена. В зависимости от знака произведения l1l2 эта особая точка будет либо седлом, либо узлом (см. рис. 1).

На каком же множестве допустимо использовать степенную замену с матрицей показателей C для того, чтобы сделать исходную особую точку элементарной?

Степенные замены переменных с невырожденной матрицей показателей определены и обратимы, вообще говоря, только во внутренности первого квадранта. Полученное после замены векторное поле может иметь на положительных координатных полуосях z и w некоторое количество особых точек помимо точки (z, w) = (0, 0), про которые мы ничего не знаем.

Вырежем квадрат P = {0 # z # e, 0 # w # e}, не содержащий других особых точек, кроме начала координат. Поскольку z = x2 / y, w = y / x, то в координатах (z, y) это множество представляет собой сектор (см. рис. 4).

Таким образом, замена переменных x = zw, y = = zw2, рассматриваемая в секторе S, соответствующем вершине диаграммы Ньютона, превращает исходное векторное поле в векторное поле с единственной элементарной особой точкой.

Как же исследовать исходное векторное поле в секторах S1 и S2 , которые являются дополнением сектора S до первого квадранта?

Поставим сектор S1 в соответствие нижнему ребру диаграммы Ньютона и сделаем в нем замену переменных

y = x2u

с матрицей Тогда в результате преобразования нижнее ребро становится вертикальным (см. рис. 3, в). Пересечение сектора S1 с некоторой окрестностью нуля превратится в прямоугольник где d - некоторое положительное число. В правой части системы, полученной из системы (3) после замены (6), выделим члены, соответствующие вертикальному отрезку - образу нижнего ребра (см. рис. 3, в). Тогда эта система запишется в виде

При замене переменных (6) векторные коэффициенты точек носителя преобразуются с помощью матрицы следующим образом:

Отсюда X1(u) = au, F1(u) = (b - 2a)u2 + h. Координаты особых точек векторного поля (7) на границе прямоугольника P1 определяются из уравнения Fl(u) = 0. Предположим, что (b - - 2a)h < 0. Тогда последнее уравнение имеет единственный положительный корень Матрица Якоби правой части (7) в точке (x, u) = (0, u0) имеет вид Ее собственные значения равны m1 = au0 и m2 = 2(b - 2a)u0 . Поэтому если a(b - 2a) > 0, то рассматриваемая особая точка - узел, а если a(b - 2a) < 0, то седло.

Аналогично исследуется сектор S2 , в нем делается замена x = uy. В результате верхнее ребро диаграммы Ньютона становится горизонтальным, сектор S2 превращается в прямоугольник а система, полученная после замены переменных, имеет вид

Пусть D = (c - d )2 - 4g(a - b) < 0, тогда система (8) не имеет в прямоугольнике P2 особых точек.

Делая в прямоугольниках, соответствующих ребрам и вершинам диаграммы Ньютона, обратные замены переменных, получим фазовый портрет в первом квадранте окрестности исходной особой точки. Для того чтобы построить фазовый портрет в любом другом квадранте, отражаем определенное в нем векторное поле V в первый квадрант и проделываем ту же процедуру. Нетрудно проверить, что в нашем примере для всех отраженных векторных полей получается та же картинка, что и для самого векторного поля V (см. рис. 4, в).

МОНОДРОМНЫЕ ОСОБЫЕ ТОЧКИ

Предположим теперь, что в нашем примере (b - - 2a)h > 0 и l1l2 < 0. Это означает, что особая точка в прямоугольнике P является седлом, а в прямоугольнике P1 нет особых точек. В этом случае траектории ведут себя так, как показано на рис. 4 зеленым цветом: они входят в каждый сектор через одну границу, а выходят через другую. Следовательно, если мы рассмотрим все траектории, начинающиеся, скажем, в точках оси OY вблизи нуля, то каждая из этих траекторий, сделав полный обход вокруг особой точки, возвращается опять на ось OY. Особые точки, в окрестности которых траектории ведут себя подобным образом, получили название монодромных особых точек. Описанный выше метод не позволяет судить о том, как поведет себя траектория в момент возвращения на ось OY : вернется ли она в ту же точку, из которой вышла, или в какую-нибудь другую.

В настоящее время доказано (Ю.С. Ильяшенко), что монодромная особая точка полиномиального векторного поля на плоскости может быть либо центром (все траектории замкнутые), либо фокусом (все траектории - спирали). Описанная выше схема не дает ответа на вопрос: является ли особая точка центром или фокусом?

Алгоритм, который позволил получить достаточные условия фокуса в случае монодромной особой точки с диаграммой Ньютона, содержащей более одного ребра, был впервые предложен Ю.С. Ильяшенко на основе методов, изложенных в мемуаре А. Дюлака "О предельных циклах". В настоящее время реализован лишь первый шаг этого алгоритма. Он состоит в следующем.

Рассмотрим гладкую кривую с началом в особой точке (трансверсаль), r - параметр на этой кривой. Преобразование D, переводящее точку r трансверсали в точку D(r) первого возвращения на трансверсаль траектории с началом в точке r, называется преобразованием монодромии монодромной особой точки. Доказано, что при подходящем выборе трансверсали преобразование монодромии имеет при x 0 асимптотику D: x cx + o(x). Если особая точка - центр, то D(r) ╞ r. Поэтому если c ? 1 или, что то же самое, ln c ? 0, то особая точка является фокусом.

Для того чтобы получить формулу для ln c, применяют метод А. Дюлака, который в случае раздутия по диаграмме Ньютона состоит в следующем. Преобразование монодромии разбивается в композицию отображений соответствия в прямоугольниках, соответствующих ребрам и вершинам диаграммы Ньютона. Отображение соответствия отображает одну сторону прямоугольника в другую вдоль траекторий определенного в нем векторного поля. На рис. 5 изображены прямоугольники, полученные после раздутия для всех четырех квадрантов. Стрелочками показаны отображения соответствия в них: f1 и f2 - отображения соответствия для прямоугольников P1 и P2 ; d - отображение соответствия для прямоугольника P. Отображения соответствия для отраженных векторных полей обозначены теми же буквами с соответствующими верхними индексами. Преобразование монодромии D представляет собой следующую суперпозицию:

Главный член отображения соответствия d в прямоугольнике P совпадает с отображением соответствия d1 для линейной системы

Вычислим это последнее отображение. Время, за которое фазовая кривая системы (10) переходит из точки (e, w) в точку (z, e), где z = d1(w), составляет

Интегрируя, получим, что

где e* - некоторая степень e,

Главный член отображения соответствия f1 в прямоугольнике P1 для системы (7) такой же, как у укороченной системы, которая получается из системы (7) отбрасыванием слагаемых с многоточием. Интегрируя уравнение

соответствующее укороченной системе, получаем, что главный член отображения f1(x) равен x c1x, где

Делая аналогичные вычисления для отраженных векторных полей, получаем, что главный член отображения dxy равен (10), а главные члены отображений dx и dy являются обратными к отображению (11). Далее главные члены отображений и линейны с коэффициентами и , причем

Аналогично для второго ребра

Вычисляя главный член суперпозиции (9), получаем, что он линеен за счет сокращения взаимнообратных степеней, а логарифм его коэффициента

ln c = u2 + lu1 .

Учитывая, что коэффициент c на самом деле не зависит от e, перейдем в (13) к пределу по e при e 0 и в силу нечетности подынтегральной функции в (12) получим, что

Итак, если в нашем примере при всех наложенных ранее условиях 2d(a - b) - b(c - d ) ? 0, то особая точка системы (3) является фокусом.

Метод Дюлака и теория нормальных форм позволяют вычислить следующие члены асимптотики преобразования монодромии, а значит, получить достаточные условия фокуса и в случае, если ln c = 0.

ЛИТЕРАТУРА

1. Мышкис А.Д. Вращение плоского векторного поля // Соросовский Образовательный Журнал. 1997. ╧ 11. С. 107-113.

2. Арнольд В.И. Обыкновенные дифференциальные уравнения. М.: Наука, 1984. 272 с.

3. Брюно А.Д. Локальный метод нелинейного анализа дифференциальных уравнений. М.: Наука, 1979. 256 с.

4. Березовская Ф.С., Медведева Н.Б. Асимптотика преобразования монодромии особой точки с фиксированной диаграммой Ньютона // Тр. сем. им. И.Г. Петровского. 1991. Вып. 15. С. 156-177.

5. Вишик М.И. Поля направлений и соответствующие им траектории // Соросовский Образовательный Журнал. 1996. ╧ 2. С. 111-117.

6. Белых В.Н. Элементарное введение в качественную теорию и теорию бифуркаций динамических систем // Там же. 1997. ╧ 1. С. 115-121.

* * *

Наталия Борисовна Медведева, кандидат физико-математических наук, доцент Челябинского государственного университета. Область научных интересов - обыкновенные дифференциальные уравнения, динамические системы. Автор 17 научных статей.