- •Следствие 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. УСТОЙЧИВОСТЬ ЛИНЕЙНЫХ СИСТЕМ_______________________________________________31
Следствие 2. Всякое решение асимптотически устойчивой линейной системы (как однородной, так и неоднородной) асимптотически устойчиво в целом.
30. Устойчивость линейной однородной системы с комплексной мат-
рицей. Рассмотрим линейную систему (1) в предположении, что элементами матрицы A(t) являются комплекснозначные функции переменной t ≥ 0. Выделим вещественную и мнимую части матрицы A: A(t) = A′(t) +i A′′(t) , где мат-
рицы A′, A′′ имеют вещественные элементы и i – мнимая единица.
Прежде всего, установим, как в рассматриваемом случае следует понимать решение системы (1). С этой целью подставим в равенство (1) комплекснозначную векторную функцию x = x(t)= u(t) +i v(t) :
ddut +i ddvt = (A′+i A′′) (u +i v) .
Приравнивая вещественные и мнимые части, получим эквивалентность
|
|
|
du(t) |
|
′ |
′′ |
||
dx(t) |
|
|
|
dt |
= A (t) u(t) −A |
(t) v(t), |
||
= A(t) x(t) |
|
|
|
|
|
|||
dt |
|
|
|
|
|
|
||
|
|
dv(t) |
|
′′ |
′ |
|
||
|
|
|
|
dt |
= A (t) u(t) + A (t) v(t), |
|||
|
|
|
|
|
|
|
|
|
которая и показывает, каким |
|
образом |
надлежит |
понимать решение |
||||
x(t) = u(t) +i v(t) |
системы (1) с комплекснозначной матрицей A(t). В частно- |
|||||||
сти, если матрица A(t) – вещественная, то получаем |
|
|||||||
|
|
|
|
|
du(t) |
= A(t) u(t), |
||
|
dx(t) |
|
|
|
|
dt |
|
|
|
= A(t) x(t) |
|
|
|
|
|||
|
dt |
|
|
|
|
|||
|
|
|
|
dv(t) |
= A(t) v(t). |
|||
|
|
|
|
|
|
dt |
||
|
|
|
|
|
|
|
|
Следует отметить, что большинство утверждений для линейных однородных систем с вещественной матрицей, известных из курса теории дифференциальных уравнений (по крайней мере, те утверждения, которые будут использоваться далее), сохраняют свою силу и для подобных систем с комплексной матрицей.
Для линейных систем с комплексной матрицей совершенно аналогично вводятся определения устойчивого (и асимптотически устойчивого) решения и системы. Внешне они ничем не отличаются от приведенных ранее определений и по этой причине здесь не приводятся. Анализ доказательств утверждений, сформулированных в данном и предыдущем параграфах, показывает, что они
ГЛАВА 1. УСТОЙЧИВОСТЬ ЛИНЕЙНЫХ СИСТЕМ_______________________________________________32
сохраняют свою справедливость (возможно, с очевидными поправками) и для систем с комплексной матрицей. В частности, например, линейная однородная система (1) с комплексной матрицей A(t) асимптотически устойчива тогда и только тогда, когда вещественная и мнимая части всякого её решения x = x(t) = u(t) +i v(t) стремятся к нулю при t → +∞.
§ 4. Устойчивость линейных систем с постоянной матрицей
Среди семейства всех дифференциальных систем класс линейных однородных систем с постоянной матрицей является наиболее простым и хорошо изученным. Как показывает следующая теорема, вопрос об устойчивости таких систем может быть сведен к решению чисто алгебраической задачи о знаках вещественных частей собственных (характеристических) значений матрицы данной системы.
10. Критерий устойчивости.
Теорема 1. Линейная однородная система
|
dx |
= A x |
(t ≥ 0) |
|
(1) |
|
dt |
|
|||
|
|
|
|
|
|
с постоянной матрицей A = (aij )n×n |
устойчива тогда и только тогда, когда |
||||
все собственные значения |
λj , |
j =1,2,..., n , матрицы |
А |
обладают |
|
неположительными вещественными частями, т.е. |
|
|
|||
|
|
Re λj |
≤ 0, j =1,2,...,n , |
|
(2) |
причём собственные значения |
λj , имеющие нулевые вещественные части, ха- |
рактеризуются тем свойством, что соответствующие им клетки Жордана сводятся к одному элементу (т.е. допускают лишь простые делители, что
равносильно выполнению равенства |
n − rang(A −λjE) = k j , где |
k j – кратность |
||
корня |
λj ). |
|
|
|
|
□ Достаточность. Обозначим |
αj +i βj , j =1,2,..., p , |
где |
i – мнимая еди- |
ница, |
собственные значения матрицы А с отрицательными вещественными |
|||
частями αj , отвечающие различным клеткам Жордана, а |
i γk , k =1,2,...,q , – |
собственные значения матрицы А с нулевыми вещественными частями, причём p +q = m – общее число клеток Жордана в нормальной форме матрицы А. То-
гда, как известно из курса дифференциальных уравнений, любое решение системы (1) может быть представлено в виде
ГЛАВА 1. УСТОЙЧИВОСТЬ ЛИНЕЙНЫХ СИСТЕМ_______________________________________________33
p |
q |
x(t)= ∑eαj t (cos(βjt) +isin (βjt)) Pj (t)+ ∑(cos(γk t) +isin (γk t)) ck , |
|
j=1 |
k=1 |
где Pj (t) – некоторые полиномиальные векторные коэффициенты, степень которых меньше кратности корня αj +i βj , а ck – постоянные вектор-столбцы.
Из выписанного представления благодаря соотношениям
eαj t P (t) |
→ 0, |
|
cos(γ |
k |
t) +i sin (γ |
k |
t) |
|
=1 |
|
|
||||||||
j |
t→+∞ |
|
|
|
|
|
|
||
|
|
|
|
|
на (t0 ,+∞). Согласно |
||||
получаем ограниченность произвольного решения x(t) |
|||||||||
теореме 1 предыдущего параграфа система (1) устойчива. |
|
|
|
|
Необходимость. Дано, что система (1) устойчива. Сначала покажем, что выполнено (2). Для этого предположим противное: найдется λs = σ + i τ – соб-
ственное значение матрицы А, для которого |
Reλs = σ > 0. Так как функция |
||||||||||||||||||||||||||||||||||
ζ = ζ(t)= eλs t c |
( |
|
|
|
c |
|
|
|
≠ 0) |
||||||||||||||||||||||||||
|
|
|
|
||||||||||||||||||||||||||||||||
является ненулевым решением системы (1), причём |
|||||||||||||||||||||||||||||||||||
|
|
|
ζ |
|
|
|
= |
|
eλs t |
|
|
|
|
|
c |
|
|
|
= eσ t |
|
|
|
c |
|
|
|
→ +∞, |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t→+∞ |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
то решение ζ неограниченно. Полученное противоречит устойчивости системы. Следовательно, (2) выполняется.
Теперь покажем, что любое собственное значение с нулевой вещественной частью обладает тем свойством, что соответствующая ему клетка Жордана сводится к одному элементу.
Пусть матрица А приведена к жордановой форме:
A = S−1 diag[J1(λ1 ),...,Jm (λm )] S,
где det(S) = |
|
S |
|
≠ 0 . Предположим |
противное: |
существует такое значение |
||||||||
|
|
|||||||||||||
λs = i μs (Re |
|
λs = 0), что ему отвечает клетка Жордана |
|
|||||||||||
|
|
|
|
|
|
|
λs |
1 |
|
0 |
0 |
|
||
|
|
|
|
|
|
|
|
0 |
λs |
|
0 |
0 |
|
|
|
|
|
|
|
(λ |
|
|
|
|
|||||
|
|
|
|
J |
|
)= |
|
|
|
|
|
|
, |
|
|
|
|
|
s |
|
s |
|
0 |
0 |
|
λs |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
0 |
0 |
|
0 |
|
|
|
|
|
|
|
|
|
|
|
λs ns ×ns |
|
ГЛАВА 1. УСТОЙЧИВОСТЬ ЛИНЕЙНЫХ СИСТЕМ_______________________________________________34
причём ns >1. Рассмотрим матрицу |
|
Θ(t)= S−1 diag[0,...,0,etJs (λs ),0,...,0] S . |
(*) |
Она является матричным решением системы (1), так как |
|
Θ& (t)= S−1 diag[0,...,0,Js (λs ) et Js (λs ),0,...,0] S =
S−1 diag[J1(λ1 ),...,Js (λs ),...,Jn (λn )] S S−1 diag[0,...,0,et Js (λs ),0,...,0] S = A Θ(t).
Из (*) следует
diag[0,...0,et Js (λs ),0,...,0]= S Θ(t) S−1 ,
откуда с использованием свойства A B ≤ AB нормы матрицы получаем
diag[0,...,0,et Js (λs ),0,...,0] = et Js (λs ) ≤ SΘ(t)S−1 .
Следовательно, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
et Js (λs ) |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Θ(t) |
|
|
|
≥ |
|
|
|
. |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
Но для всякого |
t ≥ 0 |
|
|
|
|
|
|
|
|
|
|
S |
|
|
|
S−1 |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
верно |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
|
|
|
|
|
|
|
tns −1 |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1! |
|
|
|
|
|
|
|
(ns −1)! |
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
t J |
(λ |
|
) |
|
|
λ |
t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
tns −2 |
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
e |
s |
|
s |
|
= |
e |
|
|
|
s |
|
|
0 1 |
|
|
|
|
|
|
|
|
|
|
= |
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(ns −2)! |
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
0 |
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
t |
|
|
|
|
|
|
tns −1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
1 1! |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
(ns −1)! |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
λ |
t |
|
|
|
|
|
|
|
|
|
|
|
tns −2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
σ |
t |
|
|
|
|
|
t |
|
|
|
tns −1 |
|
||||||||||||
=| e |
s |
|
| |
|
0 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= e |
|
|
s |
|
|
|
|
+ |
|
|
+... + |
|
|
|
|
|
> |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
1! |
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(ns −2)! |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(ns −1)! |
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
0 |
0 |
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(**)
(ns −1)! .