Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Понтрягин, Л. С. Обыкновенные дифференциальные уравнения учебник

.pdf
Скачиваний:
28
Добавлен:
22.10.2023
Размер:
11.65 Mб
Скачать

120

ЛИНЕЙНЫЕ УРАВНЕНИЯ С ПОСТОЯННЫМИ КОЭФФИЦИЕНТАМИ

[Гл. 2

перешел в единицу, а

вектор — в /;

тогда

вектору

S1 fti 4-

будет

соответствовать

комплексное число

t =

В

силу

этого

£

отображения фазовая траектория

(9)

перейдет

в

фазовую

траекто­

рию

на плоскости

Р*, описываемую уравнением

 

 

 

 

 

 

 

 

С=

сеи .

 

 

 

 

 

(10)

В) Ф о к у с

и

це нтр .

Перепишем

уравнение

(10)

в

полярных

координатах, положив

 

 

 

 

 

 

 

 

 

 

 

 

С= ре'*, с = Яеи .

 

 

 

 

Таким образом

получаем:

 

 

 

 

 

 

 

 

 

 

 

 

р = Я е ^ ,

= v /

- f -

 

 

 

 

это

есть уравнение движения

точки

в

плоскости

Р*.

При р 0

каждая траектория

оказывается

л о г а р и ф м и ч е с к о й

с п и р а л ь ю .

Соответствующая

картина

на

 

плоскости Р

называется

фокусом.

Если

р<^0, то

точка при

возрастании t

асимптотически

приближает­

ся к началу координат, описывая логарифмическую спираль. Это— ус­ тойчивый фокус (рис. 27, а). Если р^>0, то точка уходит от на­ чала координат в бесконечность, и мы имеем неустойчивый фокус (рис. 27, б). Если число р равно нулю, то каждая фазовая траектория, кроме положения равновесия (0, 0), замкнута, и мы имеем так назы­ ваемый центр (рис. 28).

1 161

ФАЗОВАЯ ПЛОСКОСТЬ ЛИНЕЙНОЙ СИСТЕМЫ

121

Рисунки 27 и 28 дают картину во вспомогательной фазовой пло­ скости; в плоскости Р картина аффинно искажается (см., например,

рис. 29 и 30).

 

Рис. 27.

Выше

мы рассматривали так называемые н е в ы р о ж д е н н ы е

с луча и :

корни Xj и Х4 различны и отличны от нуля. Малое

изменение элементов матрицы (а*) не меняет в этих предположениях общего характера поведения фазовых траекторий. Исключение состав­ ляет случай центра: при малом изменении элементов матрицы (а*)

