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

высшая математика

.pdf
Скачиваний:
73
Добавлен:
18.02.2016
Размер:
5.88 Mб
Скачать

 

 

d(x + λy)

= (a1 + λa2 )(x + λy) + f1(t) + λf2 (t) ,

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

интегрируя которое способом, изложенным в п. 30 Л. 48, получаем

 

 

x + λy =

e( a1 + λa2 )t {C +

[ f1(t) +

λf2 (t)]e(a1 + λa2 )t dt}.

 

(8)

 

Если уравнение (7) имеет различные вещественные корни λ

1

и λ2,

то из (8) получим два независимых интеграла системы (1), и, значит,

интегрирование системы (1) будет окончено.

 

 

 

Пример 3. Решить методом Даламбера систему:

 

 

 

 

 

 

 

dx

=

5x +

4y +

 

e

t

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

(9)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy

=

4x +

5y +

1.

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ1, 2

Решение.

Выберем λ

из

(7): 4 + 5λ = λ (5 + 4λ ), откуда

= ± 1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда по формуле (8) для случая λ = 1 будем иметь:

 

 

 

x + y = e9t {C1 + (et + 1) e9t dt}= e9t {C1 + (e8t + e9t ) dt}=

 

 

 

 

 

 

= C1e9t

1

et

 

1

.

 

(10)

 

 

 

 

 

 

 

 

 

 

 

 

9

 

 

 

 

 

 

 

 

 

 

8

 

 

 

 

 

 

Для λ = – 1,

аналогично,

получаем

 

 

 

 

 

 

x y = et {C2 + (et 1) et dt}= C2 et + tet + 1 .

 

(11)

Итак, имеем два независимых интеграла системы (9). Решая систему (10), (11) относительно х и у, получим общее решение системы (9):

 

 

 

C1

 

9t

 

 

 

C2

 

 

t

 

 

1

t

 

 

1

 

 

4

 

x =

 

 

 

e

 

+

 

 

 

 

e

 

+

 

e

t

 

 

+

 

.

2

 

2

 

 

2

8

9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C1

 

9t

 

 

C2

 

t

 

 

1

 

t

 

1

 

 

5

 

 

y =

 

 

 

e

 

 

 

 

 

e

 

 

 

e

t

+

 

 

 

 

.

 

 

2

 

 

 

2

 

2

8

9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Замечание 1. Если правые части нормальной системы ДУ имеют вид

ax + by + P(t) , где а, b – постоянные, t>0, а Р(t) – многочлен от t, то подста- t

новка t = eτ приводит эту систему к системе с постоянными коэффициентами.

3 3 1

Пример 4. Решить систему ДУ:

 

 

 

 

t

dx

= − 2x + 2y + t,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(12)

 

 

 

 

dy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

x 5y

+

 

2

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

t

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение. Сделаем замену переменного t = eτ. Тогда:

 

dx

dx

 

dτ

1 dx

 

dy

 

1 dy

 

 

=

 

 

 

 

 

 

=

 

 

 

 

,

 

 

 

 

=

 

 

 

.

 

dt

dτ

 

dt

 

t

dτ

 

dt

 

 

t

dτ

Система (12) принимает вид:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

τ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= −

2x + 2y

+

e ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dτ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy

 

 

=

 

 

x 5y +

 

e

2τ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

dτ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для решения последней системы применим метод Даламбера. Умножим второе уравнение системы на λ и сложим почленно с первым:

 

 

 

d

(x + λy)

=

(

2 − λ)x + (2 5λ)y + eτ + λe2τ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

dτ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

 

 

 

 

 

 

 

 

2 5λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x + λy) =

 

 

 

τ

 

 

2τ

 

 

 

 

 

 

 

 

 

 

(

2

− λ) x +

 

y + e

 

 

+ λe

 

.

 

 

 

(13)

 

 

 

 

2 − λ

 

 

 

 

 

 

 

 

dτ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Выберем λ так, чтобы коэффициент при у в квадратной скобке

был равен λ , т.е.

 

2

5λ

 

= λ или λ2 – 3λ + 2 = 0, откуда λ

 

