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

Конспект лекций по математике для студентов инженерно-технических специальностей. В 4 ч. Ч. 3

.pdf
Скачиваний:
1
Добавлен:
24.11.2025
Размер:
2.02 Mб
Скачать

ξ2

ξ1

0

Рис. 15.9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx1

 

= ax ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

Случай 2а.

Линейная

 

 

система

имеет

вид:

dt

 

 

 

 

 

 

корни

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx2

 

= ax2 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

характеристического

уравнения:

 

a − λ

 

0

 

= 0 ,

λ1 = λ2

= a .

 

Из уравнений

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

a − λ

 

 

 

 

 

 

 

самой

 

системы

она каноническая

– находим фазовые траектории:

 

dx2

=

x2

ln x

2

= ln x

+ ln c x

2

= cx

прямые,

заканчивающиеся или

 

 

 

 

dx1

x1

 

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

dx1

 

= x ,

 

1

Пример 1 5 . 5 3 . Линейная система: dt

 

ее характеристическое

dx2

 

= x2 ,

 

 

 

 

 

dt

 

 

уравнение:

 

1 − λ

0

 

= 0 λ1 = λ2 =1,

x1

= c1et

 

 

 

0

1 − λ

 

 

. Фазовые

 

 

 

 

x2

= c2 et x2 = cx1

траектории изображены на рисунке 15.10, особая точка – неустойчивый узел.

143

x2

x1

0

 

Рис. 15.10

 

 

 

 

 

 

Случай 3-й. Корни

характеристического

уравнения

комплексные с

ненулевой действительной

частью: λ1 = α + iβ,

λ2

= α − iβ.

В этом случае

каноническая система, полученная с помощью замены

ξ1 = Ax1

+ Bx2 ,

имеет

 

 

 

 

 

 

 

ξ2 = Ax1

Bx2 = ξ1 ,

 

 

dξ

 

 

 

 

 

d

 

 

 

 

 

 

 

 

 

 

 

ξ1 = u + iv,

 

 

 

 

 

 

ξ

 

 

 

 

 

 

 

 

 

 

вид

 

= λ1ξ1 ,

1

=

λ2 ξ1 . В канонической системе сделаем замену

 

1

 

 

 

 

 

 

 

dt

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

ξ1 = u iv,

получим уравнения:

u + iv = λ(u + iv),

и приравняем действительные и мнимые

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u iv =

λ(u iv)

 

 

 

 

 

 

 

du

= αu −βv,

 

 

 

 

 

 

 

 

 

 

части:

dt

Теперь

перейдем к полярным координатам на плоскости

 

 

 

 

 

 

 

 

 

 

dv

 

u + αv.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u0v :

 

 

 

 

 

u = r cosϕ,

v = r sin ϕ

получим

 

 

систему:

dr cosϕ −

dt

dr sin ϕ +

dt

при sin ϕ:

r sin ϕdϕ

= αr cosϕ −βr sin ϕ,

 

 

 

 

 

 

 

 

 

 

 

dt

Приравняем коэффициенты при cosϕ и

 

 

 

dϕ

 

r cosϕ

 

= αr sin ϕ + βr cosϕ.

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dr

= αr,

 

 

 

 

 

 

 

 

 

dt

Исключив t из этой системы, получим

dr

=

α

r ,

dr

=

α

dϕ,

 

 

 

 

dϕ

β

r

β

dϕ

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

144

< 0 , следовательно, x1 (t) возрастающая функция времени, а x2 (t) −

 

α

α

ϕ . Это

ln r =

ϕ + ln c . Следовательно, уравнения фазовых траекторий: r = ce β

 

β

 

 

логарифмические спирали. Точки покоя такого типа называются фокусы. При α > 0 с возрастанием ϕ возрастает r , т. е. изображающая точка удаляется по