122 ЛИНЕЙНЫЕ УРАВНЕНИЯ С ПОСТОЯННЫМИ КОЭФФИЦИЕНТАМИ [Гл. 2

равенство jj.= 0 может нарушиться, и центр перейдет в устойчивый или неустойчивый фокус. Включение этого вырожденного случая (центра) в основной текст па­ раграфа объясняется его важ­ ностью. Остальные вырожден­ ные случаи будут рассмотрены

в примерах 1 и 3.

П р и м е р ы

1.

Если матрица А системы (1) имеет лишь одно собственное значение X, то возможны два существенно различных слу­ чая, при описании которых

 

 

 

мы будем обозначать

через

А

преобразование,

соответствующее матрице А.

 

 

 

 

С л у ч а й

I.

Существует в плоскости Р базис

ft,, fts,

состоящий

из двух собственных векторов преобразования А:

 

 

 

 

 

 

Ah, — Xft,,

Ah<i — Xft2.

 

 

(11)

С л у ч а й

11.

Существует в плоскости Р такой

базис

ft„

ft2,

что

 

 

Ah, = Xft„

А!ц — X/i,2 -j- h,.

 

 

(12)

Существование базиса одного из видов (11), (12) непосредственно вытекает из теоремы 30, но здесь мы докажем этот факт непосред­ ственно. Пусть h 1 — собственный вектор преобразования А и ft2 — произвольный вектор, не коллинеарный вектору Н,. Тогда мы имеем:

 

 

 

Ah, = Xft,,

Ahi =

aft, -f-

 

Из этого

видно, что преобразование

А

имеет в базисе h b ft3

матрицу

 

 

 

 

X

а\

 

 

 

 

 

 

0

[ А / ’

 

 

так

что

его

собственными

значениями

являются X и ц, и потому

ji =

X. Если

<х= 0, то для базиса ft„

ft2

выполнены соотношения (11).

Если же а. ф 0, то, заменив вектор ft,

коллинеарным ему вектором ah,,

мы получим

базис, удовлетворяющий условию (12).

 

 

Непосредственно проверяется, что в случае I общее решение

уравнения (2) записывается

в виде:

 

 

 

 

 

 

X =

 

c*ft2ex' =

(13)

§ 16]

ФАЗОВАЯ ПЛОСКОСТЬ ЛИНЕЙНОЙ СИСТЕМЫ

123

Написанное

решение имеет начальное значение (0, х 0).

При X О

каждое решение описывает полупрямую, выходящую из начала коор­

динат. При Х<[0 движение происходит в направлении

к началу

ко­

ординат (рис. 31, а), при

Х ^>0— от начала координат

(рис. 31,

б);

относительно случая X=

0 см. пример 3.

 

 

Рис. 31.

Непосредственно проверяется также, что в случае 11 произвольное решение уравнения (2) имеет вид:

х = c'hieu - f с4 {hit - f hi) ex‘,

Разлагая это решение по базису fc„ в виде х = £’ft| -j- £4ft_> по­ лучаем уравнение траекторий в плоскости Р относительно базиса h {, Л4:

V — (с1 + cH) eu,

V2 =

cV '.

(14)

Аффинное отображение фазовой плоскости

Р,

переводящее

векторы

hi и 11г в единичные ортогональные

векторы,

направленные

по осям

координат плоскости Р*, переводит траектории плоскости Р в траек­

тории плоскости Р*, в которой

уравнения (14) дают траектории

уже

в прямоугольных координатах.

 

 

 

 

 

 

 

 

Разберем случай Xф 0 (случай X= 0 см. в примере 3). Пусть

сначала Х<^0. Рассмотрим траектории,

заполняющие

плоскость

Р*,

в этом случае. Прежде всего

из уравнений (14)

видно, что, меняя

одновременно

знаки у с1 и с2, мы получим симметрию плоскости Р*

относительно

начала координат,

при которой траектории

переходят

в траектории. Таким образом,

достаточно

рассмотреть заполнение

траекториями

верхней полуплоскости.

При

с4 = 0,

с1 ф 0

получаем

Две траектории: одну при с,

0,

другую — при с,

0.

Первая совпа­

дает с положительной полуосью

абсцисс, вторая — с

отрицательной

полуосью абсцисс; движение по обеим направлено к началу координат. Рассмотрим, далее, траекторию Cj = 0, сг 0. Мы имеем:

(15)

124

ЛИНЕЙНЫЕ УРАВНЕНИЯ С ПОСТОЯННЫМИ КОЭФФИЦИЕНТАМИ

[Гл. 2

При

t =

О получаем

точку (0, са) на оси ординат. При t, возра­

стающем

от нуля, точка сначала движется направо, затем налево,

все время

опускаясь

вниз к началу

координат, к которому она

под­

ходит

по

траектории,

касающейся

положительного направления

оси

Рис. 32.

абсцисс. При t, убывающем от нуля до — оо, точка движется налево, одновременно поднимаясь вверх, однако налево быстрее, чем вверх, так что общая тенденция ее движения — в отрицательном направлении вдоль оси абсцисс. Если в уравнениях (15) придавать константе с*

все положительные значения, то описанные таким образом траектории заполняют всю верхнюю полуплоскость (рис. 32, а). Мы имеем здесь ус­ тойчивый вырожденный узел. Если Х^>0, то траектории полу­ чаются из описанных путем зеркального отражения плоскости относи­ тельно оси ординат (рис. 32,<5), а движение по ним идет в противо­ положном направлении, т. е. от начала координат. Это — неустойчивый

§

16]

 

ФАЗОВАЯ ПЛОСКОСТЬ ЛИНЕЙНОЙ СИСТЕМЫ

125

вырожденный

узел.

На рис.

33, а, б показаны фазовые траектории

в

плоскости Р.

 

 

 

 

 

2.

Рассмотрим

линейное

однородное уравнение второго

порядка

с постоянными

коэффициентами

 

 

 

 

 

 

Л -)- ах -f- Ьх = 0.

(16)

Заменяя это уравнение нормальной системой по способу, изложенному в § 4, получаем:

х— у,

уЬх ау.

Фазовую плоскость системы (17) считают фазовой плоскостью урав­ нения (16). Непосредственно проверяется, что характеристический многочлен системы (17) совпадает с характеристическим многочленом уравнения (16), т. е. равен

Р* + ар + Ь.

(18)

Таким образом, если корни многочлена (18) комплексны, то фазовая плоскость уравнения (16) представляет собой фокус или центр. Рас­ смотрим фазовую плоскость в случае действительных различных и отличных от нуля корней многочлена (18).

Пусть X— корень многочлена (18) и h ~ ( h l, h*)— соответствую­ щий ему собственный вектор. Мы имеем тогда (учитывая вид системы

(17))

= ХА1.

Таким образом, собственное направление, соответствующее собствен­ ному значению X, определяется прямой линией, имеющей уравнение

у = Хх;

мы будем называть ее собственной прямой.

Если корни X, и Х4 отрицательны, то мы имеем устойчивый узел (см. А)). В этом случае обе собственные прямые проходят во второй и четвертой четвертях; траектории вблизи начала координат касаются той из них, которая расположена ближе к оси абсцисс (рис. 34).

Если корни Xt и Ха положительны, то мы имеем неустойчивый узел (см. А)). Обе собственные прямые проходят в первой и третьей четвертях; вблизи начала координат траектории касаются той из них, которая расположена ближе к оси абсцисс (рис. 35).

