Конспект лекций по математике для студентов инженерно-технических специальностей. В 4 ч. Ч. 3
.pdf
ξ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
|
α |
α |
ϕ . Это |
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 − 2i)α1 + 1 α2 = 0,− 4α1 + (0 − 2i)α2 = 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 zn−1 |
+... + 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 |
n−1 |
k n−1 +... + 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
