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

Лекции по дифференциальным уравнениям (Абрамов А.А

.).pdf
Скачиваний:
142
Добавлен:
12.03.2016
Размер:
3.56 Mб
Скачать

1. Основные определения и простейшие свойства

 

 

 

 

81

Существование решения задачи Коши следует из предыдущей главы.

 

Теорема: Пусть '~(t) решение системы (1), C const. Тогда

 

 

 

 

 

 

 

~(t) = '~(t + C) решение (1):

 

 

 

 

 

d~(t)

 

d'~(t + c)

 

~

 

~

 

 

~ ~

 

Доказательство:

dt

=

dt

t+c=

 

 

 

 

=

 

( )

 

 

 

 

= f '~(t) = f '~(t + c)

 

 

f

t

Теорема: Пусть '~(t) и ~(t) решения (1). Пусть '~(t1) = ~(t2). Тогда '~(t) = ~(t + t2 t1) для всех t, для которых обе части равенства определены.

Доказательство: !~(t) = ~(t+t2 t1). Тогда !~(t) решение (1), !~(t) = ~(t+t2 t1) = ~(t2) =

'~(t1).

n

!~(t) решение (1)

o^n!~(t1) = '~(t1)

) !~(t) = '~(t)o

'~(t) решение (1)

в силу единственности решения задачи Коши.

Определение: Область G называется фазовым пространством системы дифференциальных

~

уравнений (1), графики (в G) решений фазовыми траекториями, векторы f(~y) фазовыми

скоростями. Фазовое пространство вместе с фазовыми траекториями и фазовыми скоростями, называется фазовым портретом системы дифференциальных уравнений (1).

2.Положение равновесия автономной системы дифференциальных уравнений (1)

Определение: Точка A называется положением равновесия автономной системы дифференциальных уравнений (1), если ~y(t) ~a решение (1).

82 ГЛАВА 4. АВТОНОМНЫЕ СИСТЕМЫ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ

Теорема:

~a

положение равновесия

(1)

,

~

~.

 

 

 

 

 

 

 

 

 

 

d~a

 

 

f(~a) = 0

) : ~y(t) ~a

решение

(1)

,

 

 

~

. ~

~

 

dt

 

 

 

 

 

 

 

 

 

= f(~a)

f(~a) = 0

 

 

 

,

 

d~y

 

d~a

 

~

~

 

 

~ .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( : ~y(t) = ~a

 

dt =

dt

= 0 = f(~a) = f(~y)

 

 

2. Классификация положений равновесия линейной автономной системы 2 порядка.

В качестве области определения рассматриваем всю плоскость ~y.

1. Рассмотрим систему

 

 

 

 

 

 

d~y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

y2

, A = k a21

 

dt

= A~y + b;

.

A è B const. ~y =

a22

k; b =

b2

 

 

y1

 

новое

a11

a12

 

 

 

.

 

b1

 

~a положение

равновесия,

= ~y

старое

 

 

~a

 

 

 

 

~y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d~y

 

 

 

 

 

 

 

 

 

 

~ положение равновесия,

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

dt

= A~y

 

 

 

 

 

 

 

 

 

Далее 1; 2 собственные значения матрицы A.

1. 1; 2 вещественные, отличны от 0, разных знаков. Тогда сущестует вещественный базис

~ ~

h1; h2 из собственных векторов матрицы A.

1; 2 координаты в этом базисе. Наша система приводится к следующему виду:

d 1 = 1 1 dt

d 2 = 2 2 dt

1 = C1e 1t; xi2 = C2e 2t;

ãäå C1; C2 Const.

Положение равновесия, соответствующее этому случаю, называется седлом.

2. Классификация положений равновесия линейной автономной системы 2 порядка.

83

~~

2.1; 1 отличны от нуля, различны, одного знака. Тогда существует h1:h2 базис из собственных векторов A, 1; 2 координаты.

d 1 = 1 1 dt

d 2 = 2 2 dt

1 = C1e 1t; 2 = C2e 2t;

ãäå C1; C2 Const.

Мы рассмотрим два разных случая.

(a) 1, 2 отрицательны, для определ¼нности 1 < 2 < 0.

lim

1(t)

=

C1

lim

e 1t

= 0

 

 

 

 

 

2(t)

 

 

 

t!+1

 

C2 t!+1 e 2t

 

Такое положение равновесия называется устойчивым узлом.

(b)1; 2 положительны, для определ¼нности 0 < 2 < 1.

В этом случае надо поменять t ! t, и получается предыдущий случай. Полу-

чится та же самая картинка, движение будет осуществляться в противоположном направлении.

84ГЛАВА 4. АВТОНОМНЫЕ СИСТЕМЫ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ

3.1; 2 не вещественные, с вещественной частью, отличной от 0.

1 = + i

;

6= 0

2 = i

 

6= 0

Поскольку мы хотим рисовать вещественным мелом на вещественной доске, делаем следующее: составляем вектор

 

~

 

 

 

 

 

h1 = ~g + i~g2;

вектор, соответствующий 1. ~g1;~g2 вещественные.

Тогда ~

вектор, соответствующий 2. ~g1;~g2 линейно независимы, f~g1;~g2g

h2 = ~g1 i~g2

базис, 1; 2 соответствующие координаты.

В координатах 1; 2 наша система привед¼тся к такому виду:

 

 

 

d 1

= 1 + 2

 

 

 

dt

 

 

 

 

 

 

d 2

 

= 1 + 2

 

 

 

dt

1(t) = C1et sin( t + C2)2(t) = C1et cos( t + C2);

C1; C2 const.

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

Далее для определ¼нности > 0. Если оно отрицательное, поменяем местами 1; 2.

(a) < 0

Устойчивый фокус.

(b) < 0

2. Классификация положений равновесия линейной автономной системы 2 порядка.

85

Неустойчивый фокус.

Весь случай целиком называется фокус.

4. 1; 2 невещественные, с нулевой вещественной частью.

1 = i2 i

Как и в пункте 3), строим базис ~g1;~g2, получаем

