- •ГЛАВА 2. ВТОРОЙ МЕТОД ЛЯПУНОВА
- •Под вторым методом Ляпунова понимают совокупность приемов и средств исследования устойчивости решений систем дифференциальных уравнений при помощи специальных функций Ляпунова.
- •§ 1. Знакоопределенные функции
- •Рассмотрим вещественную непрерывную числовую функцию
- •Следует заметить, что в качестве W(x) из приведённого определения иногда удается использовать функцию вида
- •О п р е д е л е н и е 2. Положительно или отрицательно определенную функцию именуют знакоопределенной.
- •Пример 1. Рассмотрим функцию
- •Если |a| < 1, то она является положительно определенной, так как
- •причем
- •При a = 1 функция V1 имеет вид
- •З а м е ч а н и е 1. Если функция V = V(x) непрерывна в начале координат и V(0) = 0, то, очевидно, V допускает б.м.в.п.
- •З а м е ч а н и е 2. Для того чтобы установить, что функция V(t, x) имеет б.м.в.п., достаточно убедиться в существовании непрерывной в точке х = 0 функции W(x), удовлетворяющей следующим двум условиям:
- •Пример 2. Функция
- •на основании замечания 2 допускает б.м.в.п. в силу
- •тогда как функция
- •Упражнения
- •1) Выяснить, являются ли данные функции положительно определенными:
- •2) Определить, допускают ли указанные функции б.м.в.п.:
- •§ 2. Теорема Ляпунова об устойчивости
- •Рассмотрим систему дифференциальных уравнений
- •Будем считать, что система (1) имеет нулевое решение ξ = 0, т.е.
- •есть градиент функции V.
- •Пример 1. Пусть
- •Тогда производная функции V в силу данной системы имеет вид
- •Введем сферу
- •где произвольно выбранное ε удовлетворяет неравенствам 0 < ε < H. Очевидно, Sε – непустой компакт и функция W непрерывна на нем. Введём
- •Рассмотрим функцию
- •По условию теоремы
- •что невозможно.
- •если только ||x(t0)|| < δ. Полученное означает, что нулевое решение ξ = 0 системы (1) устойчиво по Ляпунову ■
- •Для линейной однородной системы устойчивость нулевого решения влечет ограниченность каждого решения этой системы. Поэтому имеет место следующее
- •существует положительно определенная функция V(t, x), для которой производная в силу системы неположительна, то все решения этой системы определены и ограничены при t ≥0.
- •Пример 1. Рассмотрим систему
- •и исследуем на устойчивость её нулевое решение х = 0, y = 0. Примем
- •Имеем V > 0 для всех x,y, одновременно не обращающихся в нуль, и
- •Это говорит о том, что функций, удовлетворяющих условия теоремы 1, может быть бесконечно много.
- •Проблема: построить функцию Ляпунова V(t, x), т.е. такую функцию, которая удовлетворяет условиям теоремы 1.
- •Упражнения
- •3) С помощью теоремы 1 установить устойчивость нулевого решения уравнения
- •§ 3. Теорема Ляпунова об асимптотической устойчивости
- •Вновь обратимся к системе
- •По условию теоремы
- •С учетом обозначения
- •Поэтому полученная выше оценка сверху для v(t) принимает вид
- •Отсюда при достаточно большом t следует
- •В итоге установлено равенство
- •Следовательно, предельное соотношение (#) истинно ■
- •Следствие 1. Если для линейной однородной системы
- •Заметим, что равенства нулю всех частных производных функции V вытекают из условия положительной определённости и предположения о непрерывной дифференцируемости этой функции в начале координат (в этом можно легко убедиться, предположив противное).
- •Рассмотрим систему
- •Пример 2. Исследуем на устойчивость нулевое решение уравнения
- •Если a = 0, то, очевидно, нулевое решение устойчиво, но не является асимптотически устойчивым.
- •где m+1 – четное. Следовательно, если a < 0, то функция V удовлетворяет всем требованиям последней теоремы, а значит нулевое решение асимптотически устойчиво.
- •Упражнения
- •1) С помощью надлежащей функции Ляпунова установить асимптотическую устойчивость нулевого решения уравнения
- •2) Подобрать функцию Ляпунова и с её помощью установить асимптотическую устойчивость нулевого решения системы
- •§ 4. Теоремы Ляпунова о неустойчивости
- •10. Первая теорема Ляпунова о неустойчивости.
- •Теорема 1. Пусть для системы
- •1) V(t,x) допускает б.м.в.п.;
- •Тогда нулевое решение системы (1) неустойчиво.
- •□ Согласно условию 3) имеем
- •Пусть
- •Это означает, что нулевое решение системы (1) неустойчиво по Ляпунову ■
- •20. Вторая теорема Ляпунова о неустойчивости.
- •Упражнения
- •2) Проверить, что указанные функции являются функциями Ляпунова для систем
- •имеющих неустойчивое нулевое решение.
- •§ 5. Экспоненциальная устойчивость
- •выполнено неравенство
- •Как показывает следующий результат, экспоненциальная и асимптотическая устойчивость в классе линейных однородных систем с постоянной матрицей совпадают.
- •Лемма 2. Линейная однородная система
- •с постоянной матрицей А асимптотически устойчива тогда и только тогда, когда эта система экспоненциально устойчива (т. е. экспоненциально устойчивым является всякое её решение).
- •Всякое решение линейной системы (4) может быть представлено в виде
- •Введём числа
- •Очевидно,
- •Из курса линейной алгебры известно, что
- •Далее
- •А значит,
- •с непрерывной на промежутке [0, +∞) матрицей A(t). Для таких систем нередко вводят несколько более «жёсткое» понятие экспоненциальной устойчивости, чем это было сделано выше.
- •Нетрудно понять, что из экспоненциальной устойчивости системы (5) вытекает её асимптотическая устойчивость (но не обратно). В случае, когда матрица A(t) = A постоянна, оба указанных понятия эквивалентны.
- •В следующей теореме представлен критерий экспоненциальной устойчивости системы (5).
- •Поскольку
- •Рассмотрим подынтегральное выражение из (9):
- •Используя условие 1), для него можно записать
- •В таком случае из равенства (9) вытекают неравенства
- •Аналогичным образом поступим и с правым неравенством:
- •Из последних двух неравенств легко вывести (6) ■
- •которое именуют матричным уравнением Ляпунова.
- •имела единственное решение относительно P и чтобы это решение представляло собой симметричную положительно определённую матрицу.
- •Упражнения
- •1) Проверить, являются ли экспоненциально устойчивыми нулевые решения следующих уравнений
- •§ 6. Устойчивость квазилинейных систем
- •Линейная система
- •называется системой линейного приближения системы (1).
- •то нулевое решение системы (1) асимптотически устойчиво по Ляпунову.
- •Рассмотрим функцию
- •С учётом (*) для нее получаем представление
- •Получим
- •Напомним, что x – вектор-столбец. Благодаря
- •имеем
- •если только
- •Из доказанной теоремы вместе с теоремой 1 предыдущего параграфа получаем следующее утверждение.
- •Следствие 1. Рассмотрим нелинейную систему
- •Следствие 2. Рассмотрим нелинейную автономную систему
- •и можно применить доказанную выше теорему 1 ■
- •20. Теорема Ляпунова о неустойчивости квазилинейной системы.
- •Теорема 2. Пусть дана нелинейная система
- •□ Пусть
- •Выполним замену
- •где вектор y, вообще говоря, комплексный. Имеем
- •Поэтому
- •Положим
- •и перепишем систему (*) следующим образом
- •Переходим к комплексно-сопряженной системе:
- •Очевидно,
- •Поэтому для функции
- •получаем представление
- •Полагая
- •Если теперь вернуться к исходной переменной, то получим
- •откуда
- •Пример 1. Исследовать на устойчивость нулевое решение системы второго порядка
- •где a – const. Воспользовавшись формулой Маклорена, найдем разложения правых частей системы и перепишем её в виде квазилинейной системы
- •Упражнения
- •Исследовать на асимптотическую устойчивость нулевые решения следующих систем:
- •§ 7. Теорема Зубова
- •Рассмотрим автономную дифференциальную систему
- •в которой
- •Относительно вектор-функции f будем предполагать следующее:
- •□ Ограничимся доказательством лишь части «достаточность» (доказательство второй части можно найти в [7]). По условию функция V(x) – отрицательно определена, а её производная по времени в силу системы (1) положительно определена, так как
- •С этой целью выполним замену, осуществляющую изменение переменной времени по правилу
- •где правые части уравнений непрерывны и ограничены относительно x, так как
- •вытекающим из условия 3) доказываемой теоремы ■
- •З а м е ч а н и е 1. Уравнение в частных производных
- •из последнего условия теоремы 1 называют уравнением Зубова.
- •З а м е ч а н и е 2. На основе теоремы 1 разработан метод Зубова построения оценки области асимптотической устойчивости, общая идея которого состоит в следующем. Пусть система (1) представима в виде
- •где A − постоянная матрица,
- •Существует вычислительный алгоритм построения оценок области притяжения для случая, когда в качестве V(x) выбрана положительно определенная квадратичная форма.
- •2) отрицательно определена;
- •4) удовлетворяет равенству
- •Как видим, все условия теоремы Зубова выполнены, поэтому множество A является областью притяжения нулевого решения.
- •Упражнения
- •1) Построить область асимптотической устойчивости для уравнения
- •2) Построить область асимптотической устойчивости для системы
ГЛАВА 2. ВТОРОЙ МЕТОД ЛЯПУНОВА_______________________________________________________98
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
b1 |
≤ |
W |
≤ |
b2 |
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a2 |
V |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
В таком случае из равенства (9) вытекают неравенства |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
− |
b2 |
(t−t0 ) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
− |
b1 |
(t−t0 ) |
|
|
|
|
|||||||||||
V(t |
|
, x |
|
|
|
≤ V (t; t |
|
, x |
|
) ≤ |
V(t |
|
, x |
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||
0 |
0 |
) e a1 |
|
0 |
0 |
0 |
0 |
) e |
a2 |
. |
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
Заменим V(t0, x0) |
в левом неравенстве на |
a1 |
|
|
|
x0 |
|
|
|
2 , а V1 – на |
a2 |
|
|
|
x0 |
|
|
|
2 : |
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
− |
b2 |
(t−t |
0 |
) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 . |
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
a1 || x0 ||2 e |
a1 |
≤ a2 |
|
|
x(t) |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||
Аналогичным образом поступим и с правым неравенством: |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 ≤ a2 || x0 ||2 e |
− |
b1 |
(t−t |
0 |
) |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
|
|
|
|
a1 |
|
|
|
x(t) |
|
|
|
|
a2 |
|
|
|
. |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Из последних двух неравенств легко вывести (6) ■
30. Матричное уравнение Ляпунова. Перепишем уравнение (7) с учётом
равенств V(t,x) = xT P(t) x и |
W(t,x) = xT Q(t) x : |
|
|
|
|||||||||
dV |
= |
dxT |
P(t) x +x |
T |
|
dP(t) |
x + x |
T |
P(t) |
dx |
= −x |
T |
Q(t) x . |
dt |
dt |
|
dt |
|
dt |
|
|||||||
|
|
|
|
|
|
|
|
|
|
||||
Поскольку x = x(t) |
– решение системы (5), имеем |
|
|
|
|
||||||||
dV |
= xT AT (t) P(t) + dP(t) +P(t) A(t) |
x = −xT Q(t) x . |
|||||||||||
dt |
|
|
|
|
|
dt |
|
|
|
|
|
|
|
Отсюда благодаря произвольности выбора x = x(t) следует равенство |
|
||
|
dP(t) |
+ AT (t) P(t) + P(t) A(t) = −Q(t) , |
(10) |
|
dt |
||
|
|
|
которое именуют матричным уравнением Ляпунова.
Оказывается, в частном случае, когда A – постоянная матрица, в качестве матрицы Q также можно выбрать постоянную положительно определённую матрицу1. Тогда уравнение (10) будет иметь единственное решение в виде постоянной матрицы P, причём эта матрица будет положительно определённой.
В соответствии с этим, для асимптотической устойчивости линейной однородной системы (5) с постоянной матрицей A необходимо и достаточ-
1 Напоминаем, что симметричная матрица положительно определена, если таковой же является квадратичная форма с данной матрицей.
ГЛАВА 2. ВТОРОЙ МЕТОД ЛЯПУНОВА_______________________________________________________99
но, чтобы для произвольно выбранной постоянной положительно определённой матрицы Q линейная алгебраическая система
AT P +P A = −Q
имела единственное решение относительно P и чтобы это решение представляло собой симметричную положительно определённую матрицу.
Упражнения
1)Проверить, являются ли экспоненциально устойчивыми нулевые решения следующих уравнений
i)x& = −xt 14 ,
ii)x& = −xsin2 t ,
iii)x& = −x ln(1+t)1+t .
§ 6. Устойчивость квазилинейных систем
Существует класс нелинейных дифференциальных систем, правые части которых близки к правым частям некоторых линейных систем в окрестности положения равновесия. В таком случае оказывается, что с точки зрения устойчивости в обоих классах систем положения равновесия ведут себя сходным образом. Соответствующие результаты формулируются ниже.
10. Теорема Ляпунова об устойчивости квазилинейной системы. Рас-
смотрим нелинейную дифференциальную систему
|
dx |
= A x +ϕ(t, x) |
|
|
|
|
|
|
|
|
(t ≥ 0) , |
(1) |
|
|
|||||||||||
|
dt |
|
|
|
|
|
|
|
|
|
|
|
где А – постоянная матрица размера n × n и векторная функция ϕ(t, x) |
опре- |
|||||||||||
делена и непрерывна на множестве {(t,x) | t ≥ 0, |
|
|
|
x |
|
|
|
≤ H}. |
|
|||
|
|
|
|
|
||||||||
Будем предполагать, что ϕ(t, 0) = 0 |
для всех t ≥ 0, т.е. система (1) имеет |
|||||||||||
нулевое решение (положение равновесия). |
|
|
|
|
|
|
|
|
|
|
||
Линейная система |
|
|
|
|
|
|
|
|
|
|
ddtx = A x
называется системой линейного приближения системы (1).
ГЛАВА 2. ВТОРОЙ МЕТОД ЛЯПУНОВА_______________________________________________________100
Теорема 1. Предположим, что функция ϕ(t, x) удовлетворяет условию
||ϕ( |
t, |
x) || →→ 0 |
|
|
|
равномерно по t ≥ 0 при x →0 . |
(2) |
|||||||||||||||
|
x |
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Если при этом все собственные значения λj, |
j = 1,2,...,n , матрицы А имеют |
|||||||||||||||||||||
отрицательные вещественные части, т. е. |
|
|
||||||||||||||||||||
|
|
|
|
|
|
Re λj |
< 0, j = 1,2,...,n , |
|
||||||||||||||
то нулевое решение системы (1) асимптотически устойчиво по Ляпунову. |
|
|||||||||||||||||||||
□ Через ξ(t; x) |
обозначим решение соответствующей линейной систе- |
|||||||||||||||||||||
мы: |
|
|
|
|
|
|
|
|
dξ |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
= A ξ |
|
(3) |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
dt |
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
с начальным условием |
(0,x) , т.е. |
ξ(0;x) = x Rn . |
|
|||||||||||||||||||
Пусть K(t) – нормированная фундаментальная матрица (матрицант) сис- |
||||||||||||||||||||||
темы (2), K(0) = E . Имеем |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
ξ(t; x) = K(t) x . |
|
(*) |
||||||||||||||
В силу отрицательности вещественных частей собственных значений λj |
при |
|||||||||||||||||||||
некотором N > 0 выполнено |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
||K(t)|| ≤ N e−αt |
t ≥ 0 , |
|
||||||||||||||
где max Re λj < −α < 0 . Отсюда |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
j |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
ξ(t;x) |
|
|
|
≤ N e−αt |
|
|
|
x |
|
|
|
|
t ≥ 0 . |
(**) |
||
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
Рассмотрим функцию |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
+∞ |
|
|
||||||||
|
|
|
|
|
|
V(x) = ∫ || ξ(τ; x)||2dτ. |
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
С учётом (*) для нее получаем представление |
|
|
||||||||||||||||||||
+∞ |
|
|
|
|
|
|
|
|
|
+∞ |
|
|
||||||||||
V(x) = ∫ K(τ) x, K(τ) x dτ = ∫ KT (τ) K(τ) x, x dτ = S x, x , |
|
|||||||||||||||||||||
0 |
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
ГЛАВА 2. ВТОРОЙ МЕТОД ЛЯПУНОВА_______________________________________________________101
+∞
где S := ∫ KT (τ) K(τ)dτ.
0
Таким образом, V(x) – квадратичная форма относительно переменных
x1 |
, ... , xn с вещественной симметричной матрицей S = (Sjk ) |
. |
|
|
n x n |
|
В силу (**) интеграл, определяющий V, сходится. Поэтому квадратичная |
форма V определена на Rn и конечна, причём на основании единственности решения линейной системы (2) верно
V(x) > 0 |
x ≠ 0, |
V(0) = 0 . |
|
Воспользуемся групповым свойством ξ(t;ξ(τ;x)) = ξ(t + τ;x) решений |
автономной системы (2), которое выражает тот факт, что состояние системы ξ(t+τ; x) (т.е. правая часть написанного выше равенства) к моменту времени t + τ, полученное перемещением вдоль траектории из начальной точки x = ξ(0; x), может быть достигнуто движением вдоль этой траектории из «промежуточной» точки ξ(τ;ξ(0; x)) за время t (см. рис. 2.4).
y2
ξ(τ;ξ(0; x))
ξ(t+τ; ξ(0; x)))
ξ(0; x) = x
y1
Рис. 2.4. Иллюстрация группового свойство системы.
Получим
+∞ |
+∞ |
+∞ |
V(ξ(t; x)) = ∫ || ξ(τ; ξ(t; x)) ||2dτ = ∫ || ξ(t + τ;x) ||2dτ = ∫ || ξ(τ;x)||2dτ. |
||
0 |
0 |
t |
Дифференцируя V(ξ(t; x)) |
по t в силу системы (2) в точке x, находим |