= 1,

λ

 

= 2.

 

2

λ

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При λ1 = 1 из (13) получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d(x + y)

= − 3(x + y) + eτ + e2τ ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

откуда, согласно формуле (8), имеем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x + y = e3τ[C1 + (eτ + e2τ) e3τdτ]= C1e3τ +

 

1

eτ + 1 e2τ.

 

 

(14)

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

При λ2 = 2 из (13) с использованием формулы (8) находим

 

 

 

 

 

 

 

 

 

x + 2y =

C2 e4τ +

1 eτ +

1 e2τ .

 

 

 

 

 

 

 

 

 

(15)

 

 

 

 

 

 

 

 

 

 

 

 

5

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решая систему (14)–(15) относительно х и у, получаем:

 

 

 

 

 

 

 

 

x =

 

2C1e3τ C2 e4τ + 0,3eτ +

1

e2τ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

15

 

 

 

 

 

 

 

 

 

 

 

 

3 3 2

y = − C1e3τ +

 

C2 e4τ − 0,05eτ +

2

e2τ.

 

 

 

 

 

 

 

 

 

 

 

15

 

Возвращаясь к переменному τ (eτ= t), получим общее решение

системы (12):

 

2C

 

 

 

C2

 

 

 

1

 

 

 

x =

 

1

 

+ 0,3t +

 

 

t 2 ,

 

 

3

 

4

15

 

 

 

t

 

 

 

t

 

 

 

 

y = −

C1

+

 

C2

0,05t +

 

2

t 2 .

 

3

4

15

 

 

t

 

 

t

 

 

 

 

 

! Задания для самостоятельной работы

1. Методом исключения решить следующие системы:

 

dx

=

3x + 8y,

 

dx

= − 9y,

 

dx

=

y + t,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а) dt

 

 

б) dt

 

в) dt

 

 

 

dy

= − x − 3y;

 

dy

= x;

 

dy

= x t;

 

 

 

 

dt

 

 

 

dt

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

dx

+

3x +

4y =

0,

 

 

 

 

 

г) dt

 

 

 

 

dy

+

2x +

5y =

0;

 

dt

 

 

 

 

 

4

dx

dy

+ 3x = sin t,

 

 

 

 

 

 

 

dt

д) dt

 

dx

+

y = cos t.

 

dt

 

 

 

 

 

2. Решить методом Даламбера следующие системы уравнений:

 

dx

=

5x + 4y,

 

 

dx

= 6x + y,

 

 

dx

=

2x − 4y + 1,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а) dt

 

 

 

б) dt

 

 

 

 

в) dt

 

 

 

dy

= x + 2y;

 

 

dy

 

 

= 4x + 3y;

 

 

dy

= − x + 5y;

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

dt

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

=

3x + y +

t

 

 

dx

=

2x +

4y +

 

 

 

 

 

 

 

 

 

 

e

,

 

 

 

cos t,

 

 

 

 

 

 

 

 

 

г) dt

 

 

 

д) dt

 

 

 

 

 

 

 

 

 

 

dy

=

x + 3y

et ;

 

dy

=

x

2y +

sin t.

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 3 3

Лекция 54

Устойчивость нулевого решения динамических систем второго порядка

Дано определение устойчивости по Ляпунову, приведены условия устойчивости линейных систем с постоянными коэффициентами, исследована устойчивость нелинейных систем по первому приближению.

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

 

dx

=

f (x, y),

 

 

 

 

 

dt

 

(1)

 

 

 

dy

=

g(x, y),

 

dt

 

 

 

где f (x, y) и g(x, y) – непрерывно дифференцируемые в некоторой области D плоскости Оxy (или во всей плоскости) функции. Координатную плоскость Оxy называют фазовой плоскостью.

Точками покоя (или положениями равновесия) динамической системы (1) называют такие точки (x; y), что выполняются соотноше-

ния: f (x, y) = 0, g(x, y) = 0.

Пусть g(0, 0) = f (0, 0) = 0, т.е. точка О(0; 0) (начало координат) является точкой покоя системы (1).