d 1 = 2 dt

d 2 = 2 dt

1 = C1 sin(t + C2)2 = C1 cos(t + C2)

Для определ¼нности > 0

С точки зрения выкладок это частный случай предыдущего пункта, поэтому выкладки не повторяются.

5.1 = 2 6= 0.

6.1 2 = 0.

86 ГЛАВА 4. АВТОНОМНЫЕ СИСТЕМЫ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ

Замечание: В книге Понтрягина детально разобраны все случаи.

2.Дополнительные замечания

1.Мы рисовали все картинки в том базисе, к которому мы перешли. В исходном базисе ~y1; ~y2, с уч¼том того, что была сделана замена переменной t, картники будут выглядеть так:

Рисуя вс¼ в исходной системе координат, надо тонко уч¼ть вс¼. что связано с переходом

в другую систему координат.

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

Эта классификация называется классификацией Пуанкаре линейных однородных систем дифференциальных уравнений второго порядка.

3.Все эти случаи разбиваются на два класса: 1)-3) так называемые невырожденные положения равновесия, и 4)-6) так называемые вырожденные положения равновесия.

В ч¼м смысл этого разбиения? Казалось бы, в случае 4) было получено вс¼, и выкладки являлись частным случаем случая 3). Дело в том, что каджое из невырожденных положений равновесия выделяется какими-либо условиями типа неравенства. А в определении каждого из вырожденных положений равновесия фигугирует какое-нибудь условие типа равенства.

Это важно потому что если положение равновеися невырожденное, и мы достоаточно мало меняем исходные данные, то тип положения равновесия не меняется, потому что при достаточно малом изменении исходных параметров эти значения меняются мало (непрерывные функции), а все ограничения в виде неравенств. Значит, тип невырожденного положения равновесия не меняется при достаточно малом изменении исходных параметров.

С другой стороны, для вырожденного положения равновесия это не выполняется. Центр может остаться центром, а может превратиться в фокус.

3. О нелинейных автономных системах дифференциальных уравнений

87

3. О нелинейных автономных системах дифференциальных уравнений

1.Мы рассматриваем общую нелинейную систему

 

 

 

 

 

; f~ =

:f:1: ;

 

d~y

= f~(~y); ~y =

 

y.1

 

(1)

dt

 

 

 

fn

 

 

 

 

yn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

~

 

~

 

~

 

@f

 

 

 

t; ~y; f вещественные, f(~y) определена в области G, f,

 

@~y

непрерывна в G

Åñëè

~a

положение равновесия, то

~y(t) ~a

решение

,

 

~

~.

 

 

 

 

 

 

 

 

 

 

f(~a) = 0

Часто используется следующий при¼м: ~y = ~a + ~z, j~zj мало. Тогда

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

@f

 

 

 

 

 

 

 

f~(~y) = f~(~a + ~z) = f~(~a) +

 

 

 

 

(~a)!~z + ~g(~z); j~g(~z)j = o(j~zj):

 

 

 

 

@~y

 

 

 

=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

|{z}

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

|

 

{z

 

}

 

 

 

 

 

d~z

= A~z + ~g(~z);

(2)

dt

 

 

~

A = @@~fy (~a); j~g(~z)j = o(j~zj).

(1) , (2).

d~z

= A~z

(3)

dt

 

 

Определение: (3) получено из (1) линеаризацией в окрестности положения равновесия ~a.

Один из основных при¼мов исследования поведения решений: прежде всего изучаем уравнение

(3), и выясняем нужные нам свойства на этом уравнении, затем стараемся перенести эти результаты на уравнение (2).

2.Пусть n = 2. (Линейная автономная система второго порядка).

Имеет место следующий нетривиально доказываемый факт:

Траектории уравнения (1) в окрестности положения равновесия ведут себя , если положения равновесия системы (3) невырожденное, качественно так же, как и траектории уравнения (3).

88 ГЛАВА 4. АВТОНОМНЫЕ СИСТЕМЫ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ

3. Устойчивость по Ляпунову

Определение: Решение ~y(t) ~a и само положение равновесия называется устойчивым по Ляпунову, если:

1. 9q > 0, что при j~y(0) ~aj < q решение может быть продолжено на [0; +1),

2. для любого " > 09 > 0; 6 q, такое, что если j~y(0) ~aj < , то j~y(t) ~aj < " при

0 6 t < +1.

Определение: Устойчивое по Ляпунову решение ~y(t) = ~a уравнения (1), и само положение равновесия ~a называется асимптотически устойчивым, если существует такое 0 < r 6 q

(обозначение q см. в предыдущем определении), что при ~y

 

~a

j

< r имеет место

lim ~y(t) =

j (0)

 

 

t!+1

~a.

Упражнение: Рассмотрите те типы, которые были рассмотрены в x2, и выясните, какие из них являются устойчивыми, а какие к тому же асимптотически устойчивыми.

Прежде всего, мы исследуем систему (3).

Теорема: Пусть все собственные значения матрицы A лежат в левой полуплоскости. Тогда

~

решение системы (3) ~z(t) 0 устойчиво по Ляпунову и асимптотически устойчиво.

Доказательство: ~z(t) = etA~z(0) (см. 10 главы 1).

1

; B =

 

J1 .

..

0

 

; Js =

 

1

A = CBC

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

Jm

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

10

... ......

 

1

 

 

 

 

 

 

 

 

 

0s

| {z }

ks

 

 

 

etJ1

0

 

 

 

etA = CetBC 1; etB =

...

;

 

 

 

 

0

etJm

 

 

 

3. О нелинейных автономных системах дифференциальных уравнений

89

1

etJs = e st

0

0

t!=1 : : :

... ...

...

0

tks 1

(ks 1)!

.; Re s 6 < 0

t!=1

Нам важно лишь только то, что там стоят какие-то многочлены.

< ~ < 0; jetJs j 6 je stjKe(~ )t 6 e tKe(~ )t = = Ke~t; jetAj 6 jCjjetBjjC 1j 6 Qe~t:

Итак, имеем такую оценку: j~z(t)j 6 Qe~tj~z(0)j; ~ < 0.

Центральная теорема этого параграфа. Мы возвращаемся к уравнениям (1), (2).

 

 

 

 

 

~

 

 

 

Теорема: (Ляпунов). Пусть все собственные значения матрицы

@f

(~a) лежат в левой полуплос-

@~y

кости. Тогда решение

~y(t) ~a

системы

(1)

(или, то же самое)

~

системы

(2)

устойчиво

 

 

 

~z(t) = 0

 

 

по Ляпунову и асимптотически устойчиво.

Замечание: Доказательство нужно проводить, используя уже более тщательные оценки.

Напоминание:

 

d~z

 

 

 

 

 

 

= A + ~g(~z); j~g(~z)j = o(j~zj):

(2)

 

 

 

 

 

dt

Доказательство: 80 < 9 > 0:

из j~zj 6 j j следует j~g(~z)j 6 .

 

Пока решение нашей задачи существует при 0 6 6 t, и j~g(~z)j 6 j~zj, (см. 10 главы 1)

 

 

 

 

 

 

 

 

 

~z(t) = etA~z(0) + Z0

e(t )A~g( ) d

 

Из доказательства предыдущей теоремы заимствуем оценку:

 

 

 

 

t

 

 

 

 

j~z(t)j 6 Qe~tj~z(0)j + Z0

Qe(t )~ j~z( )j d

 

 

e ~tj~z(t)j 6 Qj~z(0)j + Q Z

t

 

 

e ~j~z( )j d

 

0

По неравенству Гронуолла (см. 1 главы 3),

e ~tj~z(t)j 6 Qj~z(0)jeQ t

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

j~z(t)j 6 Qj~z(0)je(~+Q )t

90 ГЛАВА 4. АВТОНОМНЫЕ СИСТЕМЫ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ

Выберем так, чтобы ~ + Q < 0. (Чтобы отрицательное число ~ мало сдвинулось вправо.) Для этого надо взять соответствующее . Теперь взяли ~, обладающее таким свойством:

Q~ < , ~ < .

Начальные значения подчиняются условию:

j~z(0)j < ~

(4)

Если для 0 6 6 t решение ~z(t) определено, имеет место j~zj 6 и выполняется условие (4), то

j~z(t)j 6 Qj~z(0)je(~+Q)t

(5)

Сейчас мы докажем, что решение действительно определено, прич¼м для всех t, и для всех решений ~z не выходит за указанные пределы. (Иначе говоря, эти предположения выполняются).

Возьм¼м цилиндр S, такой, что

S : f0 6 t < +1; j~zj 6 g

Задача Коши с такими начальными данными имеет решение. Теперь докажем, что оно не выходит за пределы цилиндра.

 

g; 0 < T любое

S : f0 6 t 6 T; j~zj 6

Èç

(5)

получаем, что

~

.

j~z(t)j 6

.

 

 

j~z(t)j <

 

 

Теорема доказана.

Замечание: Докажите, что если хотя бы одно собственное значение матрицы лежит в правой полуплоскости, то решение не является устойчивым по Ляпунову. Подсказка: ограничиться случаем уравнения (3).