спиралям от начала координат, значит, состояние равновесия неустойчиво. При α < 0 с возрастанием ϕ убывает r , изображающая точка по фазовым траекториям асимптотически приближается к началу координат, следовательно, устойчивый фокус и устойчиво равновесное состояние системы.

 

 

4-й случай. Корни характеристического уравнения чисто мнимые (α = 0)

λ1,2

= ±iβ. Уравнения фазовых траекторий в

полярных

координатах:

dr

= 0,

dϕ

= β r = r = const, ϕ = βt + ϕ

 

т. е.

фазовые

траектории

 

dt

0

dt

 

0

 

 

 

 

 

 

 

 

 

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

соответственно равновесие является устойчивым.

 

 

 

Пример 1 5 . 5 4 .

Исследовать

устойчивость состояния

равновесия

системы: x1 = 2x1 + 5x2 ,

x2 = −2x1 4x2 .

Корни характеристического

уравнения

 

2 λ

5

 

= 0 − комплексные: λ1,2

= −1 ± i , с отрицательной действительной

 

 

 

2

4 λ

 

 

 

 

 

частью, следовательно, начало координат – точка покоя типа фокус, фазовые траектории – логарифмические спирали: r = ce−ϕ . Т. к. действительная часть корней α = −1 < 0 , то состояние равновесия системы асимптотически устойчиво. Направление движения изображающей точки по спирали – к началу координат – это легко проверить с помощью уравнений исходной системы, т. к. в 1-й четверти dxdt1 > 0 и dxdt2

убывающая. См. рис. 15.10.

145

x1

x2

0

Рис. 15.11

 

 

Пример 1 5 . 5 5 . Определить тип точки покоя системы

x1 = x2 ,

 

= −4x1.

 

x2

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

 

λ

1

 

имеет чисто мнимые корни:

 

= 0

 

 

4

λ

 

 

λ1,2 = ±2i .

dx1

dx2

Приведем три способа решения задачи: 1-й способ (см. решение примера 15.48).

Найдем фазовые траектории, исключив t из уравнений системы:

= x42x1 , 4x1dx1 = x2 dx2 , 2x12 = 12 x22 2c 4x12 + x22 = c − это эллипсы, значит,

начало координат – особая точка типа центр.

2-й способ. Решим систему уравнений, тогда общее ее решение представит параметрические уравнения фазовых траекторий в окрестности точки покоя. Формула общего решения в случае комплексных корней характеристического уравнения имеет вид:

X (t) = C

Re(Veλt ) + C

I

m

(Veλt ) .

(15.143)

1

2

 

 

 

Координаты собственного вектора V, соответствующего собственному значению λ = 2i , найдем из системы уравнений:

146

(0 − 2i1 + 1 α2 = 0,− 4α1 + (0 − 2i2 = 0,

имеющей кроме тривиального бесчисленное множество решений, т. к. ∆ = 0 ,

α

 

= 2iα

. Например, при α

 

= 1,

 

 

α

 

 

 

 

 

 

 

 

 

1

 

2

1

 

 

2

 

= 2i V =

 

 

.

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ve

2it

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

(cos 2t + isin 2t)

 

 

 

 

 

 

 

 

 

 

2i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Re(Ve

2it

) =

 

 

cos 2t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

− 2sin 2t

 

 

 

 

 

 

 

I

 

 

(Ve

2it

 

 

 

sin 2t

 

 

 

 

 

 

 

 

m

 

 

) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2cos 2t

 

 

 

 

X (t) =

C

 

 

cos 2t

 

 

+ C

 

sin 2t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

− 2sin 2t

 

 

 

2cos 2t

Следовательно, параметрические уравнения фазовых траекторий

x1 (t) = C1 cos 2t + C2 sin 2t,

x2 (t) = −2C1 sin 2t + 2C2 cos 2t.

Исключим из них t :

2

x2

2

2

2

2

2

 

2

 

x1

+

 

 

= C1

+ C2

4x1

+ x2

= C

 

.

2

 

 

 

 

 

 

 

 

 

 

 

147

3-й способ. Система представляет 4-й случай для определения типа особой точки, т. к. характеристические числа λ1,2 = ±2i чисто мнимые, следовательно, в

полярных координатах для косоугольной системы фазовые траектории имеют уравнения r = r0 , ϕ = 2t + ϕ0 это окружности, особая точка – начало координат – центр. На рис. 15.12 дан фазовый портрет системы

v

x2

 

u

x1

0

Рис. 15.12

Замечание 15.14. Если определитель матрицы системы (15.136) ∆ = 0 , λ1 = 0 , λ2 = a + d , данная линейная система имеет бесчисленное множество точек покоя, при этом:

а) если λ2 < 0 , фазовые траектории образуют семейство параллельных прямых и все точки покоя, лежащие на них, устойчивы (рис. 15.13);

x2

x1

0

Рис. 15.13

б) если λ2 > 0 , начало координат – неустойчивая точка покоя;