Будем говорить, что точка покоя х = у = 0 (или тривиальное решение x(t) y(t) 0) системы (1) устойчива по Ляпунову, если каково бы ни было ε> 0, можно найти такое δ= δ(ε) >0 , что для любого

решения col (x(t); y(t)), начальные данные которого x(0) = x0,

y(0) = y0

удовлетворяют условию

 

|x0| < δ, |y0| < δ,

(2)

выполняются неравенства

 

|x(t)| < ε, |y(t)| < ε для всех t 0.

(3)

Геометрически это означает следующее. Каким бы узким ни был цилиндр радиуса ε с осью Ot, в плоскости t = 0 найдется δ–окрест- ность точки О(0; 0; 0) такая, что все интегральные кривые x = x(t), y = y(t), выходящие из этой окрестности, для всех t 0 будут оставаться внутри этого цилиндра (рис. 1).

3 3 4

Если, кроме выполнения неравенств (3), выполняется также

условие lim

 

x(t)

 

=

lim

 

y(t)

 

= 0 , то

 

 

 

 

t → +∞

 

 

 

 

t → +∞

 

 

 

 

 

 

 

 

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

Точка покоя х = у = 0 неустойчива, если при сколь угодно малом δ> 0 хотя бы для одного решения col(x(t); y(t)) условие (3) не выполняется.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример

1.

 

Используя

 

 

 

 

Рис. 1

 

 

 

 

 

 

 

 

 

определение

устойчивости по

 

 

 

 

 

 

 

 

 

 

 

 

 

Ляпунову, показать, что решение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x(t) 0 , y(t)

0 системы

 

 

 

 

 

 

 

 

 

 

dx

=

y,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

(4)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy

=

x

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

устойчиво.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение. Любое решение системы (4), удовлетворяющее усло-

виям x(0) = x0, y(0) = y0,

имеет вид:

 

 

 

 

 

 

 

x(t) = x0cos t – y0sin t,

 

 

 

y(t) = x0sin t + y0cos t.

 

Возьмем произвольное

 

 

ε> 0 и покажем, что существует δ(ε) > 0,

такое, что при | x0| < δ,

| y0| < δ имеют место неравенства

 

 

 

x(t)

 

=

 

x0 cost y0 sin t

 

< ε,

 

 

 

y(t)

 

=

 

x0 sin t +

y0 cost

 

< ε

 

 

 

 

 

 

 

 

 

для всех t 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Это и будет означать, согласно определению, что нулевое решение x(t) 0, y(t) 0 системы (4) устойчиво по Ляпунову.

Имеем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x0 cost

y0

sin t

 

 

x0 cost

 

 

 

 

 

y0 sin t

 

 

 

 

x0

 

+

 

y0

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x0 sin t +

y0

cost

 

 

 

x0 sin t

 

 

+

 

 

y0 cost

 

 

 

x0

 

+

 

 

y0

 

 

(5)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

для всех t. Поэтому, если

 

x0

 

+

 

y0

 

< ε, то и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

для всех t .

 

x0 cost y0 sin t

 

 

 

< ε,

 

 

x0 sin t + y0 cost

 

< ε

(6)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 3 5

Следовательно, если, например, взять δ(ε) = 2ε , то при | x0| < δ, | y0| < δ

в силу (5) будут иметь место неравенства (6) для всех t 0, т.е. действительно нулевое решение системы (4) устойчиво по Ляпунову, но эта устойчивость неасимптотическая.

20. Устойчивость нулевого решения линейных систем с постоянными коэффициентами (ЛСПК). Рассмотрим частный случай системы (1) (ЛСПК):

 

dx

 

=

ax +

by,

 

 

 

 

 

 

 

dt

 

 

 

(7)

 

 

 

 

 

dy

=

cx +

dy,

 

dt

 

 

 

 

 

и пусть A =

a

b

;

det( λE

A) =

 

λ −

a

b

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

c

 

 

 

 

 

c

λ −

d

 

 

 

d

 

 

 

 

 

= λ2 (a + d)λ +

(ad

 

bc) = λ2

