- •Следствие 2. Всякое решение асимптотически устойчивой линейной системы (как однородной, так и неоднородной) асимптотически устойчиво в целом.
- •Приравнивая вещественные и мнимые части, получим эквивалентность
- •§ 4. Устойчивость линейных систем с постоянной матрицей
- •10. Критерий устойчивости.
- •Теорема 1. Линейная однородная система
- •Из выписанного представления благодаря соотношениям
- •является ненулевым решением системы (1), причём
- •Теперь покажем, что любое собственное значение с нулевой вещественной частью обладает тем свойством, что соответствующая ему клетка Жордана сводится к одному элементу.
- •Пусть матрица А приведена к жордановой форме:
- •Она является матричным решением системы (1), так как
- •Из (*) следует
- •Следовательно,
- •Поэтому неравенство (**) влечет
- •20. Критерий асимптотической устойчивости. Перейдем к решению вопроса асимптотической устойчивости линейной системы (1).
- •С помощью замены
- •перейдём к равносильной системе
- •влечёт
- •З а м е ч а н и е 2. Характеристическое уравнение для (4) имеет вид
- •В свою очередь, поведение руководителя y управляется руководителем второго ранга:
- •Соответствующее ему характеристическое уравнение имеет вид:
- •Корни этого уравнения образуют вершины правильного n-угольника на комплексной плоскости.
- •Упражнения
- •1) Исследовать на устойчивость следующие две системы и уравнение в зависимости от параметра α:
- •2) Убедиться в неустойчивости следующей линейной однородной системы четвертого порядка
- •§ 5. Полиномы Гурвица
- •О п р е д е л е н и е 1. Полином с вещественными коэффициентами
- •Теорема 1 (А. Стодола). Для того чтобы стандартный полином (1) являлся полиномом Гурвица, необходимо, чтобы все его коэффициенты были положительными.
- •□ Обозначим через
- •вещественные корни этого полинома.
- •По определению полинома Гурвица имеем
- •Из курса алгебры известно разложение на множители
- •Следующие две теоремы демонстрируют, что операция присоединения осуществляет тесную связь между стандартными полиномами Гурвица, степень которых отличается на одну единицу.
- •□ Введём семейство полиномов
- •Следовательно,
- •□ Имеем
- •Из первого равенства следует
- •Поэтому из второго равенства получаем
- •а значит
- •Пусть
- •Тогда
- •С учетом (*) имеем
- •Следствие 1. Для любого стандартного полинома
- •§ 6. Критерий Рауса-Гурвица
- •О п р е д е л е н и е 1. Матрицей Гурвица стандартного полинома
- •Благодаря сказанному имеет место
- •представляющим собой условия Гурвица, записанные для полинома g(z).
- •Рассмотрим линейную однородную систему
- •Напомним, что минор является главным, если номера занимаемых им строк совпадают с номерами занимаемых им столбцов.
- •Воспользовавшись условием положительности коэффициентов стандартного полинома как необходимым условием того, чтобы он был полиномом Гурвица, получим следующее утверждение.
- •Для асимптотической устойчивости линейной системы (2) необходимо, чтобы выполнялись неравенства
- •в частности,
- •Следствие 3. Для асимптотической устойчивости линейной системы (2) необходимо и достаточно выполнение следующих условий:
- •Пример 1. Пусть
- •Тогда
- •и условия Гурвица принимают вид
- •является асимптотически устойчивой.
- •Первое решение (на основе анализа корней характеристического полинома). Составляем характеристическое уравнение
- •Находим его корни
- •Второе решение (на основе теоремы 1 из § 5). Характеристический полином имеет вид
- •Третье решение (на основе критерия Рауса-Гурвица). Составляем матрицу Гурвица
- •Упражнения
- •Пользуясь критерием Рауса-Гурвица, исследовать на устойчивость следующие системы и уравнения
- •§ 7. Критерий Михайлова
- •□ Введем следующие обозначения
- •– его отрицательные вещественные корни;
- •Поэтому
- •Следовательно,
- •Последовательно находим:
- •Следовательно,
- •Следствие 1. Если для стандартного полинома степени n имеет место неравенство
- •то он не является полиномом Гурвица.
ГЛАВА 1. УСТОЙЧИВОСТЬ ЛИНЕЙНЫХ СИСТЕМ_______________________________________________49
A0 > 0, T1 = A1 > 0, …, Tn+1 > 0 .
Покажем, что F(z) Hn +1 . |
|
|
Согласно следствию 1 предыдущего параграфа, |
полином F(z) можно |
|
рассматривать как присоединенный к некоторому стандартному полиному |
||
f(z) = a0 +a1z +... +a n zn |
(an ≠ 0, |
a0 > 0) . |
Как и в доказательстве необходимости получаем, что главные диагональные
миноры Tk+1 матрицы Гурвица |
|
ΓF |
|
|
удовлетворяют равенствам |
|
||||||||||||||||||||||
|
|
|
|
|
T |
+1 |
= (2c)k+1 a |
0 |
|
|
|
k |
, |
k =1,2,...,n . |
|
|
||||||||||||
|
|
|
|
|
|
k |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
А так как согласно предположению |
|
|
|
Tk+1 > 0 , то k > 0 , |
k =1,2,...,n , т.е. для |
|||||||||||||||||||||||
полинома f(z) |
выполнены условия Гурвица (1). Отсюда в соответствии с ин- |
|||||||||||||||||||||||||||
дуктивным предположением получаем, что |
f(z) Hn . Таким образом, полином |
|||||||||||||||||||||||||||
F(z) |
является присоединенным к стандартному полиному Гурвица f(z) , а зна- |
|||||||||||||||||||||||||||
чит, в силу теоремы 2 из § 5, верно |
F(z) Hn+1 ■ |
|
|
|
||||||||||||||||||||||||
|
Если f(z) = a |
0 |
+a |
z +... + a |
n |
zn |
– стандартный полином Гурвица, то поли- |
|||||||||||||||||||||
ном |
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
f(z) = z |
n |
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
g |
, |
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
z |
|
|
|
|
|
|
где |
g(z) = a0zn +a1zn-1 +... +a n , |
также является стандартным полиномом Гур- |
||||||||||||||||||||||||||
вица, и обратно. |
|
|
|
z j , |
j =1,2,...,n , – |
|
|
|
|
|
|
|
|
|
и Rez j < 0 тогда и |
|||||||||||||
|
В самом деле, |
|
|
корни полинома |
f(z) |
|||||||||||||||||||||||
только тогда, когда |
|
1 |
|
– корни полинома |
g(z) , причем |
|
|
|
||||||||||||||||||||
|
z j |
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
Re |
1 |
= |
|
Rez j |
|
< 0 , |
|
|
j =1,2,...,n . |
|
|
|||||||||
|
|
|
|
|
|
|
|
z j |
|
|
z j |
|
|
2 |
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
Благодаря сказанному имеет место |
|
|
|
|
|
|
|
|
|||||||||||||||||||
вам |
Следствие 2. Условия Гурвица (1) теоремы 1 эквивалентны неравенст- |
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
~ |
0; |
~ |
1 = a n−1 > 0; |
|
~ |
= |
|
a n-1 |
|
a n |
|
> 0; … ; |
~ |
= a0 |
~ |
||||||||||||
|
|
|
|
|
||||||||||||||||||||||||
|
0 = a n > |
|
|
2 |
an-3 |
|
a n-2 |
|
n |
n-1 > 0, |
ГЛАВА 1. УСТОЙЧИВОСТЬ ЛИНЕЙНЫХ СИСТЕМ_______________________________________________50
представляющим собой условия Гурвица, записанные для полинома g(z).
Рассмотрим линейную однородную систему
|
dx |
= A x |
(t ≥ 0) |
(2) |
|
dt |
|||
|
|
|
|
|
с постоянной матрицей A = (aij ) . Пусть |
det(λE −A)= 0 – соответствующее ха- |
рактеристическое уравнение. В раскрытом виде это уравнение можно записать так:
λn − A1λn−1 + A2λn−2 +... + (−1)n An = 0 ,
где
A1 = ∑aii =SpA
|
i |
|
|
aii |
aij |
|
|
|
A2 = ∑det |
|
|
i<j |
a ji |
a jj |
……………………
Ak равно сумме всех главных миноров порядка k матрицы A
……………………
An = det(A) .
Напомним, что минор является главным, если номера занимаемых им строк совпадают с номерами занимаемых им столбцов.
Воспользовавшись условием положительности коэффициентов стандартного полинома как необходимым условием того, чтобы он был полиномом Гурвица, получим следующее утверждение.
Для асимптотической устойчивости линейной системы (2) необходимо, чтобы выполнялись неравенства
− А1 > 0; А2 > 0;...;(−1)n Аn > 0;
в частности,
SpA < 0, (−1)n det(A) > 0.
Нетрудно заметить, что для системы второго порядка (т.е. при n = 2 ) последние два условия являются и достаточными для асимптотической устойчивости.
ГЛАВА 1. УСТОЙЧИВОСТЬ ЛИНЕЙНЫХ СИСТЕМ_______________________________________________51
Следствие 3. Для асимптотической устойчивости линейной системы (2) необходимо и достаточно выполнение следующих условий:
ˆ |
|
= −A > 0; ˆ |
|
= |
|
− A1 |
1 |
|
= −A A |
+ A |
|
> 0; ... ; ˆ |
|
= (−1)n А |
|
~ |
|
> 0 . |
|
|
|
|
|
|
|
|
|||||||||||
|
1 |
1 |
2 |
|
|
− A3 |
A2 |
|
1 2 |
|
3 |
|
n |
|
n |
|
n-1 |
|
Пример 1. Пусть
f(z) =1+az + bz2 +az3 +z4 .
Тогда
|
|
a |
1 |
0 |
0 |
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
b |
|
|
a |
1 |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
a |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
Γf = |
0 |
1 |
|
|
a |
b |
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
0 |
0 |
0 |
1 |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
и условия Гурвица принимают вид |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
a |
|
1 |
0 |
|
= a 2 |
|
1 |
1 |
0 |
|
= a2 |
|
1 |
1 |
0 |
|
|
||
|
|
|
|
|
|
|
|
||||||||||||||
1 = a > 0; 2 = ab − a > 0; 3 = |
a b a |
|
|
1 b 1 |
|
|
0 b −1 1 |
|
= |
||||||||||||
|
|
0 |
|
1 |
a |
|
|
|
|
|
0 |
1 |
1 |
|
|
|
0 |
1 |
1 |
|
|
|
= a2 (b −2) > 0 , |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
что равносильно неравенствам a > 0, |
b > 2 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Пример 2. Выясним, при каких значениях параметров α |
|
и β |
система |
||||||||||||||||||
dx |
= −x +αy, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dt |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dy |
=βx − y +αz, |
|
|
|
|
|
|
|
|
|
( α,β R ). |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
dt |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dz |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= βy −z, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dt |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
является асимптотически устойчивой.
Первое решение (на основе анализа корней характеристического полинома). Составляем характеристическое уравнение
ГЛАВА 1. УСТОЙЧИВОСТЬ ЛИНЕЙНЫХ СИСТЕМ_______________________________________________52
|
λ +1 |
−α |
0 |
|
= (λ +1) [λ2 + 2λ +(1 −2αβ)]= 0 . |
|
|
||||
|
−β λ +1 −α |
|
|||
|
0 |
−β |
λ +1 |
|
|
Находим его корни
λ1 = −1, λ2,3 = −1 ± 2αβ .
Система асимптотически устойчива тогда и только тогда, когда все корни имеют отрицательные вещественные части. Следовательно, область асимптотической устойчивости в пространстве параметров определяется неравенством
αβ < 12 .
Второе решение (на основе теоремы 1 из § 5). Характеристический полином имеет вид
λ3 + 3λ2 +(3 −2αβ)λ +1 −2αβ.
Необходимое условие гурвицевости этого полинома состоит в положительно-
сти всех его коэффициентов, что вновь приводит к условию αβ < 12 . Здесь сле-
дует, однако, заметить, что в общем случае использование упомянутого необходимого условия (поскольку оно не обязательно является достаточным) может привести к более широкому множеству, чем область асимптотической устойчивости.
Третье решение (на основе критерия Рауса-Гурвица). Составляем матрицу Гурвица
3 |
−2αβ 1 −2αβ |
0 |
|
|
|
|
|
|
|
Γ = 1 |
3 |
3 − |
2αβ |
|
|
|
|
|
|
|
0 |
0 |
1 |
|
|
|
и выписываем условия |
Гурвица |
|
|
1 = 3 |
−2αβ > 0, 2 =8 −4αβ > 0 |
αβ< |
3. |
|
|
|
2 |
Полученный результат не совпадает с найденным ранее и дает более широкие границы. Почему? Дело в том, что полином, о котором идёт речь в критерии Рауса-Гурвица, должен быть стандартным. Это в данном случае означает, что