в) λ1 = λ2 = 0 , матрица система нулевая, все точки плоскости устойчивые точки покоя;

148

г) λ1 = λ2 = 0 , d = −a , bc ≤0 , начало координат – неустойчивая точка покоя.

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

 

dx1

 

= a

x

+ a

x

+... + a

x

+ X (t),

 

dt

 

11

1

12

2

1n

n

1

 

 

 

 

 

 

 

 

 

 

 

 

dx2

 

= a21 x1

+ a22 x2

+... + a2n xn

+ X 2 (t),

(15.144)

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

........................................................

 

где X1 (t), X 2 (t),..., X n (t) − функции нелинейные, разложения которых по степеням начинаются с членов 2-го порядка; эти функции отбрасывают и исследуют линеаризованную систему

x1 = a11 x1

+ a12 x2

+... + a1n xn

,

(15.145)

 

= an1 x1

+... + ann xn .

 

 

xn

 

 

 

Аналогично, при исследовании устойчивости решений автономных

нелинейных систем более общего вида:

 

 

 

 

 

 

 

dx1

=

f (x , x ,..., x ),

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

1

1

2

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx2

 

 

 

f2

(x1 , x2

,..., xn ),

 

 

 

 

 

 

=

 

(15.146)

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

.................................

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dxn

 

=

f

(x , x ,..., x ),

 

 

 

 

 

 

 

 

dt

 

 

 

n

1

2

 

n

 

 

 

 

 

 

 

 

 

 

149

используют метод первого приближения, заменяя правые части уравнений с помощью формулы Тейлора в окрестности точки покоя (x1 0, x2 0,...)

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

dx1

 

= ( f1 )x1 x1

+ ( f1 )x 2 x2

+ ... + ( f1 )xn

xn ,

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(15.147)

........................................................

 

 

 

 

 

 

 

 

 

dx

n

 

 

 

 

 

 

 

 

 

= ( f )x

+ ( f )x

+ ... + ( f )

x .

 

 

 

dt

 

n x1 1

n x2 2

n xn n

 

 

 

 

 

 

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

Теорема 15.25.

Если

действительные

части

всех

корней

λS = αS + iβS (s =1,2,...,n) характеристического

уравнения

системы

(15.147)

отрицательны, то невозмущенное движение исходной системы (15.146) устойчиво асимптотически, независимо от членов выше 1-го порядка.

Теорема 15.26. Если среди корней характеристического уравнения имеется по крайней мере один с положительной вещественной частью, то невозмущенное

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

 

Пример 1 5 . 5 6 .

Исследовать

методом

 

первого

приближения

 

 

 

3

x +

1

sin 2 y,

 

 

 

x = −

2

2

 

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

 

 

 

 

 

 

 

2x.

 

 

 

 

y = −y

 

 

 

 

 

3

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2 y)3

x ≈ −

 

x + y

 

λ

1

 

5

 

31

 

 

 

 

 

 

 

sin 2 y = 2 y

 