Sp A λ +

det A =

 

0

характеристическое уравнение матрицы А, а λ1, λ2 – корни этого уравнения.

 

Здесь Sp A – след матрицы A (сумма элементов по главной

диагонали), det A – определитель матрицы А.

 

Справедливы следующие утверждения:

 

1) для асимптотической устойчивости нулевого решения

x(t)

0, y(t) 0 системы (7) необходимо и достаточно, чтобы Re λ1 < 0 и

Re λ2

< 0;

 

2) в случае Re λ1 > 0 или Re λ2 > 0 нулевое решение системы (7)

неустойчиво;

 

3) в случае Re λ1 = Re λ2 = 0, Im λ1 = – Im λ2 0 нулевое решение

системы (7) устойчиво, но не асимптотически устойчиво;

4) в случае λ1

= λ2

= 0 нулевое решение системы может быть как

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

вание исходя из определения устойчивости).

Пример 2. Исследовать на устойчивость нулевое решение системы

 

 

 

dx

=

5x + 4y,

 

 

 

 

 

 

 

 

 

 

 

dt

 

(8)

 

 

 

 

 

 

 

dy

=

x + 2y.

 

 

 

dt

 

 

 

 

 

3 3 6

Решение. Имеем: Sp A = 5 + 2 = 7, det A = 10 – 4 = 6, характеристическое уравнение λ2 – 7λ + 6 = 0 λ 1 = 1, λ2 = 6, Re λ1 = 1 > 0. Значит, в силу 2) нулевое решение x(t) 0, y(t) 0 неустойчиво.

30. Устойчивость по первому приближению. Пусть имеем динамическую систему (1) с точкой покоя O(0; 0), где функции f (x, y) и g(x, y) дифференцируемы в окрестности начала координат достаточное число раз.

Разложим функции f и g по формуле Тейлора по х, у в окрестности начала координат:

f (x, y) = ax + by + R1(x, y), g(x, y) = cx + dy + R2(x, y),

где a = f (0,0) , b =

f (0,0)

, c =

g(0,0) , d =

g(0,0)

x

y

 

 

 

 

 

x

y

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

 

Тогда исходная система (1) примет вид:

 

 

 

dx

 

=

 

ax + by + R1(x, y),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy

=

 

cx + dy + R2 (x, y).

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

Вместо системы (1) рассмотрим систему

 

 

 

 

 

 

dx

=

ax +

by,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy

=

cx +

dy,

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

, а R1, R2 – члены

(1)

(7)

называемую системой уравнений первого приближения для системы (1).

Замечание 1. Если точка (x0; y0) – некоторое другое положение рав-

новесия системы (1) (x0

≠ 0

или y0

0) , то система первого приближения

строится так: в системе

(1)

сначала сделаем замену x = u + x

0

, y = v + y

0

и

получим функции

~

 

 

~

 

 

 

g(u, v) = g(u + x0, v

+ y0), f (u, v) = f (u + x0, v + y0),

a дальше

поступаем так же,

как и раньше с заменой х на u, а у на v.

 

 

 

 

 

Справедливы следующие утверждения:

 

 

 

 

 

1. Если все корни характеристического уравнения

 

 

 

 

 

 

λ2 Sp A λ +

det A = 0

 

 

 

(9)

имеют отрицательные вещественные части, то нулевое решение

х= у = 0, системы ( 7) и системы (1) асимптотически устойчиво.

2.Если хотя бы один корень характеристического уравнения (9) име-

3 3 7

ет положительную вещественную часть, то нулевое решение системы ( 7)

и системы (1) неустойчиво.

Говорят, что в случаях 1 и 2 возможно исследование на устойчивость по первому приближению.

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

Пример 3. Исследовать на устойчивость по первому приближе-

нию точку покоя х = 0, у = 0 системы:

 

 

 

 

 

 

x'

=

2x + y 5y 2

,

 

 

 

dx

 

 

dy

 

 

 

 

 

 

2

 

 

x

=

, y

=

(10)

 

 

 

 

 

 

 

 

 

 

.

'

=

3x +

y +