Если корни X, и Х4 имеют разные знаки,

то мы

имеем седло;

одна собственная прямая проходит во второй

и четвертой

четвертях,

а другая — в первой и третьей. В направлении

первой

из

этих пря­

мых траектории приближаются к началу координат, а в направлении второй — отходят от него (рис. 36).

126

ПИНЕЯНЫЕ УРАВНЕНИЯ С ПОСТОЯННЬШИ КОЭФФИЦИЕНТАМИ ГГл. 8

3.Рассмотрим, наконец, случай, когда хотя бы одно из собствен

ных значений матрицы А обращается в нуль.

С л у ч а й I. В нуль обращается лишь одно собственное значение:

ф О,

Х.2 = 0.

 

 

 

 

В этом случае решение можно записать

в

виде

(4), где

и S-1

даются формулами (5). Так как Х.2 =

0, то

=

const,

и движение

про­

исходит

по прямым £2 =

const в направлении к прямой

= 0 или от

нее

в зависимости от

знака числа Xt. Все точки прямой

£* =

0

являются положениями равновесия (рис. 37).

 

 

 

 

Если

же имеется единственное собственное значение

X, =

Х9 =

0,

то

могут

представиться

два случая, рассмотренные в примере

1.

 

§ 16] ФАЗОВАЯ ПЛОСКОСТЬ ЛИНЕЙНОЙ СИСТЕМЫ 127

С л у ч а й I (см. (11), А = 0). Общее решение записывается в виде

(см. (13)):

х = х л.

Этот случай имеет место, если все коэффициенты системы (1) обра­ щаются в нуль; каждая точка плоскости Р является положением равновесия.

 

С л у ч а й II (см. (12),

А= 0). Общее решение записывается в виде

(см.

(14)):

 

