- •Пример 1. Рассмотрим стандартный полином
- •З а м е ч а н и е 1 (о распределении знаков полиномов h и g). Полиномы h и g имеют следующие знаки:
- •Упражнения
- •С помощью критерия Михайлова исследовать на устойчивость следующие уравнения:
- •§ 8. Теорема Харитонова
- •О п р е д е л е н и е 1. Интервальным полиномом называется множество полиномов вида
- •О п р е д е л е н и е 2. Интервальный полином называется полиномом Гурвица, если все полиномы, образующие его, являются полиномами Гурвица.
- •являлись полиномами Гурвица.
- •□ Необходимость очевидна.
- •Достаточность. Рассмотрим годограф Михайлова
- •для полинома
- •Введем полиномы
- •Тем самым,
- •Анализ доказательства теоремы 1 дает возможность сформулировать следующее утверждение.
- •Упражнения
- •§ 9. Лемма Гронуолла-Беллмана и ее обобщения
- •Результаты данного параграфа носят вспомогательный характер; они будет использованы далее.
- •Отсюда следует неравенство
- •которое можно переписать в виде
- •Тогда
- •□ Из (3) следует
- •Применяя лемму Гронуолла-Беллмана, получаем
- •Аналогично из (3) вытекает
- •Далее
- •где m > 1 и c – некоторая положительная константа.
- •Тогда имеет место неравенство
- •□ Из неравенства (5) следует
- •Отсюда с учётом (5) получаем
- •Интегрируя полученное неравенство, находим
- •а значит
- •где благодаря (6) выражение, записанное в квадратных скобках, строго положительно. Тогда
- •§ 10. Устойчивость линейных систем с переменной матрицей
- •Теорема 1. Пусть линейная система
- •□ Воспользуемся формулой Коши
- •Отсюда
- •Применяя лемму Гронуолла-Беллмана и используя условие (3), получаем
- •Пример 1. Рассмотрим уравнение
- •Оно эквивалентно системе
- •Характеристическое уравнение для матрицы соответствующей линейной системы имеет вид
- •устойчива. Далее,
- •а значит
- •Поэтому почти линейная система (*) также является устойчивой.
- •В теореме 1 речь шла об устойчивости линейной системы с почти постоянной матрицей. Вопрос об асимптотической устойчивости указанной системы решается в следующем утверждении.
- •Теорема 2. Если линейная однородная система
- •с постоянной матрицей А асимптотически устойчива, то таковой же будет и возмущенная линейная система
- •Запишем формулу Коши для системы (5):
- •Отсюда,
- •Но, как известно [2],
- •где с = с(ε) – некоторая положительная константа. Поэтому
- •Существует T > 0, для которого
- •С учетом полученного, из (**) следует
- •с постоянными матричными коэффициентами Аm,…,А0 асимптотически устойчива, если все корни характеристического уравнения det(λE – Am) = 0 имеют отрицательные вещественные части.
- •□ Полагая
- •исходную систему приводим к виду
- •З а м е ч а н и е 2. Для линейной системы с переменной матрицей А теорема 2, вообще говоря, неверна. Это подтверждает следующий пример.
- •Пример 2. Пусть
- •20. Неравенство Важевского. Это неравенство является полезным инструментом при исследовании устойчивости решений некоторых линейных однородных дифференциальных систем с переменной матрицей.
- •Теорема 3. Для всякого решения линейной системы
- •имеет место неравенство Важевского
- •Поэтому с учётом (*) получаем
- •то система (6) асимптотически устойчива.
- •Пример 3. Исследуем на асимптотическую устойчивость систему
- •Здесь
- •Отсюда
- •а значит
- •30. Необходимые условия устойчивости и асимптотической устойчивости. Обозначим через X(t) фундаментальную матрицу линейной системы (4) и выпишем известную из курса дифференциальных уравнений формулу Лиувилля
- •Непосредственно из этой формулы при помощи рассуждений от противного вытекают следующие два утверждения.
- •Теорема 5. Для асимптотической устойчивости линейной системы (6) необходимо, чтобы
- •З а м е ч а н и е 3. Примеры показывают, что необходимые условия теорем 4 и 5 не являются достаточными.
- •Пример 4. Рассмотрим линейную нестационарную систему
- •Согласно теореме 4 отсюда следует, что исходная система неустойчива.
- •Упражнения
- •1) Исследовать на асимптотическую устойчивость систему
- •2) Установить асимптотическую устойчивость системы
- •3) Установить асимптотическую устойчивость системы
- •4) Установить неустойчивость системы
ГЛАВА 1. УСТОЙЧИВОСТЬ ЛИНЕЙНЫХ СИСТЕМ_______________________________________________67
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
u(t) ≤ c + ∫f(τ) um (τ)]d |
τ |
t ≥ t0 |
(5) |
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
t0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
t |
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∫ |
f(τ)dτ < |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t ≥ t0 , |
(6) |
|||||||
|
|
|
|
|
(m −1) cm−1 |
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
t0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где m > 1 и c – некоторая положительная константа. |
|
|||||||||||||||||||||||||||
Тогда имеет место неравенство |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
u(t) ≤ |
|
|
|
|
|
|
|
|
|
c |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t ≥ t0 . |
(7) |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
|
m−1 |
|
|
|
|
||||||
|
|
1 |
|
−(m −1) cm-1 ∫f(τ)dτ |
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
□ Из неравенства (5) следует |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
m |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∫ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
[u(t)]m ≤ c |
+ |
|
|
f(τ) um |
(τ)]dτ |
|
|
|
|
|
t ≥ t0 . |
|
||||||||||||||||
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
140 4424443 |
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
=ω(t)>0 |
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Отсюда с учётом (5) получаем |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
f(t) um (t) |
|
= |
ω′(t) |
|
|
|
≤ f(t) |
|
|
|
|
|
t ≥ t0 . |
|
||||||||||||
|
|
|
ωm (t) |
|
|
ωm(t) |
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
Интегрируя полученное неравенство, находим |
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
t |
′ |
|
|
|
|
|
|
|
ω(t) |
|
|
dω |
t |
|
|
|
|
|
||||||
|
|
|
|
∫ |
ω(τ)dτ |
= |
∫ |
|
|
|
|
∫f(τ)dτ. |
|
|||||||||||||||
|
|
|
|
|
m |
|
|
|
|
|
|
|
|
|
|
≤ |
|
|||||||||||
|
|
|
ω |
(τ) |
|
|
|
|
|
ω |
m |
|
||||||||||||||||
|
|
|
|
t0 |
|
|
|
|
ω(t0 ) |
|
|
|
|
t0 |
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
В силу m > 1 и ω(t 0 ) = c |
имеем |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
ω |
1−m |
|
ω(t) |
= ω |
1−m |
(t) − |
|
|
|
1 |
|
|
|
|
|
|
t |
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
≤ ∫f(τ)dτ, |
|
|||||||||||||
1 −m |
|
|
|
|
|
(1 −m) c |
m−1 |
|
||||||||||||||||||||
|
ω(t0 ) |
|
1 −m |
|
|
|
|
|
|
t0 |
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
а значит
ГЛАВА 1. УСТОЙЧИВОСТЬ ЛИНЕЙНЫХ СИСТЕМ_______________________________________________68
|
1 |
|
1−m |
|
|
|
m−1 |
|
t |
|
|
|
|
≥ c |
1 |
−(m −1) |
c |
|
∫f(τ)dτ , |
||
ω |
m−1 |
(t) |
|
|||||||
|
|
|
|
|
|
t |
0 |
|
||
|
|
|
|
|
|
|
|
|
|
где благодаря (6) выражение, записанное в квадратных скобках, строго положительно. Тогда
cm−1
ωm−1(t) ≤ 1 −(m −1) cm−1 ∫tt0 f(τ)dτ .
Теперь, извлекая из обеих частей данного неравенства корень положительной степени m – 1, с учётом неравенства u(t) ≤ ω(t) для всех t ≥ t0 легко прийти к
требуемому результату (7) ■
§ 10. Устойчивость линейных систем с переменной матрицей
10. Устойчивость линейных систем с почти постоянной матрицей. На основе леммы Гронуолла-Беллмана здесь будет показано, что при определённых предположениях устойчивость (а также асимптотическая устойчивость) линейной однородной дифференциальной системы, матрица которой «незначительно отличается» от некоторой постоянной матрицы A (т.е. является почти постоянной), вытекает из устойчивости (соответственно – из асимптотической устойчивости) линейной системы с матрицей A.
Теорема 1. Пусть линейная система
dx |
= A x |
(t ≥ 0) |
(1) |
|
dt |
||||
|
|
|
где A – постоянная матрица размера n ×n , устойчива. Тогда устойчивой будет и система
|
dy |
= (A + B(t)) y |
(t ≥ 0) |
(2) |
|
dt |
|||
|
|
|
|
|
где матрица B(t) C[0,+∞) такова, что |
|
|
||
|
|
+∞ |
|
|
|
|
∫|| B(t) || dt < +∞. |
|
(3) |
|
|
0 |
|
|
□ Воспользуемся формулой Коши
ГЛАВА 1. УСТОЙЧИВОСТЬ ЛИНЕЙНЫХ СИСТЕМ_______________________________________________69
|
t |
|
|
|
|
|
y(t) = et A y(0) + ∫e(t −τ) A B(τ) y(τ)dτ |
t ≥ 0 . |
|||||
Отсюда |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
|
|
|| y(t) ||≤|| et A || || y(0) || +∫|| e(t −τ) A || || B(τ) || || y(τ) ||dτ |
t ≥ 0 . |
|||||
|
0 |
|
|
|
|
|
Так как система (1) устойчива, то её фундаментальная матрица |
et A ограниче- |
|||||
на, т.е. существует константа |
K > 0, при которой || et A ||≤ K |
для всех t ≥ 0 . |
||||
Поэтому |
t |
|
|
|
|
|
|
|
|
|
|
|
|
|| y(t) ||≤ K || y(0) || +∫K || B(τ) || || y(τ) ||dτ |
t ≥ 0 . |
|||||
|
0 |
|
|
|
|
|
Применяя лемму Гронуолла-Беллмана и используя условие (3), получаем |
||||||
|
|
|
t |
|
|
|
|| y(t) ||≤ K |
|
|
|
|
≤ |
|
|| y(0) || exp K ∫|| B( τ) || dτ |
|
|||||
|
|
|
0 |
|
|
|
|
|
+∞ |
|
|
|
|
|
|
|
|
|
|
|
≤ K || y(0) || exp K ∫|| B( τ) || dτ < +∞, |
|
|||||
|
|
0 |
|
|
|
|
что означает ограниченность каждого решения y(t) системы (2) на полуоси [0;+∞) . Значит, эта система устойчива по Ляпунову ■
Как показал О. Перрон, для переменной матрицы A = A(t) доказанная теорема в общем случае не верна.
Пример 1. Рассмотрим уравнение
&& |
|
2 |
|
b |
|
(t ≥ 0). |
+(a |
|
+ |
(t +1)2 ) x = 0 |
|||
x |
|
Оно эквивалентно системе
x& = y,
y& = −(a 2 + (t +b1)2 )x,
или
|
|
|
|
|
|
|
|
|
|
|
|
|
=A |
|
|
|
=B |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
6447448 |
|
|
|||||
|
|
|
|
|
|
0 |
|
1 |
|
64748 |
|
|
0 |
0 |
|
|
|||||
|
& |
|
|
|
|
x |
0 |
|
1 x |
|
x |
|
|||||||||
x |
|
|
|
|
b |
|
|
|
|
|
|
b |
|
|
|||||||
|
|
= |
|
−a2 |
− |
|
|
0 |
|
|
= |
−a |
2 |
+ |
|
− |
0 |
. |
(*) |
||
|
|
(t +1) |
|
(t +1)2 |
|||||||||||||||||
y |
|
|
|
|
2 |
|
|
y |
|
|
0 y |
|
|
|
y |
|
|||||
|
& |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ГЛАВА 1. УСТОЙЧИВОСТЬ ЛИНЕЙНЫХ СИСТЕМ_______________________________________________70
Характеристическое уравнение для матрицы соответствующей линейной системы имеет вид
det(λΕ− Α) = |
|
λ |
−1 |
|
= λ2 |
+a 2 |
= 0 . |
|
|
||||||
|
|
a2 |
λ |
|
|
|
|
Находим его корни λ1,2
устойчива. Далее,
B =
а значит
= ±a i . Следовательно, линейная система
xy&& =А xy
|
0 |
0 |
|
| b | |
|
|
|
b |
|
|
|
||
− |
0 |
= (t +1)2 |
t ≥ 0 , |
|||
(t +1)2 |
+∞
∫B(t)dt < +∞.
0
Поэтому почти линейная система (*) также является устойчивой.
В теореме 1 речь шла об устойчивости линейной системы с почти постоянной матрицей. Вопрос об асимптотической устойчивости указанной системы решается в следующем утверждении.
Теорема 2. Если линейная однородная система
dx |
= Α x |
(t ≥ 0) |
|
dt |
|||
|
|
с постоянной матрицей А асимптотически устойчива, то таковой же будет и возмущенная линейная система
|
|
dy |
= (A + B(t)) y |
(t ≥ 0) , |
(5) |
|
|
|
|||
|
|
dt |
|
|
|
где матрица В(t) C[0,+∞) |
такова, что B(t) → 0 . |
|
|||
|
|
|
|
t→+∞ |
|
■ В силу асимптотической устойчивости системы (1) собственные значе- |
|||||
ния λj(A) |
матрицы А |
обладают отрицательными вещественными частями. |
|||
Положим |
α = max Reλj(A) < 0 и выберем ε > 0 настолько малым, чтобы было |
||||
|
j |
|
|
|
|
выполнено |
α+2ε < 0 . |
|
|
|
|