x

 

,

 

'

 

dt

'

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение. Система первого приближения

 

 

 

 

 

 

x =

 

2x +

y,

 

 

 

 

 

 

 

 

 

 

 

'

 

 

 

 

 

 

 

 

 

 

(11)

 

 

 

 

'

=

 

3x +

 

y,

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

нелинейные члены удовлетворяют нужным условиям; их порядок больше или равен двум. Составим характеристическое уравнение матрицы системы (11):

 

2

λ

1

 

= 0

или λ2

 

3λ − 1 = 0

λ1 =

3 +

 

13 , λ2

=

3

 

 

 

13

 

 

 

 

 

 

 

 

3

 

 

 

1− λ

 

 

2

2

вещественные

и

λ1 > 0.

Следовательно,

нулевое

решение

х = 0,

у = 0 системы (10)

неустойчиво.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

! Задания для самостоятельной работы

 

 

 

 

 

 

 

 

 

 

 

1.

Исследовать на устойчивость точку покоя О(0; 0)

 

 

систем:

 

 

dx

= 3x + y,

 

 

dx

=

x + 2y,

dx

 

=

x + 3y,

 

dx

 

 

= − 2x y,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

dt

 

 

dt

 

 

 

г)

dt

 

 

 

 

 

а)

dy

 

 

 

 

 

б)

 

 

 

 

 

в)

 

 

 

 

 

 

 

 

 

 

 

 

= − 2x + y;

 

 

dy

=

x + y;

 

dy

=

x + y;

 

 

dy

 

 

= 3x y.

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

dt

 

 

 

 

 

3 3 8

2. Исследовать на устойчивость по первому приближению нулевые решения х = 0, у = 0 следующих систем:

а)

в)

x'

=

x +

2y

sin y 2 ,

 

 

 

 

 

 

 

 

 

x 2

 

'

= − x − 3y

+ x

 

e

 

2

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x'

=

7x + 2 sin y

 

y 4

 

 

 

 

 

 

 

 

 

 

 

 

'

=

e

x

− 3y

− 1+

5

x

 

 

 

 

y

 

2

 

 

 

 

 

 

 

 

 

− 1 ;

,

2 ; г)

б)

x'

y'

 

'

=

− 2

+

 

sin

2

y,

x

 

 

 

= − x

− 3y

+ 4x 3 ;

y'

 

 

 

 

 

 

 

 

 

 

=

5

xe x

 

3y +

sin x 2 ,

2

 

 

 

 

 

 

 

 

 

 

 

 

 

ye

 

y 2

 

 

=

2x +

 

 

y 4 cos x.

 

2

3 3 9

Лекция 55

Фазовая плоскость

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

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

dx = f (x, y),

dt

(1)

dy = g(x, y),

dt

где f (x, y) и g (x, y) – непрерывно дифференцируемые в некоторой области D плоскости Oxy (или во всей плоскости) функции.

10. Случай линейной системы (1). Пусть соответствую-

щая (1) линейная система первого приближения (см. Л. 54) имеет вид

 

 

 

 

 

dx

 

=

ax +

by,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

(2)

 

 

 

 

 

 

 

 

 

 

 

 

 

dy

=

cx +

dy.

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

Обозначим

A =

a

b

 

 

 

 

 

 

 

 

 

. Если det A = ad cb ≠ 0, то система (2)

 

 

 

c

d

 

 

 

 

 

 

 

 

 

 

 

 

 

имеет единственное положение равновесия х = у = 0. Качественное поведение фазовых кривых (тип положения равновесия) определяется собственными числами матрицы коэффициентов этой системы, т.е. корнями уравнения

λ

2

– (a + d)λ + ad – cd = 0.

(3)

 

Если корни λ1 и λ2 уравнения (3) – действительные числа одного знака ( λ1λ2 > 0), то положение равновесия – узел, причем устойчивый, если λ1 < 0, λ2 < 0 (рис. 1), и неустойчивый, если λ1 > 0, λ2 > 0.

Если система (2) имеет вид dxdt = ax , dydt = ay , то узел называет-

3 4 0