с1 -{- cs t,

«= с*.

 

 

 

 

 

Движение происходит

равномерно по каждой из

прямых ^ = const.

Все

точки прямой £4 =

0

являются

положениями

равновесия (рис. 38).

Г Л А В А Т Р Е Т Ь Я

ЛИНЕЙНЫЕ УРАВНЕНИЯ С ПЕРЕМЕННЫМИ КОЭФФИЦИЕНТАМИ

В этой главе излагается теория линейных уравнений, сперва для нор­ мальной системы п-го порядка, а затем для одного уравнения я-го по­ рядка, причем почти все результаты, относящиеся к одному уравне­ нию, выводятся из соответствующих результатов о нормальной си­ стеме. Третий параграф главы посвящен нормальным линейным одно­ родным системам с периодическими коэффициентами. Главным результатом здесь является теорема Ляпунова о возможности линей­ ным периодическим преобразованием переменных перевести нормаль­ ную систему с периодическими коэффициентами в нормальную систе­ му с постоянными коэффициентами. В дальнейшем этот результат находит важное применение в теории устойчивости. Доказательство его очень просто, но опирается на сравнительно неэлементарную те­ орию функций от матриц. Эта теория, не являющаяся частью теории обыкновенных дифференциальных уравнений, излагается для удобства читателей в добавлении II. Таким образом, третий параграф этой главы (§ 19) является иеэлементарным благодаря используемому в нем мат­ ричному исчислению.

§ 17. Нормальная система линейных уравнений

 

Здесь будет

рассмотрена

нормальная система линейных уравнений

с переменными

коэффициентами

 

 

 

 

^

=

 

+

0 .

 

 

( 1)

 

;=1

 

 

 

 

 

 

Напомним, что

если

 

 

есть интервал

существования

и не­

прерывности коэффициентов

a1, (t) и свободных

членов b‘(t) системы

(1), то, в силу

теоремы

3,

каждое решение

может быть продолжено

на весь интервал

<С'Т* в

дальнейшем

мы будем считать,

что

каждое рассматриваемое решение задано на этом интервале и каждое рассматриваемое значение t принадлежит ему.

НОРМАЛЬНАЯ СИСТЕМА ЛИНЕЙНЫХ УРАВНЕНИЯ

129

Ф у н д а м е н т а л ь н а я с и с т е м а р е ш е н и й

В первую очередь будет рассматриваться однородная система

уравнений

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

х г =

2

a‘,{t)xJ,

 

/ = 1 ,

п,

 

 

(2)

 

 

 

 

 

 

/= |

 

 

 

 

 

 

 

 

 

или, в

векторной

записи,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X — A (t) X.

 

 

 

 

 

(3)

А)

Установим

простейшие

свойства

уравнения (3).

 

 

а)

Если

х =

<р (t)

есть

решение

 

уравнения (3),

обращающееся в

нуль при

некотором

значении t0\

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ф(*о) = 0,

 

 

 

 

 

(4)

то решение

это

тождественно равно нулю

 

 

 

 

б)

Если

 

 

ф (0 =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ф КО.

ф » ( 0 .

• ••>

фг(0

 

 

 

— решения

уравнения (3),

то

векторная

функция

 

 

 

 

 

 

 

 

ф ( 0 =

с1ф 1 ( 0 - Ь

•••

~ Ь сГфг(0>

 

 

 

где с1, ...,

сг — константы,

также

является

решением

уравнения (3).

Свойство б) проверяется непосредственно. Свойство а)

вытекает

из того,

что вектор

л; =

0,

тождественно

равный

нулю,

очевидно,

является

решением уравнения (3), а потому решение <р(£), предусмотрен­

ное в

а),

как имеющее с этим решением общее начальное условие (4),

должно с ним совпадать.

 

 

 

 

 

 

 

 

 

 

 

Б)

Пусть

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ф 1 ( 0 .

ф « ( 0 .

• • • >

Ф г ( 0

 

 

 

( 5 )

— система

решений уравнения (3). Она

 

называется

линейно зависи­

мой, если существуют такие

константы

с1,

с2,

сг,

не

Есе рав­

ные нулю,

что

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c 4 i ( 0 +

c 4

* ( 0 + . . .

+

с г ф г ( 0 =

0 .

 

 

В противном случае система (5) решений уравнения (3) называется линейно независимой. Оказывается, что если хотя бы для одного значения t = ta векторы

Ф 1< А ) . ф * ( < о ) . . . . » ф , ( * о ) ( е )

5 Понтрягин Л . С.

Соседние файлы в папке книги из ГПНТБ