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

Лекции, Простокишин В.М. / Электронный формат / ДиИУ_1_Теория_устойчивости_решений_задачи_Коши_Л_03_04_

.pdf
Скачиваний:
0
Добавлен:
21.05.2025
Размер:
3.53 Mб
Скачать

Л.03-04

1.2. Качественные исследования точки покоя простейшей системы ОДУ в случае n=2

Рассмотрим линейную систему с второго порядка постоянными коэффициентами

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

a

 

x

a

 

x

 

 

 

1

11

1

12

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

2

a

 

x

a

22

x

2

 

 

 

 

 

21 1

 

 

Точка ≡ (

,

) = (0 , 0) = O является точкой покоя этой системы.

 

 

 

1

2

 

 

 

 

 

 

 

 

 

 

В данном параграфе будем считать, что (−∞ , +∞). От системы (11) можно перейти к ОДУ первого порядка:

(11)

dx

 

a

x a

x

 

2

21

1

22

2

(12)

 

dx

 

a

x a

x

 

 

1

 

11

1

12

2

 

Фазовая траектория системы (11) является интегральной кривой уравнения (12). При этом точка ≡ ( 1 , 2) = (0 , 0) для уравнения (12) является особой точкой , так как это – точка разрыва правой части уравнения (12) и в ней нарушается условие теоремы существования и единственности решения: точка покоя системы (11) является особой точкой уравнения (12).

Исследование точки покоя O на устойчивость.

 

 

 

h

 

 

 

 

 

h

 

 

 

 

 

 

1

 

 

 

 

 

 

 

2

 

 

1

 

 

1

 

,

h

2

 

 

1

 

1) Решение определяется собственными числами 1, 2 матрицы системы и собственными векторами

h

 

1

 

 

 

h

2

.

 

 

 

 

 

 

 

 

 

 

h

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

2

2) Характеристическое уравнение

det( A E) 0

 

a11

a12

0

2 SpA det A 0

,

 

 

a21

a22

 

 

1 2 ,

 

 

 

 

 

 

3) Собственные числа могут быть:

а) действительными, при этом по знаку – положительными, отрицательными, нулевыми;

салгебраической кратностью 1 или 2;

сгеометрической кратностью 1 или 2,

б) комплексными, при этом

с действительной частью – положительной, отрицательной или нулевой.

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

Рассмотрим различные возможные случаи.

I.

1 , 2 , 1

2

и одного знака.

Пусть 2

<

1 < 0,

т.е. | 1| <

| 2|.

Общее решение системы (11):

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

X

 

1

 

 

1

t

С

 

h

2

e

t

С

 

1

 

e

t

 

 

 

 

С h e

1

2

 

2

 

1

 

1

 

 

1

 

 

 

 

 

1

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

h

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

h

 

 

 

 

 

h

 

 

 

 

 

 

1

 

 

 

 

 

 

 

2

 

 

1

 

 

1

 

,

h

2

 

 

1

 

 

h

 

1

 

 

 

h

2

 

– собственные векторы, отвечающие собственным числам 1, 2

 

 

 

 

 

 

 

 

 

h

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

2

 

асимптотически устойчиво (так как оба собственных числа отрицательны).

 

h

2

 

 

 

 

 

 

 

 

 

С

 

1

 

t

 

 

 

e

 

 

 

 

 

 

2

(13)

2

 

h

2

 

 

 

2

 

 

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

Положим С1 = 0. Тогда

x

 

2

e

t

,

С h

2

 

 

 

 

 

1

2

1

 

 

 

x

 

2

t

С h

e

 

 

 

2

2

2

2

 

параметрическое уравнение двух лучей, коллинеарных собственному

вектору

h

2

 

 

 

. Эти лучи «входят» в точку = O при

+∞ , данная точка покоя асимптотически устойчива (они

приближаются к точке покоя = O сколь угодно близко).

Назовем эти лучи фазовыми траекториями II (т.к. С2 ≠ 0,

Аналогично, если положить С2 = 0. Тогда

x

 

1

t

,

С h e

1

1

1

1

 

С1 = 0).

x

 

1

t

С h e

 

 

 

1

2

