Введение в качественный анализ динамических систем. Учебное пособие
.pdf§2. ПРИВЕДЕНИЕ ДИНАМИЧЕСКОЙ СИСТЕМЫ
ККАНОНИЧЕСКОМУ ВИДУ
При рассмотрении характера простых состояний равновесия линейные члены в системе (а ) надлежащим образом выбраны неособым линейным преобразованием
u p11 p12 , |
|
v p21 p22 |
(3.1) |
|||
(т. е. преобразованием, у которого |
|
p11 |
p12 |
|
0 ), |
приводятся к воз- |
|
|
|||||
|
|
p21 |
p22 |
|
|
|
|
|
|
|
|||
можно более простому, так называемому «каноническому», виду. Посмотрим, прежде всего, при каких условиях надлежащим преобразованием (3.1) можно привести систему (а ) к виду (невыписанные члены содержат u и v в степени не ниже второй)
du |
u |
dv |
|
2 |
v |
(3.2) |
|
|
|||||
dt |
1 |
dt |
|
|
||
|
|
|
|
|||
Подставляя в (3.2) выражения (3.1) для u и v, а затем заменяя d / dt и d / dt через их выражения из (а ), мы получаем тождества. Прирав-
ниваявэтихтождествахкоэффициентыприлинейныхчленах, получаем следующие четыре линейные однородные уравнения относительно коэффициентов p11, p12 , p21, p22 искомого линейного преобразования:
p11(a 1) p12c 0, |
p21 |
(a 2 ) p22c 0, |
|
|
(3.3) |
p11b p12 (d 1) 0, |
p21b p22 (d 2 ) 0. |
|
Эти уравнения дают для p11, p12 , p21, p22 решения, не равные тождественно нулю только в том случае, когда 1 и 2 являются корнями уравнения
|
a |
b |
|
2 (a d) (ad bc) 0 , |
(3.4) |
|
|
||||
|
c |
d |
|
|
|
которое называется характеристическим. Корни 1 и 2 |
называются |
||||
характеристическими корнями состояния равновесия (особой точки).
31
Рассмотрим следующие случаи. |
|
|
|
|
|
|
|
|
|
||
1. Корни 1 и 2 |
действительны и различны. Тогда из уравне- |
||||||||||
ния (3.3) можно найти |
p |
, p , p |
21 |
, p |
22 |
, такие, что |
|
p11 |
p12 |
|
0 , |
|
|
||||||||||
|
11 |
12 |
|
|
|
p21 |
p22 |
|
|
||
|
|
|
|
|
|
|
|
|
|
||
и, следовательно, приведение динамической системы (а), см. главу 2, к виду (3.2) возможно.
2. Корни 1 и 2 – кратные, 1 2 .
В этом случае приведение к виду (3.2), вообще говоря, невозможно. Однако в этом случае можно указать неособые преобразования, с помощью которых система приводится к виду
du |
u , |
dv |
v u |
(3.5) |
|
dt |
dt |
||||
|
|
|
(В частных случаях может быть равно нулю. Этот вид называется каноническим в случае кратных корней.)
|
3. Корни 1 |
и 2 |
– комплексно-сопряженные: 1 i , |
2 |
i ( 0 , |
может как быть, так и не быть равным нулю). |
|
В этом случае при действительных и получим комплексные сопря-
женные u и v, так что приведение к виду (3.2) невозможно. В этом случае, вводя новые переменные u1 , v1 , u u1 iv1 , v u1 iv1 , нетрудно
установить, что система может быть приведена к следующему виду:
du1 |
u |
v |
, |
dv1 |
u |
v |
, |
(3.6) |
|
||||||||
|
|
|||||||
dt |
1 |
1 |
|
dt |
1 |
1 |
|
|
|
|
|
|
|
|
|
который в этом случае считается каноническим видом. Дополнительно отметим, что если рассмотреть линейную систему
d |
a b , |
d |
c d , |
|
dt |
dt |
|||
|
|
которая получается из системы (а ) отбрасыванием нелинейных членов, то, как известно, общее решение этой системы имеет вид
c1e 1t c2e 2t ;
c1 1e 1t c2 2e 2t .
32
Здесь 1 и 2 являются корнями того же характеристического уравне-
ния (3.4).
Характеристические корни не меняются при линейной замене координат (характеристические корни являются инвариантами линейного преобразования координат), т. е. пусть дана система
dx |
ax by , |
dy |
cx dy , |
|
dt |
dt |
|||
|
|
у которой характеристические корни 1 и 2 , и после линейного преобразования
x q11u q12v , |
y q21u q22v |
получаем систему в новых координатах u и v:
du |
Au Bv , |
dv |
Cu Dv . |
|
dt |
dt |
|||
|
|
Тогда характеристические корни последней системы, т. е. корни характеристического уравнения
A |
B |
|
0 , |
|
|||
C |
D |
|
|
равны 1 и 2 .
§ 3. ВОЗМОЖНЫЙ ХАРАКТЕР ПРОСТЫХ СОСТОЯНИЙ РАВНОВЕСИЯ
ГРУБЫЕ СОСТОЯНИЯ РАВНОВЕСИЯ
Сохраняем обозначения (см. § 1): |
|
|
|
|
||
|
|
Px (x0 , y0 ) |
Py (x0 , y0 ) |
|
. |
|
|
|
|||||
|
|
|
|
, y0 ) |
|
|
|
|
Qx (x0 , y0 ) Qy (x0 |
|
|
||
Характеристическое уравнение (3.4) может быть записано в виде2 0 . Возможны следующие случаи.
33
1. 0 , 2 4 0 , корни характеристического уравнения действительны и одинаковых знаков. В этом случае все траектории, проходящие через некоторую достаточно малую окрестность состояния равновесия О, стремятся к этому равновесию:
при t , когда 1 и 2 отрицательны,
при t , когда 1 и 2 положительны.
Состояние равновесия называется устойчивым узлом, когда 1 0 ,
2 0 , и неустойчивым узлом, когда 1 0 , 2 0 .
2. 0 , корни характеристического уравнения действительны и различных знаков: 1 2 0 . Состояние равновесия является седлом.
3. 0 , 2 4 0 , 0 , корни характеристического уравнения комплексные сопряженные: 1 i , 2 i , причем действи-
тельные части этих корней отличны от нуля.
В этом случае (так же, как и в случае 1) все траектории, проходящие через некоторую достаточно малую окрестность состояния равновесия O, стремятся к такому состоянию:
при t , когда 0 , при t , когда 0 .
Состояние равновесия называется фокусом и при этом устойчивым,
когда 0 , и неустойчивым, когда 0 .
§4. ЗАМЕЧАНИЯ О ГРУБЫХ СОСТОЯНИЯХ РАВНОВЕСИЯ
Пусть O(x0 , y0 ) – состояние равновесия системы (а), см. главу 2,
авыражения для и приведены в § 3, и
Px (x0 , y0 ) |
Py (x0 |
, y0 ) |
|
2 0 |
(3.7) |
||
|
|||||||
|
, y0 ) |
|
, |
y0 ) |
|
|
|
Qx (x0 |
Qy (x0 |
|
|
|
|||
– характеристическое уравнение этого состояния равновесия.
Для простого состояния равновесия по самому его определению0 , т. е. корниуравнения(3.7) – характеристическиекорни– отличны
34
от нуля. Уравнение, определяющее направления, по которым траектории стремятся к состоянию равновесия:
|
|
(x0 |
, y0 )k |
2 |
|
|
|
|
|
|
(x0 , y0 ) 0 . |
(3.8) |
|||
|
Py |
|
Qy |
Px |
k Qx |
||||||||||
Корни 1 |
и |
2 |
характеристического уравнения (3.7) и корни k1 |
||||||||||||
и k2 уравнения (3.8) связаны соотношениями: |
|
||||||||||||||
k |
|
|
1 Px (x0 , y0 ) |
, |
|
k |
2 |
|
2 Px (x0 , y0 ) |
. |
|
||||
1 |
|
Py (x0 , y0 ) |
|
|
|
|
Py (x0 , y0 ) |
|
|||||||
|
|
|
|
|
|
|
|
|
|||||||
Очевидно, что корни k1 и k2 действительны тогда и только тогда, когда действительны 1 и 2 .
В зависимости от того, каковы характеристическиекорни состояния равновесия, система может быть в окрестности этого состояния равновесия приведена линейным преобразованием переменных к одному из следующих видов, которыеназываются каноническими (обозначения переменных сохраняются прежними).
1. Характеристические корни действительны и различны ( 1 2 ). Канонический вид системы:
x 1x (x, y) , |
y 2 y (x, y) . |
|
|
2. Характеристические корни равны ( 1 2 ). Канонический вид системы:
x x (x, y) , |
y y x (x, y) |
|
|
( может быть как равным, так и не равным нулю).
3.Характеристические корни – комплексно-сопряженные ( 1 i ,
2 i , 0 ). Канонический вид системы:
x x y (x, y) , |
y x y (x, y) . |
|
|
Таким образом, характеристические корни указывают на состояние равновесия рассматриваемой системы.
35
Глава 4
КАЧЕСТВЕННАЯ СТРУКТУРА ОКРЕСТНОСТЕЙ НЕКОТОРЫХ СЛОЖНЫХ СОСТОЯНИЙ РАВНОВЕСИЯ
§1. НАПРАВЛЕНИЯ, В КОТОРЫХ ТРАЕКТОРИИ СТРЕМЯТСЯ
КСЛОЖНОМУ СОСТОЯНИЮ РАВНОВЕСИЯ
При исследовании сложных состояний равновесия иногда бывает весьма существенно знание направлений, в которых траектории могут стремиться к этому состоянию равновесия.
Рассмотрим динамическую систему
dx |
P(x, y) , |
dy |
Q(x, y) , |
|
dt |
dt |
|||
|
|
для которой начало координат является сложным состоянием равновесия, такчто
P(0, 0) = 0, |
Q(0, 0) = 0, |
|
||||
|
|
Px (0, 0) |
Py (0, 0) |
|
0 . |
(4.1) |
|
|
|||||
|
|
|
|
|
|
|
|
|
Qx (0, 0) |
Qy (0, 0) |
|
|
|
Предположим, что разложения правых частей в ряд по степеням x, y в окрестности точек O(0, 0) имеют вид
P(x, y) Pm (x, y) (x, y) ; Q(x, y) Qm (x, y) (x, y) ,
36
где m 1 , Pm (x, y) и Qm (x, y) – однородные многочлены, состоящие
из всех членов m-го порядка соответствующих разложений, а функции(x, y) и (x, y) – ряды, состоящиеизчленовболеевысокихпорядков.
При этом полагаем, что многочлены Pm (x, y) и Qm (x, y) одновре-
менно не равны тождественно нулю (в противном случае мы бы взяли m m ).
Рассмотрим выражение
xQm (x, y) yPm (x, y) , |
(4.2) |
а также выражение
Qm (1, k) kPm (1, k) ,
которое получим из предыдущего, если поделим его на xm 1, и введем обозначения
k xy .
Рассмотрим также выражение
cos Qm (cos , sin ) sin Pm (cos , sin ) .
Имеет место следующее утверждение: всякая полутраектория системы
dxdt Pm (x, y) (x, y) ; dydt Qm (x, y) (x, y)
(Pm и Qm неравны тождественнонулю), стремящаясяк состояниюравновесия O(0, 0), либо является спиралью, стремящейся к О при t
(или t ), либо стремится к О в определенном направлении * . При этом:
1) если хоть одна из траекторий системы является спиралью, стремящейся к О при t (или t ), то все траектории, проходящие
37
через точки некоторой окрестности состояния равновесия О, являются такими же спиралями (т. е. точка есть устойчивый или неустойчивый «фокус высшей сложности»);
2) если выражение (4.2) не обращается тождественно в нуль, то наклоны k, с которыми траектории стремятся к состоянию равновесия O, удовлетворяют уравнению
Q |
(1, k*) k*P |
(1, k*) 0 |
(4.3) |
m |
m |
|
|
или, иначе, направления * , скоторымитраектории стремятсяк О, удовлетворяют уравнению
cos *Qm (cos *, sin *) sin *Pm (cos *, sin *) 0 ; |
(4.4) |
3) если
xQm (x, y) yPm (x, y) 0
и, следовательно,
Pm (x, y) xQm 1(x, y) , Qm yQm 1(x, y) ,
где Qm 1(x, y) – некоторыйнеравныйнулютождественнооднородный
многочлен степени m – 1, то какое бы направление , не удовлетворяющее уравнению
Qm 1(cos *, sin *) 0 , |
(4.5) |
мы не взяли, существует в точности одна полутраектория, стремящаяся
к O в направлении . Для особого же направления * , удовлетворяющего (4.5), может оказаться, что не существует ни одной полутраекто-
рии, стремящейсяк O вэтом направлении * , либоестьконечноечисло такихтраекторий, либо, наконец, такихтраекторийможетсуществовать бесчисленное количество.
Заметим, чтоеслисуществуеттраектория, стремящаясяксостоянию
равновесия O и имеющая определенный наклон k* , то этот наклон является действительным корнем уравнения Qm (1, k) kPm (1, k) 0 .
38
Однако если это уравнение имеет действительные корни, то это еще не означает, что существуют траектории, стремящиеся к O с этим наклоном: возможны случаи, когда при этом все траектории являются спира-
лями, или замкнутыми траекториями.
§2. СЛОЖНОЕ СОСТОЯНИЕ РАВНОВЕСИЯ (ОСОБАЯ ТОЧКА)
СНУЛЕВЫМИ ХАРАКТЕРИСТИЧЕСКИМИ КОРНЯМИ
Рассмотрим систему
dx |
P(x, y), |
dy |
Q(x, y) , |
|
dt |
dt |
|||
|
|
где P(x, y) и Q(x, y) – аналитические функции, не имеющие общего множителя, отличного от постоянного.
Пусть начало координат является сложным состоянием равновесия этой системы, т. е. мы имеем
P(0, 0) = 0, |
Q(0, 0) = 0, |
|
|
||||
|
|
Px (0, 0) |
Py (0, 0) |
|
0 |
, |
(4.6) |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
Qx (0, 0) |
Qy (0, 0) |
|
|
|
|
и, следовательно, хотя бы один из характеристических корней этого состояния равновесия равен нулю. Рассмотрим такие сложные состояния равновесия, когда в разложениях по степеням x и y функций P(x, y) и Q(x, y) хотя бы один из линейных членов не равен нулю, т. е. когда
Px (0, 0) Py (0, 0) Qx (0, 0) Qy (0, 0) 0 .
Рассмотрим наряду с величиной (0, 0) величину
Px (0, 0) Qy (0, 0) .
Среди состояний равновесия, для которых выполняется условие (4.6), естественным образом выделяются два случая в зависимости от того, чтоимеетместо (0, 0) 0 или (0, 0) 0 . Таккакхарактеристическое
уравнение имеет вид
2 0 ,
39
то, очевидно, вслучае, когда 0 , толькоодинхарактеристическийкорень равен нулю, второй же равен . Если же 0 , то оба характеристических корня равны нулю.
1. (0, 0) 0 , (0, 0) 0 . В этом случае существует неособое ли-
нейное преобразование, с помощью которого система в окрестности начала O(0, 0) может быть представлена в следующем каноническом виде:
|
dx |
P* (x, y) , |
dy |
by Q*(x, y) , |
(4.7) |
|
dt |
dt |
|
||
где b 0 , а разложения по степеням x, y функций P*(x, y) |
и Q*(x, y) |
||||
начинаются с членов, не менее чем второго порядка. |
|
||||
Введем в рассмотрение функцию |
|
||||
y (x) ,
являющуюся решением уравнения
by Q*(x, y) 0 .
Подставим функцию y (x) в P*(x, y) и введем обозначение
(x) P*(x, (x)) .
Так как предположили, что функции P(x, y) и Q(x, y) не имеют общего множителя, отличного от постоянного, то функция P*(x, (x)) не
может быть тождественно равна нулю и, следовательно, в разложении функции (x) по степеням x заведомо будут отличные от нуля члены.
Таким образом, мы можем написать, что
(x) P*(x, (x)) m xm ,
где m 2 (так как разложение P*(x, y) по степеням x и y начинается с членов не ниже второй степени и m 0 ). Число m, очевидно, характеризует кратность общей точки кривых
P*(x, y) 0 |
и by Q*(x, y) 0 . |
40