+ ...

2

 

;

 

2

 

 

= 0 , λ1,2 = −

 

±

 

i .

3!

 

 

 

1λ

4

4

 

 

 

 

 

 

2

 

 

 

 

 

y ≈ −y

2x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

150

Т. к. α = −54 , то по теореме 15.25 Ляпунова состояние равновесия заданной

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

Замечание 15.15. В “критическом случае”, если все αi < 0 за исключением хотя бы одного α =0 , надо применять другие методы исследования.

Пример 1 5 . 5 7 .

 

 

 

 

 

 

 

 

x2

+

 

 

 

2 y +

(2 y)2

 

 

(3x)2

+

(3x)4

 

x = ex+2 y − cos3x

= 1 + x +

2!

... 1 +

2!

 

...

− 1

2!

 

4!

+... ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

1

1

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2x)

 

 

 

 

 

−1 (2x)

 

 

 

 

 

 

 

y

2

 

 

 

 

 

 

 

 

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y = 4

+ 8x − 2e

y

=

2 1

 

2

 

 

+

 

 

 

 

 

...

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

− 2 1

+ y +

 

 

 

+... .

 

 

 

 

 

 

 

 

 

2!

 

 

2!

 

 

 

 

 

 

 

 

 

 

1!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x x + 2 y

,

 

1 − λ

2

 

=0

,

λ1 = 2 , λ2 = −3.

 

 

 

 

2

− 2 − λ

 

y ≈ 2x − 2 y

 

 

 

 

 

 

Т. к. λ1 >0 , по теореме 15.26 Ляпунова состояние равновесия заданной и линеаризованной систем неустойчиво.

Критерии отрицательности вещественных частей корней характеристического уравнения

Из теорем Ляпунова следует, что для суждения об устойчивости положения равновесия линейной системы или движения нелинейных систем по теории первого приближения нужно знать, является ли устойчивым полином – развернутый определитель характеристического уравнения. Устойчивым называется полином, все корни которого αk + i βk лежат слева от мнимой оси, т. е.

все αk <0 . Существуют специальные критерии, позволяющие без нахождения корней полинома, что особенно важно для полиномов степени ≥ 2 , определить:

151

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

Критерий Рауса-Гурвица

Для устойчивости полинома a

zn + a zn1

+... + a

n

необходимо и достаточно,

 

 

 

 

 

 

 

 

 

0

1

 

 

 

 

 

 

 

 

 

 

 

 

чтобы при a0 > 0

определитель Гурвица ∆n

и все его главные диагональные

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a1

a3

 

 

 

 

a1

a3

a5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

= a >

0 ,

2

=

> 0 ,

3

=

a

0

a

2

a

4

> 0 ,

 

 

1

 

 

 

 

a0

a2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

a1

a3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a1

a3

 

a5

a7

 

 

 

 

a1

 

a3

 

 

a5 .....0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a0

 

a2

 

 

a4

.....0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

=

a0

a2

 

a4

a6

 

> 0,...n =

0

 

a1

 

 

a3

.....0

 

 

> 0 .

 

 

 

 

0

a1

 

a3

a5

 

 

 

 

0

 

a0

 

 

a2 .....0

 

 

 

 

 

 

 

 

0

a0

 

a2

a4

 

 

 

 

........................

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

0

 

 

 

0.....an

 

 

 

 

Критерий Льенара-Шипара

Для уравнения

a

k n + a

n1

k n1 +... + a k + a

0

= 0 с положительными

 

n

 

1

 

коэффициентами все корни имеют отрицательную действительную часть(αk < 0),

если положительны все i с четными индексами, т. е. 2 ,4 ,... или положительны все i с нечетными индексами.

Пример 1 5 . 5 8 . Проверить устойчивость некоторой системы,

характеристический полином которой 1,25λ3 + 50,6λ2 + 675λ + 3000 .

Вычислим определители Гурвица

152

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]