2

2

 

и получаем еще два луча, коллинеарные

собственному вектору

h

1

,

также «входящие» в точку = O при +∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

которые назовем фазовыми траекториями I

(т.к. С1 ≠ 0, С2 = 0).

 

 

Пусть теперь С1 ≠ 0, С2 ≠ 0, тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

e 1t

 

 

 

2

e 2t

X x1

С h1 e 1t С h2 e 2t

С h1

С h1

 

 

1

 

1

2

 

2

 

 

 

1 1 1

 

 

 

2 2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h2

 

 

 

h2

 

x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Откуда получаем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

С h1e 1t

С

 

h2e 2t

 

 

 

 

 

 

 

 

 

 

2

 

 

1

1

2

 

2

2

2

 

 

 

 

 

 

 

 

 

 

 

dx

 

С h1e 1t

С

h2e 2t

 

 

 

 

 

 

 

 

 

 

1

 

1

1

1

 

2

2

1

 

 

 

(приближающиеся к ней сколь угодно близко),

 

 

 

 

 

1 t

 

 

2 t

 

x1

С1 1h1 e 1

С2 2 h1 e 2

,

 

 

 

 

 

 

 

(14)

 

x С h1e 1t С h2e 2t

 

 

2

1

1

2

2

2

2

(15)

Из уравнения (15) с учетом, того, что 2 1 < 0 следует, что

dx

 

 

 

1

 

 

2

e

 

С h

С h

2

1

1

2

2

2

2

 

 

 

dx

 

 

 

1

 

 

2

e

 

С h

С h

1

 

1

1

1

2

2

1

 

( )t

 

1

2

1

 

h

 

 

 

2

 

 

 

( )t

t

1

2

1

h

 

 

 

 

 

1

 

dx

 

1

 

 

 

h

 

lim

2

 

2

 

 

1

 

t dx

 

 

 

h

 

 

1

 

1

 

– фазовые траектории системы (11), отличные от траекторий I и II «входят», в точку = O при

+∞, имея общую

касательную, коллинеарную вектору

h

1

 

 

 

, т.е. касаются фазовых траекторий I, отвечающих меньшему по модулю собственному

числу 1.

В случае 2 > 1 > 0 получаем аналогичную ситуацию, только с той разницей, что точка покоя = O является неустойчивой(так как оба собственных числа положительны), и все фазовые траектории «выходят» из точки покоя (из точки сколь угодно близкой к точке покоя).

Фазовый портрет называется узлом в случаях

1, 2 ,

1 2,

1 2 > 0

действительные, неравные друг другу собственные числа одного знак

устойчивый узел, при 1 , 2 < 0

неустойчивый узел при 1 , 2 > 0 .

Рисунку соответствует | 1| < | 2|

Замечание. При действительных и неравных друг другу собственных числах существует базис из собственных векторов матрицы системы. В этом базисе матрица системы – диагональна и система (11) приводится к виду

которая имеет общее решение уравнению траекторий:

y

 

2

 

y (t) A e

t

,

1

 

 

 

 

 

1

1

 

 

 

A ( y

 

/

/ A )

2

 

1

 

 

 

 

2

1

1

 

 

 

y

t

 

A e

 

 

2

;

2

2

,

при A

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

y

 

 

1

1

1

,

 

 

 

 

 

y

2

y

2

 

 

 

2

 

 

A , A

– произвольные постоянные.

1

2

0 и y

0

при A 0 .

1

 

1

(11А)

Исключая t, приходим к

(11Б)

 

 

 

 

 

 

 

 

 

 

 

 

Оy

в точке (0; 0) если |

 

 

|

 

 

 

 

 

 

Если у собственных числе знаки совпадают, то траектории – дуги, касающиеся оси

 

1

<

 

|

 

2

| , или оси

1

 

Оy

 

в той же точке (0; 0) если |

 

 

|

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

1

>

 

|

 

2

| , при этом траекториями являются и четыре координатные полуоси и сама точка (0;

 

 

 

 

0). Такая точка покоя (точка равновесия) системы (11) (для уравнения (12) – критическая точка) называется узлом. И как мы уже отметили ранее при 1 , 2 < 0 устойчивый узел (асимптотически устойчивый),

при 1 , 2 > 0 –– неустойчивый узел.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

разных знаков (т. е.

 

 

 

 

 

 

 

 

 

0). Пусть

 

 

 

 

 

 

 

 

 

 

 

. Очевидно, что положение равновесия =

II.

 

 

 

,

 

 

 

 

 

 

,

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

<

 

 

 

<

0

<

 

 

 

 

 

1

 

2

1

2

 

1

2

2

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

O является неустойчивым (так как = 1 > 0 ). Общее решение системы (11) описывается той же формулой (13), и из тех

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

h1

 

 

2

h2

 

же соображений есть фазовые траектории I и II (лучи, коллинеарные векторам

 

h

11

,

h

 

12

, соответственно).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h2

 

 

 

h2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Но теперь, с учетом знаков собственных чисел, фазовые траектории I «выходят» из точки покоя = O (т.к. 1 > 0), а фазовые траектории II «входят» в точку покоя = O (т.к. 2 < 0).

Параметрические уравнения фазовых траекторий I и II имеют вид

 

 

 

1

t

 

 

 

2

 

t

x (t) С h e

1

x (t) С h

e

2

 

1

1

1

 

 

 

1

2

1

 

 

 

I :

 

 

1

 

t

II :

 

 

2

 

t

x

 

 

x

 

e

(t) С h e

1

(t) С h

 

2

 

2

1

2

 

 

 

2

2

2

 

 

 

Тогда с учетом (13) получаем:

 

 

 

 

 

2

e

t

0

lim (x (t) x (t)) lim С h

2

 

 

 

 

 

 

 

 

t

1

1

t

2

1

 

 

 

 

 

 

 

 

 

 

.

(16)

lim (x

(t) x

 

 

2

e

t

0

(t)) lim С h

2

 

 

 

 

 

 

 

 

t

2

2

t

2

2

 

 

 

 

 

 

 

 

 

 

поскольку 2 < 0 . Таким образом, все фазовые траектории при → +∞ имеют асимптотой фазовые траектории I. Аналогично можно показать, что при → −∞ фазовые траектории имеет асимптотой траектории II (отвечающие 2 < 0).

Для остальных фазовых траекторий (С1С2 ≠ 0) можно записать два соотношения

dx

 

 

 

1

 

 

2

e

(

)t

 

 

dx

 

1

 

С h

С h

2

1

 

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

1

1

2

2

2

2

 

 

 

,

lim

2

 

2

 

 

 

1

 

 

2

 

(

)t

 

1

dx

 

 

 

 

 

e

 

t dx

 

 

С h

С h

2

1

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1

1

1

2

2

1

 

 

 

 

 

1

 

1

dx

 

С h1e( 1 2 )t С h2

 

dx

2

 

h2

2

 

1

1

2

 

2

2

2

, lim

 

 

2

 

 

 

1

( )t

 

 

 

2

 

 

 

2

dx

 

 

 

С

h

t dx

 

h

 

С h e

1 2

 

 

 

1

 

1

1

1

 

2

2

1

1

 

1

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

(17)

(18)

Положение равновесия седло.

III.

1 , 2 ,

1

= 2 =

≠ 0 .

Возможны два случая.

 

 

 

Случай III.A Существует только один линейно независимы собственный вектор, отвечающий собственному числу

В этом случае жорданова форма матрицы систему представляет собой жорданову клетку размера 2x2, а жорданов базис

 

 

h

 

 

 

 

h

 

 

 

 

 

 

1

 

 

 

 

 

 

2

 

 

 

1

 

1

 

 

h

2

 

1

 

 

состоит из собственного вектора

h

1

и присоединенного

 

h

2

. Общее решение системы (11) имеет вид:

 

 

 

h

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2

 

 

x (t)

 

 

 

 

 

 

 

 

 

 

 

 

h

 

 

 

 

h

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

1

 

 

 

 

2

 

 

 

 

1

 

 

1

1

 

t

 

2

1

 

2

 

t

1

1

 

 

t

 

2

 

1

 

 

 

1

 

t

X (t)

 

 

 

e

 

С

(h

t h

 

)e

 

1

e

 

С

1

t

 

2

e

 

 

 

 

С h

 

 

 

 

С

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

x

(t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h

 

 

 

 

 

 

h

 

 

 

h

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

2

 

2

 

 

Рассмотрим < 0 . Тогда положение равновесия = O асимптотически устойчиво.

(19)

Положим в (19)

С2

0

. Тогда

x

 

1

t

,

x

 

1

t

С h e

 

С h e

 

1

1

1

 

 

2

1

2

 

– параметрическое уравнение двух лучей, коллинеарных

собственному вектору

h

1

 

 

 

. Эти лучи «входят» в точку покоя = O, точка покоя асимптотически устойчива. Называем эти

лучи фазовыми траекториями I.

Рассмотрим теперь в (19)

С

0

2

 

. Тогда получаем соотношение

dx

 

С h1

С (th1

h2 ) С h1

 

С h1

/ t С (h1

h2

/ t) С h1

/ t

 

h1

2

 

1

2

2

2

2

2

2

 

1

2

2

2

2

2

2

 

 

2

dx

С h1

С (th1

h2 ) С h1

С h1

/ t С (h1

h2

/ t) С h1

/ t

h1

 

 

t

1

 

1

1

2

1

1

2

1

 

1

1

2

1

1

2

1

 

 

1

Фазовые траектории системы (11), отличные от фазовых траекторий I «входят» ( < 0)в точку покоя = O, касаясь фазовых траекторий I.

В другом случае, когда > 0 , ситуация аналогична, только положение равновесия = O неустойчиво, и фазовые траектории «выходят» из точки равновесия.

Положение равновесия = O называется вырожденным узлом (устойчивым при < 0 и неустойчивым при > 0 )

Таким образом имеем одинаковые углы наклона к которым стремятся фазовых кривых в точке покоя и на бесконечности:

 

 

 

 

 

x

 

 

dx

 

 

h1

 

1) С

 

 

0

2

 

2

 

 

 

2

, t

2

 

 

 

 

 

 

 

 

 

x

 

 

dx

 

 

h1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1

 

 

1

 

 

 

 

 

 

 

 

dx

 

 

 

h1

 

2) С

 

0

 

 

2

 

 

 

2

 

 

2

 

dx

 

h1

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Случай III.Б Существует два линейно независимых собственных вектора, отвечающий собственному числу

В этом случае все ненулевые векторы – собственные, матрица системы

 

 

 

 

 

 

 

 

x

x

 

1

1

,

 

 

 

x

x

 

 

2

2

 

A

E

, система распадается на два уравнения

Общее решение системы записывается в виде

x

С e

t

,

x

С e

t

и все решения представляют собой лучи, которые

 

 

1

1

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

«выходят» из начала координат (если

 

>

0

) и «входят» в начало координат,

если

 

<

0

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Положение равновесия называют особым (или дикритическим) узлом.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

IV. 1

, 2 ,

1 = 0,

2 ≠ 0 .

Общее решение системы (11) имеет вид

 

 

 

 

 

 

 

 

x (t)

 

 

 

 

 

 

 

 

h

 

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

2

 

 

 

 

 

 

 

X (t)

 

 

1

 

 

 

 

 

1

С

h

2

e

t

С

 

1

 

С

 

 

1

 

e

t

 

 

 

 

 

 

 

 

 

 

 

С h

 

2

 

1

2

 

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

x

(t)

 

 

 

 

 

 

 

 

 

 

h

 

 

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

h

 

 

 

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1

 

 

h

2

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где собственные векторы

h

h1

и

 

h2

соответствуют собственным числам

 

1,

 

2 соответственно.

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(20)

Положив в (20)

С

0

, имеем точки покоя, лежащие на прямой, проходящей через начало координат параллельно

h1

2

 

.

 

С

0

 

 

h

2

 

Если в (20)

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

.

 

2

 

 

 

 

 

 

 

 

 

 

 

 

При этом, если 2 < 0 , то все точки покоя являются устойчивыми по Ляпунову (но не асимптотически устойчивы), а если2 > 0 , то все точки покоя неустойчивы.